Target mimics: an embedded layer of microRNA-involved gene regulatory networks in plants

The stability of transcription factor PfSPL1 participates in the response to phytoplasma stress in Paulownia fortunei

The current understanding of the pathogenesis of phytoplasma is still very limited and challenging. Here, ceRNA regulatory network and degradome sequencing identified a PfmiR156f-PfSPL regulatory module in Paulownia fortunei infected by phytoplasma, and RLM-5'RACE and dual luciferase analyses verified the relationship. The PfmiR156 cleavage site was located at 1104 nt and 1177 nt of PfSPL1 and PfSPL10, respectively. MG132 and epoxomicin, two 26S proteasome inhibitors, significantly increased the accumulation of PfSPL1. PfSPL1 was also the attack target of phytoplasma effectors (Pawb 3/9/16/37/51) after the phytoplasma invaded Paulownia. Moreover, molecular docking implied that the effectors may interact with the conserved SBP domain of the target protein PfSPL1. Basically, these results indicated that the stability of PfSPL1 was regulated by PfmiR156 cleavage activity and/or the 26S proteasome pathway at the post-translation level. The PfSPL1, which is a transcription factor, was also the one of the targets of multiple effectors attacking Paulownia. This study provides a good scope to understand the paulownia phytoplasma infecting mechanism.

Whole transcriptome analysis and construction of a ceRNA regulatory network related to leaf and petiole development in Chinese cabbage (Brassica campestris L. ssp. pekinensis)

BACKGROUND

The growth and development of leaves and petioles have a significant effect on photosynthesis. Understanding the molecular mechanisms underlying leaf and petiole development is necessary for improving photosynthetic efficiency, cultivating varieties with high photosynthetic efficiency, and improving the yield of crops of which the leaves are foodstuffs. This study aimed to identify the mRNAs, long non-coding RNAs (lncRNAs), microRNAs (miRNAs), and circular RNAs (circRNAs) related to leaf and petiole development in Chinese cabbage (Brassica campestris L. ssp. pekinensis). The data were used to construct a competitive endogenous RNA (ceRNA) network to obtain insights into the mechanisms underlying leaf and petiole development.

RESULTS

The leaves and petioles of the 'PHL' inbred line of Chinese cabbage were used as research materials for whole transcriptome sequencing. A total of 10,646 differentially expressed (DE) mRNAs, 303 DElncRNAs, 7 DEcircRNAs, and 195 DEmiRNAs were identified between leaves and petioles. Transcription factors and proteins that play important roles in leaf and petiole development were identified, including xyloglucan endotransglucosylase/hydrolase, expansion proteins and their precursors, transcription factors TCP15 and bHLH, lateral organ boundary domain protein, cellulose synthase, MOR1-like protein, and proteins related to plant hormone biosynthesis. A ceRNA regulatory network related to leaf and petiole development was constructed, and 85 pairs of ceRNA relationships were identified, including 71 DEmiRNA-DEmRNA, 12 DEmiRNA-DElncRNA, and 2 DEmiRNA-DEcircRNA pairs. Three LSH genes (BrLSH1, BrLSH2 and BrLSH3) with significant differential expression between leaves and petioles were screened from transcriptome data, and their functions were explored through subcellular localization analysis and transgenic overexpression verification. BrLSH1, BrLSH2 and BrLSH3 were nuclear proteins, and BrLSH2 inhibited the growth and development of Arabidopsis thaliana.

CONCLUSIONS

This study identifies mRNAs and non-coding RNAs that may be involved in the development of leaves and petioles in Chinese cabbage, and establishes a ceRNA regulatory network related to development of the leaves and petioles, providing valuable genomic resources for further research on the molecular mechanisms underlying leaf and petiole development in this crop species.

Expression Profiles and Characteristics of Apple lncRNAs in Roots, Phloem, Leaves, Flowers, and Fruit

LncRNAs impart crucial effects on various biological processes, including biotic stress responses, abiotic stress responses, fertility and development. The apple tree is one of the four major fruit trees in the world. However, lncRNAs's roles in different tissues of apple are unknown. We identified the lncRNAs in five tissues of apples including the roots, phloem, leaves, flowers, and fruit, and predicted the intricate regulatory networks. A total of 9440 lncRNAs were obtained. LncRNA target prediction revealed 10,628 potential lncRNA-messenger RNA (mRNA) pairs, 9410 pairs functioning in a cis-acting fashion, and 1218 acting in a trans-acting fashion. Functional enrichment analysis showed that the targets were significantly enriched in molecular functions related to photosynthesis-antenna proteins, single-organism metabolic process and glutathione metabolism. Additionally, a total of 88 lncRNAs have various functions related to microRNAs (miRNAs) as miRNA precursors. Interactions between lncRNAs and miRNAs were predicted, 1341 possible interrelations between 187 mdm-miRNAs and 174 lncRNAs (1.84%) were identified. MSTRG.121644.5, MSTRG.121644.8, MSTRG.2929.2, MSTRG.3953.2, MSTRG.63448.2, MSTRG.9870.2, and MSTRG.9870.3 could participate in the functions in roots as competing endogenous RNAs (ceRNAs). MSTRG.11457.2, MSTRG.138614.2, and MSTRG.60895.2 could adopt special functions in the fruit by working with miRNAs. A further analysis showed that different tissues formed special lncRNA-miRNA-mRNA networks. MSTRG.60895.2-mdm-miR393-MD17G1009000 may participate in the anthocyanin metabolism in the fruit. These findings provide a comprehensive view of potential functions for lncRNAs, corresponding target genes, and related lncRNA-miRNA-mRNA networks, which will increase our knowledge of the underlying development mechanism in apple.

Functional analysis of the eTM-miR171-SCL6 module regulating somatic embryogenesis in Lilium pumilum DC. Fisch

Somatic embryogenesis (SE) is of great significance in Lilium bulb production, germplasm preservation and genetic improvement. miRNAs are important regulators of plant growth and development at the transcriptional level. Previous research by our group has shown that lpu-miR171 and its target gene SCARECROW-LIKE 6 (SCL6) play an important regulatory role in lily SE, and we predicted and identified that endogenous target mimics (eTMs) can regulate lpu-miR171. However, the associated mechanism and internal regulatory network are not yet clear. In the present study, lpu-miR171 was used as an entry point to explore the regulatory network between its upstream eTMs and its downstream target gene LpSCL6, as well as to identify the mechanism of this regulatory network in Lilium SE. Tobacco transient transformation confirmed that miRNA171 significantly inhibited the expression of LpSCL6. On this basis, the Lilium stable genetic transformation system was used to demonstrate that silencing lpu-miR171a and lpu-miR171b and overexpressing LpSCL6-II and LpSCL6-I promoted starch accumulation in calli and the expression of key cell cycle genes, thus providing energy to meet preconditions for SE and accelerate the formation and development of Lilium somatic embryos. LpSCL6-II and LpSCL6-I are nuclear proteins with self-activation activity in yeast cells. In addition, we confirmed in Lilium that lpu-eTM171 is the eTM of lpu-miR171 that binds lpu-miR171 to prevent cleavage of the target gene LpSCL6, thereby promoting SE. Therefore, the present study established a new mechanism whereby the eTM-miR171-SCL6 module regulates SE in Lilium pumilum DC. Fisch. and provided new insights clarifying the mechanism of SE.

