Identification of genes associated with nitrogen-use efficiency by genome-wide transcriptional analysis of two soybean genotypes

Insights into the Early Steps of the Symbiotic Interaction between Soybean (Glycine max) and Sinorhizobium fredii Symbiosis Using Transcriptome, Small RNA, and Degradome Sequencing

Symbiotic nitrogen fixation carried out by the soybean-rhizobia symbiosis increases soybean yield and reduces the amount of nitrogen fertilizer that has been applied. MicroRNAs (miRNAs) are crucial in plant growth and development, prompting an investigation into their role in the symbiotic interaction of soybean with partner rhizobia. Through integrated small RNA, transcriptome, and degradome sequencing analysis, 1215 known miRNAs, 314 of them conserved, and 187 novel miRNAs were identified, with 44 differentially expressed miRNAs in soybean roots inoculated with Sinorhizobium fredii HH103 and a ttsI mutant. The study unveiled that the known miRNA gma-MIR398a-p5 was downregulated in the presence of the ttsI mutation, while the target gene of gma-MIR398a-p5, Glyma.06G007500, associated with nitrogen metabolism, was upregulated. The results of this study offer insights for breeding high-efficiency nitrogen-fixing soybean varieties, enhancing crop yield and quality.

The mechanism of melatonin promotion on cucumber seedling growth at different nitrogen levels

The supply level of exogenous nitrogen has a very important influence on the growth and development of cucumber. Insufficient or excessive nitrogen application will lead to metabolic disorders in the body and affect the formation of yield. Therefore, it is of great scientific and practical significance to explore the corresponding mitigation measures. Melatonin (MT) is a multi-regulatory molecule with pleiotropic effects on plant growth and development. A large number of studies have shown that the appropriate amount of melatonin supplementation is beneficial to plant growth and development by promoting root development, delaying leaf senescence, and improving fruit yield. However, the study of MT function combined with a detailed physiological analysis of nitrogen (N) absorption and metabolism in cucumber plants needs further strengthening. We performed hydroponic tests at different nitrogen levels to determine the metabolic processes associated with the enhanced tolerance to nitrogen in melatonin-treated cucumber (Cucucumis sativus L.) seedlings. Cucumber seedlings were sprayed with 100 μM melatonin or water and treated with different nitrogen in the growth chamber for 7 days. Nitrogen deficiency significantly inhibited seedling growth, and this growth inhibition was partially alleviated by melatonin. The expression analysis of related carbon and nitrogen genes showed that the genes whose expression was significantly altered by melatonin were mainly related to carbon (C) and nitrogen (N) metabolism. By enzyme activity and reactive oxygen content data analysis, melatonin-treated cucumber seedlings showed relatively stable carbon and nitrogen levels compared to untreated ones. In conclusion, MT can repair the impaired growth and development situation by regulating the nitrogen assimilation capacity and the balance between oxidation and oxidative metabolism and carbon metabolism in the cucumber under different nitrogen levels.

An integrated nitrogen utilization gene network and transcriptome analysis reveal candidate genes in response to nitrogen deficiency in Brassica napus

Nitrogen (N) is an essential factor for crop yield. Here, we characterized 605 genes from 25 gene families that form the complex gene networks of N utilization pathway in Brassica napus. We found unequal gene distribution between the A_n- and C_n-sub-genomes, and that genes derived from Brassica rapa were more retained. Transcriptome analysis indicated that N utilization pathway gene activity shifted in a spatio-temporal manner in B. napus. A low N (LN) stress RNA-seq of B. napus seedling leaves and roots was generated, which proved that most N utilization related genes were sensitive to LN stress, thereby forming co-expression network modules. Nine candidate genes in N utilization pathway were confirmed to be significantly induced under N deficiency conditions in B. napus roots, indicating their potential roles in LN stress response process. Analyses of 22 representative species confirmed that the N utilization gene networks were widely present in plants ranging from Chlorophyta to angiosperms with a rapid expansion trend. Consistent with B. napus, the genes in this pathway commonly showed a wide and conserved expression profile in response to N stress in other plants. The network, genes, and gene-regulatory modules identified here represent resources that may enhance the N utilization efficiency or the LN tolerance of B. napus.

Short-term transcriptomic analysis at organ scale reveals candidate genes involved in low N responses in NUE-contrasting tomato genotypes

BACKGROUND

Understanding the complex regulatory network underlying plant nitrogen (N) responses associated with high Nitrogen Use Efficiency (NUE) is one of the main challenges for sustainable cropping systems. Nitrate (NO₃ ^-), acting as both an N source and a signal molecule, provokes very fast transcriptome reprogramming, allowing plants to adapt to its availability. These changes are genotype- and tissue-specific; thus, the comparison between contrasting genotypes is crucial to uncovering high NUE mechanisms.

METHODS

Here, we compared, for the first time, the spatio-temporal transcriptome changes in both root and shoot of two NUE contrasting tomato genotypes, Regina Ostuni (high-NUE) and UC82 (low-NUE), in response to short-term (within 24 h) low (LN) and high (HN) NO₃ ^- resupply.

RESULTS

Using time-series transcriptome data (0, 8, and 24 h), we identified 395 and 482 N-responsive genes differentially expressed (DEGs) between RO and UC82 in shoot and root, respectively. Protein kinase signaling plant hormone signal transduction, and phenylpropanoid biosynthesis were the main enriched metabolic pathways in shoot and root, respectively, and were upregulated in RO compared to UC82. Interestingly, several N transporters belonging to NRT and NPF families, such as NRT2.3, NRT2.4, NPF1.2, and NPF8.3, were found differentially expressed between RO and UC82 genotypes, which might explain the contrasting NUE performances. Transcription factors (TFs) belonging to several families, such as ERF, LOB, GLK, NFYB, ARF, Zinc-finger, and MYB, were differentially expressed between genotypes in response to LN. A complementary Weighted Gene Co-expression Network Analysis (WGCNA) allowed the identification of LN-responsive co-expression modules in RO shoot and root. The regulatory network analysis revealed candidate genes that might have key functions in short-term LN regulation. In particular, an asparagine synthetase (ASNS), a CBL-interacting serine/threonine-protein kinase 1 (CIPK1), a cytokinin riboside 5'-monophosphate phosphoribohydrolase (LOG8), a glycosyltransferase (UGT73C4), and an ERF2 were identified in the shoot, while an LRR receptor-like serine/threonine-protein kinase (FEI1) and two TFs NF-YB5 and LOB37 were identified in the root.

DISCUSSION

Our results revealed potential candidate genes that independently and/or concurrently may regulate short-term low-N response, suggesting a key role played by cytokinin and ROS balancing in early LN regulation mechanisms adopted by the N-use efficient genotype RO.

Integrated Physiological, Transcriptomic, and Metabolomic Analyses of the Response of Peach to Nitrogen Levels during Different Growth Stages

This study performed physiological, transcriptome, and metabolite analyses of peach fruit under different nitrogen (N) conditions at different growth stages. Nitrogen management directly affected the yield, fruit quality, and metabolites of peach in different growth stages. Different fertilizing time influenced yield and leaf N concentration. RNA-Seq was used to analyze the influence of N levels at the fruit pit hardening (PH) and fruit expansion (FE) stages. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed differentially expressed genes (DEGs) related to carbon and nitrogen metabolite processes. Metabolome analysis shows that applying different nitrogen fertilizers at different growth stages of peach mainly affected metabolites related to carbon and amino acids. This research provides insight into the metabolic processes underlying different N responses during different growth stages and provides a foundation to improve the efficiency of N use in peach.

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.

Physiological, Agronomical, and Proteomic Studies Reveal Crucial Players in Rice Nitrogen Use Efficiency under Low Nitrogen Supply

Excessive use of nitrogenous fertilizers to enhance rice productivity has become a significant source of nitrogen (N) pollution and reduced sustainable agriculture. However, little information about the physiology of different growth stages, agronomic traits, and associated genetic bases of N use efficiency (NUE) are available at low-N supply. Two rice (Oryza sativa L.) cultivars were grown with optimum N (120 kg ha-1) and low N (60 kg ha-1) supply. Six growth stages were analyzed to measure the growth and physiological traits, as well as the differential proteomic profiles, of the rice cultivars. Cultivar Panvel outclassed Nagina 22 at low-N supply and exhibited improved growth and physiology at most of the growth stages and agronomic efficiency due to higher N uptake and utilization at low-N supply. On average, photosynthetic rate, chlorophyll content, plant biomass, leaf N content, and grain yield were decreased in cultivar Nagina 22 than Panvel was 8%, 11%, 21%, 19%, and 22%, respectively, under low-N supply. Furthermore, proteome analyses revealed that many proteins were upregulated and downregulated at the different growth stages under low-N supply. These proteins are associated with N and carbon metabolism and other physiological processes. This supports the genotypic differences in photosynthesis, N assimilation, energy stabilization, and rice-protein yield. Our study suggests that enhancing NUE at low-N supply demands distinct modifications in N metabolism and physiological assimilation. The NUE may be regulated by key identified differentially expressed proteins. These proteins might be the targets for improving crop NUE at low-N supply.

