Ancient rice cultivar extensively replaces phospholipids with non-phosphorus glycolipid under phosphorus deficiency

Combined Analysis of the Leaf Metabolome, Lipidome, and Candidate Gene Function: Insights into Genotypic Variation in Phosphorus Utilization Efficiency in Stylosanthes guianensis

Stylo (Stylosanthes guianensis) exhibits excellent tolerance to low-phosphate (Pi) availability, but the underlying mechanisms responsible for improving the phosphorus (P) utilization efficiency (PUE) remain unclear. This study employed metabolomics, lipidomics, and gene expression analyses to investigate the differential responses to low-Pi stress between the high-PUE genotype CF047827 and the cultivar Reyan No. 2. Results showed that CF047827 had higher expression levels of membrane lipid remodeling-related genes in its leaves compared to Reyan No. 2 under low-Pi conditions. This was accompanied by greater phospholipid degradation and non-P-containing lipid biosynthesis in the leaves of CF047827. Furthermore, the purple acid phosphatase gene SgPAP27a, which is more highly expressed in the leaves of CF047827 than in Reyan No. 2 under low-Pi conditions, was identified and functionally characterized. Its role in promoting phospholipid degradation and enhancing PUE was confirmed through heterologous expression in Arabidopsis. These findings provide insights and identify potential candidate genes for breeding high-PUE crop cultivars.

Lipidome Profiling of Phosphorus Deficiency-Tolerant Rice Cultivars Reveals Remodeling of Membrane Lipids as a Mechanism of Low P Tolerance

Plants have evolved various mechanisms for low P tolerance, one of which is changing their membrane lipid composition by remodeling phospholipids with non-phospholipids. The objective of this study was to investigate the remodeling of membrane lipids among rice cultivars under P deficiency. Rice (Oryza sativa L.) cultivars (Akamai, Kiyonishiki, Akitakomachi, Norin No. 1, Hiyadateine, Koshihikari, and Netaro) were grown in 0 (-P) and 8 (+P) mg P L^-1 solution cultures. Shoots and roots were collected 5 and 10 days after transplanting (DAT) in solution culture and subjected to lipidome profiling using liquid chromatography-mass spectrometry. Phosphatidylcholine (PC)34, PC36, phosphatidylethanolamine (PE)34, PE36, phosphatidylglycerol (PG)34, phosphatidylinositol (PI)34 were the major phospholipids and digalactosyldiacylglycerol (DGDG)34, DGDG36, 1,2-diacyl-3-O-alpha-glucuronosylglycerol (GlcADG)34, GlcADG36, monogalactosyldiacylglycerol (MGDG)34, MGDG36, sulfoquinovosyldiacylglycerol (SQDG)34 and SQDG36 were the major non-phospholipids. Phospholipids were lower in the plants that were grown under -P conditions than that in the plants that were grown under +P for all cultivars at 5 and 10 DAT. The levels of non-phospholipids were higher in -P plants than that in +P plants of all cultivars at 5 and 10 DAT. Decomposition of phospholipids in roots at 5 DAT correlated with low P tolerance. These results suggest that rice cultivars remodel membrane lipids under P deficiency, and the ability of remodeling partly contributes to low P tolerance.

Phosphorus fractions in leaves

Leaf phosphorus (P) comprises four major fractions: inorganic phosphate (P_i ), nucleic acids, phospholipids, P-containing metabolites and a residual fraction. In this review paper, we investigated whether allocation of P fractions varies among groups of terrestrial vascular plants, and is indicative of a species' strategy to use P efficiently. We found that as leaf total P concentration increases, the P_i fraction increases the most, without a plateau, while other fractions plateau. Variability of the concentrations of leaf P fractions is greatest among families > species(family) > regions > plant life forms. The percentage of total P allocated to nucleic acid-P (20-35%) and lipid-P (14-34%) varies less among families/species. High photosynthetic P-use efficiency is associated with low concentrations of all P fractions, and preferential allocation of P to metabolite-P and mesophyll cells. Sequential resorption of P from senescing leaves starts with P_i , followed by metabolite-P, and then other organic P fractions. Allocation of P to leaf P fractions varies with season. Leaf phytate concentrations vary considerably among species, associated with variation in photosynthesis and defence. Plasticity of P allocation to its fractions is important for acclimation to low soil P availability, and species-specific P allocation is needed for co-occurrence with other species.

