Integrated transcriptomics and metabolomics analyses provide insights into cold stress response in wheat

Physiological and transcriptomic analysis of tomato in response to sub-optimal temperature stress

Tomato (Solanum lycopersicum L.) is one of the most important economic crops in China. However, its quality and yield are susceptible to the adverse effects of low temperatures. In our study, two tomato cultivars, showing different tolerance to low temperatures, namely the cold-sensitive tomato cultivar (S708) and cold-tolerant tomato cultivar (T722), were grown at optimal (25/18°C) and sub-optimal (15/10°C) temperature conditions for 5 days. Our study aimed to explore the effect of sub-optimal temperature on fresh weight, chlorophyll content and chlorophyll fluorescence, soluble sugars and proline content of two tomato cultivars. Moreover, we employed RNA-Seq to analyze the transcriptomic response of tomato roots to sub-optimal temperature. The results revealed that S708 showed a more significant reduction in fresh weight, chlorophyll content, photochemical efficiency of PSII (YII), maximum quantum yield of PSII (Fv/Fm), photochemical quenching (qP) and electron transport rate (ETR) compared to T722 under the sub-optimal temperature condition. Notably, T722 maintained higher level of soluble sugars and proline in comparison to S708 uner sub-optimal temperature. RNA-seq data showed that up-regulated DEGs in both tomato cultivars were involved in "plant-pathogen interaction", "MAPK signaling pathway", "plant hormone signal transduction", and "phosphatidylinositol signaling system". Furthermore, "Amino sugar and nucleotide sugar metabolism" pathway was enriched only in T722. Moreover, under sub-optimal temperature, transcription factor genes and osmoregulation genes showed varying degrees of response in both tomato cultivars. Conclusion: In summary, our results offer detailed insights into the response characteristics of tomato to sub-optimal temperature, providing valuable references for the practical management of tomato crops under sub-optimal temperature condition.

Multi-Omics Approaches in Oil Palm Research: A Comprehensive Review of Metabolomics, Proteomics, and Transcriptomics Based on Low-Temperature Stress

Oil palm (Elaeis guineensis Jacq.) is a typical tropical oil crop with a temperature of 26-28 °C, providing approximately 35% of the total world's vegetable oil. Growth and productivity are significantly affected by low-temperature stress, resulting in inhibited growth and substantial yield losses. To comprehend the intricate molecular mechanisms underlying the response and acclimation of oil palm under low-temperature stress, multi-omics approaches, including metabolomics, proteomics, and transcriptomics, have emerged as powerful tools. This comprehensive review aims to provide an in-depth analysis of recent advancements in multi-omics studies on oil palm under low-temperature stress, including the key findings from omics-based research, highlighting changes in metabolite profiles, protein expression, and gene transcription, as well as including the potential of integrating multi-omics data to reveal novel insights into the molecular networks and regulatory pathways involved in the response to low-temperature stress. This review also emphasizes the challenges and prospects of multi-omics approaches in oil palm research, providing a roadmap for future investigations. Overall, a better understanding of the molecular basis of the response of oil palm to low-temperature stress will facilitate the development of effective breeding and biotechnological strategies to improve the crop's resilience and productivity in changing climate scenarios.

Jasmonate enhances cold acclimation in jojoba by promoting flavonol synthesis

Jojoba is an industrial oil crop planted in tropical arid areas, and its low-temperature sensitivity prevents its introduction into temperate areas. Studying the molecular mechanisms associated with cold acclimation in jojoba is advantageous for developing breeds with enhanced cold tolerance. In this study, metabolomic analysis revealed that various flavonols accumulate in jojoba during cold acclimation. Time-course transcriptomic analysis and weighted correlation network analysis (WGCNA) demonstrated that flavonol biosynthesis and jasmonates (JAs) signaling pathways played crucial roles in cold acclimation. Combining the biochemical and genetic analyses showed that ScMYB12 directly activated flavonol synthase gene (ScFLS). The interaction between ScMYB12 and transparent testa 8 (ScTT8) promoted the expression of ScFLS, but the negative regulator ScJAZ13 in the JA signaling pathway interacted with ScTT8 to attenuate the transcriptional activity of the ScTT8 and ScMYB12 complex, leading to the downregulation of ScFLS. Cold acclimation stimulated the production of JA in jojoba leaves, promoted the degradation of ScJAZ13, and activated the transcriptional activity of ScTT8 and ScMYB12 complexes, leading to the accumulation of flavonols. Our findings reveal the molecular mechanism of JA-mediated flavonol biosynthesis during cold acclimation in jojoba and highlight the JA pathway as a promising means for enhancing cold tolerance in breeding efforts.

Deacclimation and reacclimation processes in winter wheat: novel perspectives from time-series transcriptome analysis

Winter wheat achieves freezing tolerance (FT) through cold acclimation (CA) - a process which is induced by low positive temperatures in autumn. The increasing occurrences of temperature fluctuations in winter lead to deacclimation (DEA), causing premature loss of FT, and the cultivars capable of reacclimation (REA) are more likely to survive the subsequent cold spells. The genetic mechanisms of DEA and REA remain poorly understood, necessitating further research to bolster climate resilience in winter wheat. Here, we selected two winter wheat genotypes with contrasting levels of FT and conducted a ten-week-long experiment imitating low-temperature fluctuations after CA under controlled conditions. Crown and leaf tissue samples for RNA-sequencing were collected at CA, DEA, and REA time-points. It is the first transcriptomic study covering both short- and long-term responses to DEA and REA in winter wheat. The study provides novel knowledge regarding CA, DEA, and REA and discusses the gene expression patterns conferring FT under temperature fluctuations. The freezing-tolerant genotype "Lakaja DS" showed elevated photosynthetic activity in leaf tissue and upregulated cryoprotective protein-encoding genes in crowns after CA when compared to the freezing-susceptible "KWS Ferrum". "Lakaja DS" also expressed cold acclimation-associated transcripts at a significantly higher level after 1 week of DEA. Following REA, "Lakaja DS" continued to upregulate dehydrin-related genes in crowns and exhibited significantly higher expression of chitinase transcripts in leaves, when compared to "KWS Ferrum". The findings of this study shed light on the genetic mechanisms governing DEA and REA in winter wheat, thus addressing the gaps in knowledge regarding FT under low-temperature fluctuations. The identified genes should be further examined as potential molecular markers for breeding strategies focused on developing freezing-tolerant winter-type crops. Publicly available datasets generated in this study are valuable resources for further research into DEA and REA, contributing towards the enhancement of winter wheat under global climate change.

Metabolomics combined with physiology and transcriptomics reveal how Nicotiana tabacum leaves respond to cold stress

Low temperature-induced cold stress is a major threat to plant growth, development and distribution. Unraveling the responses of temperature-sensitive crops to cold stress and the mechanisms of cold acclimation are critical for food demand. In this study, combined physiological, transcriptomic, and metabolomic analyses were conducted on Nicotiana tabacum suffering short-term 4 °C cold stress. Our results showed that cold stress destroyed cellular membrane stability, decreased the chlorophyll (Chl) and carotenoid contents, and closed stomata, resulting in lipid peroxidation and photosynthesis restriction. Chl fluorescence measurements revealed that primary photochemistry, photoelectrochemical quenching and photosynthetic electron transport in Nicotiana tabacum leaves were seriously suppressed upon exposer to cold stress. Enzymatic and nonenzymatic antioxidants, including superoxide dismutase, catalase, peroxidase, reduced glutathione, proline, and soluble sugar, were all profoundly increased to trigger the cold acclimation defense against oxidative damage. A total of 178 metabolites and 16,204 genes were differentially expressed in cold-stressed Nicotiana tabacum leaves. MEturquoise and MEblue modules identified by WGCNA were highly correlated with physiological indices, and the corresponding hub genes were significantly enriched in pathways related to photosynthesis - antenna proteins and flavonoid biosynthesis. Untargeted metabolomic analysis identified specific metabolites, including sucrose, phenylalanine, glutamine, glutamate, and proline, that enhance plant cold acclimation. Combined transcriptomics and metabolomic analysis highlight the vital roles of carbohydrate and amino acid metabolism in enhancing the cold tolerance of Nicotiana tabacum. Our comprehensive investigation provides novel insights for efforts to alleviate low temperature-induced oxidative damage to Nicotiana tabacum plants and proposes a breeding target for cold stress-tolerant cultivars.

Synergistic regulation at physiological, transcriptional and metabolic levels in tomato plants subjected to a combination of salt and heat stress

With global warming and climate change, abiotic stresses often simultaneously occur. Combined salt and heat stress was a common phenomenon that was severe, particularly in arid/semi-arid lands. We aimed to reveal the systematic responsive mechanisms of tomato genotypes with different salt/heat susceptibilities to combined salt and heat stress. Morphological and physiological responses of salt-tolerant/sensitive and heat-tolerant/sensitive tomatoes at control, heat, salt and combined stress were investigated. Based on leaf Fv /Fm and H2 O2 content, samples from tolerant genotype at the four treatments for 36 h were taken for transcriptomics and metabolomics. We found that plant height, dry weight and net photosynthetic rate decreased while leaf Na+ concentration increased in all four genotypes under salt and combined stress than control. Changes in physiological indicators such as photosynthetic parameters and defence enzyme activities in tomato under combined stress were regulated by the expression of relevant genes and the accumulation of key metabolites. We screened five key pathways in tomato responding to a combination of salt and heat stress, such as oxidative phosphorylation (map00190). Synergistic regulation at morphological, physiological, transcriptional and metabolic levels in tomato plants was induced by combined stress. Heat stress was considered as a dominant stressor for tomato plants under the current combined stress. The oxidative phosphorylation pathway played a key role in tomato in response to combined stress, where tapped key genes (e.g. alternative oxidase, Aox1a) need further functional analysis. Our study will provide a valuable resource important for studying stress combination and improving tomato tolerance.

