Genome-wide identification and functional prediction of nitrogen-responsive intergenic and intronic long non-coding RNAs in maize (Zea mays L.)

Exploring the emerging role of long non-coding RNAs (lncRNAs) in plant biology: Functions, mechanisms of action, and future directions

Long non-coding RNAs (lncRNAs) are a class of RNA transcripts that surpass 200 nucleotides in length and lack discernible coding potential. LncRNAs that have been functionally characterized have pivotal functions in several plant processes, including the regulation of flowering, and development of lateral roots. It also plays a crucial role in the plant's response to abiotic stressors and exhibits vital activities in environmental adaptation. The progress in NGS (next-generation sequencing) and functional genomics technology has facilitated the discovery of lncRNA in plant species. This review is a brief explanation of lncRNA genomics, its molecular role, and the mechanism of action in plants. The review also addresses the challenges encountered in this field and highlights promising molecular and computational methodologies that can aid in the comparative and functional analysis of lncRNAs.

Transcriptome Analysis of Sesame (Sesamum indicum L.) Reveals the LncRNA and mRNA Regulatory Network Responding to Low Nitrogen Stress

Nitrogen is one of the important factors restricting the development of sesame planting and industry in China. Cultivating sesame varieties tolerant to low nitrogen is an effective way to solve the problem of crop nitrogen deficiency. To date, the mechanism of low nitrogen tolerance in sesame has not been elucidated at the transcriptional level. In this study, two sesame varieties Zhengzhi HL05 (ZZ, nitrogen efficient) and Burmese prolific (MD, nitrogen inefficient) in low nitrogen were used for RNA-sequencing. A total of 3964 DEGs (differentially expressed genes) and 221 DELs (differentially expressed lncRNAs) were identified in two sesame varieties at 3d and 9d after low nitrogen stress. Among them, 1227 genes related to low nitrogen tolerance are mainly located in amino acid metabolism, starch and sucrose metabolism and secondary metabolism, and participate in the process of transporter activity and antioxidant activity. In addition, a total of 209 pairs of lncRNA-mRNA were detected, including 21 pairs of trans and 188 cis. WGCNA (weighted gene co-expression network analysis) analysis divided the obtained genes into 29 modules; phenotypic association analysis identified three low-nitrogen response modules; through lncRNA-mRNA co-expression network, a number of hub genes and cis/trans-regulatory factors were identified in response to low-nitrogen stress including GS1-2 (glutamine synthetase 1-2), PAL (phenylalanine ammonia-lyase), CHS (chalcone synthase, CHS), CAB21 (chlorophyll a-b binding protein 21) and transcription factors MYB54, MYB88 and NAC75 and so on. As a trans regulator, lncRNA MSTRG.13854.1 affects the expression of some genes related to low nitrogen response by regulating the expression of MYB54, thus responding to low nitrogen stress. Our research is the first to provide a more comprehensive understanding of DEGs involved in the low nitrogen stress of sesame at the transcriptome level. These results may reveal insights into the molecular mechanisms of low nitrogen tolerance in sesame and provide diverse genetic resources involved in low nitrogen tolerance research.

Characterization of arbuscular mycorrhizal fungal species associating with Zea mays

Taxonomic identification of arbuscular mycorrhizal (AM) fungal spores extracted directly from the field is sometimes difficult because spores are often degraded or parasitized by other organisms. Single-spore inoculation of a suitable host plant allows for establishing monosporic cultures of AM fungi. This study aimed to propagate AM fungal spores isolated from maize soil using single spores for morphological characterization. First, trap cultures were established to trigger the sporulation of AM fungal species. Second, trap cultures were established with individual morphotypes by picking up only one spore under a dissecting microscope and transferring it to a small triangle of sterilized filter paper, which was then carefully inoculated below a root from germinated sorghum seeds in each pot and covered with a sterile substrate. All pots were placed in sunbags and maintained in a plant growth room for 120 days. Spores obtained from single spore trap cultures from each treatment, maize after oats (MO), maize after maize (MM), maize after peas (MP), and maize after soybean (MS), were extracted using the sieving method. Healthy spores were selected for morphological analysis. Direct PCR was conducted by crushing spores in RNAlater and applying three sets of primer pairs: ITS1 × ITS4, NS31 × AML2, and SSUmcf and LSUmBr. Nucleotide sequences obtained from Sanger sequencing were aligned on MEGA X. The phylogenetic tree showed that the closest neighbors of the propagated AM fungal species belonged to the genera Claroideoglomus, Funneliformis, Gigaspora, Paraglomus, and Rhizophagus. The morphological characteristics were compared to the descriptive features of described species posted on the INVAM website, and they included Acaulospora cavernata, Diversispora spurca, Funneliformis geosporus, Funneliformis mosseae, Gigaspora clarus, Gigaspora margarita, Glomus macrosporum, Paraglomus occultum, and Rhizophagus intraradices. These findings can provide a great contribution to crop productivity and sustainable management of the agricultural ecosystem. Also, the isolate analyzed could be grouped into efficient promoters of growth and mycorrhization of maize independent of their geographical location.

Long noncoding RNA TRABA suppresses β-glucosidase-encoding BGLU24 to promote salt tolerance in cotton

Salt stress severely damages the growth and yield of crops. Recently, long noncoding RNAs (lncRNAs) were demonstrated to regulate various biological processes and responses to environmental stresses. However, the regulatory mechanisms of lncRNAs in cotton (Gossypium hirsutum) response to salt stress are still poorly understood. Here, we observed that a lncRNA, trans acting of BGLU24 by lncRNA (TRABA), was highly expressed while GhBGLU24-A was weakly expressed in a salt-tolerant cotton accession (DM37) compared to a salt-sensitive accession (TM-1). Using TRABA as an effector and proGhBGLU24-A-driven GUS as a reporter, we showed that TRABA suppressed GhBGLU24-A promoter activity in double transgenic Arabidopsis (Arabidopsis thaliana), which explained why GhBGLU24-A was weakly expressed in the salt-tolerant accession compared to the salt-sensitive accession. GhBGLU24-A encodes an endoplasmic reticulum (ER)-localized β-glucosidase that responds to salt stress. Further investigation revealed that GhBGLU24-A interacted with RING-type E3 ubiquitin ligase (GhRUBL). Virus-induced gene silencing (VIGS) and transgenic Arabidopsis studies revealed that both GhBGLU24-A and GhRUBL diminish plant tolerance to salt stress and ER stress. Based on its substantial effect on ER-related degradation (ERAD)-associated gene expression, GhBGLU24-A mediates ER stress likely through the ERAD pathway. These findings provide insights into the regulatory role of the lncRNA TRABA in modulating salt and ER stresses in cotton and have potential implications for developing more resilient crops.

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.

Long noncoding RNA from Betula platyphylla, BplncSIR1, confers salt tolerance by regulating BpNAC2 to mediate reactive oxygen species scavenging and stomatal movement