RNA-Seq-Based Transcriptomics Study to Investigate the Genes Governing Nitrogen Use Efficiency in Indian Wheat Cultivars

High NUE (nitrogen use efficiency) has great practical significance for sustainable crop production. Wheat is one of the main cultivated crops worldwide for human food and nutrition. However, wheat grain productivity is dependent upon cultivars with high NUE in addition to the application of nitrogen fertilizers. In order to understand the molecular mechanisms exhibiting a high NUE response, a comparative transcriptomics study was carried out through RNA-seq analysis to investigate the gene expression that regulates NUE, in root and shoot tissue of N-efficient (PBW677) and N-inefficient (703) cultivars under optimum and nitrogen (N) stress. Differentially expressed gene analysis revealed a total of 2,406 differentially expressed genes (DEGs) present in both the contrasting cultivars under N stress. The efficient genotype PBW677 had considerably more abundant DEGs with 1,653 (903 roots +750 shoots) compared to inefficient cultivar PBW703 with 753 (96 roots +657 shoots). Gene ontology enrichment and pathway analysis of these DEGs suggested that the two cultivars differed in terms of adaptive mechanism. Gene enrichment analysis revealed that among the upregulated and downregulated genes the overrepresented and underrepresented gene categories belonged to biological processes like DNA binding, response to abiotic stimulus, photosynthesis, carbon fixation, carbohydrate metabolic process, nitrogen compound metabolic process, nitrate transport, and translation in cultivar PBW677, while the enriched biological processes were nucleosome assembly, chromatin remodeling, DNA packaging, lipid transport, sulfur compound metabolic process, protein modifications, and protein folding and refolding in N inefficient cultivar PBW703. We found several transcription factors (MYB, WRKY, RING finger protein, zinc finger protein, transporters, NRT1, amino acid transporters, sugar), protein kinases, and genes involved in N absorption, transportation, and assimilation to be highly expressed in high NUE cultivar PBW677. In our study, we report 13 potential candidate genes which showed alternate gene expression in the two contrasting cultivars under study. These genes could serve as potential targets for future breeding programs.

Global gene expression profiling under nitrogen stress identifies key genes involved in nitrogen stress adaptation in maize (Zea mays L.)

Maize is a heavy consumer of fertilizer nitrogen (N) which not only results in the high cost of cultivation but may also lead to environmental pollution. Therefore, there is a need to develop N-use efficient genotypes, a prerequisite for which is a greater understanding of N-deficiency stress adaptation. In this study, comparative transcriptome analysis was performed using leaf and root tissues from contrasting inbred lines, viz., DMI 56 (tolerant to N stress) and DMI 81 (susceptible to N stress) to delineate the differentially expressed genes (DEGs) under low-N stress. The contrasting lines were grown hydroponically in modified Hoagland solution having either sufficient- or deficient-N, followed by high-throughput RNA-sequencing. A total of 8 sequencing libraries were prepared and 88–97% of the sequenced raw reads were mapped to the reference B73 maize genome. Genes with a p value ≤ 0.05 and fold change of ≥ 2.0 or ≤ − 2 were considered as DEGs in various combinations performed between susceptible and tolerant genotypes. DEGs were further classified into different functional categories and pathways according to their putative functions. Gene Ontology based annotation of these DEGs identified three different functional categories: biological processes, molecular function, and cellular component. The KEGG and Mapman based analysis revealed that most of the DEGs fall into various metabolic pathways, biosynthesis of secondary metabolites, signal transduction, amino acid metabolism, N-assimilation and metabolism, and starch metabolism. Some of the key genes involved in N uptake (high-affinity nitrate transporter 2.2 and 2.5), N assimilation and metabolism (glutamine synthetase, asparagine synthetase), redox homeostasis (SOD, POX), and transcription factors (MYB36, AP2-EREBP) were found to be highly expressed in the tolerant genotype compared to susceptible one. The candidate genes identified in the present study might be playing a pivotal role in low-N stress adaptation in maize and hence could be useful in augmenting further research on N metabolism and development of N-deficiency tolerant maize cultivars.

The mechanisms underlying melatonin improved soybean seedling growth at different nitrogen levels

To investigate the function of melatonin (MT) on nitrogen uptake and metabolism in soybean, six groups of treatments, with and without 100μM melatonin were conducted at low, normal, and high nitrogen levels (1.5, 7.5, and 15mM, respectively). The related indexes of nitrogen metabolism and the antioxidant system of seedlings were measured and analysed. Results indicated that MT could enhance the level of nitrogen metabolism by upregulating the coding genes of enzymes related to nitrogen metabolism and increasing total nitrogen content, especially under low nitrogen levels. Under high nitrogen conditions, the addition of MT not only accelerated ammonium assimilation and utilisation by enhancing the activity of glutamine synthetase involved in ammonium assimilation, but also reduced the extent of membrane lipid peroxidation to alleviate the degree of damage by improving the activity of antioxidant enzymes. In addition, MT enhanced soybean growth with positive effects in morphological changes at different nitrogen levels, including significantly increased stem diameter, total leaf area, and root nodule number, and biomass accumulation. Finally, biomass accumulation increased under low, normal, and high nitrogen levels by 9.80%, 14.06%, and 11.44%, respectively. The results suggested that MT could enhance the soybean tolerance to low and excessive N treatments.

Integrative Transcriptomic and Proteomic Analysis Reveals an Alternative Molecular Network of Glutamine Synthetase 2 Corresponding to Nitrogen Deficiency in Rice (Oryza sativa L.)

Nitrogen (N) is an essential nutrient for plant growth and development. The root system architecture is a highly regulated morphological system, which is sensitive to the availability of nutrients, such as N. Phenotypic characterization of roots from LY9348 (a rice variety with high nitrogen use efficiency (NUE)) treated with 0.725 mM NH₄NO₃ (1/4N) was remarkable, especially primary root (PR) elongation, which was the highest. A comprehensive analysis was performed for transcriptome and proteome profiling of LY9348 roots between 1/4N and 2.9 mM NH₄NO₃ (1N) treatments. The results indicated 3908 differential expression genes (DEGs; 2569 upregulated and 1339 downregulated) and 411 differential abundance proteins (DAPs; 192 upregulated and 219 downregulated). Among all DAPs in the proteome, glutamine synthetase (GS2), a chloroplastic ammonium assimilation protein, was the most upregulated protein identified. The unexpected concentration of GS2 from the shoot to the root in the 1/4N treatment indicated that the presence of an alternative pathway of N assimilation regulated by GS2 in LY9348 corresponded to the low N signal, which was supported by GS enzyme activity and glutamine/glutamate (Gln/Glu) contents analysis. In addition, N transporters (NRT2.1, NRT2.2, NRT2.3, NRT2.4, NAR2.1, AMT1.3, AMT1.2, and putative AMT3.3) and N assimilators (NR2, GS1;1, GS1;2, GS1;3, NADH-GOGAT2, and AS2) were significantly induced during the long-term N-deficiency response at the transcription level (14 days). Moreover, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis demonstrated that phenylpropanoid biosynthesis and glutathione metabolism were significantly modulated by N deficiency. Notably, many transcription factors and plant hormones were found to participate in root morphological adaptation. In conclusion, our study provides valuable information to further understand the response of rice roots to N-deficiency stress.

Transcriptomic, proteomic, and physiological studies reveal key players in wheat nitrogen use efficiency under both high and low nitrogen supply

The effective use of available nitrogen (N) to improve crop grain yields provides an important strategy to reduce environmental N pollution and promote sustainable agriculture. However, little is known about the common genetic basis of N use efficiency (NUE) at varying N availability. Two wheat (Triticum aestivum L.) cultivars were grown in the field with high, moderate, and low N supply. Cultivar Zhoumai 27 outperformed Aikang 58 independent of the N supply and showed improved growth, canopy leaf area index, flag leaf surface area, grain number, and yield, and enhanced NUE due to both higher N uptake and utilization efficiency. Further, transcriptome and proteome analyses were performed using flag leaves that provide assimilates for grain growth. The results showed that many genes or proteins that are up- or down-regulated under all N regimes are associated with N and carbon metabolism and transport. This was reinforced by cultivar differences in photosynthesis, assimilate phloem transport, and grain protein/starch yield. Overall, our study establishes that improving NUE at both high and low N supply requires distinct adjustments in leaf metabolism and assimilate partitioning. Identified key genes/proteins may individually or concurrently regulate NUE and are promising targets for maximizing crop NUE irrespective of the N supply.

Transcriptomics reveal new insights into molecular regulation of nitrogen use efficiency in Solanum melongena

