Non-coding RNAs and Their Roles in Stress Response in Plants

Genome-Wide Discovery of miRNAs with Differential Expression Patterns in Responses to Salinity in the Two Contrasting Wheat Cultivars

Salinity is a serious environmental issue. It has a substantial effect on crop yield, as many crop species are sensitive to salinity due to climate change, and it impact is continuing to increase. Plant microRNAs (miRNAs) contribute to salinity stress response in bread wheat. However, the underlying molecular mechanisms by which miRNAs confer salt tolerance in wheat are unclear. We conducted a genome-wide discovery study using Illumina high throughput sequencing and comprehensive in silico analysis to obtain insight into the underlying mechanisms by which small RNAs confer tolerance to salinity in roots of two contrasting wheat cvv., namely Suntop (salt-tolerant) and Sunmate (salt-sensitive). A total of 191 microRNAs were identified in both cultivars, consisting of 110 known miRNAs and 81 novel miRNAs; 181 miRNAs were shared between the two cultivars. The known miRNAs belonged to 35 families consisted of 23 conserved and 12 unique families. Salinity stress induced 43 and 75 miRNAs in Suntop and Sunmate, respectively. Among them, 14 and 29 known and novel miRNAs were expressed in Suntop and 37 and 38 in Sunmate. In silico analysis revealed 861 putative target mRNAs for the 75 known miRNAs and 52 putative target mRNAs for the 15 candidate novel miRNAs. Furthermore, seven miRNAs including tae-miR156, tae-miR160, tae-miR171a-b, tae-miR319, tae-miR159a-b, tae-miR9657 and novel-mir59 that regulate auxin responsive-factor, SPL, SCL6, PCF5, R2R3 MYB, and CBL-CIPK, respectively, were predicted to contribute to salt tolerance in Suntop. This information helps further our understanding of how the molecular mechanisms of salt tolerance are mediated by miRNAs and may facilitate the genetic improvement of wheat cultivars.

Identification and Analysis of Small Interfering RNAs Associated With Heat Stress in Flowering Chinese Cabbage Using High-Throughput Sequencing

Endogenous small interfering RNAs (siRNAs) are substantial gene regulators in eukaryotes and play key functions in plant development and stress tolerance. Among environmental factors, heat is serious abiotic stress that severely influences the productivity and quality of flowering Chinese cabbage (Brassica campestris L. ssp. chinensis var. utilis Tsen et Lee). However, how siRNAs are involved in regulating gene expression during heat stress is not fully understood in flowering Chinese cabbage. Combining bioinformatical and next-generation sequencing approaches, we identified heat-responsive siRNAs in four small RNA libraries of flowering Chinese cabbage using leaves collected at 0, 1, 6, and 12 h after a 38°C heat-stress treatment; 536, 816, and 829 siRNAs exhibited substantial differential expression at 1, 6, and 12 h, respectively. Seventy-five upregulated and 69 downregulated differentially expressed siRNAs (DE-siRNAs) were common for the three time points of heat stress. We identified 795 target genes of DE-siRNAs, including serine/threonine-protein kinase SRK2I, CTR1-like, disease resistance protein RML1A-like, and RPP1, which may play a role in regulating heat tolerance. Gene ontology showed that predictive targets of DE-siRNAs may have key roles in the positive regulation of biological processes, organismal processes, responses to temperature stimulus, signaling, and growth and development. These novel results contribute to further understanding how siRNAs modulate the expression of their target genes to control heat tolerance in flowering Chinese cabbage.

Expression analysis of miRNAs and their putative target genes confirm a preponderant role of transcription factors in the early response of oil palm plants to salinity stress

BACKGROUND

Several mechanisms regulating gene expression contribute to restore and reestablish cellular homeostasis so that plants can adapt and survive in adverse situations. MicroRNAs (miRNAs) play roles important in the transcriptional and post-transcriptional regulation of gene expression, emerging as a regulatory molecule key in the responses to plant stress, such as cold, heat, drought, and salt. This work is a comprehensive and large-scale miRNA analysis performed to characterize the miRNA population present in oil palm (Elaeis guineensis Jacq.) exposed to a high level of salt stress, to identify miRNA-putative target genes in the oil palm genome, and to perform an in silico comparison of the expression profile of the miRNAs and their putative target genes.

RESULTS

A group of 79 miRNAs was found in oil palm, been 52 known miRNAs and 27 new ones. The known miRNAs found belonged to 28 families. Those miRNAs led to 229 distinct miRNA-putative target genes identified in the genome of oil palm. miRNAs and putative target genes differentially expressed under salinity stress were then selected for functional annotation analysis. The regulation of transcription, DNA-templated, and the oxidation-reduction process were the biological processes with the highest number of hits to the putative target genes, while protein binding and DNA binding were the molecular functions with the highest number of hits. Finally, the nucleus was the cellular component with the highest number of hits. The functional annotation of the putative target genes differentially expressed under salinity stress showed several ones coding for transcription factors which have already proven able to result in tolerance to salinity stress by overexpression or knockout in other plant species.

CONCLUSIONS

Our findings provide new insights into the early response of young oil palm plants to salinity stress and confirm an expected preponderant role of transcription factors - such as NF-YA3, HOX32, and GRF1 - in this response. Besides, it points out potential salt-responsive miRNAs and miRNA-putative target genes that one can utilize to develop oil palm plants tolerant to salinity stress.

Mining plant endogenous target mimics from miRNA-lncRNA interactions based on dual-path parallel ensemble pruning method

The interactions between microRNAs (miRNAs) and long non-coding RNAs (lncRNAs) play important roles in biological activities. Specially, lncRNAs as endogenous target mimics (eTMs) can bind miRNAs to regulate the expressions of target messenger RNAs (mRNAs). A growing number of studies focus on animals, but the studies on plants are scarce and many functions of plant eTMs are unknown. This study proposes a novel ensemble pruning protocol for predicting plant miRNA-lncRNA interactions at first. It adaptively prunes the base models based on dual-path parallel ensemble method to meet the challenge of cross-species prediction. Then potential eTMs are mined from predicted results. The expression levels of RNAs are identified through biological experiment to construct the lncRNA-miRNA-mRNA regulatory network, and the functions of potential eTMs are inferred through enrichment analysis. Experiment results show that the proposed protocol outperforms existing methods and state-of-the-art predictors on various plant species. A total of 17 potential eTMs are verified by biological experiment to involve in 22 regulations, and 14 potential eTMs are inferred by Gene Ontology enrichment analysis to involve in 63 functions, which is significant for further research.

Long noncoding RNA lncRNA354 functions as a competing endogenous RNA of miR160b to regulate ARF genes in response to salt stress in upland cotton

Long non-coding RNAs (lncRNAs) play important roles in response to biotic and abiotic stress through acting as competing endogenous RNAs (ceRNAs) to decoy mature miRNAs. However, whether this mechanism is involved in cotton salt stress response remains unknown. We report the characterization of an endogenous lncRNA, lncRNA354, whose expression was reduced in salt-treated cotton and was localized at the nucleus and cytoplasm. Using endogenous target mimic (eTM) analysis, we predicted that lncRNA354 had a potential binding site for miR160b. Transient expression in tobacco demonstrated that lncRNA354 was a miR160b eTM and attenuated miR160b suppression of its target genes, including auxin response factors (ARFs). Silencing or overexpressing lncRNA354 affected the expression of miR160b and target ARFs. Silencing lncRNA354 and targets GhARF17/18 resulted in taller cotton plants and enhanced the resistant to salt stress. Overexpression of lncRNA354 and targets GhARF17/18 in Arabidopsis led to dwarf plants, decreased root dry weight and reduced salt tolerance. Our results show that the lncRNA354-miR160b effect on GhARF17/18 expression may modulate auxin signalling and thus affect growth. These results also shed new light on a mechanism of lncRNA-associated responses to salt stress.