Long noncoding RNAs (lncRNAs) play an important role in abiotic stress tolerance. However, their function in conferring abiotic stress tolerance is still unclear. Herein, we characterized the function of a salt-responsive nuclear lncRNA (BplncSIR1) from Betula platyphylla (birch). Birch plants overexpressing and knocking out for BplncSIR1 were generated. BplncSIR1 was found to improve salt tolerance by inducing antioxidant activity and stomatal closure, and also accelerate plant growth. Chromatin isolation by RNA purification (ChIRP) combined with RNA sequencing indicated that BplncSIR1 binds to the promoter of BpNAC2 (encoding NAC domain-containing protein 2) to activate its expression. Plants overexpressing and knocking out for BpNAC2 were generated. Consistent with that of BplncSIR1, overexpression of BpNAC2 also accelerated plant growth and conferred salt tolerance. In addition, BpNAC2 binds to different cis-acting elements, such as G-box and 'CCAAT' sequences, to regulate the genes involved in salt tolerance, resulting in reduced ROS accumulation and decreased water loss rate by stomatal closure. Taken together, BplncSIR1 serves as the regulator of BpNAC2 to induce its expression in response to salt stress, and activated BpNAC2 accelerates plant growth and improves salt tolerance. Therefore, BplncSIR1 might be a candidate gene for molecular breeding to cultivate plants with both a high growth rate and improved salt tolerance.

Integrated analysis of genes and long non-coding RNAs in trachea transcriptome to decipher the host response during Newcastle disease challenge in different breeds of chicken

Newcastle disease is a highly infectious economically devastating disease caused by Newcastle disease Virus in Chicken (Gallus gallus). Leghorn and Fayoumi are two breeds which show differential resistance patterns towards NDV. This study aims to identify the differentially expressed genes and lncRNAs during NDV challenge which could play a potential role in this differential resistance pattern. A total of 552 genes and 1580 lncRNAs were found to be differentially expressing. Of them, 52 genes were annotated with both Immune related pathways and Gene ontologies. We found that most of these genes were upregulated in Leghorn between normal and challenged chicken but several were down regulated between different timepoints after NDV challenge, while Fayoumi showed no such downregulation. We also observed that higher number of positively correlating lncRNAs was found to be downregulated along with these genes. This shows that although Leghorn is showing higher number of differentially expressed genes in challenged than in non-challenged, most of them were downregulated during the disease between different timepoints. With this we hypothesize that the downregulation of immune related genes and co-expressing lncRNAs could play a significant role behind the Leghorn being comparatively susceptible breed than Fayoumi. The computational pipeline is available at https://github.com/Venky2804/FHSpipe.

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

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

The expression landscape and pangenome of long non-coding RNA in the fungal wheat pathogen Zymoseptoria tritici

Long non-coding RNAs (lncRNAs) are regulatory molecules interacting in a wide array of biological processes. lncRNAs in fungal pathogens can be responsive to stress and play roles in regulating growth and nutrient acquisition. Recent evidence suggests that lncRNAs may also play roles in virulence, such as regulating pathogenicity-associated enzymes and on-host reproductive cycles. Despite the importance of lncRNAs, only a few model fungi have well-documented inventories of lncRNA. In this study, we apply a recent computational pipeline to predict high-confidence lncRNA candidates in Zymoseptoria tritici, an important global pathogen of wheat impacting global food production. We analyse genomic features of lncRNAs and the most likely associated processes through analyses of expression over a host infection cycle. We find that lncRNAs are frequently expressed during early infection, before the switch to necrotrophic growth. They are mostly located in facultative heterochromatic regions, which are known to contain many genes associated with pathogenicity. Furthermore, we find that lncRNAs are frequently co-expressed with genes that may be involved in responding to host defence signals, such as oxidative stress. Finally, we assess pangenome features of lncRNAs using four additional reference-quality genomes. We find evidence that the repertoire of expressed lncRNAs varies substantially between individuals, even though lncRNA loci tend to be shared at the genomic level. Overall, this study provides a repertoire and putative functions of lncRNAs in Z. tritici enabling future molecular genetics and functional analyses in an important pathogen.

Progress of Research on the Physiology and Molecular Regulation of Sorghum Growth under Salt Stress by Gibberellin

Plant growth often encounters diverse abiotic stresses. As a global resource-based ecological problem, salinity is widely distributed and one of the major abiotic stresses affecting crop yields worldwide. Sorghum, a cereal crop with medium salt tolerance and great value for the development and utilization of salted soils, is an important source of food, brewing, energy, and forage production. However, in soils with high salt concentrations, sorghum experiences low emergence and suppressed metabolism. It has been demonstrated that the effects of salt stress on germination and seedling growth can be effectively mitigated to a certain extent by the exogenous amendment of hormonal gibberellin (GA). At present, most of the studies on sorghum salt tolerance at home and abroad focus on morphological and physiological levels, including the transcriptome analysis of the exogenous hormone on sorghum salt stress tolerance, the salt tolerance metabolism pathway, and the mining of key salt tolerance regulation genes. The high-throughput sequencing technology is increasingly widely used in the study of crop resistance, which is of great significance to the study of plant resistance gene excavation and mechanism. In this study, we aimed to review the effects of the exogenous hormone GA on leaf morphological traits of sorghum seedlings and further analyze the physiological response of sorghum seedling leaves and the regulation of sorghum growth and development. This review not only focuses on the role of GA but also explores the signal transduction pathways of GA and the performance of their responsive genes under salt stress, thus helping to further clarify the mechanism of regulating growth and production under salt stress. This will serve as a reference for the molecular discovery of key genes related to salt stress and the development of new sorghum varieties.

Uncovering the Gene Regulatory Network of Endothelial Cells in Mouse Duchenne Muscular Dystrophy: Insights from Single-Nuclei RNA Sequencing Analysis

Introduction: Duchenne muscular dystrophy (DMD) is a severe X-linked recessive disorder caused by mutations in the dystrophin gene, which leads to heart and respiratory failure. Despite the critical impact of DMD on endothelial cells (ECs), there is limited understanding of its effect on the endothelial gene network. The aim of this study was to investigate the impact of DMD on the gene regulatory network of ECs. Methods and Results: To gain insights into the role of the dystrophin muscular dystrophy gene (DMD) in ECs from Duchenne muscular dystrophy; the study utilized single-nuclei RNA sequencing (snRNA-seq) to evaluate the transcriptomic profile of ECs from skeletal muscles in DMD mutant mice (DMDmut) and wild-type control mice. The analysis showed that the DMD mutation resulted in the suppression of several genes, including SPTBN1 and the upregulation of multiple long noncoding RNAs (lncRNAs). GM48099, GM19951, and GM15564 were consistently upregulated in ECs and skeletal muscle cells from DMDmut, indicating that these dysregulated lncRNAs are conserved across different cell types. Gene ontology (GO) enrichment analysis revealed that the DMD mutation activated the following four pathways in ECs: fibrillary collagen trimer, banded collagen fibril, complex of collagen trimers, and purine nucleotide metabolism. The study also found that the metabolic pathway activity of ECs was altered. Oxidative phosphorylation (OXPHOS), fatty acid degradation, glycolysis, and pyruvate metabolism were decreased while purine metabolism, pyrimidine metabolism, and one carbon pool by folate were increased. Moreover, the study investigated the impact of the DMD mutation on ECs from skeletal muscles and found a significant decrease in their overall number, but no change in their proliferation. Conclusions: Overall, this study provides new insights into the gene regulatory program in ECs in DMD and highlights the importance of further research in this area.