Adaptive strategies of plants to conserve internal phosphorus under P deficient condition to improve P utilization efficiency

Phosphorus (P) is one of the limiting factors for plant growth and productivity due to its slow diffusion and immobilization in the soil which necessitates application of phosphatic fertilizers to meet the crop demand and obtain maximum yields. However, plants have evolved mechanisms to adapt to low P stress conditions either by increasing acquisition (alteration of belowground processes) or by internal inorganic P (Pi) utilization (cellular Pi homeostasis) or both. In this review, we have discussed the adaptive strategies that conserve the use of P and maintain cellular Pi homeostasis in the cytoplasm. These strategies involve modification in membrane lipid composition, flavanol/anthocyanin level, scavenging and reutilization of Pi adsorbed in cell wall pectin, remobilization of Pi during senescence by enzymes like RNases and purple acid phosphatases, alternative mitochondrial electron transport, and glycolytic pathways. The remobilization of Pi from senescing tissues and its internal redistribution to various cellular organelles is mediated by various Pi transporters. Although much efforts have been made to enhance P acquisition efficiency, an understanding of the physiological mechanisms conserving internal Pi and their manipulation would be useful for plants that can utilize P more efficiently to produce optimum growth per unit P uptake.

Physiological and Transcriptomic Analyses Uncover the Reason for the Inhibition of Photosynthesis by Phosphate Deficiency in Cucumis melo L

Phosphate (Pi) deficiency is a common phenomenon in agricultural production and limits plant growth. Recent work showed that long-term Pi deficiency caused the inhibition of photosynthesis and inefficient electron transport. However, the underlying mechanisms are still unknown. In this study, we used the physiological, histochemical, and transcriptomic methods to investigate the effect of low-Pi stress on photosynthetic gas exchange parameters, cell membrane lipid, chloroplast ultrastructure, and transcriptional regulation of key genes in melon seedlings. The results showed that Pi deficiency significantly downregulated the expression of aquaporin genes, induced an increase in ABA levels, and reduced the water content and free water content of melon leaves, which caused physiological drought in melon leaves. Therefore, gas exchange was disturbed. Pi deficiency also reduced the phospholipid contents in leaf cell membranes, caused the peroxidation of membrane lipids, and destroyed the ultrastructure of chloroplasts. The transcriptomic analysis showed that 822 differentially expressed genes (DEGs) were upregulated and 1254 downregulated by Pi deficiency in leaves. GO and KEGG enrichment analysis showed that DEGs significantly enriched in chloroplast thylakoid membrane composition (GO:0009535), photosynthesis-antenna proteins (map00196), and photosynthesis pathways (map00195) were downregulated by Pi deficiency. It indicated that Pi deficiency regulated photosynthesis-related genes at the transcriptional level, thereby affecting the histochemical properties and physiological functions, and consequently causing the reduced light assimilation ability and photosynthesis efficiency. It enriches the mechanism of photosynthesis inhibition by Pi deficiency.

Combining analyses of metabolite profiles and phosphorus fractions to explore high phosphorus utilization efficiency in maize

Phosphorus (P) limitation is a significant factor restricting crop production in agricultural systems, and enhancing the internal P utilization efficiency (PUE) of crops plays an important role in ensuring sustainable P use in agriculture. To better understand how P is remobilized to affect crop growth, we first screened P-efficient (B73 and GEMS50) and P-inefficient (Liao5114) maize genotypes at the same shoot P content, and then analyzed P pools and performed non-targeted metabolomic analyses to explore changes in cellular P fractions and metabolites in maize genotypes with contrasting PUE. We show that lipid P and nucleic acid P concentrations were significantly lower in lower leaves of P-efficient genotypes, and these P pools were remobilized to a major extent in P-efficient genotypes. Broad metabolic alterations were evident in leaves of P-efficient maize genotypes, particularly affecting products of phospholipid turnover and phosphorylated compounds, and the shikimate biosynthesis pathway. Taken together, our results suggest that P-efficient genotypes have a high capacity to remobilize lipid P and nucleic acid P and promote the shikimate pathway towards efficient P utilization in maize.

Leaf phosphorus fractionation in rice to understand internal phosphorus-use efficiency

BACKGROUND AND AIMS

