Identification and characterization of long non-coding RNA (lncRNA) in the developing seeds of Jatropha curcas

CcNFYB3-CcMATE35 and LncRNA CcLTCS-CcCS modules jointly regulate the efflux and synthesis of citrate to enhance aluminium tolerance in pigeon pea

Aluminium (Al) toxicity decreases crop production in acid soils in general, but many crops have evolved complex mechanisms to resist it. However, our current understanding of how plants cope with Al stress and perform Al resistance is still at the initial stage. In this study, the citrate transporter CcMATE35 was identified to be involved in Al stress response. The release of citrate was increased substantially in CcMATE35 over-expression (OE) lines under Al stress, indicating enhanced Al resistance. It was demonstrated that transcription factor CcNFYB3 regulated the expression of CcMATE35, promoting the release of citrate from roots to increase Al resistance in pigeon pea. We also found that a Long noncoding RNA Targeting Citrate Synthase (CcLTCS) is involved in Al resistance in pigeon pea. Compared with controls, overexpression of CcLTCS elevated the expression level of the Citrate Synthase gene (CcCS), leading to increases in root citrate level and citrate release, which forms another module to regulate Al resistance in pigeon pea. Simultaneous overexpression of CcNFYB3 and CcLTCS further increased Al resistance. Taken together, these findings suggest that the two modules, CcNFYB3-CcMATE35 and CcLTCS-CcCS, jointly regulate the efflux and synthesis of citrate and may play an important role in enhancing the resistance of pigeon pea under Al stress.

Genome-wide identification of long non-coding RNAs and their potential functions in radish response to salt stress

Long non-coding RNAs (lncRNAs) are increasingly recognized as cis- and trans-acting regulators of protein-coding genes in plants, particularly in response to abiotic stressors. Among these stressors, high soil salinity poses a significant challenge to crop productivity. Radish (Raphanus sativus L.) is a prominent root vegetable crop that exhibits moderate susceptibility to salt stress, particularly during the seedling stage. Nevertheless, the precise regulatory mechanisms through which lncRNAs contribute to salt response in radish remain largely unexplored. In this study, we performed genome-wide identification of lncRNAs using strand-specific RNA sequencing on radish fleshy root samples subjected to varying time points of salinity treatment. A total of 7,709 novel lncRNAs were identified, with 363 of them displaying significant differential expression in response to salt application. Furthermore, through target gene prediction, 5,006 cis- and 5,983 trans-target genes were obtained for the differentially expressed lncRNAs. The predicted target genes of these salt-responsive lncRNAs exhibited strong associations with various plant defense mechanisms, including signal perception and transduction, transcription regulation, ion homeostasis, osmoregulation, reactive oxygen species scavenging, photosynthesis, phytohormone regulation, and kinase activity. Notably, this study represents the first comprehensive genome-wide analysis of salt-responsive lncRNAs in radish, to the best of our knowledge. These findings provide a basis for future functional analysis of lncRNAs implicated in the defense response of radish against high salinity, which will aid in further understanding the regulatory mechanisms underlying radish response to salt stress.

Emerging roles of long non-coding RNAs in regulating agriculturally important seed traits

Seeds have enormous economic importance as source of calories, nutrition, edible oil, and biofuels. Therefore, seed traits like seed size and shape, weight, micronutrient content, oil content, quality, post-harvest sprouting, etc., are some of the main targets in crop improvement. Designing the strategies for their improvement benefits heavily from understanding the regulatory aspects of seed development. Recent studies indicate that long non-coding RNAs (lncRNAs) are one of the important regulators of seed development. They played a significant role in crop domestication by influencing seed traits. LncRNAs are conventionally defined as non-coding RNAs greater than 200 bp in length but lacking protein coding potential. Here we highlight the emerging pieces of evidence of lncRNA-mediated regulation of seed development through diverse mechanisms, for instance, by acting as target mimics or precursors of regulatory small RNAs or through chromatin remodeling and post-transcriptional repression. We also enumerate the insights from high-throughput transcriptomic studies from developing seeds of cereal, oilseed, biofuel, and pulse crops. We highlight the lncRNA candidates and lncRNA-mediated regulatory networks regulating seed development and related agronomic traits. Further, we discuss the potential of lncRNAs for improvement of agriculturally important seed traits through marker-assisted breeding and/or transgenic approaches.

