Genome-Wide Identification and Characterization of Long Non-Coding RNAs from Mulberry (Morus notabilis) RNA-seq Data

Emerging roles and mechanisms of lncRNAs in fruit and vegetables

With the development of genome sequencing technologies, many long non-coding RNAs (lncRNAs) have been identified in fruit and vegetables. lncRNAs are primarily transcribed and spliced by RNA polymerase II (Pol II) or plant-specific Pol IV/V, and exhibit limited evolutionary conservation. lncRNAs intricately regulate various aspects of fruit and vegetables, including pigment accumulation, reproductive tissue development, fruit ripening, and responses to biotic and abiotic stresses, through diverse mechanisms such as gene expression modulation, interaction with hormones and transcription factors, microRNA regulation, and involvement in alternative splicing. This review presents a comprehensive overview of lncRNA classification, basic characteristics, and, most importantly, recent advances in understanding their functions and regulatory mechanisms.

Long noncoding RNA transcriptome analysis reveals novel lncRNAs in Morus alba 'Yu-711' response to drought stress

Drought stress has been a key environmental factor affecting plant growth and development. The plant genome is capable of producing long noncoding RNAs (lncRNAs). To better understand white mulberry (Morus alba L.) drought response mechanism, we conducted a comparative transcriptome study comparing two treatments: drought-stressed (EG) and well-watered (CK) plants. A total of 674 differentially expressed lncRNAs (DElncRNAs) were identified. In addition, 782 differentially expressed messenger RNAs (DEmRNAs) were identified. We conducted Gene Ontology (GO) and KEGG enrichment analyses focusing on the differential lncRNAs cis-target genes. The target genes of the DElncRNAs were most significantly involved in the biosynthesis of secondary metabolites. Gene regulatory networks of the target genes involving DElncRNAs-mRNAs-DEmRNAs and DElncRNA-miRNA-DEmRNA were constructed. In the DElncRNAs-DEmRNAs network, 30 DEmRNAs involved in the biosynthesis of secondary metabolites are collocated with 46 DElncRNAs. The interaction between DElncRNAs and candidate genes was identified using LncTar. In summary, quantitative real-time polymerase chain reaction (qRT-PCR) validated nine candidate genes and seven target lncRNAs including those identified by LncTar. We predicted that the DElncRNAs-DEmRNAs might recruit microRNAs (miRNAs) to interact with gene regulatory networks under the drought stress response in mulberry. The findings will contribute to our understanding of the regulatory functions of lncRNAs under drought stress and will shed new light on the mulberry-drought stress interactions.

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.

SMRT and Illumina RNA sequencing reveal the complexity of terpenoid biosynthesis in Zanthoxylum armatum

Sichuan pepper (Zanthoxylum armatum DC) is a popular spice and is often prescribed in traditional Chinese medicine to treat vomiting, diarrhea, ascariasis and eczema, among other conditions. Volatile oils from Z. armatum leaves contain active ingredients, with terpenoids being one of the main components. In the present study, the combination of sequencing data of Z. armatum from PacBio single molecule real time (SMRT) and Illumina RNA sequencing (RNA-Seq) platforms facilitated an understanding of the gene regulatory network of terpenoid biosynthesis in pepper leaves. The leaves of three developmental stages from two Z. armatum cultivars, 'Rongchangwuci' (WC) and 'Zhuye' (ZY), were selected as test materials to construct sequencing libraries. A total of 143,122 predictions of unique coding sequences, 105,465 simple sequence repeats, 20,145 transcription factors and 4719 long non-coding RNAs (lncRNAs) were identified, and 142,829 transcripts were successfully annotated. The occurrence of alternative splicing events was verified by reverse transcription PCR, and quantitative real-time PCR was used to confirm the expression pattern of six randomly selected lncRNAs. A total of 96,931 differentially expressed genes were filtered from different samples. According to functional annotation, a total of 560 candidate genes were involved in terpenoid synthesis, of which 526 were differentially expressed genes (DEGs). To identify the key genes involved in terpenoid biosynthesis, the module genes in different samples, including structural and transcription factors genes, were analyzed using the weighted gene co-expression network method, and the co-expression network of genes was constructed. Thirty-one terpenoids were identified by gas chromatography-mass spectrometry. The correlation between 18 compounds with significantly different contents and genes with high connectivity in the module was jointly analyzed in both cultivars, yielding 12 candidate DEGs presumably involved in the regulation of terpenoid biosynthesis. Our findings showed that full-length transcriptome SMRT and Illumina RNA-Seq can play an important role in studying organisms without reference genomes and elucidating the gene regulation of a biosynthetic pathway.

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.