Nitrogen-use efficiency (NUE) is a complex trait of great interest in breeding programs because through its improvement, high crop yields can be maintained whilst N supply is reduced. In this study, we report a transcriptomic analysis of four NUE-contrasting eggplant (Solanum melongena) genotypes following short- and long-term exposure to low N, to identify key genes related to NUE in the roots and shoots. The differentially expressed genes in the high-NUE genotypes are involved in the light-harvesting complex and receptor, a ferredoxin-NADP reductase, a catalase and WRKY33. These genes were then used as bait for a co-expression gene network analysis in order to identify genes with the same trends in expression. This showed that up-regulation of WRKY33 triggered higher expression of a cluster of 21 genes and also of other genes, many of which were related to N-metabolism, that were able to improve both nitrogen uptake efficiency and nitrogen utilization efficiency, the two components of NUE. We also conducted an independent de novo experiment to validate the significantly higher expression of WRKY33 and its gene cluster in the high-NUE genotypes. Finally, examination of an Arabidopsis transgenic 35S::AtWRKY33 overexpression line showed that it had a bigger root system and was more efficient at taking up N from the soil, confirming the pivotal role of WRKY33 for NUE improvement.

Comparative Analysis Based on Transcriptomics and Metabolomics Data Reveal Differences between Emmer and Durum Wheat in Response to Nitrogen Starvation

Mounting evidence indicates the key role of nitrogen (N) on diverse processes in plant, including development and defense. Using a combined transcriptomics and metabolomics approach, we studied the response of seedlings to N starvation of two different tetraploid wheat genotypes from the two main domesticated subspecies: emmer and durum wheat. We found that durum wheat exhibits broader and stronger response in comparison to emmer as seen from the expression pattern of both genes and metabolites and gene enrichment analysis. They showed major differences in the responses to N starvation for transcription factor families, emmer showed differential reduction in the levels of primary metabolites while durum wheat exhibited increased levels of most of them to N starvation. The correlation-based networks, including the differentially expressed genes and metabolites, revealed tighter regulation of metabolism in durum wheat in comparison to emmer. We also found that glutamate and γ-aminobutyric acid (GABA) had highest values of centrality in the metabolic correlation network, suggesting their critical role in the genotype-specific response to N starvation of emmer and durum wheat, respectively. Moreover, this finding indicates that there might be contrasting strategies associated to GABA and glutamate signaling modulating shoot vs. root growth in the two different wheat subspecies.

Transcriptome Analysis for Fraxinus mandshurica Rupr. Seedlings from Different Carbon Sequestration Provenances in Response to Nitrogen Deficiency

To explore the molecular regulatory mechanism of high-carbon (C) sequestration Fraxinus mandshurica Rupr. (F. mandshurica) provenance and the expression profile of F. mandshurica during nitrogen (N) starvation, the foliage and roots of the annual Wuchang (WC) seedlings with greater C amount and Hailin (HL) seedlings with smaller C amount, which were grown in N-deficient nutrition and complete N, were used for RNA-seq and physiological determination, respectively. One thousand and fifty-seven differentially expressed genes (DEGs) between WC and HL and 8173 DEGs related to N deficiency were identified, respectively. The root of F. mandshurica responded to N deficiency more strongly than foliar. The target genes that responded to N deficiency in roots were mainly regulatory genes (transcription factors, hormones and protein kinases), and their response patterns were upregulated. The growth and N concentration in both WC and HL were reduced by the N deficiency, which might result from the decrease of the leaf Nitrate reductase (NR) and glutamine synthetase (GS) enzyme activity and ABA content, although the root-to-shoot ratio; lateral root number; lignin content; endogenous hormones content (GA, IAA and ZR); root GS and glutamate synthetase activity and transcriptional level of most of the regulatory genes were increased. The C sequestration capacity in WC was greater than that in HL, which related to the higher GS enzymes activity and transcriptional levels of regulatory genes and metabolic genes (terpenes, carbohydrates, and lipid energy). However, the C sequestration advantage of WC was significantly reduced by the N deficiency, which was due to the smaller response to N deficiency compared to HL.

Transcriptome-based approaches for clarification of nutritional responses and improvement of crop production

Genome-wide transcriptome profiling is a powerful tool for identifying key genes and pathways involved in plant development and physiological processes. This review summarizes studies that have used transcriptome profiling mainly in rice to focus on responses to macronutrients such as nitrogen, phosphorus and potassium, and spatio-temporal root profiling in relation to the regulation of root system architecture as well as nutrient uptake and transport. We also discuss strategies based on meta- and co-expression analyses with different attributed transcriptome data, which can be used for investigating the regulatory mechanisms and dynamics of nutritional responses and adaptation, and speculate on further advances in transcriptome profiling that could have potential application to crop breeding and cultivation.

Association genetics of the parameters related to nitrogen use efficiency in Brassica juncea L

Genome wide association studies allowed prediction of 17 candidate genes for association with nitrogen use efficiency. Novel information obtained may provide better understanding of genomic controls underlying germplasm variations for this trait in Indian mustard. Nitrogen use efficiency (NUE) of Indian mustard (Brassica juncea (L.) Czern & Coss.) is low and most breeding efforts to combine NUE with crop performance have not succeeded. Underlying genetics also remain unexplored. We tested 92 SNP-genotyped inbred lines for yield component traits, N uptake efficiency (NUPEFF), nitrogen utilization efficiency (NUTEFF), nitrogen harvest index (NHI) and NUE for two years at two nitrogen doses (N_o without added N and N₁₀₀ added @100 kg/ha). Genotypes IC-2489-88, M-633, MCP-632, HUJM 1080, GR-325 and DJ-65 recorded high NUE at low N. These also showed improved crop performance under high N. One determinate mustard genotype DJ-113 DT-3 revealed maximum NUTEFF. Genome wide association studies (GWAS) facilitated recognition of 17 quantitative trait loci (QTLs). Environment specificity was high. B-genome chromosomes (B02, B03, B05, B07 and B08) harbored many useful loci. We also used regional association mapping (RAM) to supplement results from GWAS. Annotation of the genomic regions around peak SNPs helped to predict several gene candidates for root architecture, N uptake, assimilation and remobilization. CAT9 (At1g05940) was consistently envisaged for both NUE and NUPEFF. Major N transporter genes, NRT1.8 and NRT3.1 were predicted for explaining variation for NUTEFF and NUPEFF, respectively. Most significant amino acid transporter gene, AAP1 appeared associated with NUE under limited N conditions. All these candidates were predicted in the regions of high linkage disequilibrium. Sequence information of the predicted candidate genes will permit development of molecular markers to aid breeding for high NUE.

Fate of nitrogen in agriculture and environment: agronomic, eco-physiological and molecular approaches to improve nitrogen use efficiency

Nitrogen is the main limiting nutrient after carbon, hydrogen and oxygen for photosynthetic process, phyto-hormonal, proteomic changes and growth-development of plants to complete its lifecycle. Excessive and inefficient use of N fertilizer results in enhanced crop production costs and atmospheric pollution. Atmospheric nitrogen (71%) in the molecular form is not available for the plants. For world's sustainable food production and atmospheric benefits, there is an urgent need to up-grade nitrogen use efficiency in agricultural farming system. The nitrogen use efficiency is the product of nitrogen uptake efficiency and nitrogen utilization efficiency, it varies from 30.2 to 53.2%. Nitrogen losses are too high, due to excess amount, low plant population, poor application methods etc., which can go up to 70% of total available nitrogen. These losses can be minimized up to 15-30% by adopting improved agronomic approaches such as optimal dosage of nitrogen, application of N by using canopy sensors, maintaining plant population, drip fertigation and legume based intercropping. A few transgenic studies have shown improvement in nitrogen uptake and even increase in biomass. Nitrate reductase, nitrite reductase, glutamine synthetase, glutamine oxoglutarate aminotransferase and asparagine synthetase enzyme have a great role in nitrogen metabolism. However, further studies on carbon-nitrogen metabolism and molecular changes at omic levels are required by using "whole genome sequencing technology" to improve nitrogen use efficiency. This review focus on nitrogen use efficiency that is the major concern of modern days to save economic resources without sacrificing farm yield as well as safety of global environment, i.e. greenhouse gas emissions, ammonium volatilization and nitrate leaching.

Genome-Wide Association Study and Identification of Candidate Genes for Nitrogen Use Efficiency in Barley (Hordeum vulgare L.)

