Metabolic and co-expression network-based analyses associated with nitrate response in rice

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.

Allantoin: A Potential Compound for the Mitigation of Adverse Effects of Abiotic Stresses in Plants

Stress-induced alterations vary with the species of plants, the intensity and duration of the exposure, and stressors availability in nature or soil. Purine catabolism acts as an inherent defensive mechanism against various abiotic stresses and plays a pivotal role in the stress acclimatisation of plants. The intermediate metabolite of purine catabolism, allantoin, compensates for soil nitrogen deficiency due to the low carbon/nitrogen ratio, thereby maintaining nitrogen homeostasis and supporting plant growth and development. Allantoin accounts for 90% of the total nitrogenous compound in legumes, while it contributes only 15% in non-leguminous plants. Moreover, studies on a variety of plant species have reported the differential accumulation of allantoin in response to abiotic stresses, endowing allantoin as a stress modulator. Allantoin functions as signalling molecule to stimulate stress-responsive genes (P5CS; pyrroline-5-carboxylase synthase) and ROS (reactive oxygen species) scavenging enzymes (antioxidant). Moreover, it regulates cross-talk between the abscisic acid and jasmonic acid pathway, and maintains ion homeostasis by increasing the accumulation of putrescine and/or spermine, consequently enhancing the tolerance against stress conditions. Further, key enzymes of purine catabolism (xanthine dehydrogenase and allantoinase) have also been explored by constructing various knockdown/knockout mutant lines to decipher their impact on ROS-mediated oxidative injury in plants. Thus, it is established that allantoin serves as a regulatory signalling metabolite in stress protection, and therefore a lower accumulation of allantoin also reduces plant stress tolerance mechanisms. This review gives an account of metabolic regulation and the possible contribution of allantoin as a photo protectant, osmoprotectant, and nitrogen recycler to reduce abiotic-stress-induced impacts on plants.

Weighted gene co-expression network analysis of nitrogen (N)-responsive genes and the putative role of G-quadruplexes in N use efficiency (NUE) in rice

Rice is an important target to improve crop nitrogen (N) use efficiency (NUE), and the identification and shortlisting of the candidate genes are still in progress. We analyzed data from 16 published N-responsive transcriptomes/microarrays to identify, eight datasets that contained the maximum number of 3020 common genes, referred to as N-responsive genes. These include different classes of transcription factors, transporters, miRNA targets, kinases and events of post-translational modifications. A Weighted gene co-expression network analysis (WGCNA) with all the 3020 N-responsive genes revealed 15 co-expression modules and their annotated biological roles. Protein-protein interaction network analysis of the main module revealed the hub genes and their functional annotation revealed their involvement in the ubiquitin process. Further, the occurrences of G-quadruplex sequences were examined, which are known to play important roles in epigenetic regulation but are hitherto unknown in N-response/NUE. Out of the 3020 N-responsive genes studied, 2298 contained G-quadruplex sequences. We compared these N-responsive genes containing G-quadruplex sequences with the 3601 genes we previously identified as NUE-related (for being both N-responsive and yield-associated). This analysis revealed 389 (17%) NUE-related genes containing G-quadruplex sequences. These genes may be involved in the epigenetic regulation of NUE, while the rest of the 83% (1811) genes may regulate NUE through genetic mechanisms and/or other epigenetic means besides G-quadruplexes. A few potentially important genes/processes identified as associated with NUE were experimentally validated in a pair of rice genotypes contrasting for NUE. The results from the WGCNA and G4 sequence analysis of N-responsive genes helped identify and shortlist six genes as candidates to improve NUE. Further, the hitherto unavailable segregation of genetic and epigenetic gene targets could aid in informed interventions through genetic and epigenetic means of crop improvement.

Transcriptome and Co-Expression Network Analysis Reveals the Molecular Mechanism of Rice Root Systems in Response to Low-Nitrogen Conditions

Nitrogen is an important nutrient for plant growth and essential metabolic processes. Roots integrally obtain nutrients from soil and are closely related to the growth and development of plants. In this study, the morphological analysis of rice root tissues collected at different time points under low-nitrogen and normal nitrogen conditions demonstrated that, compared with normal nitrogen treatment, the root growth and nitrogen use efficiency (NUE) of rice under low-nitrogen treatment were significantly improved. To better understand the molecular mechanisms of the rice root system’s response to low-nitrogen conditions, a comprehensive transcriptome analysis of rice seedling roots under low-nitrogen and control conditions was conducted in this study. As a result, 3171 differentially expressed genes (DEGs) were identified. Rice seedling roots enhance NUE and promote root development by regulating the genes related to nitrogen absorption and utilization, carbon metabolism, root growth and development, and phytohormones, thereby adapting to low-nitrogen conditions. A total of 25,377 genes were divided into 14 modules using weighted gene co-expression network analysis (WGCNA). Two modules were significantly associated with nitrogen absorption and utilization. A total of 8 core genes and 43 co-expression candidates related to nitrogen absorption and utilization were obtained in these two modules. Further studies on these genes will contribute to the understanding of low-nitrogen adaptation and nitrogen utilization mechanisms in rice.

