Identification and characterization of ncRNA-associated ceRNA networks in Arabidopsis leaf development

ETENLNC: An end to end lncRNA identification and analysis framework to facilitate construction of known and novel lncRNA regulatory networks

Long non-coding RNAs (lncRNAs) play crucial roles in the regulation of gene expression and maintenance of genomic integrity through various interactions with DNA, RNA, and proteins. The availability of large-scale sequence data from various high-throughput platforms has opened possibilities to identify, predict, and functionally annotate lncRNAs. As a result, there is a growing demand for an integrative computational framework capable of identifying known lncRNAs, predicting novel lncRNAs, and inferring the downstream regulatory interactions of lncRNAs at the genome-scale. We present ETENLNC (End-To-End-Novel-Long-NonCoding), a user-friendly, integrative, open-source, scalable, and modular computational framework for identifying and analyzing lncRNAs from raw RNA-Seq data. ETENLNC employs six stringent filtration steps to identify novel lncRNAs, performs differential expression analysis of mRNA and lncRNA transcripts, and predicts regulatory interactions between lncRNAs, mRNAs, miRNAs, and proteins. We benchmarked ETENLNC against six existing tools and optimized it for desktop workstations and high-performance computing environments using data from three different species. ETENLNC is freely available on GitHub: https://github.com/EvolOMICS-TU/ETENLNC.

Source leaves are regulated by sink strengths through non-coding RNAs and alternative polyadenylation in cucumber (Cucumis sativus L.)

BACKGROUND

The yield of major crops is generally limited by sink capacity and source strength. Cucumber is a typical raffinose family oligosaccharides (RFOs)-transporting crop. Non-coding RNAs and alternative polyadenylation (APA) play important roles in the regulation of growth process in plants. However, their roles on the sink‒source regulation have not been demonstrated in RFOs-translocating species.

RESULTS

Here, whole-transcriptome sequencing was applied to compare the leaves of cucumber under different sink strength, that is, no fruit-carrying leaves (NFNLs) and fruit-carrying leaves (FNLs) at 12th node from the bottom. The results show that 1101 differentially expressed (DE) mRNAs, 79 DE long non-coding RNAs (lncRNAs) and 23 DE miRNAs were identified, which were enriched in photosynthesis, energy production and conversion, plant hormone signal transduction, starch and carbohydrate metabolism and protein synthesis pathways. Potential co-expression networks like, DE lncRNAs-DE mRNAs/ DE miRNAs-DE mRNAs, and competing endogenous RNA (ceRNA) regulation models (DE lncRNAs-DE miRNAs-DE mRNAs) associated with sink‒source allocation, were constructed. Furthermore, 37 and 48 DE genes, which enriched in MAPK signaling and plant hormone signal transduction pathway, exist differentially APA, and SPS (CsaV3_2G033300), GBSS1 (CsaV3_5G001560), ERS1 (CsaV3_7G029600), PNO1 (CsaV3_3G003950) and Myb (CsaV3_3G022290) may be regulated by both ncRNAs and APA between FNLs and NFNLs, speculating that ncRNAs and APA are involved in the regulation of gene expression of cucumber sink‒source carbon partitioning.

CONCLUSIONS

These results reveal a comprehensive network among mRNAs, ncRNAs, and APA in cucumber sink-source relationships. Our findings also provide valuable information for further research on the molecular mechanism of ncRNA and APA to enhance cucumber yield.

LncRNA cis- and trans-regulation provides new insight into drought stress responses in wild barley

Drought is one of the most common abiotic stresses that affect barley productivity. Long noncoding RNA (lncRNA) has been reported to be widely involved in abiotic stress, however, its function in the drought stress response in wild barley remains unclear. In this study, RNA sequencing was performed to identify differentially expressed lncRNAs (DElncRNA) among two wild barley and two cultivated barley genotypes. Then, the cis-regulatory networks were according to the chromosome position and the expression level correction. The GO annotation indicates that these cis-target genes are mainly involved in "ion transport transporter activity" and "metal ion transport transporter activity". Through weighted gene co-expression network analysis (WGCNA), 10 drought-related modules were identified to contract trans-regulatory networks. The KEGG annotation demonstrated that these trans-target genes were enriched for photosynthetic physiology, brassinosteroid biosynthesis, and flavonoid metabolism. In addition, we constructed the lncRNA-mediated ceRNA regulatory network by predicting the microRNA response elements (MREs). Furthermore, the expressions of lncRNAs were verified by RT-qPCR. Functional verification of a candidate lncRNA, MSTRG.32128, demonstrated its positive role in drought response and root growth and development regulation. Hormone content analysis provided insights into the regulatory mechanisms of MSTRG.32128 in root development, revealing its involvement in auxin and ethylene signal transduction pathways. These findings advance our understanding of lncRNA-mediated regulatory mechanisms in barley under drought stress. Our results will provide new insights into the functions of lncRNAs in barley responding to drought stress.

CircRNA: An emerging star in plant research: A review

CircRNAs are a class of covalently closed non-coding RNA formed by linking the 5' terminus and the 3' terminus after reverse splicing. CircRNAs are widely found in eukaryotes, and they are highly conserved, with spatio-temporal expression specificity and stability. CircRNAs can act as miRNA sponges to regulate the expression of downstream target genes, regulating the transcription of parental genes and some can even be translated into peptides or proteins. Research on circRNAs in plants is still in its infancy compared to that in animals. With the deepening of research, the results of a variety of plant circRNAs suggest that they play an important role in growth and development, and tolerance towards abiotic stresses such as salt, drought, low temperature, high temperature and other adverse environments. In this review paper, we elaborated the molecular characteristics, mechanism of action, function and bioinformatics databases of plant circRNAs, combined with the progress of circRNA research in animals, discussed the potential mechanism of action of plant circRNAs, and proposed the unsolved problems and prospects for future application of plant circRNAs.

Identification and Characterisation of the CircRNAs Involved in the Regulation of Leaf Colour in Quercus mongolica

Circular RNAs (circRNAs) are important regulatory molecules involved in various biological processes. However, the potential function of circRNAs in the turning red process of Quercus mongolica leaves is unclear. This study used RNA-seq data to identify 6228 circRNAs in leaf samples from four different developmental stages and showed that 88 circRNAs were differentially expressed. A correlation analysis was performed between anthocyanins and the circRNAs. A total of 16 circRNAs that may be involved in regulating the colour of Mongolian oak leaves were identified. CircRNAs may affect the colour of Q. mongolica leaves by regulating auxin, cytokinin, gibberellin, ethylene, and abscisic acid. This study revealed the potential role of circRNAs in the colour change of Q. mongolica leaves.

The Construction of lncRNA/circRNA-miRNA-mRNA Networks Reveals Functional Genes Related to Growth Traits in Schima superba