Nitrogen (N) fertilizer is largely responsible for barley grain yield potential and quality, yet excessive application leads to environmental pollution and high production costs. Therefore, efficient use of N is fundamental for sustainable agriculture. In the present study, we investigated the performance of 282 barley accessions through hydroponic screening using optimal and low NH4NO3 treatments. Low-N treatment led to an average shoot dry weight reduction of 50%, but there were significant genotypic differences among the accessions. Approximately 20% of the genotypes showed high (>75%) relative shoot dry weight under low-N treatment and were classified as low-N tolerant, whereas 20% were low-N sensitive (≤55%). Low-N tolerant accessions exhibited well-developed root systems with an average increase of 60% in relative root dry weight to facilitate more N absorption. A genome-wide association study (GWAS) identified 66 significant marker trait associations (MTAs) conferring high nitrogen use efficiency, four of which were stable across experiments. These four MTAs were located on chromosomes 1H(1), 3H(1), and 7H(2) and were associated with relative shoot length, relative shoot and root dry weight. Genes corresponding to the significant MTAs were retrieved as candidate genes, including members of the asparagine synthetase gene family, several transcription factor families, protein kinases, and nitrate transporters. Most importantly, the high-affinity nitrate transporter 2.7 (HvNRT2.7) was identified as a promising candidate on 7H for root and shoot dry weight. The identified candidate genes provide new insights into our understanding of the molecular mechanisms driving nitrogen use efficiency in barley and represent potential targets for genetic improvement.

A GATA Transcription Factor from Soybean (Glycine max) Regulates Chlorophyll Biosynthesis and Suppresses Growth in the Transgenic Arabidopsis thaliana

Chlorophyll plays an essential role in photosynthetic light harvesting and energy transduction in green tissues of higher plants and is closely related to photosynthesis and crop yield. Identification of transcription factors (TFs) involved in regulating chlorophyll biosynthesis is still limited in soybean (Glycine max), and the previously identified GmGATA58 is suggested to potentially modulate chlorophyll and nitrogen metabolisms, but its complete function is still unknown. In this study, subcellular localization assay showed that GmGATA58 was localized in the nucleus. Histochemical GUS assay and qPCR assay indicated that GmGATA58 was mainly expressed in leaves and responded to nitrogen, light and phytohormone treatments. Overexpression of GmGATA58 in the Arabidopsis thaliana ortholog AtGATA21 (gnc) mutant complemented the greening defect, while overexpression in Arabidopsis wild-type led to increasing chlorophyll content in leaves through up-regulating the expression levels of the large of chlorophyll biosynthetic pathway genes, but suppressing plant growth and yield, although the net photosynthetic rate was slightly improved. Dual-luciferase reporter assay also supported that GmGATA58 activated the transcription activities of three promoters of key chlorophyll biosynthetic genes of soybean in transformed protoplast of Arabidopsis. It is concluded that GmGATA58 played an important role in regulating chlorophyll biosynthesis, but suppressed plant growth and yield in transgenic Arabidopsis.

Transcriptome Analysis of High-NUE (T29) and Low-NUE (T13) Genotypes Identified Different Responsive Patterns Involved in Nitrogen Stress in Ramie (Boehmeria nivea (L.) Gaudich)

Nitrogen-use efficiency (NUE) has significant impacts on plant growth and development. NUE in plants differs substantially in physiological resilience to nitrogen stress; however, the molecular mechanisms underlying enhanced resilience of high-NUE plants to nitrogen deficiency remains unclear. We compared transcriptome-wide gene expression between high-NUE and low-NUE ramie (Boehmeria nivea (L.) Gaudich) genotypes under nitrogen (N)-deficient and normal conditions to identify the transcriptomic expression patterns that contribute to ramie resilience to nitrogen deficiency. Two ramie genotypes with contrasting NUE were used in the study, including T29 (NUE = 46.01%) and T13 (NUE = 15.81%). Our results showed that high-NUE genotypes had higher gene expression under the control condition across 94 genes, including frontloaded genes such as GDSL esterase and lipase, gibberellin, UDP-glycosyltransferase, and omega-6 fatty acid desaturase. Seventeen stress-tolerance genes showed lower expression levels and varied little in response to N-deficiency stress in high-NUE genotypes. In contrast, 170 genes were upregulated under N deficiency in high-NUE genotypes but downregulated in low-NUE genotypes compared with the controls. Furthermore, we identified the potential key genes that enable ramie to maintain physiological resilience under N-deficiency stress, and categorized these genes into three groups based on the transcriptome and their expression patterns. The transcriptomic and clustering analysis of these nitrogen-utilization-related genes could provide insight to better understand the mechanism of linking among the three gene classes that enhance resilience in high-NUE ramie genotypes.

High-throughput sequencing reveals genes associated with high-temperature stress tolerance in sugarcane

Sugarcane (Saccharum spp.) is a major sugar crop grown in tropical and sub-tropical areas throughout the world which is vulnerable to high temperature stress due to climate change. In this present study, we have generated a transcriptome profile of sugarcane variety Co 99004 exposed to high-temperature stress (47 ^°C). The Illumina Nexseq2500 platform yielded a total of 39.28 and 13.44 million reads, corresponding to 3.9 and 1.3 gigabase pair (Gb) of the processed reads for control and high-temperature-stressed samples, respectively. Initially, the reads were de novo assembled into 118,017 unigenes with an average length of 780 bp. The longest sequence in the assembly was 21 kb. Further, these transcripts were BLASTed against GO, KEGG and COG databases to identify the novel genes/transcripts expressed due to elevated temperature conditions. The different expression analysis showed 1137 transcripts which were up-regulated during heat temperature stress when compared to control conditions. Analysis of relative gene expression showed phytepsin, ferredoxin-dependent glutamate synthase, and stress protein DDR-48 threefold increased expression during heat stress. These findings reveal novel targets for subsequent research on the genomics genetic manipulation and molecular mechanism of elevated stress tolerance in sugarcane.

Physiological characteristics and metabolomics reveal the tolerance mechanism to low nitrogen in Glycine soja leaves

To explore the regulatory mechanisms involved in the adaption to nitrogen (N) deficiency of wild soybean, the ion balance, photosynthetic characteristics, metabolic and transcriptional changes in leaves of common and low N (LN)-tolerant wild soybean seedlings under LN stress were determined. The LN-tolerant wild soybean seedlings showed a stronger ability to maintain photosynthesis and nutrient balance than common wild soybean. A total of 52 differentially accumulated metabolites, mainly related to carbon and N metabolism, were identified between the control and the LN treatment group. In general, tricarboxylic acid (TCA) cycle, shikimic acid pathway, synthetase/glutamate synthase (GS/GOGAT) cycle and accumulation of most organic acids were enhanced in LN-tolerant wild soybean, while reduced in common wild soybean under LN stress compared with their respective control group. Moreover, glycolysis, sugar and polyol and fatty acid metabolism increased in both wild soybean genotypes, and increased more in LN-tolerant wild soybean. A total of 3381 differentially expressed genes (DEGs) were identified in leaves of both wild soybean genotypes and the expressed level of DEGs associated with sugars, polyols, fatty acids and energy metabolism was significantly higher in LN-tolerant wild soybean than in common wild soybean, consistent with changes in metabolite level. Our results suggest new ideas for the study of LN tolerance of wild soybean and provide a theoretical basis for development and utilization of wild soybean resources.

Transcriptome Analysis Reveals Differences in Key Genes and Pathways Regulating Carbon and Nitrogen Metabolism in Cotton Genotypes under N Starvation and Resupply

Nitrogen (N) is the most important limiting factor for cotton production worldwide. Genotype-dependent ability to cope with N shortage has been only partially explored in cotton, and in this context, the comparison of molecular responses of cotton genotypes with different nitrogen use efficiency (NUE) is of particular interest to dissect the key molecular mechanisms underlying NUE. In this study, we employed Illumina RNA-Sequencing to determine the genotypic difference in transcriptome profile using two cotton genotypes differing in NUE (CCRI-69, N-efficient, and XLZ-30, N-inefficient) under N starvation and resupply treatments. The results showed that a large genetic variation existed in differentially expressed genes (DEGs) related to amino acid, carbon, and nitrogen metabolism between CCRI-69 and XLZ-30. Further analysis of metabolic changes in cotton genotypes under N resupply showed that nitrogen metabolism and aromatic amino acid metabolism pathways were mainly enriched in CCRI-69 by regulating carbon metabolism pathways such as starch and sucrose metabolism, glycolysis/gluconeogenesis, and pentose phosphate pathway. Additionally, we performed an expression network analysis of genes related to amino acid, carbon, and nitrogen metabolism. In total, 75 and 33 genes were identified as hub genes in shoots and roots of cotton genotypes, respectively. In summary, the identified hub genes may provide new insights into coordinating carbon and nitrogen metabolism and improving NUE in cotton.

Expression of a cucumber alanine aminotransferase2 gene improves nitrogen use efficiency in transgenic rice

Background

Rice can absorb less than 40% of applied nitrogen fertilizer, whereas the unabsorbed nitrogen fertilizer may cause environmental problems, such as algal blooms in freshwater and increased production of nitrous oxide, a greenhouse gas which is 300 times more potent than carbon dioxide. Development of nitrogen use efficient (NUE) rice is essential for more environmentally friendly rice production. Recently, NUE rice has been developed by root-specific expression of alanine aminotransferase (AlaAT) gene from barley, a monocot plant. Therefore, we tested the efficacy of AlaAT gene from cucumber in transgenic rice, aiming to provide evidence for the conservation of AlaAT gene function in monocot and dicot.

Results