Advances in Computational Methodologies for Classification and Sub-Cellular Locality Prediction of Non-Coding RNAs

Apart from protein-coding Ribonucleic acids (RNAs), there exists a variety of non-coding RNAs (ncRNAs) which regulate complex cellular and molecular processes. High-throughput sequencing technologies and bioinformatics approaches have largely promoted the exploration of ncRNAs which revealed their crucial roles in gene regulation, miRNA binding, protein interactions, and splicing. Furthermore, ncRNAs are involved in the development of complicated diseases like cancer. Categorization of ncRNAs is essential to understand the mechanisms of diseases and to develop effective treatments. Sub-cellular localization information of ncRNAs demystifies diverse functionalities of ncRNAs. To date, several computational methodologies have been proposed to precisely identify the class as well as sub-cellular localization patterns of RNAs). This paper discusses different types of ncRNAs, reviews computational approaches proposed in the last 10 years to distinguish coding-RNA from ncRNA, to identify sub-types of ncRNAs such as piwi-associated RNA, micro RNA, long ncRNA, and circular RNA, and to determine sub-cellular localization of distinct ncRNAs and RNAs. Furthermore, it summarizes diverse ncRNA classification and sub-cellular localization determination datasets along with benchmark performance to aid the development and evaluation of novel computational methodologies. It identifies research gaps, heterogeneity, and challenges in the development of computational approaches for RNA sequence analysis. We consider that our expert analysis will assist Artificial Intelligence researchers with knowing state-of-the-art performance, model selection for various tasks on one platform, dominantly used sequence descriptors, neural architectures, and interpreting inter-species and intra-species performance deviation.

Identification and characterization of early Fusarium wilt responsive mRNAs and long non-coding RNAs in banana root using high-throughput sequencing

Fusarium wilt disease, caused by Fusarium oxysporum f.sp. cubense (Foc), has been recognized as the most devastating disease to banana. The regulatory role of long non-coding RNAs (lncRNAs) in plant defense has been verified in many plant species. However, the understanding of their role during early FocTR4 (Foc tropical race 4) infection stage is very limited. In this study, lncRNA sequencing was used to reveal banana root transcriptome profile changes during early FocTR4 infection stages. Quantitative real time PCR (qRT-PCR) was performed to confirm the expression of eight differentially expressed (DE) lncRNAs (DELs) and their predicted target genes (DETs), and three DE genes (DEGs). Totally, 12,109 lncRNAs, 36,519 mRNAs and 2642 novel genes were obtained, of which 1398 (including 78 DELs, 1220 DE known genes and 100 DE novel genes) were identified as FocTR4 responsive DE transcripts. Gene function analysis revealed that most DEGs were involved in biosynthesis of secondary metabolites, plant–pathogen interaction, plant hormone signal transduction, phenylalanine metabolism, phenylpropanoid biosynthesis, alpha-linolenic acid metabolism and so on. Coincidently, many DETs have been identified as DEGs in previous transcriptome studies. Moreover, many DETs were found to be involved in ribosome, oxidative phosphorylation, lipoic acid metabolism, ubiquitin mediated proteolysis, N-glycan biosynthesis, protein processing in endoplasmic reticulum and DNA damage response pathways. QRT-PCR result showed the expression patterns of the selected transcripts were mostly consistent with our lncRNA sequencing data. Our present study showed the regulatory role of lncRNAs on known biotic and abiotic stress responsive genes and some new-found FocTR4 responsive genes, which can provide new insights into FocTR4-induced changes in the banana root transcriptome during the early pathogen infection stage.

Full-length transcriptome analysis of pecan (Carya illinoinensis) kernels

Pecan is rich in bioactive components such as fatty acids (FAs) and flavonoids and is an important nut type worldwide. Therefore, the molecular mechanisms of phytochemical biosynthesis in pecan are a focus of research. Recently, a draft genome and several transcriptomes have been published. However, the full-length mRNA transcripts remain unclear, and the regulatory mechanisms behind the quality components biosynthesis and accumulation have not been fully investigated. In this study, single-molecule long-read sequencing technology was used to obtain full-length transcripts of pecan kernels. In total, 37,504 isoforms of 16,702 genes were mapped to the reference genome. The numbers of known isoforms, new isoforms, and novel isoforms were 9013 (24.03%), 26,080 (69.54%), and 2411 (6.51%), respectively. Over 80% of the transcripts (30,751, 81.99%) had functional annotations. A total of 15,465 alternative splicing (AS) events and 65,761 alternative polyadenylation events were detected; wherein, the retained intron was the predominant type (5652, 36.55%) of AS. Furthermore, 1894 long noncoding RNAs and 1643 transcription factors were predicted using bioinformatics methods. Finally, the structural genes associated with FA and flavonoid biosynthesis were characterized. A high frequency of AS accuracy (70.31%) was observed in FA synthesis-associated genes. This study provides a full-length transcriptome data set of pecan kernels, which will significantly enhance the understanding of the regulatory basis of phytochemical biosynthesis during pecan kernel maturation.

Biochemical and Molecular Aspects of Phosphorus Limitation in Diatoms and Their Relationship with Biomolecule Accumulation

Simple Summary

Phosphorus (P) is a key nutrient involved in the transfer of energy and the synthesis of several cellular components. It has been reported that P limitation in diatoms induces the synthesis of biomolecules and the accumulation of storage compounds, such as pigments, carbohydrates and lipids, with diverse biological activities, which can be used in diverse biotechnological applications. However, the molecular and biochemical mechanisms related to how diatoms cope with P deficiency are not clear, and research into this has been limited to a few species. The integration of results obtained from omics sciences could provide a broad understanding of the response of diatoms to P limitation, and the information obtained could help to solve challenges such as biomass production, by-products yield and genetic improvement of strains.

Abstract

Diatoms are the most abundant group of phytoplankton, and their success lies in their significant adaptation ability to stress conditions, such as nutrient limitation. Phosphorus (P) is a key nutrient involved in the transfer of energy and the synthesis of several cellular components. Molecular and biochemical mechanisms related to how diatoms cope with P deficiency are not clear, and research into this has been limited to a few species. Among the molecular responses that have been reported in diatoms cultured under P deficient conditions is the upregulation of genes encoding enzymes related to the transport, assimilation, remobilization and recycling of this nutrient. Regarding biochemical responses, due to the reduction of the requirements for carbon structures for the synthesis of proteins and phospholipids, more CO₂ is fixed than is consumed by the Calvin cycle. To deal with this excess, diatoms redirect the carbon flow toward the synthesis of storage compounds such as triacylglycerides and carbohydrates, which are excreted as extracellular polymeric substances. This review aimed to gather all current knowledge regarding the biochemical and molecular mechanisms of diatoms related to managing P deficiency in order to provide a wider insight into and understanding of their responses, as well as the metabolic pathways affected by the limitation of this nutrient.

Advances in Cereal Crop Genomics for Resilience under Climate Change