Schima superba is a precious timber and fire-resistant tree species widely distributed in southern China. Currently, there is little knowledge related to its growth traits, especially with respect to molecular breeding. The lack of relevant information has delayed the development of modern breeding. The purpose is to identify probable functional genes involved in S. superba growth through whole transcriptome sequencing. In this study, a total of 32,711 mRNAs, 525 miRNAs, 54,312 lncRNAs, and 1522 circRNAs were identified from 10 S. superba individuals containing different volumes of wood. Four possible regulators, comprising three lncRNAs, one circRNA, and eleven key miRNAs, were identified from the regulatory networks of lncRNA-miRNA-mRNA and circRNA-miRNA-mRNA to supply information on ncRNAs. Several candidate genes involved in phenylpropane and cellulose biosynthesis pathways, including Ss4CL2, SsCSL1, and SsCSL2, and transcription factors, including SsDELLA2 (SsSLR), SsDELLA3 (SsSLN), SsDELLA5 (SsGAI-like2), and SsNAM1, were identified to reveal the molecular regulatory mechanisms regulating the growth traits of S. superba. The results not merely provide candidate functional genes related to S. superba growth trait and will be useful to carry out molecular breeding, but the strategy and method also provide scientists with an effective approach to revealing mechanisms behind important economic traits in other species.

An accretive detection method for in silico identification and validation of circular RNAs in wheat (Triticum aestivum L.) using RT-qPCR

BACKGROUND

Circular RNAs (circRNAs) are novel class of non-coding RNAs, which are involved in various functions at the transcriptional and post-transcriptional level in response to a fungal pathogen (Puccinia triticina), including microRNA (miRNA) sponge, RNA binding proteins sponge, regulation of parental gene and biomarkers. Detailed analysis of wheat circRNAs is essential to accelerate the regulated expression of fungal miRNAs. Therefore, we suggest a protocol to aid circRNA identification through RNA-Seq data using various algorithms based on perl script followed by validation through divergent primer designing, standard PCR, and RT-qPCR assays.

METHODS AND RESULT

The divergent primer has been widely used to detect, validate, and quantify back-spliced junction (BSJ) of circRNAs. The procedure covers index file formation, circRNA identification and BSJ detections. However, the laboratory validation of circRNA includes wheat genomic DNA isolation, RNA isolation and its cDNA conversion upto validation. In this study, we identified 28 circRNAs from RNA-Seq of S0 and R0, wherein six circRNAs are commonly present and 75% of the identified circRNAs were belongs to inter-genic, 14% were exonic and intronic category were 11%. Divergent primer designing method successfully validated the two circRNAs via RT-qPCR assay, where circRNA_2 showed less relative expression pattern than circRNA_1 in contrast with housekeeping genes.

CONCLUSION

Thus, our results of identified and validated circRNAs showed that, this protocol is quite helpful, relatively easy, reliable, and accurate for large datasets as other algorithms need various dependencies and have complex scripts with high chances of error occurrence. Additionally, analysis time will vary depending on the expertise level and the number of RNA-Seq data. This proposed protocol can also be used for a wide range of monocotyledons belonging to the Poaceae plant family.

Full-length transcriptome characterization of Platycladus orientalis based on the PacBio platform

As a unique and native conifer in China, Platycladus orientalis is widely used in soil erosion control, garden landscapes, timber, and traditional Chinese medicine. However, due to the lack of reference genome and transcriptome, it is limited to the further molecular mechanism research and gene function mining. To develop a full-length reference transcriptome, tissues from five different parts of P. orientalis and four cone developmental stages were sequenced and analyzed by single-molecule real-time (SMRT) sequencing through the PacBio platform in this study. Overall, 37,111 isoforms were detected by PacBio with an N50 length of 2,317 nt, an average length of 1,999 bp, and the GC content of 41.81%. Meanwhile, 36,120 coding sequences, 5,645 simple sequence repeats (SSRs), 1,201 non-coding RNAs (lncRNAs), and 182 alternative splicing (AS) events with five types were identified using the results obtained from the PacBio transcript isoforms. Furthermore, 1,659 transcription factors (TFs) were detected and belonged to 51 TF families. A total of 35,689 transcripts (96.17%) were annotated through the NCBI nr, KOG, Swiss-Prot and KEGG databases, and 385 transcript isoforms related to 8 types of hormones were identified incorporated into plant hormone signal transduction pathways. The assembly and revelation of the full-length transcriptome of P. orientalis offer a pioneering insight for future investigations into gene function and genetic breeding within Platycladus species.

The ncRNAs Involved in the Regulation of Abiotic Stress-Induced Anthocyanin Biosynthesis in Plants

Plants have evolved complicated defense and adaptive systems to grow in various abiotic stress environments such as drought, cold, and salinity. Anthocyanins belong to the secondary metabolites of flavonoids with strong antioxidant activity in response to various abiotic stress and enhance stress tolerance. Anthocyanin accumulation often accompanies the resistance to abiotic stress in plants to scavenge reactive oxygen species (ROS). Recent research evidence showed that many regulatory pathways such as osmoregulation, antioxidant response, plant hormone response, photosynthesis, and respiration regulation are involved in plant adaption to stress. However, the molecular regulatory mechanisms involved in controlling anthocyanin biosynthesis in relation to abiotic stress response have remained obscure. Here, we summarize the current research progress of specific regulators including small RNAs, and lncRNAs involved in the molecular regulation of abiotic stress-induced anthocyanin biosynthesis. In addition, an integrated regulatory network of anthocyanin biosynthesis controlled by microRNAs (miRNAs), long non-coding RNAs (lncRNAs), transcription factors, and stress response factors is also discussed. Understanding molecular mechanisms of anthocyanin biosynthesis for ROS scavenging in various abiotic stress responses will benefit us for resistance breeding in crop plants.

Plant long non-coding RNAs: identification and analysis to unveil their physiological functions

Eukaryotic genomes encode thousands of RNA molecules; however, only a minimal fraction is translated into proteins. Among the non-coding elements, long non-coding RNAs (lncRNAs) play important roles in diverse biological processes. LncRNAs are associated mainly with the regulation of the expression of the genome; nonetheless, their study has just scratched the surface. This is somewhat due to the lack of widespread conservation at the sequence level, in addition to their relatively low and highly tissue-specific expression patterns, which makes their exploration challenging, especially in plant genomes where only a few of these molecules have been described completely. Recently published high-quality genomes of crop plants, along with new computational tools, are considered promising resources for studying these molecules in plants. This review briefly summarizes the characteristics of plant lncRNAs, their presence and conservation, the different protocols to find these elements, and the limitations of these protocols. Likewise, it describes their roles in different plant physiological phenomena. We believe that the study of lncRNAs can help to design strategies to reduce the negative effect of biotic and abiotic stresses on the yield of crop plants and, in the future, help create fruits and vegetables with improved nutritional content, higher amounts of compounds with positive effects on human health, better organoleptic characteristics, and fruits with a longer postharvest shelf life.

Spatiotemporal miRNA and transcriptomic network dynamically regulate the developmental and senescence processes of poplar leaves

Poplar is an important afforestation and urban greening species. Poplar leaf development occurs in stages, from young to mature and then from mature to senescent; these are accompanied by various phenotypic and physiological changes. However, the associated transcriptional regulatory network is relatively unexplored. We first used principal component analysis to classify poplar leaves at different leaf positions into two stages: developmental maturity (the stage of maximum photosynthetic capacity); and the stage when photosynthetic capacity started to decline and gradually changed to senescence. The two stages were then further subdivided into five intervals by gene expression clustering analysis: young leaves, the period of cell genesis and functional differentiation (L1); young leaves, the period of development and initial formation of photosynthetic capacity (L3-L7); the period of maximum photosynthetic capacity of functional leaves (L9-L13); the period of decreasing photosynthetic capacity of functional leaves (L15-L27); and the period of senescent leaves (L29). Using a weighted co-expression gene network analysis of regulatory genes, high-resolution spatiotemporal transcriptional regulatory networks were constructed to reveal the core regulators that regulate leaf development. Spatiotemporal transcriptome data of poplar leaves revealed dynamic changes in genes and miRNAs during leaf development and identified several core regulators of leaf development, such as GRF5 and MYB5. This in-depth analysis of transcriptional regulation during leaf development provides a theoretical basis for exploring the biological basis of the transcriptional regulation of leaf development and the molecular design of breeding for delaying leaf senescence.