Whole-transcriptome analysis and construction of an anther development-related ceRNA network in Chinese cabbage (Brassica campestris L. ssp. pekinensis)

Anther development is precisely regulated by a complex gene network, which is of great significance to plant breeding. However, the molecular mechanism of anther development in Chinese cabbage is unclear. Here, we identified microRNAs (miRNAs), mRNAs, long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs) related to anther development in Chinese cabbage (Brassica campestris L. ssp. pekinensis) to construct competitive endogenous RNA (ceRNA) regulatory networks and provide valuable knowledge on anther development. Using whole-transcriptome sequencing, 9055, 585, 1344, and 165 differentially expressed mRNAs (DEmRNAs), miRNAs (DEmiRNAs), lncRNAs (DElncRNAs), and circRNAs (DEcircRNAs) were identified, respectively, in the anthers of Chinese cabbage compared with those in samples of the vegetative mass of four true leaves. An anther-related ceRNA regulatory network was constructed using miRNA targeting relationships, and 450 pairs of ceRNA relationships, including 97 DEmiRNA-DEmRNA, 281 DEmiRNA-DElncRNA, and 23 DEmiRNA-DEcircRNA interactions, were obtained. We identified important genes and their interactions with lncRNAs, circRNAs, and miRNAs involved in microsporogenesis, tapetum and callose layer development, pollen wall formation, and anther dehiscence. We analyzed the promoter activity of six predominant anther expression genes, which were expressed specifically in the anthers of Arabidopsis thaliana, indicating that they may play an important role in anther development of Chinese cabbage. This study lays the foundation for further research on the molecular mechanisms of anther growth and development in Chinese cabbage.

Bm-miR172c-5p Regulates Lignin Biosynthesis and Secondary Xylem Thickness by Altering the Ferulate 5 Hydroxylase Gene in Bacopa monnieri

MicroRNAs (miRNAs) are small non-coding, endogenous RNAs containing 20-24 nucleotides that regulate the expression of target genes involved in various plant processes. A total of 1,429 conserved miRNAs belonging to 95 conserved miRNA families and 12 novel miRNAs were identified from Bacopa monnieri using small RNA sequencing. The Bm-miRNA target transcripts related to the secondary metabolism were further selected for validation. The Bm-miRNA expression in shoot and root tissues was negatively correlated with their target transcripts. The Bm-miRNA cleavage sites were mapped within the coding or untranslated region as depicted by the modified RLM-RACE. In the present study, we validate three miRNA targets, including asparagine synthetase, cycloartenol synthase and ferulate 5 hydroxylase (F5H) and elucidate the regulatory role of Bm-miR172c-5p, which cleaves the F5H gene involved in the lignin biosynthesis. Overexpression (OE) of Bm-miR172c-5p precursor in B. monnieri suppresses F5H gene, leading to reduced lignification and secondary xylem thickness under control and drought stress. By contrast, OE of endogenous target mimics (eTMs) showed enhanced lignification and secondary xylem thickness leading to better physiological response under drought stress. Taken together, we suggest that Bm-miRNA172c-5p might be a key player in maintaining the native phenotype of B. monnieri under control and different environmental conditions.

Epigenetic regulation of nitrogen and phosphorus responses in plants

Nitrogen (N) and phosphorus (P) are two of the most important nutrients for plant growth and crop yields. In the last decade, plenty of studies have revealed the genetic factors and their regulatory networks which are involved in N and/or P uptake and utilization in different model plant species, especially in Arabidopsis and rice. However, increasing evidences have shown that epigenetic regulation also plays a vital role in modulating plant responses to nutrient availability. In this review, we make a brief summary of epigenetic regulation including histone modifications, DNA methylation, and other chromatin structure alterations in tuning N and P responses. We also give an outlook for future research directions to comprehensively dissect the involvement of epigenetic regulation in modulating nutrient response in plants.

RNA-Based Technologies for Engineering Plant Virus Resistance

In recent years, non-coding RNAs (ncRNAs) have gained unprecedented attention as new and crucial players in the regulation of numerous cellular processes and disease responses. In this review, we describe how diverse ncRNAs, including both small RNAs and long ncRNAs, may be used to engineer resistance against plant viruses. We discuss how double-stranded RNAs and small RNAs, such as artificial microRNAs and trans-acting small interfering RNAs, either produced in transgenic plants or delivered exogenously to non-transgenic plants, may constitute powerful RNA interference (RNAi)-based technology that can be exploited to control plant viruses. Additionally, we describe how RNA guided CRISPR-CAS gene-editing systems have been deployed to inhibit plant virus infections, and we provide a comparative analysis of RNAi approaches and CRISPR-Cas technology. The two main strategies for engineering virus resistance are also discussed, including direct targeting of viral DNA or RNA, or inactivation of plant host susceptibility genes. We also elaborate on the challenges that need to be overcome before such technologies can be broadly exploited for crop protection against viruses.

Molecular Manipulation of the miR399/PHO2 Expression Module Alters the Salt Stress Response of Arabidopsis thaliana

In Arabidopsis thaliana (Arabidopsis), the microRNA399 (miR399)/PHOSPHATE2 (PHO2) expression module is central to the response of Arabidopsis to phosphate (PO₄) stress. In addition, miR399 has been demonstrated to also alter in abundance in response to salt stress. We therefore used a molecular modification approach to alter miR399 abundance to investigate the requirement of altered miR399 abundance in Arabidopsis in response to salt stress. The generated transformant lines, MIM399 and MIR399 plants, with reduced and elevated miR399 abundance respectively, displayed differences in their phenotypic and physiological response to those of wild-type Arabidopsis (Col-0) plants following exposure to a 7-day period of salt stress. However, at the molecular level, elevated miR399 abundance, and therefore, altered PHO2 target gene expression in salt-stressed Col-0, MIM399 and MIR399 plants, resulted in significant changes to the expression level of the two PO₄ transporter genes, PHOSPHATE TRANSPORTER1;4 (PHT1;4) and PHT1;9. Elevated PHT1;4 and PHT1;9 PO₄ transporter levels in salt stressed Arabidopsis would enhance PO₄ translocation from the root to the shoot tissue which would supply additional levels of this precious cellular resource that could be utilized by the aerial tissues of salt stressed Arabidopsis to either maintain essential biological processes or to mount an adaptive response to salt stress.

Whole-Transcriptome RNA Sequencing Reveals the Global Molecular Responses and CeRNA Regulatory Network of mRNAs, lncRNAs, miRNAs and circRNAs in Response to Salt Stress in Sugar Beet (Beta vulgaris)

Sugar beet is an important sugar-yielding crop with some tolerance to salt, but the mechanistic basis of this tolerance is not known. In the present study, we have used whole-transcriptome RNA-seq and degradome sequencing in response to salt stress to uncover differentially expressed (DE) mRNAs, microRNAs (miRNAs), long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) in both leaves and roots. A competitive endogenous RNA (ceRNA) network was constructed with the predicted DE pairs, which revealed regulatory roles under salt stress. A functional analysis suggests that ceRNAs are implicated in copper redistribution, plasma membrane permeability, glycometabolism and energy metabolism, NAC transcription factor and the phosphoinositol signaling system. Overall, we conducted for the first time a full transcriptomic analysis of sugar beet under salt stress that involves a potential ceRNA network, thus providing a basis to study the potential functions of lncRNAs/circRNAs.

PendoTMBase: A Database for Plant Endogenous Target Mimics