AlaAT gene from cucumber (CsAlaAT2) has been successfully cloned and constructed on pCAMBIA1300 plant expression vectors under the control of tissue-specific promoter OsAnt1. Agrobacterium tumefaciens-mediated transformation of Indonesian rice cv. Fatmawati using this construct produced 14 transgenic events. Pre-screening of T1 seedlings grown in the agar medium containing low nitrogen concentration identified selected events that were superior in the root dry weight. Southern hybridization confirmed the integration of T-DNA in the selected event genomes, each of them carried 1, 2, or 3 T-DNA insertions. Efficacy assay of three lead events in the greenhouse showed that in general transgenic events had increased biomass, tiller number, nitrogen content, and grain yield compared to WT. One event, i.e., FAM13, showed an increase in yield as much as 27.9% and higher plant biomass as much as 27.4% compared to WT under the low nitrogen condition. The lead events also showed higher absorption NUE, agronomical NUE, and grain NUE as compared to WT under the low nitrogen condition.

Conclusions

The results of this study showed that root-specific expression of cucumber alanine aminotransferase2 gene improved nitrogen use efficiency in transgenic rice, which indicate the conservation of function of this gene in monocot and dicot.

Photosynthesis and yield response to elevated CO2, C4 plant foxtail millet behaves similarly to C3 species

Foxtail millet (Setaria italica) is a nutrient-rich food source traditionally grown in arid and semi-arid areas, as it is well adapted to drought climate. Yet there is limited information as how the crop responses to the changing climate. In order to investigate the response of foxtail millet to elevated [CO₂] and the underlying mechanism, the crop was grown at ambient [CO₂] (400 μmol mol^-1) and elevated [CO₂] (600 μmol mol^-1) in an open-top chamber (OTC) experimental facility in North China. The changes in leaf photosynthesis, chlorophyll fluorescence, biomass, yield and global gene expression in response to elevated [CO₂] were determined. Despite foxtail millet being a C₄ photosynthetic crop, photosynthetic rates (P_N) and intrinsic water-use efficiency (WUEi), were increased under elevated [CO₂]. Similarly, grain yield and above-ground biomass also significantly increased (P <  0.05) for the two years of experimentation under elevated [CO₂]. Increases in seeds and tiller number, spike and stem weight were the main contributors to the increased grain yield and biomass. Using transcriptomic analyses, this study further identified some genes which play a role in cell wall reinforcement, shoot initiation, stomatal conductance, carbon fixation, glycolysis / gluconeogenesis responsive to elevated [CO₂]. Changes in these genes reduced plant height, increased stem diameters, and promote CO₂ fixation. Higher photosynthetic rates at elevated [CO₂] demonstrated that foxtail millet was not photosynthetically saturated at elevated [CO₂] and its photosynthesis response to elevated [CO₂] were analogous to C₃ plants.

Genome-wide expression profiling of leaves and roots of watermelon in response to low nitrogen

BACKGROUND

Nitrogen (N) is a key macronutrient required for plant growth and development. In this study, watermelon plants were grown under hydroponic conditions at 0.2 mM N, 4.5 mM N, and 9 mM N for 14 days.

RESULTS

Dry weight and photosynthetic assimilation at low N (0.2 mM) was reduced by 29 and 74% compared with high N (9 mM). The photochemical activity (Fv/Fm) was also reduced from 0.78 at high N to 0.71 at low N. The N concentration in the leaf, stem, and root of watermelon under low N conditions was reduced by 68, 104, and 108%, respectively compared with 9 mM N treatment after 14 days of N treatment. In the leaf tissues of watermelon grown under low N conditions, 9598 genes were differentially expressed, out of which 4533 genes (47.22%) were up-regulated whereas, 5065 genes (52.78%) were down-regulated compared with high N. Similarly in the root tissues, 3956 genes were differentially expressed, out of which 1605 genes were up-regulated (40.57%) and 2351 genes were down-regulated (59.43%), compared with high N. Our results suggest that leaf tissues are more sensitive to N deficiency compared with root tissues. The gene ontology (GO) analysis showed that the availability of N significantly affected 19 biological processes, 8 cell component metabolic pathways, and 3 molecular functions in the leaves; and 13 biological processes, 12 molecular functions, and 5 cell component metabolic pathways in the roots of watermelon. The low affinity nitrate transporters, high affinity nitrate transporters, ammonium transporters, genes related with nitrogen assimilation, and chlorophyll and photosynthesis were expressed differentially in response to low N. Three nitrate transporters (Cla010066, Cla009721, Cla012765) substantially responded to low nitrate supply in the root and leaf tissues. Additionally, a large number of transcription factors (1365) were involved in adaptation to low N availability. The major transcription factor families identified in this study includes MYB, AP2-EREBP, bHLH, C2H2 and NAC.

CONCLUSION

Candidate genes identified in this study for nitrate uptake and transport can be targeted and utilized for further studies in watermelon breeding and improvement programs to improve N uptake and utilization efficiency.

De novo characterization of Phenacoccus solenopsis transcriptome and analysis of gene expression profiling during development and hormone biosynthesis

The cotton mealybug Phenacoccus solenopsis is a devastating pest of cotton causing tremendous loss in the yield of crops each year. Widespread physiological and biological studies on P. solenopsis have been carried out, but the lack of genetic information has constrained our understanding of the molecular mechanisms behind its growth and development. To understand and characterize the different developmental stages, RNA-Seq platform was used to execute de-novo transcriptome assembly and differential gene expression profiling for the eggs, first, second, third instar and adult female stages. About 182.67 million reads were assembled into 93,781 unigenes with an average length of 871.4 bp and an N50 length of 1899 bp. These unigenes sequences were annotated and classified by performing NCBI non-redundant (Nr) database, Kyoto Encyclopedia of Genes and Genomes (KEGG) and Clusters of Orthologous Groups (COG), Gene ontology (GO), the Swiss-Prot protein database (Swiss-Prot), and nearest related organism Acyrthosiphon pisum (pea aphid) database. To get more information regarding the process of metamorphosis, we performed a pairwise comparison of four developmental stages and obtained 29,415 differentially expressed genes. Some of the differentially expressed genes were associated with functional protein synthesis, anti-microbial protection, development and hormone biosynthesis. Functional pathway enrichment analysis of differentially expressed genes showed the positive correlation with specific physiological activities of each stage, and these results were confirmed by qRT-PCR experiments. This study gives a valuable genomics resource of P. solenopsis covering all its developmental stages and will promote future studies on biological processes at the molecular level.

Transcriptome and Co-Expression Network Analyses Identify Key Genes Regulating Nitrogen Use Efficiency in Brassica juncea L

Nitrate is the main source of inorganic nitrogen for plants, which also act as signaling molecule. Present study was aimed to understand nitrate regulatory mechanism in Brassica juncea cultivars, with contrasting nitrogen-use-efficiency (NUE) viz. Pusa Bold (PB, high-NUE) and Pusa Jai Kisan (PJK, low-NUE), employing RNA-seq approach. A total of 4031, 3874 and 3667 genes in PB and 2982, 2481 and 2843 genes in PJK were differentially expressed in response to early, low (0.25 mM KNO₃), medium (2 mM KNO₃) and high (4 mM KNO₃) nitrate treatments, respectively, as compared to control (0 mM KNO₃). Genes of N-uptake (NRT1.1, NRT1.8, and NRT2.1), assimilation (NR1, NR2, NiR, GS1.3, and Fd-GOGAT) and remobilization (GDH2, ASN2–3 and ALaT) were highly-upregulated in PB than in PJK in response to early nitrate treatments. We have also identified transcription factors and protein kinases that were rapidly induced in response to nitrate, suggesting their involvement in nitrate-mediated signaling. Co-expression network analysis revealed four nitrate specific modules in PB, enriched with GO terms like, “Phenylpropanoid pathway”, “Nitrogen compound metabolic process” and “Carbohydrate metabolism”. The network analysis also identified HUB transcription factors like mTERF, FHA, Orphan, bZip and FAR1, which may be the key regulators of nitrate-mediated response in B. juncea.

Transcriptome Analysis of Two Rice Varieties Contrasting for Nitrogen Use Efficiency under Chronic N Starvation Reveals Differences in Chloroplast and Starch Metabolism-Related Genes