Comprehensive insights into the regulatory mechanisms of lncRNA in alkaline-salt stress tolerance in rice

BACKGROUND

Alkaline-salt is one of the abiotic stresses that slows plant growth and developmental processes and threatens crop yield. Long non-coding RNAs (lncRNAs) are endogenous RNA found in plants that engage in a variety of cellular functions and stress responses.

METHOD

lncRNAs act as competing endogenous RNAs (ceRNA) and constitute a new set of gene control. The precise regulatory mechanism by which lncRNAs function as ceRNAs in response to alkaline-salt stress remains unclear. We identified alkaline-salt responsive lncRNAs using transcriptome-wide analysis of two varieties including alkaline-salt tolerant [WD20342 (WD)] and alkaline-salt sensitive [Caidao (CD)] rice cultivar under control and alkaline-salt stress treated [WD20342 (WDT, and Caidao (CDT)] conditions.

RESULTS

Investigating the competitive relationships between mRNAs and lncRNAs, we next built a ceRNA network involving lncRNAs based on the ceRNA hypothesis. Expression profiles revealed that a total of 65, 34, and 1549 differentially expressed (DE) lncRNAs, miRNAs, and mRNAs were identified in alkaline-salt tolerant WD (Control) vs. WDT (Treated). Similarly, 75 DE-lncRNAs, 34 DE-miRNAs, and 1725 DE-mRNAs (including up-regulated and down-regulated) were identified in alkaline-salt sensitive CD (Control) vs. CDT (Treated), respectively. An alkaline-salt stress ceRNA network discovered 321 lncRNA-miRNA-mRNA triplets in CD and CDT, with 32 lncRNAs, 121 miRNAs, and 111 mRNAs. Likewise, 217 lncRNA-miRNA-mRNA triplets in WD and WDT revealed the NONOSAT000455-osa_miR5809b-LOC_Os11g01210 triplet with the highest degree as a hub node with the most significant positive correlation in alkaline-salt stress response.

CONCLUSION

The results of our investigation indicate that osa-miR5809b is dysregulated and plays a part in regulating the defense response of rice against alkaline-salt stress. Our study highlights the regulatory functions of lncRNAs acting as ceRNAs in the mechanisms underlying alkaline-salt resistance in rice.

Role of non-coding RNAs against salinity stress in Oryza species: Strategies and challenges in analyzing miRNAs, tRFs and circRNAs

Salinity is an imbalanced concentration of mineral salts in the soil or water that causes yield loss in salt-sensitive crops. Rice plant is vulnerable to soil salinity stress at seedling and reproductive stages. Different non-coding RNAs (ncRNAs) post-transcriptionally regulate different sets of genes during different developmental stages under varying salinity tolerance levels. While microRNAs (miRNAs) are well known small endogenous ncRNAs, tRNA-derived RNA fragments (tRFs) are an emerging class of small ncRNAs derived from tRNA genes with a demonstrated regulatory role, like miRNAs, in humans but unexplored in plants. Circular RNA (circRNA), another ncRNA produced by back-splicing events, acts as target mimics by preventing miRNAs from binding with their target mRNAs, thereby reducing the miRNA's action upon its target. Same may hold true between circRNAs and tRFs. Hence, the work done on these ncRNAs was reviewed and no reports were found for circRNAs and tRFs under salinity stress in rice, either at seedling or reproductive stages. Even the reports on miRNAs are restricted to seedling stage only, in spite of severe effects on rice crop production due to salt stress during reproductive stage. Moreover, this review sheds light on strategies to predict and analyze these ncRNAs in an effective manner.

Roles of long non-coding RNAs in plant immunity

Robust plant immune systems are fine-tuned by both protein-coding genes and non-coding RNAs. Long non-coding RNAs (lncRNAs) refer to RNAs with a length of more than 200 nt and usually do not have protein-coding function and do not belong to any other well-known non-coding RNA types. The non-protein-coding, low expression, and non-conservative characteristics of lncRNAs restrict their recognition. Although studies of lncRNAs in plants are in the early stage, emerging studies have shown that plants employ lncRNAs to regulate plant immunity. Moreover, in response to stresses, numerous lncRNAs are differentially expressed, which manifests the actions of low-expressed lncRNAs and makes plant-microbe/insect interactions a convenient system to study the functions of lncRNAs. Here, we summarize the current advances in plant lncRNAs, discuss their regulatory effects in different stages of plant immunity, and highlight their roles in diverse plant-microbe/insect interactions. These insights will not only strengthen our understanding of the roles and actions of lncRNAs in plant-microbe/insect interactions but also provide novel insight into plant immune responses and a basis for further research in this field.

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.

Identification and characterization of CircRNA-associated CeRNA networks in moso bamboo under nitrogen stress

BACKGROUND

Nitrogen is a macronutrient element for plant growth and development. Circular RNAs (circRNAs) serve as pivotal regulators for the coordination between nutrient supply and plant demand. Moso bamboo (Phyllostachys edulis) is an excellent plant with fast growth, and the mechanism of the circRNA-target module in response to nitrogen remains unclear.

RESULTS

Deep small RNA sequencing results of moso bamboo seedlings under different concentrations of KNO₃ (N0 = 0 mM, N6 = 6 mM, N18 = 18 mM) were used to identify circRNAs. A total of 549 circRNAs were obtained, of which 309 were generated from corresponding parental coding genes including 66 new ones. A total of 536 circRNA-parent genes were unevenly distributed in 24 scaffolds and were associated with root growth and development. Furthermore, 52 differentially expressed circRNAs (DECs) were obtained, including 24, 33 and 15 DECs from three comparisons of N0 vs. N6, N0 vs. N18 and N6 vs. N18, respectively. Based on integrative analyses of the identified DECs, differentially expressed mRNAs (DEGs), and miRNAs (DEMs), a competitive endogenous RNA (ceRNA) network was constructed, including five DECs, eight DEMs and 32 DEGs. A regulatory module of PeSca_6:12,316,320|12,372,905-novel_miR156-PH02Gene35622 was further verified by qPCR and dual-luciferase reporter assays.

CONCLUSION

The results indicated that circRNAs could participate in multiple biological processes as miRNA sponges, including organ nitrogen compound biosynthesis and metabolic process regulation in moso bamboo. Our results provide valuable information for further study of circRNAs in moso bamboo under fluctuating nitrogen conditions.

Identification, biogenesis, function, and mechanism of action of circular RNAs in plants

Circular RNAs (circRNAs) are a class of single-stranded, closed RNA molecules with unique functions that are ubiquitously expressed in all eukaryotes. The biogenesis of circRNAs is regulated by specific cis-acting elements and trans-acting factors in humans and animals. circRNAs mainly exert their biological functions by acting as microRNA sponges, forming R-loops, interacting with RNA-binding proteins, or being translated into polypeptides or proteins in human and animal cells. Genome-wide identification of circRNAs has been performed in multiple plant species, and the results suggest that circRNAs are abundant and ubiquitously expressed in plants. There is emerging compelling evidence to suggest that circRNAs play essential roles during plant growth and development as well as in the responses to biotic and abiotic stress. However, compared with recent advances in human and animal systems, the roles of most circRNAs in plants are unclear at present. Here we review the identification, biogenesis, function, and mechanism of action of plant circRNAs, which will provide a fundamental understanding of the characteristics and complexity of circRNAs in plants.

LncRNA-mediated ceRNA regulatory network provides new insight into chlorogenic acid synthesis in sweet potato

Sweet potato (Ipomoea batatas L.) is considered a highly nutritional and economical crop due to its high contents of bioactive substances, such as anthocyanin and chlorogenic acid (CGA), especially in leaves and stems. The roles of noncoding RNAs (ncRNA), including long noncoding RNA (lncRNA) and microRNA (miRNA), in CGA synthesis, are still unknown. In this study, the differentially expressed (DE) mRNAs, miRNAs, and lncRNAs in two leafy vegetable genotypes "FS7-6-14-7" (high CGA content) and "FS7-6" (low CGA content) were identified. The cis-regulation between lncRNA and mRNA was analyzed. Then, the CGA synthesis-related modules MEBlue and MEYellow were identified to detect trans-regulation mRNA-lncRNA pairs. The GO and KEGG annotations suggested that mRNA in these two modules was significantly enriched in the secondary metabolite synthesis biosynthesis category. A competing endogenous RNAs (ceRNA) network, including 8730 miRNA-mRNA and 444 miRNA-lncRNA pairs, was constructed by DEmiRNA target prediction. Then, a CGA synthesis-related ceRNA network was obtained with lncRNA and mRNA from MEBlue and MEYellow. Finally, one relational pair, MSTRG.47662.1/mes-miR398/itb04g00990, was selected for functional validation. Overexpression of lncRNA MSTRG.47662.1 and mRNA itb04g00990 increased CGA content in both tobacco and sweet potato callus, while overexpression of miRNA mes-miR398 decreased CGA content. Meanwhile, regression analysis of the expression patterns demonstrated that MSTRG.47662.1, acting as a ceRNA, promoted itb04g00990 expression by competitively binding mes-miR398 in CGA synthesis in sweet potato. Our results provide insights into how ncRNA-mediated ceRNA regulatory networks likely contribute to CGA synthesis in leafy sweet potato.

Genome-wide circular RNA profiling and competing endogenous RNA regulatory network analysis provide new insights into the molecular mechanisms underlying early somatic embryogenesis in Dimocarpus longan Lour

Circular RNAs (circRNAs) are widely involved in plant growth and development. However, the function of circRNAs in plant somatic embryogenesis (SE) remains elusive. Here, by using high-throughput sequencing, a total of 5029 circRNAs were identified in the three stages of longan (Dimocarpus longan Lour.) early SE. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that differentially expressed (DE) circRNA host genes were enriched in the 'non-homologous end-joining' (NHEJ) and 'butanoate metabolism' pathways. In addition, the reactive oxygen species (ROS) content during longan early SE was determined. The results indicated that ROS-induced DNA double-strand breaks may not depend on the NHEJ repair pathway. Correlation analyses of the levels of related metabolites (glutamate, γ-aminobutyrate and pyruvate) and the expression levels of circRNAs and their host genes involved in butanoate metabolism were performed. The results suggested that circRNAs may act as regulators of the expression of cognate mRNAs, thereby affecting the accumulation of related compounds. A competing endogenous RNA (ceRNA) network of DE circRNAs, DE mRNAs, DE long noncoding RNAs (lncRNAs) and DE microRNAs (miRNAs) was constructed. The results showed that the putative targets of the noncoding RNA (ncRNAs) were significantly enriched in the KEGG pathways 'mitogen-activated protein kinase signaling' and 'nitrogen metabolism'. Furthermore, the expression patterns of the candidate circRNAs, lncRNAs, miRNAs and mRNAs confirmed the negative correlation between miRNAs and ceRNAs. In addition, two circRNA overexpression vectors were constructed to further verify the ceRNA network correlations in longan early SE. Our study revealed the potential role of circRNAs in longan early SE, providing new insights into the intricate regulatory mechanism underlying plant SE.

Integrative Investigation of Root-Related mRNAs, lncRNAs and circRNAs of “Muscat Hamburg” (Vitis vinifera L.) Grapevine in Response to Root Restriction through Transcriptomic Analyses

Root restriction is a physical and ecological cultivation mode which restricts plant roots into a limited container to regulate vegetative and reproduction growth by reshaping root architecture. However, little is known about related molecular mechanisms. To uncover the root-related regulatory network of endogenous RNAs under root restriction cultivation (referred to RR), transcriptome-wide analyses of mRNAs, long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs) involved in root development were performed. During root development, RR treatment had a positive effect on root weight, typically, young roots were significantly higher than conventional cultivation (referred to NR) treatment, suggesting that root architecture reconstruction under RR was attributed to the vigorous induction into lateral roots. Furthermore, a total of 26,588 mRNAs, 1971 lncRNAs, and 2615 circRNAs were identified in root of annual “Muscat Hamburg” grapevine by the transcriptomic analyses. The expression profile of mRNAs, lncRNAs and circRNA were further confirmed by the quantitative real-time PCR (RT-qPCR). Gene ontology enrichment analysis showed that a majority of the differentially expressed mRNAs, lncRNAs and circRNAs were enriched into the categories of cellular process, metabolic process, cell part, binding, and catalytic activity. In addition, the regulatory network of endogenous RNAs was then constructed by the prediction of lncRNA-miRNA-mRNA and circRNA-miRNA-mRNA network, implying that these RNAs play significant regulatory roles for root architecture shaping in response to root restriction. Our results, for the first time, the regulatory network of competitive endogenous RNAs (ceRNAs) functions of lncRNA and circRNA was integrated, and a basis for studying the potential functions of non-coding RNAs (ncRNAs) during root development of grapevine was provided.