Long Non-Coding RNAs of Plants in Response to Abiotic Stresses and Their Regulating Roles in Promoting Environmental Adaption

Abiotic stresses triggered by climate change and human activity cause substantial agricultural and environmental problems which hamper plant growth. Plants have evolved sophisticated mechanisms in response to abiotic stresses, such as stress perception, epigenetic modification, and regulation of transcription and translation. Over the past decade, a large body of literature has revealed the various regulatory roles of long non-coding RNAs (lncRNAs) in the plant response to abiotic stresses and their irreplaceable functions in environmental adaptation. LncRNAs are recognized as a class of ncRNAs that are longer than 200 nucleotides, influencing a variety of biological processes. In this review, we mainly focused on the recent progress of plant lncRNAs, outlining their features, evolution, and functions of plant lncRNAs in response to drought, low or high temperature, salt, and heavy metal stress. The approaches to characterize the function of lncRNAs and the mechanisms of how they regulate plant responses to abiotic stresses were further reviewed. Moreover, we discuss the accumulating discoveries regarding the biological functions of lncRNAs on plant stress memory as well. The present review provides updated information and directions for us to characterize the potential functions of lncRNAs in abiotic stresses in the future.

Genome-wide analysis of long non-coding RNAs (lncRNAs) in tea plants (Camellia sinensis) lateral roots in response to nitrogen application

Tea (Camellia sinensis) is one of the significant cash crops in China. As a leaf crop, nitrogen supply can not only increase the number of new shoots and leaves but also improve the tenderness of the former. However, a conundrum remains in science, which is the molecular mechanism of nitrogen use efficiency, especially long non-coding RNA (lncRNA). In this study, a total of 16,452 lncRNAs were identified through high-throughput sequencing analysis of lateral roots under nitrogen stress and control conditions, of which 9,451 were differentially expressed lncRNAs (DE-lncRNAs). To figure out the potential function of nitrogen-responsive lncRNAs, co-expression clustering was employed between lncRNAs and coding genes. KEGG enrichment analysis revealed nitrogen-responsive lncRNAs may involve in many biological processes such as plant hormone signal transduction, nitrogen metabolism and protein processing in endoplasmic reticulum. The expression abundance of 12 DE-lncRNAs were further verified by RT-PCR, and their expression trends were consistent with the results of RNA-seq. This study expands the research on lncRNAs in tea plants, provides a novel perspective for the potential regulation of lncRNAs on nitrogen stress, and valuable resources for further improving the nitrogen use efficiency of tea plants.

Regulation of seed oil accumulation by lncRNAs in Brassica napus

BACKGROUND

Studies have indicated that long non-coding RNAs (lncRNAs) play important regulatory roles in many biological processes. However, the regulation of seed oil biosynthesis by lncRNAs remains largely unknown.

RESULTS

We comprehensively identified and characterized the lncRNAs from seeds in three developing stages in two accessions of Brassica napus (B. napus), ZS11 (high oil content) and WH5557 (low oil content). Finally, 8094 expressed lncRNAs were identified. LncRNAs MSTRG.22563 and MSTRG.86004 were predicted to be related to seed oil accumulation. Experimental results show that the seed oil content is decreased by 3.1-3.9% in MSTRG.22563 overexpression plants, while increased about 2% in MSTRG.86004, compared to WT. Further study showed that most genes related to lipid metabolism had much lower expression, and the content of some metabolites in the processes of respiration and TCA (tricarboxylic acid) cycle was reduced in MSTRG.22563 transgenic seeds. The expression of genes involved in fatty acid synthesis and seed embryonic development (e.g., LEC1) was increased, but genes related to TAG assembly was decreased in MSTRG.86004 transgenic seeds.

CONCLUSION

Our results suggest that MSTRG.22563 might impact seed oil content by affecting the respiration and TCA cycle, while MSTRG.86004 plays a role in prolonging the seed developmental time to increase seed oil accumulation.

Role of long non coding RNA in plants under abiotic and biotic stresses

Evolutionary processes have evolved plants to cope with several different natural stresses. Basic physiological activities of crop plants are significantly harmed by these stresses, reducing productivity and eventually leading to death. The recent advancements in high-throughput sequencing of transcriptome and expression profiling with NGS techniques lead to the innovation of various RNAs which do not code for proteins, more specifically long non-coding RNAs (lncRNAs), undergirding regulate growth, development, and the plant defence mechanism transcriptionally under stress situations. LncRNAs are a diverse set of RNAs that play key roles in various biological processes at the level of transcription, post-transcription, and epigenetics. These are thought to serve crucial functions in plant immunity and response to changes in the environment. In plants, however, just a few lncRNAs have been functionally identified. In this review, we will address recent advancements in comprehending lncRNA regulatory functions, focusing on the expanding involvement of lncRNAs in modulating environmental stress responsiveness in plants.

Transcriptome dynamics uncovers long non-coding RNAs response to salinity stress in Chenopodium quinoa

Chenopodium quinoa is a crop with outstanding tolerance to saline soil, but long non-coding RNAs (LncRNAs) expression profile driven by salt stress in quinoa has rarely been observed yet. Based on the high-quality quinoa reference genome and high-throughput RNA sequencing (RNA-seq), genome-wide identification of LncRNAs was performed, and their dynamic response under salt stress was then investigated. In total, 153,751 high-confidence LncRNAs were discovered and dispersed intensively in chromosomes. Expression profile analysis demonstrated significant differences between LncRNAs and coding RNAs. Under salt stress conditions, 4,460 differentially expressed LncRNAs were discovered, of which only 54 were differentially expressed at all the stress time points. Besides, strongly significantly correlation was observed between salt-responsive LncRNAs and their closest neighboring genes (r = 0.346, p-value < 2.2e-16). Furthermore, a weighted co-expression network was then constructed to infer the potential biological functions of LncRNAs. Seven modules were significantly correlated with salt treatments, resulting in 210 hub genes, including 22 transcription factors and 70 LncRNAs. These results indicated that LncRNAs might interact with transcription factors to respond to salinity stress. Gene ontology enrichment of the coding genes of these modules showed that they were highly related to regulating metabolic processes, biological regulation and response to stress. This study is the genome-wide analysis of the LncRNAs responding to salt stress in quinoa. The findings will provide a solid framework for further functional research of salt responsive LncRNAs, contributing to quinoa genetic improvement.

Global Identification of White Lupin lncRNAs Reveals Their Role in Cluster Roots under Phosphorus Deficiency