Phosphorus (P) availability is often limiting for rice (Oryza sativa) production. Improving internal P-use efficiency (PUE) is crucial to sustainable food production, particularly in low-input systems. A critical aspect of PUE in plants, and one that remains poorly understood, is the investment of leaf P in different chemical P fractions (nucleic acid-P, lipid-P, inorganic-P, metabolite-P and residual-P). The overarching objective of this study was to understand how these key P fractions influence PUE.

METHODS

Three high-PUE and two low-PUE rice genotypes were grown in hydroponics with contrasting P supplies. We measured PUE, total P, P fractions, photosynthesis and biomass.

KEY RESULTS

Low investment in lipid-P was strongly associated with increased photosynthetic PUE (PPUE), achieved by reducing total leaf P concentration while maintaining rapid photosynthetic rates. All low-P plants exhibited a low investment in inorganic-P and lipid-P, but not nucleic acid-P. In addition, whole-plant PUE was strongly associated with reduced total P concentration, increased biomass and increased preferential allocation of resources to the youngest mature leaves.

CONCLUSIONS

Lipid remodelling has been shown in rice before, but we show for the first time that reduced lipid-P investment improves PUE in rice without reducing photosynthesis. This presents a novel pathway for increasing PUE by targeting varieties with reduced lipid-P investment. This will benefit rice production in low-P soils and in areas where fertilizer use is limited, improving global food security by reducing P fertilizer demands and food production costs.

Physiology of microalgae and their application to sustainable agriculture: A mini-review

Concern that depletion of fertilizer feedstocks, which are a finite mineral resource, threatens agricultural sustainability has driven the exploration of sustainable methods of soil fertilization. Given that microalgae, which are unicellular photosynthetic organisms, can take up nutrients efficiently from water systems, their application in a biological wastewater purification system followed by the use of their biomass as a fertilizer alternative has attracted attention. Such applications of microalgae would contribute to the accelerated recycling of nutrients from wastewater to farmland. Many previous reports have provided information on the physiological characteristics of microalgae that support their utility. In this review, we focus on recent achievements of studies on microalgal physiology and relevant applications and outline the prospects for the contribution of microalgae to the establishment of sustainable agricultural practices.

Specific galactolipids species correlate with rice genotypic variability for phosphate utilization efficiency

Membrane lipid remodeling helps in the efficient utilization of phosphorus (P) by replacing phospholipids with galactolipids during P deficiency. Previous studies have shown lipid remodeling in rice under P deficiency; however, main lipid classes did not show association with superior P-use-efficiency in rice genotypes. Here, diverse rice genotypes were extensively phenotyped in normal (NP) and low P (LP) conditions. Based on the phenotypic response to P deficiency, genotypes were identified as tolerant and sensitive. Further, bulks were generated differing in their physiological P-use-efficiency (PPUE) during LP condition. Shoot lipidome profiling of genotypes was performed and used to correlate the abundance of various lipid classes and their constituent species with the PPUE of the genotypes. Lipid remodeling was observed as a P-starvation-induced response in all the genotypes. However, neither total galacto- and phospholipids nor the lipid classes correlated with PPUE during P deficiency. However, the difference in PPUE in the two bulks correlated with specific lipid species of galactolipids (DGDG, MGDG). Further, DGDG34:3 had the highest Mol% among the differentially accumulated lipid species between the two bulks. Our study reveals the importance of specific galactolipids species in rice adaptation to P deficient soils and thus opens new targets for future research.

Multiple analysis of root exudates and microbiome in rice (Oryza sativa) under low P conditions

Plants release various metabolites from roots and root exudates contribute to differences in stress tolerance among plant species. Plant and soil microbes have complex interactions that are affected by biotic and abiotic factors. The purpose of this study was to examine the differences in metabolites in root exudates of rice (Oryza sativa) cultivars and their correlation with bacterial populations in the rhizosphere. Two rice cultivars (O. sativa cv. Akamai and O. sativa cv. Koshihikari) were grown in soils fertilized with 0 g P kg^-1 (- P) or 4.8 g P kg^-1 (+ P). Root exudates and root-attached soil were collected at 13 and 20 days after transplanting (DAT) and their metabolites and bacterial community structure were determined. The exudation of proline, serine, threonine, valine and 4-coumarate were increased under low P conditions in both cultivars. There was a positive correlation between the concentration of pantothenate in root exudates and the representation of members of the genera Clostridium and Sporosarcina, which were negatively correlated with root dry weight. Gracilibacter, Opitutus, Pelotomaculum, Phenylobacterium and Oxobacter were positively correlated with root dry weight and presence of allantoin, 2-aminobtyrate and GlcNac. This study provides new information about the response of plants and rhizosphere soil bacteria to low P conditions.