Whole-Transcriptome Analysis Reveals Autophagy Is Involved in Early Senescence of zj-es Mutant Rice

Senescence is a necessary stage of plant growth and development, and the early senescence of rice will lead to yield reduction and quality decline. However, the mechanisms of rice senescence remain obscure. In this study, we characterized an early-senescence rice mutant, designated zj-es (ZheJing-early senescence), which was derived from the japonica rice cultivar Zhejing22. The mutant zj-es exhibited obvious early-senescence phenotype, such as collapsed chloroplast, lesions in leaves, declined fertility, plant dwarf, and decreased agronomic traits. The ZJ-ES gene was mapped in a 458 kb-interval between the molecular markers RM5992 and RM5813 on Chromosome 3, and analysis suggested that ZJ-ES is a novel gene controlling rice early senescence. Subsequently, whole-transcriptome RNA sequencing was performed on zj-es and its wild-type rice to dissect the underlying molecular mechanism for early senescence. Totally, 10,085 differentially expressed mRNAs (DEmRNAs), 1,253 differentially expressed lncRNAs (DElncRNAs), and 614 differentially expressed miRNAs (DEmiRNAs) were identified, respectively, in different comparison groups. Based on the weighted gene co-expression network analysis (WGCNA), the co-expression turquoise module was found to be the key for the occurrence of rice early senescence. Furthermore, analysis on the competing endogenous RNA (CeRNA) network revealed that 14 lncRNAs possibly regulated 16 co-expressed mRNAs through 8 miRNAs, and enrichment analysis showed that most of the DEmRNAs and the targets of DElncRNAs and DEmiRNAs were involved in reactive oxygen species (ROS)-triggered autophagy-related pathways. Further analysis showed that, in zj-es, ROS-related enzyme activities were markedly changed, ROS were largely accumulated, autophagosomes were obviously observed, cell death was significantly detected, and lesions were notably appeared in leaves. Totally, combining our results here and the remaining research, we infer that ROS-triggered autophagy induces the programmed cell death (PCD) and its coupled early senescence in zj-es mutant rice.