With fast-evolving next-generation sequencing technology, a great amount of plant genome and transcriptome data are becoming available. Due to the availability of mature microRNA (miRNA) sequence information from the miRBase (release 21) database, it is possible to predict endogenous target mimics (eTMs) in plant by searching seed-matched target sites. We identified a total of 2669 non-redundant eTM records in 43 plant species to create a specialized web-based database platform. The platform is named PendoTMBase, which can provide details of the eTMs. Predicted pairing structure between eTMs and their target miRNA, expression levels of eTMs pairs and associated GO annotations are also stored in the database. With evaluations performed in silico, we have found that the eTMs are representative and fit for genetic studies and high value-added plants breeding. The platform also provides a function of predicting novel eTMs based on the miRNA sequences submitted by users.

Global identification and analysis revealed differentially expressed lncRNAs associated with meiosis and low fertility in autotetraploid rice

BACKGROUND

Autotetraploid rice is a useful germplasm for polyploid rice breeding. Our previous research showed that non-coding RNAs might be associated with low fertility in autotetraploid rice. However, little information is available on long non-coding RNAs (lncRNAs) involved in the low fertility of autotetraploid rice. In the present study, RNA-seq was employed to detect the differentially expressed meiosis-related lncRNAs in autotetraploid rice, and gene overexpression and knock out experiments were used to validate the potential function of candidate lncRNA.

RESULTS

A total of 444 differentially expressed lncRNAs (DEL) were detected during anther and ovary meiosis in autotetraploid rice. Of these, 328 DEL were associated with the transposable elements, which displayed low expression levels during meiosis in autotetraploid rice. We used rapid amplification of cDNA ends (RACE) assay to validate 10 DEL and found that the lncRNAs were not assembly artifacts, and six of them were conserved in tetraploid rice. Moreover, 237 and 20 lncRNAs were associated with pollen mother cell (PMC) and embryo sac mother cell (EMC) meiosis in autotetraploid rice, respectively. The differential expressions of some meiosis-related targets and its DEL regulator, including MEL1 regulated by TCONS_00068868, LOC_Os12g41350 (meiotic asynaptic mutant 1) by TCONS_00057811 in PMC, and LOC_Os12g39420 by TCONS_00144592 in EMC, were confirmed by qRT-PCR. TCONS_00057811, TCONS_00055980 and TCONS_00130461 showed anther specific expression patterns and were found to be highly expressed during meiosis. CRISPR/Cas9 editing of lncRNA57811 displayed similar morphology compared to wild type. The overexpression of lncRNA57811 resulted in low pollen fertility (29.70%) and seed setting (33%) in rice.

CONCLUSION

The differential expression levels of lncRNAs, associated with transposable elements and meiosis-regulated targets, might be endogenous noncoding regulators of pollen/embryo sac development that cause low fertility in autotetraploid rice. The results enhance our understanding about rice lncRNAs, and facilitate functional research in autotetraploid rice.

Identification and characterization of circular RNAs during wood formation of poplars in acclimation to low nitrogen availability

Circular RNA (circRNA) identification and expression profiles, and construction of circRNAs-miRNAs-mRNAs networks indicates that circRNAs are involved in wood formation of poplars in acclimation to low nitrogen availability. Circular RNAs (circRNAs) are covalently closed non-coding RNAs that play pivotal roles in various biological processes. However, circRNAs' roles in wood formation of poplars in acclimation to low nitrogen (N) availability are currently unknown. Here, we undertook a systematic identification and characterization of circRNAs in the wood of Populus × canescens exposed to either 50 (low N) or 500 (normal N) µM NH₄NO₃ using rRNA-depleted RNA-sequencing. A total of 2,509 unique circRNAs were identified, and 163 (ca. 6.5%) circRNAs were significantly differentially expressed (DE) under low N condition. We observed a positive correlation between the expression patterns of DE circRNAs and their hosting protein-coding genes. Moreover, circRNAs-miRNAs-mRNAs' networks were identified in the wood of poplars under low N availability. For instance, upregulated several circRNAs, such as circRNA1226, circRNA 1732, and circRNA392 induced increases in nuclear factor Y, subunit A1-A (NFYA1-A), NFYA1-B, and NFYA10 transcript levels via the mediation of miR169b members, which is in line with reduced xylem width and cell layers of the xylem in the wood of low N-supplied poplars. Upregulation of circRNA1006, circRNA1344, circRNA1941, circRNA901, and circRNA146 caused increased transcript level of MYB61 via the mediation of a miR5021 member, corresponding well to the higher lignin concentration in the wood of low N-treated poplars. Overall, these results indicated that DE circRNAs play an essential role in regulating gene expression via circRNAs-miRNAs-mRNAs' networks to modulate wood anatomical and chemical properties of poplars in acclimation to low N availability.

Novel Aspects of Nitrate Regulation in Arabidopsis

Nitrogen (N) is one of the most essential macronutrients for plant growth and development. Nitrate (NO₃ ^-), the major form of N that plants uptake from the soil, acts as an important signaling molecule in addition to its nutritional function. Over the past decade, significant progress has been made in identifying new components involved in NO₃ ^- regulation and starting to unravel the NO₃ ^- regulatory network. Great reviews have been made recently by scientists on the key regulators in NO₃ ^- signaling, NO₃ ^- effects on plant development, and its crosstalk with phosphorus (P), potassium (K), hormones, and calcium signaling. However, several novel aspects of NO₃ ^- regulation have not been previously reviewed in detail. Here, we mainly focused on the recent advances of post-transcriptional regulation and non-coding RNA (ncRNAs) in NO₃ ^- signaling, and NO₃ ^- regulation on leaf senescence and the circadian clock. It will help us to extend the general picture of NO₃ ^- regulation and provide a basis for further exploration of NO₃ ^- regulatory network.

Transcriptome and Phytochemical Analyses Provide New Insights Into Long Non-Coding RNAs Modulating Characteristic Secondary Metabolites of Oolong Tea (Camellia sinensis) in Solar-Withering