Alternative splicing profiling provides insights into the molecular mechanisms of peanut peg development

BACKGROUND

The cultivated peanut (Arachis hypogaea) is one of the most important oilseed crops worldwide, and the generation of pegs and formation of subterranean pods are essential processes in peanut reproductive development. However, little information has been reported about alternative splicing (AS) in peanut peg formation and development.

RESULTS

Herein, we presented a comprehensive full-length (FL) transcriptome profiling of AS isoforms during peanut peg and early pod development. We identified 1448, 1102, 832, and 902 specific spliced transcripts in aerial pegs, subterranean pegs, subterranean unswollen pegs, and early swelling pods, respectively. A total of 184 spliced transcripts related to gravity stimulation, light and mechanical response, hormone mediated signaling pathways, and calcium-dependent proteins were identified as possibly involved in peanut peg development. For aerial pegs, spliced transcripts we got were mainly involved in gravity stimulation and cell wall morphogenetic processes. The genes undergoing AS in subterranean peg were possibly involved in gravity stimulation, cell wall morphogenetic processes, and abiotic response. For subterranean unswollen pegs, spliced transcripts were predominantly related to the embryo development and root formation. The genes undergoing splice in early swelling pods were mainly related to ovule development, root hair cells enlargement, root apex division, and seed germination.

CONCLUSION

This study provides evidence that multiple genes are related to gravity stimulation, light and mechanical response, hormone mediated signaling pathways, and calcium-dependent proteins undergoing AS express development-specific spliced isoforms or exhibit an obvious isoform switch during the peanut peg development. AS isoforms in subterranean pegs and pods provides valuable sources to further understand post-transcriptional regulatory mechanisms of AS in the generation of pegs and formation of subterranean pods.

The expression of long non-coding RNAs is associated with H3Ac and H3K4me2 changes regulated by the HDA6-LDL1/2 histone modification complex in Arabidopsis

In recent years, eukaryotic long non-coding RNAs (lncRNAs) have been identified as important factors involved in a wide variety of biological processes, including histone modification, alternative splicing and transcription enhancement. The expression of lncRNAs is highly tissue-specific and is regulated by environmental stresses. Recently, a large number of plant lncRNAs have been identified, but very few of them have been studied in detail. Furthermore, the mechanism of lncRNA expression regulation remains largely unknown. Arabidopsis HISTONE DEACETYLASE 6 (HDA6) and LSD1-LIKE 1/2 (LDL1/2) can repress gene expression synergistically by regulating H3Ac/H3K4me. In this research, we performed RNA-seq and ChIP-seq analyses to further clarify the function of HDA6-LDL1/2. Our results indicated that the global expression of lncRNAs is increased in hda6/ldl1/2 and that this increased lncRNA expression is particularly associated with H3Ac/H3K4me2 changes. In addition, we found that HDA6-LDL1/2 is important for repressing lncRNAs that are non-expressed or show low-expression, which may be strongly associated with plant development. GO-enrichment analysis also revealed that the neighboring genes of the lncRNAs that are upregulated in hda6/ldl1/2 are associated with various developmental processes. Collectively, our results revealed that the expression of lncRNAs is associated with H3Ac/H3K4me2 changes regulated by the HDA6-LDL1/2 histone modification complex.

PreLnc: An Accurate Tool for Predicting lncRNAs Based on Multiple Features