Expression rewiring and methylation of non-coding RNAs involved in rhizome phenotypic variations of lotus ecotypes

Non-coding RNAs (ncRNAs), including miRNAs, lncRNAs, and circRNAs, emerge as crucial components for gene regulation. Nelumbo nucifera (lotus), a horticulturally important plant, differentiates into a temperate ecotype of enlarged rhizomes and a tropical ecotype of thin rhizomes. Nevertheless, whether and how ncRNAs can be rewired in expression and differentially methylated contributing to adaptive divergence of this storage organ in lotus ecotypes is unclear. Herein, we study the expression behaviors and DNA methylation patterns of ncRNAs in temperate and tropical lotus rhizomes. By whole transcriptome sequencing, we found both mRNAs and lncRNAs have divergent expression patterns between ecotypes, whereas miRNAs and circRNAs tended to be accession-specific or noisier in expression. The differentially expressed ncRNAs are involved in phenotypic differentiation of lotus rhizome between ecotypes, as the genes that interacted with them in the competing endogenous RNA network are enriched in functions including carbohydrate metabolism and plant hormone signaling, being critical to rhizome enlargement. Intriguingly, ncRNA-targeted genes are less prone to show positive selection or differential expression during ecotypic divergence due to constraints from ncRNA-mRNA interactions. The methylation levels of ncRNAs generally tend to be higher in temperate lotus than in tropical lotus, and differential methylation of lncRNAs also tends to have expression changes. Overall, our study of ncRNAs and their targets highlights the role of ncRNAs in rhizome growth variation between lotus ecotypes through expression rewiring and methylation modification.

The Characters of Non-Coding RNAs and Their Biological Roles in Plant Development and Abiotic Stress Response

Plant growth and development are greatly affected by the environment. Many genes have been identified to be involved in regulating plant development and adaption of abiotic stress. Apart from protein-coding genes, more and more evidence indicates that non-coding RNAs (ncRNAs), including small RNAs and long ncRNAs (lncRNAs), can target plant developmental and stress-responsive mRNAs, regulatory genes, DNA regulatory regions, and proteins to regulate the transcription of various genes at the transcriptional, posttranscriptional, and epigenetic level. Currently, the molecular regulatory mechanisms of sRNAs and lncRNAs controlling plant development and abiotic response are being deeply explored. In this review, we summarize the recent research progress of small RNAs and lncRNAs in plants, focusing on the signal factors, expression characters, targets functions, and interplay network of ncRNAs and their targets in plant development and abiotic stress responses. The complex molecular regulatory pathways among small RNAs, lncRNAs, and targets in plants are also discussed. Understanding molecular mechanisms and functional implications of ncRNAs in various abiotic stress responses and development will benefit us in regard to the use of ncRNAs as potential character-determining factors in molecular plant breeding.

LncRNA-mediated ceRNA networks provide novel potential biomarkers for peanut drought tolerance

Drought stress has been the major constraint on peanut yield and quality, and an understanding of the function of long non-coding (lncRNAs) in the peanut drought stress response is still in its infancy. In this study, two peanut varieties with contrasting drought tolerance were used to explore the functions of lncRNAs in the peanut drought response, and the results showed that the drought-tolerant variety presented greater antioxidant enzyme activity, osmotic adjustment ability, and photosynthesis under drought conditions. There were 4329 lncRNAs identified in the two varieties, of which 535 and 663 lncRNAs were differentially expressed in NH5 and FH18, respectively. The cis targets of the differentially expressed lncRNAs were putatively involved in secondary metabolite biosynthesis and other basic metabolic processes. A total of 673 competing endogenous RNA (ceRNA) pairs were selected specifically in NH5, and the associated ceRNA network revealed six lncRNAs, MSTRG.70535.2, MSTRG.86570.2, MSTRG.86570.1, MSTRG.100618.1, MSTRG.81214.2, and MSTRG.30931.1were considered as hub nodes. They were speculated to contribute to enhancing peanut drought tolerance, such as regulating transcription and plant growth processes, thereby improving the drought stress response. In this study, lncRNAs and mRNAs interaction networks were constructed to aid a comprehensive understanding of the peanut drought stress response and form a basis for future research.

Identification and characterization of circRNAs in maize seedlings under deficient nitrogen

Here, deep sequencing results of the maize transcriptome in leaves and roots were compared under high-nitrogen (HN) and low-nitrogen (LN) conditions to identify differentially expressed circRNAs (DECs). Circular RNAs (circRNAs) are covalently closed non-coding RNA with widely regulatory potency that has been identified in animals and plants. However, the understanding of circRNAs involved in responsive nitrogen deficiency remains to be elucidated. A total of 24 and 22 DECs were obtained from the leaves and roots, respectively. Ten circRNAs were validated by divergent and convergent primers, and 6 DECs showed the same expression tendency validated by reverse transcriptase-quantitative PCR. Integrating the identified differentially expressed miRNAs, 34 circRNAs could act as miRNA decoys, which might play important roles in multiple biological processes, including organonitrogen compound biosynthesis and regulation of the metabolic process. A total of 51 circRNA-parent genes located in the genome-wide association study identified loci were assessed between HN and LN conditions and were associated with root growth and development. In summary, our results provide valuable information regarding further study of maize circRNAs under nitrogen deficiency and provide new insights into screening of candidate genes as well as the improvement of maize regarding nitrogen deficiency resistance. CircRNA-miRNA-mRNA co-expression networks were constructed to explore the circRNAs that participated in biological development and nitrogen metabolism.

NGS Methodologies and Computational Algorithms for the Prediction and Analysis of Plant Circular RNAs