Oolong tea is a popular and semi-fermented beverage. During the processing of tea leaves, withering is the first indispensable process for improving flavor. However, the roles of long non-coding RNAs (lncRNAs) and the characteristic secondary metabolites during the withering of oolong tea leaves remain unknown. In this study, phytochemical analyses indicated that total polyphenols, flavonoids, catechins, epigallocatechin (EGC), catechin gallate (CG), gallocatechin gallate (GCG), epicatechin gallate (ECG), and epigallocatechin gallate (EGCG) were all less abundant in the solar-withered leaves (SW) than in the fresh leaves (FL) and indoor-withered leaves (IW). In contrast, terpenoid, jasmonic acid (JA), and methyl jasmonate (MeJA) contents were higher in the SW than in the FL and IW. By analyzing the transcriptome data, we detected 32,036 lncRNAs. On the basis of the Kyoto Encyclopedia of Genes and Genomes analysis, the flavonoid metabolic pathway, the terpenoid metabolic pathway, and the JA/MeJA biosynthesis and signal transduction pathway were enriched pathways. Additionally, 63 differentially expressed lncRNAs (DE-lncRNAs) and 23 target genes were identified related to the three pathways. A comparison of the expression profiles of the DE-lncRNAs and their target genes between the SW and IW revealed four up-regulated genes (FLS, CCR, CAD, and HCT), seven up-regulated lncRNAs, four down-regulated genes (4CL, CHI, F3H, and F3'H), and three down-regulated lncRNAs related to flavonoid metabolism; nine up-regulated genes (DXS, CMK, HDS, HDR, AACT, MVK, PMK, GGPPS, and TPS), three up-regulated lncRNAs, and six down-regulated lncRNAs related to terpenoid metabolism; as well as six up-regulated genes (LOX, AOS, AOC, OPR, ACX, and MFP2), four up-regulated lncRNAs, and three down-regulated lncRNAs related to JA/MeJA biosynthesis and signal transduction. These results suggested that the expression of DE-lncRNAs and their targets involved in the three pathways may be related to the low abundance of the total polyphenols, flavonoids, and catechins (EGC, CG, GCG, ECG, and EGCG) and the high abundance of terpenoids in the SW. Moreover, solar irradiation, high JA and MeJA contents, and the endogenous target mimic (eTM)-related regulatory mechanism in the SW were also crucial for increasing the terpenoid levels. These findings provide new insights into the greater contribution of solar-withering to the high-quality flavor of oolong tea compared with the effects of indoor-withering.

Competing Endogenous RNA Networks Underlying Anatomical and Physiological Characteristics of Poplar Wood in Acclimation to Low Nitrogen Availability

Although poplar plantations are often established on nitrogen (N)-poor soil, the physiological and molecular mechanisms underlying wood properties of poplars in acclimation to low N availability remain largely unknown. To investigate wood properties of poplars in acclimation to low N, Populus � canescens saplings were exposed to either 50 (low N) or 500 (normal N) �M NH4NO3 for 2 months. Low N resulted in decreased xylem width and cell layers of the xylem (the number of cells counted along the ray parenchyma on the stem cross section), narrower lumina of vessels and fibers, greater thickness of double fiber walls (the walls between two adjacent fiber cells), more hemicellulose and lignin deposition, and reduced cellulose accumulation in poplar wood. Consistently, concentrations of gibberellins involved in cell size determination and the abundance of various metabolites including amino acids, carbohydrates and precursors for cell wall biosynthesis were decreased in low N-supplied wood. In line with these anatomical and physiological changes, a number of mRNAs, long noncoding RNAs (lncRNAs) and microRNAs (miRNAs) were significantly differentially expressed. Competing endogenous RNA regulatory networks were identified in the wood of low N-treated poplars. Overall, these results indicate that miRNAs-lncRNAs-mRNAs networks are involved in regulating wood properties and physiological processes of poplars in acclimation to low N availability.

The long noncoding RNA T5120 regulates nitrate response and assimilation in Arabidopsis

Long noncoding RNAs (lncRNAs) are crucial regulators in many plant biological processes. However, it remains unknown whether lncRNAs can respond to nitrate or function in nitrate regulation. We detected 695 lncRNAs, 480 known and 215 novel, in Arabidopsis seedling roots; six showed altered expression in response to nitrate treatment, among which T5120 showed the highest induction. Overexpression of T5120 in Arabidopsis promoted the response to nitrate, enhanced nitrate assimilation and improved biomass and root development. Biochemical and molecular analyses revealed that NLP7, a master nitrate regulatory transcription factor, directly bound to the nitrate-responsive cis-element (NRE)-like motif of the T5120 promoter and activated T5120 transcription. In addition, T5120 partially restored the nitrate signalling and assimilation phenotypes of nlp7 mutant, suggesting that T5120 is involved in NLP7-mediated nitrate regulation. Interestingly, the expression of T5120 was regulated by the nitrate sensor NRT1.1. Therefore, T5120 is modulated by NLP7 and NRT1.1 to regulate nitrate signalling. Our work reveals a new regulatory mechanism in which lncRNA T5120 functions in nitrate regulation, providing new insights into the nitrate signalling network. Importantly, lncRNA T5120 can promote nitrate assimilation and plant growth to improve nitrogen use efficiency.

Identification of lncRNAs Responsive to Infection by Plasmodiophora brassicae in Clubroot-Susceptible and -Resistant Brassica napus Lines Carrying Resistance Introgressed from Rutabaga

Clubroot disease, caused by Plasmodiophora brassicae Woronin, is a major threat to the production of Brassica' crops. Resistance to different P. brassicae pathotypes has been reported in the A genome, chromosome A08; however, the molecular mechanism of this resistance, especially the involvement of long noncoding RNAs (lncRNAs), is not understood. We have used a strand-specific lncRNA-Seq approach to catalog lncRNAs from the roots of clubroot-susceptible and -resistant Brassica napus lines. In total, 530 differentially expressed (DE) lncRNAs were identified, including 88% of long intergenic RNAs and 11% natural antisense transcripts. Sixteen lncRNAs were identified as target mimics of the microRNAs (miRNAs) and eight were identified as the precursors of miRNAs. KEGG pathway analysis of the DE lncRNAs showed that the cis-regulated target genes mostly belong to the phenylpropanoid biosynthetic pathway (15%) and plant-pathogen interactions (15%) while the transregulated target genes mostly belong to carbon (18%) and amino acid biosynthesis pathway (19%). In all, 24 DE lncRNAs were identified from chromosome A08, which is known to harbor a quantitative trait locus conferring resistance to different P. brassicae pathotypes; however, eight of these lncRNAs showed expression only in the resistant plants. These results could form the basis for future studies aimed at delineating the roles of lncRNAs in plant-microbe interactions.

Long non-coding RNAs: Fine-tuning the developmental responses in plants

Plant developmental biology is associated with various gene regulatory pathways involved in different phases of their life cycle. In course of development, growth and differentiation of different organs in plants are regulated by specific sets of gene expression. With the advances in genomic and bioinformatic techniques, particularly high-throughput sequencing technology, many transcriptional units with no protein-coding potential have been discovered. Previously thought to be the dark matters of genome, long non-coding RNAs (lncRNAs) are gradually gaining importance as crucial players in gene regulation during different developmental phases. Some lncRNAs, showing complementarity to microRNAs (miRNAs), are used as endogenous target mimics of specific miRNA family. A number of lncRNAs can also act as natural antisense transcripts to attenuate the expression of coding genes. Although lncRNA-mediated regulations have extensively been studied in animals, plant lncRNA research is still in its initial phase. The present review highlights the regulatory mechanism and different physiological aspects of lncRNAs in plant development. In plants, lncRNAs are found to be associated with a number of plant developmental functions such as lateral root development, vernalization, photomorphogenesis, pollen development, fiber development and nodulation. Understanding these potent roles of lncRNAs in plant development can further provide novel tools for crop improvement programs in future.

ceRNAs in plants: computational approaches and associated challenges for target mimic research

The competing endogenous RNA hypothesis has gained increasing attention as a potential global regulatory mechanism of microRNAs (miRNAs), and as a powerful tool to predict the function of many noncoding RNAs, including miRNAs themselves. Most studies have been focused on animals, although target mimic (TMs) discovery as well as important computational and experimental advances has been developed in plants over the past decade. Thus, our contribution summarizes recent progresses in computational approaches for research of miRNA:TM interactions. We divided this article in three main contributions. First, a general overview of research on TMs in plants is presented with practical descriptions of the available literature, tools, data, databases and computational reports. Second, we describe a common protocol for the computational and experimental analyses of TM. Third, we provide a bioinformatics approach for the prediction of TM motifs potentially cross-targeting both members within the same or from different miRNA families, based on the identification of consensus miRNA-binding sites from known TMs across sequenced genomes, transcriptomes and known miRNAs. This computational approach is promising because, in contrast to animals, miRNA families in plants are large with identical or similar members, several of which are also highly conserved. From the three consensus TM motifs found with our approach: MIM166, MIM171 and MIM159/319, the last one has found strong support on the recent experimental work by Reichel and Millar [Specificity of plant microRNA TMs: cross-targeting of mir159 and mir319. J Plant Physiol 2015;180:45-8]. Finally, we stress the discussion on the major computational and associated experimental challenges that have to be faced in future ceRNA studies.

Functional Identification of Corynespora cassiicola-Responsive miRNAs and Their Targets in Cucumber

Target leaf spot (TLS), which is caused by Corynespora cassiicola (C. cassiicola), is one of the most important diseases in cucumber (Cucumis sativus L.). Our previous research identified several C. cassiicola-responsive miRNAs in cucumber by high-throughput sequencing, including two known miRNAs and two novel miRNAs. The target genes of these miRNAs were related to secondary metabolism. In this study, we verified the interaction between these miRNAs and target genes by histochemical staining and fluorescence quantitative assays of GUS. We transiently expressed the candidate miRNAs and target genes in cucumber cotyledons to investigate the resistance to C. cassiicola. Transient expression of miR164d, miR396b, Novel-miR1, and Novel-miR7 in cucumber resulted in decreased resistance to C. cassiicola, while transient expression of NAC (inhibited by miR164d), APE (inhibited by miR396b), 4CL (inhibited by Novel-miR1), and PAL (inhibited by Novel-miR7) led to enhanced resistance to C. cassiicola. In addition, overexpression of 4CL and PAL downregulated lignin synthesis, and overexpression of Novel-miR1 and Novel-miR7 also downregulated lignin synthesis, indicating that the regulation of 4CL and PAL by Novel-miR1 and Novel-miR7 could affect lignin content. The tobacco rattle virus (TRV) induced short tandem target mimic (STTM)-miRNA silencing vector was successfully constructed, and target miRNAs were successfully silenced. The identification of disease resistance and lignin content showed that silencing candidate miRNAs could improve cucumber resistance to C. cassiicola.

Abiotic stress induced miRNA-TF-gene regulatory network: A structural perspective

MicroRNAs (miRNAs) and transcription factors (TFs) are the largest families of trans-acting gene regulatory species, which are pivotal players in a complex regulatory network. Recently, extensive research on miRNAs and TFs in agriculture has identified these trans-acting regulatory species, as an effective tool for engineering new crop cultivars to increase yield and quality as well tolerance to environmental stresses but our knowledge of regulatory network is still not sufficient to decipher the exact mechanism. In the current work, stress-specific TF-miRNA-gene network was built for Arabidopsis under drought, cold, salt and waterlogging stress using data from reliable publically available databases; and transcriptome and degradome sequence data analysis by meta-analysis approach. Further network analysis elucidated significantly dense, scale-free, small world and hierarchical backbone of interactions. The various centrality measures highlighted several genes/TF/miRNAs as potential targets for tolerant variety cultivation. This comprehensive regulatory information will accelerate the advancement of current understanding on stress specific transcriptional and post-transcriptional regulatory mechanism and has promising utilizations for experimental biologist who are intended to improve plant crop performance under multiple Abiotic stress environments.

Phyto-miRNAs-based regulation of metabolites biosynthesis in medicinal plants

Medicinal plants, are known to produce a wide range of plant secondary metabolites (PSMs) applied as insecticides, drugs, dyes and toxins in agriculture, medicine, industry and bio-warfare plus bio-terrorism, respectively. However, production of PSMs is usually in small quantities, so we need to find novel ways to increase both quantity and quality of them. Fortunately, biotechnology suggests several options through which secondary metabolism in plants can be engineered in innovative ways to: 1) over-produce the useful metabolites, 2) down-produce the toxic metabolites, 3) produce the new metabolites. Among the ways, RNA interference (RNAi) technology which involves gene-specific regulation by small non-coding RNAs (sncRNAs) have been recently emerged as a promising tool for plant biotechnologist, not only to decipher the function of plant genes, but also for development of the plants with improved and novel traits through manipulation of both desirable and undesirable genes. Among sncRNAs, miRNAs have been recorded various regulatory roles in plants such as development, signal transduction, response to environmental stresses, metabolism. Certainly, the use of miRNAs in metabolic engineering requires identification of miRNAs involved in metabolites biosynthesis, understanding of the biosynthetic pathways, as well as the identification of key points of the pathways in which the miRNAs have their own effect. Thus, we firstly consider these three issues on metabolic engineering of medicinal plants. Our review shows, application of miRNAs can open a novel perspective to metabolic engineering of medicinal plants.