Advances and opportunities in unraveling cold-tolerance mechanisms in the world's primary staple food crops

Temperatures below or above optimal growth conditions are among the major stressors affecting productivity, end-use quality, and distribution of key staple crops including rice (Oryza sativa), wheat (Triticum aestivum), and maize (Zea mays L.). Among temperature stresses, cold stress induces cellular changes that cause oxidative stress and slowdown metabolism, limit growth, and ultimately reduce crop productivity. Perception of cold stress by plant cells leads to the activation of cold-responsive transcription factors and downstream genes, which ultimately impart cold tolerance. The response triggered in crops to cold stress includes gene expression/suppression, the accumulation of sugars upon chilling, and signaling molecules, among others. Much of the information on the effects of cold stress on perception, signal transduction, gene expression, and plant metabolism are available in the model plant Arabidopsis but somewhat lacking in major crops. Hence, a complete understanding of the molecular mechanisms by which staple crops respond to cold stress remain largely unknown. Here, we make an effort to elaborate on the molecular mechanisms employed in response to low-temperature stress. We summarize the effects of cold stress on the growth and development of these crops, the mechanism of cold perception, and the role of various sensors and transducers in cold signaling. We discuss the progress in cold tolerance research at the genome, transcriptome, proteome, and metabolome levels and highlight how these findings provide opportunities for designing cold-tolerant crops for the future.

Integrated Transcriptomics and Metabolomics Analyses Provide Insights into Qingke in Response to Cold Stress

The survival and productivity of qingke in high altitude (>4300 m, average yearly temperature <0 °C) of the Tibetan Plateau are significantly impacted by low-temperature stress. Uncovering the mechanisms underlying low-temperature stress response in cold-tolerant qingke varieties is crucial for qingke breeding. Herein, we conducted a comprehensive transcriptomic and metabolomic analysis on cold-sensitive (ZQ) and cold-tolerant (XL) qingke varieties under chilling and freezing treatments and identified lipid metabolism pathways as enriched in response to freezing treatment. Additionally, a significant positive correlation was observed between the expression of C-repeat (CRT) binding factor 10A (HvCBF10A) and Gly-Asp-Ser-Leu-motif lipase (HvGDSL) and the accumulation of multiple lipids. Functional analysis confirmed that HvCBF10A directly binds to HvGDSL, and silencing HvCBF10A resulted in a significant decrease in both HvGDSL and lipid levels, consequently impairing the cold tolerance. Overall, HvCBF10A and HvGDSL are functional units in actively regulating lipid metabolism to enhance freezing stress tolerance in qingke.

Influence of Heat Treatment on Tea Polyphenols and Their Impact on Improving Heat Tolerance in Drosophila melanogaster

This study investigated the potential mechanism of action of tea polyphenols (TPs), one of the major active ingredients in tea, to enhance heat resistance in Drosophila and the attenuating effect of heat treatment of TPs on their efficacy. The results showed that TPs were able to prolong the average survival time of Drosophila under high-temperature stress (p < 0.05), but the effect of TPs in prolonging the survival time of Drosophila melanogaster was significantly reduced (p < 0.05) with increasing TP heat-treatment time until it disappeared. The composition of TPs changed after heat treatment. It was also shown that the weakening of the effect of TPs in improving the heat tolerance of Drosophila was related to the decrease in the content of catechins and phenolic acids in their fractions as well as with the increase in the content of laccase. Transcriptomic analysis showed that the effect of TPs on heat tolerance in Drosophila melanogaster was closely related to the longevity regulation pathway, the neuroactive ligand-receptor interaction signaling pathway, and the drug metabolism-cytochrome P450 pathway. Metabolomics analysis showed that the effect of TP intervention in improving the body's heat tolerance was mainly related to amino acid metabolism and energy metabolism. However, thermal processing weakened the relevance of these transcriptomes and metabolomes. The present study reveals the mechanism of action by which heat-treated TPs affect the body's heat tolerance, which is important for the development and utilization of the heat-protection function of tea.

Photosynthesis regulation, cell membrane stabilization and methylglyoxal detoxification seems major altered pathways under cold stress as revealed by integrated multi-omics meta-analysis

Climate change has altered cold weather patterns, resulting in irregular cold weather conditions, and changing the global plant distribution pattern affecting plant development processes resulting in severe yield losses. Although molecular mechanisms and interconnections are quite well studied, a cumulative understanding of plant responses to cold stress (CS) is still lacking. Through meta-analysis, integration of data at the multi-omics level and its correlation with known physiological changes to map and understand the global changes in response to CS was made. Meta-analysis was conducted using the metafor R package program based on physiological parameters like relative electrolytic leakage, malondialdehyde, soluble sugar, proline and antioxidant enzymes activity. Proline and soluble sugars showed the highest (> 1.5 mean fold) change over control thus qualifying as global markers for studying CS. Surprisingly most up-regulated (> 15-fold) DEGs corresponded with the dehydrin family and glyoxalase superfamily proteins. Functional annotations of DEGs corresponded with photosynthesis and glycolysis pathway. Proteins responsible for cell signalling and increased soluble sugars were common in all the datasets studied thus correlating with the transcriptome and proteomic data. Proline and soluble sugars were positively regulated in all the metabolomics datasets. This study supported the earlier known players like proline and soluble sugars. Surprisingly, a new player glyoxalase seems to be contributing in CS.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01367-9.

Integrated Transcriptomic and Metabolomics Analyses Reveal Molecular Responses to Cold Stress in Coconut (Cocos nucifera L.) Seedlings

Coconut is an important tropical and subtropical fruit and oil crop severely affected by cold temperature, limiting its distribution and application. Thus, studying its low-temperature reaction mechanism is required to expand its cultivation range. We used growth morphology and physiological analyses to characterize the response of coconuts to 10, 20, and 30 d of low temperatures, combined with transcriptome and metabolome analysis. Low-temperature treatment significantly reduced the plant height and dry weight of coconut seedlings. The contents of soil and plant analyzer development (SPAD), soluble sugar (SS), soluble protein (SP), proline (Pro), and malondialdehyde (MDA) in leaves were significantly increased, along with the activities of superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT), and the endogenous hormones abscisic acid (ABA), auxin (IAA), zeatin (ZR), and gibberellin (GA) contents. A large number of differentially expressed genes (DEGs) (9968) were detected under low-temperature conditions. Most DEGs were involved in mitogen-activated protein kinase (MAPK) signaling pathway-plant, plant hormone signal transduction, plant-pathogen interaction, biosynthesis of amino acids, amino sugar and nucleotide sugar metabolism, carbon metabolism, starch and sucrose metabolism, purine metabolism, and phenylpropanoid biosynthesis pathways. Transcription factors (TFs), including WRKY, AP2/ERF, HSF, bZIP, MYB, and bHLH families, were induced to significantly differentially express under cold stress. In addition, most genes associated with major cold-tolerance pathways, such as the ICE-CBF-COR, MAPK signaling, and endogenous hormones and their signaling pathways, were significantly up-regulated. Under low temperatures, a total of 205 differentially accumulated metabolites (DAMs) were enriched; 206 DAMs were in positive-ion mode and 97 in negative-ion mode, mainly including phenylpropanoids and polyketides, lipids and lipid-like molecules, benzenoids, organoheterocyclic compounds, organic oxygen compounds, organic acids and derivatives, nucleosides, nucleotides, and analogues. Comprehensive metabolome and transcriptome analysis revealed that the related genes and metabolites were mainly enriched in amino acid, flavonoid, carbohydrate, lipid, and nucleotide metabolism pathways under cold stress. Together, the results of this study provide important insights into the response of coconuts to cold stress, which will reveal the underlying molecular mechanisms and help in coconut screening and breeding.

Multi-Omics Revealed Peanut Root Metabolism Regulated by Exogenous Calcium under Salt Stress

High salinity severely inhibits plant seedling root development and metabolism. Although plant salt tolerance can be improved by exogenous calcium supplementation, the metabolism molecular mechanisms involved remain unclear. In this study, we integrated three types of omics data (transcriptome, metabolome, and phytohormone absolute quantification) to analyze the metabolic profiles of peanut seedling roots as regulated by exogenous calcium under salt stress. (1) exogenous calcium supplementation enhanced the allocation of carbohydrates to the TCA cycle and plant cell wall biosynthesis rather than the shikimate pathway influenced by up-regulating the gene expression of antioxidant enzymes under salt stress; (2) exogenous calcium induced further ABA accumulation under salt stress by up-regulating the gene expression of ABA biosynthesis key enzymes AAO2 and AAO3 while down-regulating ABA glycosylation enzyme UGT71C5 expression; (3) exogenous calcium supplementation under salt stress restored the trans-zeatin absolute content to unstressed levels while inhibiting the root cis-zeatin biosynthesis.

Abscisic Acid Regulates Carbohydrate Metabolism, Redox Homeostasis and Hormonal Regulation to Enhance Cold Tolerance in Spring Barley

