Transcriptomic analyses of rice (Oryza sativa) genes and non-coding RNAs under nitrogen starvation using multiple omics technologies

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.

Underground communication: Long non-coding RNA signaling in the plant rhizosphere

Long non-coding RNAs (lncRNAs) have emerged as integral gene-expression regulators underlying plant growth, development, and adaptation. To adapt to the heterogeneous and dynamic rhizosphere, plants use interconnected regulatory mechanisms to optimally fine-tune gene-expression-governing interactions with soil biota, as well as nutrient acquisition and heavy metal tolerance. Recently, high-throughput sequencing has enabled the identification of plant lncRNAs responsive to rhizosphere biotic and abiotic cues. Here, we examine lncRNA biogenesis, classification, and mode of action, highlighting the functions of lncRNAs in mediating plant adaptation to diverse rhizosphere factors. We then discuss studies that reveal the significance and target genes of lncRNAs during developmental plasticity and stress responses at the rhizobium interface. A comprehensive understanding of specific lncRNAs, their regulatory targets, and the intricacies of their functional interaction networks will provide crucial insights into how these transcriptomic switches fine-tune responses to shifting rhizosphere signals. Looking ahead, we foresee that single-cell dissection of cell-type-specific lncRNA regulatory dynamics will enhance our understanding of the precise developmental modulation mechanisms that enable plant rhizosphere adaptation. Overcoming future challenges through multi-omics and genetic approaches will more fully reveal the integral roles of lncRNAs in governing plant adaptation to the belowground environment.

Switching action modes of miR408-5p mediates auxin signaling in rice

MicroRNAs (miRNAs) play fundamental roles in many developmental and physiological processes in eukaryotes. MiRNAs in plants generally regulate their targets via either mRNA cleavage or translation repression; however, which approach plays a major role and whether these two function modes can shift remains elusive. Here, we identify a miRNA, miR408-5p that regulates AUXIN/INDOLE ACETIC ACID 30 (IAA30), a critical repressor in the auxin pathway via switching action modes in rice. We find that miR408-5p usually inhibits IAA30 protein translation, but in a high auxin environment, it promotes the decay of IAA30 mRNA when it is overproduced. We further demonstrate that IDEAL PLANT ARCHITECTURE1 (IPA1), an SPL transcription factor regulated by miR156, mediates leaf inclination through association with miR408-5p precursor promoter. We finally show that the miR156-IPA1-miR408-5p-IAA30 module could be controlled by miR393, which silences auxin receptors. Together, our results define an alternative auxin transduction signaling pathway in rice that involves the switching of function modes by miR408-5p, which contributes to a better understanding of the action machinery as well as the cooperative network of miRNAs in plants.

CcNFYB3-CcMATE35 and LncRNA CcLTCS-CcCS modules jointly regulate the efflux and synthesis of citrate to enhance aluminium tolerance in pigeon pea

Aluminium (Al) toxicity decreases crop production in acid soils in general, but many crops have evolved complex mechanisms to resist it. However, our current understanding of how plants cope with Al stress and perform Al resistance is still at the initial stage. In this study, the citrate transporter CcMATE35 was identified to be involved in Al stress response. The release of citrate was increased substantially in CcMATE35 over-expression (OE) lines under Al stress, indicating enhanced Al resistance. It was demonstrated that transcription factor CcNFYB3 regulated the expression of CcMATE35, promoting the release of citrate from roots to increase Al resistance in pigeon pea. We also found that a Long noncoding RNA Targeting Citrate Synthase (CcLTCS) is involved in Al resistance in pigeon pea. Compared with controls, overexpression of CcLTCS elevated the expression level of the Citrate Synthase gene (CcCS), leading to increases in root citrate level and citrate release, which forms another module to regulate Al resistance in pigeon pea. Simultaneous overexpression of CcNFYB3 and CcLTCS further increased Al resistance. Taken together, these findings suggest that the two modules, CcNFYB3-CcMATE35 and CcLTCS-CcCS, jointly regulate the efflux and synthesis of citrate and may play an important role in enhancing the resistance of pigeon pea under Al stress.

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

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

The rice NUCLEAR FACTOR-YA5 and MICRORNA169a module promotes nitrogen utilization during nitrogen deficiency

Nitrogen (N) is essential for plant growth and development. Therefore, understanding its utilization is essential for improving crop productivity. However, much remains to be learned about plant N sensing and signaling. Here, rice (Oryza sativa) NUCLEAR FACTOR-YA5 (OsNF-YA5) expression was tightly regulated by N status and induced under N-deficient conditions. Overexpression (OE) of OsNF-YA5 in rice resulted in increased chlorophyll levels and delayed senescence compared to control plants under normal N conditions. Agronomic traits were significantly improved in OE plants and impaired in knockout mutants under N-deficient conditions. Using a dexamethasone-inducible system, we identified the putative targets of OsNF-YA5 that include amino acid, nitrate/peptide transporters, and NITRATE TRANSPORTER 1.1A (OsNRT1.1A), which functions as a key transporter in rice. OsNF-YA5 directly enhanced OsNRT1.1A expression and N uptake rate under N-deficient conditions. Besides, overexpression of OsNF-YA5 also enhanced the expression of GLUTAMINE SYNTHETASE 1/2 (GS1/2) and GLUTAMINE OXOGLUTARATE AMINOTRANSFERASE 1/2 (GOGAT1/2), increasing free amino acid contents under N-deficient conditions. Osa-miR169a expression showed an opposite pattern with OsNF-YA5 depending on N status. Further analysis revealed that osa-miR169a negatively regulates OsNF-YA5 expression and N utilization, demonstrating that an OsNF-YA5/osa-miR169a module tightly regulates rice N utilization for adaptation to N status.

Nitrate-responsive transcriptome analysis of rice RGA1 mutant reveals the role of G-protein alpha subunit in negative regulation of nitrogen-sensitivity and use efficiency

KEY MESSAGE

Nitrate-responsive transcriptomic, phenotypic and physiological analyses of rice RGA1 mutant revealed many novel RGA1-regulated genes/processes/traits related to nitrogen use efficiency, and provided robust genetic evidence of RGA1-regulation of NUE. Nitrogen (N) use efficiency (NUE) is important for sustainable agriculture. G-protein signalling was implicated in N-response/NUE in rice, but needed firm genetic characterization of the role of alpha subunit (RGA1). The knock-out mutant of RGA1 in japonica rice exhibited lesser nitrate-dose sensitivity than the wild type (WT), in yield and NUE. We, therefore, investigated its genomewide nitrate-response relative to WT. It revealed 3416 differentially expressed genes (DEGs), including 719 associated with development, grain yield and phenotypic traits for NUE. The upregulated DEGs were related to photosynthesis, chlorophyll, tetrapyrrole and porphyrin biosynthesis, while the downregulated DEGs belonged to cellular protein metabolism and transport, small GTPase signalling, cell redox homeostasis, etc. We validated 26 nitrate-responsive DEGs across functional categories by RT-qPCR. Physiological validation of nitrate-response in the mutant and the WT at 1.5 and 15 mM doses revealed higher chlorophyll and stomatal length but decreased stomatal density, conductance and transpiration. The consequent increase in photosynthesis and water use efficiency may have contributed to better yield and NUE in the mutant, whereas the WT was N-dose sensitive. The mutant was not as N-dose-responsive as the WT in shoot/root growth, productive tillers and heading date, but equally responsive as WT in total N and protein content. The RGA1 mutant was less impacted by higher N-dose or salt stress in terms of yield, protein content, photosynthetic performance, relative water content, water use efficiency and catalase activity. PPI network analyses revealed known NUE-related proteins as RGA1 interactors. Therefore, RGA1 negatively regulates N-dose sensitivity and NUE in rice.

Transcriptome analysis reveals a lncRNA-miRNA-mRNA regulatory network in OsRpp30-mediated disease resistance in rice

BACKGROUND

Long non-coding RNAs (lncRNAs) play critical roles in various biological processes in plants. Extensive studies utilizing high-throughput RNA sequencing have revealed that many lncRNAs are involved in plant disease resistance. Oryza sativa RNase P protein 30 (OsRpp30) has been identified as a positive regulator of rice immunity against fungal and bacterial pathogens. Nevertheless, the specific functions of lncRNAs in relation to OsRpp30-mediated disease resistance in rice remain elusive.

RESULTS

We conducted a comprehensive analysis of lncRNAs, miRNAs, and mRNAs expression patterns in wild type (WT), OsRpp30 overexpression (OsRpp30-OE), and OsRpp30 knockout (OsRpp30-KO) rice plants. In total, we identified 91 differentially expressed lncRNAs (DElncRNAs), 1671 differentially expressed mRNAs (DEmRNAs), and 41 differentially expressed miRNAs (DEmiRNAs) across the different rice lines. To gain further insights, we investigated the interaction between DElncRNAs and DEmRNAs, leading to the discovery of 10 trans- and 27 cis-targeting pairs specific to the OsRpp30-OE and OsRpp30-KO samples. In addition, we constructed a competing endogenous RNA (ceRNA) network comprising differentially expressed lncRNAs, miRNAs, and mRNAs to elucidate their intricate interplay in rice disease resistance. The ceRNA network analysis uncovered a set of gene targets regulated by lncRNAs and miRNAs, which were found to be involved in pathogen recognition, hormone pathways, transcription factor activation, and other biological processes related to plant immunity.