Co-expression network of heat-response transcripts: A glimpse into how splicing factors impact rice basal thermotolerance

To identify novel solutions to improve rice yield under rising temperatures, molecular components of thermotolerance must be better understood. Alternative splicing (AS) is a major post-transcriptional mechanism impacting plant tolerance against stresses, including heat stress (HS). AS is largely regulated by splicing factors (SFs) and recent studies have shown their involvement in temperature response. However, little is known about the splicing networks between SFs and AS transcripts in the HS response. To expand this knowledge, we constructed a co-expression network based on a publicly available RNA-seq dataset that explored rice basal thermotolerance over a time-course. Our analyses suggest that the HS-dependent control of the abundance of specific transcripts coding for SFs might explain the widespread, coordinated, complex, and delicate AS regulation of critical genes during a plant's inherent response to extreme temperatures. AS changes in these critical genes might affect many aspects of plant biology, from organellar functions to cell death, providing relevant regulatory candidates for future functional studies of basal thermotolerance.

Co-expression pan-network reveals genes involved in complex traits within maize pan-genome

BACKGROUND

With the advances in the high throughput next generation sequencing technologies, genome-wide association studies (GWAS) have identified a large set of variants associated with complex phenotypic traits at a very fine scale. Despite the progress in GWAS, identification of genotype-phenotype relationship remains challenging in maize due to its nature with dozens of variants controlling the same trait. As the causal variations results in the change in expression, gene expression analyses carry a pivotal role in unraveling the transcriptional regulatory mechanisms behind the phenotypes.

RESULTS

To address these challenges, we incorporated the gene expression and GWAS-driven traits to extend the knowledge of genotype-phenotype relationships and transcriptional regulatory mechanisms behind the phenotypes. We constructed a large collection of gene co-expression networks and identified more than 2 million co-expressing gene pairs in the GWAS-driven pan-network which contains all the gene-pairs in individual genomes of the nested association mapping (NAM) population. We defined four sub-categories for the pan-network: (1) core-network contains the highest represented ~ 1% of the gene-pairs, (2) near-core network contains the next highest represented 1-5% of the gene-pairs, (3) private-network contains ~ 50% of the gene pairs that are unique to individual genomes, and (4) the dispensable-network contains the remaining 50-95% of the gene-pairs in the maize pan-genome. Strikingly, the private-network contained almost all the genes in the pan-network but lacked half of the interactions. We performed gene ontology (GO) enrichment analysis for the pan-, core-, and private- networks and compared the contributions of variants overlapping with genes and promoters to the GWAS-driven pan-network.

CONCLUSIONS

Gene co-expression networks revealed meaningful information about groups of co-regulated genes that play a central role in regulatory processes. Pan-network approach enabled us to visualize the global view of the gene regulatory network for the studied system that could not be well inferred by the core-network alone.

Altering ureide transport in nodulated soybean results in whole-plant adjustments of metabolism, assimilate partitioning, and sink strength

Legumes develop a symbiotic relationship with bacteria that are housed in root nodules and fix atmospheric di-nitrogen (N2) to ammonia. In soybean (Glycine max (L.) Merr.) nodules, the final products of nitrogen (N) fixation are amino acids, and the ureides allantoin and allantoic acid that also serve as the major long-distance N transport forms. Recently, we have shown that increased expression of UPS1 (ureide permease 1) in soybean nodules results in enhanced ureide export from nodules with positive effects on N fixation and seed yield. Here, we demonstrate that changes in the ureide transport processes trigger alterations in allantoin and allantoic acid pools and partitioning throughout the transgenic plants. They further result in adjustments in amino acid availability in, and translocation to, root and shoot sinks. In addition, leaf carbon (C) capture, assimilation and allocation to sinks are improved, accommodating the increased nodule function, and root and shoot growth. Overall, we demonstrate that enhanced ureide partitioning in nodulated soybean leads to a complex rebalancing of N and C acquisition, metabolism, and transport processes with positive consequences for above- and below-ground vegetative biomass, and whole-plant N and C gains.

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.

Delineation of mechanistic approaches employed by plant growth promoting microorganisms for improving drought stress tolerance in plants

Drought stress is expected to increase in intensity, frequency, and duration in many parts of the world, with potential negative impacts on plant growth and productivity. The plants have evolved complex physiological and biochemical mechanisms to respond and adjust to water-deficient environments. The physiological and biochemical mechanisms associated with water-stress tolerance and water-use efficiency have been extensively studied. Besides these adaptive and mitigating strategies, the plant growth-promoting rhizobacteria (PGPR) play a significant role in alleviating plant drought stress. These beneficial microorganisms colonize the endo-rhizosphere/rhizosphere of plants and enhance drought tolerance. The common mechanism by which these microorganisms improve drought tolerance included the production of volatile compounds, phytohormones, siderophores, exopolysaccharides, 1-aminocyclopropane-1-carboxylate deaminase (ACC deaminase), accumulation of antioxidant, stress-induced metabolites such as osmotic solutes proline, alternation in leaf and root morphology and regulation of the stress-responsive genes. The PGPR is an easy and efficient alternative approach to genetic manipulation and crop enhancement practices because plant breeding and genetic modification are time-consuming and expensive processes for obtaining stress-tolerant varieties. In this review, we will elaborate on PGPR's mechanistic approaches in enhancing the plant stress tolerance to cope with the drought stress.