MicroRNA MIMIC binding sites: Minor flanking nucleotide alterations can strongly impact MIMIC silencing efficacy in Arabidopsis

In plants, microRNA (miRNA) target MIMICs (MIMs) have been widely used to inhibit miRNA function. They are based on the Arabidopsis INSENSITIVE TO PHOSPATE STARVATION 1 (IPS1) gene that corresponds to a non-coding RNA containing a miR399 binding site that can be modified to sequester and inhibit any miRNA of interest. However, the efficacy of miRNA inhibition of these different MIMs can vary greatly. Using MIMs that have strong efficacy (MIM159) and poor efficacy (MIM165), we investigate the underlying cause of this variation. Firstly, sequence alignments of IPS1 homologs from the Brassicaceae identified a highly conserved sequence immediately downstream of the miRNA binding site. Mutating this sequence in the context of the MIM159 attenuates its strong efficacy. This conserved flanking region contains a predicted stem-loop structure that is also predicted to be present in most modified MIMs that appear to have a strong efficacy, but not in MIM165 that has a poor efficacy. Restoring this predicted stem-loop in MIM165 via mutation of only three or five nucleotides within the conserved flanking region resulted in MIM165 variants that have very strong efficacies of miRNA inhibition. However, specifically mutating this predicted stem-loop in the MIM159 context failed to significantly reduce efficacy, and additional mutations to restore this predicted stem-loop weakened efficacy further. Although this shows there is no simple correlation between this predicted stem-loop and efficacy, these results add to the growing evidence that the sequence context of miRNA binding sites is important, and that minor nucleotide substitutions to flanking sequences of miRNA binding sites can strongly enhance or attenuate the miRNA-target interaction.

Identification of Gossypium hirsutum long non-coding RNAs (lncRNAs) under salt stress

BACKGROUND

Long non-coding RNAs (lncRNAs) represent a class of riboregulators that either directly act in long form or are processed into shorter microRNAs (miRNAs) and small interfering RNAs. Long noncoding RNAs (lncRNAs) are arbitrarily defined as RNA genes larger than 200 nt in length that have no apparent coding potential. lncRNAs have emerged as playing important roles in various biological regulatory processes and are expressed in a more tissue-specific manner than mRNA. Emerging evidence shows that lncRNAs participate in stress-responsive regulation.

RESULTS