Abscisic acid (ABA) plays a vital role in the induction of low temperature tolerance in plants. To understand the molecular basis of this phenomenon, we performed a proteomic analysis on an ABA-deficit mutant barley (Az34) and its wild type (cv Steptoe) under control conditions (25/18 °C) and after exposure to 0 °C for 24 h. Most of the differentially abundant proteins were involved in the processes of photosynthesis and metabolisms of starch, sucrose, carbon, and glutathione. The chloroplasts in Az34 leaves were more severely damaged, and the decrease in Fv/Fm was larger in Az34 plants compared with WT under low temperature. Under low temperature, Az34 plants possessed significantly higher activities of ADP-glucose pyrophosphorylase, fructokinase, monodehydroascorbate reductase, and three invertases, but lower UDP-glucose pyrophosphorylase activity than WT. In addition, concentrations of proline and soluble protein were lower, while concentration of H₂O₂ was higher in Az34 plants compared to WT under low temperature. Collectively, the results indicated that ABA deficiency induced modifications in starch and sucrose biosynthesis and sucrolytic pathway and overaccumulation of reactive oxygen species were the main reason for depressed low temperature tolerance in barley, which provide novel insights to the response of barley to low temperature under future climate change.

Combining transcriptomics and metabolomics to identify key response genes for aluminum toxicity in the root system of Brassica napus L. seedlings

By integrating QTL mapping, transcriptomics and metabolomics, 138 hub genes were identified in rapeseed root response to aluminum stress and mainly involved in metabolism of lipids, carbohydrates and secondary metabolites. Aluminum (Al) toxicity has become one of the important abiotic stress factors in areas with acid soil, which hinders the absorption of water and nutrients by roots, and consequently retards the growth of crops. A deeper understanding of the stress-response mechanism of Brassica napus may allow us to identify the tolerance gene(s) and use this information in breeding-resistant crop varieties. In this study, a population of 138 recombinant inbred lines (RILs) was subjected to aluminum stress, and QTL (quantitative trait locus) mapping was used to preliminarily locate quantitative trait loci related to aluminum stress. Root tissues from seedlings of an aluminum-resistant (R) line and an aluminum-sensitive (S) line from the RIL population were harvested for transcriptome sequencing and metabolome determination. By combining the data on quantitative trait genes (QTGs), differentially expressed genes (DEGs), and differentially accumulated metabolites (DAMs), key candidate genes related to aluminum tolerance in rapeseed were determined. The results showed that there were 3186 QTGs in the RIL population, 14,232 DEGs and 457 DAMs in the comparison between R and S lines. Lastly, 138 hub genes were selected to have a strong positive or negative correlation with 30 important metabolites (|R|≥ 0.95). These genes were mainly involved in the metabolism of lipids, carbohydrates and secondary metabolites in response to Al toxicity stress. In summary, this study provides an effective method for screening key genes by combining QTLs, transcriptome sequencing and metabolomic analysis, but also lists key genes for exploring the molecular mechanism of Al tolerance in rapeseed seedling roots.

Transcriptomics integrated with metabolomics provides a new strategy for mining key genes in response to low temperature stress in Helictotrichon virescens

H. virescens is a perennial herbaceous plant with highly tolerant to cold weather, but the key genes that respond to low temperature stress still remain unclear. Hence, RNA-seq was performed using leaves of H. virescens treated at 0 °C and 25 °C for 12 h, 36 h, and 60 h, respectively, and a total of 9416 DEGs were significantly enriched into seven KEGG pathways. The LC-QTRAP platform was performed using leaves of H. virescens leaves at 0 °C and 25 °C for 12 h, 36 h, and 60 h, respectively, and a total of 1075 metabolites were detected, which were divided into 10 categories. Additionally, 18 major metabolites, two key pathways, and six key genes were mined using a multi-omics analytical strategy. The RT-PCR results showed that with the extension of treatment time, the expression levels of key genes in the treatment group gradually increased, and the difference between the treatment group and the control group was extremely significant. Notably, the functional verification results showed that the key genes positively regulated cold tolerance of H. virescens. These results can lay a foundation for the in-depth analysis of the mechanism of response of perennial herbs to low temperature stress.

Abiotic Stress in Crop Production

The vast majority of agricultural land undergoes abiotic stress that can significantly reduce agricultural yields. Understanding the mechanisms of plant defenses against stresses and putting this knowledge into practice is, therefore, an integral part of sustainable agriculture. In this review, we focus on current findings in plant resistance to four cardinal abiotic stressors—drought, heat, salinity, and low temperatures. Apart from the description of the newly discovered mechanisms of signaling and resistance to abiotic stress, this review also focuses on the importance of primary and secondary metabolites, including carbohydrates, amino acids, phenolics, and phytohormones. A meta-analysis of transcriptomic studies concerning the model plant Arabidopsis demonstrates the long-observed phenomenon that abiotic stressors induce different signals and effects at the level of gene expression, but genes whose regulation is similar under most stressors can still be traced. The analysis further reveals the transcriptional modulation of Golgi-targeted proteins in response to heat stress. Our analysis also highlights several genes that are similarly regulated under all stress conditions. These genes support the central role of phytohormones in the abiotic stress response, and the importance of some of these in plant resistance has not yet been studied. Finally, this review provides information about the response to abiotic stress in major European crop plants—wheat, sugar beet, maize, potatoes, barley, sunflowers, grapes, rapeseed, tomatoes, and apples.

Comparative metabolomics of root-tips reveals distinct metabolic pathways conferring drought tolerance in contrasting genotypes of rice

Background

The mechanisms underlying rice root responses to drought during the early developmental stages are yet unknown.

Results

This study aimed to determine metabolic differences in IR64, a shallow-rooting, drought-susceptible genotype, and Azucena, a drought-tolerant and deep-rooting genotype under drought stress. The morphological evaluation revealed that Azucena might evade water stress by increasing the lateral root system growth, the root surface area, and length to access water. At the same time, IR64 may rely mainly on cell wall thickening to tolerate stress. Furthermore, significant differences were observed in 49 metabolites in IR64 and 80 metabolites in Azucena, for which most metabolites were implicated in secondary metabolism, amino acid metabolism, nucleotide acid metabolism and sugar and sugar alcohol metabolism. Among these metabolites, a significant positive correlation was found between allantoin, galactaric acid, gluconic acid, glucose, and drought tolerance. These metabolites may serve as markers of drought tolerance in genotype screening programs. Based on corresponding biological pathways analysis of the differentially abundant metabolites (DAMs), biosynthesis of alkaloid-derivatives of the shikimate pathway, fatty acid biosynthesis, purine metabolism, TCA cycle and amino acid biosynthesis were the most statistically enriched biological pathway in Azucena in drought response. However, in IR64, the differentially abundant metabolites of starch and sucrose metabolism were the most statistically enriched biological pathways.

Conclusion

Metabolic marker candidates for drought tolerance were identified in both genotypes. Thus, these markers that were experimentally determined in distinct metabolic pathways can be used for the development or selection of drought-tolerant rice genotypes.

Integrative Transcriptomic and Metabolomic Analysis Reveals the Molecular Mechanism of Red Maple (Acer rubrum L.) Leaf Coloring

This study employed a combination of ultraviolet spectrophotometry, LC-ESI-MS/MS system, and RNA-sequencing technology; the extracts and isolation of total RNA from the red and yellow leaf strains of red maple (Acer rubrum L.) at different developmental stages were subjected to an intercomparison of the dynamic content of chlorophyll and total anthocyanin, flavonoid metabolite fingerprinting, and gene expression. The metabonomic results indicated that one hundred and ninety-two flavonoids were identified, which could be classified into eight categories in the red maple leaves. Among them, 39% and 19% were flavones and flavonols, respectively. The metabolomic analysis identified 23, 32, 24, 24, 38, and 41 DAMs in the AR1018r vs. AR1031r comparison, the AR1018r vs. AR1119r comparison, the AR1031r vs. AR1119r comparison, the AR1018y vs. AR1031y comparison, the AR1018y vs. AR1119y comparison, and the AR1031y vs. AR1119y comparison, respectively. In total, 6003 and 8888 DEGs were identified in AR1018r vs. AR1031r comparison and in the AR1018y vs. AR1031y comparison, respectively. The GO and KEGG analyses showed that the DEGs were mainly involved in plant hormone signal transduction, flavonoid biosynthesis, and other metabolite metabolic processes. The comprehensive analysis revealed that caffeoyl-CoA 3-O-methyltransferase (Cluster-28704.45358 and Cluster-28704.50421) was up-regulated in the red strain but down-regulated in the yellow strain, while Peonidin 3-O-glucoside chloride and Pelargonidin 3-O-beta-D-glucoside were up-regulated in both the red and yellow strains. By successfully integrating the analyses on the behavior of pigment accumulation, dynamics of flavonoids, and differentially expressed genes with omics tools, the regulation mechanisms underlying leaf coloring in red maple at the transcriptomic and metabolomic levels were demonstrated, and the results provide valuable information for further research on gene function in red maple.

Integrated transcriptomics and metabolomics analysis reveals key regulatory network that response to cold stress in common Bean (Phaseolus vulgaris L.)