Accumulating evidence indicates that long non-coding RNAs (lncRNAs) have certain similarities with messenger RNAs (mRNAs) and are associated with numerous important biological processes, thereby demanding methods to distinguish them. Based on machine learning algorithms, a variety of methods are developed to identify lncRNAs, providing significant basic data support for subsequent studies. However, many tools lack certain scalability, versatility and balance, and some tools rely on genome sequence and annotation. In this paper, we propose a convenient and accurate tool "PreLnc", which uses high-confidence lncRNA and mRNA transcripts to build prediction models through feature selection and classifiers. The false discovery rate (FDR) adjusted P-value and Z-value were used for analyzing the tri-nucleotide composition of transcripts of different species. Conclusions can be drawn from the experiment that there were significant differences in RNA transcripts among plants, which may be related to evolutionary conservation and the fact that plants are under evolutionary pressure for a longer time than animals. Combining with the Pearson correlation coefficient, we use the incremental feature selection (IFS) method and the comparison of multiple classifiers to build the model. Finally, the balanced random forest was used to construct the classifier, and PreLnc obtained 91.09% accuracy for 349,186 transcripts of animals and plants. In addition, by comparing standard performance measurements, PreLnc performed better than other prediction tools.

Discovery, identification, and functional characterization of long noncoding RNAs in Arachis hypogaea L

BACKGROUND

Long noncoding RNAs (lncRNAs), which are typically > 200 nt in length, are involved in numerous biological processes. Studies on lncRNAs in the cultivated peanut (Arachis hypogaea L.) largely remain unknown.

RESULTS

A genome-wide scan of the peanut (Arachis hypogaea L.) transcriptome identified 1442 lncRNAs, which were encoded by loci distributed over every chromosome. Long intergenic noncoding RNAs accounted for 85.58% of these lncRNAs. Additionally, 189 lncRNAs were differentially abundant in the root, leaf, or seed. Generally, lncRNAs showed lower expression levels, tighter tissue-specific expression, and less splicing than mRNAs. Approximately 44.17% of the lncRNAs with an exon/intron structure were alternatively spliced; this rate was slightly lower than the splicing rate of mRNA. Transcription at the start site event was the alternative splicing (AS) event with the highest frequency (28.05%) in peanut lncRNAs, whereas the occurrence rate (30.19%) of intron retention event was the highest in mRNAs. AS changed the target gene profiles of lncRNAs and increased the diversity and flexibility of lncRNAs, which may be important for lncRNAs to execute their functions. Additionally, a substantial number of the peanut AS isoforms generated from protein-encoding genes appeared to be noncoding because they were truncated transcripts; such isoforms can be legitimately regarded as a class of lncRNAs. The predicted target genes of the lncRNAs were involved in a wide range of biological processes. Furthermore, expression pattern of several selected lncRNAs and their target genes were examined under salt stress, results showed that all of them could respond to salt stress in different manners.

CONCLUSIONS

This study provided a resource of candidate lncRNAs and expression patterns across tissues, and whether these lncRNAs are functional will be further investigated in our subsequent experiments.

Genome-Wide Identification of WOX Gene Family and Expression Analysis during Rejuvenational Rhizogenesis in Walnut (Juglans regia L.)