The nitrogen use efficiency (NUE) of crop plants is limited and enhancing it in rice, a major cereal crop, would be beneficial for farmers and the environment alike. Here we report the genome-wide transcriptome analysis of two rice genotypes, IR 64 (IR64) and Nagina 22 (N22) under optimal (N+) and chronic starvation (N−) of nitrogen (N) from 15-day-old root and shoot tissues. The two genotypes were found to be contrasting in their response to N−; IR64 root architecture and root dry weight remained almost equivalent to that under N+ conditions, while N22 showed high foraging ability but a substantial reduction in biomass under N−. Similarly, the photosynthetic pigments showed a drastic reduction in N22 under low N, while IR64 was more resilient. Nitrate reductase showed significantly low specific activity under N− in both genotypes. Glutamate synthase (GOGAT) and citrate synthase CS activity were highly reduced in N22 but not in IR64. Transcriptome analysis of these genotypes revealed nearly double the number of genes to be differentially expressed (DEGs) in roots (1016) compared to shoots (571). The response of the two genotypes to N starvation was distinctly different reflecting their morphological/biochemical response with just two and eight common DEGs in the root and shoot tissues. There were a total of 385 nitrogen-responsive DEGs (106 in shoots and 279 in roots) between the two genotypes. Fifty-two of the 89 DEGs identified as specific to N22 root tissues were also found to be differentially expressed between the two genotypes under N−. Most of these DEGs belonged to starch and chloroplast metabolism, followed by membrane and signaling proteins. Physical mapping of DEGs revealed 95 DEGs in roots and 76 in shoots to be present in quantitative trait loci (QTL) known for NUE.

Transcriptome profiling and cataloging differential gene expression in floral buds of fertile and sterile lines of cotton (Gossypium hirsutum L.)

Cytoplasmic Male Sterility is maternally inherited trait in plants, characterized by failure to produce functional pollen during anther development. Anther development is modulated through the interaction of nuclear and mitochondrial genes. In the present study, differential gene expression of floral buds at the sporogenous stage (SS) and microsporocyte stage (MS) between CGMS and its fertile maintainer line of cotton plants was studied. A total of 320 significantly differentially expressed genes, including 20 down-regulated and 37 up-regulated in CGMS comparing with its maintainer line at the SS stage, as well as and 89 down-regulated and 4 up-regulated in CGMS compared to the fertile line at MS stage. Comparing the two stages in the same line, there were 6 down-regulated differentially expressed genes only induced in CGMS and 9 up-regulated differentially expressed gene only induced in its maintainer. GO analysis revealed essential genes responsible for pollen development, and cytoskeleton category show differential expression between the fertile and CGMS lines. Validation studies by qRT-PCR shows concordance with RNA-seq result. A set of novel SSRs identified in this study can be used in evaluating genetic relationships among cultivars, QTL mapping, and marker-assisted breeding. We reported aberrant expression of genes related to pollen exine formation, and synthesis of pectin lyase, myosine heavy chain, tubulin, actin-beta, heat shock protein and myeloblastosis (MYB) protein as targets for CMS in cotton. The results of this study contribute to basic information for future screening of genes and identification of molecular portraits responsible for CMS as well as to elucidate molecular mechanisms that lead to CMS in cotton.

Stability evaluation of reference genes for gene expression analysis by RT-qPCR in soybean under different conditions

Real-time quantitative reverse transcription PCR is a sensitive and widely used technique to quantify gene expression. To achieve a reliable result, appropriate reference genes are highly required for normalization of transcripts in different samples. In this study, 9 previously published reference genes (60S, Fbox, ELF1A, ELF1B, ACT11, TUA5, UBC4, G6PD, CYP2) of soybean [Glycine max (L.) Merr.] were selected. The expression stability of the 9 genes was evaluated under conditions of biotic stress caused by infection with soybean mosaic virus, nitrogen stress, across different cultivars and developmental stages. ΔCt and geNorm algorithms were used to evaluate and rank the expression stability of the 9 reference genes. Results obtained from two algorithms showed high consistency. Moreover, results of pairwise variation showed that two reference genes were sufficient to normalize the expression levels of target genes under each experimental setting. For virus infection, ELF1A and ELF1B were the most stable reference genes for accurate normalization. For different developmental stages, Fbox and G6PD had the highest expression stability between two soybean cultivars (Tanlong No. 1 and Tanlong No. 2). ELF1B and ACT11 were identified as the most stably expressed reference genes both under nitrogen stress and among different cultivars. The results showed that none of the candidate reference genes were uniformly expressed at different conditions, and selecting appropriate reference genes was pivotal for gene expression studies with particular condition and tissue. The most stable combination of genes identified in this study will help to achieve more accurate and reliable results in a wide variety of samples in soybean.

Nitrate simultaneously enhances lipid and protein accumulation in developing yellow lupin cotyledons cultured in vitro, but not under field conditions

The research was conducted on yellow lupin (Lupinus luteus L.) mature seeds, developing cotyledons, developing pods, and seedlings. The main storage compound in yellow lupin seeds is protein, whose content may reach up to 45%. Oil content in seeds of yellow lupin is about 6%. In such protein-storing seeds there is a strong negative relationship between accumulation of storage lipid and protein. An increase in protein content causes a decrease in lipid level, and vice versa. However, simultaneous increase in lipid and protein content is possible in developing lupin cotyledons (the main storage organs of lupin seeds) cultured in vitro. Such an effect was obtained by feeding the cotyledons with nitrate (35mM). The same positive relationship in storage lipid and protein accumulation was also obtained in developing lupin pods fed with nitrate (35mM), detached from the mother plant, and maintained under quasi in vitro conditions. Fertilization of lupin plants with nitrate under field conditions (40 or 80kgNha^-1 applied before sowing, at the nodulation stage or at the flowering and pod formation stage) did not cause significant changes in lipid and protein contents in mature seeds. Experiments performed on lupin seedlings cultivated hydroponically showed that nitrate added to the medium was accumulated mainly in roots, and at a remarkably lower level in shoots. We hypothesize that the lack of stimulatory effect of nitrate on storage lipid and protein accumulation in seeds under field conditions is due to inefficient transport of nitrate from the root to developing pods in lupin plants. This causes that the level of nitrate inside the developing lupin seeds is not elevated under field conditions.

Exploring the heat-responsive chaperones and microsatellite markers associated with terminal heat stress tolerance in developing wheat

Global warming is a major threat for agriculture and food security, and in many cases the negative impacts are already apparent. Wheat is one of the most important staple food crops and is highly sensitive to the heat stress (HS) during reproductive and grain-filling stages. Here, whole transcriptome analysis of thermotolerant wheat cv. HD2985 was carried out at the post-anthesis stage under control (22 ± 3 °C) and HS-treated (42 °C, 2 h) conditions using Illumina Hiseq and Roche GS-FLX 454 platforms. We assembled ~24 million (control) and ~23 million (HS-treated) high-quality trimmed reads using different assemblers with optimal parameters. De novo assembly yielded 52,567 (control) and 59,658 (HS-treated) unigenes. We observed 785 transcripts to be upregulated and 431 transcripts to be downregulated under HS; 78 transcripts showed >10-fold upregulation such as HSPs, metabolic pathway-related genes, etc. Maximum number of upregulated genes was observed to be associated with processes such as HS-response, protein-folding, oxidation-reduction and photosynthesis. We identified 2008 and 2483 simple sequence repeats (SSRs) markers from control and HS-treated samples; 243 SSRs were observed to be overlying on stress-associated genes. Polymorphic study validated four SSRs to be heat-responsive in nature. Expression analysis of identified differentially expressed transcripts (DETs) showed very high fold increase in the expression of catalytic chaperones (HSP26, HSP17, and Rca) in contrasting wheat cvs. HD2985 and HD2329 under HS. We observed positive correlation between RNA-seq and qRT-PCR expression data. The present study culminated in greater understanding of the heat-response of tolerant genotype and has provided good candidate genes for the marker development and screening of wheat germplasm for thermotolerance.

De novo transcriptome analysis and glucosinolate profiling in watercress (Nasturtium officinale R. Br.)

Background

Watercress (Nasturtium officinale R. Br.) is an aquatic herb species that is a rich source of secondary metabolites such as glucosinolates. Among these glucosinolates, watercress contains high amounts of gluconasturtiin (2-phenethyl glucosinolate) and its hydrolysis product, 2-phennethyl isothiocyanate, which plays a role in suppressing tumor growth. However, the use of N. officinale as a source of herbal medicines is currently limited due to insufficient genomic and physiological information.

Results

To acquire precise information on glucosinolate biosynthesis in N. officinale, we performed a comprehensive analysis of the transcriptome and metabolome of different organs of N. officinale. Transcriptome analysis of N. officinale seedlings yielded 69,570,892 raw reads. These reads were assembled into 69,635 transcripts, 64,876 of which were annotated to transcripts in public databases. On the basis of the functional annotation of N. officinale, we identified 33 candidate genes encoding enzymes related to glucosinolate biosynthetic pathways and analyzed the expression of these genes in the leaves, stems, roots, flowers, and seeds of N. officinale. The expression of NoMYB28 and NoMYB29, the main regulators of aliphatic glucosinolate biosynthesis, was highest in the stems, whereas the key regulators of indolic glucosinolate biosynthesis, such as NoDof1.1, NoMYB34, NoMYB51, and NoMYB122, were strongly expressed in the roots. Most glucosinolate biosynthetic genes were highly expressed in the flowers. HPLC analysis enabled us to detect eight glucosinolates in the different organs of N. officinale. Among these glucosinolates, the level of gluconasturtiin was considerably higher than any other glucosinolate in individual organs, and the amount of total glucosinolates was highest in the flower.