Interaction between long noncoding RNA (lnc663) and microRNA (miR1128) regulates PDAT-like gene activity in bread wheat (Triticum aestivum L.)

Amylose, a starch subcomponent, can bind lipids within its helical groove and form an amylose-lipid complex, known as resistant starch type 5 (RS-5). RS contributes to lower glycaemic index of grain with health benefits. Unfortunately, genes involved in lipid biosynthesis in wheat grain remain elusive. Our study aims to characterize the lipid biosynthesis gene and its post-transcriptional regulation using the parent bread wheat variety 'C 306' and its EMS-induced mutant line 'TAC 75' varying in amylose content. Quantitative analyses of starch-bound lipids showed that 'TAC 75' has significantly higher lipid content in grains than 'C 306' variety. Furthermore, expression analyses revealed the higher expression of wheat phospholipid: diacylglycerol acyltransferase-like (PDAT-like) in the 'TAC 75' compared to the 'C 306'. Overexpression and ectopic expression of TaPDAT in yeast and tobacco leaf confirmed its ability to accumulate lipids in vivo. Enzyme activity assay showed that TaPDAT catalyzes the triacylglycerol synthesis by acylating 1,2-diacylglycerol. Interestingly, the long non-coding RNA, lnc663, was upregulated with the TaPDAT gene, while the miRNA, miR1128, downregulated in the 'TAC 75', indicating a regulatory relationship. The GFP reporter assay confirmed that the lnc663 acts as a positive regulator, and the miR1128 as a negative regulator of the TaPDAT gene, which controls lipid accumulation in wheat grain. Our findings outline TaPDAT-mediated biosynthesis of lipid accumulation and reveal the molecular mechanism of the lnc663 and miR1128 mediated regulation of the TaPDAT gene in wheat grain.

Integrated 3D genome, epigenome and transcriptome analyses reveal transcriptional coordination of circadian rhythm in rice

Photoperiods integrate with the circadian clock to coordinate gene expression rhythms and thus ensure plant fitness to the environment. Genome-wide characterization and comparison of rhythmic genes under different light conditions revealed delayed phase under constant darkness (DD) and reduced amplitude under constant light (LL) in rice. Interestingly, ChIP-seq and RNA-seq profiling of rhythmic genes exhibit synchronous circadian oscillation in H3K9ac modifications at their loci and long non-coding RNAs (lncRNAs) expression at proximal loci. To investigate how gene expression rhythm is regulated in rice, we profiled the open chromatin regions and transcription factor (TF) footprints by time-series ATAC-seq. Although open chromatin regions did not show circadian change, a significant number of TFs were identified to rhythmically associate with chromatin and drive gene expression in a time-dependent manner. Further transcriptional regulatory networks mapping uncovered significant correlation between core clock genes and transcription factors involved in light/temperature signaling. In situ Hi-C of ZT8-specific expressed genes displayed highly connected chromatin association at the same time, whereas this ZT8 chromatin connection network dissociates at ZT20, suggesting the circadian control of gene expression by dynamic spatial chromatin conformation. These findings together implicate the existence of a synchronization mechanism between circadian H3K9ac modifications, chromatin association of TF and gene expression, and provides insights into circadian dynamics of spatial chromatin conformation that associate with gene expression rhythms.

Genome-wide identification and characterization of Puccinia striiformis-responsive lncRNAs in Triticum aestivum

Wheat stripe rust (yellow rust) caused by Puccinia striiformis f. sp. tritici (Pst) is a serious biotic stress factor limiting wheat production worldwide. Emerging evidence demonstrates that long non-coding RNAs (lncRNAs) participate in various developmental processes in plants via post-transcription regulation. In this study, RNA sequencing (RNA-seq) was performed on a pair of near-isogenic lines-rust resistance line FLW29 and rust susceptible line PBW343-which differed only in the rust susceptibility trait. A total of 6,807 lncRNA transcripts were identified using bioinformatics analyses, among which 10 lncRNAs were found to be differentially expressed between resistance and susceptible lines. In order to find the target genes of the identified lncRNAs, their interactions with wheat microRNA (miRNAs) were predicted. A total of 199 lncRNAs showed interactions with 65 miRNAs, which further target 757 distinct mRNA transcripts. Moreover, detailed functional annotations of the target genes were used to identify the candidate genes, pathways, domains, families, and transcription factors that may be related to stripe rust resistance response in wheat plants. The NAC domain protein, disease resistance proteins RPP13 and RPM1, At1g58400, monodehydroascorbate reductase, NBS-LRR-like protein, rust resistance kinase Lr10-like, LRR receptor, serine/threonine-protein kinase, and cysteine proteinase are among the identified targets that are crucial for wheat stripe rust resistance. Semiquantitative PCR analysis of some of the differentially expressed lncRNAs revealed variations in expression profiles of two lncRNAs between the Pst-resistant and Pst-susceptible genotypes at least under one condition. Additionally, simple sequence repeats (SSRs) were also identified from wheat lncRNA sequences, which may be very useful for conducting targeted gene mapping studies of stripe rust resistance in wheat. These findings improved our understanding of the molecular mechanism responsible for the stripe rust disease that can be further utilized to develop wheat varieties with durable resistance to this disease.