Rejuvenation is an efficient approach used in the cuttings of trees and horticultural crops, to improve their rooting ability, especially in difficult-to-root trees. WOX gene family members are involved in cell-fate transformation through balancing the maintenance and proliferation of the stem cells. However, there are no reports about the WOX gene family in Walnut (Juglans regia L.) and its relationship between rejuvenation and adventitious roots formation (ARF). Here, a genome-wide identification of JrWOX genes and their physical and chemical properties, phylogeny, and expression profiles in different organs and during rejuvenation-induced ARF is reported. The phenotype and histology characteristics of mature and rejuvenated cuttings (Mc and Rc) are also observed. In this study, 12 genes were identified and clustered into three groups based on phylogenetics, special domains, and conserved motifs. The gene structures and conserved motifs were relatively conserved, while the 12 sequences of the JrWOXs domain were diversified. Gene expression in root, stem, leaf, female flower, immature fruit, and zygotic embryo revealed that the expression levels of JrWOX4a, JrWOX4b, JrWOX5, JrWOX11, and JrWOX13 in the root were significantly higher than those of other JrWOXs, while only the expression of JrWOX11 was exclusive to the root organ. Additionally, rejuvenation treatment significantly induced almost all JrWOX genes, except JrWOX4a, JrWOX4b, and JrWOX13 (Rc 0 vs. Mc 0). During the ARF process, the transcripts of JrWOX11 and JrWOX5 were consecutively increased on a significance level; in contrast, the transcription levels of the other JrWOXs decreased or changed insignificantly. The phenotype and histology observation indicate that rejuvenation treatment made the base of the stem expand and reduced the thickness and density of sclerenchyma between the cortex and phloem. This might provide the conditions for the formation of new meristem niches. The results provided insight into the JrWOX genes’ general characteristics and their roles in rejuvenation-induced ARF.

Single-Molecule Long-Read Sequencing Reveals the Diversity of Full-Length Transcripts in Leaves of Gnetum (Gnetales)

The limitations of RNA sequencing make it difficult to accurately predict alternative splicing (AS) and alternative polyadenylation (APA) events and long non-coding RNAs (lncRNAs), all of which reveal transcriptomic diversity and the complexity of gene regulation. Gnetum, a genus with ambiguous phylogenetic placement in seed plants, has a distinct stomatal structure and photosynthetic characteristics. In this study, a full-length transcriptome of Gnetum luofuense leaves at different developmental stages was sequenced with the latest PacBio Sequel platform. After correction by short reads generated by Illumina RNA-Seq, 80,496 full-length transcripts were obtained, of which 5269 reads were identified as isoforms of novel genes. Additionally, 1660 lncRNAs and 12,998 AS events were detected. In total, 5647 genes in the G. luofuense leaves had APA featured by at least one poly(A) site. Moreover, 67 and 30 genes from the bHLH gene family, which play an important role in stomatal development and photosynthesis, were identified from the G. luofuense genome and leaf transcripts, respectively. This leaf transcriptome supplements the reference genome of G. luofuense, and the AS events and lncRNAs detected provide valuable resources for future studies of investigating low photosynthetic capacity of Gnetum.

The Transcriptome Landscape of Walnut Interspecies Hybrid (Juglans hindsii × Juglans regia) and Regulation of Cambial Activity in Relation to Grafting

Walnuts (Juglans, Juglandaceae) are known throughout the world as economically important trees that provide fat, protein, vitamins, and minerals as a food source, and produce high-quality timber. We have amended the purpose section to say "However," the omics resources are limited, which hampered the elucidation of molecular mechanisms resulting in their economically important traits (such as yield, fertility alternation, oil synthesis, and wood formation). To enrich the omics database of walnut, there is great need for analyses of its genomic and transcriptomic characteristics. In this study, we reported for the first time of the transcriptome landscape of six important organs or tissues in walnut interspecies hybrid using next-generation sequencing technology. Over 338 million clean reads were obtained. This yielded 74,072 unigenes with an average length of 782.71 bp. To develop an understanding of gene functions and regulatory pathways, 66,355 of the unigenes were identified as homologs of annotated genes and classified into three general categories with 61 functional subcategories. 2,288 out of 2,549 unmapped unigenes had at least one BLAST hit against the public databases. A total of 1,237 transcription factor-encoding genes (TFs) and 2,297 tissue-specific unigenes were identified. Interestingly, in the new shoot between an adult seedling and a grafted tree, the expression of 9,494 unigenes were significantly different, among which 4,388 were up-regulated and 5,106 were down-regulated. Of these, 195, 177, 232, 75, 114, and 68 unigenes were related to transcription factors, cell wall, defense response, transport, plant hormone biosynthesis, and other cambial activity-related functions, respectively. The obtained sequences and putative functional data constitute a resource for future functional analyses in walnut and other woody plants. These findings will be useful in further studies addressing the molecular mechanisms underlying grafting-related cambial activity.