Characterization of Glucuronosyl-diacyl/monoacylglycerols and Discovery of Their Acylated Derivatives in Tomato Lipid Extracts by Reversed-Phase Liquid Chromatography with Electrospray Ionization and Tandem Mass Spectrometry

Glucuronic acid containing diacylglycerols (3-(O-α-d-glucuronopyranosyl)-1,2-diacyl-sn-glycerols, GlcA-DAG) are glycolipids of plant membranes especially formed under phosphate-depletion conditions. An analytical approach for the structural characterization of GlcA-DAG in red ripe tomato (Solanum lycopersicum L.) extracts, based on reversed-phase liquid chromatography (RPLC) coupled with electrospray ionization (ESI) and tandem mass spectrometry (MS/MS) using a linear ion trap, is described in this paper. At least 14 GlcA-DAG (R₁/R₂) species, including four regioisomers, containing three predominant fatty acyl chains C16:0, C18:2, and C18:3, were identified for the first time. Moreover, 29 GlcA-DAG acylated on the glucuronosyl ring (acyl-R₃ GlcA-DAG) were discovered, alongside 15 acylated lyso-forms, i.e., acylated 3-(O-α-d-glucuronosyl)monoacylglycerols, abbreviated as acyl-R₃ GlcA-MAG (R₁/0) or (0/R₂). Although many of these acylated lyso-forms were isomeric with GlcA-DAG (i.e., acyl chains with equivalent sum composition), they were successfully separated by reversed-phase liquid chromatography (RPLC) using a solid-core C18 column packed with 2.6 μm particle size. Tandem MS (and eventually MS³) data obtained from sodium adducts ([M + Na]⁺) and deprotonated molecules ([M - H]^-) were fundamental to detect diagnostic product ions related to the glucuronosyl ring and then determine the identity of all investigated glycolipids, especially to recognize the acyl chain linked to the ring. A classification of GlcA-MAG, GlcA-DAG, and acylated GlcA-DAG and GlcA-MAG was generated by an in house-built database. The discovery of acylated derivatives emphasized the already surprising heterogeneity of glucuronic acid-containing mono- and diacylglycerols in tomato plants, stimulating interesting questions on the role played by these glycolipids.

Impact of Plant Growth-Promoting Rhizobacteria Inoculation and Grafting on Tolerance of Tomato to Combined Water and Nutrient Stress Assessed via Metabolomics Analysis

In the current study, inoculation with plant growth-promoting rhizobacteria (PGPR) and grafting were tested as possible cultural practices that may enhance resilience of tomato to stress induced by combined water and nutrient shortage. The roots of tomato grown on perlite were either inoculated or not with PGPR, applying four different treatments. These were PGPR-T1, a mix of two Enterobacter sp. strains (C1.2 and C1.5); PGPR-T2, Paenibacillus sp. strain DN1.2; PGPR-T3, Enterobacter mori strain C3.1; and PGPR-T4, Lelliottia sp. strain D2.4. PGPR-treated plants were either self-grafted or grafted onto Solanum lycopersicum cv. M82 and received either full or 50% of their standard water, nitrogen, and phosphorus needs. The vegetative biomass of plants subjected to PGPR-T1 was not reduced when plants were cultivated under combined stress, while it was reduced by stress to the rest of the PGPR treatments. However, PGPR-T3 increased considerably plant biomass of non-stressed tomato plants than did all other treatments. PGPR application had no impact on fruit biomass, while grafting onto 'M82' increased fruit production than did self-grafting. Metabolomics analysis in tomato leaves revealed that combined stress affects several metabolites, most of them already described as stress-related, including trehalose, myo-inositol, and monopalmitin. PGPR inoculation with E. mori strain C3.1 affected metabolites, which are important for plant/microbe symbiosis (myo-inositol and monopalmitin). The rootstock M82 did not affect many metabolites in plant leaves, but it clearly decreased the levels of malate and D-fructose and imposed an accumulation of oleic acid. In conclusion, PGPR are capable of increasing tomato tolerance to combined stress. However, further research is required to evaluate more strains and refine protocols for their application. Metabolites that were discovered as biomarkers could be used to accelerate the screening process for traits such as stress tolerance to abiotic and/or abiotic stresses. Finally, 'M82' is a suitable rootstock for tomato, as it is capable of increasing fruit biomass production.