ICAnnoLncRNA: A Snakemake Pipeline for a Long Non-Coding-RNA Search and Annotation in Transcriptomic Sequences

Long non-coding RNAs (lncRNAs) are RNA molecules longer than 200 nucleotides that do not encode proteins. Experimental studies have shown the diversity and importance of lncRNA functions in plants. To expand knowledge about lncRNAs in other species, computational pipelines that allow for standardised data-processing steps in a mode that does not require user control up until the final result were actively developed recently. These advancements enable wider functionality for lncRNA data identification and analysis. In the present work, we propose the ICAnnoLncRNA pipeline for the automatic identification, classification and annotation of plant lncRNAs in assembled transcriptomic sequences. It uses the LncFinder software for the identification of lncRNAs and allows the adjustment of recognition parameters using genomic data for which lncRNA annotation is available. The pipeline allows the prediction of lncRNA candidates, alignment of lncRNA sequences to the reference genome, filtering of erroneous/noise transcripts and probable transposable elements, lncRNA classification by genome location, comparison with sequences from external databases and analysis of lncRNA structural features and expression. We used transcriptomic sequences from 15 maize libraries assembled by Trinity and Hisat2/StringTie to demonstrate the application of the ICAnnoLncRNA pipeline.

Long non-coding RNA and microRNA landscape of two major domesticated cotton species

Allotetraploid cotton plants Gossypium hirsutum and Gossypium barbadense have been widely cultivated for their natural, renewable textile fibres. Even though ncRNAs in domesticated cotton species have been extensively studied, systematic identification and annotation of lncRNAs and miRNAs expressed in various tissues and developmental stages under various biological contexts are limited. This influences the comprehension of their functions and future research on these cotton species. Here, we report high confidence lncRNAs and miRNA collection from G. hirsutum accession and G. barbadense accession using large-scale RNA-seq and small RNA-seq datasets incorporated into a user-friendly database, CoNCRAtlas. This database provides a wide range and depth of lncRNA and miRNA annotation based on the systematic integration of extensive annotations such as expression patterns derived from transcriptome data analysis in thousands of samples, as well as multi-omics annotations. We assume this comprehensive resource will accelerate evolutionary and functional studies in ncRNAs and inform future breeding programs for cotton improvement. CoNCRAtlas is accessible at http://www.nipgr.ac.in/CoNCRAtlas/.

Genomic survey of high-throughput RNA-Seq data implicates involvement of long intergenic non-coding RNAs (lincRNAs) in cytoplasmic male-sterility and fertility restoration in pigeon pea

BACKGROUND

Long-intergenic non-coding RNAs (lincRNAs) originate from intergenic regions and have no coding potential. LincRNAs have emerged as key players in the regulation of various biological processes in plant development. Cytoplasmic male-sterility (CMS) in association with restorer-of-fertility (Rf) systems makes it a highly reliable tool for exploring heterosis for producing commercial hybrid seeds. To date, there have been no reports of lincRNAs during pollen development in CMS and fertility restorer lines in pigeon pea.

OBJECTIVE

Identification of lincRNAs in the floral buds of cytoplasmic male-sterile (AKCMS11) and fertility restorer (AKPR303) pigeon pea lines.

METHODS

We employed a computational approach to identify lincRNAs in the floral buds of cytoplasmic male-sterile (AKCMS11) and fertility restorer (AKPR303) pigeon pea lines using RNA-Seq data.

RESULTS