CONCLUSIONS

Our study provides a comprehensive expression profiling of lncRNAs, miRNAs, and mRNAs in a collection of defense mutants in rice. To decipher the putative functional significance of lncRNAs, we constructed trans- and cis-targeting networks involving differentially expressed lncRNAs and mRNAs, as well as a ceRNA network incorporating differentially expressed lncRNAs, miRNAs, and mRNAs. Together, the findings from this study provide compelling evidence supporting the pivotal roles of lncRNAs in OsRpp30-mediated disease resistance in rice.

Evaluation of efficacy of non-coding RNA in abiotic stress management of field crops: Current status and future prospective

Abiotic stresses are responsible for the major losses in crop yield all over the world. Stresses generate harmful ROS which can impair cellular processes in plants. Therefore, plants have evolved antioxidant systems in defence against the stress-induced damages. The frequency of occurrence of abiotic stressors has increased several-fold due to the climate change experienced in recent times and projected for the future. This had particularly aggravated the risk of yield losses and threatened global food security. Non-coding RNAs are the part of eukaryotic genome that does not code for any proteins. However, they have been recently found to have a crucial role in the responses of plants to both abiotic and biotic stresses. There are different types of ncRNAs, for example, miRNAs and lncRNAs, which have the potential to regulate the expression of stress-related genes at the levels of transcription, post-transcription, and translation of proteins. The lncRNAs are also able to impart their epigenetic effects on the target genes through the alteration of the status of histone modification and organization of the chromatins. The current review attempts to deliver a comprehensive account of the role of ncRNAs in the regulation of plants' abiotic stress responses through ROS homeostasis. The potential applications ncRNAs in amelioration of abiotic stresses in field crops also have been evaluated.

Non-Targeted Metabolomic Analysis of Arabidopsis thaliana (L.) Heynh: Metabolic Adaptive Responses to Stress Caused by N Starvation

As sessile organisms, plants develop the ability to respond and survive in changing environments. Such adaptive responses maximize phenotypic and metabolic fitness, allowing plants to adjust their growth and development. In this study, we analyzed the metabolic plasticity of Arabidopsis thaliana in response to nitrate deprivation by untargeted metabolomic analysis and using wild-type (WT) genotypes and the loss-of-function nia1/nia2 double mutant. Secondary metabolites were identified using seedlings grown on a hydroponic system supplemented with optimal or limiting concentrations of N (4 or 0.2 mM, respectively) and harvested at 15 and 30 days of age. Then, spectral libraries generated from shoots and roots in both ionization modes (ESI +/-) were compared. Totals of 3407 and 4521 spectral signals (m/z_rt) were obtained in the ESI+ and ESI- modes, respectively. Of these, approximately 50 and 65% were identified as differentially synthetized/accumulated. This led to the presumptive identification of 735 KEGG codes (metabolites) belonging to 79 metabolic pathways. The metabolic responses in the shoots and roots of WT genotypes at 4 mM of N favor the synthesis/accumulation of metabolites strongly related to growth. In contrast, for the nia1/nia2 double mutant (similar as the WT genotype at 0.2 mM N), metabolites identified as differentially synthetized/accumulated help cope with stress, regulating oxidative stress and preventing programmed cell death, meaning that metabolic responses under N starvation compromise growth to prioritize a defensive response.

Comprehensive identification and characterization of lncRNAs and circRNAs reveal potential brown planthopper-responsive ceRNA networks in rice

Brown planthopper (Nilaparvata lugens Stål, BPH) is one of the most destructive pests of rice. Non-coding RNA plays an important regulatory role in various biological processes. However, comprehensive identification and characterization of long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) in BPH-infested rice have not been performed. Here, we performed a genome-wide analysis of lncRNAs and circRNAs in BPH6-transgenic (resistant, BPH6G) and Nipponbare (susceptible, NIP) rice plants before and after BPH feeding (early and late stage) via deep RNA-sequencing. A total of 310 lncRNAs and 129 circRNAs were found to be differentially expressed. To reveal the different responses of resistant and susceptible rice to BPH herbivory, the potential functions of these lncRNAs and circRNAs as competitive endogenous RNAs (ceRNAs) were predicted and investigated using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses. Dual-luciferase reporter assays revealed that miR1846c and miR530 were targeted by the lncRNAs XLOC_042442 and XLOC_028297, respectively. In responsive to BPH infestation, 39 lncRNAs and 21 circRNAs were predicted to combine with 133 common miRNAs and compete for miRNA binding sites with 834 mRNAs. These mRNAs predictably participated in cell wall organization or biogenesis, developmental growth, single-organism cellular process, and the response to stress. This study comprehensively identified and characterized lncRNAs and circRNAs, and integrated their potential ceRNA functions, to reveal the rice BPH-resistance network. These results lay a foundation for further study on the functions of lncRNAs and circRNAs in the rice-BPH interaction, and enriched our understanding of the BPH-resistance response in rice.

Pokkali: A Naturally Evolved Salt-Tolerant Rice Shows a Distinguished Set of lncRNAs Possibly Contributing to the Tolerant Phenotype

Pokkali is a strong representation of how stress-tolerant genotypes have evolved due to natural selection pressure. Numerous omics-based investigations have indicated different categories of stress-related genes and proteins, possibly contributing to salinity tolerance in this wild rice. However, a comprehensive study towards understanding the role of long-noncoding RNAs (lncRNAs) in the salinity response of Pokkali has not been done to date. We have identified salt-responsive lncRNAs from contrasting rice genotypes IR64 and Pokkali. A total of 63 and 81 salinity-responsive lncRNAs were differentially expressed in IR64 and Pokkali, respectively. Molecular characterization of lncRNAs and lncRNA-miRNA-mRNA interaction networks helps to explore the role of lncRNAs in the stress response. Functional annotation revealed that identified lncRNAs modulate various cellular processes, including transcriptional regulation, ion homeostasis, and secondary metabolite production. Additionally, lncRNAs were predicted to bind stress-responsive transcription factors, namely ERF, DOF, and WRKY. In addition to salinity, expression profiling was also performed under other abiotic stresses and phytohormone treatments. A positive modulation in TCONS_00035411, TCONS_00059828, and TCONS_00096512 under both abiotic stress and phytohormone treatments could be considered as being of potential interest for the further functional characterization of IncRNA. Thus, extensive analysis of lncRNAs under various treatments helps to delineate stress tolerance mechanisms and possible cross-talk.

Transcriptomic and Physiological Analyses of Two Rice Restorer Lines under Different Nitrogen Supplies Provide Novel Insights into Hybrid Rice Breeding

Improving plant nitrogen-use efficiency (NUE) has great significance for various crops, particularly in hybrid breeding. Reducing nitrogen inputs is key to achieving sustainable rice production and mitigating environmental problems. In this study, we analyzed the transcriptomic and physiological changes in two indica restorer lines (Nanhui511 [NH511] and Minghui23 [MH23]) under high nitrogen (HN) and low nitrogen (LN) conditions. Compared to MH23, NH511 was more sensitive to different nitrogen supplies and exhibited higher nitrogen uptake and NUE under HN conditions by increasing lateral root and tiller numbers in the seedling and maturation stages, respectively. NH511 also exhibited a lower survival rate than MH23 when planted in a chlorate-containing hydroponic solution, indicating its HN uptake ability under different nitrogen-supply conditions. Transcriptomic analysis showed that NH511 has 2456 differentially expressed genes, whereas MH23 had only 266. Furthermore, these genes related to nitrogen utilization showed differential expression in NH511 under HN conditions, while the opposite was observed in MH23. Our findings revealed that NH511 could be regarded as elite rice and used for breeding high-NUE restorer lines by regulating and integrating nitrogen-utilization genes, which provides novel insights for the cultivation of high-NUE hybrid rice.

Expression levels of nitrogen assimilation-related genes, physiological responses, and morphological adaptations of three indica rice (Oryza sativa L. ssp. indica) genotypes subjected to nitrogen starvation conditions

