Genome-Wide Identification and Functional Prediction of Novel Drought-Responsive lncRNAs in Pyrus betulifolia

Whole-Transcriptome Sequencing Reveals the Global Molecular Responses and NAC Transcription Factors Involved in Drought Stress in Dendrobium catenatum

Dendrobium catenatum is a highly drought-tolerant herb, which usually grows on cliffs or in the branches of trees, yet the underlying molecular mechanisms for its tolerance remain poorly understood. We conducted a comprehensive study utilizing whole-transcriptome sequencing approaches to investigate the molecular response to extreme drought stress in D. catenatum. A large number of differentially expressed mRNAs, lncRNAs, and circRNAs have been identified, and the NAC transcription factor family was highly enriched. Meanwhile, 46 genes were significantly up-regulated in the ABA-activated signaling pathway. In addition to the 89 NAC family members accurately identified in this study, 32 members were found to have different expressions between the CK and extreme drought treatment. They may regulate drought stress through both ABA-dependent and ABA-independent pathways. Moreover, the 32 analyzed differentially expressed DcNACs were found to be predominantly expressed in the floral organs and roots. The ceRNA regulatory network showed that DcNAC87 is at the core of the ceRNA network and is regulated by miR169, miR393, and four lncRNAs. These investigations provided valuable information on the role of NAC transcription factors in D. catenatum's response to drought stress.

Long non-coding RNA-mediated epigenetic response for abiotic stress tolerance in plants

Plants perceive environmental fluctuations as stress and confront several stresses throughout their life cycle individually or in combination. Plants have evolved their sensing and signaling mechanisms to perceive and respond to a variety of stresses. Epigenetic regulation plays a critical role in the regulation of genes, spatiotemporal expression of genes under stress conditions and imparts a stress memory to encounter future stress responses. It is quintessential to integrate our understanding of genetics and epigenetics to maintain plant fitness, achieve desired genetic gains with no trade-offs, and durable long-term stress tolerance. The long non-coding RNA >200 nts having no coding potential (or very low) play several roles in epigenetic memory, contributing to the regulation of gene expression and the maintenance of cellular identity which include chromatin remodeling, imprinting (dosage compensation), stable silencing, facilitating nuclear organization, regulation of enhancer-promoter interactions, response to environmental signals and epigenetic switching. The lncRNAs are involved in a myriad of stress responses by activation or repression of target genes and hence are potential candidates for deploying in climate-resilient breeding programs. This review puts forward the significant roles of long non-coding RNA as an epigenetic response during abiotic stresses in plants and the prospects of deploying lncRNAs for designing climate-resilient plants.

Analysis of lncRNAs in Lupinus mutabilis (Tarwi) and Their Potential Role in Drought Response

Lupinus mutabilis is a legume with high agronomic potential and available transcriptomic data for which lncRNAs have not been studied. Therefore, our objective was to identify, characterize, and validate the drought-responsive lncRNAs in L. mutabilis. To achieve this, we used a multilevel approach based on lncRNA prediction, annotation, subcellular location, thermodynamic characterization, structural conservation, and validation. Thus, 590 lncRNAs were identified by at least two algorithms of lncRNA identification. Annotation with the PLncDB database showed 571 lncRNAs unique to tarwi and 19 lncRNAs with homology in 28 botanical families including Solanaceae (19), Fabaceae (17), Brassicaceae (17), Rutaceae (17), Rosaceae (16), and Malvaceae (16), among others. In total, 12 lncRNAs had homology in more than 40 species. A total of 67% of lncRNAs were located in the cytoplasm and 33% in exosomes. Thermodynamic characterization of S03 showed a stable secondary structure with -105.67 kcal/mol. This structure included three regions, with a multibranch loop containing a hairpin with a SECIS-like element. Evaluation of the structural conservation by CROSSalign revealed partial similarities between L. mutabilis (S03) and S. lycopersicum (Solyc04r022210.1). RT-PCR validation demonstrated that S03 was upregulated in a drought-tolerant accession of L. mutabilis. Finally, these results highlighted the importance of lncRNAs in tarwi improvement under drought conditions.