Conclusions

This study has enhanced our understanding of functional genomics of N. officinale, including the glucosinolate biosynthetic pathways of this plant. Ultimately, our data will be helpful for further research on watercress bio-engineering and better strategies for exploiting its anti-carcinogenic properties.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3792-5) contains supplementary material, which is available to authorized users.

RNA-Seq analysis of nodule development at five different developmental stages of soybean (Glycine max) inoculated with Bradyrhizobium japonicum strain 113-2

Nodule development directly affects nitrogen fixation efficiency during soybean growth. Although abundant genome-based information related to nodule development has been released and some studies have reported the molecular mechanisms that regulate nodule development, information on the way nodule genes operate in nodule development at different developmental stages of soybean is limited. In this report, notably different nodulation phenotypes in soybean roots inoculated with Bradyrhizobium japonicum strain 113-2 at five developmental stages (branching stage, flowering stage, fruiting stage, pod stage and harvest stage) were shown, and the expression of nodule genes at these five stages was assessed quantitatively using RNA-Seq. Ten comparisons were made between these developmental periods, and their differentially expressed genes were analysed. Some important genes were identified, primarily encoding symbiotic nitrogen fixation-related proteins, cysteine proteases, cystatins and cysteine-rich proteins, as well as proteins involving plant-pathogen interactions. There were no significant shifts in the distribution of most GO functional annotation terms and KEGG pathway enrichment terms between these five development stages. A cystatin Glyma18g12240 was firstly identified from our RNA-seq, and was likely to promote nodulation and delay nodule senescence. This study provides molecular material for further investigations into the mechanisms of nitrogen fixation at different soybean developmental stages.

MicroRNA and Transcription Factor: Key Players in Plant Regulatory Network

Recent achievements in plant microRNA (miRNA), a large class of small and non-coding RNAs, are very exciting. A wide array of techniques involving forward genetic, molecular cloning, bioinformatic analysis, and the latest technology, deep sequencing have greatly advanced miRNA discovery. A tiny miRNA sequence has the ability to target single/multiple mRNA targets. Most of the miRNA targets are transcription factors (TFs) which have paramount importance in regulating the plant growth and development. Various families of TFs, which have regulated a range of regulatory networks, may assist plants to grow under normal and stress environmental conditions. This present review focuses on the regulatory relationships between miRNAs and different families of TFs like; NF-Y, MYB, AP2, TCP, WRKY, NAC, GRF, and SPL. For instance NF-Y play important role during drought tolerance and flower development, MYB are involved in signal transduction and biosynthesis of secondary metabolites, AP2 regulate the floral development and nodule formation, TCP direct leaf development and growth hormones signaling. WRKY have known roles in multiple stress tolerances, NAC regulate lateral root formation, GRF are involved in root growth, flower, and seed development, and SPL regulate plant transition from juvenile to adult. We also studied the relation between miRNAs and TFs by consolidating the research findings from different plant species which will help plant scientists in understanding the mechanism of action and interaction between these regulators in the plant growth and development under normal and stress environmental conditions.

Physiological, genomic and transcriptional diversity in responses to boron deficiency in rapeseed genotypes

Allotetraploid rapeseed (Brassica napus L. A_nA_nC_nC_n, 2n=4x=38) is highly susceptible to boron (B) deficiency, a widespread limiting factor that causes severe losses in seed yield. The genetic variation in the sensitivity to B deficiency found in rapeseed genotypes emphasizes the complex response architecture. In this research, a B-inefficient genotype, 'Westar 10' ('W10'), responded to B deficiencies during vegetative and reproductive development with an over-accumulation of reactive oxygen species, severe lipid peroxidation, evident plasmolysis, abnormal floral organogenesis, and widespread sterility compared to a B-efficient genotype, 'Qingyou 10' ('QY10'). Whole-genome re-sequencing (WGS) of 'QY10' and 'W10' revealed a total of 1 605 747 single nucleotide polymorphisms and 218 755 insertions/deletions unevenly distributed across the allotetraploid rapeseed genome (~1130Mb). Digital gene expression (DGE) profiling identified more genes related to B transporters, antioxidant enzymes, and the maintenance of cell walls and membranes with higher transcript levels in the roots of 'QY10' than in 'W10' under B deficiency. Furthermore, based on WGS and bulked segregant analysis of the doubled haploid (DH) line pools derived from 'QY10' and 'W10', two significant quantitative trait loci (QTLs) for B efficiency were characterized on chromosome C2, and DGE-assisted QTL-seq analyses then identified a nodulin 26-like intrinsic protein gene and an ATP-binding cassette (ABC) transporter gene as the corresponding candidates regulating B efficiency. This research facilitates a more comprehensive understanding of the differential physiological and transcriptional responses to B deficiency and abundant genetic diversity in rapeseed genotypes, and the DGE-assisted QTL-seq analyses provide novel insights regarding the rapid dissection of quantitative trait genes in plant species with complex genomes.

De novo transcriptome sequencing and comparative analysis to discover genes related to floral development in Cymbidium faberi Rolfe

Cymbidium faberi is a traditional orchid flower in China that is highly appreciated for its fragrant aroma from its zygomorphic flowers. One bottleneck of the commercial production of C. faberi is the long vegetative growth phase of the orchid and the difficulty of the regulation of its flowering time. Moreover, its flower size, shape and color are often targeting traits for orchid breeders. Understanding the molecular mechanisms of floral development in C. faberi will ultimately benefit the genetic improvement of this orchid plant. The goal of this study is to identify potential genes and regulatory networks related to the floral development in C. faberi by using transcriptome sequencing, de novo assembly and computational analyses. The vegetative and flower buds of C. faberi were sampled for such comparisons. The RNA-seq yielded about 189,300 contigs that were assembled into 172,959 unigenes. Furthermore, a total of 13,484 differentially expressed unigenes (DEGs) were identified between the vegetative and flower buds. There were 7683 down-regulated and 5801 up-regulated DEGs in the flower buds compared to those in the vegetative buds, among which 3430 and 6556 DEGs were specifically enriched in the flower or vegetative buds, respectively. A total of 173 DEGs orthologous to known genes associated with the floral organ development, floral symmetry and flowering time were identified, including 12 TCP transcription factors, 34 MADS-box genes and 28 flowering time related genes. Furthermore, expression levels of ten genes potentially involved in floral development and flowering time were verified by quantitative real-time PCR. The identified DEGs will facilitate the functional genetic studies for further understanding the flower development of C. faberi.

Transcriptomics-assisted quantitative trait locus fine mapping for the rapid identification of a nodulin 26-like intrinsic protein gene regulating boron efficiency in allotetraploid rapeseed

Allotetraploid rapeseed (Brassica napus L., An An Cn Cn , 2n = 4x = 38) is extraordinarily susceptible to boron (B) deficiency, a ubiquitous problem causing severe losses in seed yield. The breeding of B-efficient rapeseed germ plasm is a cost-effective and environmentally friendly strategy for the agricultural industry; however, genes regulating B efficiency in allotetraploid rapeseed have not yet been isolated. In this research, quantitative trait locus (QTL) fine mapping and digital gene expression (DGE) profiling were combined to identify the candidate genes underlying the major-effect QTL qBEC-A3a, which regulates B efficiency. Comparative phenotype analyses of the near-isogenic lines (NILs) indicated that qBEC-A3a plays a significant role in improving B efficiency under B deficiency. Exploiting QTL fine mapping and DGE analyses revealed a nodulin 26-like intrinsic protein (NIP) gene, which encodes a likely boric acid channel. The gene co-expression network for putative B transporters also highlighted its central role in the efficiency of B uptake. An integration of whole-genome re-sequencing (WGS) with bulked segregant analysis (BSA) authenticated the emerging availability of QTL-seq for the QTL analyses in allotetraploid rapeseed. Transcriptomics-assisted QTL mapping and comparative genomics provided novel insights into the rapid identification of quantitative trait genes (QTGs) in plant species with complex genomes.

Identification and Comparative Analysis of CBS Domain-Containing Proteins in Soybean (Glycine max) and the Primary Function of GmCBS21 in Enhanced Tolerance to Low Nitrogen Stress

Nitrogen is an important macronutrient required for plant growth, and is a limiting factor for crop productivity. Improving the nitrogen use efficiency (NUE) is therefore crucial. At present, the NUE mechanism is unclear and information on the genes associated with NUE in soybeans is lacking. cystathionine beta synthase (CBS) domain-containing proteins (CDCPs) may be implicated in abiotic stress tolerance in plants. We identified and classified a CBS domain-containing protein superfamily in soybean. A candidate gene for NUE, GmCBS21, was identified. GmCBS21 gene characteristics, the temporal expression pattern of the GmCBS21 gene, and the phenotype of GmCBS21 overexpression in transgenic Arabidopsis thaliana under low nitrogen stress were analyzed. The phenotypes suggested that the transgenic Arabidopsis thaliana seedlings performed better under the nitrogen-deficient condition. GmCBS21-overexpressing transgenic plants exhibit higher low nitrogen stress tolerance than WT plants, and this suggests its role in low nitrogen stress tolerance in plants. We conclude that GmCBS21 may serve as an excellent candidate for breeding crops with enhanced NUE and better yield.