Adapting to climate change, providing sufficient human food and nutritional needs, and securing sufficient energy supplies will call for a radical transformation from the current conventional adaptation approaches to more broad-based and transformative alternatives. This entails diversifying the agricultural system and boosting productivity of major cereal crops through development of climate-resilient cultivars that can sustainably maintain higher yields under climate change conditions, expanding our focus to crop wild relatives, and better exploitation of underutilized crop species. This is facilitated by the recent developments in plant genomics, such as advances in genome sequencing, assembly, and annotation, as well as gene editing technologies, which have increased the availability of high-quality reference genomes for various model and non-model plant species. This has necessitated genomics-assisted breeding of crops, including underutilized species, consequently broadening genetic variation of the available germplasm; improving the discovery of novel alleles controlling important agronomic traits; and enhancing creation of new crop cultivars with improved tolerance to biotic and abiotic stresses and superior nutritive quality. Here, therefore, we summarize these recent developments in plant genomics and their application, with particular reference to cereal crops (including underutilized species). Particularly, we discuss genome sequencing approaches, quantitative trait loci (QTL) mapping and genome-wide association (GWAS) studies, directed mutagenesis, plant non-coding RNAs, precise gene editing technologies such as CRISPR-Cas9, and complementation of crop genotyping by crop phenotyping. We then conclude by providing an outlook that, as we step into the future, high-throughput phenotyping, pan-genomics, transposable elements analysis, and machine learning hold much promise for crop improvements related to climate resilience and nutritional superiority.

PlncRNA-HDeep: plant long noncoding RNA prediction using hybrid deep learning based on two encoding styles

Background

Long noncoding RNAs (lncRNAs) play an important role in regulating biological activities and their prediction is significant for exploring biological processes. Long short-term memory (LSTM) and convolutional neural network (CNN) can automatically extract and learn the abstract information from the encoded RNA sequences to avoid complex feature engineering. An ensemble model learns the information from multiple perspectives and shows better performance than a single model. It is feasible and interesting that the RNA sequence is considered as sentence and image to train LSTM and CNN respectively, and then the trained models are hybridized to predict lncRNAs. Up to present, there are various predictors for lncRNAs, but few of them are proposed for plant. A reliable and powerful predictor for plant lncRNAs is necessary.

Results

To boost the performance of predicting lncRNAs, this paper proposes a hybrid deep learning model based on two encoding styles (PlncRNA-HDeep), which does not require prior knowledge and only uses RNA sequences to train the models for predicting plant lncRNAs. It not only learns the diversified information from RNA sequences encoded by p-nucleotide and one-hot encodings, but also takes advantages of lncRNA-LSTM proposed in our previous study and CNN. The parameters are adjusted and three hybrid strategies are tested to maximize its performance. Experiment results show that PlncRNA-HDeep is more effective than lncRNA-LSTM and CNN and obtains 97.9% sensitivity, 95.1% precision, 96.5% accuracy and 96.5% F1 score on Zea mays dataset which are better than those of several shallow machine learning methods (support vector machine, random forest, k-nearest neighbor, decision tree, naive Bayes and logistic regression) and some existing tools (CNCI, PLEK, CPC2, LncADeep and lncRNAnet).

Conclusions

PlncRNA-HDeep is feasible and obtains the credible predictive results. It may also provide valuable references for other related research.

Role of non-coding RNAs in plant immunity

Crops are exposed to attacks by various pathogens that cause substantial yield losses and severely threaten food security. To cope with pathogenic infection, crops have elaborated strategies to enhance resistance against pathogens. In addition to the role of protein-coding genes as key regulators in plant immunity, accumulating evidence has demonstrated the importance of non-coding RNAs (ncRNAs) in the plant immune response. Here, we summarize the roles and molecular mechanisms of endogenous ncRNAs, especially microRNAs (miRNAs), long ncRNAs (lncRNAs), and circular RNAs (circRNAs), in plant immunity. We discuss the coordination between miRNAs and small interfering RNAs (siRNAs), between lncRNAs and miRNAs or siRNAs, and between circRNAs and miRNAs in the regulation of plant immune responses. We also address the role of cross-kingdom mobile small RNAs in plant-pathogen interactions. These insights improve our understanding of the mechanisms by which ncRNAs regulate plant immunity and can promote the development of better approaches for breeding disease-resistant crops.

Ensemble Deep Learning Based on Multi-level Information Enhancement and Greedy Fuzzy Decision for Plant miRNA-lncRNA Interaction Prediction

MicroRNAs (miRNAs) and long non-coding RNAs (lncRNAs) are both non-coding RNAs (ncRNAs) and their interactions play important roles in biological processes. Computational methods, such as machine learning and various bioinformatics tools, can predict potential miRNA-lncRNA interactions, which is significant for studying their mechanisms and biological functions. A growing number of RNA interaction predictors for animal have been reported, but they are unreliable for plant due to the differences of ncRNAs in animal and plant. It is urgent to build a reliable plant predictor, especially for cross-species. This paper proposes an ensemble deep learning model based on multi-level information enhancement and greedy fuzzy decision (PmliPEMG) for plant miRNA-lncRNA interaction prediction. The fusion complex features, multi-scale convolutional long short-term memory networks, and attention mechanism are adopted to enhance the sample information at the feature, scale, and model levels, respectively. An ensemble deep learning model is built based on a novel method (greedy fuzzy decision) which greatly improves the efficiency. The multi-level information enhancement and greedy fuzzy decision are verified to have the positive effects on prediction performance. PmliPEMG can be applied to the cross-species prediction. It shows better performance and stronger generalization ability than state-of-the-art predictors and may provide valuable references for related research.

RNA sequencing-based exploration of the effects of far-red light on lncRNAs involved in the shade-avoidance response of D. officinale

Dendrobium officinale (D. officinale) is a valuable medicinal plant with a low natural survival rate, and its shade-avoidance response to far-red light is as an important strategy used by the plant to improve its production efficiency. However, the lncRNAs that play roles in the shade-avoidance response of D. officinale have not yet been investigated. This study found that an appropriate proportion of far-red light can have several effects, including increasing the leaf area and accelerating stem elongation, in D. officinale. The effects of different far-red light treatments on D. officinale were analysed by RNA sequencing technology, and a total of 69 and 78 lncRNAs were differentially expressed in experimental group 1 (FR1) versus the control group (CK) (FR1-CK) and in experimental group 4 (FR4) versus the CK (FR4-CK), respectively. According to GO and KEGG analyses, most of the differentially expressed lncRNA targets are involved in the membrane, some metabolic pathways, hormone signal transduction, and O-methyltransferase activity, among other functions. Physiological and biochemical analyses showed that far-red light promoted the accumulation of flavonoids, alkaloids, carotenoids and polysaccharides in D. officinale. The effect of far-red light on D. officinalemight be closely related to the cell membrane and Ca2+ transduction. Based on a Cytoscape analysis and previous research, this study also found that MSTRG.38867.1, MSTRG.69319.1, and MSTRG.66273.1, among other components, might participate in the far-red light signalling network through their targets and thus regulate the shade-avoidance response of D. officinale. These findings will provide new insights into the shade-avoidance response of D. officinale.

Contribution of Small RNA Pathway to Interactions of Rice with Pathogens and Insect Pests

Small RNAs (sRNAs) are mainly classified into microRNAs (miRNAs) and small interfering RNAs (siRNAs) according to their origin. miRNAs originate from single-stranded RNA precursors, whereas siRNAs originate from double-stranded RNA precursors that are synthesized by RNA-dependent RNA polymerases. Both of single-stranded and double-stranded RNA precursors are processed into sRNAs by Dicer-like proteins. Then, the sRNAs are loaded into ARGONAUTE proteins, forming RNA-induced silencing complexes (RISCs). The RISCs repress the expression of target genes with sequences complementary to the sRNAs through the cleavage of transcripts, the inhibition of translation or DNA methylation. Here, we summarize the recent progress of sRNA pathway in the interactions of rice with various parasitic organisms, including fungi, viruses, bacteria, as well as insects. Besides, we also discuss the hormone signal in sRNA pathway, and the emerging roles of circular RNAs and long non-coding RNAs in rice immunity. Obviously, small RNA pathway may act as a part of rice innate immunity to coordinate with growth and development.

