- 1College of Life Sciences and Agronomy, Zhoukou Normal University, Zhoukou, China
- 2School of Mechanical and Electrical Engineering, Zhoukou Normal University, Zhoukou, China
Many soils face dual challenges of cadmium (Cd) contamination and salinization. However, the response of crops, especially wheat, to combined Cd and salinity stress is not understood. Here, wheat was grown in a hydroponic model for 14 days under single and combined Cd and NaCl stresses. Growth parameters, tissue Cd2+ and Na+ contents, and leaf chlorophyll (Chl), O2•−, and MDA levels were determined. Comparative transcriptomic and metabolomic analyses of the leaves were performed. The results showed that combined stress had a greater inhibitory effect on Chl contents and generated more O2•− and MDA, resulting in more severe wheat growth retardation than those under Cd or NaCl stress. Stress-induced decrease in Chl levels may be attributed to the inhibition of Chl biosynthesis, activation of Chl degradation, or a decline in glutamate content. Cd addition weakened the promotional effect of NaCl on SOS1 gene expression, thereby increasing the Na+ content. Contrastingly, NaCl supplementation downregulated the Nramp and ZIP gene expressions related to Cd uptake and transport, thereby impeding Cd2+ accumulation. All stresses enhanced tryptophan content via promoting tryptophan biosynthesis. Meanwhile, Cd and NaCl stresses activated phenylpropanoid biosynthesis and purine metabolism, respectively, thereby increasing the levels of caffeic acid, fumaric acid, and uric acid. Activating the TCA cycle was important in the wheat’s response to combined stress. Additionally, NaCl and combined stresses affected starch and sucrose metabolism, resulting in sucrose and trehalose accumulation. Our findings provide a comprehensive understanding of the response of wheat to the combined Cd and salinity stress.
1 Introduction
Cadmium (Cd) is a non-essential heavy metal in organisms. Owing to mining and smelting activities, electroplating, industrial solid waste discharge, sewage irrigation, and overuse of pesticides and fertilizers, large amounts of Cd are released into the environment, especially into the soil (Uchimiya et al., 2020). According to a national soil pollution survey report announced by China’s Ministry of Environmental Protection, soil Cd has the highest over-standard rate (7.0%) among heavy metals contamination. Cd concentrations in soils (including agricultural soil) of the area around a Pb-zinc minein Liuzhou, China and a lead-zinc smelter in Zhuzhou, China ranged from 6.56–17.93 and 0.11–23.3 mg·kg-1, respectively (Yu et al., 2021; Zhou et al., 2022). In addition to China, other countries worldwide such as Japan, Thailand, Korea, Pakistan, India, Vietnam, and Peru also face soil Cd contamination (Arao et al., 2010; Sriprachote et al., 2012; Ahn et al., 2017; Rehman et al., 2017; Kumar et al., 2019; Nguyen et al., 2020; Thomas et al., 2023). As a highly toxic and non-biodegradable heavy metal, Cd damages crop growth and easily accumulates in crops, thereby endangering human health through the food chain (Hu et al., 2022). Therefore, Cd contamination has attracted considerable attention worldwide.
In addition to Cd contamination, soil salinization is a non-negligible abiotic stress for crops. Saline soils are estimated to cover 833 million hectares worldwide (https://www.fao.org). Even worse, saline agricultural soils are increasing at a rate of 10% every year due to climate change, brackish water irrigation, and industrial pollution (Li et al., 2021). High salinity often results in crop growth inhibition and yield loss via ionic toxicity, osmotic stress, and oxidative stress, causing economic losses exceeding 27.3 billion dollars annually (Qadir et al., 2014; Yang and Guo, 2018; Zhu, 2016).
Increasing evidences revealed that many soils were suffering from the dual challenges of Cd contamination and salinity stress, particularly in arid and semi-arid regions and coastal areas (Cheng et al., 2019; Lutts and Lefèvre, 2015). For instance, high levels of Cd (4.17 mg·kg-1) and salinity (EC 4.67 dS·m-1) were detected simultaneously in the soil from Murcia and its surroundings in Spain (Acosta et al., 2011). Furthermore, high soil Cd (2.0 mg·kg-1) and salinity (EC 6.55 dS·m-1) were also observed around an abandoned factory in Northern China (Wu et al., 2022). Therefore, crops may be subjected to dual stress from Cd and salinity.
Wheat (Triticum aestivum L.) is the most cultivated edible cereal globally, providing 20% of the daily calories and protein requirements for humans (Tricker et al., 2018). In 2021, the total wheat cultivation area was 220.8 million hectares, producing 770.9 million tons (https://www.fao.org/faostat). Since wheat is a moderately salinity-tolerant crop with strong Cd enrichment ability (Li and Zhou, 2019; Quan et al., 2021), it is of great significance to study the response of wheat to combined Cd and salinity stress to ensure food security and cultivate Cd- and salinity-tolerant varieties. However, studies on Cd and salinity combination in wheat are limited, mainly focusing on antioxidant properties (Abbas et al., 2018; Mühling and Läuchli, 2003; Shafi et al., 2009). Little is known about the transcriptomic and metabolomic mechanisms underlying the wheat response to a combination of Cd and salinity stresses.
In this study, Zhoumai 36, a winter wheat cultivar planted in the Huang-Huai-Hai region, China, was exposed to single and combined stresses of Cd and salinity for 14 days using a hydroponic model. The growth parameters were recorded, and leaf chlorophyll (Chl), superoxide anion (O2•−), malondialdehyde (MDA), tissue Na+, and tissue Cd2+ contents were measured. Comparative transcriptomic and metabolomic analyses of the leaves were conducted to reveal the underlying response mechanisms of wheat to single and combined stresses. The findings of this study will improve our understanding of wheat response and adaptation to combined Cd and salinity stress.
2 Materials and methods
2.1 Plant culture and exposure experiment
Wheat cultivar Zhoumai 36 seeds were surface-sterilized and germinated, as described in our previous study (Yue et al., 2022). Ten well-germinated seeds were transferred to each glass tissue culture bottle containing 325 mL of 1/2 Hoagland solution. The bottles were randomly assigned to four groups, all in quintuplicates: (1) Control group: no addition of NaCl and Cd2+; (2) NaCl group: 150 mM NaCl; (3) Cd group: 100 μM Cd2+; and (4) Combined NaCl+Cd group: 150 mM NaCl and 100 μM Cd2+. To avoid osmotic shock, NaCl (2.925 g) is added into the nutrient solution in increments of 50 mM every 24 h to reach the final concentration of 150 mM. Cd2+ was added as CdCl2·2.5 H2O. The NaCl concentration was adopted from our previous study (Yue et al., 2022), and the Cd concentration was selected based on the Cd tolerance experiment of Zhoumai 36 (Supplementary Figure S1) and soil Cd content in China (Yu et al., 2021; Zhou et al., 2022). The hydroponics experiment was performed for 14 days under 25 ± 0.5°C, 60 ± 3% relative humidity, and 14/10 h of light/dark period. The nutrient solution was supplemented daily to compensate for water loss caused by evaporation and plant transpiration and to ensure plants received as consistent nutrition and water as possible throughout the experiment.
2.2 Sample collection and growth parameters measurement
After 14 days of treatment, the wheat seedlings were harvested and washed thrice with deionized water. The primary root length and shoot height were measured and recorded. Subsequently, the roots and leaves were separated. A portion of the roots and leaves were dried for 48 h at 75°C in an oven and used to determine the elemental content after weighing the dry weight. The remaining leaves were immediately frozen in liquid nitrogen and stored at -80°C until biochemical, transcriptomic and metabolomic analyses.
2.3 Measurement of Chl content
Chl was extracted as previously described (Hiscox and Israelstam, 1979). Detailed procedures can be found in the Supplementary Methods. The absorbance of the extract solutions was measured at 645 and 663 nm using a SpectraMax i3x microplate reader (Molecular Devices, San Jose, USA). Finally, the total Chl, Chl a, and Chl b contents were calculated according to established equations (Arnon, 1949).
2.4 Measurement of tissues Na+ and Cd2+ content
Oven-dried root and shoot samples were digested in glass digestion tubes containing 8 mL concentrated HNO3 and 2 mL H2O2. The digestion was performed in a SPH108 digestion instrument (ALVA Instrument Co., Ltd., Jinan, China) under the following conditions: 12 min at 130°C, and 50–60 min at 180°C. Subsequently, the digested samples were diluted to 45 mL with deionized water and filtered using a 0.45 μm filter membrane. The Na+ and Cd2+ contents were determined using an A3AFG flame atomic absorption spectrophotometer (Purkinje General Instrument Co., Ltd., Beijing, China).
2.5 Determination of O2•− and MDA contents
O2•− content was determined using a superoxide anion content assay kit (Boxbio, Beijing, China). MDA content was determined using an MDA colorimetric assay kit (Elabscience, Wuhan, China). Briefly, 0.1 g fresh leaf tissue was homogenized with 1 mL PBS buffer (0.01M, pH7.4) and centrifuged for 10,000× g at 4°C for 20 min. The resulting supernatant was collected, and the O2•− and MDA contents were measured according to the manufacturer’s instructions. The supernatant protein concentration was determined by using a BCA protein colorimetric assay kit (Elabscience).
2.6 Leaf transcriptome sequencing, bioinformatics analysis, and quantitative real-time PCR validation
Total RNA was extracted from the leaf samples (n = 3 per group), and RNA quality assessment was performed as previously described (Yue et al., 2022). High-quality total RNA (1 μg) was used to construct RNA-seq transcriptome libraries using Stranded mRNA Prep Ligation (Illumina, San Diego, USA). Detailed library construction procedures are provided in the Supplementary Methods. After quantification by Qubit 4.0 Fluorometer (Thermo Fisher, Waltham, USA), the paired-end libraries were sequenced in an Illumina Novaseq 6000 system by Majorbio Biotech Co., Ltd, Shanghai, China.
Raw reads were processed using fastp software to remove the adapter, low-quality, higher N-ratio (> 10%), and <20 bp reads. The resulting clean reads were aligned with the wheat reference genome (version IWGSC) using HISAT2 software and assembled using StringTie software. Gene expression levels were calculated using the transcript per million reads method in RSEM software. Differentially expressed genes (DEGs) between different treatment groups were identified using DESeq2 software, and the screening criteria were p adjust < 0.01 and |Log2 (fold change [FC])| ≥ 1.0. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of these DEGs were conducted at p adjust < 0.05 using Goatools and KOBAS software, respectively. Besides, the protein-protein interaction (PPI) network of DEGs in the key pathways was analyzed based on the STRING database and visualized using the NetworkX package. The raw data were uploaded to the NCBI Sequence Read Archive database under the project accession number PRJNA1020113.
Six DEGs (3 up-regulated and 3 down-regulated) were randomly selected to verify the accuracy of the transcriptome data using qPCR, as previously described (Yue et al., 2021). All primer sequences are listed in Supplementary Table S1.
2.7 Leaf non-targeted metabolomics and data analysis
Briefly, 50 mg frozen leaf samples (n = 6 per group) were homogenized with 400 µL 80% methanol (containing 0.02 mg·mL-1 L-2-chlorophenylalanine) for 6 min at -10°C in a Wonbio-96c tissue grinder (Wonbio, Shanghai, China). The homogenized solutions were ultrasonicated for 30 min at 5°C by an SBL-10DT ultrasonic cleaner (Scientz, Ningbo, China). They were then maintained at -20°C for 30 min and centrifuged at 13,000× g for 15 min at 4°C. Finally, the supernatant was transferred to an injection vial for ultra-high performance liquid chromatography tandem mass spectrometry (UPLC-MS/MS) analysis. Detailed UPLC-MS/MS procedures are provided in the Supplementary Methods and Supplementary Tables S2, S3.
The raw UPLC-MS/MS data were pretreated using Progenesis QI 2.3 software (Waters, Milford, USA), and a data matrix containing the retention time, mass-to-charge (m/z) ratio, and peak intensity was generated. Subsequently, the data matrix was uploaded to the Majorbio Cloud platform (https://cloud.majorbio.com) for data analysis. Principal component analysis (PCA) and orthogonal projections to latent structure-discriminant analysis (OPLS-DA) were performed using the ropls (Version 1.6.2) R package. Differentially abundant metabolites (DAMs) between different treatments were identified based on variable importance in the projection (VIP) > 1 (obtained by the OPLS-DA model) and p < 0.05 (generated by Student’s test). Pathway enrichment analysis of DAMs was conducted using the KEGG database. The raw data were uploaded to the EMBL-EBI MetaboLights database with the accession number MTBLS9697.
2.8 Statistical analysis
Data are expressed as means ± standard deviation (SD). Significant differences between treatment groups were assessed by one-way analysis of variance (ANOVA) with Tukey’s multiple comparisons test using SPSS statistics software 19.0 (IBM, Armonk, USA). In addition, Pearson’s correlation analysis was also performed using SPSS 19.0. p < 0.05 and p < 0.01 were considered statistically significant and extremely significant, respectively.
3 Results
3.1 Effect of Cd, NaCl, and combined stresses on wheat seedling growth parameters
The growth of wheat seedlings was significantly inhibited under Cd, NaCl, and combined stresses compared with that in the control group (Figure 1A). Specifically, primary root length (PRL) significantly decreased in the presence of Cd and NaCl+Cd (Figure 1B). Shoot height (SH) and plant dry weight (PDW) significantly decreased under all three stress treatments (Figures 1C, D). Notably, combined stress showed more severe inhibition of PRL (except for Cd vs. NaCl+Cd), SH, and PDW than Cd or NaCl stress (Figures 1B-D).
Figure 1. Adverse effects of Cd, NaCl, and NaCl+Cd stresses on wheat seedling growth. (A) Photograph of wheat seedlings; (B) Primary root length of wheat plants; (C) Shoot height of wheat plants; (D) Dry weights of wheat plants. Data are expressed as the means ± SD (n = 15 per group). * and ** represent p < 0.05 and p < 0.01, respectively.
3.2 Leaf transcriptomic analysis under Cd, NaCl, and combined stresses
To analyze the responses of global gene expression in leaves to Cd, NaCl, and combined stresses, RNA-seq was performed; an overview of the RNA-seq results is presented in Supplementary Tables S4-S6. In total, 7,178 (5,076 upregulated and 2,102 downregulated), 17,745 (10,616 upregulated and 7,129 downregulated), and 12,000 (7,833 upregulated and 4,167 downregulated) DEGs were identified in the control vs. Cd, control vs. NaCl, and control vs. NaCl+Cd comparison groups, respectively (Figure 2A; Supplementary Figure S2). Among these DEGs, 2,551 genes (10.80%) were common to all stress groups, whereas 2,945 (12.47%), 7,811 (33.07%), and 2,108 (8.93%) genes were unique to the Cd, NaCl, and combined stress groups, respectively (Figure 2B). qPCR analysis showed similar results to those obtained from the RNA sequencing, indicating that the transcriptome data is reliable (Supplementary Figure S3).
Figure 2. Transcriptomic analysis of wheat leaf under Cd, NaCl, and NaCl+Cd stresses. (A) Statistical graph of DEGs under Cd, NaCl, and combined stresses compared to those in the Control group. Red bars represent upregulated genes, and blue bars represent downregulated genes. (B) UpSet plot of DEGs. (C-E) GO enrichment analysis of DEGs in response to Cd (C), NaCl (D), and combined (E) stresses. The size of the dots represents the number of DEGs enriched in the GO term, and p adjust < 0.05 represents statistical significance. (F) Venn diagram of the top 20 significantly enriched KEGG pathways among the three stress treatments.
GO enrichment analysis revealed that these DEGs in response to Cd, NaCl, and combined stresses were significantly enriched in 842, 1,371, and 1,110 GO terms, respectively (Supplementary Tables S7, S8). The top 20 GO terms are shown in Figures 2C-E. Notably, the GO terms related to Na+ transport, including “sodium ion transport”, “sodium:proton antiporter activity”, “sodium ion transmembrane transporter activity”, “sodium ion import across plasma membrane”, and “sodium ion transmembrane transport” were significantly enriched under NaCl and combined stresses. In addition, the GO terms “response to metal ion”, “metal ion transport”, “metal ion binding”, and “metal ion transmembrane transporter activity” were significantly enriched under three stresses. KEGG enrichment analysis revealed that the DEGs in response to Cd, NaCl, and combined stresses were significantly enriched in 44, 45, and 46 pathways, respectively (Supplementary Tables S9, S10). Among the top 20 pathways, 7 were common to all three stresses, while 9, 8, and 8 were unique to Cd, NaCl, and combined stresses, respectively (Figure 2F; Supplementary Figure S4). A detailed description is provided in the Supplementary Results.
3.3 Effect of Cd, NaCl, and combined stresses on Na+ and Cd2+ uptake in wheat seedlings
The Na+ content in the roots and shoots significantly increased under NaCl and combined stresses compared with those in the control group (Figures 3A, B). Notably, Na+ content in the roots under combined stress was 2.17-fold higher than that under NaCl stress. Cd2+ was not detected in the control and NaCl groups. The Cd2+ content in the roots and shoots significantly decreased under combined stress compared with that under Cd stress (Figures 3C, D).
Figure 3. Root Na+ (A), shoot Na+ (B), root Cd2+ (C), and shoot Cd2+ (D) contents in wheat seedlings under Cd, NaCl, and NaCl+Cd stresses. Data are expressed as the means ± SD (n = 3 per group). ** represent p < 0.01.
A total of 25 and 20 DEGs were significantly enriched in GO terms related to Na+ transport under NaCl and combined stresses, respectively. Among these DEGs, some encoding salt overly sensitive 1 (SOS1), SOS2, and SOS3 were significantly upregulated under both stress treatments (Table 1). Remarkably, the relative expression level of these upregulated SOS1 genes under combined stress was only 0.74-fold that of NaCl stress alone. In addition, 7 and 3 DEGs encoding Na+/H+ antiporters (NHXs) were significantly upregulated under NaCl and combined stresses, respectively.
Cd and combined stresses affected many DEGs related to Cd uptake and transport (Table 2). Under Cd stress, 5 DEGs encoding natural resistance-associated macrophage protein (Nramp) and 8 DEGs encoding zinc transporter (ZIP) were significantly upregulated. Under combined stress, 4 DEGs encoding Nramp (2 upregulated and 2 downregulated) and 20 DEGs (7 upregulated and 13 downregulated) encoding ZIP were significantly altered.
3.4 Leaf metabolic analysis under Cd, NaCl, and combined stresses
PCA revealed that the samples from the control, Cd, NaCl, and NaCl+Cd groups were separated, indicating that the metabolic profiles in the four groups were significantly different (Figures 4A, B). DAMs analysis based on VIP > 1.0 and p < 0.05 showed that 556 (375 increase and 181 decrease), 384 (220 increase and 177 decrease), and 397 (169 increase and 215 decrease) DAMs were identified under Cd, NaCl, and combined stresses, respectively, compared with those in the control group (Figure 4C). Among these DAMs, 142 were observed under all stress treatments, and 269, 62, and 54 were unique to Cd, NaCl and combined stresses, respectively (Figure 4D). KEGG pathway enrichment analysis of DAMs demonstrated that 25, 22, and 22 pathways were significantly enriched under Cd, NaCl, and combined stresses, respectively (Figure 4E). Among these, 12 pathways, including “Alanine, aspartate and glutamate metabolism”, “Phenylalanine, tyrosine and tryptophan biosynthesis”, “Lysine biosynthesis”, and “Glycine, serine and threonine metabolism”, were shared under all three stresses (Figures 4F-H). In addition to these common pathways, Cd stress specifically affected 11 other pathways, including “Phenylpropanoid biosynthesis”, “Riboflavin metabolism”, “Arginine biosynthesis”, and “Cyanoamino acid metabolism” (Figure 4F). NaCl stress specifically interfered with “Galactose metabolism” and “Cysteine and methionine metabolism” (Figure 4G). Unlike single stress, “Citrate cycle (TCA cycle)” and “Carbon fixation in photosynthetic organisms” were the unique pathways under combined stress (Figure 4H).
Figure 4. Metabolic profiling of wheat leaf subjected to Cd, NaCl, and NaCl+Cd stresses. (A, B) PCA score plot of metabolites obtained by positive and negative ion modes. (C) Statistical graph of DAMs under Cd, NaCl, and combined stresses compared with those in the control group. Red bars represent upregulated metabolites, and blue bars represent downregulated metabolites. (D) UpSet plot of DAMs. (E) Venn diagram of significantly enriched KEGG pathways among Cd, NaCl, and NaCl+Cd stressed groups. (F-H) KEGG enrichment analysis of DAMs in response to Cd (F), NaCl (G), and NaCl+Cd (H) stresses. The size of the dots represents the number of DAMs enriched in the pathway, and p < 0.05 represents statistical significance.
3.5 Effect of Cd, NaCl, and combined stresses on Chl content, Chl pathway, and photosynthetic proteins
The total Chl, Chl a, and Chl b contents in the leaves were significantly reduced under Cd, NaCl, and combined stresses, compared with those in the control group (Figures 5A-C). Notably, wheat plants under combined stress showed higher inhibition of these indicators than those under Cd or NaCl stress. Pearson’s correlation analysis showed that the total Chl, Chl a, and Chl b contents were positively correlated with PDW (Supplementary Table S11).
Figure 5. Chlorophyll content, chlorophyll-related gene expression, and photosynthesis-related gene expression in wheat seedlings under Cd, NaCl, and NaCl+Cd stresses. (A) Total Chl content. (B) Chl a content. (C) Ch b content. (D) Chl biosynthesis and degradation pathway. (E) Relative abundance of glutamate. (F) Photosynthesis pathway. Red represents upregulated DEGs, blue represents downregulated DEGs, and green represents up- and downregulated DEGs. Data are expressed as the means ± SD (n = 3 per group). * and ** represent p < 0.05 and p < 0.01, respectively.
At the transcriptional level, NaCl stress significantly downregulated the expression of DEGs encoding glutamate-tRNA ligase (GluRS), glutamate-tRNA reductase (GluTR), Mg-chelatase (MgCh), NADPH-protochlorophyllide oxidoreductase (POR), and divinyl chlorophyllide a 8-vinyl-reductase (DVR) in the Chl biosynthesis pathway and upregulated the expression of DEGs encoding Chl b reductase (CBR), chlorophyllase (CLH), pheophytinase (PPH), pheophorbide a oxygenase (PAO), RCC reductase (RCCR), and TIC55 in the Chl degradation pathway (Figure 5D). NaCl+Cd stress significantly upregulated CBR, CLH, PPH, PAO, and TIC55 expressions in the Chl degradation pathway. Cd stress significantly upregulated the CLH expression but significantly downregulated the CBR, PAO, and RCCR expressions, and significantly reduced L-glutamate content at the metabolic level (Figure 5E). Moreover, 22 DEGs encoding 16 subunit proteins, 64 DEGs encoding 26 subunit proteins, and 22 DEGs encoding 11 subunit proteins related to the photosynthesis pathway were affected by Cd, NaCl, and combined stresses, respectively (Figure 5F). Specifically, 11 subunits related to photosystem II (PSII), PSI, and F-type ATPase were significantly downregulated under Cd stress. 21 subunits related to PSII, PSI, photosynthetic electron transport, and F-type ATPase were significantly downregulated under NaCl stress. 6 subunits related to PSII, PSI, photosynthetic electron transport, and F-type ATPase were significantly downregulated under combined stress.
3.6 Effect of Cd, NaCl and combined stresses on O2•− and MDA contents
O2•− and MDA contents in the leaf were significantly increased under Cd, NaCl, and combined stresses compared with those in the control group (Figures 6A, B). Meanwhile, wheat plants under combined stress showed higher O2•− and MDA contents than those under Cd or NaCl stress. Pearson’s correlation analysis showed that O2•− and MDA contents were negatively correlated with PDW (Supplementary Table S11).
Figure 6. Contents of leaf superoxide anion (A) and malondialdehyde (B) in wheat seedlings under Cd, NaCl, and NaCl+Cd stresses. Data are expressed as the means ± SD (n=3 per group). * and ** represent p < 0.05 and p < 0.01, respectively.
3.7 Integrative analysis of transcriptomics and metabolomics data
Integrated transcriptome and metabolome analyses were performed to identify the critical pathways in which the DEGs and DAMs participate. There were 12, 14, and 11 KEGG pathways that simultaneously enriched DEGs and DAMs in control vs. Cd, control vs. NaCl, and control vs. NaCl+Cd comparison groups, respectively (Supplementary Figure S5). In each comparison group, 4 common pathways were detected, with the “Phenylalanine, tyrosine and tryptophan biosynthesis” pathway exhibiting the most significant degree of enrichment (Figure 7). Furthermore, 7, 5, and 1 specific pathways were detected in control vs. Cd, control vs. NaCl, and control vs. NaCl+Cd comparisons, respectively (Figure 7). Based on the p < 0.01 and number of enriched DAMs and DEGs, “Phenylpropanoid biosynthesis” was the most representative unique pathway in the control vs. Cd comparison group (Figure 7A). “Purine metabolism” was the most representative unique pathway in the control vs. NaCl comparison group (Figure 7B). “TCA cycle” was the only unique pathway in the control vs. NaCl+Cd comparison group (Figure 7C). In addition, 5 common pathways were found between the NaCl stress and combined stress groups, indicating that the response of wheat to combined stress was highly similar to that under NaCl stress. Among these common pathways, “Starch and sucrose metabolism” was the most representative and enriched the most significant number of DAMs and DEGs (Figures 7B, C).
Figure 7. Integrative analysis of transcriptomics and metabolomics. (A) KEGG pathway simultaneously enriched DEGs and DAMs in Control vs. Cd. (B) KEGG pathway simultaneously enriched DEGs and DAMs in Control vs. NaCl. (C) KEGG pathway simultaneously enriched DEGs and DAMs in Control vs. NaCl+Cd. Red font represents common pathways, and green font represents unique pathways.
3.8 Expression patterns of DEGs and DAMs associated with phenylalanine, tyrosine and tryptophan biosynthesis and starch and sucrose metabolism
Under Cd, NaCl and combined stresses, many genes were significantly altered in the phenylalanine, tyrosine and tryptophan (Trp) biosynthesis pathway, especially those involved in Trp biosynthesis (Figure 8A; Supplementary Table S12). Specifically, the DEGs encoding AroF, AroK, and TrpF were significantly upregulated under all three stress treatments. Other DEGs encoding TrpB, TrpC, TrpD, TrpE/G, AroA, AroB, and AroC were significantly upregulated under one or two stress treatments. At the metabolic level, the Trp content increased significantly under all three stresses. The 3-dehydro-shikimate, shikimate, and 3-dehydroquinate contents were significantly increased under two stresses.
Figure 8. Effect of Cd, NaCl, and NaCl+Cd stresses on the genes and metabolites associated with phenylalanine, tyrosine and tryptophan biosynthesis, and starch and sucrose metabolism. (A) Phenylalanine, tyrosine and tryptophan biosynthesis pathway during Cd, NaCl, and NaCl+Cd stresses. (B) Starch and sucrose metabolism pathway during NaCl and NaCl+Cd stresses. Red represents upregulated DEGs or increased DAMs. Blue represents downregulated DEGs or decreased DAMs. Green represents up- and downregulated DEGs.
As shown in Figure 8B and Supplementary Table S13, the expressions of some major genes involved in starch and sucrose metabolism were significantly altered by NaCl and combined stresses. Notably, most DEGs showed consistent expression patterns under these two treatments, such as sucrose synthase (SUS), sucrose-6-phosphatase (SPP), trehalose-6-phosphate synthase (TPS), trehalose-6-phosphate phosphatase (TPP), and trehalose (TREH). At the metabolic level, sucrose and trehalose contents significantly increased, while dextran content significantly decreased.
3.9 Analysis of unique pathways under Cd, NaCl, and combined stresses
In the phenylpropanoid biosynthesis pathway, 134 and 11 DEGs were significantly upregulated and downregulated, respectively, during Cd stress (Figure 9A; Supplementary Table S14). Specifically, the expression of genes encoding phenylalanine/tyrosine ammonia-lyase (PTAL), trans-cinnamate 4-monooxygenase (CYP73A), 4-coumarate-CoA ligase (4CL), 5-O-(4-coumaroyl)-D-quinate 3’-monooxygenase (C3’H), caffeoyl-CoA O-methyltransferase (CCoAOMT), and cinnamyl-alcohol dehydrogenase (CAD) was significantly upregulated. Meanwhile, the expression of genes encoding phenylalanine ammonia-lyase (PAL, 12 upregulated and 4 downregulated), caffeate O-methyltransferase (COMT, 6 upregulated and 1 downregulated), ferulate-5-hydroxylase (F5H, 7 upregulated and 1 downregulated), cinnamoyl-CoA reductase (CCR, 6 upregulated and 1 downregulated), shikimate O-hydroxycinnamoyltransferase (HCT, 1 upregulated and 1 downregulated), and peroxidase (POD, 76 upregulated and 2 downregulated) was also significantly affected. In addition, the levels of some pivotal metabolites, such as caffeic acid, ferulic acid, and sinapyl alcohol, were significantly increased.
Figure 9. Expression profiles of genes and metabolites associated with phenylpropanoid biosynthesis, purine metabolism, and TCA cycle. (A) Phenylpropanoid biosynthesis pathway under Cd stress. (B) Purine metabolism pathway under NaCl stress. (C) TCA cycle pathway under NaCl+Cd stress. Red represents upregulated DEGs or increased DAMs. Blue represents downregulated DEGs or decreased DAMs. Green represents up- and downregulated DEGs.
Under NaCl stress, 74 DEGs and 6 DAMs were significantly altered in the purine metabolism pathway (Figure 9B; Supplementary Table S14). The levels of all DAMs were significantly increased, including GMP, guanosine, inosine, adenine, deoxyadenosine and uric acid. Furthermore, there was a significant upregulation of the expression of genes involved in these DAMs, such as nucleoside-diphosphate kinase (NDK), ribonucleoside-diphosphate reductase (RNR), 5’-nucleotidase (NT5C), AMP deaminase (AMPD), and xanthine dehydrogenase/oxidase (XDH).
Unlike under single stress treatment, combined stress explicitly affected the TCA cycle (Figure 9C; Supplementary Table S14). At the transcriptional level, the expression of genes encoding citrate synthase (CS), aconitate hydratase (ACO), 2-oxoglutarate dehydrogenase (2-OGDH), dihydrolipoyllysine-residue succinyltransferase (DLST), dihydrolipoyl dehydrogenase (DLD), and phosphoenolpyruvate carboxykinase (PEPCK), was significantly upregulated. Furthermore, the expression of genes encoding succinate dehydrogenase (SDH, 10 upregulated and 2 downregulated) and pyruvate dehydrogenase E1 component (PDH, 5 upregulated and 1 downregulated) was remarkably altered. At the metabolic level, isocitrate content was significantly increased, while cis-aconitate and malate contents were significantly decreased.
In the PPI network, COMT (TraesCS3B02G448600), CAD (TraesCS6D02G162800, TraesCS6B02G201600), and 4CL (TraesCS6B02G179900, TraesCS6A02G151700, TraesCS2D02G150400, TraesCS6D02G141700) in the phenylpropanoid biosynthesis pathway, AMPD (TraesCS2D02G084300, TraesCS2A02G087000, TraesCS2B02G101700) and ADK (TraesCS4B02G324400, TraesCS7A02G382500, TraesCS2B02G579200) in the purine metabolism pathway, and CS (TraesCS6B02G194700, TraesCS6D02G155800), SDH (TraesCS5A02G189800, TraesCS7A02G314000, TraesCS7D02G311200), and DLST (TraesCS2D02G039500) in the TCA cycle had the highest scores for betweenness centrality, thus reflecting their functional importance in the network (Supplementary Figure S6).
4 Discussion
4.1 Combined stress results in more severe growth inhibition and higher O2•− and MDA levels than Cd or NaCl stress
Globally, soil salinity stress and Cd pollution have become non-negligible abiotic stresses. The retardation of crop growth due to Cd and salinity has been observed in wheat, maize, rice, and other crops (Huang et al., 2021; Sun et al., 2022; Zhang et al., 2023). Consistent with the results, in the present study, we found significant decreases in PRL, SH, and PDW in Cd- and NaCl-stressed wheat seedlings. Wheat plants under combined stress exhibited lower PRL, SH and PDW than those under single stress, indicating that combined stress caused more severe damage to wheat growth. Similar results have been observed in other wheat cultivars, quinoa, and Spartina alterniflora exposed to combined Cd and salinity stress (Abdal et al., 2023; Chai et al., 2013; Shafi et al., 2009). Generally, abiotic stress causes the accumulation of large amounts of reactive oxygen species (ROS) in plants, resulting in severe oxidative damage to proteins, DNA, and lipids (Das and Roychoudhury, 2014). Our study observed significant increases in O2•− and MDA contents under all stress treatments, with combined stress having the highest contents. O2•− is usually the first formed ROS, and can be converted into more reactive and toxic OH• and 1O2, leading to membrane lipid peroxidation (Halliwell, 2006). Therefore, these results indicated that combined stress had more serious oxidative damage than those of NaCl or Cd stress, which may be one of the reasons for the more severe growth inhibition. The positive correlation among PDW, O2•−, and MDA further confirmed this.
4.2 Cd increases Na+ content, whereas NaCl reduces Cd2+ content
Salinity stress usually results in Na+ accumulation in plants. In this study, high concentrations of Na+ were observed in NaCl-treated wheat plants. Plants have evolved a series of salinity-resistance strategies to adapt to salinity stress. Among these, the SOS pathway is one of the most important, and involves three proteins: plasma membrane Na+/H+ antiporter SOS1, calcineurin B-like (CBL)-interacting protein kinase SOS2, and calcium-binding protein SOS3 (Zhu, 2000). When plants sense excess Na+, SOS3 activates and recruits SOS2 to the plasma membrane. SOS2 phosphorylates and activates SOS1 to exclude excessive Na+ from the cytosol (Ali et al., 2023; Ji et al., 2013). In this study, the SOS1, SOS2, and SOS3 transcriptional levels were significantly upregulated under NaCl and combined stresses, indicating that wheat plants activated the SOS pathway under these two stresses. Notably, the SOS1 upregulation under the combined stress was only 0.74-fold that under salinity stress alone, indicating that Cd may weaken the SOS1 expression, thereby inhibiting Na+ efflux and increasing Na+ content under combined stress. In addition to SOS-mediated Na+ efflux, plants transport Na+ into the vacuoles for sequestration by the vacuole-localized NHXs under salt stress (Xiao and Zhou, 2023). In this study, some DEGs encoding vacuolar NHXs were significantly upregulated under NaCl and combined stresses, indicating that wheat plants initiated vacuolar compartmentalization of Na+.
In contrast to the Na+ content, the Cd2+ content was significantly lower under combined stress than that under Cd stress. Plants contain many proteins involved in Cd uptake and transport. Among them, Nramp and ZIP are responsible for Cd uptake and translocation to aboveground parts (Sasaki et al., 2012; Hou et al., 2023). In this study, 5 DEGs encoding Nramp and 8 DEGs encoding ZIPs were significantly upregulated under Cd stress, whereas 2 DEGs encoding Nramp and most DEGs encoding ZIP were significantly downregulated under combined stress. Therefore, the low Cd accumulation under combined stress may be attributed to the partial downregulation of Nramp and ZIP gene expression by NaCl.
4.3 Cd, NaCl, and combined stresses inhibit photosynthetic pigments and photosynthetic proteins
Chl is essential for photosynthesis in most plants, where Chl a converts light energy into chemical energy, and Chl b captures and transmits light energy (Trukhachev et al., 2022). In this study, Cd, NaCl, and combined stresses sharply suppressed the total Chl, Chl a, and Chl b contents. In plants, Chl biosynthesis commences with the precursor 5-aminolevulinic acid (ALA), which is generated from glutamate under catalysis by GluRS, GluTR and glutamate-1-semialdehyde transaminase (Tanaka et al., 2011). ALA is then gradually converted into protoporphyrin IX (Proto IX), which is further catalyzed by the enzymes MgCh, POR, DVR, Chl synthetase (ChlS), and chlorophyllide an oxygenase (CAO) to produce Chl a and Chl b (Masuda and Takamiya, 2004; Liu and Zheng, 2008). In this study, we observed that many DEGs involved in the Chl biosynthesis pathway under NaCl stress were significantly downregulated, including GluRS, GluTR, MgCh, POR, and DVR, indicating that NaCl stress inhibited Chl biosynthesis. Contrastingly, many DEGs involved in Chl degradation, such as CBR, CLH, PPH, PAO, and TIC55, were significantly upregulated under NaCl and NaCl+Cd stresses. In the Chl degradation pathway, Chl b is first converted to Chl a through a two-step reaction catalyzed by CBR and hydroxymethyl chlorophyll a reductase (HCAR) (Horie et al., 2009; Meguro et al., 2011). Chl a is converted into Pheide a by removing phytol and Mg2+ catalyzed by Mg-dechelatase (MDCase) and PPH or CLH enzymes (Lin and Charng, 2021). Pheide a undergoes ring opening catalyzed by PAO to generate linear red chlorophyll metabolites (RCC), which are further converted into colorless pFCC and hpFCC under the catalysis of RCC reductase (RCCR) and TIC55 (Bhattacharyya et al., 2015; Kuai et al., 2018). Therefore, we suggested that NaCl stress reduced the Chl content by inhibiting Chl biosynthesis and activating Chl degradation, whereas NaCl+Cd reduced the Chl content by activating Chl degradation. Unlike NaCl and NaCl+Cd stresses, Cd stress significantly downregulated the expression of some DEGs related to Chl degradation, such as CBR, PAO, and RCCR. Notably, we found that the L-glutamate content, the initial substrate for Chl biosynthesis, was significantly decreased under Cd stress, which may partially explain why Cd stress reduced Chl content.
Photosynthesis plays a decisive role in plant growth and development, mainly involving four pigment-protein complexes: photosystem I (PSI), photosystem II (PSII), cytochrome b6f complex, and ATP synthase complex (Allen et al., 2011). PSI and PSII are multi-subunit membrane protein complexes, where PSII catalyzes water splitting, oxygen evolution, and plastoquinone reduction, whereas PSI generates reducing power for the reduction of NADP+ to NADPH (Gao et al., 2018). ATP synthase participates in ATP synthesis, mediated by transmembrane proton motive force (Zhang et al., 2015). In this study, multiple genes encoding PSII, PSI, and F-type ATPase were significantly downregulated under Cd, NaCl, and combined stresses, indicating the inhibition of photosynthetic proteins, which may further lead to a decline in photosynthesis. In summary, the stress-induced inhibition of photosynthetic pigments and proteins in this study may affect photosynthesis and carbon fixation, ultimately resulting in growth inhibition in wheat seedlings.
4.4 Cd, NaCl, and combined stresses activate Trp biosynthesis and increase Trp content
Under all three stresses, there was a significant upregulation of DEGs related to Trp biosynthesis, indicating that Cd, NaCl, and combined stresses activated Trp biosynthesis. The observed increase in Trp content under all three stresses further supported our hypothesis. Trp is an essential amino acid for protein synthesis in living organisms and has recently been shown to play an important role in protecting plants against abiotic stress (Godoy et al., 2021; Zemanová et al., 2014). For instance, exogenous Trp improved Cd tolerance in broccoli (Brassica oleracea var. italica) and Arabidopsis (Arabidopsis thaliana) (Jiang et al., 2022; Li and Qi, 2023). The application of Trp alleviated the detrimental effects of salt stress in red pepper (Capsicum annuum L.) and sugar beet (Beta vulgaris L.) (Hozayn et al., 2020; Jamil et al., 2018). Therefore, activating Trp biosynthesis may improve wheat adaptation to Cd, NaCl, and combined stresses.
4.5 Starch and sucrose metabolism plays an important role under NaCl and combined stresses
In this study, numerous DEGs involved in starch and sucrose metabolism were significantly altered, resulting in increased sucrose and trehalose contents under NaCl and combined stresses. Similarly, increased levels of soluble sugars (sucrose, trehalose, glucose, and fructose) have been observed in many plants under salinity, drought, water, chilling, and nutrient deficiency stresses (Martínez-Noël and Tognetti, 2018; Saddhe et al., 2021). Accumulated soluble sugars act as osmo-protectants in maintaining osmotic homeostasis, stabilizing membrane structure, and scavenging toxic ROS against different stresses (Amist and Singh, 2020; Singh et al., 2015). Exogenous application of sucrose and trehalose improved plant tolerance to salinity and heat stresses (Gong and Chen, 2021; Yang et al., 2014). Therefore, the accumulation of sucrose and trehalose observed in this study may help wheat better cope with salinity and combined stresses.
4.6 Phenylpropanoid biosynthesis plays a unique role under Cd stress
For phenylpropanoid biosynthesis, multitudinous DEGs encoding PTAL, CYP73A, 4CL, C3’H, CCoAOMT, CAD, COMT, F5H, CCR, and POD were significantly upregulated under Cd stress, indicating that Cd stress activated the pathway. As one of the most important metabolic pathways in plants, the phenylpropanoid biosynthesis pathway produces many secondary metabolites, including phenolic acids possessing redox properties (Zhang and Liu, 2015; Sharma et al., 2019). In this study, there were significant increases in caffeic acid (CA) and fumaric acid (FA) under Cd stress. CA is synthesized from phenylalanine or tyrosine under the sequential catalysis of PAL/PTAL, CYP73A, and C3’H, and then transformed into FA via COMT catalysis (Riaz et al., 2019). They are well-known as potential antioxidants that improve abiotic stress resistance in plants by scavenging ROS (Chen et al., 2020; Shu et al., 2022; Wan et al., 2015). Therefore, wheat may increase CA and FA levels to cope with Cd-induced oxidative damage by activating the phenylpropanoid biosynthesis pathway.
4.7 Purine metabolism plays a unique role under NaCl stress
Under NaCl stress, all DAMs, including uric acid, and nearly all DEGs in the purine metabolism pathway were markedly increased or upregulated in the present study, suggesting that the purine metabolism pathway plays a positive role in the response to NaCl stress. Our results were consistent with those of previous studies showing that the activation of enzymes and increased levels in intermediary metabolites are associated with purine metabolism in stressed plants (Alamillo et al., 2010; Liu et al., 2020; Rose et al., 2012; Sagi et al., 1998). Purine metabolism not only recycles nitrogen by disintegrating the purine ring, but also produces uric acid which functions to scavenge ROS (Ma et al., 2016; Soltabayeva et al., 2018; Werner and Witte, 2011). Sun et al. (2021) reported that exogenous application of uric acid elevated the resistance of apple plants to salinity stress by enhancing ROS scavenging. Therefore, we suggested that wheat plants may increase intermediate metabolites, including uric acid, by activating the purine metabolism pathway to protect plants from oxidative damage due to NaCl stress.
4.8 TCA cycle plays a unique role under combined stress
Under combined stress, the “TCA cycle” was the only specific pathway observed. As the epicenter of cell metabolism, the TCA cycle is essential in energy and material metabolism (Zhang and Fernie, 2023). In this study, critical genes related to the TCA cycle, such as CS, ACO, 2-OGDH, DLST, DLD, and SDH, were consistently upregulated when subjected to combined stress. CS initiates the first step of the TCA cycle and catalyzes the condensation of acetyl-CoA and oxaloacetate to form citrate, which is then converted to isocitrate under the action of ACO (Lehmann et al., 2009; Schmidtmann et al., 2014). 2-OGDH, DLST, and DLD are components of the 2-OGDH complex, and convert 2-oxoglutarate into succinyl CoA and generate further CO2 and NADH (Araújo et al., 2012). SDH participates in the oxidation of succinate to fumarate with the simultaneous reduction of flavin adenine dinucleotide (FAD) to FADH2 (Huang and Millar, 2013). NADH and FADH2 molecules pass their electrons into the electron transport chain and undergo oxidative phosphorylation to generate most of the energy (ATP) required for plant growth, cell division, and other physiological processes (Ghifari et al., 2023). Therefore, the upregulation of genes involved in the TCA cycle may produce more energy provision for wheat plants to counteract the combined stress.
5 Conclusions
Our study showed that combined stress caused more severe adverse effects on wheat growth and Chl content and produced more O2•− and MDA contents than those under NaCl or Cd stress. Stress-induced decrease in Chl levels may be related to the inhibition of Chl biosynthesis, activation of Chl degradation, or a decline in glutamate content. Integrated transcriptomics and metabolomics analyses revealed similarities and differences in the responses of wheat to single and combined stresses. All three stresses activated Trp biosynthesis, increasing Trp levels. Meanwhile, Cd stress specifically activated phenylpropanoid biosynthesis and increased CA and FA contents. NaCl stress specifically activated purine metabolism and enhanced intermediate metabolite levels. Combined stress explicitly activated the TCA cycle, leading to increased energy provision. In addition, NaCl and combined stresses affected starch and sucrose metabolism, resulting in the accumulation of sucrose and trehalose.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
ZY: Conceptualization, Writing – original draft, Writing – review & editing. YL: Investigation, Methodology, Writing – review & editing. LZ: Investigation, Methodology, Writing – review & editing. QZ: Investigation, Methodology, Writing – review & editing. YW: Investigation, Methodology, Writing – review & editing. YH: Investigation, Methodology, Writing – review & editing. MZ: Data curation, Formal analysis, Writing – review & editing. YC: Data curation, Visualization, Writing – review & editing. ZW: Supervision, Visualization, Writing – review & editing. LH: Funding acquisition, Supervision, Writing – review & editing. KM: Conceptualization, Writing – original draft, Writing – review & editing, Supervision.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the Henan Provincial Department of Science and Technology (No.242102111088; 222102320365).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2024.1378226/full#supplementary-material
References
Abbas, T., Rizwan, M., Ali, S., Adrees, M., Zia-Ur-Rehman, M., Qayyum, M. F., et al. (2018). Effect of biochar on alleviation of cadmium toxicity in wheat (Triticum aestivum L.) grown on Cd-contaminated saline soil. Environ. Sci. pollut. Res. Int. 25, 25668–25680. doi: 10.1007/s11356-017-8987-4
Abdal, N., Abbas, G., Asad, S. A., Ghfar, A. A., Shah, G. M., Rizwan, M., et al. (2023). Salinity mitigates cadmium-induced phytotoxicity in quinoa (Chenopodium quinoa Willd.) by limiting the Cd uptake and improved responses to oxidative stress: implications for phytoremediation. Environ. Geochem Health 45, 171–185. doi: 10.1007/s10653-021-01082-y
Acosta, J. A., Jansen, B., Kalbitz, K., Faz, A., Martinez-Martinez, S. (2011). Salinity increases mobility of heavy metals in soils. Chemosphere 85, 1318–1324. doi: 10.1016/j.chemosphere.2011.07.046
Ahn, S. C., Chang, J. Y., Lee, J. S., Yu, H. Y., Jung, A. R., Kim, J. Y., et al. (2017). Exposure factors of cadmium for residents in an abandoned metal mine area in Korea. Environ. Geochem Health 39, 1059–1070. doi: 10.1007/s10653-016-9872-7
Alamillo, J. M., Dí Az-Leal, J. L., Sánchez-Moran, M. A. V., Pineda, M. (2010). Molecular analysis of ureide accumulation under drought stress in Phaseolus vulgaris L. Plant Cell Environ. 33, 1828–1837. doi: 10.1111/j.1365-3040.2010.02187.x
Ali, A., Petrov, V., Yun, D. J., Gechev, T. (2023). Revisiting plant salt tolerance: novel components of the SOS pathway. Trends Plant Sci. 28, 1060–1069. doi: 10.1016/j.tplants.2023.04.003
Allen, J. F., de Paula, W. B. M., Puthiyaveetil, S., Nield, J. (2011). A structural phylogenetic map for chloroplast photosynthesis. Trends Plant Sci. 16, 645–655. doi: 10.1016/j.tplants.2011.10.004
Amist, N., Singh, N. (2020). “The role of sugars in the regulation of environmental stress,” in Plant life under changing environment (Salt Lake City: Academic Press), 497–512.
Arao, T., Ishikawa, S., Murakami, M., Abe, K., Maejima, Y., Makino, T. (2010). Heavy metal contamination of agricultural soil and countermeasures in Japan. Paddy Water Environ. 8, 247–257. doi: 10.1007/s10333-010-0205-7
Araújo, W. L., Tohge, T., Osorio, S., Lohse, M., Balbo, I., Krahnert, I., et al. (2012). Antisense inhibition of the 2-oxoglutarate dehydrogenase complex in tomato demonstrates its importance for plant respiration and during leaf senescence and fruit maturation. Plant Cell 24, 2328–2351. doi: 10.1105/tpc.112.099002
Arnon, D. (1949). Copper enzymes in isolated chloroplasts. phenoloxidase in Beta vulgaris. Plant Physiol. 24, 1–15. doi: 10.1104/pp.24.1.1
Bhattacharyya, D., Gnanasekaran, P., Kumar, R. K., Kushwaha, N. K., Sharma, V. K., Yusuf, M. A., et al. (2015). A geminivirus betasatellite damages the structural and functional integrity of chloroplasts leading to symptom formation and inhibition of photosynthesis. J. Exp. Bot. 66, 5881–5895. doi: 10.1093/jxb/erv299
Chai, M. W., Shi, F. C., Li, R. L., Liu, F. C., Qiu, G. Y., Liu, L. M. (2013). Effect of NaCl on growth and Cd accumulation of halophyte Spartina alterniflora under CdCl2 stress. S Afr J. Bot. 85, 63–69. doi: 10.1016/j.sajb.2012.12.004
Chen, S., Lin, R., Lu, H., Wang, Q., Yang, J., Liu, J., et al. (2020). Effects of phenolic acids on free radical scavenging and heavy metal bioavailability in Kandelia obovata under cadmium and zinc stress. Chemosphere 249, 126341. doi: 10.1016/j.chemosphere.2020.126341
Cheng, M., Kopittke, P. M., Wang, A., Tang, C. (2019). Salinity decreases Cd translocation by altering Cd speciation in the halophytic Cd-accumulator Carpobrotus rossii. Ann. Bot. 123, 121–132. doi: 10.1093/aob/mcy148
Das, K., Roychoudhury, A. (2014). Reactive oxygen species (ROS) and response of antioxidants as ROS-scavengers during environmental stress in plants. Front. Environ. Sci. 2, 53. doi: 10.3389/fenvs.2014.00053
Gao, J., Wang, H., Yuan, Q., Feng, Y. (2018). Structure and function of the photosystem supercomplexes. Front. Plant Sci. 9, 357. doi: 10.3389/fpls.2018.00357
Ghifari, A. S., Saha, S., Murcha, M. W. (2023). The biogenesis and regulation of the plant oxidative phosphorylation system. Plant Physiol. 192, 728–747. doi: 10.1093/plphys/kiad108
Godoy, F., Olivos-Hernández, K., Stange, C., Handford, M. (2021). Abiotic stress in crop species: improving tolerance by applying plant metabolites. Plants 10, 186. doi: 10.3390/plants10020186
Gong, H. L., Chen, Q. Q. (2021). Exogenous sucrose protects potato seedlings against heat stress by enhancing the antioxidant defense system. J. Soil Sci. Plant Nutr. 21, 1511–1519. doi: 10.1007/s42729-021-00457-5
Halliwell, B. (2006). Reactive species and antioxidants. Redox biology is a fundamental theme of aerobic life. Plant Physiol. 141, 312–322. doi: 10.1104/pp.106.077073
Hiscox, J. D., Israelstam, G. F. (1979). A method for the extraction of chlorophyll from leaf tissue without maceration. Can. J. Bot. 57, 1332–1334. doi: 10.1139/b79-163
Horie, Y., Ito, H., Kusaba, M., Tanaka, R., Tanaka, A. (2009). Participation of chlorophyll b reductase in the initial step of the degradation of light-harvesting chlorophyll a/b-Protein complexes in Arabidopsis. J. Biol. Chem. 284, 17449–17456. doi: 10.1074/jbc.M109.008912
Hou, L., Ji, S., Zhang, Y., Wu, X., Zhang, L., Liu, P. (2023). The mechanism of silicon on alleviating cadmium toxicity in plants: A review. Front. Plant Sci. 14, 1141138. doi: 10.3389/fpls.2023.1141138
Hozayn, M., Abd-Elmonem, A. A., Samaha, G. M. (2020). The physiological effect of pre-soaking with tryptophan on sugar beet (Beta vulgaris L.) productivity under different levels of salinity stresses. Bull. Natl. Res. Cent 44, 65. doi: 10.1186/s42269-020-00324-w
Hu, J. H., Chen, G. L., Xu, K., Wang, J. (2022). Cadmium in cereal crops: uptake and transport mechanisms and minimizing strategies. J. Agric. Food Chem. 70, 5961–5974. doi: 10.1021/acs.jafc.1c07896
Huang, H., Li, M., Rizwan, M., Dai, Z., Yuan, Y., Hossain, M. M., et al. (2021). Synergistic effect of silicon and selenium on the alleviation of cadmium toxicity in rice plants. J. Hazard Mater 401, 123393. doi: 10.1016/j.jhazmat.2020.123393
Huang, S., Millar, A. H. (2013). Succinate dehydrogenase: the complex roles of a simple enzyme. Curr. Opin. Plant Biol. 16, 344–349. doi: 10.1016/j.pbi.2013.02.007
Jamil, M., Kharal, M. A., Ahmad, M., Abbasi, G. H., Nazli, F., Hussain, A., et al. (2018). Inducing salinity tolerance in red pepper (Capsicum annuum L.) through exogenous application of proline and L-tryptophan. Soil Environ. 37, 160–168. doi: 10.25252/SE/18/31052
Ji, H., Pardo, J. M., Batelli, G., van Ooste, M. J., Bressan, R. A., Li, X. (2013). The salt overly sensitive (SOS) pathway: established and emerging roles. Mol. Plant 6, 275–286. doi: 10.1093/mp/sst017
Jiang, J., Wang, Z., Kong, X., Chen, Y., Li, J. (2022). Exogenous tryptophan application improves cadmium tolerance and inhibits cadmium upward transport in broccoli (Brassica oleracea var. italica). Front. Plant Sci. 13, 969675. doi: 10.3389/fpls.2022.969675
Kuai, B., Chen, J., Hörtensteiner, S. (2018). The biochemistry and molecular biology of chlorophyll breakdown. J. Exp. Bot. 69, 751–767. doi: 10.1093/jxb/erx322
Kumar, V., Sharma, A., Kaur, P., Singh Sidhu, G. P., Bali, A. S., Bhardwaj, R., et al. (2019). Pollution assessment of heavy metals in soils of India and ecological risk assessment: A state-of-the-art. Chemosphere 216, 449–462. doi: 10.1016/j.chemosphere.2018.10.066
Lehmann, M., Schwarzländer, M., Obata, T., Sirikantaramas, S., Burow, M., Olsen, C. E., et al. (2009). The metabolic response of Arabidopsis roots to oxidative stress is distinct from that of heterotrophic cells in culture and highlights a complex relationship between the levels of transcripts, metabolites, and flux. Mol. Plant 2, 390–406. doi: 10.1093/mp/ssn080
Li, H., La, S., Zhang, X., Gao, L., Tian, Y. (2021). Salt-induced recruitment of specific root-associated bacterial consortium capable of enhancing plant adaptability to salt stress. ISME J. 15, 2865–2882. doi: 10.1038/s41396-021-00974-2
Li, Y., Qi, X. (2023). Tryptophan pretreatment adjusts transcriptome and metabolome profiles to alleviate cadmium toxicity in Arabidopsis. J. Hazard Mater 2452, 131226. doi: 10.1016/j.jhazmat.2023.131226
Li, X., Zhou, D. (2019). A meta-analysis on phenotypic variation in cadmium accumulation of rice and wheat: Implications for food cadmium risk control. Pedosphere 29, 545–553. doi: 10.1016/S1002-0160(19)60828-3
Lin, Y. P., Charng, Y. Y. (2021). Chlorophyll dephytylation in chlorophyll metabolism: a simple reaction catalyzed by various enzymes. Plant Sci. 302, 110682. doi: 10.1016/j.plantsci.2020.110682
Liu, L., Wang, B., Liu, D., Zou, C., Wu, P., Wang, Z., et al. (2020). Transcriptomic and metabolomic analyses reveal mechanisms of adaptation to salinity in which carbon and nitrogen metabolism is altered in sugar beet roots. BMC Plant Biol. 20, 138. doi: 10.1186/s12870-020-02349-9
Liu, H. H., Zheng, C. C. (2008). MPEC: an important gene in the chlorophyll biosynthesis pathway in photosynthetic organisms. Photosynthetica 46, 321–328. doi: 10.1007/s11099-008-0060-6
Lutts, S., Lefèvre, I. (2015). How can we take advantage of halophyte properties to cope with heavy metal toxicity in salt-affected areas? Ann. Bot. 115, 509–528. doi: 10.1093/aob/mcu264
Ma, X., Wang, W., Bittner, F., Schmidt, N., Berkey, R., Zhang, L., et al. (2016). Dual and opposing roles of xanthine dehydrogenase in defense-associated reactive oxygen species metabolism in Arabidopsis. Plant Cell 28, 1108–1126. doi: 10.1105/tpc.15.00880
Martínez-Noël, G. M. A., Tognetti, J. A. (2018). “Sugar signaling under abiotic stress in plants,” in Plant metabolites and regulation under environmental stress (Salt Lake City: Academic Press), 397–406.
Masuda, T., Takamiya, K. (2004). Novel insights into the enzymology, regulation and physiological functions of light-dependent protochlorophyllide oxidoreductase in angiosperms. Photosynth Res. 81, 1–29. doi: 10.1023/B:PRES.0000028392.80354.7c
Meguro, M., Ito, H., Takabayashi, A., Tanaka, R., Tanaka, A. (2011). Identification of the 7-hydroxymethyl chlorophyll a reductase of the chlorophyll cycle in Arabidopsis. Plant Cell 23, 3442–3453. doi: 10.1105/tpc.111.089714
Mühling, K. H., Läuchli, A. (2003). Interaction of NaCl and Cd stress on compartmentation pattern of cations, antioxidant enzymes and proteins in leaves of two wheat genotypes differing in salt tolerance. Plant Soil 253, 219–231. doi: 10.1023/A:1024517919764
Nguyen, T. H., Hoang, H. N. T., Bien, N. Q., Tuyen, L. H., Kim, K. W. (2020). Contamination of heavy metals in paddy soil in the vicinity of Nui Phao multi-metal mine, North Vietnam. Environ. Geochem Health 42, 4141–4158. doi: 10.1007/s10653-020-00611-5
Qadir, M., Quillérou, E., Nangia, V., Murtaza, G., Singh, M., Thomas, R. J., et al. (2014). Economics of salt-induced land degradation and restoration. Nat. Resour Forum 38, 282–295. doi: 10.1111/narf.2014.38.issue-4
Quan, X., Liang, X., Li, H., Xie, C., He, W., Qin, Y. (2021). Identification and characterization of wheat germplasm for salt tolerance. Plants 10, 268. doi: 10.3390/plants10020268
Rehman, Z. U., Khan, S., Brusseau, M. L., Shah, M. T. (2017). Lead and cadmium contamination and exposure risk assessment via consumption of vegetables grown in agricultural soils of five-selected regions of Pakistan. Chemosphere 168, 1589–1596. doi: 10.1016/j.chemosphere.2016.11.152
Riaz, U., Kharal, M. A., Murtaza, G., Zaman, Q., Javaid, S., Malik, H. A., et al. (2019). Prospective roles and mechanisms of caffeic acid in counter plant stress: A mini review. Pakistan J. Agric. Res. 32, 8–19. doi: 10.17582/journal.pjar/2019/32.1.8.19
Rose, M. T., Rose, T. J., Pariasca-Tanaka, J., Yoshihashi, T., Neuweger, H., Goesmann, A., et al. (2012). Root metabolic response of rice (Oryza sativa L.) genotypes with contrasting tolerance to zinc deficiency and bicarbonate excess. Planta 236, 959–973. doi: 10.1007/s00425-012-1648-4
Saddhe, A. A., Manuka, R., Penna, S. (2021). Plant sugars: Homeostasis and transport under abiotic stress in plants. Physiol. Plant 171, 739–755. doi: 10.1111/ppl.v171.4
Sagi, M., Omarov, R. T., Lips, S. H. (1998). The Mo-hydroxylases xanthine dehydrogenase and aldehyde oxidase in ryegrass as affected by nitrogen and salinity. Plant Sci. 135, 125–135. doi: 10.1016/S0168-9452(98)00075-2
Sasaki, A., Yamaji, N., Yokosho, K., Ma, J. F. (2012). Nramp5 is a major transporter responsible for manganese and cadmium uptake in rice. Plant Cell 24, 2155–2167. doi: 10.1105/tpc.112.096925
Schmidtmann, E., König, A. C., Orwat, A., Leister, D., Hartl, M., Finkemeier, I. (2014). Redox regulation of Arabidopsis mitochondrial citrate synthase. Mol. Plant 7, 156–169. doi: 10.1093/mp/sst144
Shafi, M., Bakht, J., Hassan, M. J., Raziuddin, M., Zhang, G. (2009). Effect of cadmium and salinity stresses on growth and antioxidant enzyme activities of wheat (Triticum aestivum L.). Bull. Environ. Contam Toxicol. 82, 772–776. doi: 10.1007/s00128-009-9707-7
Sharma, A., Shahzad, B., Rehman, A., Bhardwaj, R., Landi, M., Zheng, B. (2019). Response of phenylpropanoid pathway and the role of polyphenols in plants under abiotic stress. Molecules 24, 2452. doi: 10.3390/molecules24132452
Shu, P., Li, Y., Li, Z., Xiang, L., Sheng, J., Shen, L. (2022). Ferulic acid enhances chilling tolerance in tomato fruit by up-regulating the gene expression of CBF transcriptional pathway in MAPK3-dependent manner. Postharvest Biol. Tec 185, 111775. doi: 10.1016/j.postharvbio.2021.111775
Singh, M., Kumar, J., Singh, S., Singh, V., Prasad, S. (2015). Roles of osmoprotectants in improving salinity and drought tolerance in plants: a review. Rev. Environ. Sci. Bio/Technol 14, 407–426. doi: 10.1007/s11157-015-9372-8
Soltabayeva, A., Srivastava, S., Kurmanbayeva, A., Bekturova, A., Fluhr, R., Sagi, M. (2018). Early senescence in older leaves of low nitrate-grown Atxdh1 uncovers a role for purine catabolism in N Supply. Plant Physiol. 178, 1027–1044. doi: 10.1104/pp.18.00795
Sriprachote, A., Kanyawongha, P., Ochiai, K., Matoh, T. (2012). Current situation of cadmium-polluted paddy soil, rice and soybean in the Mae Sot District, Tak Province, Thailand. Soil Sci. Plant Nutr. 58, 349–359. doi: 10.1080/00380768.2012.686435
Sun, T., Pei, T., Yang, L., Zhang, Z., Li, M., Liu, Y., et al. (2021). Exogenous application of xanthine and uric acid and nucleobase-ascorbate transporter MdNAT7 expression regulate salinity tolerance in apple. BMC Plant Biol. 21, 52. doi: 10.1186/s12870-021-02831-y
Sun, L., Zhang, X., Ouyang, W., Yang, E., Cao, Y., Sun, R. (2022). Lowered Cd toxicity, uptake and expression of metal transporter genes in maize plant by ACC deaminase-producing bacteria Achromobacter sp. J. Hazard Mater 423, 127036. doi: 10.1016/j.jhazmat.2021.127036
Tanaka, R., Kobayashi, K., Masuda, T. (2011). Tetrapyrrole metabolism in Arabidopsis thaliana. Arabidopsis Book 9, e0145. doi: 10.1199/tab.0145
Thomas, E., Atkinson, R., Zavaleta, D., Rodriguez, C., Lastra, S., Yovera, F., et al. (2023). The distribution of cadmium in soil and cacao beans in Peru. Sci. Total Environ. 881, 163372. doi: 10.1016/j.scitotenv.2023.163372
Tricker, P. J., ElHabti, A., Schmidt, J., Fleury, D. (2018). The physiological and genetic basis of combined drought and heat tolerance in wheat. J. Exp. Bot. 69, 3195–3210. doi: 10.1093/jxb/ery081
Trukhachev, V. I., Seregina, I. I., Belopukhov, S. L., Dmitrevskaya, I. I., Fomina, T. I., Zharkikh, O. A., et al. (2022). The effect of stressful ecological conditions on chlorophyll content in the leaves of spring wheat plants. IOP Conf Ser. Earth Environ. Sci. 981, 32093. doi: 10.1088/1755-1315/981/3/032093
Uchimiya, M., Bannon, D., Nakanishi, H., McBride, M. B., Williams, M. A., Yoshihara, T. (2020). Chemical speciation, plant uptake, and toxicity of heavy metals in agricultural soils. J. Agric. Food Chem. 68, 12856–12869. doi: 10.1021/acs.jafc.0c00183
Wan, Y., Zhang, Y., Zhang, L., Zhou, Z., Li, X., Shi, Q., et al. (2015). Caffeic acid protects cucumber against chilling stress by regulating antioxidant enzyme activity and proline and soluble sugar contents. Acta Physiol. Plant 37, 1706. doi: 10.1007/s11738-014-1706-6
Werner, A. K., Witte, C. P. (2011). The biochemistry of nitrogen mobilization: purine ring catabolism. Trends Plant Sci. 16, 381–387. doi: 10.1016/j.tplants.2011.03.012
Wu, L., Yue, W., Zheng, N., Guo, M., Teng, Y. (2022). Assessing the impact of different salinities on the desorption of Cd, Cu and Zn in soils with combined pollution. Sci. Total Environ. 836, 155725. doi: 10.1016/j.scitotenv.2022.155725
Xiao, F., Zhou, H. (2023). Plant salt response: Perception, signaling, and tolerance. Front. Plant Sci. 13, 1053699. doi: 10.3389/fpls.2022.1053699
Yang, Y., Guo, Y. (2018). Elucidating the molecular mechanisms mediating plant salt-stress responses. N Phytol. 217, 523–539. doi: 10.1111/nph.2018.217.issue-2
Yang, L., Zhao, X., Zhu, H., Paul, M., Zu, Y., Tang, Z. (2014). Exogenous trehalose largely alleviates ionic unbalance, ROS burst, and PCD occurrence induced by high salinity in Arabidopsis seedlings. Front. Plant Sci. 5, 570. doi: 10.3389/fpls.2014.00570
Yu, F., Yao, Y., Feng, J., Wang, X., Ma, J., Liu, K., et al. (2021). Enterobacter sp. FM-1 inoculation influenced heavy metal-induced oxidative stress in pakchoi (Brassica campestris L. ssp. Chinensis Makino) and water spinach (Ipomoea aquatic F.) cultivated in cadmium and lead co-contaminated soils. Plant Soil 459, 155–171. doi: 10.1007/s11104-020-04752-5
Yue, Z., Chen, Y., Chen, C., Ma, K., Tian, E., Wang, Y., et al. (2021). Endophytic Bacillus altitudinis WR10 alleviates Cu toxicity in wheat by augmenting reactive oxygen species scavenging and phenylpropanoid biosynthesis. J. Hazard Mater 405, 124272. doi: 10.1016/j.jhazmat.2020.124272
Yue, Z., Chen, Y., Wang, Y., Zheng, L., Zhang, Q., Liu, Y., et al. (2022). Halotolerant Bacillus altitudinis WR10 improves salt tolerance in wheat via a multi-level mechanism. Front. Plant Sci. 13, 941388. doi: 10.3389/fpls.2022.941388
Zemanová, V., Pavlík, M., Pavlíková, D., Tlusto, P. (2014). The significance of methionine, histidine and tryptophan in plant responses and adaptation to cadmium stress. Plant Soil Environ. 60, 426–432. doi: 10.17221/544/2014-PSE
Zhang, Y., Fernie, A. R. (2023). The role of TCA cycle enzymes in plants. Adv. Biol. 7, e2200238. doi: 10.1002/adbi.202200238
Zhang, X., Liu, C. J. (2015). Multifaceted regulations of gateway enzyme phenylalanine ammonia-lyase in the biosynthesis of phenylpropanoids. Mol. Plant 8, 17–27. doi: 10.1016/j.molp.2014.11.001
Zhang, D., Liu, J., Zhang, Y., Wang, H., Wei, S., Zhang, X., et al. (2023). Morphophysiological, proteomic and metabolomic analyses reveal cadmium tolerance mechanism in common wheat (Triticum aestivum L.). J. Hazard Mater 445, 130499. doi: 10.1016/j.jhazmat.2022.130499
Zhang, X. C., Liu, M., Zhao, Y. (2015). How does transmembrane electrochemical potential drive the rotation of Fo motor in an ATP synthase? Protein Cell 6, 784–791. doi: 10.1007/s13238-015-0217-6
Zhou, Y., Jiang, D., Ding, D., Wu, Y., Wei, J., Kong, L., et al. (2022). Ecological-health risks assessment and source apportionment of heavy metals in agricultural soils around a super-sized lead-zinc smelter with a long production history, in China. Environ. pollut. 307, 119487. doi: 10.1016/j.envpol.2022.119487
Zhu, J. K. (2000). Genetic analysis of plant salt tolerance using Arabidopsis. Plant Physiol. 124, 941–948. doi: 10.1104/pp.124.3.941
Keywords: Triticum aestivum, Cadmium, NaCl, multiple stresses, omics
Citation: Yue Z, Liu Y, Zheng L, Zhang Q, Wang Y, Hao Y, Zhang M, Chen Y, Wang Z, He L and Ma K (2024) Integrated transcriptomics, metabolomics and physiological analyses reveal differential response mechanisms of wheat to cadmium and/or salinity stress. Front. Plant Sci. 15:1378226. doi: 10.3389/fpls.2024.1378226
Received: 29 January 2024; Accepted: 13 September 2024;
Published: 01 October 2024.
Edited by:
Jayakumar Bose, Western Sydney University, AustraliaReviewed by:
Shikha Verma, Ben-Gurion University of the Negev, IsraelYingpeng Hua, Hunan Agricultural University, China
Copyright © 2024 Yue, Liu, Zheng, Zhang, Wang, Hao, Zhang, Chen, Wang, He and Ma. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Keshi Ma, emtudW1rc0AxNjMuY29t