Long noncoding RNAs and miRNAs regulating terpene and tartaric acid biosynthesis in rose-scented geranium

This study aimed to explore the noncoding RNAs, which have emerged as key regulatory molecules in biological processes, in rose-scented geranium. We analyzed RNA-seq data revealing 26 784 long noncoding RNAs (lncRNAs) and 871 miRNAs in rose-scented geranium. A total of 466 lncRNAs were annotated using different plant lncRNA public databases. Furthermore, 372 lncRNAs and 99 miRNAs were detected that target terpene and tartarate biosynthetic pathways. An interactome, comprising of lncRNAs, miRNAs, and mRNAs, was constructed that represents a noncoding RNA regulatory network of the target mRNAs. Real-time quantitative PCR expression validation was done for selected lncRNAs involved in the regulation of terpene and tartaric acid pathways. This study provides the first insights into the regulatory functioning of noncoding RNAs in rose-scented geranium.

Plant lncRNAs are enriched in and move systemically through the phloem in response to phosphate deficiency

In response to phosphate (Pi) deficiency, it has been shown that micro-RNAs (miRNAs) and mRNAs are transported through the phloem for delivery to sink tissues. Growing evidence also indicates that long non-coding RNAs (lncRNAs) are critical regulators of Pi homeostasis in plants. However, whether lncRNAs are present in and move through the phloem, in response to Pi deficiency, remains to be established. Here, using cucumber as a model plant, we show that lncRNAs are enriched in the phloem translocation stream and respond, systemically, to an imposed Pi-stress. A well-known lncRNA, IPS1, the target mimic (TM) of miRNA399, accumulates to a high level in the phloem, but is not responsive to early Pi deficiency. An additional 24 miRNA TMs were also detected in the phloem translocation stream; among them miRNA171 TMs and miR166 TMs were induced in response to an imposed Pi stress. Grafting studies identified 22 lncRNAs which move systemically into developing leaves and root tips. A CU-rich PTB motif was further identified in these mobile lncRNAs. Our findings revealed that lncRNAs respond to Pi deficiency, non-cell-autonomously, and may act as systemic signaling agents to coordinate early Pi deficiency signaling, at the whole-plant level.

Long Non-coding RNAs Coordinate Developmental Transitions and Other Key Biological Processes in Grapevine

Long non-coding RNAs (lncRNAs) are transcripts >200 nucleotides that have prominently surfaced as dynamic regulatory molecules. Using computational approaches, we identified and characterized 56,441 lncRNAs in grapevine (Vitis vinifera) by harnessing RNA-seq data from 10 developmental stages of leaf, inflorescence, and berry tissues. We conducted differential expression analysis and determined tissue- and developmental stage-specificity of lncRNAs in grapevine, which indicated their spatiotemporal regulation. Functional annotation using co-expression analysis revealed their involvement in regulation of developmental transitions in sync with transcription factors (TFs). Further, pathway enrichment analysis revealed lncRNAs associated with biosynthetic and secondary metabolic pathways. Additionally, we identified 115, 560, and 133 lncRNAs as putative miRNA precursors, targets, and endogenous target mimics, respectively, which provided an insight into the interplay of regulatory RNAs. We also explored lncRNA-mediated regulation of extra-chromosomal genes–i.e., mitochondrial and chloroplast coding sequences and observed their involvement in key biological processes like ‘photosynthesis’ and ‘oxidative phosphorylation’. In brief, these transcripts coordinate important biological functions via interactions with both coding and non-coding RNAs as well as TFs in grapevine. Our study would facilitate future experiments in unraveling regulatory mechanisms of development in this fruit crop of economic importance.

Three vital RNA functions and interactions in the process of silk gland apoptosis in silkworm Bombyx mori