Whole-Genome DNA Methylation Analysis in Hydrogen Peroxide Overproducing Transgenic Tobacco Resistant to Biotic and Abiotic Stresses

Epigenetic regulation is a key component of stress responses, acclimatization and adaptation processes in plants. DNA methylation is a stable mark plausible for the inheritance of epigenetic traits, such that it is a potential scheme for plant breeding. However, the effect of modulators of stress responses, as hydrogen peroxide (H₂O₂), in the methylome status has not been elucidated. A transgenic tobacco model to the CchGLP gene displayed high H₂O₂ endogen levels correlated with biotic and abiotic stresses resistance. The present study aimed to determine the DNA methylation status changes in the transgenic model to obtain more information about the molecular mechanism involved in resistance phenotypes. The Whole-genome bisulfite sequencing analysis revealed a minimal impact of overall levels and distribution of methylation. A total of 9432 differential methylated sites were identified in distinct genome regions, most of them in CHG context, with a trend to hypomethylation. Of these, 1117 sites corresponded to genes, from which 83 were also differentially expressed in the plants. Several genes were associated with respiration, energy, and calcium signaling. The data obtained highlighted the relevance of the H₂O₂ in the homeostasis of the system in stress conditions, affecting at methylation level and suggesting an association of the H₂O₂ in the physiological adaptation to stress functional linkages may be regulated in part by DNA methylation.

Long Non-Coding RNAs, the Dark Matter: An Emerging Regulatory Component in Plants

Long non-coding RNAs (lncRNAs) are pervasive transcripts of longer than 200 nucleotides and indiscernible coding potential. lncRNAs are implicated as key regulatory molecules in various fundamental biological processes at transcriptional, post-transcriptional, and epigenetic levels. Advances in computational and experimental approaches have identified numerous lncRNAs in plants. lncRNAs have been found to act as prime mediators in plant growth, development, and tolerance to stresses. This review summarizes the current research status of lncRNAs in planta, their classification based on genomic context, their mechanism of action, and specific bioinformatics tools and resources for their identification and characterization. Our overarching goal is to summarize recent progress on understanding the regulatory role of lncRNAs in plant developmental processes such as flowering time, reproductive growth, and abiotic stresses. We also review the role of lncRNA in nutrient stress and the ability to improve biotic stress tolerance in plants. Given the pivotal role of lncRNAs in various biological processes, their functional characterization in agriculturally essential crop plants is crucial for bridging the gap between phenotype and genotype.

RNAi Mediated Hypoxia Stress Tolerance in Plants

Small RNAs regulate various biological process involved in genome stability, development, and adaptive responses to biotic or abiotic stresses. Small RNAs include microRNAs (miRNAs) and small interfering RNAs (siRNAs). MicroRNAs (miRNAs) are regulators of gene expression that affect the transcriptional and post-transcriptional regulation in plants and animals through RNA interference (RNAi). miRNAs are endogenous small RNAs that originate from the processing of non-coding primary miRNA transcripts folding into hairpin-like structures. The mature miRNAs are incorporated into the RNA-induced silencing complex (RISC) and drive the Argonaute (AGO) proteins towards their mRNA targets. siRNAs are generated from a double-stranded RNA (dsRNA) of cellular or exogenous origin. siRNAs are also involved in the adaptive response to biotic or abiotic stresses. The response of plants to hypoxia includes a genome-wide transcription reprogramming. However, little is known about the involvement of RNA signaling in gene regulation under low oxygen availability. Interestingly, miRNAs have been shown to play a role in the responses to hypoxia in animals, and recent evidence suggests that hypoxia modulates the expression of various miRNAs in plant systems. In this review, we describe recent discoveries on the impact of RNAi on plant responses to hypoxic stress in plants.

Interactions and links among the noncoding RNAs in plants under stresses

The complex interplay among sRNAs, lncRNAs and circRNAs has been implicated in plants under biotic and abiotic stresses. Here, we review current advances in our understanding of ncRNA interactions and links, which have considerable potential for improving the agronomic traits and the environmental adaptability of plants. Plants can respond to biotic or abiotic stresses. To cope with various conditions, numerous intricate molecular regulatory mechanisms have evolved in plants. Noncoding RNAs (ncRNAs) can be divided into small noncoding RNAs (sRNAs), long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs). Emerging evidence has demonstrated that interplay among the ncRNAs acts as a novel layer in the regulatory mechanisms, which has attracted substantial interest. Links between sRNAs can affect plant immune responses and development in synergistic or antagonistic manners. Additionally, multiple interactions between lncRNAs and sRNAs are involved in crop breeding, disease resistance and high tolerance to environmental stresses. Here, we summarize current knowledge of the interactions and links among the ncRNAs in plant responses to stresses and the methods for identifying ncRNA interactions. Furthermore, challenges and prospects for further progress in elucidating ncRNA interactions and links are highlighted.

Plant Non-Coding RNAs: Origin, Biogenesis, Mode of Action and Their Roles in Abiotic Stress

As sessile species, plants have to deal with the rapidly changing environment. In response to these environmental conditions, plants employ a plethora of response mechanisms that provide broad phenotypic plasticity to allow the fine-tuning of the external cues related reactions. Molecular biology has been transformed by the major breakthroughs in high-throughput transcriptome sequencing and expression analysis using next-generation sequencing (NGS) technologies. These innovations have provided substantial progress in the identification of genomic regions as well as underlying basis influencing transcriptional and post-transcriptional regulation of abiotic stress response. Non-coding RNAs (ncRNAs), particularly microRNAs (miRNAs), short interfering RNAs (siRNAs), and long non-coding RNAs (lncRNAs), have emerged as essential regulators of plants abiotic stress response. However, shared traits in the biogenesis of ncRNAs and the coordinated cross-talk among ncRNAs mechanisms contribute to the complexity of these molecules and might play an essential part in regulating stress responses. Herein, we highlight the current knowledge of plant microRNAs, siRNAs, and lncRNAs, focusing on their origin, biogenesis, modes of action, and fundamental roles in plant response to abiotic stresses.

Genome-wide analysis of long non-coding RNAs responsive to multiple nutrient stresses in Arabidopsis thaliana

Nutrient stress is the most important environmental stress that limits plant growth and development. Although recent evidence highlights the vital functions of long non-coding RNAs (lncRNA) in response to single nutrient stress in some model plants, a comprehensive investigation of the effect of lncRNAs in response to nutrient stress has not been performed in Arabidopsis thaliana. Here, we presented the identification and characterization of lncRNAs under seven nutrient stress conditions. The expression pattern analysis revealed that aberrant expression of lncRNAs is a stress-specific manner under nutrient stress conditions and that lncRNAs are more sensitive to nutrient stress than protein-coding genes (PCGs). Moreover, competing endogenous RNA (ceRNA) network and lncRNA-mRNA co-expression network (CEN) were constructed to explore the potential function of these lncRNAs under nutrient stress conditions. We further combined different expressed lncRNAs with ceRNA network and CEN to select key lncRNAs in response to nutrient stress. Together, our study provides important information for further insights into the role of lncRNAs in response to stress in plants.

Genome Size Evolution Mediated by Gypsy Retrotransposons in Brassicaceae