Identification and characterization of long noncoding RNAs involved in the aluminum stress response in Medicago truncatula via genome-wide analysis

Numerous studies have shown that plant long noncoding RNAs (lncRNAs) play an important regulatory role in the plant response to environmental stress. However, there are no reports on lncRNAs regulating and enhancing aluminum (Al) stress tolerance in legumes. This study analyzed the role of lncRNAs in response to Al stress in the legume model plant Medicago truncatula. A total of 219.49 Gb clean data were generated: 3,284 lncRNA genes were identified, of which 515 were differentially expressed, and 1,254 new genes were functionally annotated through database alignment. We further predicted and classified putative targets of these lncRNAs and found that they were enriched in biological processes and metabolic pathways such as plant hormone signal transduction, cell wall modification and the tricarboxylic acid (TCA) cycle. Finally, we characterized the functions of 2 Al-activated-malate-transporter-related lncRNAs in yeast. The recombinant plasmids of MSTRG.12506.5 and MSTRG.34338.20 were transformed into yeast, and these yeast exhibited better growth than those carrying empty vectors on medium supplemented with 10 μM AlCl₃ and showed that they have biological functions affording Al stress tolerance. These findings suggest that lncRNAs are involved in regulating plant responses to Al stress. Our findings help to understand the role of lncRNAs in the response to Al stress in legumes and provide candidate lncRNAs for further studies.

Genome-wide identification of long intergenic non-coding RNAs for Ralstonia solanacearum resistance in tomato (Solanum lycopersicum)

There is growing evidences indicating that long intergenic ncRNAs (lincRNAs) play key roles in plant development and stress responses. To research tomato lincRNA functions during the interaction between tomato and Ralstonia solanacearum, RNA-seq data of tomato plants inoculated with R. solanacearum was analyzed. In this study, 315 possible lincRNAs were identified from RNA-seq data. Then 23 differentially expressed lincRNAs between tomato plants inoculated with R. solanacearum and control were identified and a total of 171 possible target genes for these differentially expressed lincRNAs were predicted. Through GO and KEGG analysis, we found that lincRNA might be involved in jasmonic acid and ethylene signaling pathways to respond to tomato bacterial wilt infection. Furthermore, lincRNA may also be involved in regulating the expression of AGO protein. Subsequently, analysis of expression patterns between differentially expressed lincRNAs and adjacent mRNAs by qRT-PCR revealed that part of lincRNAs and their possible target genes exhibited positive correlation. Taken together, these results suggest that lincRNAs play potential roles in tomato against R. solanacearum infection and will provide fundamental information about the lincRNA-based plant defense mechanisms.

Long non-coding RNA transcriptome landscape of anthers at different developmental stages in response to drought stress in tomato

Pollen development is particularly susceptible to drought stress. Long non-coding RNAs (lncRNAs) are known to play a role in plant development and responses to drought stress. However, the expression profile and putative function of lncRNAs in drought-induced male sterility remain largely unknown. In this study, we investigated the lncRNA transcriptome landscape of tomato anthers at early and late stages-tetrad-vacuolated microspore (TED-VUM) and binucleate-mature pollen (BIN-MP) anthers, respectively-in response to drought stress using RNA-sequencing. In total, we identified 67,770 lncRNAs, of which 3053 lncRNAs were drought responsive. Interestingly, there were more differentially expressed (DE) lncRNAs in TED-VUM (2879) than in BIN-MP (174) anthers, which was consistent with the TED-VUM anthers being more drought sensitive. Functional enrichment analysis revealed that the target genes of DE lncRNAs were significantly enriched in diverse metabolic processes, including in carbohydrate metabolism and hormone synthesis. Co-expression analysis also identified 1407 lncRNAs that strongly co-expressed with 8 target genes that are involved in hormone (abscisic acid and jasmonic acid) and carbohydrate (sucrose and starch) metabolisms and tapetum development, highlighting the potential of lncRNA-target-gene modulation of anther development under drought stress. Our results serve as a baseline for future investigations of the potential function of lncRNAs in plant reproductive development under drought stress.

Drought tolerance improvement in Solanum lycopersicum: an insight into "OMICS" approaches and genome editing