Nitrogen (N) is an essential nutrient available to the plants in form of nitrate and ammonium. It is a macronutrient important for the plant growth and development, especially in cereal crops, which consume it for the production of amino acids, proteins/enzymes, nucleic acids, cell wall complexes, plant hormones, and vitamins. In rice production, 17 kg N uptake is required to produce 1 ton of rice. Considering this, many techniques have been developed to evaluate leaf greenness or SPAD value for assessing the amount of N application in the rice cultivar to maximize the grain yield. The aim of the present study was to investigate the morpho-physiological characteristics and relative expression level of N assimilation in three different rice genotypes (MT2, RD31, KDML105) under 1.00 × (full N), 0.50 × , 0.25 × (N depletion), and 0.00 × (N deficiency) at seedling stage and the morpho-physiological traits and the grain yield attributes under 1.00 × (full N) and 0.25 × (N depletion) were compared. Leaf chlorosis and growth inhibition in rice seedlings under N deficiency were evidently observed. Shoot height, number of leaves, shoot fresh weight, shoot dry weight, and root fresh weight in KDML105 under N deficiency were decreased by 27.65%, 42.11%, 65.44%, 47.90%, and 54.09% over the control (full N). Likewise, leaf greenness was lowest in KDML105 under N deficiency (78.57% reduction over the full N), leading to low photosynthetic abilities. In addition, expression of nitrogen assimilation-related genes, OsNR1, OsGln1;1, and OsGln2, in KDML105 under N depletion were increased within 3 h and then declined after the long incubation period, whereas those were unchanged in cvs. MT2 and RD31. Similarly, relative expression level of OsNADH-GOGAT, OsFd-GOGAT, and OsAspAt1 in KDML105 was peaked when subjected to 0.50 × N for 6 h and then declined after the long incubation period. Moreover, overall growth characters and physiological changes in cv. RD31 at vegetative stage under 0.25 × N were retained better than those in cvs. KDML105 and MT2, resulting in high yield at the harvesting process. In summary, N assimilated-related genes in rice seedlings under N depletion were rapidly regulated within 3-6 h, especially cv. KDML105 and MT2, then downregulated, resulting in physiological changes, growth inhibition, and yield reduction.

Comparative transcriptome analyses under individual and combined nutrient starvations provide insights into N/P/K interactions in rice

Crops often suffer from simultaneous limitations of multiple nutrients in soils, including nitrogen (N), phosphorus (P) and potassium (K), which are three major macronutrients essential for ensuring growth and yield. Although plant responses to individual N, P, and K deficiency have been well documented, our understanding of the responses to combined nutrient deficiencies and the crosstalk between nutrient starvation responses is still limited. Here, we compared the physiological responses in rice under seven kinds of single and multiple low nutrient stress of N, P and K, and used RNA sequencing approaches to compare their transcriptome changes. A total of 13,000 genes were found to be differentially expressed under all these single and multiple low N/P/K stresses, and 66 and 174 of them were shared by all these stresses in roots and shoots, respectively. Functional enrichment analyses of the DEGs showed that a group of biological and metabolic processes were shared by these low N/P/K stresses. Comparative analyses indicated that DEGs under multiple low nutrient stress was not the simple summation of single nutrient stress. N was found to be the predominant factor affecting the transcriptome under combined nutrient stress. N, P, or K availability exhibited massive influences on the transcriptomic responses to starvation of other nutrients. Many genes involved in nutrient transport, hormone signaling, and transcriptional regulation were commonly responsive to low N/P/K stresses. Some transcription factors were predicted to regulate the expression of genes that are commonly responsive to N, P, and K starvations. These results revealed the interactions between N, P, and K starvation responses, and will be helpful for further elucidation of the molecular mechanisms underlying nutrient interactions.

Long term nitrogen deficiency alters expression of miRNAs and alters nitrogen metabolism and root architecture in Indian dwarf wheat (Triticum sphaerococcum Perc.) genotypes

The important roles of plant microRNAs (miRNAs) in adaptation to nitrogen (N) deficiency in different crop species especially cereals (rice, wheat, maize) have been under discussion since last decade with little focus on potential wild relatives and landraces. Indian dwarf wheat (Triticum sphaerococcum Percival) is an important landrace native to the Indian subcontinent. Several unique features, especially high protein content and resistance to drought and yellow rust, make it a very potent landrace for breeding. Our aim in this study is to identify the contrasting Indian dwarf wheat genotypes based on nitrogen use efficiency (NUE) and nitrogen deficiency tolerance (NDT) traits and the associated miRNAs differentially expressed under N deficiency in selected genotypes. Eleven Indian dwarf wheat genotypes and a high NUE bread wheat genotype (for comparison) were evaluated for NUE under control and N deficit field conditions. Based on NUE, selected genotypes were further evaluated under hydroponics and miRNome was compared by miRNAseq under control and N deficit conditions. Among the identified, differentially expressed miRNAs in control and N starved seedlings, the target gene functions were associated with N metabolism, root development, secondary metabolism and cell-cycle associated pathways. The key findings on miRNA expression, changes in root architecture, root auxin abundance and changes in N metabolism reveal new information on the N deficiency response of Indian dwarf wheat and targets for genetic improvement of NUE.

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.

Transcriptome Profiling and Chlorophyll Metabolic Pathway Analysis Reveal the Response of Nitraria tangutorum to Increased Nitrogen

To identify genes that respond to increased nitrogen and assess the involvement of the chlorophyll metabolic pathway and associated regulatory mechanisms in these responses, Nitraria tangutorum seedlings were subjected to four nitrogen concentrations (N0, N6, N36, and N60: 0, 6, 36, and 60 mmol·L^-1 nitrogen, respectively). The N. tangutorum seedling leaf transcriptome was analyzed by high-throughput sequencing (Illumina HiSeq 4000), and 332,420 transcripts and 276,423 unigenes were identified. The numbers of differentially expressed genes (DEGs) were 4052 in N0 vs. N6, 6181 in N0 vs. N36, and 3937 in N0 vs. N60. Comparing N0 and N6, N0 and N36, and N0 and N60, we found 1101, 2222, and 1234 annotated DEGs in 113, 121, and 114 metabolic pathways, respectively, classified in the Kyoto Encyclopedia of Genes and Genomes database. Metabolic pathways with considerable accumulation were involved mainly in anthocyanin biosynthesis, carotenoid biosynthesis, porphyrin and chlorophyll metabolism, flavonoid biosynthesis, and amino acid metabolism. N36 increased δ-amino levulinic acid synthesis and upregulated expression of the magnesium chelatase H subunit, which promoted chlorophyll a synthesis. Hence, N36 stimulated chlorophyll synthesis rather than heme synthesis. These findings enrich our understanding of the N. tangutorum transcriptome and help us to research desert xerophytes' responses to increased nitrogen in the future.

Genome-wide analysis of long non-coding RNAs in sugar beet (Beta vulgaris L.) under drought stress

Drought stress is one of the most severe abiotic stresses that restrict global crop production. Long non-coding RNAs (lncRNAs) have been proved to play a key role in response to drought stress. However, genome-wide identification and characterization of drought-responsive lncRNAs in sugar beet is still lacking. Thus, the present study focused on analyzing lncRNAs in sugar beet under drought stress. We identified 32017 reliable lncRNAs in sugar beet by strand-specific high-throughput sequencing. A total of 386 differentially expressed lncRNAs (DElncRNAs) were found under drought stress. The most significantly upregulated and downregulated lncRNAs were TCONS_00055787 (upregulated by more than 6000 fold) and TCONS_00038334 (downregulated by more than 18000 fold), respectively. Quantitative real-time PCR results exhibited a high concordance with RNA sequencing data, which conformed that the expression patterns of lncRNAs based on RNA sequencing were highly reliable. In addition, we predicted 2353 and 9041 transcripts that were estimated to be the cis- and trans-target genes of the drought-responsive lncRNAs. As revealed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, the target genes of DElncRNAs were significantly enriched in organelle subcompartment, thylakoid, endopeptidase activity, catalytic activity, developmental process, lipid metabolic process, RNA polymerase activity, transferase activity, flavonoid biosynthesis and several other terms associated with abiotic stress tolerance. Moreover, 42 DElncRNAs were predicted as potential miRNA target mimics. LncRNAs have important effects on plant adaptation to drought conditions through the interaction with protein-encoding genes. The present study leads to greater insights into lncRNA biology and offers candidate regulators for improving the drought tolerance of sugar beet cultivars at the genetic level.

Contribution of RdDM to the ecotype-specific differential methylation on conserved as well as highly variable regions between Arabidopsis ecotypes

BACKGROUND

Several studies showed genome-wide DNA methylation during Arabidopsis embryogenesis and germination. Although it has been known that the change of DNA methylation mainly occurs at CHH context mediated by small RNA-directed DNA methylation pathway during seed ripening and germination, the causality of the methylation difference exhibited in natural Arabidopsis ecotypes has not been thoroughly studied.

RESULTS