Sugar and Nitrate Sensing: A Multi-Billion-Year Story

Sugars and nitrate play a major role in providing carbon and nitrogen in plants. Understanding how plants sense these nutrients is crucial, most notably for crop improvement. The mechanisms underlying sugar and nitrate sensing are complex and involve moonlighting proteins such as the nitrate transporter NRT1.1/NFP6.3 or the glycolytic enzyme HXK1. Major components of nutrient signaling, such as SnRK1, TOR, and HXK1, are relatively well conserved across eukaryotes, and the diversification of components such as the NRT1 family and the SWEET sugar transporters correlates with plant terrestrialization. In plants, Tre6P plays a hormone-like role in plant development. In addition, nutrient signaling has evolved to interact with the more recent hormone signaling, allowing fine-tuning of physiological and developmental responses.

Identification of gene co-expression networks and key genes regulating flavonoid accumulation in apple (Malus × domestica) fruit skin

Anthocyanin provides a red color for apple and health benefit for human. To better understand the molecular mechanisms of regulating apple color formation, we analyzed 27 transcriptomes of fruit skin from three cultivars 'Huashuo' (red-skinned), 'Hongcuibao' (red-skinned), and 'Golden Delicious' (yellow-skinned) at 0, 2, and 6 days after bag removal. Using pairwise comparisons and weighted gene co-expression network analyses (WGCNA), we constructed 17 co-expression modules. Among them, a specific module was negatively correlated to anthocyanin accumulation. The genes in the module are enriched in flavonoid biosynthesis pathways. These pathway genes were used to construct gene co-expression network of anthocyanin accumulation. Finally, a R2R3-MYB repressor designated MdMYB28 was identified as a key hub gene in the anthocyanin metabolism network. During the anthocyanin accumulation of apple fruit skin reaching a peak, MdMYB28 expression level was negatively correlated with the anthocyanin content. MdMYB28 was shown to directly bind to the promoter of MdMYB10 in yeast one-hybrid analyses. Over-expression of MdMYB28 decreased the anthocyanin biosynthesis in tobacco flower petals, suggesting that MdMYB28 acts as a negatively regulator of anthocyanin biosynthesis.

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.

Meta-Analysis of Yield-Related and N-Responsive Genes Reveals Chromosomal Hotspots, Key Processes and Candidate Genes for Nitrogen-Use Efficiency in Rice

Nitrogen-use efficiency (NUE) is a function of N-response and yield that is controlled by many genes and phenotypic parameters that are poorly characterized. This study compiled all known yield-related genes in rice and mined them from the N-responsive microarray data to find 1,064 NUE-related genes. Many of them are novel genes hitherto unreported as related to NUE, including 80 transporters, 235 transcription factors (TFs), 44 MicroRNAs (miRNAs), 91 kinases, and 8 phosphatases. They were further shortlisted to 62 NUE-candidate genes following hierarchical methods, including quantitative trait locus (QTL) co-localization, functional evaluation in the literature, and protein-protein interactions (PPIs). They were localized to chromosomes 1, 3, 5, and 9, of which chromosome 1 with 26 genes emerged as a hotspot for NUE spanning 81% of the chromosomes. Further, co-localization of the NUE genes on NUE-QTLs resolved differences in the earlier studies that relied mainly on N-responsive genes regardless of their role in yield. Functional annotations and PPIs for all the 1,064 NUE-related genes and also the shortlisted 62 candidates revealed transcription, redox, phosphorylation, transport, development, metabolism, photosynthesis, water deprivation, and hormonal and stomatal function among the prominent processes. In silico expression analysis confirmed differential expression of the 62 NUE-candidate genes in a tissue/stage-specific manner. Experimental validation in two contrasting genotypes revealed that high NUE rice shows better photosynthetic performance, transpiration efficiency and internal water-use efficiency in comparison to low NUE rice. Feature Selection Analysis independently identified one-third of the common genes at every stage of hierarchical shortlisting, offering 6 priority targets to validate for improving the crop NUE.

Genome- wide structural and functional variant discovery of rice landraces using genotyping by sequencing