The Silkworm Bombyx mori is an important insect in terms of economics and a model organism with a complete metamorphosis. The economic importance of silkworms is dependent on the functions of the silkgland, a specialized organ that synthesizes silk proteins. The silk gland undergoes massive degeneration during the larval to pupal stage, which involves in cell apoptosis. In this paper, high throughput sequencing was used to detect the expression of messenger RNA (mRNA), long noncoding RNA (lncRNA), and microRNA (miRNA) from silk glands of Day 3 in the fifth instar larvae (L5D3) and the spinning 36h (sp36h). We analyzed the Gene Ontology (GO) functions of target genes of the differentially expressed lncRNAs and miRNAs. We investigated the regulations of mRNA, lncRNA, and miRNA on silk gland apoptosis in L5D3 and sp36h. In total, 10,947 lncRNAs were detected in the silk gland and the index number TCONS-00021360 lncRNA may be involved in the process of apoptosis. In addition, 344 miRNAs targeted 285 mRNAs were related to the death process under GO entry. The results indicated that miRNAs play an important role in the molecular regulation of the silk gland apoptosis compared with that of lncRNAs. Finally, we screened 746 lncRNAs and 20 miRNAs that might interact with BmDredd, and drew an interaction network among them.

Transcriptional regulation of osmotic stress tolerance in wheat (Triticum aestivum L.)

The current review provides an updated, new insights into the regulation of transcription mediated underlying mechanisms of wheat plants to osmotic stress perturbations. Osmotic stress tolerance mechanisms being complex are governed by multiple factors at physiological, biochemical and at the molecular level, hence approaches like "OMICS" that can underpin mechanisms behind osmotic tolerance in wheat is of paramount importance. The transcription factors (TFs) are a class of molecular proteins, which are involved in regulation, modulation and orchestrating the responses of plants to a variety of environmental stresses. Recent reports have provided novel insights on the role of TFs in osmotic stress tolerance via direct molecular links. However, our knowledge on the regulatory role TFs during osmotic stress tolerance in wheat remains limited. The present review in its first part sheds light on the importance of studying the role of osmotic stress tolerance in wheat plants and second aims to decipher molecular mechanisms of TFs belonging to several classes, including DREB, NAC, MYB, WRKY and bHLH, which have been reported to engage in osmotic stress mediated gene expression in wheat and third part covers the systems biology approaches to understand the transcriptional regulation of osmotic stress and the role of long non-coding RNAs in response to osmotic stress with special emphasis on wheat. The current concept may lead to an understanding in molecular regulation and signalling interaction of TFs under osmotic stress to clarify challenges and problems for devising potential strategies to improve complex regulatory events involved in plant tolerance to osmotic stress adaptive pathways in wheat.

Genome-wide identification, characterization and expression analysis of long non-coding RNAs in different tissues of apple

Recently, the long non-coding RNAs (lncRNAs), which play important roles in various complex biological processes, have received more attention in plants. However, little information is available on lncRNAs in woody fruit trees and their potential regulatory roles remain poorly understood, especially in apple. Here, a total of 1726 high-confidence lncRNAs (hc-lncRNAs) were identified in different tissues including young fruits, shoot tips, stem phloem and root tips using high-throughput sequencing. These lncRNAs are distributed across all 17 apple chromosomes, and >85% come from intergenic regions. The apple lncRNAs have longer transcript lengths and greater exon numbers than protein-coding genes. Additionally, among the 1726 hc-lncRNA, 850 are predicted to have target genes. These target genes are involved in many processes including hormone signaling, sugar metabolism, and the cell cycle and stress responses. Furthermore, 57, 74, 168 and 78 lncRNAs specifically expressed in root tips, shoot tips, young fruits and stem phloem were analyzed using the COG (cluster of orthologous group)and GO (gene ontology) databases. Young fruits contain the most unique lncRNAs, which are involved in biological processes such as energy production and conversion, carbohydrate transport and metabolism, posttranslational modification and protein turnover. Quantitative real-time PCR (qRT-PCR) is employed to confirm the different expression levels among the four tissues. Moreover, the expression levels of eight fruit-related lncRNAs are investigated during different fruit development stages, which indicates they have important roles in fruit ripening and sugar metabolism. Overall, our genome-wide research on lncRNAs in different apple tissues provides valuable clues and information that can help elucidate the potential roles of lncRNAs in the growth and development of apple, as well as in other fruit trees.