Cold temperatures can be detrimental to crop survival and productivity. Breeding progress can be improved by understanding the molecular basis of low temperature tolerance. We investigated the key routes and critical metabolites related to low temperature resistance in cold-tolerant and -sensitive common bean cultivars 120 and 093, respectively. Many potential genes and metabolites implicated in major metabolic pathways during the chilling stress response were identified through transcriptomics and metabolomics research. Under chilling stress, the expression of many genes involved in lipid, amino acid, and flavonoid metabolism, as well as metabolite accumulation increased in the two bean types. Malondialdehyde (MDA) content was lower in 120 than in 093. Regarding amino acid metabolism, 120 had a higher concentration of acidic amino acids than 093, whereas 093 had a higher concentration of basic amino acids. Methionine accumulation was clearly higher in 120 than in 093. In addition, 120 had a higher concentration of many types of flavonoids than 093. Flavonoids, methionine and malondialdehyde could be used as biomarkers of plant chilling injury. Transcriptome analysis of hormone metabolism revealed considerably greater, expression of abscisic acid (ABA), gibberellin (GA), and jasmonic acid (JA) in 093 than in 120 during chilling stress, indicating that hormone regulation modes in 093 and 120 were different. Thus, chilling stress tolerance is different between 093 and 120 possibly due to transcriptional and metabolic regulation.

Integrated Metabolomic and Transcriptomic Analysis Reveals That Amino Acid Biosynthesis May Determine Differences in Cold-Tolerant and Cold-Sensitive Tea Cultivars

Cold stress is one of the major abiotic stresses limiting tea production. The planting of cold-resistant tea cultivars is one of the most effective measures to prevent chilling injury. However, the differences in cold resistance between tea cultivars remain unclear. In the present study, we perform a transcriptomic and metabolomic profiling of Camellia sinensis var. "Shuchazao" (cold-tolerant, SCZ) and C. sinensis var. assamica "Yinghong 9" (cold-sensitive, YH9) during cold acclimation and analyze the correlation between gene expression and metabolite biosynthesis. Our results show that there were 51 differentially accumulated metabolites only up-regulated in SCZ in cold-acclimation (CA) and de-acclimation (DA) stages, of which amino acids accounted for 18%. The accumulation of L-arginine and lysine in SCZ in the CA stage was higher than that in YH9. A comparative transcriptomic analysis showed an enrichment of the amino acid biosynthesis pathway in SCZ in the CA stage, especially "arginine biosynthesis" pathways. In combining transcriptomic and metabolomic analyses, it was found that genes and metabolites associated with amino acid biosynthesis were significantly enriched in the CA stage of SCZ compared to CA stage of YH9. Under cold stress, arginine may improve the cold resistance of tea plants by activating the polyamine synthesis pathway and CBF (C-repeat-binding factor)-COR (cold-regulated genes) regulation pathway. Our results show that amino acid biosynthesis may play a positive regulatory role in the cold resistance of tea plants and assist in understanding the cold resistance mechanism differences among tea varieties.

Transcriptomic profiling of wheat stem during meiosis in response to freezing stress

Low temperature injury in spring has seriously destabilized the production and grain quality of common wheat. However, the molecular mechanisms underlying spring frost tolerance remain elusive. In this study, we investigated the response of a frost-tolerant wheat variety Zhongmai8444 to freezing stress at the meiotic stage. Transcriptome profiles over a time course were subsequently generated by high-throughput sequencing. Our results revealed that the prolonged freezing temperature led to the significant reductions in plant height and seed setting rate. Cell wall thickening in the vascular tissue was also observed in the stems. RNA-seq analyses demonstrated the identification of 1010 up-regulated and 230 down-regulated genes shared by all time points of freezing treatment. Enrichment analysis revealed that gene activity related to hormone signal transduction and cell wall biosynthesis was significantly modulated under freezing. In addition, among the identified differentially expressed genes, 111 transcription factors belonging to multiple gene families exhibited dynamic expression pattern. This study provided valuable gene resources beneficial for the breeding of wheat varieties with improved spring frost tolerance.

Transcriptomics integrated with widely targeted metabolomics reveals the cold resistance mechanism in Hevea brasiliensis

The rubber tree is the primary source of natural rubber and is mainly cultivated in Southeast Asian countries. Low temperature is the major abiotic stress affecting the yield of the rubber tree. Therefore, uncovering the cold resistance mechanism in the rubber tree is necessary. The present study used RNA-sequencing technology and ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) to analyze the transcriptomic and metabolomic changes in two rubber tree clones with different cold resistance capacities (temperature-sensitive Reyan 8-79 and cold-resistant Yunyan 77-4) at 0 h, 2 h, 6 h, and 20 h of exposure to 4°C. Independent analysis of the transcriptome and metabolitome showed that under prolonged low-temperature treatment, Yunyan 77-4 expressed more genes involved in regulating enzyme activity, changing cell permeability, and synthesizing significant metabolites, such as flavonoids and amino acids, than Reyan 8-79. The KEGG annotation and enrichment analysis identified arginine metabolism and biosynthesis of flavonoids as the major pathway associated with cold resistance. Integrated transcriptome and metabolome analysis showed that the increase in the expression of genes modulated flavonoid biosynthesis, arginine biosynthesis, and anthocyanins biosynthesis, resulting in higher levels of metabolites, such as naringenin chalcone, apigenin, dihydroquercetin, cyanidin 3-glucoside, L-arginosuccinate, N-acetyl-ornithine, ornithine, and N-acetyl-glutamate, in Yunyan 77-4 than in Reyan 8-79 after prolonged low-temperature treatment. Phylogenetic analysis identified the genes, such as CHS (gene356) and F3H (gene33147) of flavonoid biosynthesis and NAGS (gene16028, gene33765), ArgC (gene2487), and ASS (gene6161) of arginine biosynthesis were the key genes involved in the cold resistant of rubber tree. Thus, the present study provides novel insights into how rubber clones resist cold and is a valuable reference for cold-resistance breeding.

Cloning of TaeRF1 gene from Caucasian clover and its functional analysis responding to low-temperature stress

Low temperature (LT) is an important threat to the normal growth of plants. In this study, based on the full-length transcriptome sequencing results, the cold resistance genes were cloned from Caucasian clover with strong cold resistance. We cloned the CDS of TaeRF1, which is 1311 bp in length and encodes 436 amino acids. The molecular weight of the protein is 48.97 kDa, which had no transmembrane structure, and its isoelectric point (pI) was 5.42. We predicted the structure of TaeRF1 and found 29 phosphorylation sites. Subcellular localization showed that TaeRF1 was localized and expressed in cell membrane and chloroplasts. The TaeRF1 gene was induced by stress due to cold, salt, alkali and drought and its expression level was higher in roots and it was more sensitive to LT. Analysis of transgenic A. thaliana plants before and after LT treatment showed that the TaeRF1 gene enhanced the removal of excess H₂O₂, and increased the activity of antioxidant enzymes, thus improving the plant's ability to resist stress. Additionally, the OE lines showed increased cold tolerance by upregulating the transcription level of cold-responsive genes (CBF1, CBF2, COR15B, COR47, ICE1, and RD29A). This study demonstrates that TaeRF1 is actively involved in the responses of plants to LT stress. We also provide a theoretical basis for breeding and a potential mechanism underlying the responses of Caucasian clover to abiotic stress.

Genome-Wide Identification of bHLH Transcription Factor in Medicago sativa in Response to Cold Stress

Alfalfa represents one of the most important legume forages, and it is also applied as an organic fertilizer to improve soil quality. However, this perennial plant is native to warmer temperate regions, and its valuable cold-acclimation-related regulatory mechanisms are still less known. In higher plants, the bHLH transcription factors play pleiotropic regulatory roles in response to abiotic stresses. The recently released whole genome sequencing data of alfalfa allowed us to identify 469 MsbHLHs by multi-step homolog search. Herein, we primarily identified 65 MsbHLH genes that significantly upregulated under cold stress, and such bHLHs were classified into six clades according to their expression patterns. Interestingly, the phylogenetic analysis and conserved motif screening of the cold-induced MsbHLHs showed that the expression pattern is relatively varied in each bHLH subfamily, this result indicating that the 65 MsbHLHs may be involved in a complex cold-responsive regulatory network. Hence, we analyzed the TFBSs at promoter regions that unraveled a relatively conserved TFBS distribution with genes exhibiting similar expression patterns. Eventually, to verify the core components involved in long-term cold acclimation, we examined transcriptome data from a freezing-tolerant species (cv. Zhaodong) in the field and compared the expression of cold-sensitive/tolerant subspecies of alfalfa, giving 11 bHLH as candidates, which could be important for further cold-tolerance enhancement and molecular breeding through genetic engineering in alfalfa.

Physiology and Molecular Breeding in Sustaining Wheat Grain Setting and Quality under Spring Cold Stress

Spring cold stress (SCS) compromises the reproductive growth of wheat, being a major constraint in achieving high grain yield and quality in winter wheat. To sustain wheat productivity in SCS conditions, breeding cultivars conferring cold tolerance is key. In this review, we examine how grain setting and quality traits are affected by SCS, which may occur at the pre-anthesis stage. We have investigated the physiological and molecular mechanisms involved in floret and spikelet SCS tolerance. It includes the protective enzymes scavenging reactive oxygen species (ROS), hormonal adjustment, and carbohydrate metabolism. Lastly, we explored quantitative trait loci (QTLs) that regulate SCS for identifying candidate genes for breeding. The existing cultivars for SCS tolerance were primarily bred on agronomic and morphophysiological traits and lacked in molecular investigations. Therefore, breeding novel wheat cultivars based on QTLs and associated genes underlying the fundamental resistance mechanism is urgently needed to sustain grain setting and quality under SCS.

Comparative Analysis of Antioxidant Accumulation under Cold Acclimation, Deacclimation and Reacclimation in Winter Wheat