In this study, in order to develop a comprehensive catalogue of lncRNAs in upland cotton under salt stress, we performed whole-transcriptome strand-specific RNA sequencing for three-leaf stage cotton seedlings treated with salt stress (S_NaCl) and controls (S_CK). In total we identified 1117 unique lncRNAs in this study and 44 differentially expressed RNAs were identified as potential non-coding RNAs. For the differentially expressed lncRNAs that were identified as intergenic lncRNAs (lincRNA), we analysed the gene ontology enrichment of cis targets and found that cis target protein-coding genes were mainly enriched in stress-related categories. Real-time quantitative PCR confirmed that all selected lincRNAs responsive to salt stress. We found lnc_388 was likely as regulator of Gh_A09G1182. And lnc_883 may participate in regulating tolerance to salt stress by modulating the expression of Gh_D03G0339 MS_channel. We then predicted the target mimics for miRNA in Gossypium. six miRNAs were identified, and the result of RT-qPCR with lncRNA and miRNA suggested that lnc_973 and lnc_253 may regulate the expression of ghr-miR399 and ghr-156e as a target mimic under salt stress.

CONCLUSIONS

We identified 44 lincRNAs that were differentially expressed under salt stress. These lincRNAs may target protein-coding genes via cis-acting regulation. We also discovered that specifically-expressed lincRNAs under salt stress may act as endogenous target mimics for conserved miRNAs. These findings extend the current view on lincRNAs as ubiquitous regulators under stress stress.

Long noncoding RNAs in the model species Brachypodium distachyon

Eukaryotic genomes are pervasively transcribed and only a small portion of the transcribed sequences belongs to protein coding genes. High-throughput sequencing technology contributed to consolidate this perspective, allowing the identification of numerous noncoding RNAs with key roles in biological processes. Long noncoding RNAs (lncRNAs) are transcripts longer than 200 nt with limited phylogenetic conservation, expressed at low levels and characterized by tissue/organ specific expression profiles. Although a large set of lncRNAs has been identified, the functional roles of lncRNAs are only beginning to be recognized and the molecular mechanism of lncRNA-mediated gene regulation remains largely unexplored, particularly in plants where their annotation and characterization are still incomplete. Using public and proprietary poly-(A)+ RNA-seq data as well as a collection of full length ESTs from several organs, developmental stages and stress conditions in three Brachypodium distachyon inbred lines, we describe the identification and the main features of thousands lncRNAs. Here we provide a genome-wide characterization of lncRNAs, highlighting their intraspecies conservation and describing their expression patterns among several organs/tissues and stress conditions. This work represents a fundamental resource to deepen our knowledge on long noncoding RNAs in C3 cereals, allowing the Brachypodium community to exploit these results in future research programs.

Regulation of rice root development by a retrotransposon acting as a microRNA sponge

It is well documented that transposable elements (TEs) can regulate the expression of neighbouring genes. However, their ability to act in trans and influence ectopic loci has been reported rarely. We searched in rice transcriptomes for tissue-specific expression of TEs and found them to be regulated developmentally. They often shared sequence homology with co-expressed genes and contained potential microRNA-binding sites, which suggested possible contributions to gene regulation. In fact, we have identified a retrotransposon that is highly transcribed in roots and whose spliced transcript constitutes a target mimic for miR171. miR171 destabilizes mRNAs encoding the root-specific family of SCARECROW-Like transcription factors. We demonstrate that retrotransposon-derived transcripts act as decoys for miR171, triggering its degradation and thus results in the root-specific accumulation of SCARECROW-Like mRNAs. Such transposon-mediated post-transcriptional control of miR171 levels is conserved in diverse rice species.

Regulation of rice root development by a retrotransposon acting as a microRNA sponge

It is well documented that transposable elements (TEs) can regulate the expression of neighbouring genes. However, their ability to act in trans and influence ectopic loci has been reported rarely. We searched in rice transcriptomes for tissue-specific expression of TEs and found them to be regulated developmentally. They often shared sequence homology with co-expressed genes and contained potential microRNA-binding sites, which suggested possible contributions to gene regulation. In fact, we have identified a retrotransposon that is highly transcribed in roots and whose spliced transcript constitutes a target mimic for miR171. miR171 destabilizes mRNAs encoding the root-specific family of SCARECROW-Like transcription factors. We demonstrate that retrotransposon-derived transcripts act as decoys for miR171, triggering its degradation and thus results in the root-specific accumulation of SCARECROW-Like mRNAs. Such transposon-mediated post-transcriptional control of miR171 levels is conserved in diverse rice species.

Improving the Quality of Positive Datasets for the Establishment of Machine Learning Models for pre-microRNA Detection

MicroRNAs (miRNAs) are involved in the post-transcriptional regulation of protein abundance and thus have a great impact on the resulting phenotype. It is, therefore, no wonder that they have been implicated in many diseases ranging from virus infections to cancer. This impact on the phenotype leads to a great interest in establishing the miRNAs of an organism. Experimental methods are complicated which led to the development of computational methods for pre-miRNA detection. Such methods generally employ machine learning to establish models for the discrimination between miRNAs and other sequences. Positive training data for model establishment, for the most part, stems from miRBase, the miRNA registry. The quality of the entries in miRBase has been questioned, though. This unknown quality led to the development of filtering strategies in attempts to produce high quality positive datasets which can lead to a scarcity of positive data. To analyze the quality of filtered data we developed a machine learning model and found it is well able to establish data quality based on intrinsic measures. Additionally, we analyzed which features describing pre-miRNAs could discriminate between low and high quality data. Both models are applicable to data from miRBase and can be used for establishing high quality positive data. This will facilitate the development of better miRNA detection tools which will make the prediction of miRNAs in disease states more accurate. Finally, we applied both models to all miRBase data and provide the list of high quality hairpins.

Genome-wide Discovery of Circular RNAs in the Leaf and Seedling Tissues of Arabidopsis Thaliana

BACKGROUND

Recently, identification and functional studies of circular RNAs, a type of non-coding RNAs arising from a ligation of 3' and 5' ends of a linear RNA molecule, were conducted in mammalian cells with the development of RNA-seq technology.

METHOD

Since compared with animals, studies on circular RNAs in plants are less thorough, a genome-wide identification of circular RNA candidates in Arabidopsis was conducted with our own developed bioinformatics tool to several existing RNA-seq datasets specifically for non-coding RNAs.

RESULTS

A total of 164 circular RNA candidates were identified from RNA-seq data, and 4 circular RNA transcripts, including both exonic and intronic circular RNAs, were experimentally validated. Interestingly, our results show that circular RNA transcripts are enriched in the photosynthesis system for the leaf tissue and correlated to the higher expression levels of their parent genes. Sixteen out of all 40 genes that have circular RNA candidates are related to the photosynthesis system, and out of the total 146 exonic circular RNA candidates, 63 are found in chloroplast.

Discovery and Annotation of Plant Endogenous Target Mimicry Sequences from Public Transcriptome Libraries: A Case Study of Prunus persica

Novel transcript discovery through RNA sequencing has substantially improved our understanding of the transcriptome dynamics of biological systems. Endogenous target mimicry (eTM) transcripts, a novel class of regulatory molecules, bind to their target microRNAs (miRNAs) by base pairing and block their biological activity. The objective of this study was to provide a computational analysis framework for the prediction of putative eTM sequences in plants, and as an example, to discover previously un-annotated eTMs in Prunus persica (peach) transcriptome. Therefore, two public peach transcriptome libraries downloaded from Sequence Read Archive (SRA) and a previously published set of long non-coding RNAs (lncRNAs) were investigated with multi-step analysis pipeline, and 44 putative eTMs were found. Additionally, an eTM-miRNA-mRNA regulatory network module associated with peach fruit organ development was built via integration of the miRNA target information and predicted eTM-miRNA interactions. My findings suggest that one of the most widely expressed miRNA families among diverse plant species, miR156, might be potentially sponged by seven putative eTMs. Besides, the study indicates eTMs potentially play roles in the regulation of development processes in peach fruit via targeting specific miRNAs. In conclusion, by following the step-by step instructions provided in this study, novel eTMs can be identified and annotated effectively in public plant transcriptome libraries.

