Changes in the Transcriptomic Profiles of Allohexaploid Wheat in Response to Low Fe Stress

doi:10.21203/rs.3.rs-1249808/v1

Download PDF

Research Article

Changes in the Transcriptomic Profiles of Allohexaploid Wheat in Response to Low Fe Stress

https://doi.org/10.21203/rs.3.rs-1249808/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Background: Around the world, especially in aerobic soils, because the solubility of iron (Fe) is very low, plants are often stressed by low Fe tress, and the growth, yield, and quality of crops will be inhibited in the case of severe Fe deficiency. The metabolism of Fe in plants is controlled by a series of complex transport, storage, and regulatory mechanisms to maintain the homeostasis of Fe in cells. Allohexaploid wheat (Triticum aestivum L.) is an important food crop and is more sensitive to low Fe tress. Although some studies have been conducted on the low Fe tress response of different plant species, these mechanisms are still unclear in wheat.

Results: According to the results of transmission electron microscope, and paraffin section, under low Fe stress, leaf chlorosis was due to damage to the chloroplast structure, and the effect on root structure was mediated by reducing the rate of cell division in meristems and reducing cell elongation in the elongation zone. ICP-MS showed that low Fe stress significantly limited the absorption of essential elements, including N, Pi, K, Ca, Mg, Fe, Mn, Cu, Zn, and B nutrients. RNA sequencing revealed the transcriptomic changes of wheat under low Fe stress. The results showed that 378 and 2,619 differentially expressed genes (DEGs) were identified in the shoots and roots, respectively. These DEGs were mainly involved in the synthesis of Fe chelating agents, ion transport, photosynthesis, amino acid metabolism, protein synthesis, and other processes. Next, to find the core genes responding to low Fe stress, the gene co-expression network diagram was constructed. The results indicated that TaIRT1b-4A, TaNAS2-6D, TaNAS1a-6A, TaNAS1-6B, and TaNAAT1b-1D might play a key role in dealing with low-Fe stress.

Conclusions: These research results might help to fully understand the morphological and molecular responses of plants to low Fe stress, and provide excellent genetic resources for the genetic modification of Fe deficiency crops.

Ion homeostasis

Low Fe stress

Morpho-physiologic response

Transcriptomic

Triticum aestivum L

Iron (Fe) plays an important role in maintaining the normal growth and development of higher plants, such as participating in chlorophyll synthesis, respiration, redox reactions, and electron transfer [1]. Although the Fe content in the soil is very high, it is often in the form of Fe³⁺ which is hardly soluble in water due to the influence of soil pH and oxygen partial pressure. The total amount of soluble Fe is less than 10⁻¹⁰ mol·L⁻¹ and the Fe concentration required by plants to maintain growth is about 10 ⁸ mol·L⁻¹, so it is difficult for plants to use [1, 2], especially in calcareous soils with higher pH, low Fe is more serious [3]. Studies have shown that the lack of Fe will not only affect the growth and development of plants but also affect the absorption of Fe by animals and humans [2]. According to investigations, more than 40% of the world’s soils are severely Fe deficient [4]. If there is no active regulation mechanism to allow plants to obtain sufficient Fe, most plants will show Fe-deficiency symptoms. Therefore, low Fe and chlorosis in plants have become a universal concern around the world, and it is also been one of the hot spots in plant nutrition research since the beginning of this century.

Faced with different Fe nutritional statuses in the environment, higher plants have formed two different strategies in the long-term evolution process to adapt to the changes in Fe content in the external environment: reduction strategy (strategy I) and chelation strategy (strategy II) [5]. Monocotyledonous non-graminaceous plants and all dicotyledonous plants adopt strategy I to deal with a low Fe environment. When these plants are lacking Fe, their root epidermis adopts three processes of acidification, reduction, and absorption to absorb rhizosphere Fe [3]. First, a proton ATPase (H⁺-ATPase) is expressed in large quantities, pumping a large number of protons (H⁺) to the root environment to acidify the rhizosphere environment, thereby improving the availability of Fe in the rhizosphere environment. Second, up-regulate the expression of a gene encoding Fe reductase (FRO), and then reduce the Fe³⁺ adsorbed on the root epidermis to Fe²⁺. At the same time, the iron-regulated transporter (IRT) gene encoding specific uptake of Fe²⁺ is also up-regulated to transport Fe²⁺ into root cells [3]. Strategy II adopted by monocotyledonous gramineous plants consists of two parts [6]. First, by synthesizing a plant Fe carrier (phytosiderophore, PS) that can bind to Fe³⁺ with high affinity, it is then secreted into the rhizosphere environment and combined with Fe³⁺ to form a Fe³⁺-PS complex. At the same time, under the action of the YSL (YS1-like protein) family transporter, the Fe³⁺-PS complex near the root epidermal cells is transported across the membrane to the inside of the plant, thereby completing the absorption of Fe in monocotyledonous grasses [7]. The proposal of strategy I and strategy II opened up a new field of research on the specific response of plant nutrient stress, enabled people to have a deeper understanding of the mechanism of plant low Fe tress, and also provides an important theoretical basis for solving the problem of low Fe in plants.

Triticum aestivum L. is the most important crop in the world and the second-largest crop in my country after rice. Its yield and quality directly affect the improvement of the human diet [8]. Fe, as an essential micronutrient element for Triticum aestivum, has important biological functions. Low Fe tress will cause distortion or destruction of the cell membrane system and intracellular organelles [9], and the phenotype will be chlorosis and chlorosis [10]. Low Fe tress can also induce the reduction of photosynthetic pigments in the shoots, weaken the capture and transmission of light energy, damage the photosystem II reaction center, and lead to the reduction of excitation energy capture efficiency, electron transfer efficiency, and photosynthetic efficiency. Thus inhibiting the growth and development of seedlings leads to a decrease in biomass [11, 12]. Under low Fe tress, the activity of Fe reductase in the rhizosphere is significantly increased, thereby increasing the ability of plants to reduce Fe and increasing the Fe use efficiency of Triticum aestivum roots [13]. Previous studies have shown that the Triticum aestivum genome contains two homologous genes TaVIT1 and TaVIT2 of vacuolar Fe transporter (VIT), which have different expression patterns, but the expression in the endosperm is lower. Overexpression of TaVIT2 can enhance the transport of vacuolar Fe in the endosperm, thereby increasing the transport of Fe to seeds [14]. Although many scholars have studied Triticum aestivum roots with low Fe stress, most of them mainly focus on photosynthetic characteristics, root characteristics, seed germination, and plant physiological and biochemical effects. There is still a lack of in-depth research on the adaptability of Triticum aestivum to low Fe tress. Therefore, studying the resistance mechanism of Triticum aestivum under low Fe stress has very important theoretical and practical guiding significance for improving the quality of crops.

The nutritional quality of crops, especially the effect of microelements represented by Fe on human health, has attracted more and more attention. In recent years, important progress has been made in the molecular mechanisms of Fe absorption, transport, and metabolism in higher plants. Although some key genes, such as fer-like iron deficiency-induced transcription factor (FIT), have been successfully cloned, their relationship with other regulatory factors in the molecular regulatory network has become a hot spot for future research. In addition, the lack of Fe is also a major problem faced by the growing population. Increasing the Fe content in grains has become a hot issue for enhanced biological breeding in the future. Therefore, the research on the molecular mechanism of plants' efficient absorption and utilization of Fe has become an important direction for future plant nutrition research, and at the same time, it provides a basis for the cultivation of high-quality new varieties of Fe with high-efficiency utilization.

Plant material and low Fe tress treatment

The bread wheat (Triticum aestivum L.) cultivar zhengmai 1860 cultivated by Henan Academy of Agricultural Sciences (Zhengzhou, Henan Province, China) was used in this research. We germinated the seeds of zhengmai 1860 on wet gauze for 7 days and then selected the seedlings with the same growing vigor for hydroponic cultivation in a black plastic bottle containing 10 L of Hoagland nutrient solution (pH 5.8). Half of the seedlings were cultivated in a nutrient solution containing sufficient Fe, and half were cultivated in a nutrient solution containing low Fe stress. Hydroponic solutions using basic nutrition solution contained 1.0 mM KH₂PO₄, 5.0 mM KNO₃, 5.0 mM Ca(NO₃)₂·4H₂O, 2.0 mM MgSO₄·7H₂O, 9.0 μM MnCl₂·4H₂O, 0.80 μM ZnSO₄·7H₂O, 0.30 μM CuSO₄·5H₂O, 0.10 μM Na₂MoO₄·2H₂O, and 46 μM H₃BO₃ with different Fe(III)-EDTA concentrations (Fe-sufficient condition: 50 μM; low Fe stress condition: 2 μM). The rapeseed seedlings were cultivated in an illuminated chamber with the following growth condition: the light intensity of 200 μmol m^-2 s^-1, the temperature of 25°C daytime/22°C night, the light period of 16 h photoperiod/8 h dark, the relative humidity of 70%. During hydroponic culture, the nutrient solutions were refreshed every 5 days. The pH of the nutrient solution was adjusted to 5.8 with NaOH.