Low temperature during cold acclimation (CA) leads to the accumulation of detrimental reactive oxygen species (ROS) in plant tissues, which are scavenged by antioxidants, such as ascorbate and glutathione. However, there is a lack of studies examining the dynamics of antioxidants throughout CA, deacclimation (DEA), and reacclimation (REA) in winter wheat. Six winter wheat genotypes were selected to assess the effect of CA, DEA, and REA on the concentrations of ascorbate and glutathione in leaf and crown tissues under two CA temperature treatments. Higher levels of total, reduced, and oxidised ascorbate were determined in leaves, whereas crowns accumulated higher concentrations of nicotinamide adenine dinucleotide (NAD⁺) after CA, DEA, and REA. Constant low temperature (CLT) during CA led to higher contents of ascorbate and glutathione in both tissues at all stages of acclimation, in comparison with prolonged higher low temperature (PHLT). The concentrations of antioxidants increased after CA, tended to decrease after DEA, and returned to CA levels after REA. Significant positive correlations between freezing tolerance (FT) and antioxidants were only determined under the CA at CLT treatment, thus, affirming the negative effect of PHLT during CA on the FT of winter wheat.

The effect of cold acclimation, deacclimation and reacclimation on metabolite profiles and freezing tolerance in winter wheat

Global climate change will cause longer and warmer autumns, thus negatively affecting the quality of cold acclimation (CA) and reducing the freezing tolerance (FT) of winter wheat. Insufficient FT and fluctuating temperatures during winter can accelerate the deacclimation (DEA) process, whereas reacclimation (REA) is possible only while the vernalization requirement is unfulfilled. Six winter wheat genotypes with different winter hardiness profiles were used to evaluate the impact of constant low-temperature (2°C) and prolonged higher low-temperature (28 days at 10°C followed by 2°C until day 49) on shoot biomass and metabolite accumulation patterns in leaf and crown tissues throughout 49 days of CA, 7 days of DEA, and 14 days of REA. The FT of winter wheat was determined as LT₃₀ values by conducting freezing tests after CA, DEA, and REA. Shoot biomass accumulation, projected as the green leaf area (GLA), was investigated by non-destructive RGB imaging-based phenotyping. Dynamics of carbohydrates, hexose phosphates, organic acids, proteins, and amino acids were assessed in leaf and crown tissues. Results revealed that exposure to higher low-temperature induced higher accumulation of shoot biomass and had a negative impact on FT of winter wheat. Prolonged higher low-temperature negatively affected the accumulation of soluble carbohydrates, protein content and amino acids, and had a positive effect on starch accumulation in leaf and crown tissues after CA, in comparison with the constant low-temperature treatment. DEA resulted in significantly reduced FT. Lower concentrations of glucose-6-phosphate, sucrose and proline, as well as higher concentrations of starch in leaves and crowns were found after DEA. The majority of the genotypes regained FT after REA; higher concentrations of glucose and malate in leaves, and sucrose in crown tissue were observed, whereas starch accumulation was decreased in both tissues. Negative correlations were determined between FT and starch concentration in leaves and crowns, while proline and proteins, accumulated in crowns, showed positive correlations with FT. This study broadens the knowledge regarding the effect of different low-temperature regimes on the dynamics of metabolite accumulation in winter wheat throughout CA, DEA, and REA, and its relationship to biomass accumulation and FT.

Transcription-associated metabolomic profiling reveals the critical role of frost tolerance in wheat

BACKGROUND

Low temperature is a crucial stress factor of wheat (Triticum aestivum L.) and adversely impacts on plant growth and grain yield. Multi-million tons of grain production are lost annually because crops lack the resistance to survive in winter. Particularlly, winter wheat yields was severely damaged under extreme cold conditions. However, studies about the transcriptional and metabolic mechanisms underlying cold stresses in wheat are limited so far.

RESULTS

In this study, 14,466 differentially expressed genes (DEGs) were obtained between wild-type and cold-sensitive mutants, of which 5278 DEGs were acquired after cold treatment. 88 differential accumulated metabolites (DAMs) were detected, including P-coumaroyl putrescine of alkaloids, D-proline betaine of mino acids and derivativ, Chlorogenic acid of the Phenolic acids. The comprehensive analysis of metabolomics and transcriptome showed that the cold resistance of wheat was closely related to 13 metabolites and 14 key enzymes in the flavonol biosynthesis pathway. The 7 enhanced energy metabolites and 8 up-regulation key enzymes were also compactly involved in the sucrose and amino acid biosynthesis pathway. Moreover, quantitative real-time PCR (qRT-PCR) revealed that twelve key genes were differentially expressed under cold, indicating that candidate genes POD, Tacr7, UGTs, and GSTU6 which were related to cold resistance of wheat.

CONCLUSIONS

In this study, we obtained the differentially expressed genes and differential accumulated metabolites in wheat under cold stress. Using the DEGs and DAMs, we plotted regulatory pathway maps of the flavonol biosynthesis pathway, sucrose and amino acid biosynthesis pathway related to cold resistance of wheat. It was found that candidate genes POD, Tacr7, UGTs and GSTU6 are related to cold resistance of wheat. This study provided valuable molecular information and new genetic engineering clues for the further study on plant resistance to cold stress.

Combined Analysis of BSA-Seq Based Mapping, RNA-Seq, and Metabolomic Unraveled Candidate Genes Associated with Panicle Grain Number in Rice (Oryza sativa L.)

Rice grain yield is a complex and highly variable quantitative trait consisting of several key components, including the grain weight, the effective panicles per unit area, and the grain number per panicle (GNPP). The GNPP is a significant contributor to grain yield controlled by multiple genes (QTL) and is crucial for improvement. Attempts have been made to find genes for this trait, which has always been a challenging and arduous task through conventional methods. We combined a BSA analysis, RNA profiling, and a metabolome analysis in the present study to identify new candidate genes involved in the GNPP. The F₂ population from crossing R4233 (high GNPP) and Ce679 (low GNPP) revealed a frequency distribution fitting two segregated genes. Three pools, including low, middle, and high GNPP, were constructed and a BSA analysis revealed six candidate regions spanning 5.38 Mb, containing 739 annotated genes. Further, a conjunctive analysis of BSA-Seq and RNA-Seq showed 31 differentially expressed genes (DEGs) in the candidate intervals. Subsequently, a metabolome analysis showed 1024 metabolites, with 71 significantly enriched, including 44 up and 27 downregulated in Ce679 vs. R4233. A KEGG enrichment analysis of these 31 DEGs and 71 differentially enriched metabolites (DEMs) showed two genes, Os12g0102100 and Os01g0580500, significantly enriched in the metabolic pathways’ biosynthesis of secondary metabolites, cysteine and methionine metabolism, and fatty acid biosynthesis. Os12g0102100, which encodes for the alcohol dehydrogenase superfamily and a zinc-containing protein, is a novel gene whose contribution to the GNPP is not yet elucidated. This gene coding for mitochondrial trans-2-enoyl-CoA reductase is involved in the biosynthesis of myristic acid, also known as tetradecanoic acid. The Os01g0580500 coding for the enzyme 1-aminoclopropane-1-carboxylate oxidase (OsACO7) is responsible for the final step of the ethylene biosynthesis pathway through the conversion of 1-aminocyclopropane-1-carboxylic acid (ACC) into ethylene. Unlike Os12g0102100, this gene was significantly upregulated in R4233, downregulated in Ce679, and significantly enriched in two of the three metabolite pathways. This result pointed out that these two genes are responsible for the difference in the GNPP in the two cultivars, which has never been identified. Further validation studies may disclose the physiological mechanisms through which they regulate the GNPP in rice.

Transcriptome analysis of the winter wheat Dn1 in response to cold stress

BACKGROUND

Heilongjiang Province has a long and cold winter season (the minimum temperature can reach -30 ℃), and few winter wheat varieties can safely overwinter. Dongnongdongmai1 (Dn1) is the first winter wheat variety that can safely overwinter in Heilongjiang Province. This variety fills the gap for winter wheat cultivation in the frigid region of China and greatly increases the land utilization rate. To understand the molecular mechanism of the cold response, we conducted RNA-sequencing analysis of Dn1 under cold stress.

RESULTS

Approximately 120,000 genes were detected in Dn1 under cold stress. The numbers of differentially expressed genes (DEGs) in the six comparison groups (0 ℃ vs. 5 ℃, -5 ℃ vs. 5 ℃, -10 ℃ vs. 5 ℃, -15 ℃ vs. 5 ℃, -20 ℃ vs. 5 ℃ and -25 ℃ vs. 5 ℃) were 11,313, 8313, 15,636, 13,671, 14,294 and 13,979, respectively. Gene Ontology functional annotation suggested that the DEGs under cold stress mainly had "binding", "protein kinase" and "catalytic" activities and were involved in "oxidation-reduction", "protein phosphorylation" and "carbohydrate metabolic" processes. Kyoto Encyclopedia of Genes and Genomes enrichment analysis indicated that the DEGs performed important functions in cold signal transduction and carbohydrate metabolism. In addition, major transcription factors (AP2/ERF, bZIP, NAC, WRKY, bHLH and MYB) participating in the Dn1 cold stress response were activated by low temperature.

CONCLUSION

This is the first study to explore the Dn1 transcriptome under cold stress. Our study comprehensively analysed the key genes involved in cold signal transduction and carbohydrate metabolism in Dn1 under cold stress. The results obtained by transcriptome analysis could help to further explore the cold resistance mechanism of Dn1 and provide basis for breeding of cold-resistant crops.

Advances in Metabolomics-Driven Diagnostic Breeding and Crop Improvement