In this study we compared DNA methylation difference using comparative pairwise multi-omics dynamics in Columbia-0 (Col) and Cape Verde Island (Cvi) ecotypes. Arabidopsis genome was divided into two regions, common regions in both ecotypes and Col-specific regions, depending on the reads mapping of whole genome bisulfite sequencing libraries from both ecotypes. Ecotype comparison was conducted within common regions and the levels of DNA methylation on common regions and Col-specific regions were also compared. we confirmed transcriptome were relatively dynamic in stage-wise whereas the DNA methylome and small RNAome were more ecotype-dependent. While the global CG methylation remains steady during maturation and germination, we found genic CG methylation differs the most between the two accessions. We also found that ecotype-specific differentially methylated regions (eDMR) are positively correlated with ecotype-specifically expressed 24-nt small RNA clusters. In addition, we discovered that Col-specific regions enriched with transposable elements (TEs) and structural variants that tend to become hypermethylated, and TEs in Col-specific regions were longer in size, more pericentromeric, and more hypermethylated than those in the common regions. Through the analysis of RdDM machinery mutants, we confirmed methylation on Col-specific region as well as on eDMRs in common region are contributed by RdDM pathway. Lastly, we demonstrated that highly variable sequences between ecotypes (HOT regions) were also affected by RdDM-mediated regulation.

CONCLUSIONS

Through ecotype comparison, we revealed differences and similarities of their transcriptome, methylome and small RNAome both in global and local regions. We validated the contribution of RdDM causing differential methylation of common regions. Hypermethylated ecotype-specific regions contributed by RNA-directed DNA methylation pathway largely depend on the presence of TEs and copy-gain structural variations. These ecotype-specific regions are frequently associated with HOT regions, providing evolutionary insights into the epigenome dynamics within a species.

Molecular and physiological mechanisms of tea (Camellia sinensis (L.) O. Kuntze) leaf and root in response to nitrogen deficiency

BACKGROUND

As an economically important crop, tea is strongly nitrogen (N)-dependent. However, the physiological and molecular mechanisms underlying the response of N deficiency in tea are not fully understood. Tea cultivar "Chunlv2" [Camellia sinensis (L.) O. Kuntze] were cultured with a nutrient solution with 0 mM [N-deficiency] or 3 mM (Control) NH₄NO₃ in 6 L pottery pots containing clean river sands.

RESULTS

N deficiency significantly decreased N content, dry weight, chlorophyll (Chl) content, L-theanine and the activities of N metabolism-related enzymes, but increased the content of total flavonoids and polyphenols in tea leaves. N deficiency delayed the sprouting time of tea buds. By using the RNA-seq technique and subsequent bioinformatics analysis, 3050 up-regulated and 2688 down-regulated differentially expressed genes (DEGs) were isolated in tea leaves in response to N deficiency. However, only 1025 genes were up-regulated and 744 down-regulated in roots. Gene ontology (GO) term enrichment analysis showed that 205 DEGs in tea leaves were enriched in seven GO terms and 152 DEGs in tea roots were enriched in 11 GO items based on P < 0.05. In tea leaves, most GO-enriched DEGs were involved in chlorophyll a/b binding activities, photosynthetic performance, and transport activities. But most of the DEGs in tea roots were involved in the metabolism of carbohydrates and plant hormones with regard to the GO terms of biological processes. N deficiency significantly increased the expression level of phosphate transporter genes, which indicated that N deficiency might impair phosphorus metabolism in tea leaves. Furthermore, some DEGs, such as probable anion transporter 3 and high-affinity nitrate transporter 2.7, might be of great potential in improving the tolerance of N deficiency in tea plants and further study could work on this area in the future.

CONCLUSIONS

Our results indicated N deficiency inhibited the growth of tea plant, which might be due to altered N metabolism and expression levels of DEGs involved in the photosynthetic performance, transport activity and oxidation-reduction processes.

Ectopic Expression of Sugarcane ScAMT1.1 Has the Potential to Improve Ammonium Assimilation and Grain Yield in Transgenic Rice under Low Nitrogen Stress

In China, nitrogen (N) fertilizer is excessively used in sugarcane planting areas, while the nitrogen use efficiency (NUE) of sugarcane is relatively low. Mining and identifying the key genes in response to low N stress in sugarcane can provide useful gene elements and a theoretical basis for developing sugarcane varieties with high NUE. In our study, RNA-Seq combined with qRT-PCR analysis revealed that the ScAMT1.1 gene responded positively to low N stress, resulting in the stronger low N tolerance and high NUE ability of sugarcane cultivar ROC22. Then, ScAMT1.1 was cloned from sugarcane. The full-length cDNA of the ScAMT1.1 gene is 1868 bp, containing a 1491 bp open reading frame (ORF), and encoding 496 amino acids. ScAMT1.1 belongs to the AMT superfamily and shares 91.57% homologies with AMT1.1 from Oryza sativa. Furthermore, it was stably overexpressed in rice (O. sativa). Under low N treatment, the plant height and the fresh weight of the ScAMT1.1-overexpressed transgenic rice were 36.48% and 51.55% higher than that of the wild-type, respectively. Both the activity of ammonium assimilation key enzymes GS and GDH, and the expression level of ammonium assimilation key genes, including GS1.1, GS1.2, GDH, Fd-GOGAT, and NADH-GOGAT2 in the transgenic plants, were significantly higher than that of the wild-type. The grain number and grain yield per plant in the transgenic rice were 6.44% and 9.52% higher than that of the wild-type in the pot experiments, respectively. Taken together, the sugarcane ScAMT1.1 gene has the potential to improve ammonium assimilation ability and the yield of transgenic rice under low N fertilizer conditions. This study provided an important functional gene for improving sugarcane varieties with high NUE.

MicroRNAs modulating nutrient homeostasis: a sustainable approach for developing biofortified crops

During their lifespan, sessile plants have to cope with bioavailability of the suboptimal nutrient concentration and have to constantly sense/evolve the connecting web of signal cascades for efficient nutrient uptake, storage, and translocation for proper growth and metabolism. However, environmental fluctuations and escalating anthropogenic activities are making it a formidable challenge for plants. This is adding to (micro)nutrient-deficient crops and nutritional insecurity. Biofortification is emerging as a sustainable and efficacious approach which can be utilized to combat the micronutrient malnutrition. A biofortified crop has an enriched level of desired nutrients developed using conventional breeding, agronomic practices, or advanced biotechnological tools. Nutrient homeostasis gets hampered under nutrient stress, which involves disturbance in short-distance and long-distance cell-cell/cell-organ communications involving multiple cellular and molecular components. Advanced sequencing platforms coupled with bioinformatics pipelines and databases have suggested the potential roles of tiny signaling molecules and post-transcriptional regulators, the microRNAs (miRNAs) in key plant phenomena including nutrient homeostasis. miRNAs are seen as emerging targets for biotechnology-based biofortification programs. Thus, understanding the mechanistic insights and regulatory role of miRNAs could open new windows for exploring them in developing nutrient-efficient biofortified crops. This review discusses significance and roles of miRNAs in plant nutrition and nutrient homeostasis and how they play key roles in plant responses to nutrient imbalances/deficiencies/toxicities covering major nutrients-nitrogen (N), phosphorus (P), sulfur (S), magnesium (Mg), iron (Fe), and zinc (Zn). A perspective view has been given on developing miRNA-engineered biofortified crops with recent success stories. Current challenges and future strategies have also been discussed.

Long non-coding RNAs as the regulatory hubs in rice response to salt stress

Salinity seriously constrains growth and fertility of rice worldwide. Long non-coding RNAs (lncRNAs) play crucial roles in plant abiotic stress response. However, salt responsive lncRNAs are poorly understood in rice. Herein, salt responsive lncRNAs (DE-lncRNAs) were identified in FL478 (salt tolerant) compared to its susceptible parent (IR29) using RNA-seq in root tissues at seedling stage. In FL478 and IR29, 8724 and 9235 transcripts with length of > 200 bp were nominated as potential lncRNAs, respectively. Rigorous filtering left four (in FL478) and nine (in IR29) DE-lncRNAs with only 2 DE-lncRNAs in common. ATAC-seq data showed that the genomic regions of all four lncRNAs in FL478 and 6/9 in IR29 are significantly accessible for transcription. Weighted correlation network analysis (WGCNA) revealed that lncRNA.2-FL was highly correlated with 173 mRNAs as trans-targets and a gene encoding pentatricopeptide repeat (PPR) protein was predicted as cis-target of lncRNA.2-FL. In silico mutagenesis analysis proposed the same transcription factor binding sites (TFBSs) in vicinity of the trans- and cis-regulatory target genes of lncRNA.2-FL, which significantly affect their transcription start site (TSS). This study provides new insights into involvement of the DE-lncRNAs in rice response to salt stress. Among them, lncRNA.2-FL may play a significant regulatory role in the salt stress tolerance of FL478.

Identification and Functional Analysis of Long Non-Coding RNA (lncRNA) in Response to Seed Aging in Rice