Rice landraces are vital genetic resources for agronomic and quality traits but the undeniable collection of Kerala landraces remains poorly delineated. To effectively conserve, manage, and use these resources, understanding the genomic structure of germplasm is essential. Genotyping by sequencing (GBS) enables identification of an immense number of single nucleotide polymorphism (SNP) and insertion deletion (InDel) from 96 rice germplasm. In the present study, a total of 16.9 × 10⁷ reads were generated, and among that 16.3 × 10⁷ reads were mapped to the indica reference genome. Exploring GBS data unfolded a wide genomic variations including 82,59,639 SNPs and 1,07,140 Indels. Both neighbor-joining tree and principal coordinate analysis with InDel markers revealed the selected germplasm in this study as highly diverse in structure. We assembled unmapped reads which were further employed for gene ontology analysis. These unmapped sequences that are generally expelled from subsequent studies of GBS data analysis may exist as an unexplored resort for several novel significant biological findings. The discovery of SNPs from the haplotyping results of GS3 and GIF1 genes provided insight into marker- assisted selection based on grain size and yield and can be utilized for rice yield improvement. To our knowledge, this is the first report on structural variation analysis using the GBS platform in rice landraces collected from Kerala. Genomic information from this study endows with valuable resources for perceptive rice landrace structure and can also facilitate sequencing-based molecular breeding.

Gene regulatory network and its constituent transcription factors that control nitrogen-deficiency responses in rice

Increase in the nitrogen (N)-use efficiency and optimization of N response in crop species are urgently needed. Although transcription factor-based genetic engineering is a promising approach for achieving these goals, transcription factors that play key roles in the response to N deficiency have not been studied extensively. Here, we performed RNA-seq analysis of root samples of 20 Asian rice (Oryza sativa) accessions with differential nutrient uptake. Data obtained from plants exposed to N-replete and N-deficient conditions were subjected to coexpression analysis and machine learning-based pathway inference to dissect the gene regulatory network required for the response to N deficiency. Four transcription factors, including members of the G2-like and bZIP families, were predicted to function as key regulators of gene transcription within the network in response to N deficiency. Cotransfection assays validated inferred novel regulatory pathways, and further analyses using genome-edited knockout lines suggested that these transcription factors are important for N-deficiency responses in planta. Many of the N deficiency-responsive genes, including those encoding key regulators within the network, were coordinately regulated by transcription factors belonging to different families. Transcription factors identified in this study could be valuable for the modification of N response and metabolism.

Role of ureides in source-to-sink transport of photoassimilates in non-fixing soybean

Nitrogen (N)-fixing soybean plants use the ureides allantoin and allantoic acid as major long-distance transport forms of N, but in non-fixing, non-nodulated plants amino acids mainly serve in source-to-sink N allocation. However, some ureides are still synthesized in roots of non-fixing soybean, and our study addresses the role of ureide transport processes in those plants. In previous work, legume ureide permeases (UPSs) were identified that are involved in cellular import of allantoin and allantoic acid. Here, UPS1 from common bean was expressed in the soybean phloem, which resulted in enhanced source-to-sink transport of ureides in the transgenic plants. This was accompanied by increased ureide synthesis and elevated allantoin and allantoic acid root-to-sink transport. Interestingly, amino acid assimilation, xylem transport, and phloem partitioning to sinks were also strongly up-regulated. In addition, photosynthesis and sucrose phloem transport were improved in the transgenic plants. These combined changes in source physiology and assimilate partitioning resulted in increased vegetative growth and improved seed numbers. Overall, the results support that ureide transport processes in non-fixing plants affect source N and carbon acquisition and assimilation as well as source-to-sink translocation of N and carbon assimilates with consequences for plant growth and seed development.

Transcriptomic and network analyses reveal distinct nitrate responses in light and dark in rice leaves (Oryza sativa Indica var. Panvel1)

Nitrate (N) response is modulated by light, but not understood from a genome-wide perspective. Comparative transcriptomic analyses of nitrate response in light-grown and etiolated rice leaves revealed 303 and 249 differentially expressed genes (DEGs) respectively. A majority of them were exclusive to light (270) or dark (216) condition, whereas 33 DEGs were common. The latter may constitute response to N signaling regardless of light. Functional annotation and pathway enrichment analyses of the DEGs showed that nitrate primarily modulates conserved N signaling and metabolism in light, whereas oxidation–reduction processes, pentose-phosphate shunt, starch-, sucrose- and glycerolipid-metabolisms in the dark. Differential N-regulation of these pathways by light could be attributed to the involvement of distinctive sets of transporters, transcription factors, enriched cis-acting motifs in the promoters of DEGs as well as differential modulation of N-responsive transcriptional regulatory networks in light and dark. Sub-clustering of DEGs-associated protein–protein interaction network constructed using experimentally validated interactors revealed that nitrate regulates a molecular complex consisting of nitrite reductase, ferredoxin-NADP reductase and ferredoxin. This complex is associated with flowering time, revealing a meeting point for N-regulation of N-response and N-use efficiency. Together, our results provide novel insights into distinct pathways of N-signaling in light and dark conditions.

Metabolic signatures provide novel insights to Picrorhiza kurroa adaptation along the altitude in Himalayan region