Circular RNAs (circRNAs) are a class of single-stranded RNAs derived from exonic, intronic, and intergenic regions from precursor messenger RNAs (pre-mRNA), where a noncanonical back-splicing event occurs, in which the 5' and 3' ends are attached by covalent bond. CircRNAs participate in the regulation of gene expression at the transcriptional and posttranscriptional level primarily as miRNA and RNA-binding protein (RBP) sponges, but also involved in the regulation of alternative RNA splicing and transcription. CircRNAs are widespread and abundant in plants where they have been involved in stress responses and development. Through the analysis of all publications in this field in the last five years, we can summarize that the identification of these molecules is carried out through next generation sequencing studies, where samples have been previously treated to eliminate DNA, rRNA, and linear RNAs as a means to enrich circRNAs. Once libraries are prepared, they are sequenced and subsequently studied from a bioinformatics point of view. Among the different tools for identifying circRNAs, we can highlight CIRI as the most used (in 60% of the published studies), as well as CIRCExplorer (20%) and find_circ (20%). Although it is recommended to use more than one program in combination, and preferably developed specifically to treat with plant samples, this is not always the case. It should also be noted that after identifying these circular RNAs, most of the authors validate their findings in the laboratory in order to obtain bona fide results.

Circular RNAs acting as ceRNAs mediated by miRNAs may be involved in the synthesis of soybean fatty acids

Soybean oil is composed of fatty acids and glycerol. The content and composition of fatty acids partly determine the quality of soybean seeds. Circular RNAs (circRNAs) are endogenous non-coding RNAs that competitively bind to microRNAs (miRNAs) through miRNA recognition elements, thereby acting as sponges to regulate the expression of target genes. Although circRNAs have been identified previously in soybean, only their expression has been investigated without exploration of the competitive endogenous RNAs (ceRNAs) network of circRNAs-miRNAs-mRNAs. In this study, circRNAs in immature pods of a low linolenic acid soybean Mutant 72' (MT72) and the wild-type control 'Jinong 18' (JN18) were systematically identified and analyzed at 30 and 40 days after flowering using high-throughput sequencing technology. We identified 6377 circRNAs, of which 114 were differentially expressed. Gene ontology and KEGG pathway analyses of targeted mRNAs in the ceRNAs network indicated that the differentially expressed circRNAs may be involved in fatty acid transport, suggesting that circRNAs may play a post-transcriptional regulatory role in soybean oil synthesis. This study provides a foundation for future exploration of the function of circRNAs in soybean and presents novel insights to guide further studies of plant circRNAs.

Genome-Wide Identification and Characterization of Potato Long Non-coding RNAs Associated With Phytophthora infestans Resistance

Long non-coding RNA (lncRNA) is a crucial regulatory mechanism in the plant response to biotic and abiotic stress. However, their roles in potato (Solanum tuberosum L.) resistance to Phytophthora infestans (P. infestans) largely remain unknown. In this study, we identify 2857 lncRNAs and 33,150 mRNAs of the potato from large-scale published RNA sequencing data. Characteristic analysis indicates a similar distribution pattern of lncRNAs and mRNAs on the potato chromosomes, and the mRNAs were longer and had more exons than lncRNAs. Identification of alternative splicing (AS) shows that there were a total of 2491 lncRNAs generated from AS and the highest frequency (46.49%) of alternative acceptors (AA). We performed R package TCseq to cluster 133 specific differentially expressed lncRNAs from resistance lines and found that the lncRNAs of cluster 2 were upregulated. The lncRNA targets were subject to KEGG pathway enrichment analysis, and the interactive network between lncRNAs and mRNAs was constructed by using GENIE3, a random forest machine learning algorithm. Transient overexpression of StLNC0004 in Nicotiana benthamiana significantly suppresses P. infestans growth compared with a control, and the expression of extensin (NbEXT), the ortholog of the StLNC0004 target gene, was significantly upregulated in the overexpression line. Together, these results suggest that lncRNAs play potential functional roles in the potato response to P. infestans infection.

Analysis of Rice Transcriptome Reveals the LncRNA/CircRNA Regulation in Tissue Development

BACKGROUND

Long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) can play important roles in many biological processes. However, no study of the influence of epigenetics factors or the 3D structure of the genome in their regulation is available in plants.

RESULTS

In the current analysis, we identified a total of 15,122 lncRNAs and 7902 circRNAs in three tissues (root, leaf and panicle) in the rice varieties Minghui 63, Zhenshan 97 and their hybrid Shanyou 63. More than 73% of these lncRNAs and parental genes of circRNAs (P-circRNAs) are shared among Oryza sativa with high expression specificity. We found that, compared with protein-coding genes, the loci of these lncRNAs have higher methylation levels and the loci of circRNAs tend to locate in the middle of genes with high CG and CHG methylation. Meanwhile, the activated lncRNAs and P-circRNAs are mainly transcribed from demethylated regions containing CHH methylation. In addition, ~ 53% lncRNAs and ~ 15% P-circRNAs are associated with transposable elements (TEs), especially miniature inverted-repeat transposable elements and RC/Helitron. We didn't find correlation between the expression of lncRNAs and histone modifications; however, we found that the binding strength and interaction of RNAPII significantly affects lncRNA expression. Interestingly, P-circRNAs tend to combine active histone modifications. Finally, we found that lncRNAs and circRNAs acting as competing-endogenous RNAs have the potential to regulate the expression of genes, such as osa-156 l-5p (related to yield) and osa-miR444a-3p (related to N/P metabolism) confirmed through dual-luciferase reporter assays, with important roles in the growth and development of rice, laying a foundation for future rice breeding analyses.

CONCLUSIONS

In conclusion, our study comprehensively analyzed the important regulatory roles of lncRNA/circRNA in the tissue development of Indica rice from multiple perspectives.

Identification of circRNAs in the Liver of Whitespotted Bamboo Shark (Chiloscyllium plagiosum)

Whitespotted bamboo shark (Chiloscyllium plagiosum), a member of the cartilaginous fish family, has an extremely large liver and demonstrates a strong regeneration ability and immune regulation. Circular RNAs (circRNAs) is an important class of non-coding RNAs. Increasing evidences suggest that circRNAs are a kind of potential regulators. Recently, researchers have isolated and identified different circRNAs from various species, while few reports were on the circRNAs of C. plagiosum. In this study, we have identified a total of 4,558 circRNAs in the liver of C. plagiosum. This finding suggests that circRNAs are not evenly distributed in the chromosomes and follow the GT-AG rule during cyclization. Alternative back-splicing might exist in shark circRNAs as shown by the authenticity identification of predicted circRNAs. The binding strength of circRNAs (<2,000 bp) and the detected miRNAs in shark liver were simultaneously analyzed to construct an mRNA–miRNA–circRNA network for the Glutathione S-transferase P1 gene, and the circRNA authenticity was simultaneously verified. Our data provide not only novel insights into the rich existence of circRNAs in marine animals, but also a basis for characterizing functions of identified circRNAs in the liver homeostasis of C. plagiosum.

Integrated Transcriptome, microRNA, and Phytochemical Analyses Reveal Roles of Phytohormone Signal Transduction and ABC Transporters in Flavor Formation of Oolong Tea (Camellia sinensis) during Solar Withering

The unique aroma and flavor of oolong tea develop during the withering stage of postharvest processing. We explored the roles of miRNA-related regulatory networks during tea withering and their effects on oolong tea quality. We conducted transcriptome and miRNA analyses to identify differentially expressed (DE) miRNAs and target genes among fresh leaves, indoor-withered leaves, and solar-withered leaves. We identified 32 DE-miRNAs and 41 target genes involved in phytohormone signal transduction and ABC transporters. Further analyses indicated that these two pathways regulated the accumulation of flavor-related metabolites during tea withering. Flavonoid accumulation was correlated with the miR167d_1-ARF-GH3, miR845-ABCC1-3/ABCC2, miR166d-5p_1-ABCC1-2, and miR319c_3-PIF-ARF modules. Terpenoid content was correlated with the miR171b-3p_2-DELLA-MYC2 and miR166d-5p_1-ABCG2-MYC2 modules. These modules inhibited flavonoid biosynthesis and enhanced terpenoid biosynthesis in solar-withered leaves. Low auxin and gibberellic acid contents and circRNA-related regulatory networks also regulated the accumulation of flavor compounds in solar-withered leaves. Our analyses reveal how solar withering produces high-quality oolong tea.