Identification and validation of reference genes for qRT-PCR analysis in mulberry (Morus alba L.)

Mulberry (Morus alba L.) is an important economic tree species in many countries. Quantitative real time PCR (qRT-PCR) has become a widely used method for gene expression studies in plants. A suitable reference gene is essential to ensure accurate and reliable results for qRT-PCR analyses. However, no reports describing the selection of reference genes have been published for mulberry. In this work, we evaluated the stability of twenty candidate reference genes in different plant tissues and under different stress conditions by qRT-PCR in mulberry using algorithms in two programs—geNorm and NormFinder. The results revealed that TUB2, UBI4, ACTIN3 and RPL4 were ranked as the most stable reference genes in the samples subsets, whereas EF1α4 and TUB3showed the least stability with both algorithms. To further validate the stability of the reference genes, the expression patterns of six genes of mulberry were analyzed by normalization with the selected reference genes. Our study will benefit future analyses of gene expression in mulberry.

Genome-wide analysis of long non-coding RNAs in Pichia pastoris during stress by RNA sequencing

Long non-coding RNAs play significant roles in many biological processes. The roles of lncRNAs in Pichia pastoris remain unclear. In this work, we focused on the identification of lncRNAs in P. pastoris and exploration of their potential roles in stress response to PLA₂ overexpression and methanol induction. By strand specific RNA sequencing, 208 novel long non-coding RNAs were identified and analyzed. Bioinformatic analysis showed potential trans-target genes and cis-regulated genes of 39 differential lncRNAs. Functional annotation and sequence motif analysis indicated that lncRNAs participate in pathways related to methanol degradation and production of the recombinant protein. The differential expression of lncRNAs was validated by qRT-PCR. Lastly, the potential functions of three lncRNAs were evaluated by knockdown of their expression and analysis of the expression levels of target genes. Our study identifies novel lncRNAs in P. pastoris induced during use as a bioreactor, facilitating future functional research.

Non-coding RNAs and plant male sterility: current knowledge and future prospects

Latest outcomes assign functional role to non-coding (nc) RNA molecules in regulatory networks that confer male sterility to plants. Male sterility in plants offers great opportunity for improving crop performance through application of hybrid technology. In this respect, cytoplasmic male sterility (CMS) and sterility induced by photoperiod (PGMS)/temperature (TGMS) have greatly facilitated development of high-yielding hybrids in crops. Participation of non-coding (nc) RNA molecules in plant reproductive development is increasingly becoming evident. Recent breakthroughs in rice definitively associate ncRNAs with PGMS and TGMS. In case of CMS, the exact mechanism through which the mitochondrial ORFs exert influence on the development of male gametophyte remains obscure in several crops. High-throughput sequencing has enabled genome-wide discovery and validation of these regulatory molecules and their target genes, describing their potential roles performed in relation to CMS. Discovery of ncRNA localized in plant mtDNA with its possible implication in CMS induction is intriguing in this respect. Still, conclusive evidences linking ncRNA with CMS phenotypes are currently unavailable, demanding complementing genetic approaches like transgenics to substantiate the preliminary findings. Here, we review the recent literature on the contribution of ncRNAs in conferring male sterility to plants, with an emphasis on microRNAs. Also, we present a perspective on improved understanding about ncRNA-mediated regulatory pathways that control male sterility in plants. A refined understanding of plant male sterility would strengthen crop hybrid industry to deliver hybrids with improved performance.

A Novel LncRNA, MuLnc1, Associated With Environmental Stress in Mulberry (Morus multicaulis)