Phosphorus (P) deficiency heterogeneously affected plant nutritional status and physiological performance, ultimately leading to a severe yield reduction. A few putative long non-coding RNAs (lncRNAs) responding to P-starvation in the model crops Arabidopsis thaliana and Oryza sativa have been characterized. White lupin (Lupinus albus) is of prime importance, and is a legume with increasing agronomic value as a protein crop as it exhibits extreme tolerance to nutrient deficiency, particularly P deficiency. Despite its adapted nature to P deficiency, nothing is known about low P-induced lncRNAs in white lupin roots. To address this issue, we identified 39,840 mRNA and 2028 lncRNAs in the eight developmental stages of white lupin root (S0–S7 and lateral root, LR) grown under P deficiency. From these 2028 lncRNAs, 1564 were intergenic and 464 natural antisense intergenic transcript (NAT) lncRNAs. We further predicted six potential targets of miRNAs with twelve lncRNAs, which may regulate P-deficiency-related processes. Moreover, the weighted gene co-expression network analysis (WGCNA) revealed seven modules that were correlated with the expression pattern of lncRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed 606 GO terms and 27 different pathways including signal transduction, energy synthesis, detoxification, and Pi transport. In addition, we screened 13 putative lncRNAs that showed a distinct expression pattern in each root, indicating their role in the P deficiency regulatory network. Therefore, white lupin may be a reference legume to characterize P-deficiency-responsive novel lncRNAs, which would highlight the role of lncRNAs in the regulation of plant responses to P deficiency.

Identification and characterization of heat-responsive lncRNAs in maize inbred line CM1

BACKGROUND

Frequent occurrence of extreme high temperature is a major threat to crop production. Increasing evidence demonstrates that long non-coding RNAs (lncRNAs) have important biological functions in the regulation of the response to heat stress. However, the regulatory mechanism of lncRNAs involved in heat response requires further exploration and the regulatory network remains poorly understood in maize.

RESULTS

In this research, high-throughput sequencing was adopted to systematically identify lncRNAs in maize inbred line CM1. In total, 53,249 lncRNAs (259 known lncRNAs and 52,990 novel lncRNAs) were detected, of which 993 lncRNAs showed significantly differential expression (DElncRNAs) under heat stress. By predicting the target genes, 953 common targets shared by cis- and trans-regulation of the DElncRNAs were identified, which exhibited differential expression between the control and the heat stress treatments. Functional annotation indicated that a number of important biological processes and pathways, including photosynthesis, metabolism, translation, stress response, hormone signal transduction, and spliceosome, were enriched for the common targets, suggesting that they play important roles in heat response. A lncRNA-mediated regulatory network was constructed to visualize the molecular response mechanism in response to heat stress, which represented the direct regulatory relationships of DElncRNAs, differentially expressed miRNAs, target genes, and functional annotations.

CONCLUSIONS

This study lays a foundation for further elucidation of the regulatory mechanism for the response to heat stress in the maize inbred line CM1. The findings provide important information for identification of heat-responsive genes, which will be beneficial for the molecular breeding in the cultivation of heat-tolerant maize germplasm.

Whole-transcriptome analysis and construction of an anther development-related ceRNA network in Chinese cabbage (Brassica campestris L. ssp. pekinensis)

Anther development is precisely regulated by a complex gene network, which is of great significance to plant breeding. However, the molecular mechanism of anther development in Chinese cabbage is unclear. Here, we identified microRNAs (miRNAs), mRNAs, long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs) related to anther development in Chinese cabbage (Brassica campestris L. ssp. pekinensis) to construct competitive endogenous RNA (ceRNA) regulatory networks and provide valuable knowledge on anther development. Using whole-transcriptome sequencing, 9055, 585, 1344, and 165 differentially expressed mRNAs (DEmRNAs), miRNAs (DEmiRNAs), lncRNAs (DElncRNAs), and circRNAs (DEcircRNAs) were identified, respectively, in the anthers of Chinese cabbage compared with those in samples of the vegetative mass of four true leaves. An anther-related ceRNA regulatory network was constructed using miRNA targeting relationships, and 450 pairs of ceRNA relationships, including 97 DEmiRNA-DEmRNA, 281 DEmiRNA-DElncRNA, and 23 DEmiRNA-DEcircRNA interactions, were obtained. We identified important genes and their interactions with lncRNAs, circRNAs, and miRNAs involved in microsporogenesis, tapetum and callose layer development, pollen wall formation, and anther dehiscence. We analyzed the promoter activity of six predominant anther expression genes, which were expressed specifically in the anthers of Arabidopsis thaliana, indicating that they may play an important role in anther development of Chinese cabbage. This study lays the foundation for further research on the molecular mechanisms of anther growth and development in Chinese cabbage.

Integrated analysis of long non-coding RNAs and mRNAs reveals the regulatory network of maize seedling root responding to salt stress

Background

Long non-coding RNAs (lncRNAs) play important roles in response to abiotic stresses in plants, by acting as cis- or trans-acting regulators of protein-coding genes. As a widely cultivated crop worldwide, maize is sensitive to salt stress particularly at the seedling stage. However, it is unclear how the expressions of protein-coding genes are affected by non-coding RNAs in maize responding to salt tolerance.

Results

The whole transcriptome sequencing was employed to investigate the differential lncRNAs and target transcripts responding to salt stress between two maize inbred lines with contrasting salt tolerance. We developed a flexible, user-friendly, and modular RNA analysis workflow, which facilitated the identification of lncRNAs and novel mRNAs from whole transcriptome data. Using the workflow, 12,817 lncRNAs and 8,320 novel mRNAs in maize seedling roots were identified and characterized. A total of 742 lncRNAs and 7,835 mRNAs were identified as salt stress-responsive transcripts. Moreover, we obtained 41 cis- and 81 trans-target mRNA for 88 of the lncRNAs. Among these target transcripts, 11 belonged to 7 transcription factor (TF) families including bHLH, C2H2, Hap3/NF-YB, HAS, MYB, WD40, and WRKY. The above 8,577 salt stress-responsive transcripts were further classified into 28 modules by weighted gene co-expression network analysis. In the salt-tolerant module, we constructed an interaction network containing 79 nodes and 3081 edges, which included 5 lncRNAs, 18 TFs and 56 functional transcripts (FTs). As a trans-acting regulator, the lncRNA MSTRG.8888.1 affected the expressions of some salt tolerance-relative FTs, including protein-serine/threonine phosphatase 2C and galactinol synthase 1, by regulating the expression of the bHLH TF.

Conclusions

The contrasting genetic backgrounds of the two inbred lines generated considerable variations in the expression abundance of lncRNAs and protein-coding transcripts. In the co-expression networks responding to salt stress, some TFs were targeted by the lncRNAs, which further regulated the salt tolerance-related functional transcripts. We constructed a regulatory pathway of maize seedlings to salt stress, which was mediated by the hub lncRNA MSTRG.8888.1 and participated by the bHLH TF and its downstream target transcripts. Future work will be focused on the functional revelation of the regulatory pathway.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08286-7.

Regulatory long non-coding RNAs in root growth and development