Determination of physiological indexes

10 days after the treatment, selected 5 seedlings with uniform growth under different treatments, took samples from the shoots and roots, and measured their fresh weight. Then, dried at 65ºC to constant weight, weighed the dry weight, and calculated the root-to-shoot ratio and relative water content of the shoots. According to the previous research [48], the shoots were cut into silk and mixed, and 0.2 g of the samples were weighed and placed in a test tube containing 5 ml of 80% acetone solution and soaked in a dark place at room temperature for 24 h. 80% acetone solution was used as the reference solution, and the DU®640 UV-Vis spectrophotometer (Beckman) was used to measure the absorbance at wavelengths of 470 nm, 645 nm, and 663 nm to determine the chlorophyll content.

We took fresh root samples of the whole plant in 5 mL 0.4 M mannitol at 25°C, 100 r/min shaker for 4 hours, and passed the conductivity meter (DDSJ-318, Shanghai, China) to measure the conductivity of the solution and record it as the initial leakage rate. Then put the samples in a water bath at 85°C for 20 minutes, measured the leakage rate of the solution, and recorded it as the total leakage rate. Ion leakage rate = initial leakage rate/total leakage rate [49].

The root scanner (Microtek Scan Maker I800, WinRHIZO Pro) was used to measure the total root length, primary root length, root surface area, root volume, number of root tips, lateral root length, and average root diameter of each plant, and 5 biological replicates were made for each treatment. The Suzhou Mengxi Biomedical Technology Kit was used to determine the root activity of wheat under different treatments by the tetrazolium (TTC) reduction method, and 5 biological replicates were performed for the different treatments.

The Keming biological kit was used to determine the content of proline (Pro), superoxide anion (OFR), hydrogen peroxide (H₂O₂), and malondialdehyde (MDA) in the shoots and roots of wheat under different treatments, and 5 biological replicates were made for each treatment.

Microscopy analysis

10 days after treatment of the wheat plants, the roots were taken and placed in a glass container filled with ddH₂O. The morphology of the root hairs was observed with a stereoscopic microscope (Leica M125C).

Paraffin sections were used to study the cell size in wheat root tips. First, put fresh wheat root tips of about 0.5 cm into a 5 ml reagent tube, added 2-3 ml FAA fixative (formaldehyde: glacial acetic acid: 70% ethanol=1:1:18), and performed aspiration treatment on the fixed material in time. It was then embedded after dehydration, transparency, and wax dipping. When using conventional paraffin section technology to slice the fixed material, the material needed to be stained by the safranine-fast green double staining method and then sliced with a microtome (Leica RM2235). The slice thickness was 8-10 μm, and a typical section was selected for microscopy, and finally took pictures under a stereoscopic microscope (Leica M125C) [50].

The shoots of wheat under different treatments were cleaned with distilled water, the excess water was absorbed by the absorbent paper, and they were fixed in fresh 4% glutaraldehyde fixative for 12 h. Then rinsed with 0.1 mol/L phosphate-buffered solution (PBS) of pH 6.8 for 4-5 times, and put the rinsed material into 1% osmium acid solution for fixation at 4 ℃ for 1 h, and then dehydrated with 100% ethanol. Infiltration, finally resin embedding and polymerization, the tissue block was taken out and stored in a desiccator. The sample was cut into ultra-thin sections of 50-70 nm, double-stained with uranyl acetate-lead citrate, observed and photographed by transmission electron microscope (H-7650; Hitachi, Tokyo, Japan) [51, 52].

Detection of ROS in wheat plants under different treatments

ROS accumulation can be estimated based on the histochemical staining intensity of total H₂O₂ and O₂^¯ in the roots. 3’3’-Diaminobenzidine (DAB) can react with H₂O₂ in the presence of peroxidase to form a dark brown polymer product, thereby detecting the H₂O₂ content in plant tissues [53]. Nitro blue tetrazole (NBT) can react with oxygen-free radical O₂^¯ to produce dark blue insoluble formazan, which can be used to detect the O₂^¯ content in plant tissues [54]. Finally, a stereoscopic microscope (Leica M125C) was used to take pictures.

Ionoqmic analysis

The Triticum aestivum seedlings after 7-day seed germination were hydroponically grown under Fe-sufficient and low-Fe solution for 10 days until sampling. Each sample contained 5 independent biological replicates for the ionomic analyses under low Fe stress. The over-dried shoot and root tissues were subsequently transferred to an HNO₃/HClO₄ mixture (4:1, v/v) at 200℃ until the digestion was completed. The diluted supernatant was determined to quantify the concentrations of mineral elements using an inductively coupled plasma mass spectrometry (ICP-MS; NexIONTM 350X, PerkinElmer).

Transcriptome Analysis

For the low Fe treatment, the 7-d-old uniform seedlings after seed germination were hydroponically grown under Fe-sufficient and low-Fe solution for 10 days until sampling. The shoots and roots of the aforementioned Triticum aestivum plants were harvested, and three independent biological replicates were used for each treatment. Pre-chilled Trizol (Takara Bio Inc, Kusatsu, Shiga, Japan) was used to isolate total RNA, following which the RNA integrity number (RIN) was assessed. Samples of RNA with the RIN values > 8.0 were obtained to construct strand-specific cDNA libraries, which were further employed for the paired-end sequencing (read length = 150 bp) on a lane of an Illumina Hiseq 4000 platform. The fragments per kilobase of exon model per million mapped reads (FPKM) values were normalized to quantify the gene expression abundances, and both false discovery rate and P values < 0.05 were used to identify the differentially expressed genes (DEGs) [55]. PANTHER (http://www.pantherdb.org/data/) [56] and Kyoto Encyclopedia of Genes and Genomes (KEGG) (http://www.kegg.jp/) [57], respectively, were used to perform gene ontology (GO) and pathway enrichment analysis of the DEGs. Heat maps showing differential gene expression were delineated using a Multi experiment Viewer (http://www.tm4.org/mev.html) [58].

Reverse transcription-quantitative polymerase chain reaction assays

Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) assays were performed to validate the accuracy of transcriptome sequencing data. After removing genomic DNA from RNA samples with RNase-free DNase I, total RNA was used as RT templates for cDNA synthesis using the PrimeScriptTM RT Reagent Kit Eraser (Perfect Real Time; TaKaRa, Shiga, Japan). The RT-qPCR assays were performed to detect relative gene expression using SYBR®Premix Ex TaqTM II (TliRNaseH Plus) (TaKaRa, Shiga, Japan) using a Bio-Rad C1000 Touch Thermal Cycler of CFX96TM Real-time PCR detection system.

The RT-qPCR program was as follows: 95ºC for 3 min, 40 cycles of 95ºC for 10 s, and 60ºC for 30 s. The melting curve was plotted as follows to analyze the primer gene-specificity: 95ºC for 15 s, 60ºC for 1 min, and 60-95ºC for 15 s (+0.3ºC/cycle). The expression data of target genes were normalized using BnaEF1-α [59], and the relative gene expression levels were calculated according to the 2^-ΔΔCT method [60]. Each sample contained three independent biological replicates.

Gene co-expression network analysis