The dynamic activity of transposable elements (TEs) contributes to the vast diversity of genome size and architecture among plants. Here, we examined the genomic distribution and transposition activity of long terminal repeat retrotransposons (LTR-RTs) in Arabidopsis thaliana (Ath) and three of its relatives, Arabidopsis lyrata (Aly), Eutrema salsugineum (Esa), and Schrenkiella parvula (Spa), in Brassicaceae. Our analyses revealed the distinct evolutionary dynamics of Gypsyretrotransposons, which reflects the different patterns of genome size changes of the four species over the past million years. The rate of Gypsy transposition in Aly is approximately five times more rapid than that of Ath and Esa, suggesting an expanding Aly genome. Gypsy insertions in Esa are strictly confined to pericentromeric heterochromatin and associated with dramatic centromere expansion. In contrast, Gypsy insertions in Spa have been largely suppressed over the last million years, likely as a result of a combination of an inherent molecular mechanism of preferential DNA removal and purifying selection at Gypsy elements. Additionally, species-specific clades of Gypsy elements shaped the distinct genome architectures of Aly and Esa.

PmliPred: a method based on hybrid model and fuzzy decision for plant miRNA-lncRNA interaction prediction

MOTIVATION

The studies have indicated that not only microRNAs (miRNAs) or long non-coding RNAs (lncRNAs) play important roles in biological activities, but also their interactions affect the biological process. A growing number of studies focus on the miRNA-lncRNA interactions, while few of them are proposed for plant. The prediction of interactions is significant for understanding the mechanism of interaction between miRNA and lncRNA in plant.

RESULTS

This article proposes a new method for fulfilling plant miRNA-lncRNA interaction prediction (PmliPred). The deep learning model and shallow machine learning model are trained using raw sequence and manually extracted features, respectively. Then they are hybridized based on fuzzy decision for prediction. PmliPred shows better performance and generalization ability compared with the existing methods. Several new miRNA-lncRNA interactions in Solanum lycopersicum are successfully identified using quantitative real time-polymerase chain reaction from the candidates predicted by PmliPred, which further verifies its effectiveness.

AVAILABILITY AND IMPLEMENTATION

The source code of PmliPred is freely available at http://bis.zju.edu.cn/PmliPred/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Long non-coding RNAs in plants: emerging modulators of gene activity in development and stress responses

MAIN CONCLUSION

Long non-coding RNAs modulate gene activity in plant development and stress responses by various molecular mechanisms. Long non-coding RNAs (lncRNAs) are transcripts larger than 200 nucleotides without protein coding potential. Computational approaches have identified numerous lncRNAs in different plant species. Research in the past decade has unveiled that plant lncRNAs participate in a wide range of biological processes, including regulation of flowering time and morphogenesis of reproductive organs, as well as abiotic and biotic stress responses. LncRNAs execute their functions by interacting with DNA, RNA and protein molecules, and by modulating the expression level of their targets through epigenetic, transcriptional, post-transcriptional or translational regulation. In this review, we summarize characteristics of plant lncRNAs, discuss recent progress on understanding of lncRNA functions, and propose an experimental framework for functional characterization.

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

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

Versatile interactions and bioinformatics analysis of noncoding RNAs

Advances in RNA sequencing technologies and computational methodologies have provided a huge impetus to noncoding RNA (ncRNA) study. Once regarded as inconsequential results of transcriptional promiscuity, ncRNAs were later found to exert great roles in various aspects of biological functions. They are emerging as key players in gene regulatory networks by interacting with other biomolecules (DNA, RNA or protein). Here, we provide an overview of ncRNA repertoire and highlight recent discoveries of their versatile interactions. To better investigate the ncRNA-mediated regulation, it is necessary to make full use of innovative sequencing techniques and computational tools. We further describe a comprehensive workflow for in silico ncRNA analysis, providing up-to-date platforms, databases and tools dedicated to ncRNA identification and functional annotation.

Systematic Identification and Analysis of Light-Responsive Circular RNA and Co-expression Networks in Lettuce (Lactuca sativa)

Circular RNA (circRNA) is a covalently-closed single-stranded RNA molecule that plays an important role in transcriptional regulation of gene expression in a variety of species. Light intensity is a pivotal environmental factor affecting plant growth and development. However, little is known regarding photoresponsive plant circRNAs. Here, we aimed to investigate the expression and function of circRNAs in lettuce leaves in response to different light intensity treatments. We performed RNA sequencing (RNA-Seq) on leaves of lettuce (Lactuca sativa) to determine circRNA expression profiles and reverse-transcription polymerase chain reaction (PCR) to validate the candidate circRNA molecules. We then combined bioinformatics approach to explore the function of the parental genes of circRNA, including network, Gene Ontology, and Kyoto Encyclopedia of Genes and Genomes pathway analysis. We identified 1650 circRNAs in lettuce, of which 1508 (86.40%) were derived from exons. Using real-time PCR, we characterized 10 validated differentially expressed circRNAs and their parental genes, all of which showed expression patterns consistent with RNA-Seq data. Interestingly, the expression of circRNA was, in some cases, inversely correlated with the expression of the parental gene. Furthermore, analysis of the circRNA-microRNA-mRNA network suggests that circRNAs may be involved in plant hormone signaling and chlorophyll metabolism during photoreactivity. These findings provide an essential reference basis for studying circRNAs' biological mechanisms in light-treated plants.

Characterization and Analysis of the Full-Length Transcriptomes of Multiple Organs in Pseudotaxus chienii (W.C.Cheng) W.C.Cheng

Pseudotaxus chienii, a rare tertiary relict species with economic and ecological value, is a representative of the monotypic genus Pseudotaxus that is endemic to China. P. chienii can adapt well to habitat isolation and ecological heterogeneity under a variety of climate and soil conditions, and is able to survive in harsh environments. However, little is known about the molecular and genetic resources of this long-lived conifer. Herein, we sequenced the transcriptomes of four organs of P. chienii using the PacBio Isoform Sequencing and Illumina RNA Sequencing platforms. Based on the PacBio Iso-Seq data, we obtained 44,896, 58,082, 50,485, and 67,638 full-length unigenes from the root, stem, leaf, and strobilus, respectively, with a mean length of 2692 bp, and a mean N50 length of 3010.75 bp. We then comprehensively annotated these unigenes. The number of organ-specific expressed unigenes ranged from 4393 in leaf to 9124 in strobilus, suggesting their special roles in physiological processes, organ development, and adaptability in the different four organs. A total of 16,562 differentially expressed genes (DEGs) were identified among the four organs and clustered into six subclusters. The gene families related to biotic/abiotic factors, including the TPS, CYP450, and HSP families, were characterized. The expression levels of most DEGs in the phenylpropanoid biosynthesis pathway and plant–pathogen interactions were higher in the root than in the three other organs, suggesting that root constitutes the main organ of defensive compound synthesis and accumulation and has a stronger ability to respond to stress. The sequences were analyzed to predict transcription factors, long non-coding RNAs, and alternative splicing events. The expression levels of most DEGs of C2H2, C3H, bHLH, and bZIP families in the root and stem were higher than those in the leaf and strobilus, indicating that these TFs may play a crucial role in the survival of the root and stem. These results comprise the first comprehensive gene expression profiles obtained for different organs of P. chienii. Our findings will facilitate further studies on the functional genomics, adaptive evolution, and phylogeny of P. chienii, and lay the foundation for the development of conservation strategies for this endangered conifer.

Mechanisms of Plant Responses and Adaptation to Soil Salinity

Soil salinity is a major environmental stress that restricts the growth and yield of crops. Understanding the physiological, metabolic, and biochemical responses of plants to salt stress and mining the salt tolerance-associated genetic resource in nature will be extremely important for us to cultivate salt-tolerant crops. In this review, we provide a comprehensive summary of the mechanisms of salt stress responses in plants, including salt stress-triggered physiological responses, oxidative stress, salt stress sensing and signaling pathways, organellar stress, ion homeostasis, hormonal and gene expression regulation, metabolic changes, as well as salt tolerance mechanisms in halophytes. Important questions regarding salt tolerance that need to be addressed in the future are discussed.