Climate change continues to threaten global crop output by reducing annual productivity. As a result, global food security is now considered as one of the most important challenges facing humanity. To address this challenge, modern crop breeding approaches are required to create plants that can cope with increased abiotic/biotic stress. Metabolomics is rapidly gaining traction in plant breeding by predicting the metabolic marker for plant performance under a stressful environment and has emerged as a powerful tool for guiding crop improvement. The advent of more sensitive, automated, and high-throughput analytical tools combined with advanced bioinformatics and other omics techniques has laid the foundation to broadly characterize the genetic traits for crop improvement. Progress in metabolomics allows scientists to rapidly map specific metabolites to the genes that encode their metabolic pathways and offer plant scientists an excellent opportunity to fully explore and rationally harness the wealth of metabolites that plants biosynthesize. Here, we outline the current application of advanced metabolomics tools integrated with other OMICS techniques that can be used to: dissect the details of plant genotype–metabolite–phenotype interactions facilitating metabolomics-assisted plant breeding for probing the stress-responsive metabolic markers, explore the hidden metabolic networks associated with abiotic/biotic stress resistance, facilitate screening and selection of climate-smart crops at the metabolite level, and enable accurate risk-assessment and characterization of gene edited/transgenic plants to assist the regulatory process. The basic concept behind metabolic editing is to identify specific genes that govern the crucial metabolic pathways followed by the editing of one or more genes associated with those pathways. Thus, metabolomics provides a superb platform for not only rapid assessment and commercialization of future genome-edited crops, but also for accelerated metabolomics-assisted plant breeding. Furthermore, metabolomics can be a useful tool to expedite the crop research if integrated with speed breeding in future.

Increase Crop Resilience to Heat Stress Using Omic Strategies

Varieties of various crops with high resilience are urgently needed to feed the increased population in climate change conditions. Human activities and climate change have led to frequent and strong weather fluctuation, which cause various abiotic stresses to crops. The understanding of crops' responses to abiotic stresses in different aspects including genes, RNAs, proteins, metabolites, and phenotypes can facilitate crop breeding. Using multi-omics methods, mainly genomics, transcriptomics, proteomics, metabolomics, and phenomics, to study crops' responses to abiotic stresses will generate a better, deeper, and more comprehensive understanding. More importantly, multi-omics can provide multiple layers of information on biological data to understand plant biology, which will open windows for new opportunities to improve crop resilience and tolerance. However, the opportunities and challenges coexist. Interpretation of the multidimensional data from multi-omics and translation of the data into biological meaningful context remained a challenge. More reasonable experimental designs starting from sowing seed, cultivating the plant, and collecting and extracting samples were necessary for a multi-omics study as the first step. The normalization, transformation, and scaling of single-omics data should consider the integration of multi-omics. This review reports the current study of crops at abiotic stresses in particular heat stress using omics, which will help to accelerate crop improvement to better tolerate and adapt to climate change.

Proteomic analysis reveals the molecular mechanism underlying the cold acclimation and freezing tolerance of wheat (Triticum aestivum L.)

Cold acclimation (CA) is an important evolutionary adaptive mechanism for wheat freezing resistence. To clarify the molecular basis of wheat CA and freezing tolerance, the effects of CA (4 °C) and non-CA (20 °C) treatments and freezing stress (-5 °C) on the proteins in the wheat crown were characterized via an iTRAQ-based proteomic analysis. A total of 669 differentially accumulated proteins (DAPs) were identified after the CA, of which seven were also DAPs in the CA plants exposed to freezing stress. Additionally, the 15 DAPs in the CA group and the 23 DAPs in the non-CA group after the freezing treatment differed substantially. Functional analyses indicated that CA enhanced freezing tolerance by regulating proteins involved in signal transduction, carbohydrate metabolism, stress and defense responses, and phenylpropanoid biosynthesis. An integrated transcriptomic, proteomic, and metabolomic analysis revealed significant changes in various components of the glutathione metabolic pathway. The overexpression and silencing of Wdhn13 in Arabidopsis and wheat resulted in increased tolerance and sensitivity to freezing stress, respectively, suggesting Wdhn13 promotes freezing tolerance. Overall, our study offers insights into the regulatory network underlying the CA and freezing tolerance of wheat, which may be useful for elucidating wheat freezing resistance.

Genomic Survey and Cold-Induced Expression Patterns of bHLH Transcription Factors in Liriodendron chinense (Hemsl) Sarg

bHLH transcription factors play an animated role in the plant kingdom during growth and development, and responses to various abiotic stress. In this current study, we conducted, the genome-wide survey of bHLH transcription factors in Liriodendron chinense (Hemsl) Sarg., 91 LcbHLH family members were identified. Identified LcbHLH gene family members were grouped into 19 different subfamilies based on the conserved motifs and phylogenetic analysis. Our results showed that LcbHLH genes clustered in the same subfamily exhibited a similar conservative exon-intron pattern. Hydrophilicity value analysis showed that all LcbHLH proteins were hydrophilic. The Molecular weight (Mw) of LcbHLH proteins ranged from 10.19 kD (LcbHLH15) to 88.40 kD (LcbHLH50). A greater proportion, ~63%, of LcbHLH proteins had a theoretical isoelectric point (pI) less than seven. Additional analysis on the collinear relationships within species and among dissimilar species illustrated that tandem and fragment duplication are the foremost factors of amplification of this family in the evolution process, and they are all purified and selected. RNA-seq and real-time quantitative PCR analysis of LcbHLH members showed that the expression of LcbHLH35, 55, and 86 are up-regulated, and the expression of LcbHLH9, 20, 39, 54, 56, and 69 is down-regulated during cold stress treatments while the expression of LcbHLH24 was up-regulated in the short term and then later down-regulated. From our results, we concluded that LcbHLH genes might participate in cold-responsive processes of L. chinense. These findings provide the basic information of bHLH gene in L. chinense and their regulatory roles in plant development and cold stress response.

Transcriptome and Metabolome Analyses Reveal Molecular Responses of Two Pepper (Capsicum annuum L.) Cultivars to Cold Stress

Low temperature is a significant factor affecting field-grown pepper. The molecular mechanisms behind peppers' response to cold stress remain unknown. Transcriptomic and metabolomic analyses were used to investigate the responses of two pepper cultivars, XS (cold-sensitive) and GZ (cold-resistant), to cold stress; these were screened from 45 pepper materials. In this study, compared with the control group (0 h), we identified 10,931 differentially expressed genes (DEGs) in XS and GZ, 657 differentially expressed metabolites (DEMs) in the positive ion mode, and 390 DEMs in the negative ion mode. Most DEGs were involved in amino acid biosynthesis, plant hormone signal transduction, and the mitogen-activated protein kinase (MAPK) signaling pathway. Furthermore, metabolomic analysis revealed that the content of free polyamines (PAs), plant hormones, and osmolytes, mainly contained increased putrescine, spermine, spermidine, abscisic acid (ABA), jasmonic acid (JA), raffinose, and proline, in response to cold stress. Importantly, the regulation of the ICE (inducer of CBF expression)-CBF (C repeat binding factors)-COR (cold regulated) pathway by Ca²⁺ signaling, MAPK signaling, and reactive oxygen species (ROS) signaling plays a key role in regulating responses of peppers to cold stress. Above all, the results of the present study provide important insights into the response of peppers to cold stress, which will reveal the potential molecular mechanisms and contribute to pepper screening and breeding in the future.

Integrated Transcriptomics and Widely Targeted Metabolomics Analyses Provide Insights Into Flavonoid Biosynthesis in the Rhizomes of Golden Buckwheat (Fagopyrum cymosum)

Golden buckwheat (Fagopyrum cymosum) is used in Traditional Chinese Medicine. It has received attention because of the high value of its various medicinal and nutritional metabolites, especially flavonoids (catechin and epicatechin). However, the metabolites and their encoding genes in golden buckwheat have not yet been identified in the global landscape. This study performed transcriptomics and widely targeted metabolomics analyses for the first time on rhizomes of golden buckwheat. As a result, 10,191 differentially expressed genes (DEGs) and 297 differentially regulated metabolites (DRMs) were identified, among which the flavonoid biosynthesis pathway was enriched in both transcriptome and metabolome. The integration analyses of the transcriptome and the metabolome revealed a network related to catechin, in which four metabolites and 14 genes interacted with each other. Subsequently, an SG5 R2R3-MYB transcription factor, named FcMYB1, was identified as a transcriptional activator in catechin biosynthesis, as it was positively correlated to eight flavonoid biosynthesis genes in their expression patterns and was directly bound to the promoters of FcLAR2 and FcF3'H1 by yeast one hybrid analysis. Finally, a flavonoid biosynthesis pathway was proposed in the rhizomes of golden buckwheat, including 13 metabolites, 11 genes encoding 9 enzymes, and 1 MYB transcription factor. The expression of 12 DEGs were validated by qRT-PCR, resulting in a good agreement with the Pearson R ranging from 0.83 to 1. The study provided a comprehensive flavonoid biosynthesis and regulatory network of golden buckwheat.

Comparative transcriptome profiling of a resistant vs susceptible bread wheat (Triticum aestivum L.) cultivar in response to water deficit and cold stress