INTRODUCTION

Along the altitude, environmental conditions vary significantly that might influence plant performance and distribution. Adaptation to these changing conditions is a complex biological process that involves reprogramming of genes, proteins and metabolites. The metabolic response of medicinal plants along the altitude has been less explored yet.

OBJECTIVES

In the present study, we investigated the adaptation strategies of Picrorhiza kurroa Royle ex Benth. along the altitude in organ specific manner using metabolomic approach.

METHODS

Picrorhiza kurroa plants at flowering stage were randomly sampled from three altitudes viz. 3400, 3800 and 4100 masl in the Himalayan region. Leaf, root and rhizome were used for LC-MS based non-targeted metabolite profiling and targeted analysis of sugars, amino acids, picrosides and their corresponding phenolic acids.

RESULTS

A total of 220, primary and secondary metabolites (SMs) were identified (p < 0.05) representing an extensive inventory of metabolites and their spatial distribution in P. kurroa. Differential accumulation of metabolites suggests source-sink carbon partitioning, occurrence of partial TCA cycle, ascorbate metabolism, purine catabolism and salvage route, pyrimidine synthesis, lipid alteration besides gibberellins and cytokinin inhibition might be an adaptive strategy to alpine environmental stress along the altitude. Further, marked differences of organ and altitude specific SMs reflect alteration in secondary metabolic pathways. Significant accumulation of picrosides suggests their probable role in P. kurroa adaptation.

CONCLUSION

This study provides a platform that would be useful in deciphering the role of metabolites considered to be involved in plant adaptation.

Specific adaptations are selected in opposite sun exposed Antarctic cryptoendolithic communities as revealed by untargeted metabolomics

Antarctic cryptoendolithic communities are self-supporting borderline ecosystems spreading across the extreme conditions of the Antarctic desert and represent the predominant life-form in the ice-free areas of McMurdo Dry Valleys, accounted as the closest terrestrial Martian analogue. Components of these communities are highly adapted extremophiles and extreme-tolerant microorganisms, among the most resistant known to date. Recently, studies investigated biodiversity and community composition in these ecosystems but the metabolic activity of the metacommunity has never been investigated. Using an untargeted metabolomics, we explored stress-response of communities spreading in two sites of the same location, subjected to increasing environmental pressure due to opposite sun exposure, accounted as main factor influencing the diversity and composition of these ecosystems. Overall, 331 altered metabolites (206 and 125 unique for north and south, respectively), distinguished the two differently exposed communities. We also selected 10 metabolites and performed two-stage Receiver Operating Characteristic (ROC) analysis to test them as potential biomarkers. We further focused on melanin and allantoin as protective substances; their concentration was highly different in the community in the shadow or in the sun. These results clearly indicate that opposite insolation selected organisms in the communities with different adaptation strategies in terms of key metabolites produced.

Novel low-nitrogen stress-responsive long non-coding RNAs (lncRNA) in barley landrace B968 (Liuzhutouzidamai) at seedling stage

BACKGROUND

Reducing the dependence of crop production on chemical fertilizer with its associated costs, carbon footprint and other environmental problems is a challenge for agriculture. New solutions are required to solve this problem, and crop breeding for high nitrogen use efficiency or tolerance of low nitrogen availability has been widely considered to be a promising approach. However, the molecular mechanisms of high nitrogen use efficiency or low-nitrogen tolerance in crop plants are still to be elucidated, including the role of long non-coding RNAs (lncRNAs).

RESULTS

In this study, we identified 498 lncRNAs in barley (Hordeum vulgare) landrace B968 (Liuzhutouzidamai), of which 487 were novel, and characterised 56 that were responsive to low-nitrogen stress. For functional analysis of differentially-expressed lncRNAs, the gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment of co-expressed and co-located protein-coding genes were analyzed, and interactions with annotated co-expressed protein coding genes or micro RNAs (miRNAs) were further predicted. Target mimicry prediction between differentially-expressed lncRNAs and miRNAs identified 40 putative target mimics of lncRNAs and 58 target miRNAs. Six differentially-expressed lncRNAs were further validated by qPCR, and one in particular showed consistent differential expression using both techniques. Expression levels of most of the lncRNAs were found to be very low, and this may be the reason for the apparent inconsistency between RNA-seq and qPCR data.

CONCLUSIONS

The analysis of lncRNAs that are differentially-expressed under low-nitrogen stress, as well as their co-expressed or co-located protein coding genes and target mimics, could elucidate complex and hitherto uncharacterised mechanisms involved in the adaptation to low-nitrogen stress in barley and other crop plants.

Identification of Regulatory Networks and Hub Genes Controlling Nitrogen Uptake in Tea Plants [Camellia sinensis (L.) O. Kuntze]