Many lncRNAs have been shown to play a vital role in aging processes. However, how lncRNAs regulate seed aging remains unknown. In this study, we performed whole transcriptome strand-specific RNA sequencing of samples from rice embryos, analyzed the differences in expression of rice seed lncRNAs before and after artificial aging treatment (AAT), and systematically screened 6002 rice lncRNAs. During the AAT period, the expression levels of most lncRNAs (454) were downregulated and only four were upregulated among the 458 differentially expressed lncRNAs (DELs). Cis- or trans-regulated target genes of the four upregulated lncRNAs were mainly related to base repair, while 454 downregulated lncRNAs were related to plant-pathogen interaction, plant hormones, energy metabolism, and secondary metabolism. The pathways of DEL target genes were similar with those of differentially expressed mRNAs (DEGs). A competing endogenous RNA (ceRNA) network composed of 34 lncRNAs, 24 microRNAs (miRNA), and 161 mRNAs was obtained. The cDNA sequence of lncRNA LNC_037529 was obtained by rapid amplification of cDNA ends (RACE) cloning with a total length of 1325 bp, a conserved 5' end, and a non-conserved 3' end. Together, our findings indicate that genome-wide selection for lncRNA downregulation was an important mechanism for rice seed aging. LncRNAs can be used as markers of seed aging in rice. These findings provide a future path to decipher the underlying mechanism associated with lncRNAs in seed aging.

A long noncoding RNA HILinc1 enhances pear thermotolerance by stabilizing PbHILT1 transcripts through complementary base pairing

As global warming intensifies, heat stress has become a major environmental constraint threatening crop production and quality worldwide. Here, we characterize Heat-induced long intergenic noncoding RNA 1 (HILinc1), a cytoplasm-enriched lincRNA that plays a key role in thermotolerance regulation of pear (Pyrus spp.). HILinc1 Target 1 (PbHILT1) which is the target transcript of HILinc1, was stabilized via complementary base pairing to upregulate its expression. PbHILT1 could bind to Heat shock transcription factor A1b (PbHSFA1b) to enhance its transcriptional activity, leading to the upregulation of a major downstream transcriptional regulator, Multiprotein bridging factor 1c (PbMBF1c), during heat response. Transient overexpressing of either HILinc1 or PbHILT1 increases thermotolerance in pear, while transient silencing of HILinc1 or PbHILT1 makes pear plants more heat sensitive. These findings provide evidences for a new regulatory mechanism by which HILinc1 facilitates PbHSFA1b activity and enhances pear thermotolerance through stabilizing PbHILT1 transcripts.

Heat stress in pear cultivar results in upregulation of long non-coding RNA HILinc1, which binds to and stabilizes PbHILT1 mRNA, which codes for a protein that interacts with heat shock factor A1b, improving thermotolerance.

Identification of long non-coding RNAs involved in floral scent of Rosa hybrida

Long non-coding RNAs (lncRNAs) were found to play important roles in transcriptional, post-transcriptional, and epigenetic gene regulation in various biological processes. However, lncRNAs and their regulatory roles remain poorly studied in horticultural plants. Rose is economically important not only for their wide use as garden and cut flowers but also as important sources of natural fragrance for perfume and cosmetics industry, but presently little was known about the regulatory mechanism of the floral scent production. In this paper, a RNA-Seq analysis with strand-specific libraries, was performed to rose flowers in different flowering stages. The scented variety 'Tianmidemeng' (Rosa hybrida) was used as plant material. A total of 13,957 lncRNAs were identified by mining the RNA-Seq data, including 10,887 annotated lncRNAs and 3070 novel lncRNAs. Among them, 10,075 lncRNAs were predicted to possess a total of 29,622 target genes, including 54 synthase genes and 24 transcription factors related to floral scent synthesis. 425 lncRNAs were differentially expressed during the flowering process, among which 19 were differentially expressed among all the three flowering stages. Using weighted correlation network analysis (WGCNA), we correlate the differentially-expressed lncRNAs to synthesis of individual floral scent compounds. Furthermore, regulatory function of one of candidate lncRNAs for floral scent synthesis was verified using VIGS method in the rose. In this study, we were able to show that lncRNAs may play important roles in floral scent production in the rose. This study also improves our understanding of how plants regulate their secondary metabolism by lncRNAs.

B, Iquebal MA, Bhattacharya RC, Rai A, Kumar D. Genome wide identification of lncRNAs and circRNAs having regulatory role in fruit shelf life in health crop cucumber (Cucumis sativus L.)

Cucumber is an extremely perishable vegetable; however, under room conditions, the fruits become unfit for consumption 2-3 days after harvesting. One natural variant, DC-48 with an extended shelf-life was identified, fruits of which can be stored up to 10-15 days under room temperature. The genes involved in this economically important trait are regulated by non-coding RNAs. The study aims to identify the long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) by taking two contrasting genotypes, DC-48 and DC-83, at two different fruit developmental stages. The upper epidermis of the fruits was collected at 5 days and 10 days after pollination (DAP) for high throughput RNA sequencing. The differential expression analysis was performed to identify differentially expressed (DE) lncRNAs and circRNAs along with the network analysis of lncRNA, miRNA, circRNA, and mRNA interactions. A total of 97 DElncRNAs were identified where 18 were common under both the developmental stages (8 down regulated and 10 upregulated). Based on the back-spliced reads, 238 circRNAs were found to be distributed uniformly throughout the cucumber genomes with the highest numbers (71) in chromosome 4. The majority of the circRNAs (49%) were exonic in origin followed by inter-genic (47%) and intronic (4%) origin. The genes related to fruit firmness, namely, polygalacturonase, expansin, pectate lyase, and xyloglucan glycosyltransferase were present in the target sites and co-localized networks indicating the role of the lncRNA and circRNAs in their regulation. Genes related to fruit ripening, namely, trehalose-6-phosphate synthase, squamosa promoter binding protein, WRKY domain transcription factors, MADS box proteins, abscisic stress ripening inhibitors, and different classes of heat shock proteins (HSPs) were also found to be regulated by the identified lncRNA and circRNAs. Besides, ethylene biosynthesis and chlorophyll metabolisms were also found to be regulated by DElncRNAs and circRNAs. A total of 17 transcripts were also successfully validated through RT PCR data. These results would help the breeders to identify the complex molecular network and regulatory role of the lncRNAs and circRNAs in determining the shelf-life of cucumbers.

Genome-wide identification and characterization of long noncoding RNAs during peach (Prunus persica) fruit development and ripening

LncRNAs represent a class of RNA transcripts of more than 200 nucleotides (nt) in length without discernible protein-coding potential. The expression levels of lncRNAs are significantly affected by stress or developmental cues. Recent studies have shown that lncRNAs participate in fruit development and ripening processes in tomato and strawberry; however, in other fleshy fruits, the association between lncRNAs and fruit ripening remains largely elusive. Here, we constructed 9 ssRNA-Seq libraries from three different peach (Prunus persica) fruit developmental stages comprising the first and second exponential stages and the fruit-ripening stage. In total, 1500 confident lncRNAs from 887 loci were obtained according to the bioinformatics analysis. The lncRNAs identified in peach fruits showed distinct characteristics compared with protein-coding mRNAs, including lower expression levels, lower complexity of alternative splicing, shorter isoforms and smaller numbers of exons. Expression analysis identified 575 differentially expressed lncRNAs (DELs) classified into 6 clusters, among which members of Clusters 1, 2, 4 and 5 were putatively associated with fruit development and ripening processes. Quantitative real-time PCR revealed that the DELs indeed had stage-specific expression patterns in peach fruits. GO and KEGG enrichment analysis revealed that DELs might be associated with fruit-ripening-related physiological and metabolic changes, such as flavonoid biosynthesis, fruit texture softening, chlorophyll breakdown and aroma compound accumulation. Finally, the similarity analysis of lncRNAs within different plant species indicated the low sequence conservation of lncRNAs. Our study reports a large number of fruit-expressed lncRNAs and identifies fruit development phase-specific expressed lncRNA members, which highlights their potential functions in fruit development and ripening processes and lays the foundations for future functional research.

Insights into nitrogen metabolism in the wild and cultivated lettuce as revealed by transcriptome and weighted gene co-expression network analysis

Large amounts of nitrogen fertilizers applied during lettuce (Lactuca sativa L.) production are lost due to leaching or volatilization, causing severe environmental pollution and increased costs of production. Developing lettuce varieties with high nitrogen use efficiency (NUE) is the eco-friendly solution to reduce nitrogen pollution. Hence, in-depth knowledge of nitrogen metabolism and assimilation genes and their regulation is critical for developing high NUE varieties. In this study, we performed comparative transcriptomic analysis of the cultivated lettuce (L. sativa L.) and its wild progenitor (L. serriola) under high and low nitrogen conditions. A total of 2,704 differentially expressed genes were identified. Key enriched biological processes included photosynthesis, oxidation-reduction process, chlorophyll biosynthetic process, and cell redox homeostasis. The transcription factors (TFs) belonging to the ethylene responsive factor family and basic helix-loop-helix family were among the top differentially expressed TFs. Using weighted gene co-expression network analysis we constructed nine co-expression modules. Among these, two modules were further investigated because of their significant association with total nitrogen content and photosynthetic efficiency of photosystem II. Three highly correlated clusters were identified which included hub genes for nitrogen metabolism, secondary metabolites, and carbon assimilation, and were regulated by cluster specific TFs. We found that the expression of nitrogen transportation and assimilation genes varied significantly between the two lettuce species thereby providing the opportunity of introgressing wild alleles into the cultivated germplasm for developing lettuce cultivars with more efficient use of nitrogen.