We predicted a total of 2145 potential lincRNAs of which 966 were observed to be differentially expressed between the sterile and fertile pollen. We identified, 927 cis-regulated and 383 trans-regulated target genes of the lincRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the target genes revealed that these genes were specifically enriched in pathways like pollen and pollen tube development, oxidative phosphorylation, etc. We detected 23 lincRNAs that were co-expressed with 17 pollen-related genes with known functions. Fifty-nine lincRNAs were predicted to be endogenous target mimics (eTMs) for 25 miRNAs, and found to be associated with pollen development. The, lincRNA regulatory networks revealed that different lincRNA-miRNA-mRNA networks might be associated with CMS and fertility restoration.

CONCLUSION

Thus, this study provides valuable information by highlighting the functions of lincRNAs as regulators during pollen development in pigeon pea and utilization in hybrid seed production.

Genome-wide identification and characterization of tissue-specific non-coding RNAs in black pepper (Piper nigrum L.)

Long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) are the two classes of non-coding RNAs (ncRNAs) present predominantly in plant cells and have various gene regulatory functions at pre- and post-transcriptional levels. Previously deemed as "junk", these ncRNAs have now been reported to be an important player in gene expression regulation, especially in stress conditions in many plant species. Black pepper, scientifically known as Piper nigrum L., despite being one of the most economically important spice crops, lacks studies related to these ncRNAs. From a panel of 53 RNA-Seq datasets of black pepper from six tissues, namely, flower, fruit, leaf, panicle, root, and stem of six black pepper cultivars, covering eight BioProjects across four countries, we identified and characterized a total of 6406 lncRNAs. Further downstream analysis inferred that these lncRNAs regulated 781 black pepper genes/gene products via miRNA-lncRNA-mRNA network interactions, thus working as competitive endogenous RNAs (ceRNAs). The interactions may be various mechanisms like miRNA-mediated gene silencing or lncRNAs acting as endogenous target mimics (eTMs) of the miRNAs. A total of 35 lncRNAs were also identified to be potential precursors of 94 miRNAs after being acted upon by endonucleases like Drosha and Dicer. Tissue-wise transcriptome analysis revealed 4621 circRNAs. Further, miRNA-circRNA-mRNA network analysis showed 432 circRNAs combining with 619 miRNAs and competing for the binding sites on 744 mRNAs in different black pepper tissues. These findings can help researchers to get a better insight to the yield regulation and responses to stress in black pepper in endeavor of higher production and improved breeding programs in black pepper varieties.

Genome-wide screening and characterization of long noncoding RNAs involved in flowering/bolting of Lactuca sativa

BACKGROUND

Lettuce (Lactuca sativa L.) is considered the most important vegetable in the leafy vegetable group. However, bolting affects quality, gives it a bitter taste, and as a result makes it inedible. Bolting is an event induced by the coordinated effects of various environmental factors and endogenous genetic components. Although bolting/flowering responsive genes have been identified in most sensitive and non-sensitive species, non-coding RNA molecules like long non-coding RNAs (lncRNAs) have not been investigated in lettuce. Hence, in this study, potential long non-coding RNAs that regulate flowering /bolting were investigated in two lettuce strains S24 (resistant strain) and S39 (susceptible strain) in different flowering times to better understand the regulation of lettuce bolting mechanism. For this purpose, we used two RNA-seq datasets to discover the lncRNA transcriptome profile during the transition from vegetative to reproductive phase.

RESULTS

For identifying unannotated transcripts in these datasets, a 7-step pipeline was employed to filter out these transcripts and terminate with 293 novel lncRNAs predicted by PLncPRO and CREMA. These transcripts were then utilized to predict cis and trans flowering-associated targets and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Computational predictions of target gene function showed the involvement of putative flowering-related genes and enrichment of the floral regulators FLC, CO, FT, and SOC1 in both datasets. Finally, 17 and 18 lncRNAs were proposed as competing endogenous target mimics (eTMs) for novel and known lncRNA miRNAs, respectively.

CONCLUSION

Overall, this study provides new insights into lncRNAs that control the flowering time of plants known for bolting, such as lettuce, and opens new windows for further study.

Identification of a Novel Long Non-coding RNA, lnc-ATMIN-4:2, and its Clinicopathological and Prognostic Significance in Advanced Gastric Cancer

BACKGROUND/AIM

Long non-coding RNAs (lncRNAs) are emerging as significant regulators of gene expression and a novel promising biomarker for cancer diagnosis and prognosis. This study identified a novel, differentially expressed lncRNA in advanced gastric cancer (AGC), Inc-ATMIN-4:2, and evaluated its clinicopathological and prognostic significance.