CircPlant: An Integrated Tool for circRNA Detection and Functional Prediction in Plants

The recent discovery of circular RNAs (circRNAs) and characterization of their functional roles have opened a new avenue for understanding the biology of genomes. circRNAs have been implicated to play important roles in a variety of biological processes, but their precise functions remain largely elusive. Currently, a few approaches are available for novel circRNA prediction, but almost all these methods are intended for animal genomes. Considering that the major differences between the organization of plant and mammal genomes cannot be neglected, a plant-specific method is needed to enhance the validity of plant circRNA identification. In this study, we present CircPlant, an integrated tool for the exploration of plant circRNAs, potentially acting as competing endogenous RNAs (ceRNAs), and their potential functions. With the incorporation of several unique plant-specific criteria, CircPlant can accurately detect plant circRNAs from high-throughput RNA-seq data. Based on comparison tests on simulated and real RNA-seq datasets from Arabidopsis thaliana and Oryza sativa, we show that CircPlant outperforms all evaluated competing tools in both accuracy and efficiency. CircPlant is freely available at http://bis.zju.edu.cn/circplant.

Splicing conservation signals in plant long noncoding RNAs

Long noncoding RNAs (lncRNAs) have recently emerged as prominent regulators of gene expression in eukaryotes. LncRNAs often drive the modification and maintenance of gene activation or gene silencing states via chromatin conformation rearrangements. In plants, lncRNAs have been shown to participate in gene regulation, and are essential to processes such as vernalization and photomorphogenesis. Despite their prominent functions, only over a dozen lncRNAs have been experimentally and functionally characterized. Similar to its animal counterparts, the rates of sequence divergence are much higher in plant lncRNAs than in protein coding mRNAs, making it difficult to identify lncRNA conservation using traditional sequence comparison methods. Beyond this, little is known about the evolutionary patterns of lncRNAs in plants. Here, we characterized the splicing conservation of lncRNAs in Brassicaceae. We generated a whole-genome alignment of 16 Brassica species and used it to identify synthenic lncRNA orthologs. Using a scoring system trained on transcriptomes from A. thaliana and B. oleracea, we identified splice sites across the whole alignment and measured their conservation. Our analysis revealed that 17.9% (112/627) of all intergenic lncRNAs display splicing conservation in at least one exon, an estimate that is substantially higher than previous estimates of lncRNA conservation in this group. Our findings agree with similar studies in vertebrates, demonstrating that splicing conservation can be evidence of stabilizing selection. We provide conclusive evidence for the existence of evolutionary deeply conserved lncRNAs in plants and describe a generally applicable computational workflow to identify functional lncRNAs in plants.

Characterization and Function of Circular RNAs in Plants

CircRNAs are covalently closed-loop single-stranded RNA molecules ubiquitously expressing in eukaryotes. As an important member of the endogenous ncRNA family, circRNAs are associated with diverse biological processes and can regulate transcription, modulate alternative splicing, and interact with miRNAs or proteins. Compared to abundant advances in animals, studies of circRNAs in plants are rapidly emerging. The databases and analysis tools for plant circRNAs are constantly being developed. Large numbers of circRNAs have been identified and characterized in plants and proved to play regulatory roles in plant growth, development, and stress responses. Here, we review the biogenesis, characteristics, bioinformatics resources, and biological functions of plant circRNAs, and summarize the distinct circularization features and differentially expression patterns comparison with animal-related results.

Bioinformatics Resources for Plant Abiotic Stress Responses: State of the Art and Opportunities in the Fast Evolving -Omics Era

Abiotic stresses are among the principal limiting factors for productivity in agriculture. In the current era of continuous climate changes, the understanding of the molecular aspects involved in abiotic stress response in plants is a priority. The rise of -omics approaches provides key strategies to promote effective research in the field, facilitating the investigations from reference models to an increasing number of species, tolerant and sensitive genotypes. Integrated multilevel approaches, based on molecular investigations at genomics, transcriptomics, proteomics and metabolomics levels, are now feasible, expanding the opportunities to clarify key molecular aspects involved in responses to abiotic stresses. To this aim, bioinformatics has become fundamental for data production, mining and integration, and necessary for extracting valuable information and for comparative efforts, paving the way to the modeling of the involved processes. We provide here an overview of bioinformatics resources for research on plant abiotic stresses, describing collections from -omics efforts in the field, ranging from raw data to complete databases or platforms, highlighting opportunities and still open challenges in abiotic stress research based on -omics technologies.

Comparative Transcriptome Analysis Reveals New lncRNAs Responding to Salt Stress in Sweet Sorghum

Long non-coding RNAs (lncRNAs) can enhance plant stress resistance by regulating the expression of functional genes. Sweet sorghum is a salt-tolerant energy crop. However, little is known about how lncRNAs in sweet sorghum respond to salt stress. In this study, we identified 126 and 133 differentially expressed lncRNAs in the salt-tolerant M-81E and the salt-sensitive Roma strains, respectively. Salt stress induced three new lncRNAs in M-81E and inhibited two new lncRNAs in Roma. These lncRNAs included lncRNA13472, lncRNA11310, lncRNA2846, lncRNA26929, and lncRNA14798, which potentially function as competitive endogenous RNAs (ceRNAs) that influence plant responses to salt stress by regulating the expression of target genes related to ion transport, protein modification, transcriptional regulation, and material synthesis and transport. Additionally, M-81E had a more complex ceRNA network than Roma. This study provides new information regarding lncRNAs and the complex regulatory network underlying salt-stress responses in sweet sorghum.

Genome-wide identification of long non-coding RNAs and circular RNAs reveal their ceRNA networks in response to cucumber green mottle mosaic virus infection in watermelon

Long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) play vital roles in plant defense responses against viral infections. However, there is no systematic understanding of lncRNAs and circRNAs and their competing endogenous RNA (ceRNA) networks in watermelon under cucumber green mottle mosaic virus (CGMMV) stress. Here, we present the characterization and expression profiles of lncRNAs and circRNAs in watermelon leaves 48-h post-inoculation (48 hpi) with CGMMV, with mock inoculation as a control. Deep sequencing analysis revealed 2373 lncRNAs and 606 circRNAs in the two libraries. Among them, 67 lncRNAs (40 upregulated and 27 downregulated) and 548 circRNAs (277 upregulated and 271 downregulated) were differentially expressed (DE) in the 48 hpi library compared with the control library. Furthermore, 263 cis-acting matched lncRNA-mRNA pairs were detected for 49 of the DE-lncRNAs. KEGG pathway analysis of the cis target genes of the DE-lncRNAs revealed significant associations with phenylalanine metabolism, the citrate cycle (TCA cycle), and endocytosis. Additionally, 30 DE-lncRNAs were identified as putative target mimics of 33 microRNAs (miRNAs), and 153 DE-circRNAs were identified as putative target mimics of 88 miRNAs. Furthermore, ceRNA networks of lncRNA/circRNA-miRNA-mRNA in response to CGMMV infection are described, with 12 DE-lncRNAs and 65 DE-circRNAs combining with 22 miRNAs and competing for the miRNA binding sites on 29 mRNAs. The qRT-PCR validation of selected lncRNAs and circRNAs showed a general correlation with the high-throughput sequencing results. This study provides a valuable resource of lncRNAs and circRNAs involved in the response to CGMMV infection in watermelon.