The gene co-expression networks were constructed to identify gene interactions and locate core genes that connect most neighboring genes involved in the response to low Fe stress in Triticum aestivum. For each pair of genes, the threshold of the Pearson correlation value was set according to the default settings (http://plantgrn.noble.org/DeGNServer/Analysis.jsp), and then gene co-expression networks were constructed by CYTOSCAPE 3.2.1 (Donnelly Centre for Cellular and Biomolecular Research, University of Toronto, Toronto, Canada) (http://www.cytoscape.org/) [61].

Statistical analysis

Statistical Productions and Service Solutions 17.0 (SPSS, Chicago, IL, USA) was used to perform statistical tests. One-way analysis of variance followed by Tukey’s honestly significant difference (P < 0.05, P < 0.01, and P < 0.001) multiple comparison tests were used to determine the significant differences.

The morphological response of wheat to low Fe stress

The typical symptom of low Fe in plants is the thin green veins on the pale green leaves. As the symptoms worsen, the leaves become thinner and white (or pale yellow), and brown spots and bad tissue appear on the leaves. After the leaves mature, the boundary between the green veins and the light green or light yellow mesophyll becomes more obvious. Our results showed that under low Fe stress, the growth and development of the root system of wheat were inhibited, and the root hair density became higher (Fig. 1A, C). Leaves were pale yellow or whitish, young leaves were chlorosis, and there was obvious visible inter-vein yellowing on the leaves (Fig. 1A, B).

The biomass of wheat was measured and found that under low Fe stress, the fresh weight and dry weight of the shoots were reduced by about 1.2 times, and the fresh weight and dry weight of the roots were reduced by about 1.3 times and 1 time, respectively (Fig. 2A, B). However, the relative water content of shoots showed no significant difference under different treatments (Figure. 2C). Determination of photosynthetic pigment content in wheat found that under low Fe stress, the concentrations of total chlorophylls (including chlorophyll a and b) and carotenoids (including xanthophyll and carotene) were significantly reduced (Fig. 2D, F). The ratio of the concentration of chlorophyll a to chlorophyll b was increased significantly under low Fe stress, which was also a sign of plant adaptation to adversity stress (Fig. 2E). In addition, low Fe stress reduced the ratio of total chlorophyll concentration to carotenoid concentration (Fig. 2E). Then we measured the ion leakage of the root system and found that there was no difference in ion leakage rate under different treatment conditions (Fig. 2G).

The damage of organelles was observed by transmission electron microscope, and the intracellular ultrastructure of the morphological difference between the control and the Fe-deficiency stress was studied (Fig. 3A). Our results showed that under normal Fe supply condition, chloroplast was well developed, complete in structure, convex lens-shaped, rich in content, many thylakoids, clear and complete basal lamella, and chloroplast membrane, and multiple basal lamellas, neatly arranged. However, under the low Fe condition, the chloroplast deformed and became smaller and the chloroplast membrane was blurred; the lamella system could not be stacked to form basal particles, and most of the basal lamellas disappeared and became blurred.

Root growth is determined by the balance between cell division and cell elongation [15]. To study the role of Fe in the growth and development of the root system, we observed the cell size of the meristematic zone and elongation zone in the root system grown under normal or low Fe stress (Fig. 3B). The meristem size was calculated as the area where isodiametric cells extend from the quiescent center (QC) to the cell, which was twice the length of the immediately preceding cell [16]. The boundary of the transition zone is different in each cell type, so in all the analyses performed here, the cortical cell file was used to define the boundary [15]. Our results showed that compared with the control, the number of cells in the meristematic zone of the root system under Fe-deficient condition was significantly reduced, and the length of the cells in the elongation zone was significantly shorter (Fig. 3B-D).

Under the low Fe treatment, the root-to-shoot ratio was reduced but the difference was not significant. It could be inferred that the low Fe stress inhibited the growth of the shoots higher than that of the roots (Fig. 4A). Subsequently, the specific effects of low Fe stress on the growth of the wheat root system were analyzed. Low Fe stress was significantly reduced root-related parameters, including total root length, main root length, root surface area, root volume, number of root tips, average root diameter, and lateral root length (Fig. 4B-H). Then we measured the root activity of wheat and found that the root activity of wheat seedlings decreased significantly under low Fe stress, indicating that low Fe stress reduced the ability of roots to absorb Fe (Fig. 4I). These results indicated that the growth and development of wheat seedlings under low Fe stress were severely inhibited.

Physiological response of wheat to low Fe stress

To further understand the physiological response of wheat to low Fe stress, the contents of some osmotic adjustment substances that may participate in the regulation of low Fe stress resistance were tested. Under low Fe stress, the Pro concentration in the shoots and roots was increased significantly, but the concentration of MAD was only increased significantly in the shoots (Fig. 5A, D). The content of MDA is a sign of plasma membrane peroxidation. Our results showed that compared with the control, there was no significant difference in the concentration of MAD in roots under the low Fe stress, which indicated that the low Fe stress did not cause peroxidative damage to the plasma membrane (Fig. 5D). Under the stress of low Fe, the concentrations of OFR and H₂O₂ in the shoots and roots were significantly higher than those of the control condition (Fig. 5B, C). The results showed that low Fe stress accumulated more ROS.

Fe deprivation caused ROS accumulation

When plants are stressed by Fe deficiency, they will produce a large amount of ROS, which will cause oxidative damage. To analyze the effect of low Fe on the accumulation of ROS in wheat plants, DAB and NBT staining were performed on the shoots and roots of wheat after different treatments to observe the accumulation of H₂O₂ and O₂^¯ (Fig. 6A, B). The results showed that compared with normal Fe treatment, the DAB staining and NBT staining colors of wheat plant shoots and roots under low Fe stress were dark brown or dark blue, respectively, indicating that the shoots under low Fe stress, the concentrations of H₂O₂ and O₂^¯ in the roots and roots were significantly higher than those of the control (Fig. 6A, B). These research results showed that wheat plants produced more ROS after low Fe treatment. This result was consistent with the conclusion of Fig. 5B and C, that was, low Fe stress increased the accumulation of ROS.

Effect of low Fe stress on mineral nutrient elements in wheat

Inductively coupled plasma mass spectrometry (ICP-MS) is used to determine the ion expression profile of control and low Fe stress. ICP-MS quantitative analysis showed that the Fe⁺ concentration in the shoots and roots of wheat was significantly reduced under the low Fe stress treatment (Fig. 7B). Subsequently, we tested the concentrations of other metal cations, including Cu²⁺, Mn²⁺, Zn²⁺, Mg²⁺, Na⁺, Ca²⁺, Cd²⁺, and K⁺ (Fig. 7). In general, the response of low Fe stress to Cu²⁺, Mg²⁺, Na⁺, Cd^2+, and K⁺ was similar. Under low Fe stress, the concentrations of Cu²⁺, Mg²⁺, Na⁺, Cd²⁺, and K⁺ in the shoots did not change significantly; on the contrary, the concentrations of these cations increased significantly in the roots (Fig. 7A, E, F, H, I). Under low Fe stress, the Mn²⁺ concentration was increased significantly in the shoots and roots (Fig. 7C). Unlike the above cations, the concentrations of Zn²⁺ and Ca²⁺ did not change significantly in the shoots and roots (Fig. 7D, G).

Transcriptional response of wheat to low Fe stress

To understand the molecular basis of Fe tolerance in wheat seedlings under low Fe stress, an RNA-seq library was established using the shoots and roots of low Fe seedlings as materials. The correlation between gene expression levels in samples is an important indicator for testing the reliability of experiments. In this study, under the same treatment, the Pearson correlation coefficient between each pair of biological replications was mostly higher than 0.90, indicating that the similarity between samples in the transcriptome was very high (Figure S1). RNA-Seq generated more than 1,145,14 million Raw reads (Table 1). Among these reads, the GC content of libraries was about 54.00%. After quality control, 1,134,550,000 Clean reads were obtained, the average error rate of sequencing bases was less than 0.0249%, the Q20 value was more than 98.03%, and the Q30 value was more than 94.06%. Among them, 7,592-10.23 million Clean reads were located in the wheat genome (Table 1). Principal component analysis showed that there were significant differences in the expression patterns between different treatments and different wheat tissues (Fig. 8A). 2,861 genes were identified as differentially expressed (FDR <0.05, fold change ≥ 2). There were a total of 137 differentially expressed genes (DEGs) in the shoots and roots (Fig. 8B). The number of DEGs in the shoots was relatively small (378), under Fe-deficient stress, 313 genes were up-regulated and 65 genes were down-regulated (Fig. 8B, C). There were a large number of DEGs in the roots (2,619), under low Fe stress, there were 1,589 up-regulated genes and 1,030 down-regulated genes (Fig. 8B, D).

To obtain the functional information of DEGs, GO annotation analysis was performed on DEGs, and DEGs were classified into three types: biological process (BP), cell component (CC), and molecular function (MF). GO entries with p-value < 0.001 were considered to be significantly enriched. In the shoots and roots under the stress of low Fe, BP was mainly enriched in the metabolic process and cellular process (Fig. 8E); CC was mainly enriched in the cell part and membrane part (Fig. 8E); in the MF, binding and catalytic activity were the two most abundant GO terms (Fig. 8E). The KEGG annotation divided the pathways that DEGs participate in into five categories: Cellular Processes, Environmental Information Processing, Genetic Information Processing, Metabolism, and Organismal Systems. Most DEGs in Cellular Processes were annotated to Endocytosis pathway, most DEGs in Environmental Information Processing were annotated to MAPK signaling pathway–plant pathway, most DEGs in Genetic Information Processing were annotated to Protein processing in endoplasmic reticulum pathway, most of Metabolism DEGs were annotated to the Phenylpropanoid biosynthesis pathway, and most DEGs in Organic Systems were annotated to the Plant-pathogen interaction pathway (Fig. 8F).

To further deepen our understanding of the functions of these DEGs, we conducted GO enrichment analysis (Fig. 8G). Among them, GO terms with p-value < 0.001 were considered to be significantly enriched. We found that most DEGs were related to plant metabolism and biosynthesis processes, such astricarboxylic acid biosynthetic process, nicotianamine synthase activity, nicotianamine metabolic process, nicotianamine biosynthetic process, nicotianamine aminotransferase activity, and L−methionine biosynthetic process, etc. These results indicated that the synthesis of Fe chelator and Fe transport was regulated to maintain Fe homeostasis. Significantly enriched GO terms related to ion transport and ion homeostasis included iron ion transmembrane transport, intracellular sequestering of iron ion, sequestering of metal ion and sequestering of iron ion, etc. The above results indicated that the transporter played an important role during the low Fe period. In addition, the GO terms related to photosynthesis, such as photosynthesis, light harvesting in photosystem I, were also significantly enriched, indicating that low Fe affected the photosynthesis of wheat, which was closely related to the decrease of chlorophyll content in the shoots (Fig. 2D, E, F). These results emphasized the importance of membrane and/or membrane-localized metal ion transporters and regulatory and metabolic proteins in the stress of low Fe.

The KEGG database was used to further determine the pathways involved in the response of wheat to low Fe stress (Fig. 8H). Studies have shown that in plants exposed to abiotic stress, the accumulation of amino acids was believed to have a beneficial effect in the process of plant stress adaptation [17-19]. Most of the DEGs were involved in amino acid metabolism pathways, such as Alanine, aspartate and glutamate metabolism, Tyrosine metabolism, Cysteine, and methionine metabolism and Phenylalanine metabolism, and other KEGG pathways, which indicated that amino acid metabolism and protein synthesis played an important role in the response of wheat to Fendeficiency stress. DEGs were also involved in pathways related to photosynthesis, such as Photosynthesis − antenna proteins and Carotenoid biosynthesis, which indicated that low Fe had a significant impact on the photosynthesis of wheat. In addition, Phenylpropanoid biosynthesis, MAPK signaling pathway-plant, and other KEGG pathways were also significantly enriched (Fig. 8H).

Transcriptional response of Fe transport-related genes under low Fe stress

GO enrichment indicated that the synthesis of Fe chelator and Fe transport played an important role during low Fe stress. Among many DEGs, genes related to Fe ion homeostasis are the key genes for wheat tolerance to low Fe. Therefore, we analyzed the expression of genes related to Fe uptake and transport under Fe-deficiency stress. Figure S2 showed a molecular model of genes involved in Fe absorption and transport in plant roots (Figure S2A), chloroplasts, mitochondria, and vacuoles (Figure S2B). Transcription results analysis showed that most of the key genes involved in strategy I uptake of Fe, such as Fe³⁺ chelate reductase (FRO), iron-regulated transporter (IRT), natural resistance-associated macrophage protein (NRAPM), and Yellow stripe-like (YSL), were significantly up-regulated under Fe-deficiency stress. Among them, TaFRO2-2A (Figure S3A) and TaIRT1a-4A (Figure S3E), TaNRAMP2-4A (Figure S3B), and TaYSL15a-6D (Figure S3C) were up-regulated by 3 times, 6 times, 5 times, and 70 times, respectively, which might play a decisive role in the process of Fe uptake by plants.

In strategy II, the expression abundance of genes involved in plant Fe vector synthesis and Fe uptake and transport-related genes was significantly up-regulated. The process of synthesis and secretion of plant Fe siderophores is very complicated and requires the participation of 13 enzymes, but there are mainly four key enzymes: methionine synthetase (SAM), Nicotinamide synthase (NAS), Nicotinamide aminotransferase (NAAT), and Deoxyergate synthase (DMAS). Transcription analysis showed that these four key enzyme genes were significantly induced under low Fe stress. Among them, TaSAM3-6B (Figure S3D), TaNAS1c-6D (Figure S3I), TaNAAT1b-1B (Figure S3G), and TaDMAS1-4A (Figure S3F) had the most differential expression fold changes, which were up-regulated by 3 times, 3105 times, and 2718, and 24 times, respectively. The GO enrichment results also found that the GO terms related to nicotinamide (NA) synthesis were significantly enriched under low Fe stress (Fig. 8G). It might be related to the up-regulation of the expression abundance of the four key enzyme genes mentioned above, SAM, NAS, NAAT, and DMAS. And the synthesis of plant Fe carrier is related to the methionine salvage pathway because SAM is the key precursor for the synthesis of malic acid substances (MAs), which are produced by methionine. Transcriptome KEGG enrichment analysis also showed that methionine metabolism, L-methionine biosynthetic process, and L-methionine salvage from the S-adenosylmethionine pathway were significantly enriched (Fig. 8H). The above results indicated that amino acid metabolism played an important role in wheat response to low Fe. In strategy II, genes involved in Fe³⁺ transport such as transporter of mugnetic acid (TOM) (Figure S3J) and multidrug and toxin efflux family (MATE) (Figure S3K) were significantly up-regulated under low Fe stress. Among them, the differential expression of TaTOM-2B and TaMATE-4A changed the most, up-regulated by 2 times and 20 times, respectively. It was speculated that it might play a key role in the transport of iron from the xylem to neighboring cells. Oligopeptide transporter (OPT) is a polypeptide transporter located on the plasma membrane of plant epidermal cells. Transcriptome analysis showed that it was significantly induced by low Fe, and the differential expression of TaOPT3-5B was significantly up-regulated by low Fe by 44 times (Figure S3M), which might play a key role in the process of transmembrane transport of Fe³⁺-PS complexes near epidermal cells to plants.

The process of vacuolar iron transporter ginseng (VIT) for vacuolar membrane positioning and vacuole loading of Fe into seed embryos not only promotes the development of seeds but is also crucial for seed germination. NRAMP is located on the vacuole membrane, mediates the transport of Fe from the vacuole of the seed embryo to the cytoplasm, and participates in the supply of Fe during seed germination. Transcriptome results showed that the expression abundance of most VIT genes was significantly down-regulated under low Fe stress, especially TaVIT2-5B was down-regulated by 3 times (Figure S3R), which might play a key role in promoting seed development. NRAMP3/4 was significantly induced by low Fe stress, and TaNRAMP3-7D was up-regulated and expressed 3 times, which might play an important role in the transport of Fe from the vacuole of the seed embryo to the cytoplasm (Figure S3B). Mitochondrial m-type thioredoxin in chloroplast (ATM) is mainly involved in the process of exporting Fe-S to mitochondria, thereby maintaining the balance of Fe ion in mitochondria. Transcriptome analysis found that TraesCS5A02G110300 had the highest expression abundance in ATM, which might play an important role in maintaining Fe ion balance in mitochondria (Figure S3H). Permease in chioroplasts (PIC) is the first protein identified to participate in chloroplast Fe transport. Its expression is not regulated by Fe and is a constitutively expressed protein. Transcriptome analysis showed that PIC was significantly induced by low Fe in the shoots, and TraesCS1A02G165700 had the highest expression abundance (Figure S3L). Iron efflux transporter ferroportin (FPN) located in the chloroplast can regulate the intracellular Fe content by participating in the transport of NA or Fe-NA complex into the chloroplast. Transcriptome results showed that FPN was significantly induced by low Fe in roots, and TaFPN1-7A was significantly up-regulated by 60 times (Figure S3P), which might play a key role in maintaining Fe homeostasis in chloroplasts under low Fe stress. Nonintrinsic ABC protein (NAP) is located on the inner chloroplast membrane. It encodes a non-plasma membrane-localized nucleotide-binding domain subunit of the ABC transporter. Transcriptome analysis showed that NAP was significantly induced by low Fe in the shoots, and TraesCS3A02G047200 was significantly up-regulated. This result indicated that it might be part of the chloroplast Fe-ABC transport complex and played an important role in regulating the dynamic balance of Fe in the chloroplast (Figure S3Q). Another contributor to promoting the transfer of Fe ion into the chloroplast is mitoferrin-like, which contains chloroplast transit peptides and is mainly located in rosette leaves. Transcriptome results showed that the expression abundance of most mitoferrin or mitoferrin-like under low Fe stress decreased significantly, and the expression levels of TraesCS6D02G228200 and TraesCSU02G183100 were significantly down-regulated. It was speculated that they might play an important role in maintaining Fe content in chloroplasts (Figure S3N, O). The above results emphasized the importance of Fe ion uptake and transport proteins and regulatory and metabolic proteins in the stress of low Fe.

Involved in the transcriptional response of photosynthesis-related genes under low Fe stress

The GO enrichment analysis showed that GO terms related to photosynthesis, such as photosynthesis, light harvesting in photosystem I, were significantly enriched (Fig. 8G). And the KEGG enrichment analysis also showed that Photosynthesis − antenna proteins and carotenoid biosynthesis KEGG pathways were significantly enriched under low Fe stress (Fig. 8H), which indicated that low Fe significantly affected the photosynthesis of wheat plants. Many genes related to photosynthesis, such as photosynthetic antenna protein and key genes in the biosynthesis of carotenoids were identified (Figure S4). Through KEGG pathway analysis, a total of 33 genes encoding photosynthesis-antenna proteins (Figure S4A) and 39 genes encoding carotenoid biosynthesis (Figure S4B) were screened, and their expression levels were shown in Figure S4. The results showed that the expression levels of most genes involved in the photosynthetic antenna protein pathway and carotenoid biosynthetic pathway were significantly induced by low Fe (Figure S4A, B). The ratio of the concentration of chlorophyll a to chlorophyll b under low Fe stress was significantly higher than that under the control condition (Fig. 2E), which was a characteristic of plants adapting to adversity stress. Among them, chlorophyll an oxygenase (CAO), chlorophyll (ide) b reductase (CBR), chlorophyll synthase (CS), and chlorophyllase play an important role in the conversion between chlorophyll a and chlorophyll b (Figure S4C). Among the DEGs involved in photosynthesis, we found that the genes encoding these enzymes showed an increased expression pattern under low Fe stress (Figure S4C). There are two synthetic pathways for plant carotenoid substances: MEP pathway and MVA pathway. These two pathways generate geranylgeranyl diphosphate (GGPP), the direct precursor substance for carotenoid synthesis. GGPP generates the first carotenoid substance phytoene. After dehydrogenation, cyclization, hydroxylation, epoxidation, etc., it is transformed into other carotenoids. Transcriptome analysis showed that the expression of key enzyme genes in the process of carotenoid syntheses such as phytoene synthase (PSY), β-carotene hydroxylase (β-OHase), zeaxanthin epoxidase (ZEP), 9-cis-epoxycarotenoid dioxygenase (NCED), LUT5, and violaxanthin de-epoxidase (VDE) was significantly inhibited under low Fe stress (Figure S4D). And we found that Chlorophyll a/b binding proteins (Figure S4E), Photosynthesis II reaction proteins (Figure S4F), RuBisCo subunit binding proteins (Figure S2G), Mg²⁺ chelatase (Figure S4H), and Ribulose bisphosphate carboxylase/oxygenase activase (Figure S4I) were significantly down-regulated under low Fe stress. These results all indicated that the severe degradation of wheat chlorophyll and photosynthesis was inhibited under low Fe stress.

Transcriptional response of other ion transporters to low Fe stress

Among the numerous DEGs, genes related to ion homeostasis are the key genes for wheat resistance to Fe-deficiency stress. Fig. 9A showed a molecular model of the key genes responsible for regulating K⁺, Na⁺/Cl^-, and Ca²⁺ transport. Transcription analysis showed that most of the K⁺ transporter genes, including the K⁺ efflux transporter gene KEA located on the chloroplast, the vacuolar K⁺ inflow transporter gene KCO, the plasma membrane-located K⁺ inflow transporter gene AKT/KAT and HKT, and the K⁺ efflux gene SKOR was up-regulated under low Fe stress (Fig. 9B). The expression of Na⁺/H⁺ antiporter gene NHX, especially NHX2 involved in the Na⁺ compartment in the vacuole, and SOS1/NHX7 regulated the outflow of Na⁺ in the cell, were significantly up-regulated. To avoided cell damage due to excessive accumulation of Na⁺ in the cell (Fig. 9C). NHD could regulate the efflux of Na⁺, and its expression level was significantly up-regulated under low Fe stress, especially the expression of TaNHD1-5B was significantly up-regulated by 2 times, which might help reduce chloroplast damage caused by excessive Na⁺ (Fig. 9C). The chloroplast-localized bile acid: Na⁺ cotransporter (BASS) could regulate Na⁺ influx, and its expression level was significantly inhibited by low Fe stress, especially the expression level of TaBASS1D-1D was significantly down-regulated by 3 times (Fig. 9C). In addition, low Fe stress also induced the expression of most ALMT genes, which were related to the transport of Cl^- in vacuoles, especially the expression of TaALMT12-1B was significantly up-regulated by a factor of 2 under low Fe stress (Fig. 9C). Subsequently, the gene expression related to Na⁺/Ca²⁺ transport was studied, including CCX, CAX, ANXD, GLR, and CNGC. The expression levels of most genes changed significantly under low Fe stress. This result indicated that these genes might be involved in ion homeostasis under low Fe stress (Fig. 9D). For example, under low Fe stress, TaCCX expression was increased, TaCAX expression was decreased, and TaANXD expression was increased. This might be an important reason for reducing cytoplasmic Na⁺ and increasing cytoplasmic Ca²⁺ concentration (Fig. 9D). Most of the plasma membrane-localized NSCC, including GLR and CNGC, were up-regulated under low Fe stress (Fig. 9D). This might be an important way for Ca²⁺ and Na⁺ to enter the cytoplasm under low Fe stress (Fig. 9D). In addition, under low Fe stress, the expression levels of genes PHT2;1 (Fig. 9E) and COPT (Fig. 9F) involved in the uptake and transport of Pi and Cu²⁺ by roots were significantly down-regulated, while genes involved in Mg²⁺ uptake and transport were significantly induced by low Fe tress (Fig. 9G). The above results all indicated that the maintenance of ion homeostasis in cells was of great significance for plants to respond to low Fe tress.

Expression profile of cell cycle-related genes and ROS metabolism-related genes under low Fe stress

To study whether the inhibitory effect of Fe deficiency on root growth is related to the degree of cell division and differentiation, we combined the transcriptome to analyze the rich expression of cell cycle-related genes (Figure S5A). Cyclin and cyclin-dependent kinase belong to cell cycle control proteins, which play a key role in the process of mitosis [20].LRP (lateral root primordia) is specifically expressed in adventitious roots and lateral root primordia. It is a transcription factor that regulates root elongation and has strong tissue specificity [20]. Transcriptome analysis showed that under low Fe tress, the expression abundance of these genes was significantly down-regulated (Figure S5A). It indicated that the shortening of roots under low Fe stress might be due to the decrease in the division frequency of cells in the meristematic zone and the shortened cell elongation in the elongation zone.

According to the results in Fig. 5 and Fig. 6, we combined the transcriptome to analyze the expression of key enzyme genes in the ROS metabolic pathway (Figure S5B, C). The respiratory burst oxidase homologous gene (RBOH) encodes NADPH oxidase and plays a key role in the synthesis of ROS [21]. Under low Fe stress, 11 RBOH differentially expressed genes were significantly up-regulated in shoots or roots (Figure S5B). The expression of SOD, CAT, and PDH genes in most shoots and roots increased under low Fe stress, which might be necessary for ROS clearance under low Fe stress (Figure S5C).

Gene co-expression network analysis

In allohexaploid wheat, multiple copies of genes are ubiquitous. Therefore, the identification of core genes is an important prerequisite for understanding the molecular mechanisms of important agronomic traits. The systematic analysis of the transcriptional response of genes related to Fe uptake and transport under low Fe tress will help us to fully understand the adaptation mechanism of plants to Fe-deficiency stress. To identify the core members of genes related to Fe uptake and transport, we constructed a gene co-expression network (Figure S6A). The results showed that some Fe²⁺ absorption and transport related genes TaIRT1b-4A and plant siderophores related genes such as TaNAS2-6D, TaNAS1a-6A, TaNAS1-6B and TaNAAT1b-1D were identified as core target genes, which might play a key role in plant response to low Fe stress (Figure S6A).

Worldwide, low Fe tress severely inhibits plant growth, yield, and crop quality. Allohexaploid wheat is a major food crop and is highly sensitive to the stress of low Fe. Analyzing the physiological and molecular mechanisms of wheat tolerance to low Fe tress is an effective way to cultivate low Fe wheat varieties.

Transcriptomics-assisted analysis of morphophysiological responses of wheat to low Fe tress

In higher plants, low Fe is characterized by chlorosis between the veins of young leaves, while mature leaves remain green due to Fe immobility [22]. Our research results showed that under low Fe stress, wheat seedlings had obvious plant damage, including yellowing of leaves (Fig. 1) and inhibition of root growth (Fig. 1, 3). Although Fe is not a component of chlorophyll, it is indispensable in the synthesis of chlorophyll precursors. It is also an important indicator for studying wheat growth characteristics, physiological changes, and Fe nutritional status [23]. The degree of leaves xanthosis caused by low Fe in crops is mainly expressed by chlorophyll content because chlorophyll content is an effective indicator for judging whether crops are lacking Fe [24, 25]. As a functional indicator of photosynthetic pigment and light response [26], the concentration ratio of chlorophyll a and chlorophyll b under low Fe tress was significantly lower than the control condition (Fig. 2D), indicating that low Fe tress had a more obvious inhibitory effect on chlorophyll a. In addition, the concentration of total chlorophyll and carotenoids were decreased under low Fe tress (Fig. 2D, F), which also indicated that chlorophyll was degraded or damaged photosynthesis under low Fe stress [27]. Whether the stress of plant nutrient deficiency or physiological nutrient deficiency can cause the metabolic disorder of its leaves, severe nutrient deficiency will affect the organ structure of the plant, resulting in changes in the cell and subcellular morphology of plant leaves. Previous studies in rice [28] and Arabidopsis thaliana [29] found that the reduction or number of chloroplasts would inhibit the photosynthetic rate, and the destruction of the chloroplast structure would limit the photosynthetic rate. In this study, low Fe tress caused chloroplast separation and chlorophyll degradation, which might further lead to photosynthesis inhibition (Fig. 2D-F, Fig. 3A). The analysis of DEGs showed that the photosynthesis-related KEGG pathway was highly enriched (Fig. 8H), and the genes related to the chlorophyll biosynthesis pathway were significantly down-regulated under low Fe tress (Figure S4). The root system is an organ for wheat to absorb essential mineral elements for growth and development. Because root growth is related to cell division and elongation rate, our analysis of the root system of wheat showed that under low Fe stress, cell division was reduced, resulting in a reduction in the elongation of meristem cells, which was related to the observed reduction in growth (Fig. 2B, C).

Some research results show that under the stress of low Fe, the content of ROS and osmotic adjustment substances change significantly [30]. When adversity stress forces the production of a large amount of ROS in plant cells, it will cause damage to macromolecular substances and other components in the cell, hindering the normal metabolism and growth of plants, and even death [31]. Our results showed that under the condition of low Fe, the concentration of OFR and H₂O₂ were increased significantly (Fig. 5B, C), which initiated membrane lipid peroxidation and damages the structure of the cell membrane system. However, the ion leakage rate (Fig. 2G) and MAD content (Fig. 5D) in the root system were not different from those of the control, indicating that the accumulation of ROS did not cause great damage to the cell membrane system in the wheat root system. The content of osmotic regulators such as Pro and MAD accumulated in the cells in a large amount, thereby regulating the osmotic potential of the cells and further maintaining the normal metabolism of the cells (Fig. 5A, D).

Transcriptomics-assisted analysis of ion response of wheat to low Fe tress

Fe homeostasis in plants is jointly participated by a variety of Fe ligands, transporters, and regulatory factors. Usually, Fe is combined with ligand for transportation or storage to reduce the risk of free Fe producing toxic ROS [32]. NA and Deoxysarcosine (DMA) are Fe chelating agents and are responsible for the acquisition and transport of Fe. NA has long-distance and cellular Fe transport functions in higher plants [33–35]. DMA is mainly used for the long-distance transportation of Fe in the xylem and phloem [36, 37]. We also found that GO terms related to NA and DAM synthesis were significantly enriched in our transcriptome (Fig. 8G). In addition, expressions of key enzymes in NA and DAM synthesis pathways, such as SAM, NAS, NAAT, and DMAS, were significantly up-regulated (Figure S3D, I, G, F), which was consistent with previous results in wheat [38–40] and rice [33–35]. Previous studies have shown that the expression abundance of genes related to Fe uptake and transport in mitochondria [41, 42], chloroplasts [43, 44], and vacuoles [45] are significantly up-regulated under Fe-deficiency stress. This conclusion was consistent with our transcriptome results (10H-R). The increased expression of these genes in young leaves under low Fe tress s might be to unload more Fe into mesophyll cells to meet the normal absorption and utilization of Fe by young leaves. But when the Fe content in the cell exceeds the limit, the plant can transport Fe²⁺-chelates into the organelles to store Fe such as vacuoles, avoiding cell toxicity. Studies have shown that when AtVTL1, AtVTL2, or AtVTL5 is expressed under the 35S promoter in the background of nramp3/nramp4 or vit1-1, in the two mutants [46], AtVTL1, AtVTL2 or AtVTL5 can restore root growth. In yeast, the VIT homologous Ca²⁺-Sensitive-Cross-Complementer (CCC1) protein is a vacuolar Fe transporter that transports Fe to the vacuole under condition of excess Fe [45]. The pastatin promoter is used on cassava for overexpansion, and the expression of atvit1 gene is positively correlated with the increase of Fe concentration in storage roots. This result shows that the VIT gene is involved in mediating the vacuolar retention of Fe [47]. The above results all indicate that the dynamic balance of Fe in cells is essential for the normal development of plants.

Fe not only plays an important role in plant physiological functions but also plays an important role in maintaining various cell processes. In the past few decades, more and more progress has been made in understanding how plants adapt to low Fe in the soil. However, how to observe the Fe dynamics in plants is still a significant challenge. Although there is a lot of information indicating that many genes are responsible for the absorption of Fe from the soil, the transportation from roots to shoots, the storage in cells, and even the regulation of transcription and post-transcriptional levels. However, further research is still needed to reveal the further connection and integration between Fe-deficient signaling pathways and developmental and physiological networks. This study analyzed the morphological and physiological changes of the shoots and roots under low Fe tress and combined the transcriptome to analyze the expression profiles of genes, which related to Fe uptake and transport, photosynthesis, and other ion transporter-related genes. Finally, key genes affecting plant response to low Fe stress were identified. In summary, these results might provide excellent genetic resources for the apparent regulatory mechanism of wheat resistance to low Fe tress and crop genetic improvement.

ATM: Mitochondrial m-type thioredoxin in chloroplast; BASS: chloroplast-localized bile acid; BP: biological process; CAO: chlorophyll an oxygenase; CBR: chlorophyll (ide) b reductase; CC: cell component; CCC1: Ca²⁺-Sensitive cross-Cross-Complementer; CS: chlorophyll synthase; CV: Coefficient of Variation; DAB: 3’3’-Diaminobenzidine; DEGs: differentially expressed genes; DMA: Deoxysarcosine; DMAS: Deoxyergate synthase; Fe: iron; FIT: fer-like iron deficiency-induced transcription factor; FPKM: fragments per kilobase of exon model per million mapped reads; FPN: Iron efflux transporter ferroportin; FRO: Fe³⁺ chelate reductase; GGPP: geranylgeranyl diphosphate; GO: Gene Ontology; H⁺-ATPase: proton ATPase; H₂O₂: hydrogen peroxide; ICP-MS: Inductively coupled plasma mass spectrometry; IRT: Iron-regulated transporter; KEGG: Kyoto Encyclopedia of Genes and Genomes; LRP: lateral root primordia; MATE: Multidrug and toxin efflux family; MDA: malondialdehyde; MF: molecular function; NAAT: Nicotinamide aminotransferase; NAP: Nonintrinsic ABC protein; NAS: Nicotinamide synthase; NBT: Nitro blue tetrazole; NCED: 9-cis-epoxycarotenoid dioxygenase; NRAMP: Natural resistance-associated macrophage protein; OFR: superoxide anion; OPT: Oligopeptide transporter; PIC: Permease in chioroplasts; Pro: proline; PS: phytosiderophore; PSY: phytoene synthase; QC: quiescent center; RBOH: respiratory burst oxidase homologous gene; RIN: RNA integrity number; RT-qPCR: Reverse transcription-quantitative polymerase chain reaction; SAM: Methionine synthetase; TOM: Transporter of mugnetic acid; TTC: tetrazolium; VDE: violaxanthin de-epoxidase; VIT: Vacuolar iron transporter ginseng; YSL: Yellow stripe-like; ZEP: zeaxanthin epoxidase; β-OHase: β-carotene hydroxylase.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Availability of data and materials

All data generated or analysed during this study are included in this published article [and its supplementary information files]. All the data and materials that are required to reproduce these findings can be shared by contacting Dr. Ying-peng Hua ([email protected]).

Plant Ethics

We declare that the collecting of plant material is allowed by the Convention on the Trade in Endangered Species of Wild Fauna and Flora: https://www.cites.org/, and plant material was collected with permission in accordance with national and international guidelines.

In this study, all the bread wheat (Triticum aestivum L.) cultivar zhengmai 1860 were obtained from the Henan Academy of Agricultural Sciences (Zhengzhou, Henan Province, China).

Competing interests

The authors declare that they have no competing interests.

Funding

This study was financially supported by the National Natural Science Foundation of China (31801923), Major Collaborative Innovation Project of Zhengzhou City (Key Discipline Construction Project of Zhengzhou University) (NO. xkzdjc201905), and Youth Innovation Project of Key discipline of Zhengzhou University (NO. XKZDQN202002).

Authors' contributions

Ying-peng Hua, Jin-yong Huang, Ting Zhou, and Jin-yong Huang planned and designed the research and analyzed data. Ying-peng Hua and Yue Wang performed the experiments. Yue Wang wrote the original manuscript. Ying-peng Hua reviewed and edited the manuscript. All authors read and approved the final manuscript.

Acknowledgments

We thanked all the colleagues in our laboratory for providing useful discussions and technical assistance. We are very grateful to the editor and reviewers for critically evaluating the manuscript and providing constructive comments for its improvement.

Author information

Affiliations

School of Agricultural Sciences, Zhengzhou University, Zhengzhou, 450001, China

Email address for each author:

Yue Wang: [email protected]

Ying-peng Hua: [email protected]

Ting Zhou: [email protected]

Jin-yong Huang: [email protected]n

Cai-peng Yue: [email protected]

Pii Y, Penn A, Terzano R, Crecchio C, Mimmo T, Cesco S. Plant-microorganism-soil interactions influence the Fe availability in the rhizosphere of cucumber plants. Plant Physiol Biochem. 2015;87:45–52.
Guerinot ML, Yi Y. Iron: utritious, noxious, and not readily available. Plant Physiol. 1994;104:815–820.
Kim SA, Guerinot ML. Mining iron: iron uptake and transport in plants. FEBS Lett. 2007;581:2273–80.
Briat JF, Curie C, Gaymard F. Iron utilization and metabolism in plants. Curr Opin Plant Biol. 2007;10:276–82.
Briat JF, Curie C, Gaymard F. Iron utilization and metabolism in plants. Curr Opin Plant Biol. 2007;10:276–82.
Römheld V, Marschner H. Evidence for a specific uptake system for iron phytosiderophores in roots of grasses. Plant Physiol. 1986;80:175–80.
Curie C, Panaviene Z, Loulergue C, Dellaporta SL, Briat JF, Walker EL. Maize yellow stripe1 encodes a membrane protein directly involved in Fe(III) uptake. Nature. 2001;409:346–9.
Dong K, Zhen S, Cheng Z, Cao H, Ge P, Yan Y. Proteomic analysis reveals key proteins and phosphoproteins upon seed germination of wheat (Triticum aestivum L.). Front Plant Sci. 2015;6:1017.
Spiller S, Terry N. Limiting factors in photosynthesis: II. Iron stress reduces photochemical capacity by reducing the number of photosynthetic units. Plant Physiol. 1980;65:121–5.
Garnica M, Bacaicoa E, Mora V, San Francisco S, Baigorri R, Zamarreño AM, Garcia-Mina JM. Shoot iron status and auxin are involved in iron deficiency-induced phytosiderophores release in wheat. BMC Plant Biol. 2018;18:105.
Khan MA, Castro-Guerrero NA, McInturf SA, Nguyen NT, Dame AN, Wang J, Bindbeutel RK, Joshi T, Jurisson SS, Nusinow DA, Mendoza-Cozatl DG. Changes in iron availability in Arabidopsis are rapidly sensed in the leaf vasculature and impaired sensing leads to opposite transcriptional programs in leaves and roots. Plant Cell Environ. 2018;41:2263–2276.
Salomé PA, Oliva M, Weigel D, Krämer U. Circadian clock adjustment to plant iron status depends on chloroplast and phytochrome function. EMBO J. 2013;32:511–23.
Oburger E, Gruber B, Schindlegger Y, Schenkeveld WD, Hann S, Kraemer SM, Wenzel WW, Puschenreiter M. Root exudation of phytosiderophores from soil-grown wheat. New Phytol. 2014;203:1161–74.
Connorton JM, Jones ER, Rodríguez-Ramiro I, Fairweather-Tait S, Uauy C, Balk J. Wheat vacuolar iron transporter TaVIT2 transports Fe and Mn and is effective for biofortification. Plant Physiol. 2017;174:2434–2444.
Hetherington FM, Kakkar M, Topping JF, Lindsey K. Gibberellin signaling mediates lateral root inhibition in response to K⁺-deprivation. Plant Physiol. 2021;185:1198–1215.
González-García MP, Vilarrasa-Blasi J, Zhiponova M, Divol F, Mora-García S, Russinova E, Caño-Delgado AI. Brassinosteroids control meristem size by promoting cell cycle progression in Arabidopsis roots. Development. 2011;138:849–59.
Brosché M, Vinocur B, Alatalo ER, Lamminmäki A, Teichmann T, Ottow EA, Djilianov D, Afif D, Bogeat-Triboulot MB, Altman A, Polle A, Dreyer E, Rudd S, Paulin L, Auvinen P, Kangasjärvi J. Gene expression and metabolite profiling of populus euphratica growing in the negev desert. Genome Biol. 2005;6:R101.
Usadel B, Bläsing OE, Gibon Y, Poree F, Höhne M, Günter M, Trethewey R, Kamlage B, Poorter H, Stitt M. Multilevel genomic analysis of the response of transcripts, enzyme activities and metabolites in Arabidopsis rosettes to a progressive decrease of temperature in the non-freezing range. Plant Cell Environ. 2008;31:518–47.
Lugan R, Niogret MF, Leport L, Guégan JP, Larher FR, Savouré A, Kopka J, Bouchereau A. Metabolome and water homeostasis analysis of Thellungiella salsuginea suggests that dehydration tolerance is a key response to osmotic stress in this halophyte. Plant J. 2010;64:215–29.
Glab N, Labidi B, Qin LX, Trehin C, Bergounioux C, Meijer L. Olomoucine, an inhibitor of the cdc2/cdk2 kinases activity, blocks plant cells at the G1 to S and G2 to M cell cycle transitions. FEBS Lett. 1994;353:207–11.
Kaya H, Takeda S, Kobayashi MJ, Kimura S, Iizuka A, Imai A, Hishinuma H, Kawarazaki T, Mori K, Yamamoto Y, Murakami Y, Nakauchi A, Abe M, Kuchitsu K. Comparative analysis of the reactive oxygen species-producing enzymatic activity of Arabidopsis NADPH oxidases. Plant J. 2019;98:291–300.
Li Y, Wang N, Zhao F, Song X, Yin Z, Huang R, Zhang C. Changes in the transcriptomic profiles of maize roots in response to iron-deficiency stress. Plant Mol Biol. 2014;85:349–63.
Hu X, Page MT, Sumida A, Tanaka A, Terry MJ, Tanaka R. The iron-sulfur cluster biosynthesis protein SUFB is required for chlorophyll synthesis, but not phytochrome signaling. Plant J. 2017;89:1184–1194.
Kalaji HM, Bąba W, Gediga K, Goltsev V, Samborska IA, Cetner MD, Dimitrova S, Piszcz U, Bielecki K, Karmowska K, Dankov K, Kompała-Bąba A. Chlorophyll fluorescence as a tool for nutrient status identification in rapeseed plants. Photosynth Res. 2018;136:329–343.
Gupta OP, Pandey V, Saini R, Narwal S, Malik VK, Khandale T, Ram S, Singh GP. Identifying transcripts associated with efficient transport and accumulation of Fe and Zn in hexaploid wheat (T. aestivum L.). J Biotechnol. 2020;316:46–55.
Aman R, Carle R, Conrad J, Beifuss U, Schieber A. Isolation of carotenoids from plant materials and dietary supplements by high-speed counter-current chromatography. J Chromatogr A. 2005;1074:99–105.
Guo A, Hu Y, Shi M, Wang H, Wu Y, Wang Y. Effects of iron deficiency and exogenous sucrose on the intermediates of chlorophyll biosynthesis in Malus halliana. PLoS One. 2020;15:e0232694.
Li Y, Ren B, Yang X, Xu G, Shen Q, Guo S. Chloroplast downsizing under nitrate nutrition restrained mesophyll conductance and photosynthesis in rice (Oryza sativa L.) under drought conditions. Plant Cell Physiol. 2012;53:892–900.
Weise SE, Carr DJ, Bourke AM, Hanson DT, Swarthout D, Sharkey TD. The arc mutants of Arabidopsis with fewer large chloroplasts have a lower mesophyll conductance. Photosynth Res. 2015;124:117–26.
Chu Q, Sha Z, Maruyama H, Yang L, Pan G, Xue L, Watanabe T. Metabolic reprogramming in nodules, roots, and leaves of symbiotic soybean in response to iron deficiency. Plant Cell Environ. 2019;42:3027–3043.
Nadarajah KK. ROS homeostasis in abiotic stress tolerance in plants. Int J Mol Sci. 2020;21:5208.
Alvarez-Fernández A, Díaz-Benito P, Abadía A, López-Millán AF, Abadía J. Metal species involved in long distance metal transport in plants. Front Plant Sci. 2014;5:105.
Bashir K, Inoue H, Nagasaka S, Takahashi M, Nakanishi H, Mori S, Nishizawa NK. Cloning and characterization of deoxymugineic acid synthase genes from graminaceous plants. J Biol Chem. 2006;281:32395–402.
Inoue H, Higuchi K, Takahashi M, Nakanishi H, Mori S, Nishizawa NK. Three rice nicotianamine synthase genes, OsNAS1, OsNAS2, and OsNAS3 are expressed in cells involved in long-distance transport of iron and differentially regulated by iron. Plant J. 2003;36:366–81.
Inoue H, Takahashi M, Kobayashi T, Suzuki M, Nakanishi H, Mori S, Nishizawa NK. Identification and localisation of the rice nicotianamine aminotransferase gene OsNAAT1 expression suggests the site of phytosiderophore synthesis in rice. Plant Mol Biol. 2008;66:193–203.
Kakei Y, Yamaguchi I, Kobayashi T, Takahashi M, Nakanishi H, Yamakawa T, Nishizawa NK. A highly sensitive, quick and simple quantification method for nicotianamine and 2'-deoxymugineic acid from minimum samples using LC/ESI-TOF-MS achieves functional analysis of these components in plants. Plant Cell Physiol. 2009;50:1988–93.
Nishiyama R, Kato M, Nagata S, Yanagisawa S, Yoneyama T. Identification of Zn-nicotianamine and Fe-2'-Deoxymugineic acid in the phloem sap from rice plants (Oryza sativa L.). Plant Cell Physiol. 2012;53:381–90.
Kaur G, Shukla V, Kumar A, Kaur M, Goel P, Singh P, Shukla A, Meena V, Kaur J, Singh J, Mantri S, Rouached H, Pandey AK. Integrative analysis of hexaploid wheat roots identifies signature components during iron starvation. J Exp Bot. 2019;70:6141–6161.
Bonneau J, Baumann U, Beasley J, Li Y, Johnson AA. Identification and molecular characterization of the nicotianamine synthase gene family in bread wheat. Plant Biotechnol J. 2016;14:2228–2239.
Beasley JT, Bonneau JP, Johnson AAT. Characterisation of the nicotianamine aminotransferase and deoxymugineic acid synthase genes essential to Strategy II iron uptake in bread wheat (Triticum aestivum L.). PLoS One. 2017;12:e0177061.
Jain A, Connolly EL. Mitochondrial iron transport and homeostasis in plants. Front Plant Sci. 2013;4:348.
Bashir K, Ishimaru Y, Nishizawa NK. Identification and characterization of the major mitochondrial Fe transporter in rice. Plant Signal Behav. 2011;6:1591–3.
Müller B, Kovács K, Pham HD, Kavak Y, Pechoušek J, Machala L, Zbořil R, Szenthe K, Abadía J, Fodor F, Klencsár Z, Solti Á. Chloroplasts preferentially take up ferric-citrate over iron-nicotianamine complexes in Brassica napus. Planta. 2019;249:751–763.
Han J, Song X, Li P, Yang H, Yin L. Maize ZmFDR3 localized in chloroplasts is involved in iron transport. Sci China C Life Sci. 2009;52:864–71.
Ram H, Sardar S, Gandass N. Vacuolar iron transporter (Like) proteins: regulators of cellular iron accumulation in plants. Physiol Plant. 2021;171:823–832.
Gollhofer J, Timofeev R, Lan P, Schmidt W, Buckhout TJ. Vacuolar-Iron-Transporter1-Like proteins mediate iron homeostasis in Arabidopsis. PLoS One. 2014;9:e110468.
Narayanan N, Beyene G, Chauhan RD, Gaitán-Solis E, Grusak MA, Taylor N, Anderson P. Overexpression of Arabidopsis VIT1 increases accumulation of iron in cassava roots and stems. Plant Sci. 2015;240:170–81.
Arnon DI. Copper enzymes in isolated chloroplasts. Polyphenol oxidase in beta vulgaris. Plant Physiol. 1949;24:1–15.
Sakuraba Y, Park SY, Kim YS, Wang SH, Yoo SC, Hörtensteiner S, Paek NC. Arabidopsis STAY-GREEN2 is a negative regulator of chlorophyll degradation during leaf senescence. Mol Plant. 2014;7:1288–1302.
Persson DP, Chen A, Aarts MG, Salt DE, Schjoerring JK, Husted S. Multi-element bioimaging of Arabidopsis thaliana roots. Plant Physiol. 2016;172:835–847.
Yuan J, Wang X, Zhou H, Li Y, Zhang J, Yu S, Wang M, Hao M, Zhao Q, Liu L, Li M, Li J. Comparison of Sample Preparation Techniques for Inspection of Leaf Epidermises Using Light Microscopy and Scanning Electronic Microscopy. Front Plant Sci. 2020;11:133.
Kong L, Wang F, Feng B, Li S, Si J, Zhang B. The structural and photosynthetic characteristics of the exposed peduncle of wheat (Triticum aestivum L.): an important photosynthate source for grain-filling. BMC Plant Biol. 2010;10:141.
Jambunathan N. Determination and detection of reactive oxygen species (ROS), lipid peroxidation, and electrolyte leakage in plants. Methods Mol Biol. 2010;639:292–8.
Shaw AK, Hossain Z. Impact of nano-CuO stress on rice (Oryza sativa L.) seedlings. Chemosphere. 2013;93:906–15.
Wan HP, Chen LL, Guo JB, Li Q, Wen J, Yi B, Ma CZ, Tu JX, Fu TD, Shen JX. Genome-wide association study reveals the genetic architecture underlying salt tolerance-related traits in rapeseed (Brassica napus L.). Front Plant Sci. 2017;8:593.
Mi H, Lazareva-Ulitsky B, Loo R, et al. The PANTHER database of protein families, subfamilies, functions and pathways. Nucleic Acids Res. 2005;33:D284-D288.
Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36:D480-D484.
Eisen MB, Spellman PT, Brown PO, Botstein D. Cluster analysis and display of genome-wide expression patterns. P Natl Acad Sci USA. 1998;95:14863–68.
Maillard A, Etienne P, Diquélou S, Trouverie J, Billard V, Yvin JC, Ourry A. Nutrient deficiencies modify the ionomic composition of plant tissues: a focus on cross-talk between molybdenum and other nutrients in Brassica napus. J Exp Bot. 2016;67:5631–5641.
Chang S, Chen W, Yang J. Another formula for calculating the gene change rate in real-time RT-PCR. Mol Biol Rep. 2009;36:2165–8.
Li M, Li D, Tang Y, Wu F, Wang J. CytoCluster: A Cytoscape Plugin for Cluster Analysis and Visualization of Biological Networks. Int J Mol Sci. 2017;18:1880.

Table 1 Statistics of wheat transcriptome sequencing data under low Fe stress
Sample	Raw reads	Clean reads	Mapped Reads	Error rate(%)	Q20(%)	Q30(%)	GC content(%)
HIR1	103,367,072	102,485,512	93,941,511 (91.66%)	0.0244	98.24	94.79	55.13
HIR2	92,962,154	91,950,192	82,671,080 (89.91%)	0.0246	98.16	94.60	54.20
HIR3	93,650,802	92,663,246	83,454,929 (90.06%)	0.0246	98.17	94.63	54.67
HIS1	81,831,040	81,022,372	75,929,551 (93.71%)	0.0249	98.03	94.26	55.55
HIS2	97,354,592	96,467,486	90,270,260 (93.58%)	0.0244	98.21	94.75	56.63
HIS3	109,751,026	108,750,122	102,330,323 (94.10%)	0.0244	98.24	94.80	55.66
LIR1	88,487,912	87,710,632	81,340,207 (92.74%)	0.0247	98.12	94.44	55.03
LIR2	84,443,470	83,702,758	76,906,608 (91.88%)	0.0248	98.08	94.38	54.53
LIR3	108,537,716	107,600,262	100,077,545 (93.01%)	0.0245	98.19	94.64	55.20
LIS1	85,080,798	84,380,564	79,656,284 (94.40%)	0.0241	98.33	95.06	57.10
LIS2	106,018,904	105,038,702	98,685,159 (93.95%)	0.0245	98.19	94.67	55.62
LIS3	93,658,786	92,779,306	87,098,185 (93.88%)	0.0247	98.10	94.43	55.92
Note: HIR (Hight Iron Root), HIS (Hight Iron Shoot), LIR (Low Iron Root), LIS (Low Iron Shoot). 1, 2, and 3 respectively represent different biological repetitions. Error rate (%) is the average error rate of sequencing bases corresponding to the quality control data. Q20% is the proportion of the nucleotide quality value larger than 20; Q30% is the proportion of the nucleotide quality value larger than 30; GC content (%) is proportion of guanidine and cytosine nucleotides among the total nucleotides.

No competing interests reported.

Download PDF

Editorial decision: Major revision
04 Mar, 2022
Reviews received at journal
19 Feb, 2022
Reviewers agreed at journal
09 Feb, 2022
Reviewers agreed at journal
02 Feb, 2022
Reviewers invited by journal
02 Feb, 2022
Editor assigned by journal
02 Feb, 2022
Editor invited by journal
01 Feb, 2022
Submission checks completed at journal
01 Feb, 2022
First submitted to journal
11 Jan, 2022

You are reading this latest preprint version

Changes in the Transcriptomic Profiles of Allohexaploid Wheat in Response to Low Fe Stress

Status:

Version 1

Abstract

Figures

Background

Materials And Methods

Results

Discussion

Conclusions

Abbreviations

Declarations

References

Tables

Additional Declarations

Supplementary Files

Status:

Version 1