An Integrated Regulatory Network of mRNAs, microRNAs, and lncRNAs Involved in Nitrogen Metabolism of Moso Bamboo

Nitrogen is a key macronutrient essential for plant growth and development, and its availability has a strong influence on biological processes. Nitrogen fertilizer has been widely applied in bamboo forests in recent decades; however, the mechanism of nitrogen metabolism in bamboo is not fully elucidated. Here, we characterized the morphological, physiological, and transcriptome changes of moso bamboo in response to different schemes for nitrogen addition to illuminate the regulation mechanism of nitrogen metabolism. The appropriate addition of nitrogen improved the chlorophyll content and Pn (net photosynthetic rate) of leaves, the nitrogen and ammonium contents of the seedling roots, the biomass of the whole seedling, the number of lateral roots, and the activity of enzymes involved in nitrogen metabolism in the roots. Based on the whole transcriptome data of the roots, a total of 8,632 differentially expressed mRNAs (DEGs) were identified under different nitrogen additions, such as 52 nitrate transporter genes, 6 nitrate reductase genes, 2 nitrite reductase genes, 2 glutamine synthase genes, 2 glutamate synthase genes (GOGAT), 3 glutamate dehydrogenase genes, and 431 TFs belonging to 23 families. Meanwhile, 123 differentially expressed miRNAs (DEMs) and 396 differentially expressed lncRNAs (DELs) were characterized as nitrogen responsive, respectively. Furthermore, 94 DEM-DEG pairs and 23 DEL-DEG pairs involved in nitrogen metabolism were identified. Finally, a predicted regulatory network of nitrogen metabolism was initially constructed, which included 17 nitrogen metabolic pathway genes, 15 TFs, 4 miRNAs, and 10 lncRNAs by conjoint analysis of DEGs, DEMs, and DELs and their regulatory relationships, which was supported by RNA-seq data and qPCR results. The lncRNA-miRNA-mRNA network provides new insights into the regulation mechanism of nitrogen metabolism in bamboo, which facilitates further genetic improvement for bamboo to adapt to the fluctuating nitrogen environment.

Nitrate-dependent regulation of miR444-OsMADS27 signalling cascade controls root development in rice

Nitrate is an important nutrient and a key signalling molecule for plant development. A number of transcription factors involved in the response to nitrate and their regulatory mechanisms have been identified. However, little is known about the transcription factors involved in nitrate sensing and their regulatory mechanisms among crop plants. In this study, we identified functions of a nitrate-responsive miR444:MADS-box transcription factor OsMADS27 module and its downstream targets mediating rice root growth and stress responses. Transgenic rice plants expressing miR444 target mimic improved rice root growth. Although miR444 has the potential to target multiple genes, we identified OsMADS27 as the major miR444 target that regulates the expression of nitrate transporters, as well as several key genes including expansins, and those associated with auxin signalling, to promote root growth. In agreement with this, overexpression of miRNA-resistant OsMADS27 improved root development and tolerance to abiotic stresses, while its silencing suppressed root growth. OsMADS27 mediated robust stress tolerance in plants through its ability to bind to the promoters of specific stress regulators, as observed in ChIP-seq analysis. Our results provide evidence of a nitrate-dependent miR444-OsMADS27 signalling cascade involved in the regulation of rice root growth, as well as its surprising role in stress responses.

Regulatory role of microRNAs (miRNAs) in the recent development of abiotic stress tolerance of plants

MicroRNAs (miRNAs) are a distinct groups of single-stranded non-coding, tiny regulatory RNAs approximately 20-24 nucleotides in length. miRNAs negatively influence gene expression at the post-transcriptional level and have evolved considerably in the development of abiotic stress tolerance in a number of model plants and economically important crop species. The present review aims to deliver the information on miRNA-mediated regulation of the expression of major genes or Transcription Factors (TFs), as well as genetic and regulatory pathways. Also, the information on adaptive mechanisms involved in plant abiotic stress responses, prediction, and validation of targets, computational tools, and databases available for plant miRNAs, specifically focus on their exploration for engineering abiotic stress tolerance in plants. The regulatory function of miRNAs in plant growth, development, and abiotic stresses consider in this review, which uses high-throughput sequencing (HTS) technologies to generate large-scale libraries of small RNAs (sRNAs) for conventional screening of known and novel abiotic stress-responsive miRNAs adds complexity to regulatory networks in plants. The discoveries of miRNA-mediated tolerance to multiple abiotic stresses, including salinity, drought, cold, heat stress, nutritional deficiency, UV-radiation, oxidative stress, hypoxia, and heavy metal toxicity, are highlighted and discussed in this review.

Recent Advances in Agronomic and Physio-Molecular Approaches for Improving Nitrogen Use Efficiency in Crop Plants

The efficiency with which plants use nutrients to create biomass and/or grain is determined by the interaction of environmental and plant intrinsic factors. The major macronutrients, especially nitrogen (N), limit plant growth and development (1.5-2% of dry biomass) and have a direct impact on global food supply, fertilizer demand, and concern with environmental health. In the present time, the global consumption of N fertilizer is nearly 120 MT (million tons), and the N efficiency ranges from 25 to 50% of applied N. The dynamic range of ideal internal N concentrations is extremely large, necessitating stringent management to ensure that its requirements are met across various categories of developmental and environmental situations. Furthermore, approximately 60 percent of arable land is mineral deficient and/or mineral toxic around the world. The use of chemical fertilizers adds to the cost of production for the farmers and also increases environmental pollution. Therefore, the present study focused on the advancement in fertilizer approaches, comprising the use of biochar, zeolite, and customized nano and bio-fertilizers which had shown to be effective in improving nitrogen use efficiency (NUE) with lower soil degradation. Consequently, adopting precision farming, crop modeling, and the use of remote sensing technologies such as chlorophyll meters, leaf color charts, etc. assist in reducing the application of N fertilizer. This study also discussed the role of crucial plant attributes such as root structure architecture in improving the uptake and transport of N efficiency. The crosstalk of N with other soil nutrients plays a crucial role in nutrient homeostasis, which is also discussed thoroughly in this analysis. At the end, this review highlights the more efficient and accurate molecular strategies and techniques such as N transporters, transgenes, and omics, which are opening up intriguing possibilities for the detailed investigation of the molecular components that contribute to nitrogen utilization efficiency, thus expanding our knowledge of plant nutrition for future global food security.

Nitrate-responsive transcriptome analysis reveals additional genes/processes and associated traits viz. height, tillering, heading date, stomatal density and yield in japonica rice

Our transcriptomic analysis expanded the repertoire of nitrate-responsive genes/processes in rice and revealed their phenotypic association with root/shoot, stomata, tiller, panicle/flowering and yield, with agronomic implications for nitrogen use efficiency. Nitrogen use efficiency (NUE) is a multigenic quantitative trait, involving many N-responsive genes/processes that are yet to be fully characterized. Microarray analysis of early nitrate response in excised leaves of japonica rice revealed 6688 differentially expressed genes (DEGs), including 2640 hitherto unreported across multiple functional categories. They include transporters, enzymes involved in primary/secondary metabolism, transcription factors (TFs), EF-hand containing calcium binding proteins, hormone metabolism/signaling and methytransferases. Some DEGs belonged to hitherto unreported processes viz. alcohol, lipid and trehalose metabolism, mitochondrial membrane organization, protein targeting and stomatal opening. 1158 DEGs were associated with growth physiology and grain yield or phenotypic traits for NUE. We identified seven DEGs for shoot apical meristem, 66 for leaf/culm/root, 31 for tiller, 70 for heading date/inflorescence/spikelet/panicle, 144 for seed and 78 for yield. RT-qPCR validated nitrate regulation of 31 DEGs belonging to various important functional categories/traits. Physiological validation of N-dose responsive changes in plant development revealed that relative to 1.5 mM, 15 mM nitrate significantly increased stomatal density, stomatal conductance and transpiration rate. Further, root/shoot growth, number of tillers and grain yield declined and panicle emergence/heading date delayed, despite increased photosynthetic rate. We report the binding sites of diverse classes of TFs such as WRKY, MYB, HMG etc., in the 1 kb up-stream regions of 6676 nitrate-responsive DEGs indicating their role in regulating nitrate response/NUE. Together, these findings expand the repertoire of genes and processes involved in genomewide nitrate response in rice and reveal their physiological, phenotypic and agronomic implications for NUE.