Solanum lycopersicum (tomato) is an internationally acclaimed vegetable crop that is grown worldwide. However, drought stress is one of the most critical challenges for tomato production, and it is a crucial task for agricultural biotechnology to produce drought-resistant cultivars. Although breeders have done a lot of work on the tomato to boost quality and quantity of production and enhance resistance to biotic and abiotic stresses, conventional tomato breeding approaches have been limited to improving drought tolerance because of the intricacy of drought traits. Many efforts have been made to better understand the mechanisms involved in adaptation and tolerance to drought stress in tomatoes throughout the years. "Omics" techniques, such as genomics, transcriptomics, proteomics, and metabolomics in combination with modern sequencing technologies, have tremendously aided the discovery of drought-responsive genes. In addition, the availability of biotechnological tools, such as plant transformation and the recently developed genome editing system for tomatoes, has opened up wider opportunities for validating the function of drought-responsive genes and the generation of drought-tolerant varieties. This review highlighted the recent progresses for tomatoes improvement against drought stress through "omics" and "multi-omics" technologies including genetic engineering. We have also discussed the roles of non-coding RNAs and genome editing techniques for drought stress tolerance improvement in tomatoes.

Identification and functional deciphering suggested the regulatory roles of long intergenic ncRNAs (lincRNAs) in increasing grafting pepper resistance to Phytophthora capsici

BACKGROUND

As a popular and valuable technique, grafting is widely used to protect against soil-borne diseases and nematodes in vegetable production. Growing evidences have revealed that long intergenic ncRNAs (lincRNAs) are strictly regulated and play essential roles in plants development and stress responses. Nevertheless, genome-wide identification and function deciphering of pepper lincRNAs, especially for their roles in improving grafting pepper resistance to Phytophthora capsici is largely unknown.

RESULTS

In this study, RNA-seq data of grafting and control pepper plants with or without P. capsici inoculation were used to identify lincRNAs. In total, 2,388 reliable lincRNAs were identified. They were relatively longer and contained few exons than protein-coding genes. Similar to coding genes, lincRNAs had higher densities in euchromatin regions; and longer chromosome transcribed more lincRNAs. Expression pattern profiling suggested that lincRNAs commonly had lower expression than mRNAs. Totally, 607 differentially expressed lincRNAs (DE-lincRANs) were identified, of which 172 were found between P. capsici resistance grafting pepper sample GR and susceptible sample LDS. The neighboring genes of DE-lincRNAs and miRNAs competitively sponged by DE-lincRNAs were identified. Subsequently, the expression level of DE-lincRNAs was further confirmed by qRT-PCR and regulation patterns between DE-lincRNAs and neighboring mRNAs were also validated. Function annotation revealed that DE-lincRNAs increased the resistance of grafting prepper to P. capsici by modulating the expression of disease-defense related genes through cis-regulating and/or lincRNA-miRNA-mRNA interaction networks.

CONCLUSIONS

This study identified pepper lincRNAs and suggested their potential roles in increasing the resistance level of grafting pepper to P. capsici.

Non-Coding RNAs in Response to Drought Stress

Drought stress causes changes in the morphological, physiological, biochemical and molecular characteristics of plants. The response to drought in different plants may vary from avoidance, tolerance and escape to recovery from stress. This response is genetically programmed and regulated in a very complex yet synchronized manner. The crucial genetic regulations mediated by non-coding RNAs (ncRNAs) have emerged as game-changers in modulating the plant responses to drought and other abiotic stresses. The ncRNAs interact with their targets to form potentially subtle regulatory networks that control multiple genes to determine the overall response of plants. Many long and small drought-responsive ncRNAs have been identified and characterized in different plant varieties. The miRNA-based research is better documented, while lncRNA and transposon-derived RNAs are relatively new, and their cellular role is beginning to be understood. In this review, we have compiled the information on the categorization of non-coding RNAs based on their biogenesis and function. We also discuss the available literature on the role of long and small non-coding RNAs in mitigating drought stress in plants.

Comparative Transcriptomic Analyses Provide Insights into the Enzymatic Browning Mechanism of Fresh-Cut Sand Pear Fruit