Effect of phosphorus supply on root traits of two Brassica oleracea L. genotypes

BACKGROUND

Phosphorus (P) deficiency limits crop production worldwide. Crops differ in their ability to acquire and utilise the P available. The aim of this study was to determine root traits (root exudates, root system architecture (RSA), tissue-specific allocation of P, and gene expression in roots) that (a) play a role in P-use efficiency and (b) contribute to large shoot zinc (Zn) concentration in Brassica oleracea.

RESULTS

Two B. oleracea accessions (var. sabellica C6, a kale, and var. italica F103, a broccoli) were grown in a hydroponic system or in a high-throughput-root phenotyping (HTRP) system where they received Low P (0.025 mM) or High P (0.25 mM) supply for 2 weeks. In hydroponics, root and shoot P and Zn concentrations were measured, root exudates were profiled using both Fourier-Transform-Infrared spectroscopy and gas-chromatography-mass spectrometry and previously published RNAseq data from roots was re-examined. In HTRP experiments, RSA (main and lateral root number and lateral root length) was assessed and the tissue-specific distribution of P was determined using micro-particle-induced-X-ray emission. The C6 accession had greater root and shoot biomass than the F103 accession, but the latter had a larger shoot P concentration than the C6 accession, regardless of the P supply in the hydroponic system. The F103 accession had a larger shoot Zn concentration than the C6 accession in the High P treatment. Although the F103 accession had a larger number of lateral roots, which were also longer than in the C6 accession, the C6 accession released a larger quantity and number of polar compounds than the F103 accession. A larger number of P-responsive genes were found in the Low P treatment in roots of the F103 accession than in roots of the C6 accession. Expression of genes linked with "phosphate starvation" was up-regulated, while those linked with iron homeostasis were down-regulated in the Low P treatment.

CONCLUSIONS

The results illustrate large within-species variability in root acclimatory responses to P supply in the composition of root exudates, RSA and gene expression, but not in P distribution in root cross sections, enabling P sufficiency in the two B. oleracea accessions studied.

Insights Into Oxidized Lipid Modification in Barley Roots as an Adaptation Mechanism to Salinity Stress

Lipidomics is an emerging technology, which aims at the global characterization and quantification of lipids within biological matrices including biofluids, cells, whole organs and tissues. The changes in individual lipid molecular species in stress treated plant species and different cultivars can indicate the functions of genes affecting lipid metabolism or lipid signaling. Mass spectrometry-based lipid profiling has been used to track the changes of lipid levels and related metabolites in response to salinity stress. We have developed a comprehensive lipidomics platform for the identification and direct qualification and/or quantification of individual lipid species, including oxidized lipids, which enables a more systematic investigation of peroxidation of individual lipid species in barley roots under salinity stress. This new lipidomics approach has improved with an advantage of analyzing the composition of acyl chains at the molecular level, which facilitates to profile precisely the 18:3-containing diacyl-glycerophosphates and allowed individual comparison of lipids across varieties. Our findings revealed a general decrease in most of the galactolipids in plastid membranes, and an increase of glycerophospholipids and acylated steryl glycosides, which indicate that plastidial and extraplastidial membranes in barley roots ubiquitously tend to form a hexagonal II (HII) phase under salinity stress. In addition, salt-tolerant and salt-sensitive cultivars showed contrasting changes in the levels of oxidized membrane lipids. These results support the hypothesis that salt-induced oxidative damage to membrane lipids can be used as an indication of salt stress tolerance in barley.

Early Responses of Brassica oleracea Roots to Zinc Supply Under Sufficient and Sub-Optimal Phosphorus Supply