Genome-wide identification and functional prediction of cold and/or drought-responsive lncRNAs in cassava

Cold and drought stresses seriously affect cassava (Manihot esculenta) plant growth and yield. Recently, long noncoding RNAs (lncRNAs) have emerged as key regulators of diverse cellular processes in mammals and plants. To date, no systematic screening of lncRNAs under abiotic stress and their regulatory roles in cassava has been reported. In this study, we present the first reference catalog of 682 high-confidence lncRNAs based on analysis of strand-specific RNA-seq data from cassava shoot apices and young leaves under cold, drought stress and control conditions. Among them, 16 lncRNAs were identified as putative target mimics of cassava known miRNAs. Additionally, by comparing with small RNA-seq data, we found 42 lncNATs and sense gene pairs can generate nat-siRNAs. We identified 318 lncRNAs responsive to cold and/or drought stress, which were typically co-expressed concordantly or discordantly with their neighboring genes. Trans-regulatory network analysis suggested that many lncRNAs were associated with hormone signal transduction, secondary metabolites biosynthesis, and sucrose metabolism pathway. The study provides an opportunity for future computational and experimental studies to uncover the functions of lncRNAs in cassava.

Root transcriptome of two contrasting indica rice cultivars uncovers regulators of root development and physiological responses

The huge variation in root system architecture (RSA) among different rice (Oryza sativa) cultivars is conferred by their genetic makeup and different growth or climatic conditions. Unlike model plant Arabidopsis, the molecular basis of such variation in RSA is very poorly understood in rice. Cultivars with stable variation are valuable resources for identification of genes involved in RSA and related physiological traits. We have screened for RSA and identified two such indica rice cultivars, IR-64 (OsAS83) and IET-16348 (OsAS84), with stable contrasting RSA. OsAS84 produces robust RSA with more crown roots, lateral roots and root hairs than OsAS83. Using comparative root transcriptome analysis of these cultivars, we identified genes related to root development and different physiological responses like abiotic stress responses, hormone signaling, and nutrient acquisition or transport. The two cultivars differ in their response to salinity/dehydration stresses, phosphate/nitrogen deficiency, and different phytohormones. Differential expression of genes involved in salinity or dehydration response, nitrogen (N) transport, phosphate (Pi) starvation signaling, hormone signaling and root development underlies more resistance of OsAS84 towards abiotic stresses, Pi or N deficiency and its robust RSA. Thus our study uncovers gene-network involved in root development and abiotic stress responses in rice.

MicroRNA858 Is a Potential Regulator of Phenylpropanoid Pathway and Plant Development

MicroRNAs (miRNAs) are endogenous, noncoding small RNAs that function as critical regulators of gene expression. In plants, miRNAs have shown their potential as regulators of growth, development, signal transduction, and stress tolerance. Although the miRNA-mediated regulation of several processes is known, the involvement of miRNAs in regulating secondary plant product biosynthesis is poorly understood. In this study, we functionally characterized Arabidopsis (Arabidopsis thaliana) miR858a, which putatively targets R2R3-MYB transcription factors involved in flavonoid biosynthesis. Overexpression of miR858a in Arabidopsis led to the down-regulation of several MYB transcription factors regulating flavonoid biosynthesis. In contrast to the robust growth and early flowering of miR858OX plants, reduction of plant growth and delayed flowering were observed in Arabidopsis transgenic lines expressing an artificial miRNA target mimic (MIM858). Genome-wide expression analysis using transgenic lines suggested that miR858a targets a number of regulatory factors that modulate the expression of downstream genes involved in plant development and hormonal and stress responses. Furthermore, higher expression of MYBs in MIM858 lines leads to redirection of the metabolic flux towards the synthesis of flavonoids at the cost of lignin synthesis. Altogether, our study has established the potential role of light-regulated miR858a in flavonoid biosynthesis and plant growth and development.

The use of high-throughput sequencing methods for plant microRNA research

MicroRNA (miRNA) acts as a critical regulator of gene expression at post-transcriptional and occasionally transcriptional levels in plants. Identification of reliable miRNA genes, monitoring the procedures of transcription, processing and maturation of the miRNAs, quantification of the accumulation levels of the miRNAs in specific biological samples, and validation of miRNA-target interactions become the basis for thoroughly understanding of the miRNA-mediated regulatory networks and the underlying mechanisms. Great progresses have been achieved for sequencing technology. Based on the high degree of sequencing depth and coverage, the high-throughput sequencing (HTS, also called next-generation sequencing) technology provides unprecedentedly efficient way for genome-wide or transcriptome-wide studies. In this review, we will introduce several HTS platform-based methods useful for plant miRNA research, including RNA-seq (RNA sequencing), RNA-PET-seq (paired end tag sequencing of RNAs), sRNA-seq (small RNA sequencing), dsRNA-seq (double-stranded RNA sequencing), ssRNA-seq (single-stranded RNA sequencing) and degradome-seq (degradome sequencing). In particular, we will provide some special cases to illustrate the novel use of HTS methods for investigation of the processing modes of the miRNA precursors, identification of the RNA editing sites on miRNA precursors, mature miRNAs and target transcripts, re-examination of the current miRNA registries, and discovery of novel miRNA species and novel miRNA-target interactions. Summarily, we opinioned that integrative use of the above mentioned HTS methods could make the studies on miRNAs more efficient.

Genome-wide identification and functional analysis of lincRNAs acting as miRNA targets or decoys in maize

Background

Long intergenic noncoding RNAs (lincRNAs) are endogenous non-coding RNAs (ncRNAs) that are transcribed from ‘intergenic’ regions of the genome and may play critical roles in regulating gene expression through multiple RNA-mediated mechanisms. MicroRNAs (miRNAs) are single-stranded small ncRNAs of approximately 21–24 nucleotide (nt) that are involved in transcriptional and post-transcriptional gene regulation. While miRNAs functioning as mRNA repressors have been studied in detail, the influence of miRNAs on lincRNAs has seldom been investigated in plants.

Methods

LincRNAs as miRNA targets or decoys were predicted via GSTAr.pl script with a set of rules, and lincRNAs as miRNA targets were validated by degradome data. Conservation analysis of lincRNAs as miRNA targets or decoys were conducted using BLASTN and MAFFT. The function of lincRNAs as miRNA targets were predicted via a lincRNA-mRNA co-expression network, and the function of lincRNAs as miRNA decoys were predicted according to the competing endogenous RNA (ceRNA) hypothesis.

Results

In this work, we developed a computational method and systematically predicted 466 lincRNAs as 165 miRNA targets and 86 lincRNAs as 58 miRNA decoys in maize (Zea mays L.). Furthermore, 34 lincRNAs predicted as 33 miRNA targets were validated based on degradome data. We found that lincRNAs acting as miRNA targets or decoys are a common phenomenon, which indicates that the regulated networks of miRNAs also involve lincRNAs. To elucidate the function of lincRNAs, we reconstructed a miRNA-regulated network involving 78 miRNAs, 117 lincRNAs and 8834 mRNAs. Based on the lincRNA-mRNA co-expression network and the competing endogenous RNA hypothesis, we predicted that 34 lincRNAs that function as miRNA targets and 86 lincRNAs that function as miRNA decoys participate in cellular and metabolic processes, and play role in catalytic activity and molecular binding functions.