Pear (Pyrus spp.) is one of the most commonly consumed temperate fruits, having considerable economic and health importance. Fresh-cut or processed pear fruits are prone to browning because of the abundant phenolic compounds; however, little is known about the molecular mechanisms underlying enzymatic browning of fresh-cut sand pear fruit. In this study, fruits of two sand pear genotypes (low browning cultivar ‘Eli No.2′ and high browning cultivar ‘Weiningdahuangli’) were used to analyze the molecular mechanism of enzymatic browning by SMRT-seq and RNA-seq. The results generated 69,122 consensus isoforms, 21,336 new transcripts, 7105 alternative splicing events, and 254 long non-coding RNAs (lncRNAs). Furthermore, five genes related to enzymatic browning were predicted to be targets of six lncRNAs, and 9930 differentially expressed genes (DEGs) were identified between two different flesh browning cultivars. Meanwhile, most DEGs (e.g., PAL, 4CL, CAD, CCR, CHS, and LAR) involved in the phenylpropanoid biosynthesis pathway were up-regulated, and the expression of PPO and POD were highly expressed in the high-browning cultivar. Interestingly, the transcript level of PbrPPO4 (Pbr000321.4) was significantly higher than other PPO and POD genes, and a high level of total polyphenol and PPO activity were observed in the high browning cultivar. We found that the expression of lncRNA PB.156.1 was significantly positively correlated with the target gene PbrPPO4 (Pbr000321.4). The results suggest that PbrPPO4 might act as a major contributor and a key enzyme encoding gene in regulating fresh-cut sand pear fruit enzymatic browning; the expression of PbrPPO4 was probably regulated by lncRNA PB.156.1. Altogether, the transcriptomic and physiological analyses expand the knowledge of sand pear flesh enzymatic browning at the molecular level and provide a foundation for germplasm resources for molecular breeding of high polyphenol and low browning cultivars in sand pears.

The transcriptomic response to heat stress of a jujube (Ziziphus jujuba Mill.) cultivar is featured with changed expression of long noncoding RNAs

Long non-coding RNA (lncRNA) of plant species undergoes dynamic regulation and acts in developmental and stress regulation. Presently, there is little information regarding the identification of lncRNAs in jujube (Ziziphus jujuba Mill.), and it is uncertain whether the lncRNAs could respond to heat stress (HS) or not. In our previous study, a cultivar (Hqing1-HR) of Z. jujuba were treated by HS (45°C) for 0, 1, 3, 5 and 7 days, and it was found that HS globally changed the gene expression by RNA sequencing (RNA-seq) experiments and informatics analyses. In the current study, 8260 lncRNAs were identified successfully from the previous RNA-seq data, and it indicated that lncRNAs expression was also altered globally, suggesting that the lncRNAs might play vital roles in response to HS. Furthermore, bioinformatics analyses of potential target mRNAs of lncRNAs with cis-acting mechanism were performed, and it showed that multiple differentially expressed (DE) mRNAs co-located with DElncRNAs were highly enriched in pathways associated with response to stress and regulation of metabolic process. Taken together, these findings not only provide a comprehensive identification of lncRNAs but also useful clues for molecular mechanism response to HS in jujube.

Single-Molecule Long-Read Sequencing of Purslane (Portulaca oleracea) and Differential Gene Expression Related with Biosynthesis of Unsaturated Fatty Acids

This study aimed to obtain the full-length transcriptome of purslane (Portulaca oleracea); assorted plant samples were used for single-molecule real-time (SMRT) sequencing. Based on SMRT, functional annotation of transcripts, transcript factors (TFs) analysis, simple sequence repeat analysis and long non-coding RNAs (LncRNAs) prediction were accomplished. Total 15.33-GB reads were produced; with 9,350,222 subreads and the average length of subreads, 1640 bp was counted. With 99.99% accuracy, after clustering, 132,536 transcripts and 78,559 genes were detected. All unique SMART transcripts were annotated in seven functional databases. 4180 TFs (including transcript regulators) and 7289 LncRNAs were predicted. The results of RNA-seq were confirmed with qRT–PCR analysis. Illumina sequencing of leaves and roots of two purslane genotypes was carried out. Amounts of differential expression genes and related KEGG pathways were found. The expression profiles of related genes in the biosynthesis of unsaturated fatty acids pathway in leaves and roots of two genotypes of purslane were analyzed. Differential expression of genes in this pathway built the foundation of ω-3 fatty acid accumulation in different organs and genotypes of purslane. The aforementioned results provide sequence information and may be a valuable resource for whole-genome sequencing of purslane in the future.