PATIENTS AND METHODS

Whole transcriptome sequencing was performed to identify differentially expressed lncRNAs in AGC tissue samples. We also analyzed lnc-ATMIN-4:2 expression in 317 patients with AGC using RNA in situ hybridization.

RESULTS

High (>30 dots) lnc-ATMIN-4:2 expression significantly correlated with younger age, poorly differentiated histology, diffuse type, deeper invasion depth, perineural invasion, lymph node metastasis, and higher stage group. In addition, high lnc-ATMIN-4:2 expression was significantly associated with worse overall survival in patients with AGC.

CONCLUSION

This study elucidated the significance of lncRNAs in AGC and indicated the value of lnc-ATMIN-4:2 expression as a predictive biomarker for the overall survival of patients with AGC.

Genome-wide analysis suggests the potential role of lncRNAs during seed development and seed size/weight determination in chickpea

The integrated transcriptome data analyses suggested the plausible roles of lncRNAs during seed development in chickpea. The candidate lncRNAs associated with QTLs and those involved in miRNA-mediated seed size/weight determination in chickpea have been identified. Long non-coding RNAs (lncRNAs) are important regulators of various biological processes. Here, we identified lncRNAs at seven successive stages of seed development in small-seeded and large-seeded chickpea cultivars. In total, 4751 lncRNAs implicated in diverse biological processes were identified. Most of lncRNAs were conserved between the two cultivars, whereas only a few of them were conserved in other plants, suggesting their species-specificity. A large number of lncRNAs differentially expressed between the two chickpea cultivars associated with seed development-related processes were identified. The lncRNAs acting as precursors of miRNAs and those mimicking target protein-coding genes of miRNAs involved in seed size/weight determination, including HAIKU1, BIG SEEDS1, and SHB1, were also revealed. Further, lncRNAs located within seed size/weight associated quantitative trait loci were also detected. Overall, we present a comprehensive resource and identified candidate lncRNAs that may play important roles during seed development and seed size/weight determination in chickpea.

Long Non-Coding RNAs: New Players in Plants

During the process of growth and development, plants are prone to various biotic and abiotic stresses. They have evolved a variety of strategies to resist the adverse effects of these stresses. lncRNAs (long non-coding RNAs) are a type of less conserved RNA molecules of more than 200 nt (nucleotides) in length. lncRNAs do not code for any protein, but interact with DNA, RNA, and protein to affect transcriptional, posttranscriptional, and epigenetic modulation events. As a new regulatory element, lncRNAs play a critical role in coping with environmental pressure during plant growth and development. This article presents a comprehensive review on the types of plant lncRNAs, the role and mechanism of lncRNAs at different molecular levels, the coordination between lncRNA and miRNA (microRNA) in plant immune responses, the latest research progress of lncRNAs in plant growth and development, and their response to biotic and abiotic stresses. We conclude with a discussion on future direction for the elaboration of the function and mechanism of lncRNAs.

Exosomal long noncoding RNAs - the lead thespian behind the regulation, cause and cure of autophagy-related diseases

Recent advances in exosome biology have revealed significant roles of exosome and their contents in intercellular communication. Among various exosomal content, long non-coding RNAs (lncRNAs), which have a large size (˃ 200 nt) and lack protein coding potential, are known to play key roles in intercellular communication and novel biomarkers of various metabolic disorders. Moreover, long non-coding RNAs are often involved in the regulation of various cellular processes such as autophagy, apoptosis, cell proliferation. On the other hand, autophagy is the central regulating point that controls the various metabolic functions of the body. This process is known to prevent diseases and promote longevity. Therefore, the present review discusses the relationship between diseases and autophagy, and also look into the biological functions of exosome-associated lncRNAs in regulating autophagy. Furthermore, this review will summarize some of the studies that provide novel insights into the pathogenesis of autophagy-related diseases followed by the non-canonical roles played by autophagy and related proteins in the development of exosome biogenesis.

Comprehensive Analysis of Long Non-coding RNA Modulates Axillary Bud Development in Tobacco (Nicotiana tabacum L.)