Conclusions

This work provides a comprehensive view of miRNA-regulated networks and indicates that lincRNAs can participate in a layer of regulatory interactions as miRNA targets or decoys in plants, which will enable in-depth functional analysis of lincRNAs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2024-0) contains supplementary material, which is available to authorized users.

Regulation of Nicotine Biosynthesis by an Endogenous Target Mimicry of MicroRNA in Tobacco

The interaction between noncoding endogenous target mimicry (eTM) and its corresponding microRNA (miRNA) is a newly discovered regulatory mechanism and plays pivotal roles in various biological processes in plants. Tobacco (Nicotiana tabacum) is a model plant for studying secondary metabolite alkaloids, of which nicotine accounts for approximately 90%. In this work, we identified four unique tobacco-specific miRNAs that were predicted to target key genes of the nicotine biosynthesis and catabolism pathways and an eTM, novel tobacco miRNA (nta)-eTMX27, for nta-miRX27 that targets QUINOLINATE PHOSPHORIBOSYLTRANSFERASE2 (QPT2) encoding a quinolinate phosphoribosyltransferase. The expression level of nta-miRX27 was significantly down-regulated, while that of QPT2 and nta-eTMX27 was significantly up-regulated after topping, and consequently, nicotine content increased in the topping-treated plants. The topping-induced down-regulation of nta-miRX27 and up-regulation of QPT2 were only observed in plants with a functional nta-eTMX27 but not in transgenic plants containing an RNA interference construct targeting nta-eTMX27. Our results demonstrated that enhanced nicotine biosynthesis in the topping-treated tobacco plants is achieved by nta-eTMX27-mediated inhibition of the expression and functions of nta-miRX27. To our knowledge, this is the first report about regulation of secondary metabolite biosynthesis by an miRNA-eTM regulatory module in plants.

Long non-coding RNAs and their biological roles in plants

With the development of genomics and bioinformatics, especially the extensive applications of high-throughput sequencing technology, more transcriptional units with little or no protein-coding potential have been discovered. Such RNA molecules are called non-protein-coding RNAs (npcRNAs or ncRNAs). Among them, long npcRNAs or ncRNAs (lnpcRNAs or lncRNAs) represent diverse classes of transcripts longer than 200 nucleotides. In recent years, the lncRNAs have been considered as important regulators in many essential biological processes. In plants, although a large number of lncRNA transcripts have been predicted and identified in few species, our current knowledge of their biological functions is still limited. Here, we have summarized recent studies on their identification, characteristics, classification, bioinformatics, resources, and current exploration of their biological functions in plants.

Long noncoding RNA transcriptome of plants

Since their discovery more than two decades ago, animal long noncoding RNAs (lncRNAs) have emerged as important regulators of many biological processes. Recently, a large number of lncRNAs have also been identified in higher plants, and here, we review their identification, classification and known regulatory functions in various developmental events and stress responses. Knowledge gained from a deeper understanding of this special group of noncoding RNAs may lead to biotechnological improvement of crops. Some possible examples in this direction are discussed.

Motif content comparison between monocot and dicot species

While a number of DNA sequence motifs have been functionally characterized, the full repertoire of motifs in an organism (the motifome) is yet to be characterized. The present study wishes to widen the scope of motif content analysis in different monocot and dicot species that include both rice species, Brachypodium, corn, wheat as monocots and Arabidopsis, Lotus japonica, Medicago truncatula, and Populus tremula as dicots. All possible existing motifs were analyzed in different regions of genomes such as were found in different sets of sequences in these species: the whole genome, core proximal and distal promoters, 5′ and 3′ UTRs, and the 1st introns. Due to the increased number of species involved in this study compared to previous works, species relationships were analyzed based on the similarity of common motif content. Certain secondary structure elements were inferred in the genomes of these species as well as new unknown motifs. The distribution of 20 motifs common to the studied species were found to have a significantly larger occurrence within the promoters and 3′ UTRs of genes, both being regulatory regions. Motifs common to the promoter regions of japonica rice, Brachypodium, and corn were also found in a number of orthologous and paralogous genes. Some of our motifs were found to be complementary to miRNA elements in Brachypodium distachyon and japonica rice.

Genome-wide identification and functional prediction of novel and drought-responsive lincRNAs in Populus trichocarpa

Protein-coding genes are considered to be a dominant component of the eukaryotic transcriptome; however, many studies have shown that intergenic, non-coding transcripts also play an important role. Long intergenic non-coding RNAs (lincRNAs) were found to play a vital role in human and Arabidopsis. However, lincRNAs and their regulatory roles remain poorly characterized in woody plants, especially Populus trichocarpa (P. trichocarpa). A large set of Populus RNA-Seq data were examined with high sequencing depth under control and drought conditions and a total of 2542 lincRNA candidates were identified. In total, 51 lincRNAs and 20 lincRNAs were identified as putative targets and target mimics of known Populus miRNAs, respectively. A total of 504 lincRNAs were found to be drought responsive, eight of which were confirmed by RT-qPCR. These findings provide a comprehensive view of Populus lincRNAs, which will enable in-depth functional analysis.

Trans-acting small interfering RNA4: key to nutraceutical synthesis in grape development?

The facility and versatility of microRNAs (miRNAs) to evolve and change likely underlies how they have become dominant constituents of eukaryotic genomes. In this opinion article I propose that trans-acting small interfering RNA gene 4 (TAS4) evolution may be important for biosynthesis of polyphenolics, arbuscular symbiosis, and bacterial pathogen etiologies. Expression-based and phylogenetic evidence shows that TAS4 targets two novel grape (Vitis vinifera L.) MYB transcription factors (VvMYBA6, VvMYBA7) that spawn phased small interfering RNAs (siRNAs) which probably function in nutraceutical bioflavonoid biosynthesis and fruit development. Characterization of the molecular mechanisms of TAS4 control of plant development and integration into biotic and abiotic stress- and nutrient-signaling regulatory networks has applicability to molecular breeding and the development of strategies for engineering healthier foods.

Widespread long noncoding RNAs as endogenous target mimics for microRNAs in plants

Target mimicry is a recently identified regulatory mechanism for microRNA (miRNA) functions in plants in which the decoy RNAs bind to miRNAs via complementary sequences and therefore block the interaction between miRNAs and their authentic targets. Both endogenous decoy RNAs (miRNA target mimics) and engineered artificial RNAs can induce target mimicry effects. Yet until now, only the Induced by Phosphate Starvation1 RNA has been proven to be a functional endogenous microRNA target mimic (eTM). In this work, we developed a computational method and systematically identified intergenic or noncoding gene-originated eTMs for 20 conserved miRNAs in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa). The predicted miRNA binding sites were well conserved among eTMs of the same miRNA, whereas sequences outside of the binding sites varied a lot. We proved that the eTMs of miR160 and miR166 are functional target mimics and identified their roles in the regulation of plant development. The effectiveness of eTMs for three other miRNAs was also confirmed by transient agroinfiltration assay.