Epigenomic landscape and epigenetic regulation in maize

Epigenetic regulation has been implicated in the control of multiple agronomic traits in maize. Here, we review current advances in our understanding of epigenetic regulation, which has great potential for improving agronomic traits and the environmental adaptability of crops. Epigenetic regulation plays vital role in the control of complex agronomic traits. Epigenetic variation could contribute to phenotypic diversity and can be used to improve the quality and productivity of crops. Maize (Zea mays L.), one of the most widely cultivated crops for human food, animal feed, and ethanol biofuel, is a model plant for genetic studies. Recent advances in high-throughput sequencing technology have made possible the study of epigenetic regulation in maize on a genome-wide scale. In this review, we discuss recent epigenetic studies in maize many achieved by Chinese research groups. These studies have explored the roles of DNA methylation, posttranslational modifications of histones, chromatin remodeling, and noncoding RNAs in the regulation of gene expression in plant development and environment response. We also provide our future prospects for manipulating epigenetic regulation to improve crops.

Distinct non-coding RNAs confer root-dependent sense transgene-induced post-transcriptional gene silencing and nitrogen-dependent post-transcriptional regulation to AtAMT1;1 transcripts in Arabidopsis roots

High-affinity ammonium uptake in roots mediate by AMT1-type ammonium transporters, which are tightly controlled at multiple regulatory levels for adapting various nitrogen availability. For Arabidopsis AtAMT1;1 gene, in addition to the transcriptional and post-translational controls, an organ-dependent and N-dependent post-transcriptional regulation was suggested as an additional regulatory step for fine tuning ammonium uptake, but the underlying mechanisms remain to be elucidated. Here, we showed that degradation of AtAMT1;1 transcript in roots of Pro35s:AtAMT1;1-transformed atamt1;1-1 Arabidopsis plants resulted from RDR6-dependent sense transgene-induced post-transcriptional gene silencing (S-PTGS). The siRNAs for S-PTGS may derive from the aberrant RNA, of which the production was co-determined by sequence feature and excessive expression of AtAMT1;1. Switching to the expression of AtAMT1;1 driven by ProAtUBQ10 or of AtAMT1;1 mutated at two siRNA-targeted hotspots reduced AtAMT1;1-specific siRNAs and overcame S-PTGS in roots. In roots of these lines, however, the steady-state transcript levels of AtAMT1;1 still significantly decreased under conditions of N-sufficiency compared with N-deficiency, confirming a N-dependent post-transcriptional regulatory manner. A crucial role of the 207-bp 3'-end sequence of AtAMT1;1 was further demonstrated by N-dependent accumulation of chimeric-AtAMT1;1 transcript in T-DNA insertion lines and of GFP-tagged chimeric-AtAMT1;1 transcript in transgenic lines. A novel non-coding RNA (ncRNA), which was highly abundant in N-sufficient roots, may target the above-identified 3'-end region for the degrading AtAMT1;1 transcript. This degradation could be prevented by a mutation on the AtAMT1;1 transcript at a potential cleavage site (+1458). These results suggested two distinct mechanisms of regulating AtAMT1;1 mRNA turnover by ncRNA for strictly control of ammonium uptake in roots.

A point mutation in LTT1 enhances cold tolerance at the booting stage in rice

The cold tolerance of rice at the booting stage is a main factor determining sustainability and regional adaptability. However, relatively few cold tolerance genes have been identified that can be effectively used in breeding programmes. Here, we show that a point mutation in the low-temperature tolerance 1 (LTT1) gene improves cold tolerance by maintaining tapetum degradation and pollen development, by activation of systems that metabolize reactive oxygen species (ROS). Cold-induced ROS accumulation is therefore prevented in the anthers of the ltt1 mutants allowing correct development. In contrast, exposure to cold stress dramatically increases ROS accumulation in the wild type anthers, together with the expression of genes encoding proteins associated with programmed cell death and with the accelerated degradation of the tapetum that ultimately leads to pollen abortion. These results demonstrate that appropriate ROS management is critical for the cold tolerance of rice at the booting stage. Hence, the ltt1 mutation can significantly improve the seed setting ability of cold-sensitive rice varieties under low-temperature stress conditions, with little yield penalty under optimal temperature conditions. This study highlights the importance of a valuable genetic resource that may be applied in rice breeding programmes to enhance cold tolerance.

A genome-wide identification, characterization and functional analysis of salt-related long non-coding RNAs in non-model plant Pistacia vera L. using transcriptome high throughput sequencing

Long non-coding RNAs (lncRNAs) play crucial roles in regulating gene expression in response to plant stresses. Given the importance regulatory roles of lncRNAs, providing methods for predicting the function of these molecules, especially in non-model plants, is strongly demanded by researchers. Here, we constructed a reference sequence for lncRNAs in P. vera (Pistacia vera L.) with 53220 transcripts. In total, we identified 1909 and 2802 salt responsive lncRNAs in Ghazvini, a salt tolerant cultivar, after 6 and 24 h salt treatment, respectively and 1820 lncRNAs in Sarakhs, a salt sensitive cultivar, after 6 h salt treatment. Functional analysis of these lncRNAs by several hybrid methods, revealed that salt responsive NAT-related lncRNAs associated with transcription factors, CERK1, LEA, Laccase genes and several genes involved in the hormone signaling pathways. Moreover, gene ontology (GO) enrichment analysis of salt responsive target genes related to top five selected lncRNAs showed their involvement in the regulation of ATPase, cation transporter, kinase and UDP-glycosyltransferases genes. Quantitative real-time PCR (qRT-PCR) experiment results of lncRNAs, pre-miRNAs and mature miRNAs were in accordance with our RNA-seq analysis. In the present study, a comparative analysis of differentially expressed lncRNAs and microRNA precursors between salt tolerant and sensitive pistachio cultivars provides valuable knowledge on gene expression regulation under salt stress condition.

A Comprehensive Transcriptomics Analysis Reveals Long Non-Coding RNA to be Involved in the Key Metabolic Pathway in Response to Waterlogging Stress in Maize

Waterlogging stress (WS) in a dynamic environment seriously limits plant growth, development, and yield. The regulatory mechanism underlying WS conditions at an early stage in maize seedlings is largely unknown. In the present study, the primary root tips of B73 seedlings were sampled before (0 h) and after (2 h, 4 h, 6 h, 8 h, 10 h, and 12 h) WS and then subjected to transcriptome sequencing, resulting in the identification of differentially expressed protein-coding genes (DEpcGs) and long non-coding RNAs (DElncRs) in response to WS. These DEpcGs were classified into nine clusters, which were significantly enriched in several metabolic pathways, such as glycolysis and methionine metabolism. Several transcription factor families, including AP2-EREBP, bZIP, NAC, bHLH, and MYB, were also significantly enriched. In total, 6099 lncRNAs were identified, of which 3190 were DElncRs. A co-expression analysis revealed lncRNAs to be involved in 11 transcription modules, 10 of which were significantly associated with WS. The DEpcGs in the four modules were enriched in the hypoxia response pathways, including phenylpropanoid biosynthesis, MAPK signaling, and carotenoid biosynthesis, in which 137 DElncRs were also co-expressed. Most of the co-expressed DElncRs were co-localized with previously identified quantitative trait loci associated with waterlogging tolerance. A quantitative reverse transcription-polymerase chain reaction analysis of DEpcG and DElncR expression among the 32 maize genotypes after 4 h of WS verified significant expression correlations between them as well as significant correlation with the phenotype of waterlogging tolerance. Moreover, the high proportion of hypoxia response elements in the promoter region increased the reliability of the DElncRs identified in this study. These results provide a comprehensive transcriptome in response to WS at an early stage of maize seedlings and expand our understanding of the regulatory network involved in hypoxia in plants.