Environmental stresses are major constraints that limit the leaf productivity and quality of mulberry. LncRNAs have emerged as important regulators in response to biotic and abiotic stresses in plants. However, the functions and mechanisms of most lncRNAs remain largely unknown. A novel lncRNA designated as MuLnc1 was found to be cleaved by mul-miR3954 and produce secondary siRNAs in a 21 nt phase in mulberry. It was demonstrated that one of the siRNAs produced, si161579, can silence the expression of the calmodulin-like protein gene CML27 of mulberry (MuCML27). When MuCML27 was heterologously expressed in Arabidopsis, the transgenic plants exhibited enhanced resistance to Botrytis cinerea and Pseudomonas syringae pv tomato DC3000. In addition, the transgenic MuCML27-overexpressing Arabidopsis plants are more tolerant to salt and drought stresses. Furthermore, the network of mul-miR3954-MuLnc1-siRNAs-mRNAs was modeled to elucidate the interaction between lncRNAs and sRNAs with mRNAs. All of these, taken together, suggest that MuLnc1 was associated with environmental stress in mulberry and may be considered as a potential genetic improvement target gene of mulberry. The information provided may shed light on the complicated gene expression regulatory mechanisms in mulberry stress responses.

Present Scenario of Long Non-Coding RNAs in Plants

Small non-coding RNAs have been extensively studied in plants over the last decade. In contrast, genome-wide identification of plant long non-coding RNAs (lncRNAs) has recently gained momentum. LncRNAs are now being recognized as important players in gene regulation, and their potent regulatory roles are being studied comprehensively in eukaryotes. LncRNAs were first reported in humans in 1992. Since then, research in animals, particularly in humans, has rapidly progressed, and a vast amount of data has been generated, collected, and organized using computational approaches. Additionally, numerous studies have been conducted to understand the roles of these long RNA species in several diseases. However, the status of lncRNA investigation in plants lags behind that in animals (especially humans). Efforts are being made in this direction using computational tools and high-throughput sequencing technologies, such as the lncRNA microarray technique, RNA-sequencing (RNA-seq), RNA capture sequencing, (RNA CaptureSeq), etc. Given the current scenario, significant amounts of data have been produced regarding plant lncRNAs, and this amount is likely to increase in the subsequent years. In this review we have documented brief information about lncRNAs and their status of research in plants, along with the plant-specific resources/databases for information retrieval on lncRNAs.

Genome-wide analysis of long non-coding RNAs at the mature stage of sea buckthorn (Hippophae rhamnoides Linn) fruit

Long non-coding RNAs (lncRNAs), which are >200nt longer transcripts, potentially play important roles in almost all biological processes in plants and mammals. However, the functions and profiles of lncRNAs in fruit is less understood. Therefore, it is urgent and necessary to identify and analyze the functions of lncRNAs in sea buckthorns. Using RNA-sequencing, we synthetically identified lncRNAs in mature fruit from the red and yellow sea buckthorn. We obtained 567,778,938 clean reads from six samples and identified 3428 lncRNAs in mature fruit, including 2498 intergenic lncRNAs, 593 anti-sense lncRNAs, and 337 intronic lncRNAs. We also identified 3819 and 2295 circular RNAs in red and yellow sea buckthorn Fruit. In the aspects of gene architecture and expression, our results showed significant differences among the three lncRNA subtypes. We also investigated the effect of lncRNAs on its cis and trans target genes. Based on target genes analysis, we obtained 61 different expression lncRNAs (DE-lncRNAs) between these two sea buckthorns, including 23 special expression lncRNAs in red fruit and 22 special expression lncRNAs in yellow fruit. Importantly, we found a few DE-lncRNAs play cis and trans roles for genes in the Carotenoid biosynthesis, ascorbate and aldarate metabolism and fatty acid metabolism pathways. Our study provides a resource for lncRNA studies in mature fruit. It probably encourages researchers to deeply study fruit-coloring. It expands our knowledge about lncRNA biology and the annotation of the sea buckthorn genome.