As sessile organisms, plants have evolved sophisticated mechanisms of gene regulation to cope with changing environments. Among them, long non-coding RNAs (lncRNAs) are a class of RNAs regulating gene expression at both transcriptional and post-transcriptional levels. They are highly responsive to environmental cues or developmental processes and are generally involved in fine-tuning plant responses to these signals. Roots, in addition to anchoring the plant to the soil, allow it to absorb the major part of its mineral nutrients and water. Furthermore, roots directly sense environmental constraints such as mineral nutrient availability and abiotic or biotic stresses and dynamically adapt their growth and architecture. Here, we review the role of lncRNAs in the control of root growth and development. In particular, we highlight their action in fine-tuning primary root growth and the development of root lateral organs, such as lateral roots and symbiotic nodules. Lastly, we report their involvement in plant response to stresses and the regulation of nutrient assimilation and homeostasis, two processes leading to the modification of root architecture. LncRNAs could become interesting targets in plant breeding programs to subtly acclimate crops to coming environmental changes.

Nitrogen use efficiency (NUE): elucidated mechanisms, mapped genes and gene networks in maize (Zea mays L.)

Nitrogen, the vital primary plant growth nutrient at deficit soil conditions, drastically affects the growth and yield of a crop. Over the years, excess use of inorganic nitrogenous fertilizers resulted in pollution, eutrophication and thereby demanding the reduction in the use of chemical fertilizers. Being a C₄ plant with fibrous root system and high NUE, maize can be deployed to be the best candidate for better N uptake and utilization in nitrogen deficient soils. The maize germplasm sources has enormous genetic variation for better nitrogen uptake contributing traits. Adoption of single cross maize hybrids as well as inherent property of high NUE has helped maize cultivars to achieve the highest growth rate among the cereals during last decade. Further, considering the high cost of nitrogenous fertilizers, adverse effects on soil health and environmental impact, maize improvement demands better utilization of existing genetic variation for NUE via introgression of novel allelic combinations in existing cultivars. Marker assisted breeding efforts need to be supplemented with introgression of genes/QTLs related to NUE in ruling varieties and thereby enhancing the overall productivity of maize in a sustainable manner. To achieve this, we need mapped genes and network of interacting genes and proteins to be elucidated. Identified genes may be used in screening ideal maize genotypes in terms of better physiological functionality exhibiting high NUE. Future genome editing may help in developing lines with increased productivity under low N conditions in an environment of optimum agronomic practices.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01113-z.

Identification of Long Non-Coding RNAs Associated with Tomato Fruit Expansion and Ripening by Strand-Specific Paired-End RNA Sequencing

As emerging essential regulators in plant development, long non-coding RNAs (lncRNAs) have been extensively investigated in multiple horticultural crops, as well as in different tissues of plants. Tomato fruits are an indispensable part of people’s diet and are consumed as fruits and vegetables. Meanwhile, tomato is widely used as a model to study the ripening mechanism in fleshy fruit. Although increasing evidence shows that lncRNAs are involved in lots of biological processes in tomato plants, the comprehensive identification of lncRNAs in tomato fruit during its expansion and ripening and their functions are partially known. Here, we performed strand-specific paired-end RNA sequencing (ssRNA-seq) of tomato Heinz1706 fruits at five different developmental stages, as well as flowers and leaves. We identified 17,674 putative lncRNAs by referencing the recently released SL4.0 and annotation ITAG4.0 in tomato plants. Many lncRNAs show different expression patterns in fleshy fruit at different developmental stages compared with leaves or flowers. Our results indicate that lncRNAs play an important role in the regulation of tomato fruit expansion and ripening, providing informative lncRNA candidates for further studies in tomato fruits. In addition, we also summarize the recent advanced progress in lncRNAs mediated regulation on horticultural fruits. Hence, our study updates the understanding of lncRNAs in horticultural plants and provides resources for future studies relating to the expansion and ripening of tomato fruits.

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.

Expression Patterns and Regulation of Non-Coding RNAs during Synthesis of Cellulose in Eucalyptus grandis Hill

Cellulose, an essential structural component in the plant cell wall and a renewable biomass resource, plays a significant role in nature. Eucalyptus’s excellent timber tree species (including Eucalyptus grandis Hill) provide many raw materials for the paper and wood industries. The synthesis of cellulose is a very complex process involving multiple genes and regulated by various biological networks. However, research on regulating associated genes and non-coding RNAs during cellulose synthesis in E. grandis remains lacking. In this study, the wood anatomical characteristics and chemical indexes of E. grandis were analyzed by taking three different parts (diameter at breast height (DBH), middle and upper part of the trunk) from the main stem of E. grandis as raw materials. The role of non-coding RNAs (Long non-coding RNA, lncRNA; Micro RNA, miRNA; Circle RNA, circRNA) on regulating candidate genes was presented, and the network map of ceRNA (Competing endogenous RNA) regulation during wood cellulose biosynthesis of E. grandis was constructed. The transcriptome sequencing of nine samples obtained from the trunk of the immature xylem in E. grandis at DBH, middle and upper parts had a 95.81 G clean reading, 57,480 transcripts, 7365 lncRNAs, and 5180 circRNAs. Each sample had 172–306 known miRNAs and 1644–3508 new miRNAs. A total of 190 DE-lncRNAs (Differentially expressed long non-coding RNAs), 174 DE-miRNAs (Differentially expressed micro RNAs), and 270 DE-circRNAs (Differentially expressed circle RNAs) were obtained by comparing transcript expression levels. Four lncRNAs and nine miRNAs were screened out, and the ceRNA regulatory network was constructed. LncRNA1 and lncRNA4 regulated the genes responsible for cellulose synthesis in E. grandis, which were overexpressed in 84K (Populus Alba × Populus glandulosa) poplar. The cellulose and lignin content in lncRNA4-oe were significantly higher than wild type 84K poplar and lncRNA1-oe. The average plant height, middle and basal part of the stem diameter in lncRNA4-oe were significantly higher than the wild type. However, there was no significant difference between the growth of lncRNA1-oe and the wild type. Further studies are warranted to explore the molecular regulatory mechanism of cellulose biosynthesis in Eucalyptus species.

Bacterial and archaeal diversity in oil fields and reservoirs and their potential role in hydrocarbon recovery and bioprospecting

Hydrocarbon is a primary source of energy in the current urbanized society. Considering the increasing demand, worldwide oil productions are declining due to maturity of oil fields and because of difficulty in discovering new oil fields to substitute the exploited ones. To meet current and future energy demands, further exploitation of oil resources is highly required. Microorganisms inhabiting in these areas exhibit highly diverse catabolic activities to degrade, transform, or accumulate various hydrocarbons. Enrichment of hydrocarbon-utilizing bacteria in oil basin is caused by continuous long duration and low molecular weight hydrocarbon microseepage which plays a very important role as an indicator for petroleum prospecting. The important microbial metabolic processes in most of the oil reservoir are sulfate reduction, fermentation, acetogenesis, methanogenesis, NO₃^- reduction, and Fe (III) and Mn (IV) reduction. The microorganisms residing in these sites have critical control on petroleum composition, recovery, and production methods. Physical characteristics of heavy oil are altered by microbial biotransformation and biosurfactant production. Considering oil to be one of the most vital energy resources, it is important to have a comprehensive understanding of petroleum microbiology. This manuscript reviews the recent research work referring to the diversity of bacteria in oil field and reservoir sites and their applications for enhancing oil transformation in the target reservoir and geomicrobial prospecting scope for petroleum exploration.