Nitrogen (N) uptake, as the first step of N metabolism, is a key limiting factor for plant growth. To understand the gene expression networks that control N absorption and metabolism in tea plants, we analyzed transcriptomes in the young roots of two groups of tea plants with significantly different growth rates under different N treatments (0, 0.2, and 2 mM). Using pairwise comparisons and weighted gene co-expression network analyses (WGCNA), we successfully constructed 16 co-expression modules. Among them, a specific module (turquoise) that substantially responded to the low N treatment was identified. Based on KEGG analysis, the relative genes that enriched in the "N metabolism" pathways were used to construct gene co-expression networks of N metabolism. Finally, a high-affinity ammonium (NH₄⁺) transporter designated CsAMT1.2 was identified as a hub gene in the N metabolism network in tea plant roots and the gene expression could be highly induced by N resupply. The gene functional analysis revealed that CsAMT1.2 could make functional complementation of MEP1, MEP2, and MEP3 genes in 31019b yeast cells and improve NH₄⁺ uptake rate in 31019b at low NH₄⁺ level. Thus, CsAMT1.2 was a key gene controlling N uptake in tea plants and might play a vital role in promoting NH₄⁺ uptake from the environment in tea roots. This study provided a useful foundation for improving the NUE in tea plantations.

A Tripartite Interaction among the Basidiomycete Rhodotorula mucilaginosa, N2-Fixing Endobacteria, and Rice Improves Plant Nitrogen Nutrition

Nitrogen (N) limits crop yield, and improvement of N nutrition remains a key goal for crop research; one approach to improve N nutrition is identifying plant-interacting, N2-fixing microbes. Rhodotorula mucilaginosa JGTA-S1 is a basidiomycetous yeast endophyte of narrowleaf cattail (Typha angustifolia). JGTA-S1 could not convert nitrate or nitrite to ammonium but harbors diazotrophic (N2-fixing) endobacteria (Pseudomonas stutzeri) that allow JGTA-S1 to fix N2 and grow in a N-free environment; moreover, P. stutzeri dinitrogen reductase was transcribed in JGTA-S1 even under adequate N. Endobacteria-deficient JGTA-S1 had reduced fitness, which was restored by reintroducing P. stutzeri JGTA-S1 colonizes rice (Oryza sativa), significantly improving its growth, N content, and relative N-use efficiency. Endofungal P. stutzeri plays a significant role in increasing the biomass and ammonium content of rice treated with JGTA-S1; also, JGTA-S1 has better N2-fixing ability than free-living P. stutzeri and provides fixed N to the plant. Genes involved in N metabolism, N transporters, and NODULE INCEPTION-like transcription factors were upregulated in rice roots within 24 h of JGTA-S1 treatment. In association with rice, JGTA-S1 has a filamentous phase and P. stutzeri only penetrated filamentous JGTA-S1. Together, these results demonstrate an interkingdom interaction that improves rice N nutrition.

croFGD: Catharanthus roseus Functional Genomics Database