Shoot zinc (Zn) concentration in Brassica oleracea is affected by soil Zn and phosphorus (P) supply. Most problematic is the negative impact of P fertilizers on Zn concentrations in crops, which makes balancing yield and mineral quality challenging. To evaluate early molecular mechanisms involved in the accumulation of large shoot Zn concentrations regardless of the P supply, two B. oleracea accessions differing in root architecture and root exudates were grown hydroponically for two weeks with different combinations of P and Zn supply. Ionome profiling and deep RNA sequencing of roots revealed interactions of P and Zn in planta, without apparent phenotypic effects. In addition, increasing P supply did not reduce tissue Zn concentration. Substantial changes in gene expression in response to different P and/or Zn supplies in roots of both accessions ensured nutritionally sufficient P and Zn uptake. Numerous genes were differentially expressed after changing Zn or P supply and most of them were unique to only one accession, highlighting their different strategies in achieving nutrient sufficiency. Thus, different gene networks responded to the changing P and Zn supply in the two accessions. Additionally, enrichment analysis of gene ontology classes revealed that genes involved in lipid metabolism, response to starvation, and anion transport mechanisms were most responsive to differences in P and Zn supply in both accessions. The results agreed with previously studies demonstrating alterations in P and Zn transport and phospholipid metabolism in response to reduced P and Zn supply. It is anticipated that improved knowledge of genes responsive to P or Zn supply will help illuminate the roles in uptake and accumulation of P and Zn and might identify candidate genes for breeding high-yield-high-Zn brassicas.

Adaptation of Synechococcus sp. PCC 7942 to phosphate starvation by glycolipid accumulation and membrane lipid remodeling

Organisms use various adaptive strategies against phosphate stress, including lipid remodeling. Here, the response of major membrane lipids to phosphate stress was analyzed in Synechococcus sp. PCC 7942. Unlike plants and eukaryotic microalgae, no significant increases in neutral lipids were found, whereas glycolipids content increased to as high as 6.13% (of dry cell weight, DCW) and phospholipids decreased to 0.34% (of DCW) after 16 days of cultivation without phosphate. Glycolipids accumulation were mainly attributed to the significant increase of digalactosyldiacylglycerol (DGDG) by 50% and sulfoquinovosyldiaclglycerol (SQDG) by 90%, both of which acted as complementary lipids for phosphatidylglycerol (PG) in the cyanobacterial membrane. Also, a notable increase in content (by 48%) of C18 fatty acids (especially C18:1) was observed in all glycolipids at the expense of C12 and C14 (72%). These changes may contribute to membrane fluidity and photosynthetic activity for basic cell metabolism and phosphate stress adaptation. Lipidomic analyses showed the reduction of PG 18:1/16: 0 (by 52%) with the increase of DGDG 18:1/16:0 (133%) and SQDG 18:1/16:0 (245%), strongly suggesting a direct conversion of PG to DGDG and SQDG. Moreover, the decreasing amount of monogalactosyldiacylglycerol (MGDG) 16:1/16:0 (22%) was consistent with the increase of free fatty acids (125%) on day 2 of phosphate absence, which suggested that MGDG is more likely to provide a pool of fatty acids for de novo synthesis of glycolipids. This study provides valuable insight into cyanobacteria adaptation strategies to phosphate stress by membrane lipid remodeling and unveils the underlying acyl chain fluxes into glycolipids.

Molecular mechanisms underpinning phosphorus-use efficiency in rice

Orthophosphate (H₂ PO₄^- , Pi) is an essential macronutrient integral to energy metabolism as well as a component of membrane lipids, nucleic acids, including ribosomal RNA, and therefore essential for protein synthesis. The Pi concentration in the solution of most soils worldwide is usually far too low for maximum growth of crops, including rice. This has prompted the massive use of inefficient, polluting, and nonrenewable phosphorus (P) fertilizers in agriculture. We urgently need alternative and more sustainable approaches to decrease agriculture's dependence on Pi fertilizers. These include manipulating crops by (a) enhancing the ability of their roots to acquire limiting Pi from the soil (i.e. increased P-acquisition efficiency) and/or (b) increasing the total biomass/yield produced per molecule of Pi acquired from the soil (i.e. increased P-use efficiency). Improved P-use efficiency may be achieved by producing high-yielding plants with lower P concentrations or by improving the remobilization of acquired P within the plant so as to maximize growth and biomass allocation to developing organs. Membrane lipid remodelling coupled with hydrolysis of RNA and smaller P-esters in senescing organs fuels P remobilization in rice, the world's most important cereal crop.