Noncoding-RNA-Mediated Regulation in Response to Macronutrient Stress in Plants

Macronutrient elements including nitrogen (N), phosphorus (P), potassium (K), calcium (Ca), magnesium (Mg), and sulfur (S) are required in relatively large and steady amounts for plant growth and development. Deficient or excessive supply of macronutrients from external environments may trigger a series of plant responses at phenotypic and molecular levels during the entire life cycle. Among the intertwined molecular networks underlying plant responses to macronutrient stress, noncoding RNAs (ncRNAs), mainly microRNAs (miRNAs) and long ncRNAs (lncRNAs), may serve as pivotal regulators for the coordination between nutrient supply and plant demand, while the responsive ncRNA-target module and the interactive mechanism vary among elements and species. Towards a comprehensive identification and functional characterization of nutrient-responsive ncRNAs and their downstream molecules, high-throughput sequencing has produced massive omics data for comparative expression profiling as a first step. In this review, we highlight the recent findings of ncRNA-mediated regulation in response to macronutrient stress, with special emphasis on the large-scale sequencing efforts for screening out candidate nutrient-responsive ncRNAs in plants, and discuss potential improvements in theoretical study to provide better guidance for crop breeding practices.

Comparative Transcriptome Analysis Reveals Sex-Based Differences during the Development of the Adult Parasitic Wasp Cotesia vestalis (Hymenoptera: Braconidae)

The endoparasitic wasp Cotesia vestalis is an important biological agent for controlling the population of Plutella xylostella, a major pest of cruciferous crops worldwide. Though the genome of C. vestalis has recently been reported, molecular mechanisms associated with sexual development have not been comprehensively studied. Here, we combined PacBio Iso-Seq and Illumina RNA-Seq to perform genome-wide profiling of pharate adult and adult development of male and female C. vestalis. Taking advantage of Iso-Seq full-length reads, we identified 14,466 novel transcripts as well as 8770 lncRNAs, with many lncRNAs showing a sex- and stage-specific expression pattern. The differentially expressed gene (DEG) analyses showed 2125 stage-specific and 326 sex-specific expressed genes. We also found that 4819 genes showed 11,856 alternative splicing events through combining the Iso-Seq and RNA-Seq data. The results of comparative analyses showed that most genes were alternatively spliced across developmental stages, and alternative splicing (AS) events were more prevalent in females than in males. Furthermore, we identified six sex-determining genes in this parasitic wasp and verified their sex-specific alternative splicing profiles. Specifically, the characterization of feminizer and doublesex splicing between male and female implies a conserved regulation mechanism of sexual development in parasitic wasps.

Comparative analysis of long noncoding RNAs in angiosperms and characterization of long noncoding RNAs in response to heat stress in Chinese cabbage

Long noncoding RNAs (lncRNAs) are widely present in different species and play critical roles in response to abiotic stresses. However, the functions of lncRNAs in Chinese cabbage under heat stress remain unknown. Here, we first conducted a global comparative analysis of 247,242 lncRNAs among 37 species. The results indicated that lncRNAs were poorly conserved among different species, and only 960 lncRNAs were homologous to 524 miRNA precursors. We then carried out lncRNA sequencing for a genome-wide analysis of lncRNAs and their target genes in Chinese cabbage at different stages of heat treatment. In total, 18,253 lncRNAs were identified, of which 1229 differentially expressed (DE) lncRNAs were characterized as being heat-responsive. The ceRNA network revealed that 38 lncRNAs, 16 miRNAs, and 167 mRNAs were involved in the heat response in Chinese cabbage. Combined analysis of the cis- and trans-regulated genes indicated that the targets of DE lncRNAs were significantly enriched in the "protein processing in endoplasmic reticulum" and "plant hormone signal transduction" pathways. Furthermore, the majority of HSP and PYL genes involved in these two pathways exhibited similar expression patterns and responded to heat stress rapidly. Based on the networks of DE lncRNA-mRNAs, 29 and 22 lncRNAs were found to interact with HSP and PYL genes, respectively. Finally, the expression of several critical lncRNAs and their targets was verified by qRT-PCR. Overall, we conducted a comparative analysis of lncRNAs among 37 species and performed a comprehensive analysis of lncRNAs in Chinese cabbage. Our findings expand the knowledge of lncRNAs involved in the heat stress response in Chinese cabbage, and the identified lncRNAs provide an abundance of resources for future comparative and functional studies.

PLncDB V2.0: a comprehensive encyclopedia of plant long noncoding RNAs