Catharanthus roseus is a medicinal plant, which can produce monoterpene indole alkaloid (MIA) metabolites with biological activity and is rich in vinblastine and vincristine. With release of the scaffolded genome sequence of C. roseus, it is necessary to annotate gene functions on the whole-genome level. Recently, 53 RNA-seq datasets are available in public with different tissues (flower, root, leaf, seedling, and shoot) and different treatments (MeJA, PnWB infection and yeast elicitor). We used in-house data process pipeline with the combination of PCC and MR algorithms to construct a co-expression network exploring multi-dimensional gene expression (global, tissue preferential, and treat response) through multi-layered approaches. In the meanwhile, we added miRNA-target pairs, predicted PPI pairs into the network and provided several tools such as gene set enrichment analysis, functional module enrichment analysis, and motif analysis for functional prediction of the co-expression genes. Finally, we have constructed an online croFGD database (http://bioinformatics.cau.edu.cn/croFGD/). We hope croFGD can help the communities to study the C. roseus functional genomics and make novel discoveries about key genes involved in some important biological processes.

Opposite fates of the purine metabolite allantoin under water and nitrogen limitations in bread wheat

Degradation of nitrogen-rich purines is tightly and oppositely regulated under drought and low nitrogen supply in bread wheat. Allantoin is a key target metabolite for improving nitrogen homeostasis under stress. The metabolite allantoin is an intermediate of the catabolism of purines (components of nucleotides) and is known for its housekeeping role in nitrogen (N) recycling and also for its function in N transport and storage in nodulated legumes. Allantoin was also shown to differentially accumulate upon abiotic stress in a range of plant species but little is known about its role in cereals. To address this, purine catabolic pathway genes were identified in hexaploid bread wheat and their chromosomal location was experimentally validated. A comparative study of two Australian bread wheat genotypes revealed a highly significant increase of allantoin (up to 29-fold) under drought. In contrast, allantoin significantly decreased (up to 22-fold) in response to N deficiency. The observed changes were accompanied by transcriptional adjustment of key purine catabolic genes, suggesting that the recycling of purine-derived N is tightly regulated under stress. We propose opposite fates of allantoin in plants under stress: the accumulation of allantoin under drought circumvents its degradation to ammonium (NH₄⁺) thereby preventing N losses. On the other hand, under N deficiency, increasing the NH₄⁺ liberated via allantoin catabolism contributes towards the maintenance of N homeostasis.

Overexpression of miR169o, an Overlapping MicroRNA in Response to Both Nitrogen Limitation and Bacterial Infection, Promotes Nitrogen Use Efficiency and Susceptibility to Bacterial Blight in Rice

Limiting nitrogen (N) supply contributes to improved resistance to bacterial blight (BB) caused by Xanthomonas oryzae pv. oryzae (Xoo) in susceptible rice (Oryza sativa). To understand the regulatory roles of microRNAs (miRNAs) in this phenomenon, 63 differentially expressed overlapping miRNAs in response to Xoo infection and N limitation stress in rice were identified through deep RNA sequencing and stem-loop quantitative real-time PCR. Among these, miR169o was further assessed as a typical overlapping miRNA through the overexpression of the miR169o primary gene. Osa-miR169o-OX plants were taller, and had more biomass accumulation with significantly increased nitrate and total amino acid contents in roots than the wild type (WT). Transcript level assays showed that under different N supply conditions, miR169o oppositely regulated NRT2, and this is reduced under normal N supply conditions but remarkably induced under N-limiting stress. On the other hand, osa-miR169o-OX plants also displayed increased disease lesion lengths and reduced transcriptional levels of defense gene (PR1b, PR10a, PR10b and PAL) compared with the WT after inoculation with Xoo. In addition, miR169o impeded Xoo-mediated NRT transcription. Therefore, the overlapping miR169o contributes to increase N use efficiency and negatively regulates the resistance to BB in rice. Consistently, transient expression of NF-YA genes in rice protoplasts promoted the transcripts of PR genes and NRT2 genes, while it reduced the transcripts of NRT1 genes. Our results provide novel and additional insights into the co ordinated regulatory mechanisms of cross-talk between Xoo infection and N deficiency responses in rice.

Transcriptome analysis for identifying possible gene regulations during maize root emergence and formation at the initial growth stage

The root plays an important role during plant development and growth, i.e., the plant body maintenance, nutrient storage, absorption of water, oxygen and nutrient from the soil, and storage of water and carbohydrates, etc. The objective of this study was attempted to determine root-specific genes at the initial developmental stages of maize by using network-based transcriptome analysis. The raw data obtained using RNA-seq were filtered for quality control of the reads with the FASTQC tool, and the filtered reads were pre-proceed using the TRIMMOMATIC tool. The enriched BINs of the DEGs were detected using PageMan analysis with the ORA_FISHER statistical test, and genes were assigned to metabolic pathways by using the MapMan tool, which was also used for detecting transcription factors (TFs). For reconstruction of the co-expression network, we used the algorithm for the reconstruction of accurate cellular networks (ARACNE) in the R package, and then the reconstructed co-expression network was visualized using the Cytoscape tool. RNA-seq. was performed using maize shoots and roots at different developmental stages of root emergence (6-10 days after planting, VE) and 1 week after plant emergence (V2). A total of 1286 differentially expressed genes (DEGs) were detected in both tissues. Many DEGs involved in metabolic pathways exhibited altered mRNA levels between VE and V2. In addition, we observed gene expression changes for 113 transcription factors and found five enriched cis-regulatory elements in the 1-kb upstream regions of both DEGs. The network-based transcriptome analysis showed two modules as co-expressed gene clusters differentially expressed between the shoots and roots during plant development. The DEGs of one module exhibited gene expressional coherence in the maize root tips, suggesting that their functional relationships are associated with the initial developmental stage of the maize root. Finally, we confirmed reliable mRNA levels of the hub genes in the potential sub-network related to initial root development at the different developmental stages of VE, V2, and 2 weeks after plant emergence.

Trehalose 6-Phosphate Regulates Photosynthesis and Assimilate Partitioning in Reproductive Tissue

Transgenic maize (Zea mays) that expresses rice (Oryza sativa) TREHALOSE PHOSPHATE PHOSPHATASE1 (TPP1) from the rice MADS6 promoter, which is active over the flowering period, produces higher yields than wild type. This yield increase occurs with or without drought conditions during flowering. To understand the mechanistic basis of the increased yield, we characterized gene expression and metabolite profiles in leaves and developing female reproductive tissue, comprising florets, node, pith, and shank, over the flowering period with and without drought. The MADS6 promoter was most active in the vasculature, particularly phloem companion cells in florets and pith, consistent with the largest decreases in trehalose 6-phosphate (T6P) levels (2- to 3-fold) being found in pith and florets. Low T6P led to decreased gene expression for primary metabolism and increased gene expression for secondary metabolism, particularly lipid-related pathways. Despite similar changes in gene expression, the pith and floret displayed opposing assimilate profiles: sugars, sugar phosphates, amino acids, and lipids increased in florets, but decreased in pith. Possibly explaining this assimilate distribution, seven SWEET genes were found to be up-regulated in the transgenic plants. SnRK1 activity and the expression of the gene for the SnRK1 beta subunit, expression of SnRK1 marker genes, and endogenous trehalose pathway genes were also altered. Furthermore, leaves of the transgenic maize maintained a higher photosynthetic rate for a longer period compared to wild type. In conclusion, we found that decreasing T6P in reproductive tissues down-regulates primary metabolism and up-regulates secondary metabolism, resulting in different metabolite profiles in component tissues. Our data implicate T6P/ SnRK1 as a major regulator of whole-plant resource allocation for crop yield improvement.

Identification and Analysis of Rice Yield-Related Candidate Genes by Walking on the Functional Network

Rice (Oryza sativa L.) is one of the most important staple foods in the world. It is possible to identify candidate genes associated with rice yield using the model of random walk with restart on a functional similarity network. We demonstrated the high performance of this approach by a five-fold cross-validation experiment, as well as the robustness of the parameter r. We also assessed the strength of associations between known seeds and candidate genes in the light of the results scores. The candidates ranking at the top of the results list were considered to be the most relevant rice yield-related genes. This study provides a valuable alternative for rice breeding and biology research. The relevant dataset and script can be downloaded at the website: http://lab.malab.cn/jj/rice.htm.

Comparative transcriptome analyses of flower development in four species of Achimenes (Gesneriaceae)

BACKGROUND

Flowers have an amazingly diverse display of colors and shapes, and these characteristics often vary significantly among closely related species. The evolution of diverse floral form can be thought of as an adaptive response to pollination and reproduction, but it can also be seen through the lens of morphological and developmental constraints. To explore these interactions, we use RNA-seq across species and development to investigate gene expression and sequence evolution as they relate to the evolution of the diverse flowers in a group of Neotropical plants native to Mexico-magic flowers (Achimenes, Gesneriaceae).

RESULTS

The assembled transcriptomes contain between 29,000 and 42,000 genes expressed during development. We combine sequence orthology and coexpression clustering with analyses of protein evolution to identify candidate genes for roles in floral form evolution. Over 25% of transcripts captured were distinctive to Achimenes and overrepresented by genes involved in transcription factor activity. Using a model-based clustering approach we find dynamic, temporal patterns of gene expression among species. Selection tests provide evidence of positive selection in several genes with roles in pigment production, flowering time, and morphology. Combining these approaches to explore genes related to flower color and flower shape, we find distinct patterns that correspond to transitions of floral form among Achimenes species.

CONCLUSIONS

The floral transcriptomes developed from four species of Achimenes provide insight into the mechanisms involved in the evolution of diverse floral form among closely related species with different pollinators. We identified several candidate genes that will serve as an important and useful resource for future research. High conservation of sequence structure, patterns of gene coexpression, and detection of positive selection acting on few genes suggests that large phenotypic differences in floral form may be caused by genetic differences in a small set of genes. Our characterized floral transcriptomes provided here should facilitate further analyses into the genomics of flower development and the mechanisms underlying the evolution of diverse flowers in Achimenes and other Neotropical Gesneriaceae.

Allantoin, a stress-related purine metabolite, can activate jasmonate signaling in a MYC2-regulated and abscisic acid-dependent manner

Allantoin is a metabolic intermediate of purine catabolism that often accumulates in stressed plants. Recently, we used Arabidopsis knockout mutants (aln) of ALLANTOINASE to show that this purine metabolite activates abscisic acid (ABA) production, thereby stimulating stress-related gene expression and enhancing seedling tolerance to abiotic stress. A detailed re-examination of the microarray data of an aln mutant (aln-1) confirmed the increased expression of ABA-related genes and also revealed altered expression of genes involved in jasmonic acid (JA) responses, probably under the control of MYC2, a master switch in the JA signaling pathway. Consistent with the transcriptome profiles, the aln-1 mutant displayed increased JA levels and enhanced responses to mechanical wounding and exogenous JA. Moreover, aln mutants demonstrated modestly increased susceptibility to Pseudomonas syringae and Pectobacterium carotovorum, probably reflecting the antagonistic action of MYC2 on the defense against these bacterial phytopathogens. Exogenously administered allantoin elicited the expression of JA-responsive genes, including MYC2, in wild-type plants, supporting the idea that allantoin might be responsible for the observed JA-related phenotypes of aln mutants. However, mutants deficient in bioactive JA (jar1-1), insensitive to JA (myc2-3), or deficient in ABA (aba2-1 and bglu18) suppressed the effect of exogenous allantoin. The suppression was further confirmed in aln-1 jar1-1 and aln-1 bglu18 double mutants. These results indicate that allantoin can activate the MYC2-regulated JA signaling pathway through ABA production. Overall, this study suggests a possible connection of purine catabolism with stress hormone homeostasis and signaling, and highlights the potential importance of allantoin in these interactions.