Genome-Wide Identification of lncRNAs and Analysis of ceRNA Networks During Tomato Resistance to Phytophthora infestans

Our previous studies have revealed the function of long noncoding RNAs (lncRNAs) and microRNAs (miRNAs) in tomato in response to Phytophthora infestans infection. However, the interaction relationships between lncRNAs and miRNAs during tomato resistance to P. infestans infection are unknown. In this study, 9,011 lncRNAs were identified from tomato plants, including 115 upregulated and 81 downregulated lncRNAs. Among these, 148 were found to be differentially expressed and might affect the expression of 771 genes, which are composed of 887 matched lncRNA-mRNA pairs. In total, 88 lncRNAs were identified as endogenous RNAs (ceRNAs) and predicted to decoy 46 miRNAs. Degradome sequencing revealed that 11 miRNAs that were decoyed by 20 lncRNAs could target 30 genes. These lncRNAs, miRNAs, and target genes were predicted to form 10 regulatory modules. Among them, lncRNA42705/lncRNA08711, lncRNA39896, and lncRNA11265/lncRNA15816 might modulate MYB, HD-Zip, and NAC transcription factors by decoying miR159, miR166b, and miR164a-5p, respectively. Upon P. infestans infection, the expression levels of lncRNA42705 and lncRNA08711 displayed a negative correlation with the expression level of miR159 and a positive correlation with the expression levels of MYB genes. Tomato plants in which lncRNA42705 and lncRNA08711 were silenced displayed increased levels of miR159 and decreased levels of MYB, respectively. The result demonstrated that lncRNAs might function as ceRNAs to decoy miRNAs and affect their target genes in tomato plants, increasing resistance to disease.

Non-Coding RNAs Participate in the Regulation of CRY-DASH in the Growth and Early Development of Saccharina japonica (Laminariales, Phaeophyceae)

CRY-DASH, a new cryptochrome blue light receptor, can repair damaged DNA and regulate secondary metabolism and development of fungus. However, its role in regulation during the growth of Saccharina japonica is still unclear. After cloning the full-length of CRY-DASH from S. japonica (sjCRY-DASH), we deduced that its open reading frame was 1779 bp long and encoded 592 amino acids. sjCRY-DASH transcription was rapidly upregulated within 30 min in response to blue light and exhibited 24 h periodicity with different photoperiods. Moreover, sjCRY-DASH maintained the same periodicity in suitable growth temperature, suggesting a close relationship between this periodicity and circadian rhythm regulation. Novel-m3234-5p, which was targeted to sjCRY-DASH, decreased with increasing sjCRY-DASH transcription, acting as a negative modulator of sjCRY-DASH. Six long non-coding RNAs classified as long intergenic non-coding RNAs (lincRNAs) exhibited co-expression with sjCRY-DASH. A miRNA sjCRY DASH lincRNA network was consequently identified. By predicting the endogenous competing mRNAs of novel-m3234-5p, we found that sjCRY-DASH indirectly participated in the regulation of DNA damage repair, protein synthesis and processing, and actin transport. In conclusion, our results revealed that non-coding RNAs participate in the regulation of sjCRY-DASH, which played vital roles in the growth and early development of S. japonica.

Whole-transcriptome RNA sequencing reveals the global molecular responses and ceRNA regulatory network of mRNAs, lncRNAs, miRNAs and circRNAs in response to copper toxicity in Ziyang Xiangcheng (Citrus junos Sieb. Ex Tanaka)

BACKGROUND

Copper (Cu) toxicity has become a potential threat for citrus production, but little is known about related mechanisms. This study aims to uncover the global landscape of mRNAs, long non-coding RNAs (lncRNAs), circular RNAs (circRNAs) and microRNAs (miRNAs) in response to Cu toxicity so as to construct a regulatory network of competing endogenous RNAs (ceRNAs) and to provide valuable knowledge pertinent to Cu response in citrus.

RESULTS

Tolerance of four commonly used rootstocks to Cu toxicity was evaluated, and 'Ziyang Xiangcheng' (Citrus junos) was found to be the most tolerant genotype. Then the roots and leaves sampled from 'Ziyang Xiangcheng' with or without Cu treatment were used for whole-transcriptome sequencing. In total, 5734 and 222 mRNAs, 164 and 5 lncRNAs, 45 and 17 circRNAs, and 147 and 130 miRNAs were identified to be differentially expressed (DE) in Cu-treated roots and leaves, respectively, in comparison with the control. Gene ontology enrichment analysis showed that most of the DEmRNAs and targets of DElncRNAs and DEmiRNAs were annotated to the categories of 'oxidation-reduction', 'phosphorylation', 'membrane', and 'ion binding'. The ceRNA network was then constructed with the predicted pairs of DEmRNAs-DEmiRNAs and DElncRNAs-DEmiRNAs, which further revealed regulatory roles of these DERNAs in Cu toxicity.

CONCLUSIONS

A large number of mRNAs, lncRNAs, circRNAs, and miRNAs in 'Ziyang Xiangcheng' were altered in response to Cu toxicity, which may play crucial roles in mitigation of Cu toxicity through the ceRNA regulatory network in this Cu-tolerant rootstock.

Identification and characterization of CircRNAs involved in the regulation of wheat root length

BACKGROUND

Recent studies indicate that circular RNAs (circRNAs) may play important roles in the regulation of plant growth and development. Plant roots are the main organs of nutrient and water uptake. However, whether circRNAs involved in the regulation of plant root growth remains to be elucidated.

METHODS

LH9, XN979 and YN29 are three Chinese wheat varieties with contrasting root lengths. Here, the root circRNA expression profiles of LH9, XN979 and YN29 were examined by using high-throughput sequencing technology.

RESULTS

Thirty-three and twenty-two differentially expressed circRNAs (DECs) were identified in the YN29-LH9 comparison and YN29-XN979 comparison, respectively. Among them, ten DECs coexisted in both comparisons. As the roots of both LH9 and XN979 were significantly larger and deeper than YN29, the ten DECs coexisting in the two comparisons were highly likely to be involved in the regulation of wheat root length. Moreover, three of the ten DECs have potential miRNA binding sites. Real-time PCR analysis showed that the expression levels of the potential binding miRNAs exhibited significant differences between the long root plants and the short root plants.

CONCLUSIONS

The expression levels of some circRNAs exhibited significant differences in wheat varieties with contrasting root phenotypes. Ten DECs involved in the regulation of wheat root length were successfully identified in which three of them have potential miRNAs binding sites. The expression levels of putative circRNA-binding miRNAs were correlated with their corresponding circRNAs. Our results provide new clues for studying the potential roles of circRNAs in the regulation of wheat root length.