The Arabidopsis Hypoxia Inducible AtR8 Long Non-Coding RNA also Contributes to Plant Defense and Root Elongation Coordinating with WRKY Genes under Low Levels of Salicylic Acid

AtR8 lncRNA was previously identified in the flowering plant Arabidopsis thaliana as an abundant Pol III-transcribed long non-coding RNA (lncRNA) of approximately 260 nt. AtR8 lncRNA accumulation is responsive to hypoxic stress and salicylic acid (SA) treatment in roots, but its function has not yet been identified. In this study, microarray analysis of an atr8 mutant and wild-type Arabidopsis indicated a strong association of AtR8 lncRNA with the defense response. AtR8 accumulation exhibited an inverse correlation with an accumulation of two WRKY genes (WRKY53/WRKY70) when plants were exposed to exogenous low SA concentrations (20 µM), infected with Pseudomonas syringae, or in the early stage of development. The highest AtR8 accumulation was observed 5 days after germination, at which time no WRKY53 or WRKY70 mRNA was detectable. The presence of low levels of SA resulted in a significant reduction of root length in atr8 seedlings, whereas wrky53 and wrky70 mutants exhibited the opposite phenotype. Taken together, AtR8 lncRNA participates in Pathogenesis-Related Proteins 1 (PR-1)-independent defense and root elongation, which are related to the SA response. The mutual regulation of AtR8 lncRNA and WRKY53/WRKY70 is mediated by Nonexpressor of Pathogenesis-Related Gene 1 (NPR1).

Insight into the Role of Epigenetic Processes in Abiotic and Biotic Stress Response in Wheat and Barley

Environmental stresses such as salinity, drought, heat, freezing, heavy metal and even pathogen infections seriously threaten the growth and yield of important cereal crops including wheat and barley. There is growing evidence indicating that plants employ sophisticated epigenetic mechanisms to fine-tune their responses to environmental stresses. Here, we provide an overview of recent developments in understanding the epigenetic processes and elements—such as DNA methylation, histone modification, chromatin remodeling, and non-coding RNAs—involved in plant responses to abiotic and biotic stresses in wheat and barley. Potentials of exploiting epigenetic variation for the improvement of wheat and barley are discussed.

Machine learning improves our knowledge about miRNA functions towards plant abiotic stresses

During the last two decades, human has increased his knowledge about the role of miRNAs and their target genes in plant stress response. Biotic and abiotic stresses result in simultaneous tissue-specific up/down-regulation of several miRNAs. In this study, for the first time, feature selection algorithms have been used to investigate the contribution of individual plant miRNAs in Arabidopsis thaliana response towards different levels of several abiotic stresses including drought, salinity, cold, and heat. Results of information theory-based feature selection revealed that miRNA-169, miRNA-159, miRNA-396, and miRNA-393 had the highest contributions to plant response towards drought, salinity, cold, and heat, respectively. Furthermore, regression models, i.e., decision tree (DT), support vector machines (SVMs), and Naïve Bayes (NB) were used to predict the plant stress by having the plant miRNAs' concentration. SVM with Gaussian kernel was capable of predicting plant stress (R2 = 0.96) considering miRNA concentrations as input features. Findings of this study prove the performance of machine learning as a promising tool to investigate some aspects of miRNAs' contribution to plant stress responses that have been undiscovered until today.

Identification of a novel anthocyanin synthesis pathway in the fungus Aspergillus sydowii H-1

BACKGROUND

Anthocyanins are common substances with many agro-food industrial applications. However, anthocyanins are generally considered to be found only in natural plants. Our previous study isolated and purified the fungus Aspergillus sydowii H-1, which can produce purple pigments during fermentation. To understand the characteristics of this strain, a transcriptomic and metabolomic comparative analysis was performed with A. sydowii H-1 from the second and eighth days of fermentation, which confer different pigment production.

RESULTS

We found five anthocyanins with remarkably different production in A. sydowii H-1 on the eighth day of fermentation compared to the second day of fermentation. LC-MS/MS combined with other characteristics of anthocyanins suggested that the purple pigment contained anthocyanins. A total of 28 transcripts related to the anthocyanin biosynthesis pathway was identified in A. sydowii H-1, and almost all of the identified genes displayed high correlations with the metabolome. Among them, the chalcone synthase gene (CHS) and cinnamate-4-hydroxylase gene (C4H) were only found using the de novo assembly method. Interestingly, the best hits of these two genes belonged to plant species. Finally, we also identified 530 lncRNAs in our datasets, and among them, three lncRNAs targeted the genes related to anthocyanin biosynthesis via cis-regulation, which provided clues for understanding the underlying mechanism of anthocyanin production in fungi.

CONCLUSION

We first reported that anthocyanin can be produced in fungus, A. sydowii H-1. Totally, 31 candidate transcripts were identified involved in anthocyanin biosynthesis, in which CHS and C4H, known as the key genes in anthocyanin biosynthesis, were only found in strain H1, which indicated that these two genes may contribute to anthocyanins producing in H-1. This discovery expanded our knowledges of the biosynthesis of anthocyanins and provided a direction for the production of anthocyanin.

Integrated Analysis of Long Non-coding RNAs (lncRNAs) and mRNAs Reveals the Regulatory Role of lncRNAs Associated With Salt Resistance in Camellia sinensis

Tea plant (Camellia sinensis), an important economic crop, is seriously affected by various abiotic stresses, including salt stress, which severely diminishes its widespread planting. However, little is known about the roles of long non-coding RNAs (lncRNAs) in transcriptional regulation under salt stress. In this study, high-throughput sequencing of tea shoots under salt-stress and control conditions was performed. Through sequencing analysis, 16,452 unique lncRNAs were identified, including 172 differentially expressed lncRNAs (DE-lncRNAs). The results of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of their cis- and trans-target genes showed that these DE-lncRNAs play important roles in many pathways such as the galactinol synthase (GOLS), calcium signaling pathway, and interact with transcription factors (TFs) under salt stress. The data from the gene-specific antisense oligodeoxynucleotide-mediated reduction in the lncRNA MSTRG.139242.1 and its predicted interacting gene, TEA027212.1 (Ca²⁺-ATPase 13), in tea leaves revealed that MSTRG.139242.1 may function in the response of tea plants to high salinity. In addition, 12 lncRNAs were predicted to be target mimics of 17 known mature miRNAs, such as miR156, that are related to the salt-stress response in C. sinensis. Our results provide new insights into lncRNAs as ubiquitous regulators in response to salt stress in tea plants.

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.

Genome-Wide Identification of Long Non-coding RNA in Trifoliate Orange (Poncirus trifoliata (L.) Raf) Leaves in Response to Boron Deficiency