Genome-Wide Identification and Characterization of Long Noncoding RNAs in Populus × canescens Roots Treated With Different Nitrogen Fertilizers

Nitrate (NO₃ ^-) and ammonium (NH₄ ⁺) are the primary forms of inorganic nitrogen acquired by plant roots. LncRNAs, as key regulators of gene expression, are a class of non-coding RNAs larger than 200 bp. However, knowledge about the regulatory role of lncRNAs in response to different nitrogen forms remains limited, particularly in woody plants. Here, we performed strand-specific RNA-sequencing of P. × canescens roots under three different nitrogen fertilization treatments. In total, 324 lncRNAs and 6,112 mRNAs were identified as showing significantly differential expression between the NO₃ ^- and NH₄NO₃ treatments. Moreover, 333 lncRNAs and 6,007 mRNAs showed significantly differential expression between the NH₄ ⁺ and NH₄NO₃ treatments. Further analysis suggested that these lncRNAs and mRNAs have different response mechanisms for different nitrogen forms. In addition, functional annotation of cis and trans target mRNAs of differentially expressed lncRNAs indicated that 60 lncRNAs corresponding to 49 differentially expressed cis and trans target mRNAs were involved in plant nitrogen metabolism and amino acid biosynthesis and metabolism. Furthermore, 42 lncRNAs were identified as putative precursors of 63 miRNAs, and 28 differentially expressed lncRNAs were potential endogenous target mimics targeted by 96 miRNAs. Moreover, ceRNA regulation networks were constructed. MSTRG.6097.1, MSTRG.13550.1, MSTRG.2693.1, and MSTRG.12899.1, as hub lncRNAs in the ceRNA networks, are potential candidate lncRNAs for studying the regulatory mechanism in poplar roots under different nitrogen fertilization treatments. The results provide a basis for obtaining insight into the molecular mechanisms of lncRNA responses to different nitrogen forms in woody plants.

miRNAomic Approach to Plant Nitrogen Starvation

Nitrogen (N) is one of the indispensable nutrients required by plants for their growth, development, and survival. Being a limited nutrient, it is mostly supplied exogenously to the plants, to maintain quality and productivity. The increased use of N fertilizers is associated with high-cost inputs and negative environmental consequences, which necessitates the development of nitrogen-use-efficient plants for sustainable agriculture. Understanding the regulatory mechanisms underlying N metabolism in plants under low N is one of the prerequisites for the development of nitrogen-use-efficient plants. One of the important and recently discovered groups of regulatory molecules acting at the posttranscriptional and translational levels are microRNAs (miRNAs). miRNAs are known to play critical roles in the regulation of gene expression in plants under different stress conditions including N stress. Several classes of miRNAs associated with N metabolism have been identified so far. These nitrogen-responsive miRNAs may provide a platform for a better understanding of the regulation of N metabolism and pave a way for the development of genotypes for better N utilization. The current review presents a brief outline of miRNAs and their regulatory role in N metabolism.

Transcriptome Analysis Reveals Different Responsive Patterns to Nitrogen Deficiency in Two Wheat Near-Isogenic Lines Contrasting for Nitrogen Use Efficiency

Simple Summary

Nitrogen (N) limitation is the key factor for wheat production worldwide. Therefore, the development of genotypes with improved nitrogen use efficiency (NUE) is a prerequisite for sustainable and productive agriculture. Exploring the molecular mechanisms of low N stress tolerance is significant for breeding wheat cultivars with high NUE. To clarify the underlying molecular mechanisms of enhanced resilience to low N in high-NUE wheat, we performed an RNA sequencing (RNA-seq) analysis. In the current research, two wheat near-isogenic lines (NILs) differing dramatically in NUE were used to measure gene expression differences under different N treatments. There was a dramatic difference between two wheat NILs in response to N deficiency at the transcriptional level, and the classification of identified candidate genes may provide new valuable insights into the resilience mechanism of wheat.

Abstract

The development of crop cultivars with high nitrogen use efficiency (NUE) under low-N fertilizer inputs is imperative for sustainable agriculture. However, there has been little research on the molecular mechanisms underlying enhanced resilience to low N in high-NUE plants. The comparison of the transcriptional responses of genotypes contrasting for NUE will facilitate an understanding of the key molecular mechanism of wheat resilience to low-N stress. In the current study, the RNA sequencing (RNA-seq) technique was employed to investigate the genotypic difference in response to N deficiency between two wheat NILs (1Y, high-NUE, and 1W, low-NUE). In our research, high- and low-NUE wheat NILs showed different patterns of gene expression under N-deficient conditions, and these N-responsive genes were classified into two major classes, including “frontloaded genes” and “relatively upregulated genes”. In total, 103 and 45 genes were identified as frontloaded genes in high-NUE and low-NUE wheat, respectively. In summary, our study might provide potential directions for further understanding the molecular mechanism of high-NUE genotypes adapting to low-N stress.

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.

Genome-wide in silico analysis of long intergenic non-coding RNAs from rice peduncles at the heading stage

Long intergenic non-coding RNAs (lincRNAs) belong to the category of long non-coding RNAs (lncRNAs), originated from intergenic regions, which do not code for proteins. LincRNAs perform prominent role in regulation of gene expression during plant development and stress response by directly interacting with DNA, RNA, or proteins, or triggering production of small RNA regulatory molecules. Here, we identified 2973 lincRNAs and investigated their expression dynamics during peduncle elongation in two Indian rice cultivars, Pokkali and Swarna, at the time of heading. Differential expression analysis revealed common and cultivar-specific expression patterns, which we utilized to infer the lincRNA candidates with potential involvement in peduncle elongation and panicle exsertion. Their putative targets were identified using in silico prediction methods followed by pathway mapping and literature-survey based functional analysis. Further, to infer the mechanism of action, we identified the lincRNAs which potentially act as miRNA precursors or target mimics.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01059-2.

LncRNA TCONS_00021861 is functionally associated with drought tolerance in rice (Oryza sativa L.) via competing endogenous RNA regulation

BACKGROUND

Water deficit is an abiotic stress that retards plant growth and destabilizes crop production. Long non coding RNAs (lncRNAs) are a class of non-coding endogenous RNAs that participate in diverse cellular processes and stress responses in plants. lncRNAs could function as competing endogenous RNAs (ceRNA) and represent a novel layer of gene regulation. However, the regulatory mechanism of lncRNAs as ceRNA in drought stress response is yet unclear.

RESULTS

In this study, we performed transcriptome-wide identification of drought-responsive lncRNAs in rice. Thereafter, we constructed a lncRNA-mediated ceRNA network by analyzing competing relationships between mRNAs and lncRNAs based on ceRNA hypothesis. A drought responsive ceRNA network with 40 lncRNAs, 23 miRNAs and 103 mRNAs was obtained. Network analysis revealed TCONS_00021861/miR528-3p/YUCCA7 regulatory axis as a hub involved in drought response. The miRNA-target expression and interaction were validated by RT-qPCR and RLM-5'RACE. TCONS_00021861 showed significant positive correlation (r = 0.7102) with YUCCA7 and negative correlation with miR528-3p (r = -0.7483). Overexpression of TCONS_00021861 attenuated the repression of miR528-3p on YUCCA7, leading to increased IAA (Indole-3-acetic acid) content and auxin overproduction phenotypes.

CONCLUSIONS

TCONS_00021861 could regulate YUCCA7 by sponging miR528-3p, which in turn activates IAA biosynthetic pathway and confer resistance to drought stress. Our findings provide a new perspective of the regulatory roles of lncRNAs as ceRNAs in drought resistance of rice.

The N uptake-associated physiological processes at late growth stage in wheat (Triticum aestivum) under N deprivation combined with deficit irrigation condition

Elucidating physiological mechanisms underlying the plant N uptake benefits breeding of high N use efficiency (NUE) crop cultivars. In this study, we investigated the growth and N uptake-associated processes in wheat under N deprivation and deficit irrigation, using two contrasting NUE cultivars. Compared with sufficient-N (SN), deficient-N (DN) treatment reduced plant biomass, N accumulation, and yields in two cultivars (high NUE Shinong 086 and N deprivation-sensitive Jimai 585), suggesting that N deprivation negatively regulates plant growth and N uptake. Shinong 086 was better on growth and N uptake-associated traits than Jimai 585 due to the improved root biomass across soil profile, which was consistent with the decrease of available N contents in soil layers. These results suggested that the improved root system architecture (RAS) enhances plant acquirement for soil N under N- and water-deprivation condition, contributing to the plant N uptake and yield formation capacities. Transcriptome investigation revealed that numerous genes were differentially expressed (DE) in the N-deprived Shinong 086 plants, which involve the regulation of complicate biochemical pathways. These results suggested that the modified RAS and N uptake in high NUE plants are accomplished underlying the regulation of numerous DE genes. TaWRKY20, a gene in ZFP transcription factor family, was functionally characterized for the role in mediating plant N uptake. Overexpression of it conferred plants improved growth and N uptake under DN due to its regulation on TaNRT2.1 and TaNRT2.2, two nitrate transporter genes. Our investigation provides insights in high NUE mechanisms in wheat under N deprivation.