Long non-coding RNAs (lncRNAs) regulate gene expression and are crucial for plant growth and development. However, the mechanisms underlying the effects of activated lncRNAs on axillary bud development remain largely unknown. By lncRNA transcriptomes of axillary buds in topped and untopped tobacco plants, we identified a total of 13,694 lncRNAs. LncRNA analysis indicated that the promoted growth of axillary bud by topping might be partially ascribed to the genes related to hormone signal transduction and glycometabolism, trans-regulated by differentially expressed lncRNAs, such as MSTRG.52498.1, MSTRG.60026.1, MSTRG.17770.1, and MSTRG.32431.1. Metabolite profiling indicated that auxin, abscisic acid and gibberellin were decreased in axillary buds of topped tobacco lines, while cytokinin was increased, consistent with the expression levels of related lncRNAs. MSTRG.52498.1, MSTRG.60026.1, MSTRG.17770.1, and MSTRG.32431.1 were shown to be influenced by hormones and sucrose treatments, and were associated with changes of axillary bud growth in the overexpression of NtCCD8 plants (with reduced axillary buds) and RNA interference of NtTB1 plants (with increased axillary buds). Moreover, MSTRG.28151.1 was identified as the antisense lncRNA of NtTB1. Silencing of MSTRG.28151.1 in tobacco significantly attenuated the expression of NtTB1 and resulted in larger axillary buds, suggesting the vital function of MSTRG.28151.1 axillary bud developmen by NtTB1. Our findings shed light on lncRNA-mRNA interactions and their functional roles in axillary bud growth, which would improve our understanding of lncRNAs as important regulators of axillary bud development and plant architecture.

Identification of the Complex Interplay Between Nematode-Related lncRNAs and Their Target Genes in Glycine max L

Soybean (Glycine max) is a major plant protein source and oilseed crop. However, plant-parasitic nematodes (PPNs) affect its annual yield. In the current study, in order to better understand the regulation of defense mechanism against PPNs in soybean, we investigated the role of long non-coding RNAs (lncRNAs) in response to two nematode species, Heterodera glycines (SCN: soybean cyst nematode) and Rotylenchulus reniformis (reniform). To this end, two publicly available RNA-seq data sets (SCN data set and RAD: reniform-associated data set) were employed to discover the lncRNAome profile of soybean under SCN and reniform infection, respectively. Upon identification of unannotated transcripts in these data sets, a seven-step pipeline was utilized to sieve these transcripts, which ended up in 384 and 283 potential lncRNAs in SCN data set and RAD, respectively. These transcripts were then used to predict cis and trans nematode-related targets in soybean genome. Computational prediction of target genes function, some of which were also among differentially expressed genes, revealed the involvement of putative nematode-responsive genes as well as enrichment of multiple stress responses in both data sets. Finally, 15 and six lncRNAs were proposed to be involved in microRNA-mediated regulation of gene expression in soybean in response to SNC and reniform infection, respectively. Collectively, this study provides a novel insight into the signaling and regulatory network of soybean-pathogen interactions and opens a new window for further research.

Biologia Futura: progress and future perspectives of long non-coding RNAs in forest trees

Forest trees are affected by climate change, anthropogenic pressure, as well as abiotic and biotic stresses. Conventional tree breeding has so far been limited to enhance overall productivity, and our understanding of the genetic basis of quantitative traits is still inadequate. Quantum leaps in next-generation sequencing technologies and bioinformatics have permitted the exploration and identification of various non-coding regions of the genome other than protein coding genes. These genomic regions produce various types of non-coding RNAs and regulate myriads of biological functions at epigenetic, transcriptional and translational levels. Recently, long non-coding RNAs (lncRNAs) which act as molecular switch have been identified to be pivotal molecules in forest trees. This review focuses on progress made in regulatory mechanisms in various developmental phases like wood formation, adventitious rooting and flowering and stress responses. It was predicted that complex regulatory interactions among lncRNA, miRNA and gene exist. LncRNAs can function as a sponge for miRNAs, reducing the suppressive effect of miRNAs on target mRNAs and perhaps adding a new layer of regulatory interactions among non-coding RNA classes in trees. Furthermore, network analysis revealed the interactions of lncRNA and genes during the expression of several important genes. The insights generated about lncRNAs in forest trees would enable improvement of economically important traits including the devastating abiotic and biotic stresses. In addition, solid understanding on the wide range of regulatory functions of lncRNAs on traits influencing biomass productivity and adaptation would aid the applications of biotechnology in genetic improvement of forest trees.