Long non-coding RNAs (lncRNAs) play important roles in plant growth and stress responses. As a dominant abiotic stress factor in soil, boron (B) deficiency stress has impacted the growth and development of citrus in the red soil region of southern China. In the present work, we performed a genome-wide identification and characterization of lncRNAs in response to B deficiency stress in the leaves of trifoliate orange (Poncirus trifoliata), an important rootstock of citrus. A total of 2101 unique lncRNAs and 24,534 mRNAs were predicted. Quantitative real-time polymerase chain reaction (qRT-PCR) experiments were performed for a total of 16 random mRNAs and lncRNAs to validate their existence and expression patterns. Expression profiling of the leaves of trifoliate orange under B deficiency stress identified 729 up-regulated and 721 down-regulated lncRNAs, and 8419 up-regulated and 8395 down-regulated mRNAs. Further analysis showed that a total of 84 differentially expressed lncRNAs (DELs) were up-regulated and 31 were down-regulated, where the number of up-regulated DELs was 2.71-fold that of down-regulated. A similar trend was also observed in differentially expressed mRNAs (DEMs, 4.21-fold). Functional annotation of these DEMs was performed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, and the results demonstrated an enrichment of the categories of the biosynthesis of secondary metabolites (including phenylpropanoid biosynthesis/lignin biosynthesis), plant hormone signal transduction and the calcium signaling pathway. LncRNA target gene enrichment identified several target genes that were involved in plant hormones, and the expression of lncRNAs and their target genes was significantly influenced. Therefore, our results suggest that lncRNAs can regulate the metabolism and signal transduction of plant hormones, which play an important role in the responses of citrus plants to B deficiency stress. Co-expression network analysis indicated that 468 significantly differentially expressed genes may be potential targets of 90 lncRNAs, and a total of 838 matched lncRNA-mRNA pairs were identified. In summary, our data provides a rich resource of candidate lncRNAs and mRNAs, as well as their related pathways, thereby improving our understanding of the role of lncRNAs in response to B deficiency stress, and in symptom formation caused by B deficiency in the leaves of trifoliate orange.

Identification of conserved and novel miRNAs responsive to heat stress in flowering Chinese cabbage using high-throughput sequencing

Plant microRNAs (miRNAs) are noncoding and endogenous key regulators that play significant functions in regulating plant responses to stress, and plant growth and development. Heat stress is a critical abiotic stress that reduces the yield and quality of flowering Chinese cabbage (Brassica campestris L. ssp. chinensis var. utilis Tsen et Lee). However, limited information is available on whether miRNAs are involved in the regulation of heat stress in B. campestris. A high-throughput sequencing approach was used to identify novel and conserved heat-responsive miRNAs in four small RNA libraries of flowering Chinese cabbage using leaves collected at 0 h, 1 h, 6 h and 12 h after a 38 °C heat-stress treatment. The analysis identified 41 conserved miRNAs (belonging to 19 MIR families), of which MIR156, MIR159, MIR168, MIR171 and MIR1885 had the most abundant molecules. Prediction and evaluation of novel miRNAs using the unannotated reads resulted in 18 candidate miRNAs. Differential expression analysis showed that most of the identified miRNAs were downregulated in heat-treated groups. To better understand functional importance, bioinformatic analysis predicted 432 unique putative target miRNAs involved in cells, cell parts, catalytic activity, cellular processes and abiotic stress responses. Furthermore, the Kyoto Encyclopedia of Genes and Genomes maps of flowering Chinese cabbage identified the significant role of miRNAs in stress adaptation and stress tolerance, and in several mitogen-activated protein kinases signaling pathways including cell death. This work presents a comprehensive study of the miRNAs for understanding the regulatory mechanisms and their participation in the heat stress of flowering Chinese cabbage.

Transcription Factors Associated with Abiotic and Biotic Stress Tolerance and Their Potential for Crops Improvement

In field conditions, crops are adversely affected by a wide range of abiotic stresses including drought, cold, salt, and heat, as well as biotic stresses including pests and pathogens. These stresses can have a marked effect on crop yield. The present and future effects of climate change necessitate the improvement of crop stress tolerance. Plants have evolved sophisticated stress response strategies, and genes that encode transcription factors (TFs) that are master regulators of stress-responsive genes are excellent candidates for crop improvement. Related examples in recent studies include TF gene modulation and overexpression approaches in crop species to enhance stress tolerance. However, much remains to be discovered about the diverse plant TFs. Of the >80 TF families, only a few, such as NAC, MYB, WRKY, bZIP, and ERF/DREB, with vital roles in abiotic and biotic stress responses have been intensively studied. Moreover, although significant progress has been made in deciphering the roles of TFs in important cereal crops, fewer TF genes have been elucidated in sorghum. As a model drought-tolerant crop, sorghum research warrants further focus. This review summarizes recent progress on major TF families associated with abiotic and biotic stress tolerance and their potential for crop improvement, particularly in sorghum. Other TF families and non-coding RNAs that regulate gene expression are discussed briefly. Despite the emphasis on sorghum, numerous examples from wheat, rice, maize, and barley are included. Collectively, the aim of this review is to illustrate the potential application of TF genes for stress tolerance improvement and the engineering of resistant crops, with an emphasis on sorghum.

Small RNA Mobility: Spread of RNA Silencing Effectors and its Effect on Developmental Processes and Stress Adaptation in Plants

Plants are exposed every day to multiple environmental cues, and tight transcriptome reprogramming is necessary to control the balance between responses to stress and processes of plant growth. In this context, the silencing phenomena mediated by small RNAs can drive transcriptional and epigenetic regulatory modifications, in turn shaping plant development and adaptation to the surrounding environment. Mounting experimental evidence has recently pointed to small noncoding RNAs as fundamental players in molecular signalling cascades activated upon exposure to abiotic and biotic stresses. Although, in the last decade, studies on stress responsive small RNAs increased significantly in many plant species, the physiological responses triggered by these molecules in the presence of environmental stresses need to be further explored. It is noteworthy that small RNAs can move either cell-to-cell or systemically, thus acting as mobile silencing effectors within the plant. This aspect has great importance when physiological changes, as well as epigenetic regulatory marks, are inspected in light of plant environmental adaptation. In this review, we provide an overview of the categories of mobile small RNAs in plants, particularly focusing on the biological implications of non-cell autonomous RNA silencing in the stress adaptive response and epigenetic modifications.

Non-Coding RNAs and their Integrated Networks

Eukaryotic genomes are pervasively transcribed. Besides protein-coding RNAs, there are different types of non-coding RNAs that modulate complex molecular and cellular processes. RNA sequencing technologies and bioinformatics methods greatly promoted the study of ncRNAs, which revealed ncRNAs' essential roles in diverse aspects of biological functions. As important key players in gene regulatory networks, ncRNAs work with other biomolecules, including coding and non-coding RNAs, DNAs and proteins. In this review, we discuss the distinct types of ncRNAs, including housekeeping ncRNAs and regulatory ncRNAs, their versatile functions and interactions, transcription, translation, and modification. Moreover, we summarize the integrated networks of ncRNA interactions, providing a comprehensive landscape of ncRNAs regulatory roles.

Non-coding RNAs and transposable elements in plant genomes: emergence, regulatory mechanisms and roles in plant development and stress responses

This review will provide evidence for the indispensable function of these elements in regulating plant development and resistance to biotic and abiotic stresses, as well as their evolutionary role in facilitating plant adaptation. Over millions of years of evolution, plant genomes have acquired a complex constitution. Plant genomes consist not only of protein coding sequences, but also contain large proportions of non-coding sequences. These include introns of protein-coding genes, and intergenic sequences such as non-coding RNA, repeat sequences and transposable elements. These non-coding sequences help to regulate gene expression, and are increasingly being recognized as playing an important role in genome organization and function. In this review, we summarize the known molecular mechanisms by which gene expression is regulated by several species of non-coding RNAs (microRNAs, long non-coding RNAs, and circular RNAs) and by transposable elements. We further discuss how these non-coding RNAs and transposable elements evolve and emerge in the genome, and the potential influence and importance of these non-coding RNAs and transposable elements in plant development and in stress responses.

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

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