Long non-coding RNAs: emerging players regulating plant abiotic stress response and adaptation

BACKGROUND

The immobile nature of plants means that they can be frequently confronted by various biotic and abiotic stresses during their lifecycle. Among the various abiotic stresses, water stress, temperature extremities, salinity, and heavy metal toxicity are the major abiotic stresses challenging overall plant growth. Plants have evolved complex molecular mechanisms to adapt under the given abiotic stresses. Long non-coding RNAs (lncRNAs)-a diverse class of RNAs that contain > 200 nucleotides(nt)-play an essential role in plant adaptation to various abiotic stresses.

RESULTS

LncRNAs play a significant role as 'biological regulators' for various developmental processes and biotic and abiotic stress responses in animals and plants at the transcription, post-transcription, and epigenetic level, targeting various stress-responsive mRNAs, regulatory gene(s) encoding transcription factors, and numerous microRNAs (miRNAs) that regulate the expression of different genes. However, the mechanistic role of lncRNAs at the molecular level, and possible target gene(s) contributing to plant abiotic stress response and adaptation, remain largely unknown. Here, we review various types of lncRNAs found in different plant species, with a focus on understanding the complex molecular mechanisms that contribute to abiotic stress tolerance in plants. We start by discussing the biogenesis, type and function, phylogenetic relationships, and sequence conservation of lncRNAs. Next, we review the role of lncRNAs controlling various abiotic stresses, including drought, heat, cold, heavy metal toxicity, and nutrient deficiency, with relevant examples from various plant species. Lastly, we briefly discuss the various lncRNA databases and the role of bioinformatics for predicting the structural and functional annotation of novel lncRNAs.

CONCLUSIONS

Understanding the intricate molecular mechanisms of stress-responsive lncRNAs is in its infancy. The availability of a comprehensive atlas of lncRNAs across whole genomes in crop plants, coupled with a comprehensive understanding of the complex molecular mechanisms that regulate various abiotic stress responses, will enable us to use lncRNAs as potential biomarkers for tailoring abiotic stress-tolerant plants in the future.

Genome-wide analysis of long non-coding RNAs (lncRNAs) in two contrasting rapeseed (Brassica napus L.) genotypes subjected to drought stress and re-watering

BACKGROUND

Drought stress is a major abiotic factor that affects rapeseed (Brassica napus L.) productivity. Though previous studies indicated that long non-coding RNAs (lncRNAs) play a key role in response to drought stress, a scheme for genome-wide identification and characterization of lncRNAs' response to drought stress is still lacking, especially in the case of B. napus. In order to further understand the molecular mechanism of the response of B. napus to drought stress, we compared changes in the transcriptome between Q2 (a drought-tolerant genotype) and Qinyou8 (a drought-sensitive genotype) responding drought stress and rehydration treatment at the seedling stage.

RESULTS

A total of 5546 down-regulated and 6997 up-regulated mRNAs were detected in Q2 compared with 7824 and 10,251 in Qinyou8, respectively; 369 down-regulated and 108 up- regulated lncRNAs were detected in Q2 compared with 449 and 257 in Qinyou8, respectively. LncRNA-mRNA interaction network analysis indicated that the co-expression network of Q2 was composed of 145 network nodes and 5175 connections, while the co-expression network of Qinyou8 was composed of 305 network nodes and 22,327 connections. We further identified 34 transcription factors (TFs) corresponding to 126 differentially expressed lncRNAs in Q2, and 45 TFs corresponding to 359 differentially expressed lncRNAs in Qinyou8. Differential expression analysis of lncRNAs indicated that up- and down-regulated mRNAs co-expressed with lncRNAs participated in different metabolic pathways and were involved in different regulatory mechanisms in the two genotypes. Notably, some lncRNAs were co-expressed with BnaC07g44670D, which are associated with plant hormone signal transduction. Additionally, some mRNAs co-located with XLOC_052298, XLOC_094954 and XLOC_012868 were mainly categorized as signal transport and defense/stress response.