Bread wheat (Triticum aestivum L.) is one of the most important agricultural plants wearing abiotic stresses, such as water deficit and cold, that cause its productivity reduction. Since resistance to abiotic factors is a multigenic trait, therefore modern genome-wide approaches can help to involve various genetic material in breeding. One technique is full transcriptome analysis that reveals groups of stress response genes serving marker-assisted selection markers. Comparing transcriptome profiles of the same genetic material under several stresses is essential and makes the whole picture. Here, we addressed this by studying the transcriptomic response to water deficit and cold stress for two evolutionarily distant bread wheat varieties: stress-resistant cv. Saratovskaya 29 (S29) and stress-sensitive cv. Yanetzkis Probat (YP). For the first time, transcriptomes for these cultivars grown under abiotic stress conditions were obtained using Illumina based MACE technology. We identified groups of genes involved in response to cold and water deficiency stresses, including responses to each stress factor and both factors simultaneously that may be candidates for resistance genes. We discovered a core group of genes that have a similar pattern of stress-induced expression changes. The particular expression pattern was revealed not only for the studied varieties but also for the published transcriptomic data on cv. Jing 411 and cv. Fielder. Comparative transcriptome profiling of cv. S29 and cv. YP in response to water deficit and cold stress confirmed the hypothesis that stress-induced expression change is unequal within a homeologous gene group. As a rule, at least one changed significantly while the others had a relatively lower expression. Also, we found several SNPs distributed throughout the genomes of cv. S29 and cv. YP and distinguished the studied varieties from each other and the reference cv. Chinese Spring. Our results provide new data for genomics-assisted breeding of stress-tolerant wheat cultivars.

Comparative Multi-Omics of Tender Shoots from a Novel Evergrowing Tea Cultivar Provide Insight into the Winter Adaptation Mechanism

Tea (Camellia sinensis [L.] O. Kuntze) tree is a perennial plant in which winter dormancy is an important biological adaptation to environmental changes. We discovered and reported a novel tea tree cultivar that can generate tender shoots in winter several years ago, but the molecular mechanism for this unique phenotype remains unknown . Here, we conducted comparative transcriptomics, proteomics and metabolomics along with phytohormone quantitation between the winter and spring tender shoots to investigate the physiological basis and putative regulatory mechanisms of its evergrowing character during winter. Our multi-omics study has led to the following findings. Gibberellin (GA) levels and key enzymes for GA biosynthesis and the signal transduction pathway were increased in the winter shoots, causing the ABA/GA content ratio to decrease, which might play a key regulatory role in maintaining normal growth during winter. The abundance of proteins, genes and metabolites involved in energy metabolism was all increased in winter shoots, indicating that energy is critical for continuous growth under the relatively weak-light and low-temperature environment. Abiotic resistance-related proteins and free amino acids were also increased in abundance in the winter shoots, which possibly represents an adaptation response to winter conditions. These results allowed us to hypothesize a novel molecular mechanism of adaptation for this unique tender shoot evergrowing in winter.

Maize metabolome and proteome responses to controlled cold stress partly mimic early-sowing effects in the field and differ from those of Arabidopsis

In Northern Europe, sowing maize one-month earlier than current agricultural practices may lead to moderate chilling damage. However, studies of the metabolic responses to low, non-freezing, temperatures remain scarce. Here, genetically-diverse maize hybrids (Zea mays, dent inbred lines crossed with a flint inbred line) were cultivated in a growth chamber at optimal temperature and then three decreasing temperatures for 2 days each, as well as in the field. Leaf metabolomic and proteomic profiles were determined. In the growth chamber, 50% of metabolites and 18% of proteins changed between 20 and 16°C. These maize responses, partly differing from those of Arabidopsis to short-term chilling, were mapped on genome-wide metabolic maps. Several metabolites and proteins showed similar variation for all temperature decreases: seven MS-based metabolite signatures and two proteins involved in photosynthesis decreased continuously. Several increasing metabolites or proteins in the growth-chamber chilling conditions showed similar trends in the early-sowing field experiment, including trans-aconitate, three hydroxycinnamate derivatives, a benzoxazinoid, a sucrose synthase, lethal leaf-spot 1 protein, an allene oxide synthase, several glutathione transferases and peroxidases. Hybrid groups based on field biomass were used to search for the metabolite or protein responses differentiating them in growth-chamber conditions, which could be of interest for breeding.

Integrated transcriptome and metabolome analyses of biochar-induced pathways in response to Fusarium wilt infestation in pepper

The present study used soils contaminated with Fusarium oxysporum f. sp. capsici (CCS) and CCS amended with bamboo biochar (CCS + BC) to grow the pepper variety Qujiao No.1. The physiological performance, and transcriptome and metabolome profiling in leaf (L) and fruit (F) of Qujiao No.1 were conducted. Application of biochar improved soil properties, pepper plant nutrition and increased activities of enzymes related to pest/disease resistance, leading to superior physiological performance and lesser F. wilt disease incidence than plants from CCS. Most of the differentially expressed genes (DEGs) and differentially accumulated metabolites (DAMs) were involved in protein processing in endoplasmic reticulum (fruit), plant pathogen interaction (fruit), photosynthesis (leaf), phenylpropanoid biosynthesis (both tissues) and metabolic pathways (both tissues). Biochar improved plant photosynthesis, enhanced the immune system, energy production and increased stress signaling pathways. Overall, our results provide evidence of a number of pathways induced by biochar in pepper regulating its response to F. wilt disease.

Can omics deliver temperature resilient ready-to-grow crops?

Plants are extensively well-thought-out as the main source for nourishing natural life on earth. In the natural environment, plants have to face several stresses, mainly heat stress (HS), chilling stress (CS) and freezing stress (FS) due to adverse climate fluctuations. These stresses are considered as a major threat for sustainable agriculture by hindering plant growth and development, causing damage, ultimately leading to yield losses worldwide and counteracting to achieve the goal of "zero hunger" proposed by the Food and Agricultural Organization (FAO) of the United Nations. Notably, this is primarily because of the numerous inequities happening at the cellular, molecular and/or physiological levels, especially during plant developmental stages under temperature stress. Plants counter to temperature stress via a complex phenomenon including variations at different developmental stages that comprise modifications in physiological and biochemical processes, gene expression and differences in the levels of metabolites and proteins. During the last decade, omics approaches have revolutionized how plant biologists explore stress-responsive mechanisms and pathways, driven by current scientific developments. However, investigations are still required to explore numerous features of temperature stress responses in plants to create a complete idea in the arena of stress signaling. Therefore, this review highlights the recent advances in the utilization of omics approaches to understand stress adaptation and tolerance mechanisms. Additionally, how to overcome persisting knowledge gaps. Shortly, the combination of integrated omics, genome editing, and speed breeding can revolutionize modern agricultural production to feed millions worldwide in order to accomplish the goal of "zero hunger."

Joint Transcriptomic and Metabolic Analysis of Flavonoids in Cyclocarya paliurus Leaves

Flavonoids are a class of commonly occurring natural compounds in the plant kingdom with various biological activities. This study compares the content of flavonoids in Cyclocarya paliurus at different developmental stages to better inform the selection of the optimal picking period. Thus, we analyzed the transcriptome and metabolome of C. paliurus at different developmental stages. The transcriptome analysis revealed 44 genes involved in the biosynthesis of flavonoids in C. paliurus, with 10 differentially expressed genes across the four different developmental stages. The metabolites were separated and identified by a combination of chromatography and mass spectrometry, followed by multi-reaction monitoring mode analysis of triple quadrupole mass spectrometry for complete metabolite quantification. In the flavonoid synthesis pathway, a total of 137 differential flavonoids were detected. The joint transcriptome and metabolome analysis showed that the expression trends in differential metabolites and genes were significantly related. Four MYB transcription factors and two bHLH transcription factors that are closely related to flavonoid biosynthesis were identified. The regulation network of flavonoid biosynthesis in C. paliurus was thus established, providing guidance for follow-up research.

Genome-Wide Identification and Low Temperature Responsive Pattern of Actin Depolymerizing Factor (ADF) Gene Family in Wheat (Triticum aestivum L.)

The actin depolymerizing factor (ADF) gene family, which is conserved in eukaryotes, is important for plant development, growth, and stress responses. Cold stress restricts wheat growth, development, and distribution. However, genome-wide identification and functional analysis of the ADF family in wheat is limited. Further, because of the promising role of ADF genes in cold response, there is need for an understanding of the function of this family on wheat under cold stress. In this study, 25 ADF genes (TaADFs) were identified in the wheat genome and they are distributed on 15 chromosomes. The TaADF gene structures, duplication events, encoded conversed motifs, and cis-acting elements were investigated. Expression profiles derived from RNA-seq data and real-time quantitative PCR analysis revealed the tissue- and temporal-specific TaADF expression patterns. In addition, the expression levels of TaADF13/16/17/18/20/21/22 were significantly affected by cold acclimation or freezing conditions. Overexpression of TaADF16 increased the freezing tolerance of transgenic Arabidopsis, possibly because of enhanced ROS scavenging and changes to the osmotic regulation in cells. The expression levels of seven cold-responsive genes were up-regulated in the transgenic Arabidopsis plants, regardless of whether the plants were exposed to low temperature. These findings provide fundamental information about the wheat ADF genes and may help to elucidate the regulatory effects of the encoded proteins on plant development and responses to low-temperature stress.

Past and Future of Plant Stress Detection: An Overview From Remote Sensing to Positron Emission Tomography