Harnessing Next Generation Sequencing in Climate Change: RNA-Seq Analysis of Heat Stress-Responsive Genes in Wheat (Triticum aestivum L.)

Wheat is a staple food worldwide and provides 40% of the calories in the diet. Climate change and global warming pose a threat to wheat production, however, and demand a deeper understanding of how heat stress might impact wheat production and wheat biology. However, it is difficult to identify novel heat stress associated genes when the genomic information is not available. Wheat has a very large and complex genome that is about 37 times the size of the rice genome. The present study sequenced the whole transcriptome of the wheat cv. HD2329 at the flowering stage, under control (22°±3°C) and heat stress (42°C, 2 h) conditions using Illumina HiSeq and Roche GS-FLX 454 platforms. We assembled more than 26.3 and 25.6 million high-quality reads from the control and HS-treated tissues transcriptome sequences respectively. About 76,556 (control) and 54,033 (HS-treated) contigs were assembled and annotated de novo using different assemblers and a total of 21,529 unigenes were obtained. Gene expression profile showed significant differential expression of 1525 transcripts under heat stress, of which 27 transcripts showed very high (>10) fold upregulation. Cellular processes such as metabolic processes, protein phosphorylation, oxidations-reductions, among others were highly influenced by heat stress. In summary, these observations significantly enrich the transcript dataset of wheat available on public domain and show a de novo approach to discover the heat-responsive transcripts of wheat, which can accelerate the progress of wheat stress-genomics as well as the course of wheat breeding programs in the era of climate change.

Development of chloroplast genomic resources for Cynara

In this study, new chloroplast (cp) resources were developed for the genus Cynara, using whole cp genomes from 20 genotypes, by means of high-throughput sequencing technologies. Our target species included seven globe artichokes, two cultivated cardoons, eight wild artichokes, and three other wild Cynara species (C. baetica, C. cornigera and C. syriaca). One complete cp genome was isolated using short reads from a whole-genome sequencing project, while the others were obtained by means of long-range PCR, for which primer pairs are provided here. A de novo assembly strategy combined with a reference-based assembly allowed us to reconstruct each cp genome. Comparative analyses among the newly sequenced genotypes and two additional Cynara cp genomes ('Brindisino' artichoke and C. humilis) retrieved from public databases revealed 126 parsimony informative characters and 258 singletons in Cynara, for a total of 384 variable characters. Thirty-nine SSR loci and 34 other INDEL events were detected. After data analysis, 37 primer pairs for SSR amplification were designed, and these molecular markers were subsequently validated in our Cynara genotypes. Phylogenetic analysis based on all cp variable characters provided the best resolution when compared to what was observed using only parsimony informative characters, or only short 'variable' cp regions. The evaluation of the molecular resources obtained from this study led us to support the 'super-barcode' theory and consider the total cp sequence of Cynara as a reliable and valuable molecular marker for exploring species diversity and examining variation below the species level.

Identification of MicroRNAs in Response to Different Day Lengths in Soybean Using High-Throughput Sequencing and qRT-PCR

MicroRNAs (miRNAs) are short, non-coding single-strand RNA molecules that play important roles in plant growth, development and stress responses. Flowering time affects the seed yield and quality of soybean. However, the miRNAs involved in the regulation of flowering time in soybean have not been reported until recently. Here, high-throughput sequencing and qRT-PCR were used to identify miRNAs involved in soybean photoperiodic pathways. The first trifoliate leaves of soybean that receive the signal of light treatment were used to construct six libraries (0, 8, and 16 h under short-day (SD) treatment and 0, 8, and 16 h under long-day (LD) treatment). The libraries were sequenced using Illumina Solexa. A total of 318 known plant miRNAs belonging to 163 miRNA families and 81 novel predicted miRNAs were identified. Among these, 23 miRNAs at 0 h, 65 miRNAs at 8 h and 83 miRNAs at 16 h, including six novel predicted miRNAs at 8 h and six novel predicted miRNAs at 16 h, showed differences in abundance between LD and SD treatments. Furthermore, the results of GO and KEGG analyses indicated that most of the miRNA targets were transcription factors. Seven miRNAs at 0 h, 23 miRNAs (including four novel predicted miRNAs) at 8 h, 16 miRNAs (including one novel predicted miRNA) at 16 h and miRNA targets were selected for qRT-PCR analysis to assess the accuracy of the sequencing and target prediction. The results indicated that the expression patterns of the selected miRNAs and miRNA targets showed no differences between the qRT-PCR and sequencing results. In addition, 23 miRNAs at 0 h, 65 miRNAs at 8 h and 83 miRNAs at 16 h responded to day length changes in soybean, including six novel predicted miRNAs at 8 h and six novel predicted miRNAs at 16 h. These results provided an important molecular basis to understand the regulation of flowering time through photoperiodic pathways in soybean.

Comparative Transcriptome Analysis between the Cytoplasmic Male Sterile Line NJCMS1A and Its Maintainer NJCMS1B in Soybean (Glycine max (L.) Merr.)

BACKGROUND

The utilization of soybean heterosis is probably one of the potential approaches in future yield breakthrough as was the situation in rice breeding in China. Cytoplasmic male sterility (CMS) plays an important role in the production of hybrid seeds. However, the molecular mechanism of CMS in soybean remains unclear.

RESULTS

The comparative transcriptome analysis between cytoplasmic male sterile line NJCMS1A and its near-isogenic maintainer NJCMS1B in soybean was conducted using Illumina sequencing technology. A total of 88,643 transcripts were produced in Illumina sequencing. Then 56,044 genes were obtained matching soybean reference genome. Three hundred and sixty five differentially expressed genes (DEGs) between NJCMS1A and NJCMS1B were screened by threshold, among which, 339 down-regulated and 26 up-regulated in NJCMS1A compared to in NJCMS1B. Gene Ontology (GO) annotation showed that 242 DEGs were annotated to 19 functional categories. Clusters of Orthologous Groups of proteins (COG) annotation showed that 265 DEGs were classified into 19 categories. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that 46 DEGs were assigned to 33 metabolic pathways. According to functional and metabolic pathway analysis combined with reported literatures, the relations between some key DEGs and the male sterility of NJCMS1A were discussed. qRT-PCR analysis validated that the gene expression pattern in RNA-Seq was reliable. Finally, enzyme activity assay showed that energy supply was decreased in NJCMS1A compared to in NJCMS1B.

CONCLUSIONS

We concluded that the male sterility of NJCMS1A might be related to the disturbed functions and metabolism pathways of some key DEGs, such as DEGs involved in carbohydrate and energy metabolism, transcription factors, regulation of pollen development, elimination of reactive oxygen species (ROS), cellular signal transduction, and programmed cell death (PCD) etc. Future research will focus on cloning and transgenic function validation of possible candidate genes associated with soybean CMS.

Genome-wide survey of the soybean GATA transcription factor gene family and expression analysis under low nitrogen stress

GATA transcription factors are transcriptional regulatory proteins that contain a characteristic type-IV zinc finger DNA-binding domain and recognize the conserved GATA motif in the promoter sequence of target genes. Previous studies demonstrated that plant GATA factors possess critical functions in developmental control and responses to the environment. To date, the GATA factors in soybean (Glycine max) have yet to be characterized. Thus, this study identified 64 putative GATA factors from the entire soybean genomic sequence. The chromosomal distributions, gene structures, duplication patterns, phylogenetic tree, tissue expression patterns, and response to low nitrogen stress of the 64 GATA factors in soybean were analyzed to further investigate the functions of these factors. Results indicated that segmental duplication predominantly contributed to the expansion of the GATA factor gene family in soybean. These GATA proteins were phylogenetically clustered into four distinct subfamilies, wherein their gene structure and motif compositions were considerably conserved. A comparative phylogenetic analysis of the GATA factor zinc finger domain sequences in soybean, Arabidopsis (Arabidopsis thaliana), and rice (Oryza sativa) revealed four major classes. The GATA factors in soybean exhibited expression diversity among different tissues; some of these factors showed tissue-specific expression patterns. Numerous GATA factors displayed upregulation or downregulation in soybean leaf in response to low nitrogen stress, and two GATA factors GATA44 and GATA58 were likely to be involved in the regulation of nitrogen metabolism in soybean. Overexpression of GmGATA44 complemented the reduced chlorophyll phenotype of the Arabidopsis ortholog AtGATA21 mutant, implying that GmGATA44 played an important role in modulating chlorophyll biosynthesis. Overall, our study provides useful information for the further analysis of the biological functions of GATA factors in soybean and other crops.