Low nitrogen availability inhibits the phosphorus starvation response in maize (Zea mays ssp. mays L.)

BACKGROUND

Nitrogen (N) and phosphorus (P) are macronutrients essential for crop growth and productivity. In cultivated fields, N and P levels are rarely sufficient, contributing to the gap between realized and potential production. Fertilizer application increases nutrient availability, but is not available to all farmers, nor are current rates of application sustainable or environmentally desirable. Transcriptomic studies of cereal crops have revealed dramatic responses to either low N or low P single stress treatments. In the field, however, levels of both N and P may be suboptimal. The interaction between N and P starvation responses remains to be fully characterized.

RESULTS

We characterized growth and root and leaf transcriptomes of young maize plants under nutrient replete, low N, low P or combined low NP conditions. We identified 1555 genes to respond to our nutrient treatments, in one or both tissues. A large group of genes, including many classical P starvation response genes, were regulated antagonistically between low N and P conditions. An additional experiment over a range of N availability indicated that a mild reduction in N levels was sufficient to repress the low P induction of P starvation genes. Although expression of P transporter genes was repressed under low N or low NP, we confirmed earlier reports of P hyper accumulation under N limitation.

CONCLUSIONS

Transcriptional responses to low N or P were distinct, with few genes responding in a similar way to the two single stress treatments. In combined NP stress, the low N response dominated, and the P starvation response was largely suppressed. A mild reduction in N availability was sufficient to repress the induction of P starvation associated genes. We conclude that activation of the transcriptional response to P starvation in maize is contingent on N availability.

MicroRNAs: Potential Targets for Developing Stress-Tolerant Crops

Crop yield is challenged every year worldwide by changing climatic conditions. The forecasted climatic scenario urgently demands stress-tolerant crop varieties to feed the ever-increasing global population. Molecular breeding and genetic engineering approaches have been frequently exploited for developing crops with desired agronomic traits. Recently, microRNAs (miRNAs) have emerged as powerful molecules, which potentially serve as expression markers during stress conditions. The miRNAs are small non-coding endogenous RNAs, usually 20-24 nucleotides long, which mediate post-transcriptional gene silencing and fine-tune the regulation of many abiotic- and biotic-stress responsive genes in plants. The miRNAs usually function by specifically pairing with the target mRNAs, inducing their cleavage or repressing their translation. This review focuses on the exploration of the functional role of miRNAs in regulating plant responses to abiotic and biotic stresses. Moreover, a methodology is also discussed to mine stress-responsive miRNAs from the enormous amount of transcriptome data available in the public domain generated using next-generation sequencing (NGS). Considering the functional role of miRNAs in mediating stress responses, these molecules may be explored as novel targets for engineering stress-tolerant crop varieties.

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.

Down regulation of transcripts involved in selective metabolic pathways as an acclimation strategy in nitrogen use efficient genotypes of rice under low nitrogen

To understand the molecular mechanism of nitrogen use efficiency (NUE) in rice, two nitrogen (N) use efficient genotypes and two non-efficient genotypes were characterized using transcriptome analyses. The four genotypes were evaluated for 3 years under low and recommended N field conditions for 12 traits/parameters of yield, straw, nitrogen content along with NUE indices and 2 promising donors for rice NUE were identified. Using the transcriptome data generated from GS FLX 454 Roche and Illumina HiSeq 2000 of two efficient and two non-efficient genotypes grown under field conditions of low N and recommended N and their de novo assembly, differentially expressed transcripts and pathways during the panicle development were identified. Down regulation was observed in 30% of metabolic pathways in efficient genotypes and is being proposed as an acclimation strategy to low N. Ten sub metabolic pathways significantly enriched with additional transcripts either in the direction of the common expression or contra-regulated to the common expression were found to be critical for NUE in rice. Among the up-regulated transcripts in efficient genotypes, a hypothetical protein OsI_17904 with 2 alternative forms suggested the role of alternative splicing in NUE of rice and a potassium channel SKOR transcript (LOC_Os06g14030) has shown a positive correlation (0.62) with single plant yield under low N in a set of 16 rice genotypes. From the present study, we propose that the efficient genotypes appear to down regulate several not so critical metabolic pathways and divert the thus conserved energy to produce seed/yield under long-term N starvation.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-020-02631-5.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Heterotrimeric G-protein α subunit (RGA1) regulates tiller development, yield, cell wall, nitrogen response and biotic stress in rice

G-proteins are implicated in plant productivity, but their genome-wide roles in regulating agronomically important traits remain uncharacterized. Transcriptomic analyses of rice G-protein alpha subunit mutant (rga1) revealed 2270 differentially expressed genes (DEGs) including those involved in C/N and lipid metabolism, cell wall, hormones and stress. Many DEGs were associated with root, leaf, culm, inflorescence, panicle, grain yield and heading date. The mutant performed better in total weight of filled grains, ratio of filled to unfilled grains and tillers per plant. Protein–protein interaction (PPI) network analysis using experimentally validated interactors revealed many RGA1-responsive genes involved in tiller development. qPCR validated the differential expression of genes involved in strigolactone-mediated tiller formation and grain development. Further, the mutant growth and biomass were unaffected by submergence indicating its role in submergence response. Transcription factor network analysis revealed the importance of RGA1 in nitrogen signaling with DEGs such as Nin-like, WRKY, NAC, bHLH families, nitrite reductase, glutamine synthetase, OsCIPK23 and urea transporter. Sub-clustering of DEGs-associated PPI network revealed that RGA1 regulates metabolism, stress and gene regulation among others. Predicted rice G-protein networks mapped DEGs and revealed potential effectors. Thus, this study expands the roles of RGA1 to agronomically important traits and reveals their underlying processes.

Integrated analysis of lncRNA and mRNA transcriptomes reveals the potential regulatory role of lncRNA in kiwifruit ripening and softening

Kiwifruit has gained increasing attention worldwide for its unique flavor and high nutritional value. Rapid softening after harvest greatly shortens its shelf-life and reduces the commercial value. Therefore, it is imperative and urgent to identify and clarify its softening mechanism. This study aimed to analyze and compare the long noncoding RNA (lncRNA) and mRNA expression patterns in ABA-treated (ABA) and room temperature (RT)-stored fruits with those in freshly harvested fruits (CK) as control. A total of 697 differentially expressed genes (DEGs) and 81 differentially expressed lncRNAs (DELs) were identified while comparing ABA with CK, and 458 DEGs and 143 DELs were detected while comparing RT with CK. The Kyoto Encyclopedia of Genes and Genomes analysis of the identified DEGs and the target genes of DELs revealed that genes involved in starch and sucrose metabolism, brassinosteroid biosynthesis, plant hormone signal transduction, and flavonoid biosynthesis accounted for a large part. The co-localization networks, including 38 DEGs and 31 DELs in ABA vs. CK, and 25 DEGs and 25 DELs in RT vs. CK, were also performed. Genes related to fruit ripening, such as genes encoding β-galactosidase, mannan endo-1,4-β-mannosidase, pectinesterase/pectinesterase inhibitor, and NAC transcription factor, were present in the co-localization network, suggesting that lncRNAs were involved in regulating kiwifruit ripening. Notably, several ethylene biosynthesis- and signaling-related genes, including one 1-aminocyclopropane-1-carboxylic acid oxidase gene and three ethylene response factor genes, were found in the co-localization network of ABA vs. CK, suggesting that the promoting effect of ABA on ethylene biosynthesis and fruit softening might be embodied by increasing the expression of these lncRNAs. These results may help understand the regulatory mechanism of lncRNAs in ripening and ABA-induced fruit softening of kiwifruit.

Noncoding RNA: An Insight into Chloroplast and Mitochondrial Gene Expressions

Regulation of gene expression in any biological system is a complex process with many checkpoints at the transcriptional, post-transcriptional and translational levels. The control mechanism is mediated by various protein factors, secondary metabolites and a newly included regulatory member, i.e., noncoding RNAs (ncRNAs). It is known that ncRNAs modulate the mRNA or protein profiles of the cell depending on the degree of complementary and context of the microenvironment. In plants, ncRNAs are essential for growth and development in normal conditions by controlling various gene expressions and have emerged as a key player to guard plants during adverse conditions. In order to have smooth functioning of the plants under any environmental pressure, two very important DNA-harboring semi-autonomous organelles, namely, chloroplasts and mitochondria, are considered as main players. These organelles conduct the most crucial metabolic pathways that are required to maintain cell homeostasis. Thus, it is imperative to explore and envisage the molecular machineries responsible for gene regulation within the organelles and their coordination with nuclear transcripts. Therefore, the present review mainly focuses on ncRNAs origination and their gene regulation in chloroplasts and plant mitochondria.