Zhiyin Jiao1†
Yannan Shi1†Jinping Wang1Zhifang Wang1Xing Zhang1Xinyue Jia1Qi Du1Jingtian Niu1Bocheng Liu1Ruiheng Du1Guisu Ji1
Junfeng Cao2*Peng Lv1*- 1Institute of Millet Crops, Hebei Academy of Agriculture and Forestry Sciences/ Hebei Branch of National Sorghum Improvement center/ Key Laboratory of Genetic Improvement and Utilization for Featured Coarse Cereals (Co-construction by Ministry and Province), Ministry of Agriculture and Rural Affairs/ Key Laboratory of Minor Cereal Crops of Hebei Province, Shijiazhuang, China
- 2Key Laboratory of Urban Agriculture (South), Ministry of Agriculture, Plant Biotechnology Research Center, Fudan-SJTU-Nottingham Plant Biotechnology R&D Center, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai, China
Cadmium (Cd) pollution is a serious threat to plant growth and human health. Although the mechanisms controlling the Cd response have been elucidated in other species, they remain unknown in Sorghum (Sorghum bicolor (L.) Moench), an important C4 cereal crop. Here, one-week-old sorghum seedlings were exposed to different concentrations (0, 10, 20, 50, 100, and 150 μM) of CdCl2 and the effects of these different concentrations on morphological responses were evaluated. Cd stress significantly decreased the activities of the enzymes peroxidase (POD), superoxide dismutase (SOD), glutathione S-transferase (GST) and catalase (CAT), and increased malondialdehyde (MDA) levels, leading to inhibition of plant height, decreases in lateral root density and plant biomass production. Based on these results, 10 μM Cd concentration was chosen for further transcription and metabolic analyses. A total of 2683 genes and 160 metabolites were found to have significant differential abundances between the control and Cd-treated groups. Multi-omics integrative analysis revealed that the flavonoid biosynthesis pathway plays a critical role in regulating Cd stress responses in sorghum. These results provide new insights into the mechanism underlying the response of sorghum to Cd.
Introduction
Environmental pollution with heavy metals pollution, such as cadmium (Cd), lead (Pb), copper (Cu), and zinc (Zn), has rapidly increased over recent decades. This pollution not only impairs growth and development in both plants and animals, but also threatens human health, and has become a serious environmental problem (He et al., 2013; Xue et al., 2014; Mourad et al., 2021; Meng et al., 2022). Cd is highly toxic and greatly inhibits the growth and development of plant roots, stems and other organs, causing a decrease in plant yield and quality (Feng et al., 2018; Mahmud et al., 2018; Santoro et al., 2022). Remediation of cadmium-polluted soil is therefore urgent.
Hyperaccumulators are plants capable of growing in soils with high levels of contamination and able to concentrate high levels of toxins in their tissues. Cd hyperaccumulators have been used worldwide tophyto-extract the metals from the contaminated soils in an eco-friendly and green approach (Guan et al., 2019; Bai et al., 2020).
Terrestrial plants uptake Cd from the soil mainly through their roots, and during this process, the root system initially suffers from Cd stress, triggering signal transduction (Li et al., 2012; Duan et al., 2015; Abozeid et al., 2017; Meng et al., 2022). Cd toxicity in plants can lead indirectly to the formation of reactive oxygen species (ROS), which in turn leads to cellular redox imbalance, lipid peroxidation, enzyme inactivation, disruption of nutrient absorption and assimilation, and eventually leads to the death of plant cells (Gill and Tuteja, 2010; Petrov et al., 2015; Han et al., 2016). Certain plants have therefore evolved intricate Cd detoxification strategies (He et al., 2013; Yang et al., 2021). Cd can cause changes in microstructural level, subsequently altering physiological and molecular processes. In response to Cd exposure, gene expression in certain plants exhibit can change, the antioxidative system can be activated, protein abundance can be altered and certain metabolic pathways involved in Cd uptake, transport and detoxification can be activated (He et al., 2015; Corso et al., 2018; Leng et al., 2020). Such Cd response processes have been explored largely in plant species with low biomass, such as Arabidopsis thaliana (L.) Heynh. (Arabidopsis), Sedum aalfredii Hance (S.alfredii) and S olanum nigrum L(S.nigrum) (Han et al., 2016; Luo et al., 2016; Wang et al., 2022). However, a low biomass constrains the effectiveness of these plant species in phytoremediation.
Sorghum (Sorghum bicolor (L.) Moench) is a C4 model plant, and is widely grown as a food, feed and energy crop (Xin et al., 2008; Xiao et al., 2021; Shi et al., 2022). Previous studies have proposed sorghum for phytoremediation due to its high production of biomass, high tolerance to abiotic stresses and the relatively high Cd accumulation observed in some genotypes (Jia et al., 2017; Feng et al., 2018; Yuan et al., 2019; Xiao et al., 2021). Sorghum is able to remedicate Cd-contaminated soil, and the polluted plant tissues can be harvested and used in the production of ethanol gas, thereby preventing Cd pollution from re-entering the food chain. This means that sorghum could be a connection between phytoremediation and bioethanol production, achieving a win-win effect (Jia et al., 2017; Feng et al., 2018; Mishra et al., 2020). Thus, it is essential to understand Cd toxicity in this species, and its response to Cd. Hydroponics and pot experiments have been conducted to study the physiological responses of different sorghum varieties to Cd stress. Cd was found to have accumulated primarily in the roots of sorghum plants and to have been partially transferred to the above ground parts. Chlorophyll content and the activities of antioxidative enzymes were affected by Cd accumulation, and led to growth reduction (Bonfranceschi et al., 2009; Soudek et al., 2014; Jawad Hassan et al., 2020). Comparative transcription analysis experiments demonstrated that higher Cd accumulation depends on multilevel coordination of efficient Cd uptake and transport (Feng et al., 2018). However, there are only a few preliminary studies into the effects of Cd stress on the growth of sorghum, and the physiological and molecular mechanisms underlying the sorghum response to Cd stress remain elusive.
In recent years, multi-omics analysis technologies have been widely utilized to study the mechanisms underlying plant responses to abiotic stress. Transcriptomics and metabolomics represent powerful tools to comprehensively analyze the biological processes and their metabolic regulation in plants growing in stressful environments (Mishra et al., 2020; Feng et al., 2021). In this study, an integrated morpho-physiological analysis combined with transcriptomics and metabolomics was performed to gain insights into the molecular mechanisms underlying sorghum responses to Cd stress. We aimed to reveal the key responsive candidate genes and pathways that respond to Cd stress in sorghum and to provide new insights to inform the genetic engineering of new plants for use in phytoremediation.
Materials and methods
Plant material
Sorghum (Sorghum bicolor) inbred line “407B” was used in this study. Healthy seeds were immersed in distilled water and germinated in petri dishes on filter papers. Then, seedlings with uniform growth were selected and cultivated in hydroponic boxes with 1 L 1/2 Murashige and Skoog (MS) solution (pH 5.8) (Adhikari et al., 2020). All plants were cultivated in a growth room maintained at 28 ± 2°C, with a photoperiod of 16/8 (day/night). The seedlings were grown in hydroponic medium for 7 days, and then these seedlings were subjected to different concentrations of CdCl2 (0, 10, 20, 50, 100 and 150 μM) (AR, AladdinBio-Chem Technology Co., Ltd, Shanghai, China) for 3 days. After treatment, sorghum plants were collected, and the plant height, root length and lateral density were measured. All samples were ground in liquid nitrogen and stored at -80°C.
Measurement of growth index and determination of Cd concentration
The plant height, root length and the number of lateral roots were immediately measured after 3 days of growth in a Cd stress environment (Jia et al., 2017). The shoots and roots were detached and used for fresh and dry weight measurements. The shoots and roots immediately weighed to obtain the fresh weight and then dried to constant weight at 80°C to obtain the dry weight.
The Cd content in the sorghum seedlings was assayed according to a method described previously with slight modifications (Xiao et al., 2021). Briefly, roots and aboveground parts of the sorghum plants were sampled and rinsed with distilled water. Samples were then dried in an oven at 80°C until constant weight and then pulverized to a fine powder by passing through a mesh (100 mesh fractions). Afterwards, the obtained ash residues were digested using the nitric-perchloric acid (9:1, v/v) wet digestion method. The amount of Cd was measured using an atomic absorption spectrophotometer (GGX-810, Haiguang Instrument Co., Ltd, Beijing, China).
Measurement of MDA, and H2O2 levels and activities of antioxidant enzyme
Fresh roots (0.2 g) were detached from the control and Cd treatment plants. Root concentrations of the enzymes superoxide dismutase (SOD), peroxidase (POD), catalase (CAT), glutathione S-transferase (GST), and of malondialdehyde (MDA) and hydrogen peroxide (H2O2) were analyzed using a SOD assay kit (G0101W), POD assay kit (G0107W), CAT assay kit (G0105W), GST assay kit (G0208W), MDA assay kit (G0109W) and H2O2 assay kit (G0112W), respectively, following the manufacturer’s instructions (Suzhou Grace Biotechnology Co., Ltd., Suzhou, China). The assays were performed as previously described (Jia et al., 2022). Three replicates of each treatment were used in each experiment.
Total RNA extraction, preparation of cDNA library and Illumina RNA-seq
Based on the physiological results, samples from the 0 and 10 μM CdCl2 treatments were selected for transcription and metabolic analysis. The roots of the treated and control sorghum seedlings were washed with sterilized distilled water, and were immediately frozen in liquid nitrogen. Total RNA was extracted using the NEBNext®UltraTM RNA Library Prep Kit for Illumina®(NEB, USA) following manufacturer’s recommendations. The purity, concentration and integrity of RNA were assessed using a NanoPhotometer® spectrophotometer (IMPLEN, CA, USA), Qubit® RNA Assay Kit in Qubit®2.0 Flurometer (Life Technologies, CA, USA) and RNA Nano 6000 Assay Kit of the Bioanalyzer 2100 system (Agilent Technologies, CA, USA), respectively (Ren et al., 2022). One μg high quality RNA per sample was then used to construct the cDNA library, and was sequenced using an Illumina HiSeq sequencing platform by Metware Biotechnology Co., Ltd. (Wuhan, China). The software fastp v 0.19.3 was used to clean and trim the original data to obtain high-quality clean reads (reads with adapters were filtered out, as were paired reads were the N content exceeded 10% of the base number of the reads, or when the number of low-quality (Q ≤ 20) bases contained in the read exceeds 50% of the bases). Clean data were then mapped to the sorghum reference genome (Sbicolor_454_v3.0.1) using the HISAT v2.1.0 software. The expression abundance of reads was quantified using the fragments per kilobase of transcript per million base pairs (FPKM) value. The differentially expressed genes (DEGs) in the control and treatment groups were screened using DESeq2 v1.22.1 (Ross Ihaka, University of Auckland, New Zealand) with a threshold of |log2 Fold Change | ≥1 and False Discovery Rate (FDR)< 0.05.
Widely-targeted metabolomics analysis
The samples described above were also used in the metabolomics analysis. The extract analysis, and the identification and quantification of metabolites were performed by MetWare (Wuhan, China) using standard procedures and based on a previously published protocol (Sun et al., 2020; Ren et al., 2022). Briefly, the freeze-dried roots were crushed using a grinder (MM 400, Retsch) at 30Hz for 1.5 min. The powder was placed in 1.2 ml of 70% aqueous methanol solution (4°C) and vortexed six times, then left over night. The extracts were filtrated through a microporous membrane (0.22 μm pore size) after being centrifuged at 12,000 g for 3 min. The extracts were analyzed using an ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) system. Analysis of metabolite data was conducted with the Analyst 1.6.3 software. Metabolites with a fold change (FC) ≥ 2 and FC ≤ 0.5, and variable importance in projection (VIP) scores > 1 were considered to be differentially expressed metabolites (DEMs).
Statistical analysis
All the experimental data were expressed as mean ± standard error (SE). The SPSS 26.0 Program (SPSS Inc., Chicago, IL, USA) was used for one-way ANOVA analysis at a significance level of p< 0.05. Origin 8 (OriginLab, Northampton, Massachusetts, USA) was used for chart processing.
Results
Growth characteristics of plants in responding to Cd stress
To determine the potential effect of Cd stress on the growth of sorghum seedlings, the phenotypic and physiological characteristics of sorghum under Cd stress were analyzed and compared with a control group without Cd. The growth of the sorghum seedlings was markedly influenced by different concentrations of Cd, and the seedlings displayed visible poisoning symptoms including chlorosis both in leaves and roots, especially in the 150 μM Cd group (Figure 1A). With increasing Cd concentrations, the seedling height of the Cd-treated group was significantly reduced (14.92-47.73%) compared to that of the control group (Figure 1B). Seedlings treated with 10, 20, 50, 100 and 150 μM Cd had 12.85%, 53.28%, 63.87%, 86.13% and 88.32% fewer lateral roots, respectively (Figure 1C). Nevertheless, the length of the roots was not significantly different between these groups (Figure 1D). In addition, we monitored the plant fresh weight and dry weight of sorghum seedlings shoots and roots. Compared with the control, the fresh and dry weight of the aboveground and underground parts of sorghum seedlings under different cadmium concentrations were significantly reduced (Supplementary Figure S1). The seedlings treated with 10 μM Cd showed the minimum reduction in plant growth, and this concentration was chosen for further analysis. We calculated and compared the Cd concentrations in the shoot and root organs. We found that the Cd levels in the roots were more than 8.17 times higher than those in the shoots (Figures 1E, F). Taken together, these results suggest that Cd stress can negatively affect plant height and lateral root density in sorghum.
Figure 1 Effect of different Cd concentrations on (A) Phenotype, (B) Plant height, (C) The number of lateral roots and (D) Root length. (E) Cadmium concentrations in sorghum roots and (F) sorghum shoots under 10 μM Cd treatments. Asterisks denote significant differences (*P< 0.05; **P< 0.01). Bars = 1 cm.
Effect of Cd on the antioxidative enzymes in sorghum root
In order to further explore the effect of Cd treatment on sorghum, we analyzed the activities of antioxidative enzymes in the roots of sorghum seedlings. Under Cd stress, the activities of the tested antioxidative enzymes in the sorghum were all inhibited by varying degrees. The activity of CAT in sorghum seedlings following Cd treatment group was 37.57% lower than of the control group (Figure 2A). Similarly, following Cd treatment, the SOD, POD and GST activities decreased by 21.56%, 10.99% and 12.59%, respectively, compared with those in the control group (Figures 2B–D). We also determined the MDA and H2O2 content in the sorghum roots. The level of MDA in Cd-treated plants was 1.33-fold that in control plants (Figure 2E), However, the H2O2 content in the roots was not significantly different, from that in the control, although it exhibited a slight increase with increasing Cd levels (Figure 2F). In general, Cd exposure causes a severe oxidative injury in plant tissues.
Figure 2 Effect of different Cd concentrations on (A) CAT activity, (B) SOD activity, (C) POD activity, (D) GST activity, (E) MDA content and (F) H2O2 content under 0 and 10 μM Cd treatments. Data are means ± SE (n = 3). Asterisks denote significant differences (*P< 0.05; **P< 0.01).
Transcriptome analysis of sorghum roots under Cd stress
We then wanted to investigate the potential molecular mechanism underlying the phenotypes observed in sorghum seedlings under Cd stress. Cd is a kind of trace heavy metal, which can produce toxic effects on plants and human body at low doses (Haider et al., 2021). Previous studies have focused on the harm of high concentrations to plants (Leng et al., 2020; Yang et al., 2021; Wang et al., 2022), while the low dose of Cd in soil is also a common problem to be solved in human life and crop production (Kubier et al., 2019; Wang et al., 2019). Combined with our physiological data and the results of the study on cadmium tolerance in sorghum (Jia et al., 2017; Feng et al., 2018), 10μM is a suitable concentration choice, so we chose this concentration for further research Samples from the 0 and 10 μM Cd treatment groups were selected as study materials for the transcriptome and metabolome analyses. In total, 6 libraries (control and Cd treated plants, each with three biological replicates) were generated and 407,987,342 raw data were sequenced. After removing low-quality reads, a total of 397,496,580 clean reads remained. The Q20, Q30 and GC content values of the clean reads were 97.93 – 98.32%, 94.14 – 95.02% and 54.18 – 55.42%, respectively (Supplementary Table S1), indicating that the transcriptome sequencing results were reliable and could be used for further DEG analysis.
The percentage of reads that could be mapped to the sorghum reference genome was 91.86 – 93.62%, among which 89.02 – 90.76% were uniquely mapped (Supplementary Table S1), indicating that the selected reference genome was suitable. The correlation analysis revealed that treatment groups had appropriate intra - group biological repeatability (Supplementary Figure S2A). Principal component analysis (PCA) analysis demonstrated that the samples within the control and treatment groups clustered together in their respective groups, and the expression of gene clusters between the control and Cd treatment groups was clearly distinguished (Supplementary Figure S2B). Based on standard cut-offs (|fold change| > 2 and corrected p-value< 0.05), we identified the differentially expressed genes (DEGs) in sorghum roots in response to Cd stress. A total of 2683 DEGs (1812 up- and 871 down-regulated) were characterized in the control and Cd treatments (Figure 3A; Supplementary Table S2).
Figure 3 Transcriptomic analysis of sorghum roots under 0 and 10 μM Cd treatments. (A) Volcano plots of DEGs. (B) GO classification of DEGs. (C) KEGG enrichment of DEGs.
GO and KEGG pathway analysis of DEGs
To assess the biological functions of DEGs in sorghum roots under Cd treatment, Gene Ontology (GO) and KEGG enrichment analyses of DEGs were carried out. GO classification analysis showed that the DEGs could be predominantly classified into three ontologies concerning biological process (BP), molecular function (MF) and cellular component (CC), as shown in Figure 3B. Theses DEGs could be further divided into 32 functional GO terms. In the BP group, the most abundant terms were “cellular process”, “metabolic process” and “response to stimulus”. In the MF category, “binding”, “catalytic activity”, “transcription regulator activity” and “transporter activity” were the most highly represented terms; and in the CC processes, the most abundant term was “cellular anatomical entity” (Figure 3B; Supplementary Table S3). Additionally, KEGG analysis assigned these DEGs to 120 KEGG pathways (Supplementary Table S4), of which the top 50 are given in Figure 3C. In particular, “Biosynthesis of secondary metabolites” (ko01110), “Glutathione (GSH) metabolism” (ko00480), “Flavonoid biosynthesis” (ko00941), “ATP-binding cassette (ABC) transporters” (ko02010), “Phenylpropanoid biosynthesis” (ko00940), “Metabolic pathways” (ko01100) and “Flavone and flavonol biosynthesis” (ko00944) were significantly enriched. These results indicate that sorghum roots can effectively enhance GO terms such as metabolic process, response to stimulus, catalytic activity or transporter activity to improve the Cd tolerance of plants.
Metabolomic analysis of sorghum root under Cd stress
To further elucidate the responses of sorghum roots to Cd stress, we conducted widely targeted metabolomics analyses using an UPLC-MS/MS system. PCA and orthogonal projections to latent structures-discriminant analysis (OPLS-DA) score plot results showed that different samples under the same treatment conditions clustered together, whereas samples from different treatment groups were significantly separated (Supplementary Figures S3A, B). The correlation analysis of metabolites in different samples showed biological repeatability between groups (Supplementary Figure S3C). In addition, these DEMs fell into 12 categories (Supplementary Figure S3D). The most abundant classes were flavonoids (20.84%), followed by phenolic acids (17.04%), amino acids and their derivatives (12.6%), alkaloids (10.87%), lipids (10.48%), others (9.13%), organic acids (6.56%), nucleotides and their derivatives (4.69%), lignans and coumarins (3.47%), terpenoids (3.02%), quinones (1.09%) and tannins (0.19%). A total of 1555 metabolites was detected (Supplementary Table S5). We screened the DEMs based on the selection criteria listed in the methods section above. The abundance of 160 metabolites changed differentially between the two treatment groups, and included 35 down-regulated and 125 up-regulated metabolites (Figure 4A). Detailed information regarding the 160 metabolites is given in Supplementary Table S6.
Figure 4 Metabolite analysis of sorghum roots under 0 and 10 μM Cd treatments. (A) Volcano plots of DEMs. (B) KEGG classification of DEMs. (C) Heat map of DEGs in flavonoid biosynthesis pathway.
KEGG enrichment analysis of DEMs
To further confirm the key metabolic pathways related to sorghum root responses to Cd toxicity, we conducted KEGG enrichment pathways analysis. As shown in Figure 4B and Supplementary Table S7, KEGG pathway enrichment analysis indicated that these DEMs could be catagorized into 47 species, the most significant of which were metabolic pathways (38), biosynthesis of secondary metabolites (26), followed by flavonoid biosynthesis (13), ABC transporters (9) and others. The results of this analysis suggest that these metabolites could play a vital role in the response to Cd stress in sorghum roots.
We then constructed a heat map of flavonoid biosynthesis metabolite expression to better visualize the expression patterns of related metabolites in the control vs Cd-treated comparison groups. Thirteen DEMs including hesperetin-7-O-glucoside, butin; 7,3’,4’-trihydroxyflavanone, naringenin chalcone; 2’,4,4’,6’-tetrahydroxychalcone, 5,4’-dihydroxy-7-methoxyflavanone (sakuranetin), taxifolin (dihydroquercetin), phloretin, eriodictyol (5,7,3’,4’-tetrahydroxyflavanone), luteolin (5,7,3’,4’-tetrahydroxyflavone), 3,5,7-trihydroxyflavanone (pinobanksin), galangin (3,5,7-trihydroxyflavone), naringenin (5,7,4’-trihydroxyflavanone) and apigenin; 4’,5,7-trihydroxyflavone, were found to be related to flavonoid biosynthesis (Figure 4C). These results were essentially consistent with the transcription results.
Integrated analysis of the transcriptome and metabolome analyses
To understand the mechanisms underlying the response of sorghum to Cd stress, we then conducted an analysis of the relationship between DEMs and DEGs. The results of KEGG pathway enrichment analysis gave the number of differential metabolites and differential genes enriched in the flavonoid biosynthesis pathway, ABC transporters, flavone and flavonol biosynthesis pathway and other pathways (Figure 5A; Supplementary Table S8). As shown in Figure 5B and Supplementary Table S9, the co-expression network analysis of DEGs and DEMs in control and Cd-treated groups showed that the terms were mainly enriched in flavonoid biosynthesis-sorghum bicolor (sbi00941), isoflavonoid biosynthesis (map00943), and phenylpropanoid biosynthesis-sorghum bicolor (sbi00940 and sbi00260). The results suggest that the flavonoid biosynthesis pathway might have a crucial role in the responses of sorghum roots to Cd stress.
Figure 5 Transcriptomic and metabolomic analysis of sorghum roots under 0 and 10 μM Cd treatments. (A) KEGG enrichment analysis histogram of the DEGs and DEMs of sorghum roots. (B) KGML network of DEGs and DEMs. The square represents the gene or gene product (red: up-regulated, green: down-regulated), and the diamond represents the pathway. (C) Parallel displays of DEGs and DEMs involved in flavonoid biosynthesis in sorghum roots in response to Cd stress. The triangle and rectangle represent the metabolites and genes regulated under Cd treatment, respectively (red, up-regulation; blue, down-regulation). The types of metabolites are distinguished by long rectangular boxes with different colors (CHS, chalcone synthase; CHI, chalcone-flavonone isomerase; FLS, flavonol synthase; F3M, flavonoid 3’-monooxygenase; F3’5’H, flavonoid 3′,5′-hydroxylase; NOMT, naringenin 7-O-methyltransferase; FLAS, flavonol synthase II).
To further investigate the effect of Cd stress on of genes and metabolites related to flavonoid biosynthesis in sorghum roots, the DEGs and DEMs related to flavonoid biosynthesis pathway were analyzed. The key genes in the flavonoid biosynthesis pathway, including NOMT (naringenin 7-O-methyltransferase), F3’5’H (flavonoid 3′,5′-hydroxylase), CHS (chalcone synthase), CHI (chalcone-flavonone isomerase) and F3M (flavonoid 3’-monooxygenase), were found to be upregulated in the Cd-treated group. The increase in naringenin coincided with the upregulation of genes encoding CHS (chalcone synthase) and CHI (chalcone-flavonone isomerase). Moreover, metabolites associated with flavonoid biosynthesis, including hesperetin-7-O-glucoside, accumulated significantly in sorghum under Cd treatment (Figure 5C; Supplementary Table S10). The integrated analysis results demonstrated that the DEGs and DEMs related to the flavonoid biosynthesis pathway were a response to Cd stress in sorghum roots. Sorghum therefore has a complex network of gene and metabolite regulation to improve its resistance to Cd.
Discussion
Effect of Cd stress on the growth phenotype and physiology of sorghum
Cd, a toxic and nonessential heavy metal, can alter plant physiological and biochemical parameters, inhibit plant growth and hinder plant yield (Li et al., 2019; Wang et al., 2022). This study explored the effects of different concentrations of Cd on the growth of sorghum seedlings. The experimental results demonstrated the negative impacts of Cd on sorghum root growth and on the enzyme activities in the roots. Sorghum leaves curled and the tips became brownish under Cd stress, and consistent with previous studies, Cd induced chlorosis and a dwarf phenotype in these plants (Roy et al., 2016). In addition, a dose-dependent decline in shoot height, number of lateral root and plant biomass production, with high levels of Cd accumulating in both roots and shoots (Figure 1). Previous studies have demonstrated that exposure of Pb in T. pratense and maize seedlings under VI conditions has toxic effects in a concentration-dependent manner, and that Pb is more toxic at high concentrations (Adhikari et al., 2020; Meng et al., 2022). We found that the growth inhibition of sorghum seedlings treated with high Cd concentrations was greater than that under lower Cd concentrations (Figures 1B, C). A similar phenomenon has been repeatedly observed in many plant species under Cd stress (Chen et al., 2020; Khan et al., 2020). Soudek et al. (2014) found that 100-500 μM Cd treated sorghum plants after 7 days, 14 days, 28 days already have a reduction in the growth of the root system. However, our findings differ from this to some degree. Different Cd concentrations higher than 10 μM significantly reduced the number of lateral roots in sorghum seedlings in our experiments, but had no effect on root length (Figure 1D). It has been suggested that Cd affects plant root development in a time and concentration-dependent manner (Abozeid et al., 2017; Zhan et al., 2017; Cao et al., 2019), so we speculate that the Cd treatment in our experiment does not affect sorghum root length due to the short treatment time. We also observed a poisoned phenotype in the leaves and roots of young sorghum seedlings under Cd stress, especially at elevated concentrations (Figure 1A). Published studies have demonstrated that Cd has a toxic effect on plant leaves because it damages the ultrastructure of chloroplasts, so that the plant leaves appear yellowish (Liu et al., 2019; Chen et al., 2021). Furthermore, the accumulation of Cd in plants also damages the mitochondria and other organelles, significantly reducing the efficiency of photosynthesis and therefore the energy supply, further inhibiting plant growth (Minkina et al., 2018; Chen et al., 2020; Yadav et al., 2021; Wang et al., 2022). In other words, the toxic effect of Cd on plants is concentration-dependent, and Cd has a greater impact on plant growth at elevated concentrations.
Heavy metal stress leads to the production of reactive oxygen species (ROS) in plants (Adhikari et al., 2020). Plants are able to protect themselves from oxidative damage caused by stress through enzymatic systems, with the POD and CAT enzymes, for example, responsible for the conversion of H2O2 to water. Another ROS scavenging enzyme, SOD, also plays an essential role in mitigating oxidative stress. SOD can modulate the amount of superoxide radicals , which involves the dismutation of into H2O2 and O2 under stress conditions (Liu et al., 2008; Zornoza et al., 2010). The main function of GST is considered to be detoxification, and this enzyme also plays a major role in ROS scavenging (Dixon et al., 2002; Zagorchev et al., 2013). In the current work, the activities of SOD, POD, CAT and GST in the roots of sorghum all decreased following treatment of the plants with 10 μM Cd (Figure 2). This suggests that the activities of these antioxidant enzymes in sorghum roots are suppressed under 10 μM Cd stress. Previous studies have verified that Cd-induced oxidative stress of membrane lipids, harms the thylakoid membrane structure (He et al., 2017; Saleem et al., 2020). Our results show that levels of MDA in the group treated with Cd were significantly higher than those in the control group (Figure 2E). We also found that there was no significant difference in H2O2 content between the control and Cd-treated groups, although there was a slight increase. Previous studies have shown that the accumulation of H2O2 is time-and concentration-dependent (Lv et al., 2017), so we speculate that this may be the reason why H2O2 has not changed significantly in the current experiment.(Figure 2F). From our physiological data, 10µM Cd caused damage to the growth of the sorghum seedlings. Cd stress thus caused the sorghum seedlings to suffer from oxidative stress, which then impacted the growth of the seedlings.
The transcriptomic and metabolomic analyses of sorghum roots in response to Cd toxicity
To further investigate the physiological and molecular mechanisms of the response in sorghum to Cd, the metabolomic and transcriptomic analyses of sorghum roots exposed to Cd were adopted. From the transcriptomic data, we identified 2683 DEGs among the treatment groups (Figure 3A, Supplementary Table S2). These DEGs were involved in certain biological pathways, including metabolic process, response to stimulus, catalytic activity, and transporter activity (Figure 3B), indicating that the sorghum roots could efficiently activate the defense system and induce the expression of stress-related genes to when subjected to Cd toxicity. Previous studies have shown that Cd stress enhances the levels of various metabolites related to amino acid, carbohydrate, and nucleotide biosynthesis pathways via regulation of the expression of DEGs (Singh et al., 2016; Wang et al., 2022). We identified 160 DEMs (Figure 4A). In addition, we found that some metabolic pathways were induced markedly in sorghum roots under Cd stress, including flavonoid biosynthesis, ABC transporters, flavone and flavonol biosynthesis, and other pathways (Figure 4B; Supplementary Table S7). Transcriptome and cDNA-microarray studies have found that ABC transporters were involved in the tolerance and detoxification of heavy metals in plants, and have identified a large number of ABC transporter genes that are induced or inhibited under Cd treatments (Bovet et al., 2005; Li et al., 2022a). Moreover, ABC transporters in transgenic plants have higher Cd tolerance compared to the wild type, and are thought to increase their tolerance of Cd stress by modifying the transport and distribution of Cd, for example by exporting Cd from the root cells or by forming Cd conjugates and subsequently sequestrated Cd in vacuoles via ABC transporters for detoxification (Wojas et al., 2009; Park et al., 2012; Brunetti et al., 2015; Wang et al., 2017; Fu et al., 2019). Winter et al. (2015) reported that the ABC transporter can assist the transport of Cd to the xylem through the ectoplasmic pathway to transport and fix Cd. Our results confirm the function of ABC transporters in cadmium stress, and also verify the reliability of our omic data. Previous studies have also shown that phenylpropanoid biosynthesis may play an important role in Cd resistance (Feng et al., 2018; Quan et al., 2020; Zhu et al., 2020). Furthermore, KEGG pathway analysis found that the phenylpropanoid biosynthesis pathway showed significant differences between two sorghum genotypes (one high-Cd accumulation and one low-Cd accumulation) upon Cd stress (Feng et al., 2018). In our study, the phenylpropanoid biosynthesis pathway was enriched in the KEGG pathways of the DEGs, and in the co-expression network analysis of DEGs and DEMs. However, the genotype of Cd accumulation in 407B could not be determined by our study, and needs further research.
By combining the transcriptomic and metabolic results, we identified that the flavonoid biosynthesis pathway was a crucial pathway in Cd stress responses in sorghum seedlings, and that there were several DEGs and DEMs in the sorghum roots activated in the response to Cd stress. We have provided a more complete molecular outline of sorghum resistance to Cd stress (Figure 5C, Supplementary Table S10). Flavonoids are a ubiquitous group of polyphenol compounds, belonging to non-enzymatic antioxidants and composed of flavones, flavonols, flavanones, isoflavones and anthocyanins (Nijveldt et al., 2001; Sun et al., 2020). It has been demonstrated that flavonoids play an important role in tolerance to environmental stresses, such as heavy metals and drought in plants (Meng et al., 2022). The accumulation of flavonoids is a hallmark of plant stress, and the flavonoids synthesis can be activated under oxidative injury. Recent studies have shown that flavonoids not only have apparent roles in stress protection, but also participate in regulating auxin transport and altering plant growth (Li et al., 2022b; Winkel et al., 2022). In our study, certain flavonoid metabolites accumulated under Cd stress (Figure 5). CHS catalyzes the condensation of one molecule of 4- coumaroyl-CoA and three molecules of malonyl-CoA to naringenin chalcone, which is the substrate for CHI and which is converted to naringenin, CHS and CHI are the key rate-limiting enzymes in the flavonoid biosynthesis pathway (Rani et al., 2012; Liu et al., 2021). Wang et al. (2021) have shown that there were devoid of tricin-lignin in the cell walls and depleted of flavones in CHS- and CHI-deficient rice mutants. The biosynthesis of lignin affected the growth of plant seedlings and altered plant morpho-physiological under Cd treatment (Feng et al., 2018). We screened the CHS and CHI genes, which were up-regulated following Cd treatment. We observed an accumulation of the substances involved in the flavonoid biosynthesis pathway, including naringenin chalcone and naringenin (Figure 5; Supplementary Table S10). Some secondary metabolites, including flavones, are mainly used in plants for chelating metal ions to confer Cd resistance in plants (Zhu et al., 2020; Li et al., 2022; Sharma et al., 2022). In this study, flavones, including galangin, apigenin and luteolin, accumulated during Cd stress, and the genes encoding these key enzymes in this pathway, such as FLS, F3M and FLAS, were upregulated following Cd treatment (Figure 5). Increased levels of flavonoids in plants can reduce Cd poisoning through chelation and passivation (Li et al., 2015; Cai et al., 2017). This is consistent with our experimental data. Consequently, we infer that the up-regulation of genes involved in the flavonoid synthesis pathway together with the increase in flavonoid metabolites are closely related to sorghum resistance to Cd stress. These findings confirm the importance of the flavonoid metabolic pathway, which is involved in the response of sorghum roots to Cd stress.
Conclusion
In summary, this study investigated the effects of different concentrations of Cd on the growth and physiology of sorghum seedlings. Cd toxicity significantly inhibited plant height and lateral root development, especially at high Cd concentrations. Cd stress also altered the activities of antioxidant enzymes (SOD, POD, GST, CAT) and the concentrations of MDA. Moreover, transcriptome and metabolomics analyses revealed the mechanism underlying the response of sorghum to low concentrations (10 μM) of Cd. Transcriptomic and metabolomic analyses revealed that Cd treatment induced numerous pivotal DEGs and altered certain key metabolites enriched in the flavonoid biosynthesis pathway. Our findings provide new insights into the physiological and molecular mechanisms of sorghum responses to Cd stress. In this current study, we used only a single sorghum variety and a single Cd concentration. Future research should also investigate proteomics to further explore the subtle molecular mechanisms underlying sorghum responses to Cd at different concentration.
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 below: https://ngdc.cncb.ac.cn/gsa. accession number CRA009573.
Author contributions
ZJ, YS, PL, and JC conceived, supervised the experiment, revised the manuscript, and coordinated the study. ZJ, JW, ZW, XZ, XJ, and QD performed research. BL, JN, RD, and GJ assisted in writing the original manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the China Agriculture Reseach System (CARS-06-14.5-B5), Modern Agriculture Research System of Hebei Province (HBCT2018070204), HAAFS Basic Science and Technology Contract Project (HBNKY-BGZ-02) and Natural Science Foundation of Hebei Province for Youths, China (No. C2022301063).
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.2023.1144265/full#supplementary-material
Supplementary Figure 1 | Effect of different Cd concentrations on (A) Aboveground fresh weight, (B) Aboveground dry weight, (C) Underground fresh weight, (D) Underground dry weight.
Supplementary Figure 2 | Transcriptomic differences between control and Cd-treated groups. (A) PCC analysis and (B) PCA score plot.
Supplementary Figure 3 | Metabolomic differences between control and Cd-treated groups. (A) PCA score plot, (B) OPLS-DA score plot, (C) PCC analysis, (D) Classification of secondary metabolites.
Supplementary Table 1 | Statistics regarding transcriptome sequencing results.
Supplementary Table 2 | Summary information of DEGs.
Supplementary Table 3 | GO analysis of DEGs.
Supplementary Table 4 | KEGG enrichment of DEGs.
Supplementary Table 5 | Information about all metabolites.
Supplementary Table 6 | Information about DEMs.
Supplementary Table 7 | KEGG classification of DEMs.
Supplementary Table 8 | KEGG enrichment of DEGs and DEMs.
Supplementary Table 9 | Co-expression network of DEGs and DEMs.
Supplementary Table 10 | List of DEGs and DEMs involved in the flavonoid biosynthesis pathway.
References
Abozeid, A., Ying, Z., Lin, Y., Liu, J., Zhang, Z., Tang, Z. (2017). Ethylene improves root system development under cadmium stress by modulating superoxide anion concentration in Arabidopsis thaliana. Front. Plant Sci. 8, 253. doi: 10.3389/fpls.2017.00253
Adhikari, A., Adhikari, S., Ghosh, S., Azahar, I., Shaw, A. K., Roy, D., et al. (2020). Imbalance of redox homeostasis and antioxidant defense status in maize under chromium (VI) stress. Environ. Exp. Bot. 169, 103873. doi: 10.1016/j.envexpbot.2019.103873
Bai, Z. Q., Zhu, L., Chang, H. X., Wu, J. W., Luo, Z. B. (2020). Enhancement of cadmium accumulation in sweet sorghum as affected by nitrate. Plant Biol. 23, 66–73. doi: 10.1111/plb.13186
Bonfranceschi, B. A., Flocco, C. G., Donati, E. R. (2009). Study of the heavy metal phytoextraction capacity of two forage species growing in an hydroponic environment. J. Haz. Mat. 165, 366–371. doi: 10.1016/j.jhazmat.2008.10.024
Bovet, L., Feller, U., Martinoia, E. (2005). Possible involvement of plant ABC transporters in cadmium detoxification: a cDNA sub-microarray approach. Environ. Int. 31, 263–267. doi: 10.1016/j.envint.2004.10.011
Brunetti, P., Zanella, L., De Paolis, A., Di Litta, D., Cecchetti, V., Falasca, G., et al. (2015). Cadmium-inducible expression of the ABC-type transporter AtABCC3 increases phytochelatin-mediated cadmium tolerance in arabidopsis. J. Exp. Bot. 66, 3815–3829. doi: 10.1093/jxb/erv185
Cai, S.-Y., Zhang, Y., Xu, Y.-P., Qi, Z.-Y., Li, M.-Q., Ahammed, G. J., et al. (2017). HsfA1a upregulates melatonin biosynthesis to confer cadmium tolerance in tomato plants. J. Pineal Res. 62, e12387. doi: 10.1111/jpi.12387
Cao, J., Ye, C., Hao, G., Dabney-Smith, C., Hunt, A. G., Li, Q. Q. (2019). Root hair single cell type specific profiles of gene expression and alternative polyadenylation under cadmium stress. Front. Plant Sci. 10, 589. doi: 10.3389/fpls.2019.00589
Chen, P., Li, Z., Luo, D., Jia, R., Lu, H., Tang, M., et al. (2021). Comparative transcriptomic analysis reveals key genes and pathways in two different cadmium tolerance kenaf (Hibiscus cannabinus l.) cultivars. Chemosphere 263, 128211. doi: 10.1016/j.chemosphere.2020.128211
Chen, H.-C., Zhang, S.-L., Wu, K.-J., Li, R., He, X.-R., He, D.-N., et al. (2020). The effects of exogenous organic acids on the growth, photosynthesis and cellular ultrastructure of Salix variegata franch. under cd stress. Ecotoxicol. Environ. Saf. 187, 109790. doi: 10.1016/j.ecoenv.2019.109790
Corso, M., Schvartzman, M. S., Guzzo, F., Souard, F., Malkowski, E., Hanikenne, M., et al. (2018). Contrasting cadmium resistance strategies in two metallicolous populations of arabidopsis halleri. New Phytol. 218, 283–297. doi: 10.1111/nph.14948
Dixon, D. P., Lapthorn, A., Edwards, R. (2002). Plant glutathione transferases. Genome Biol. 3, 10. doi: 10.1186/gb-2002-3-3-reviews3004
Duan, X., Li, X., Ding, F., Zhao, J., Guo, A., Zhang, L., et al. (2015). Interaction of nitric oxide and reactive oxygen species and associated regulation of root growth in wheat seedlings under zinc stress. Ecotoxicol. Environ. Saf. 113, 95–102. doi: 10.1016/j.ecoenv.2014.11.030
Feng, Z., Ji, S., Ping, J., Cui, D. (2021). Recent advances in metabolomics for studying heavy metal stress in plants. TrAC Trends Anal. Chem. 143, 116402. doi: 10.1016/j.trac.2021.116402
Feng, J., Jia, W., Lv, S., Bao, H., Miao, F., Zhang, X., et al. (2018). Comparative transcriptome combined with morpho-physiological analyses revealed key factors for differential cadmium accumulation in two contrasting sweet sorghum genotypes. Plant Biotechnol. J. 16, 558–571. doi: 10.1111/pbi.12795
Fu, S., Lu, Y., Zhang, X., Yang, G., Chao, D., Wang, Z., et al. (2019). The ABC transporter ABCG36 is required for cadmium tolerance in rice. J. Exp. Bot. 70, 5909–5918. doi: 10.1093/jxb/erz335
Gill, S. S., Tuteja, N. (2010). Reactive oxygen species and antioxidant machinery in abiotic stress tolerance in crop plants. Plant Physiol. Biochem. 48, 909–930. doi: 10.1016/j.plaphy.2010.08.016
Guan, M. Y., Zhu, Y. X., Liu, X. X., Jin, C. W. (2019). Induction of s-nitrosoglutathione reductase reduces root cadmium uptake by inhibiting iron-regulated transporter 1. Plant Soil 438, 251–262. doi: 10.1007/s11104-019-04014-z
Haider, F. U., Liqun, C., Coulter, J. A., Cheema, S. A., Wu, J., Zhang, R., et al. (2021). Cadmium toxicity in plants: Impacts and remediation strategies. Ecotoxicol. Environ. Saf. 211, 111887. doi: 10.1016/j.ecoenv.2020.111887
Han, X., Yin, H., Song, X., Zhang, Y., Liu, M., Sang, J., et al. (2016). Integration of small RNAs, degradome and transcriptome sequencing in hyperaccumulatorSedum alfrediiuncovers a complex regulatory network and provides insights into cadmium phytoremediation. Plant Biotechnol. J. 14, 1470–1483. doi: 10.1111/pbi.12512
He, J., Li, H., Luo, J., Ma, C., Li, S., Qu, L., et al. (2013). A transcriptomic network underlies microstructural and physiological responses to cadmium in Populus × canescens. Plant Physiol. 162, 424–439. doi: 10.1104/pp.113.215681
He, J., Li, H., Ma, C., Zhang, Y., Polle, A., Rennenberg, H., et al. (2015). Overexpression of bacterial γ-glutamylcysteine synthetase mediates changes in cadmium influx, allocation and detoxification in poplar. New Phytol. 205, 240–254. doi: 10.1111/nph.13013
He, S., Yang, X., He, Z., Baligar, V. C. (2017). Morphological and physiological responses of plants to cadmium toxicity: A review. Pedosphere 27, 421–438. doi: 10.1016/S1002-0160(17)60339-4
Jawad Hassan, M., Ali Raza, M., Ur Rehman, S., Ansar, M., Gitari, H., Khan, I., et al. (2020). Effect of cadmium toxicity on growth, oxidative damage, antioxidant defense system and cadmium accumulation in two sorghum cultivars. Plants 9, 1575. doi: 10.3390/plants9111575
Jia, Y., Liu, H., Wang, H., Zou, D., Qu, Z., Wang, J., et al. (2022). Effects of root characteristics on panicle formation in japonica rice under low temperature water stress at the reproductive stage. Field Crops Res. 277, 108395. doi: 10.1016/j.fcr.2021.108395
Jia, W., Miao, F., Lv, S., Feng, J., Zhou, S., Zhang, X., et al. (2017). Identification for the capability of cd-tolerance, accumulation and translocation of 96 sorghum genotypes. Ecotoxicol. Environ. Saf. 145, 391–397. doi: 10.1016/j.ecoenv.2017.07.002
Khan, K. Y., Ali, B., Stoffella, P. J., Cui, X., Yang, X., Guo, Y. (2020). Study amino acid contents, plant growth variables and cell ultrastructural changes induced by cadmium stress between two contrasting cadmium accumulating cultivars of brassica rapa ssp. chinensis l. (pak choi). Ecotoxicol. Environ. Saf. 200, 110748. doi: 10.1016/j.ecoenv.2020.110748
Kubier, A., Wilkin, R. T., Pichler, T. (2019). Cadmium in soils and groundwater: A review. Appl. Geochemistry 19, 104388. doi: 10.1016/j.apgeochem.2019.104388
Leng, Y., Li, Y., Wen, Y., Zhao, H., Wang, Q., Li, S.-W. (2020). Transcriptome analysis provides molecular evidences for growth and adaptation of plant roots in cadimium-contaminated environments. Ecotoxicol. Environ. Saf. 204, 111098. doi: 10.1016/j.ecoenv.2020.111098
Li, X., Cui, X., Zhang, X., Liu, W., Cui, Z. (2019). Combined toxicity and detoxification of lead, cadmium and arsenic in Solanum nigrum l. J. Haz. Mat. 389, 121874. doi: 10.1016/j.jhazmat.2019.121874
Li, N., Liu, J., Yang, L., Kang, Y., Cao, Y., Chen, K., et al. (2022b). Synergistic effect of antioxidant systems enhance cadmium phytoextraction and translocation in Amaranthus hypochondriacus under rutin application. South Afr. J. Bot. 149, 582–590. doi: 10.1016/j.sajb.2022.06.053
Li, J., Lu, H., Liu, J., Hong, H., Yan, C. (2015). The influence of flavonoid amendment on the absorption of cadmium in Avicennia marina roots. Ecotoxicol. Environ. Saf. 120, 1–6. doi: 10.1016/j.ecoenv.2015.05.004
Li, L., Yan, X., Li, J., Wu, X., Wang, X. (2022a). Metabolome and transcriptome association analysis revealed key factors involved in melatonin mediated cadmium-stress tolerance in cotton. Front. Plant Sci. 13, 995205. doi: 10.3389/fpls.2022.995205
Li, X., Yang, Y., Zhang, J., Jia, L., Li, Q., Zhang, T., et al. (2012). Zinc induced phytotoxicity mechanism involved in root growth of Triticum aestivum l. Ecotoxicol. Environ. Saf. 86, 198–203. doi: 10.1016/j.ecoenv.2012.09.021
Liu, S., Ali, S., Yang, R., Tao, J., Ren, B. (2019). A newly discovered cd-hyperaccumulator lantana camara l. J. Haz. Mat. 371, 233–242. doi: 10.1016/j.jhazmat.2019.03.016
Liu, W., Feng, Y., Yu, S., Fan, Z., Li, X., Li, J., et al. (2021). The flavonoid biosynthesis network in plants. Int. J. Mol. Sci. 22, 12824. doi: 10.3390/ijms222312824
Liu, J.-n., Zhou, Q.-x., Sun, T., Ma, L. Q., Wang, S. (2008). Growth responses of three ornamental plants to cd and cd–Pb stress and their metal accumulation characteristics. J. Haz. Mat. 151, 261–267. doi: 10.1016/j.jhazmat.2007.08.016
Luo, Z.-B., He, J., Polle, A., Rennenberg, H. (2016). Heavy metal accumulation and signal transduction in herbaceous and woody plants: Paving the way for enhancing phytoremediation efficiency. Biotechnol. Adv. 34, 1131–1148. doi: 10.1016/j.biotechadv.2016.07.003
Lv, W., Yang, L., Xu, C., Shi, Z., Shao, J., Xian, M., et al. (2017). Cadmium disrupts the balance between hydrogen peroxide and superoxide radical by regulating endogenous hydrogen sulfide in the root tip of brassica rapa. Front. Plant Sci. 8, 985298. doi: 10.3389/fpls.2017.00232
Mahmud, J. A., Hasanuzzaman, M., Nahar, K., Bhuyan, M. H. M. B., Fujita, M. (2018). Insights into citric acid-induced cadmium tolerance and phytoremediation in Brassica juncea l.: Coordinated functions of metal chelation, antioxidant defense and glyoxalase systems. Ecotoxicol. Environ. Saf. 147, 990–1001. doi: 10.1016/j.ecoenv.2017.09.045
Meng, L., Yang, Y., Ma, Z., Jiang, J., Zhang, X., Chen, Z., et al. (2022). Integrated physiological, transcriptomic and metabolomic analysis of the response of Trifolium pratense l. to Pb toxicity. J. Haz. Mat. 436, 129128. doi: 10.1016/j.jhazmat.2022.129128
Minkina, T., Fedorenko, G., Nevidomskaya, D., Fedorenko, A., Chaplygin, V., Mandzhieva, S. (2018). Morphological and anatomical changes of Phragmites australis cav. due to the uptake and accumulation of heavy metals from polluted soils. Sci. Total Environ. 636, 392–401. doi: 10.1016/j.scitotenv.2018.04.306
Mishra, D., Kumar, S., Mishra, B. N. (2020). An overview of morpho-physiological, biochemical, and molecular responses of Sorghum towards heavy metal stress. Rev. Environ. Contaminat. Toxicol. 256, 155–177. doi: 10.1007/398_2020_61
Mourad, A. M. I., Amin, A.E.-E.A.Z., Dawood, M. F. A. (2021). Genetic variation in kernel traits under lead and tin stresses in spring wheat diverse collection. Environ. Exp. Bot. 192, 104646. doi: 10.1016/j.envexpbot.2021.104646
Nijveldt, R. J., Nood, E. v., Hoorn, D. E. C. v., Boelens, P. G., Norren, K. v., Leeuwen, P. (2001). Flavonoids: a review of probable mechanisms of action and potential applications. Am. J. Clin. Nutr. 7, 418–425. doi: 10.1093/ajcn/74.4.418
Park, J., Song, W. Y., Ko, D., Eom, Y., Hansen, T. H., Schiller, M., et al. (2012). The phytochelatin transporters AtABCC1 and AtABCC2 mediate tolerance to cadmium and mercury. Plant J. 69, 278–288. doi: 10.1111/j.1365-313X.2011.04789.x
Petrov, V., Hille, J., Mueller-Roeber, B., Gechev, T. S. (2015). ROS-mediated abiotic stress-induced programmed cell death in plants. Front. Plant Sci. 6, 69. doi: 10.3389/fpls.2015.00069
Quan, M., Liu, X., Xiao, L., Chen., P., Song., F., Lu., W., et al. (2020). Transcriptome analysis and association mapping reveal the genetic regulatory network response to cadmium stress in Populus tomentosa. J. Exp. Bot. 72, 576–591. doi: 10.1093/jxb/eraa434
Rani, A., Singh, K., Ahuja, P. S., Kumar, S. (2012). Molecular regulation of catechins biosynthesis in tea [Camellia sinensis (L.) o. kuntze]. Gene 495, 205–210. doi: 10.1016/j.gene.2011.12.029
Ren, G., Yang, P., Cui, J., Gao, Y., Yin, C., Bai, Y., et al. (2022). Multiomics analyses of two sorghum cultivars reveal the molecular mechanism of salt tolerance. Front. Plant Sci. 13, 886805. doi: 10.3389/fpls.2022.886805
Roy, S. K., Cho, S.-W., Kwon, S. J., Kama, A. H. M., Kim, S.-W., Oh, M.-W., et al. (2016). Morpho-physiological and proteome level responses to cadmium stress in sorghum. PloS One 11, e0150431. doi: 10.1371/journal.pone.0150431
Saleem, M. H., Ali, S., Rehman, M., Rana, M. S., Rizwan, M., Kamran, M., et al. (2020). Influence of phosphorus on copper phytoextraction via modulating cellular organelles in two jute (Corchorus capsularis l.) varieties grown in a copper mining soil of hubei province, China. Chemosphere 248, 126032. doi: 10.1016/j.chemosphere.2020.126032
Santoro, D. F., Sicilia, A., Testa, G., Cosentino, S. L., Lo Piero, A. R. (2022). Global leaf and root transcriptome in response to cadmium reveals tolerance mechanisms in Arundo donax l. BMC Genomics 23, 427. doi: 10.1186/s12864-022-08605-6
Sharma, P., Bakshi, P., Kaur, R., Sharma, A., Bhardwaj, R., El−Sheikh, M., et al. (2022). Inoculation of plant-growth-promoting rhizobacteria and earthworms in the rhizosphere reinstates photosynthetic attributes and secondary metabolites in Brassica juncea l. under chromium toxicity. Plant Soil 21, 1–15. doi: 10.1007/s11104-022-05765-y
Shi, Y., Wang, X., Wang, J., Niu, J., Du, R., Ji, G., et al. (2022). Systematical characterization of GRF gene family in sorghum, and their potential functions in aphid resistance. Gene 836, 146669. doi: 10.1016/j.gene.2022.146669
Singh, S., Parihar, P., Singh, R., Singh, V. P., Prasad, S. M. (2016). Heavy metal tolerance in plants: Role of transcriptomics, proteomics, metabolomics, and ionomics. Front. Plant Sci. 6, 1143. doi: 10.3389/fpls.2015.01143
Soudek, P., Petrová, Š., Vaňková, R., Song, J., Vaněk, T. (2014). Accumulation of heavy metals using sorghum sp. Chemosphere 104, 15–24. doi: 10.1016/j.chemosphere.2013.09.079
Sun, J., Qiu, C., Ding, Y., Wang, Y., Sun, L., Fan, K., et al. (2020). Fulvic acid ameliorates drought stress-induced damage in tea plants by regulating the ascorbate metabolism and flavonoids biosynthesis. BMC Genomics 21, 411. doi: 10.1186/s12864-020-06815-4
Wang, J., Chen, X., Chu, S., You, Y., Chi, Y., Wang, R., et al. (2022). Comparative cytology combined with transcriptomic and metabolomic analyses of Solanum nigrum l. in response to cd toxicity. J. Haz. Mat. 423, 127168. doi: 10.1016/j.jhazmat.2021.127168
Wang, P., Chen, H., Kopittke, P. M., Zhao, F. J. (2019). Cadmium contamination in agricultural soils of China and the impact on food safety. Environ. pollut. 249, 1038–1048. doi: 10.1016/j.envpol.2019.03.063
Wang, L., Lui, A. C. W., Liu, H., Yuri, T.-K., Chen, M.-X., Zhu, F.-Y., et al. (2021). Deficiency in flavonoid biosynthesis genes CHS, CHI, and CHIL alters rice flavonoid and lignin profiles. Plant Physiol. 188, 1993–2011. doi: 10.1093/plphys/kiab606
Wang, Y., Wang, X., Wang, C., Peng, F., Wang, R., Xiao, X., et al. (2017). Transcriptomic profiles reveal the interactions of Cd/Zn in dwarf polish wheat (Triticum polonicum l.) roots. Front. Physiol. 8, 168. doi: 10.3389/fphys.2017.00168
Winkel-Shirle, B. (2002). Biosynthesis of flavonoids and effects of stress. Curr. Opin. Plant Biol. 5, 218–223. doi: 10.1016/S1369-5266(02)00256-X
Winter, G., Todd, C. D., Trovato, M., Forlani, G., Funck, D. (2015). Physiological implications of arginine metabolism in plants. Front. Plant Sci. 6. doi: 10.3389/fpls.2015.00534
Wojas, S., Hennig, J., Plaza, S., Geisler, M., Siemianowski, O., Skłodowska, A., et al. (2009). Ectopic expression of arabidopsis ABC transporter MRP7 modifies cadmium root-to-shoot transport and accumulation. Environ. pollut. 157, 2781–2789. doi: 10.1016/j.envpol.2009.04.024
Xiao, M.-Z., Sun, R., Du, Z.-Y., Yang, W.-B., Sun, Z., Yuan, T.-Q. (2021). A sustainable agricultural strategy integrating cd-contaminated soils remediation and bioethanol production using sorghum cultivars. Ind. Crops Product. 162, 113299. doi: 10.1016/j.indcrop.2021.113299
Xin, Z., Li Wang, M., Barkley, N. A., Burow, G., Franks, C., Pederson, G., et al. (2008). Applying genotyping (TILLING) and phenotyping analyses to elucidate gene function in a chemically induced sorghum mutant population. BMC Plant Biol. 8, 103. doi: 10.1186/1471-2229-8-103
Xue, D., Jiang, H., Deng, X., Zhang, X., Wang, H., Xu, X., et al. (2014). Comparative proteomic analysis provides new insights into cadmium accumulation in rice grain under cadmium stress. J. Haz. Mat. 280, 269–278. doi: 10.1016/j.jhazmat.2014.08.010
Yadav, V., Arif, N., Kováč, J., Singh, V. P., Tripathi, D. K., Chauhan, D. K., et al. (2021). Structural modifications of plant organs and tissues by metals and metalloids in the environment: A review. Plant Physiol. Biochem. 159, 100–112. doi: 10.1016/j.plaphy.2020.11.047
Yang, Z., Yang, F., Liu, J.-L., Wu, H.-T., Yang, H., Shi, Y., et al. (2021). Heavy metal transporters: Functional mechanisms, regulation, and application in phytoremediation. Sci. Total Environ. 809, 151099. doi: 10.1016/j.scitotenv.2021.151099
Yuan, X., Xiong, T., Yao, S., Liu, C., Yin, Y., Li, H., et al. (2019). A real filed phytoremediation of multi-metals contaminated soils by selected hybrid sweet sorghum with high biomass and high accumulation ability. Chemosphere 237, 124536. doi: 10.1016/j.chemosphere.2019.124536
Zagorchev, L., Seal, C., Kranner, I., Odjakova, M. (2013). A central role for thiols in plant tolerance to abiotic stress. Int. J. Mol. Sci. 14, 7405–7432. doi: 10.3390/ijms14047405
Zhan, Y.-H., Zhang, C.-H., Zheng, Q.-X., Huang, Z.-A., Yu, C.-L. (2017). Cadmium stress inhibits the growth of primary roots by interfering auxin homeostasis in Sorghum bicolor seedlings. J. Plant Biol. 60, 593–603. doi: 10.1007/s12374-017-0024-0
Zhu, H., Ai, H., Hu, Z., Du, D., Sun, J., Chen, K., et al. (2020). Comparative transcriptome combined with metabolome analyses revealed key factors involved in nitric oxide (NO)-regulated cadmium stress adaptation in tall fescue. BioMed. Cent. 21, 601. doi: 10.1186/s12864-020-07017-8
Keywords: cadmium stress, sorghum bicolor, physiological, transcriptome, metabolome, flavonoid biosynthesis
Citation: Jiao Z, Shi Y, Wang J, Wang Z, Zhang X, Jia X, Du Q, Niu J, Liu B, Du R, Ji G, Cao J and Lv P (2023) Integration of transcriptome and metabolome analyses reveals sorghum roots responding to cadmium stress through regulation of the flavonoid biosynthesis pathway. Front. Plant Sci. 14:1144265. doi: 10.3389/fpls.2023.1144265
Received: 14 January 2023; Accepted: 10 February 2023;
Published: 23 February 2023.
Edited by:
Ting Zhao, Zhejiang University, ChinaReviewed by:
Ashutosh Sharma, DAV University, IndiaMingyang Quan, Beijing Forestry University, China
Jinjuan Tan, Zhejiang Academy of Agricultural Sciences, China
Jian Chen, Jiangsu Academy of Agricultural Sciences, China
Copyright © 2023 Jiao, Shi, Wang, Wang, Zhang, Jia, Du, Niu, Liu, Du, Ji, Cao and Lv. 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: Junfeng Cao, caojunfeng@sibs.ac.cn; Peng Lv, PENGL001@126.com
†These authors have contributed equally to this work