Plant stress detection is considered one of the most critical areas for the improvement of crop yield in the compelling worldwide scenario, dictated by both the climate change and the geopolitical consequences of the Covid-19 epidemics. A complicated interconnection of biotic and abiotic stressors affect plant growth, including water, salt, temperature, light exposure, nutrients availability, agrochemicals, air and soil pollutants, pests and diseases. In facing this extended panorama, the technology choice is manifold. On the one hand, quantitative methods, such as metabolomics, provide very sensitive indicators of most of the stressors, with the drawback of a disruptive approach, which prevents follow up and dynamical studies. On the other hand qualitative methods, such as fluorescence, thermography and VIS/NIR reflectance, provide a non-disruptive view of the action of the stressors in plants, even across large fields, with the drawback of a poor accuracy. When looking at the spatial scale, the effect of stress may imply modifications from DNA level (nanometers) up to cell (micrometers), full plant (millimeters to meters), and entire field (kilometers). While quantitative techniques are sensitive to the smallest scales, only qualitative approaches can be used for the larger ones. Emerging technologies from nuclear and medical physics, such as computed tomography, magnetic resonance imaging and positron emission tomography, are expected to bridge the gap of quantitative non-disruptive morphologic and functional measurements at larger scale. In this review we analyze the landscape of the different technologies nowadays available, showing the benefits of each approach in plant stress detection, with a particular focus on the gaps, which will be filled in the nearby future by the emerging nuclear physics approaches to agriculture.

Integrated Analysis of Metabolome and Transcriptome Reveals Insights for Cold Tolerance in Rapeseed (Brassica napus L.)

Rapeseed (Brassica napus L.) is an important oilseed crop in the world. Its productivity is significantly influenced by numerous abiotic stresses, including cold stress (CS). Consequently, enhancement in CS tolerance is becoming an important area for agricultural investigation and crop improvement. Therefore, the current study aimed to identify the stress-responsive genes, metabolites, and metabolic pathways based on a combined transcriptome and metabolome analysis to understand the CS responses and tolerance mechanisms in the cold-tolerant (C18) and cold-sensitive (C6) rapeseed varieties. Based on the metabolome analysis, 31 differentially accumulated metabolites (DAMs) were identified between different comparisons of both varieties at the same time points. From the transcriptome analysis, 2,845, 3,358, and 2,819 differentially expressed genes (DEGs) were detected from the comparison of C6-0 vs. C18-0, C6-1 vs. C18-1, and C6-7 vs. C18-7. By combining the transcriptome and metabolome data sets, we found that numerous DAMs were strongly correlated with several differentially expressed genes (DEGs). A functional enrichment analysis of the DAMs and the correlated DEGs specified that most DEGs and DAMs were mainly enriched in diverse carbohydrates and amino acid metabolisms. Among them, starch and sucrose metabolism and phenylalanine metabolism were significantly enriched and played a vital role in the CS adaption of rapeseed. Six candidate genes were selected from the two pathways for controlling the adaption to low temperature. In a further validation, the T-DNA insertion mutants of their Arabidopsis homologous, including 4cl3, cel5, fruct4, ugp1, axs1, and bam2/9, were characterized and six lines differed significantly in levels of freezing tolerance. The outcome of the current study provided new prospects for the understanding of the molecular basis of CS responses and tolerance mechanisms in rapeseed and present a set of candidate genes for use in improving CS adaptability in the same plant.

Impacts of Nitrogen Deficiency on Wheat (Triticum aestivum L.) Grain During the Medium Filling Stage: Transcriptomic and Metabolomic Comparisons

Nitrogen (N) supplementation is essential to the yield and quality of bread wheat (Triticum aestivum L.). The impact of N-deficiency on wheat at the seedling stage has been previously reported, but the impact of distinct N regimes applied at the seedling stage with continuous application on filling and maturing wheat grains is lesser known, despite the filling stage being critical for final grain yield and flour quality. Here, we compared phenotype characteristics such as grain yield, grain protein and sugar quality, plant growth, leaf photosynthesis of wheat under N-deficient and N-sufficient conditions imposed prior to sowing (120 kg/hm²) and in the jointing stage (120 kg/hm²), and then evaluated the effects of this continued stress through RNA-seq and GC-MS metabolomics profiling of grain at the mid-filling stage. The results showed that except for an increase in grain size and weight, and in the content of total sugar, starch, and fiber in bran fraction and white flour, the other metrics were all decreased under N-deficiency conditions. A total of 761 differentially expressed genes (DEGs) and 77 differentially accumulated metabolites (DAMs) were identified. Under N-deficiency, 51 down-regulated DEGs were involved in the process of impeding chlorophyll synthesis, chloroplast development, light harvesting, and electron transfer functions of photosystem, which resulted in the SPAD and Pn value decreased by 32 and 15.2% compared with N-sufficiency, inhibited photosynthesis. Twenty-four DEGs implicated the inhibition of amino acids synthesis and protein transport, in agreement with a 17-42% reduction in ornithine, cysteine, aspartate, and tyrosine from metabolome, and an 18.6% reduction in grain protein content. However, 14 DEGs were implicated in promoting sugar accumulation in the cell wall and another six DEGs also enhanced cell wall synthesis, which significantly increased fiber content in the endosperm and likely contributed to increasing the thousands-grain weight (TGW). Moreover, RNA-seq profiling suggested that wheat grain can improve the capacity of DNA repair, iron uptake, disease and abiotic stress resistance, and oxidative stress scavenging through increasing the content levels of anthocyanin, flavonoid, GABA, galactose, and glucose under N-deficiency condition. This study identified candidate genes and metabolites related to low N adaption and tolerance that may provide new insights into a comprehensive understanding of the genotype-specific differences in performance under N-deficiency conditions.

Integrated Transcriptomics and Metabolomics Analysis Reveal Key Metabolism Pathways Contributing to Cold Tolerance in Peanut

Low temperature (non-freezing) is one of the major limiting factors in peanut (Arachis hypogaea L.) growth, yield, and geographic distribution. Due to the complexity of cold-resistance trait in peanut, the molecular mechanism of cold tolerance and related gene networks were largely unknown. In this study, metabolomic analysis of two peanut cultivars subjected to chilling stress obtained a set of cold-responsive metabolites, including several carbohydrates and polyamines. These substances showed a higher accumulation pattern in cold-tolerant variety SLH than cold-susceptible variety ZH12 under cold stress, indicating their importance in protecting peanut from chilling injuries. In addition, 3,620 cold tolerance genes (CTGs) were identified by transcriptome sequencing, and the CTGs were most significantly enriched in the "phenylpropanoid biosynthesis" pathway. Two vital modules and several novel hub genes were obtained by weighted gene co-expression network analysis (WGCNA). Several key genes involved in soluble sugar, polyamine, and G-lignin biosynthetic pathways were substantially higher and/or responded more quickly in SLH (cold tolerant) than ZH12 (cold susceptible) under low temperature, suggesting they might be crucial contributors during the adaptation of peanut to low temperature. These findings will not only provide valuable resources for study of cold resistance in peanut but also lay a foundation for genetic modification of cold regulators to enhance stress tolerance in crops.

Abscisic acid-mimicking ligand AMF4 induced cold tolerance in wheat by altering the activities of key carbohydrate metabolism enzymes

AMF4, a recently synthetic ABA-mimicking ligand, has been reported to have long-lasting effects in inducing the expression of stress-responsive genes, hence conferring abiotic tolerance in plants. To test the effect of AMF4 on cold tolerance induction, the wheat plants were firstly foliar sprayed with AMF4 (10 μM), then the AMF4 treated and the control plants were exposed to a 24-h low temperature treatment (2/0 °C). Under low temperature stress, the AMF plants possessed significantly higher relative water content, membrane stability index and ATP concentration in leaf, while lower leaf H₂O₂ concentration, compared with the control plants. The AMF4 treatment significantly increased the activities of APX, Ca²⁺-ATPase in the chloroplasts, while decreased SOD activity under low temperature, in relation to the control plants. In addition, the AMF plants showed significantly higher activities of key carbohydrate metabolism enzymes, such as UDP-glucose pyrophorylase, hexokinase, fructokinase, ADP-Glucose pyrophosphorylase, phosphoglucomutase, glucose-6-phosphate dehydrogenase and phosphoglucoisomerase, in relation to the control plants under low temperature. These results demonstrate that AMF4 could be used to induce cold tolerance in wheat, and provide novel insights into the potential ways to enhance abiotic stress tolerance using the synthetic ABA-mimicking ligands.

Multi-Omics Analyses Reveal the Molecular Mechanisms Underlying the Adaptation of Wheat (Triticum aestivum L.) to Potassium Deprivation

Potassium (K) is essential for regulating plant growth and mediating abiotic stress responses. Elucidating the biological mechanism underlying plant responses to K-deficiency is crucial for breeding new cultivars with improved K uptake and K utilization efficiency. In this study, we evaluated the extent of the genetic variation among 543 wheat accessions differing in K-deficiency tolerance at the seedling and adult plant stages. Two accessions, KN9204 and BN207, were identified as extremely tolerant and sensitive to K-deficiency, respectively. The accessions were exposed to normal and K-deficient conditions, after which their roots underwent ionomic, transcriptomic, and metabolomic analyses. Under K-deficient conditions, KN9204 exhibited stronger root growth and maintained higher K concentrations than BN207. Moreover, 19,440 transcripts and 162 metabolites were differentially abundant in the roots of both accessions according to transcriptomic and metabolomic analyses. An integrated analysis of gene expression and metabolite profiles revealed that substantially more genes, including those related to ion homeostasis, cellular reactive oxygen species homeostasis, and the glutamate metabolic pathway, were up-regulated in KN9204 than in BN207 in response to low-K stress. Accordingly, these candidate genes have unique regulatory roles affecting plant K-starvation tolerance. These findings may be useful for further clarifying the molecular changes underlying wheat root adaptations to K deprivation.