Long noncoding RNAs (lncRNAs) are transcripts longer than 200 nucleotides with little or no protein coding potential. The expanding list of lncRNAs and accumulating evidence of their functions in plants have necessitated the creation of a comprehensive database for lncRNA research. However, currently available plant lncRNA databases have some deficiencies, including the lack of lncRNA data from some model plants, uneven annotation standards, a lack of visualization for expression patterns, and the absence of epigenetic information. To overcome these problems, we upgraded our Plant Long noncoding RNA Database (PLncDB, http://plncdb.tobaccodb.org/), which was based on a uniform annotation pipeline. PLncDB V2.0 currently contains 1 246 372 lncRNAs for 80 plant species based on 13 834 RNA-Seq datasets, integrating lncRNA information from four other resources including EVLncRNAs, RNAcentral and etc. Expression patterns and epigenetic signals can be visualized using multiple tools (JBrowse, eFP Browser and EPexplorer). Targets and regulatory networks for lncRNAs are also provided for function exploration. In addition, PLncDB V2.0 is hierarchical and user-friendly and has five built-in search engines. We believe PLncDB V2.0 is useful for the plant lncRNA community and data mining studies and provides a comprehensive resource for data-driven lncRNA research in plants.

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.

Roles of Non-Coding RNAs in Response to Nitrogen Availability in Plants

Nitrogen (N) is an essential nutrient for plant growth and development; therefore, N deficiency is a major limiting factor in crop production. Plants have evolved mechanisms to cope with N deficiency, and the role of protein-coding genes in these mechanisms has been well studied. In the last decades, regulatory non-coding RNAs (ncRNAs), such as microRNAs (miRNAs), small interfering RNAs (siRNAs), and long ncRNAs (lncRNAs), have emerged as important regulators of gene expression in diverse biological processes. Recent advances in technologies for transcriptome analysis have enabled identification of N-responsive ncRNAs on a genome-wide scale. Characterization of these ncRNAs is expected to improve our understanding of the gene regulatory mechanisms of N response. In this review, we highlight recent progress in identification and characterization of N-responsive ncRNAs in Arabidopsis thaliana and several other plant species including maize, rice, and Populus.

The PacBio Full-Length Transcriptome of the Tea Aphid as a Reference Resource

The tea aphid, Aphis aurantii, has become one of the destructive pests in tea plantations in the tropics and subtropics. Very few functional studies have so far focused on the developmental and reproductive biology at a molecular level, because of the lack of comprehensive genetic information. Full-length transcriptomes represent a very highly efficient approach to obtain reference gene sequences in non-model insects. In the present study, the transcriptome of A. aurantii was comprehensively sequenced using PacBio Iso-Seq technology. A total of 46.8 Gb nucleotides and 15,938 non-redundant full-length transcripts were obtained, 13,498 (84.69%) of which were annotated into seven databases. Of these transcripts, 2,029 alternative splicing events and 15,223 simple sequence repeats were detected. Among these transcripts, 4,571 (28.68%) and 11,367 (71.32%) were long non-coding RNAs (lncRNAs) and protein-coding genes, respectively. Five hundred and ninety transcription factors were detected. The first full-length transcriptome represents a significant increase in the known genetic information of A. aurantii. It will assist the future functional study of genes involved in its development and reproduction.

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.

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.

Large-Scale Discovery of Non-conventional Peptides in Maize and Arabidopsis through an Integrated Peptidogenomic Pipeline

Non-conventional peptides (NCPs), which include small open reading frame-encoded peptides, play critical roles in fundamental biological processes. In this study, we developed an integrated peptidogenomic pipeline using high-throughput mass spectra to probe a customized six-frame translation database and applied it to large-scale identification of NCPs in plants.A total of 1993 and 1860 NCPs were unambiguously identified in maize and Arabidopsis, respectively. These NCPs showed distinct characteristics compared with conventional peptides and were derived from introns, 3' UTRs, 5' UTRs, junctions, and intergenic regions. Furthermore, our results showed that translation events in unannotated transcripts occur more broadly than previously thought. In addition, we found that dozens of maize NCPs are enriched within regions associated with phenotypic variations and domestication selection, indicating that they potentially are involved in genetic regulation of complex traits and domestication in maize. Taken together, our study developed an integrated peptidogenomic pipeline for large-scale identification of NCPs in plants, which would facilitate global characterization of NCPs from other plants. The identification of large-scale NCPs in both monocot (maize) and dicot (Arabidopsis) plants indicates that a large portion of plant genome can be translated into biologically functional molecules, which has important implications for functional genomic studies.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Identification and characterization of long non-coding RNAs as competing endogenous RNAs in the cold stress response of Triticum aestivum

Long non-coding RNAs (lncRNAs) play important roles in plant development and stress responses. MicroRNAs (miRNAs) are involved in transcriptional and post-transcriptional gene regulation. It is not clear how lncRNA-mediated plant responses to cold stress and how lncRNAs, miRNAs and target mRNAs cooperate subject to the competing endogenous RNA (ceRNA). We interpreted the function of lncRNAs in the winter wheat cultivar Dongnongdongmai 1 (Dn1). A total of 9970 putative lncRNAs were initially identified from three Dn1 lncRNA libraries (5 °C, -10 °C and -25 °C) using high-throughput sequencing. Among the 14,626 genes detected via weighted gene co-expression network analysis, 7435 lncRNAs were co-expressed with 7191 mRNAs. We found six modules related to cold resistance in the lncRNA-mRNA weighted co-expression network, and the functions of mRNAs were similar in each module. Antioxidant systems and hormones played important roles in low-temperature responses. RNA sequencing analysis revealed that interactions between the 384 lncRNAs and 70 miRNAs were required for ceRNA activity. According to ceRNA activity, 225 lncRNAs, 60 miRNAs and 621 target mRNAs were involved in the regulatory networks of the cold stress response. Notably, a conserved region was found in the complementary regions of lncRNAs and miR164/408 but had reverse expression trends in the ceRNA network. Our results reveal possible roles of lncRNAs-mRNAs in the regulatory networks associated with tolerance to low temperature and provide useful information for more strategic use of genomic resources in wheat breeding.

High-resolution profile of transcriptomes reveals a role of alternative splicing for modulating response to nitrogen in maize

Background

The fluctuation of nitrogen (N) contents profoundly affects the root growth and architecture in maize by altering the expression of thousands of genes. The differentially expressed genes (DEGs) in response to N have been extensively reported. However, information about the effects of N variation on the alternative splicing in genes is limited.

Results

To reveal the effects of N on the transcriptome comprehensively, we studied the N-starved roots of B73 in response to nitrate treatment, using a combination of short-read sequencing (RNA-seq) and long-read sequencing (PacBio-sequencing) techniques. Samples were collected before and 30 min after nitrate supply. RNA-seq analysis revealed that the DEGs in response to N treatment were mainly associated with N metabolism and signal transduction. In addition, we developed a workflow that utilizes the RNA-seq data to improve the quality of long reads, increasing the number of high-quality long reads to about 2.5 times. Using this workflow, we identified thousands of novel isoforms; most of them encoded the known functional domains and were supported by the RNA-seq data. Moreover, we found more than 1000 genes that experienced AS events specifically in the N-treated samples, most of them were not differentially expressed after nitrate supply-these genes mainly related to immunity, molecular modification, and transportation. Notably, we found a transcription factor ZmNLP6, a homolog of AtNLP7-a well-known regulator for N-response and root growth-generates several isoforms varied in capacities of activating downstream targets specifically after nitrate supply. We found that one of its isoforms has an increased ability to activate downstream genes. Overlaying DEGs and DAP-seq results revealed that many putative targets of ZmNLP6 are involved in regulating N metabolism, suggesting the involvement of ZmNLP6 in the N-response.

Conclusions

Our study shows that many genes, including the transcription factor ZmNLP6, are involved in modulating early N-responses in maize through the mechanism of AS rather than altering the transcriptional abundance. Thus, AS plays an important role in maize to adapt N fluctuation.

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.

GFF Utilities: GffRead and GffCompare

GTF (Gene Transfer Format) and GFF (General Feature Format) are popular file formats used by bioinformatics programs to represent and exchange information about various genomic features, such as gene and transcript locations and structure. GffRead and GffCompare are open source programs that provide extensive and efficient solutions to manipulate files in a GTF or GFF format. While GffRead can convert, sort, filter, transform, or cluster genomic features, GffCompare can be used to compare and merge different gene annotations. Availability and implementation: GFF utilities are implemented in C++ for Linux and OS X and released as open source under an MIT license  ( https://github.com/gpertea/gffread, https://github.com/gpertea/gffcompare).

GFF Utilities: GffRead and GffCompare

Summary: GTF (Gene Transfer Format) and GFF (General Feature Format) are popular file formats used by bioinformatics programs to represent and exchange information about various genomic features, such as gene and transcript locations and structure. GffRead and GffCompare are open source programs that provide extensive and efficient solutions to manipulate files in a GTF or GFF format. While GffRead can convert, sort, filter, transform, or cluster genomic features, GffCompare can be used to compare and merge different gene annotations.

Availability and implementation: GFF utilities are implemented in C++ for Linux and OS X and released as open source under an MIT license  ( https://github.com/gpertea/gffread, https://github.com/gpertea/gffcompare).

The Emerging Role of Long Non-Coding RNAs in Plant Defense Against Fungal Stress

Growing interest and recent evidence have identified long non-coding RNA (lncRNA) as the potential regulatory elements for eukaryotes. LncRNAs can activate various transcriptional and post-transcriptional events that impact cellular functions though multiple regulatory functions. Recently, a large number of lncRNAs have also been identified in higher plants, and an understanding of their functional role in plant resistance to infection is just emerging. Here, we focus on their identification in crop plant, and discuss their potential regulatory functions and lncRNA-miRNA-mRNA network in plant pathogen stress responses, referring to possible examples in a model plant. The knowledge gained from a deeper understanding of this colossal special group of plant lncRNAs will help in the biotechnological improvement of crops.

Cross-Species Root Transcriptional Network Analysis Highlights Conserved Modules in Response to Nitrate between Maize and Sorghum

Plants have evolved complex mechanisms to respond to the fluctuation of available nitrogen (N) in soil, but the genetic mechanisms underlying the N response in crops are not well-documented. In this study, we generated a time series of NO₃⁻-mediated transcriptional profiles in roots of maize and sorghum, respectively. Using weighted gene co-expression network analysis, we identified modules of co-expressed genes that related to NO₃⁻ treatments. A cross-species comparison revealed 22 conserved modules, of which four were related to hormone signaling, suggesting that hormones participate in the early nitrate response. Three other modules are composed of genes that are mainly upregulated by NO₃⁻ and involved in nitrogen and carbohydrate metabolism, including NRT, NIR, NIA, FNR, and G6PD2. Two G2-like transcription factors (ZmNIGT1 and SbNIGT1), induced by NO₃⁻ stimulation, were identified as hub transcription factors (TFs) in the modules. Transient assays demonstrated that ZmNIGT1 and SbNIGT1 are transcriptional repressors. We identified the target genes of ZmNIGT1 by DNA affinity-purification sequencing (DAP-Seq) and found that they were significantly enriched in catalytic activity, including carbon, nitrogen, and other nutrient metabolism. A set of ZmNIGT1 targets encode transcription factors (ERF, ARF, and AGL) that are involved in hormone signaling and root development. We propose that ZmNIGT1 and SbNIGT1 are negative regulators of nitrate responses that play an important role in optimizing nutrition metabolism and root morphogenesis. Together with conserved N responsive modules, our study indicated that, to encounter N variation in soil, maize and sorghum have evolved an NO₃⁻-regulatory network containing a set of conserved modules and transcription factors.

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.

GFF Utilities: GffRead and GffCompare

GTF (Gene Transfer Format) and GFF (General Feature Format) are popular file formats used by bioinformatics programs to represent and exchange information about various genomic features, such as gene and transcript locations and structure. GffRead and GffCompare are open source programs that provide extensive and efficient solutions to manipulate files in a GTF or GFF format. While GffRead can convert, sort, filter, transform, or cluster genomic features, GffCompare can be used to compare and merge different gene annotations. Availability and implementation: GFF utilities are implemented in C++ for Linux and OS X and released as open source under an MIT license  ( https://github.com/gpertea/gffread, https://github.com/gpertea/gffcompare).

Identification, characterization, and functional prediction of circular RNAs in maize

Circular RNAs (circRNAs) are a new type of intracellular regulator that have been widely identified in animals and plants by high-throughput sequencing. However, there are still few functional studies on circRNAs in plants. To better understand maize circRNAs and their potential functions, we identified 1199 circRNAs in maize from RiboMinus RNA-Seq transcriptome data, and found distinct features of splicing site selection bias, longer flanking introns, and miniature inverted-repeat transposable element (MITE) insertions in flanking introns in maize circRNAs compared to other plant circRNAs. In total, 31 and 36 orthologous circRNAs were identified in rice and maize, respectively, but the orthologous parental genes could not produce orthologous circRNAs, mostly because of long-sequence insertions/deletions at flanking introns and approximately 24.3% of them contained MITE sequences. The majority of maize circRNAs showed high diversity of expression under different treatments and/or in different genetic backgrounds, implying that circRNAs could be involved in various regulatory networks. Twenty-six ecircRNAs were predicted to contain one or more target mimics, and 229 circRNAs had high coding potential, indicating that circRNAs could perform peptide-encoding functions in plants. These results will broaden understanding of the roles of circRNAs in plants and support further functional work on maize.

Long Non-coding RNA in Plants in the Era of Reference Sequences

The discovery of non-coding RNAs (ncRNAs), and the subsequent elucidation of their functional roles, was largely delayed due to the misidentification of non-protein-coding parts of DNA as "junk DNA," which forced ncRNAs into the shadows of their protein-coding counterparts. However, over the past decade, insight into the important regulatory roles of ncRNAs has led to rapid progress in their identification and characterization. Of the different types of ncRNAs, long non-coding RNAs (lncRNAs), has attracted considerable attention due to their mRNA-like structures and gene regulatory functions in plant stress responses. While RNA sequencing has been commonly used for mining lncRNAs, a lack of widespread conservation at the sequence level in addition to relatively low and highly tissue-specific expression patterns challenges high-throughput in silico identification approaches. The complex folding characteristics of lncRNA molecules also complicate target predictions, as the knowledge about the interaction interfaces between lncRNAs and potential targets is insufficient. Progress in characterizing lncRNAs and their targets from different species may hold the key to efficient identification of this class of ncRNAs from transcriptomic and potentially genomic resources. In wheat and barley, two of the most important crops, the knowledge about lncRNAs is very limited. However, recently published high-quality genomes of these crops are considered as promising resources for the identification of not only lncRNAs, but any class of molecules. Considering the increasing demand for food, these resources should be used efficiently to discover molecular mechanisms lying behind development and a/biotic stress responses. As our understanding of lncRNAs expands, interactions among ncRNA classes, as well as interactions with the coding sequences, will likely define novel functional networks that may be modulated for crop improvement.

Novel Aspects of Nitrate Regulation in Arabidopsis

Nitrogen (N) is one of the most essential macronutrients for plant growth and development. Nitrate (NO₃ ^-), the major form of N that plants uptake from the soil, acts as an important signaling molecule in addition to its nutritional function. Over the past decade, significant progress has been made in identifying new components involved in NO₃ ^- regulation and starting to unravel the NO₃ ^- regulatory network. Great reviews have been made recently by scientists on the key regulators in NO₃ ^- signaling, NO₃ ^- effects on plant development, and its crosstalk with phosphorus (P), potassium (K), hormones, and calcium signaling. However, several novel aspects of NO₃ ^- regulation have not been previously reviewed in detail. Here, we mainly focused on the recent advances of post-transcriptional regulation and non-coding RNA (ncRNAs) in NO₃ ^- signaling, and NO₃ ^- regulation on leaf senescence and the circadian clock. It will help us to extend the general picture of NO₃ ^- regulation and provide a basis for further exploration of NO₃ ^- regulatory network.