CONCLUSIONS

The results of this study increased our understanding of expression characterization of rapeseed lncRNAs in response to drought stress and re-watering, which would be useful to provide a reference for the further study of the function and action mechanisms of lncRNAs under drought stress and re-watering.

Spatio-Temporal Transcriptional Dynamics of Maize Long Non-Coding RNAs Responsive to Drought Stress

Long non-coding RNAs (lncRNAs) have emerged as important regulators in plant stress response. Here, we report a genome-wide lncRNA transcriptional analysis in response to drought stress using an expanded series of maize samples collected from three distinct tissues spanning four developmental stages. In total, 3488 high-confidence lncRNAs were identified, among which 1535 were characterized as drought responsive. By characterizing the genomic structure and expression pattern, we found that lncRNA structures were less complex than protein-coding genes, showing shorter transcripts and fewer exons. Moreover, drought-responsive lncRNAs exhibited higher tissue- and development-specificity than protein-coding genes. By exploring the temporal expression patterns of drought-responsive lncRNAs at different developmental stages, we discovered that the reproductive stage R1 was the most sensitive growth stage with more lncRNAs showing altered expression upon drought stress. Furthermore, lncRNA target prediction revealed 653 potential lncRNA-messenger RNA (mRNA) pairs, among which 124 pairs function in cis-acting mode and 529 in trans. Functional enrichment analysis showed that the targets were significantly enriched in molecular functions related to oxidoreductase activity, water binding, and electron carrier activity. Multiple promising targets of drought-responsive lncRNAs were discovered, including the V-ATPase encoding gene, vpp4. These findings extend our knowledge of lncRNAs as important regulators in maize drought response.

Long Non-Coding RNAs: Rising Regulators of Plant Reproductive Development

Long non-coding RNAs (lncRNAs) exert a multitude of functions in regulating numerous biological processes. Recent studies have uncovered a growing number of lncRNAs within the plant genome. These molecules show striking tissue-specific expression patterns, suggesting that they exert regulatory functions in the growth and development processes of plants. Plant reproductive development is tightly regulated by both environmental and endogenous factors. As plant reproductive development is a crucial aspect of crop breeding, lncRNAs that modulate reproductive development are now particularly worth regarding. Here, we summarize findings that implicate lncRNAs in the control of plant reproductive development, especially in flowering control. Additionally, we elaborate on the regulation mechanisms of lncRNAs, tools for research on their function and mechanism, and potential directions of future research.

Functional identification of lncRNAs in sweet cherry (Prunus avium) pollen tubes via transcriptome analysis using single-molecule long-read sequencing

Sweet cherry (Prunus avium) is a popular fruit with high nutritional value and excellent flavor. Although pollen plays an important role in the double fertilization and subsequent fruit production of this species, little is known about its pollen tube transcriptome. In this study, we identified 16,409 transcripts using single-molecule sequencing. After filtering 292 transposable elements, we conducted further analyses including mRNA classification, gene function prediction, alternative splicing (AS) analysis, and long noncoding RNA (lncRNA) identification to gain insight into the pollen transcriptome. The filtered transcripts could be matched with 3,438 coding region sequences from the sweet cherry genome. GO and KEGG analyses revealed complex biological processes during pollen tube elongation. A total of 2043 AS events were predicted, 7 of which were identified in different organs, such as the leaf, pistil and pollen tube. Using BLASTnt and the Coding-Potential Assessment Tool (CPAT), we distinguished a total of 284 lncRNAs, among which 154 qualified as natural antisense transcripts (NATs). As the NATs could be the reverse complements of coding mRNA sequences, they might bind to coding sequences. Antisense transfection assays showed that the NATs could regulate the expression levels of their complementary sequences and even affect the growth conditions of pollen tubes. In summary, this research characterizes the transcripts of P. avium pollen and lays the foundation for elucidating the physiological and biochemical mechanisms underlying sexual reproduction in the male gametes of this species.