Zhanling Zhang1,2†
Lele Zhong3,4†
Wanting Xiao1,2
Yaping Du1,2
Guiqi Han1,2,5Zhuyun Yan1,2Dongmei He1,2*Chuan Zheng1,2 - 1School of Pharmacy, Chengdu University of Traditional Chinese Medicine, Chengdu, Sichuan, China
- 2State Key Laboratory of Characteristic Chinese Medicine Resources in Southwest China, Chengdu University of Traditional Chinese Medicine, Chengdu, Sichuan, China
- 3Evaluation and Utilization of Strategic Rare Metals and Rare Earth Resource Key Laboratory of Sichuan Province, Chengdu, Sichuan, China
- 4Chengdu Analytical & Testing Center, Sichuan Bureau of Geology & Mineral Resources, Chengdu, Sichuan, China
- 5College of Medical Technology, Chengdu University of Traditional Chinese Medicine, Chengdu, Sichuan, China
Introduction: Ligusticum chuanxiong Hort. is a widely used medicinal plant, but its growth and quality can be negatively affected by contamination with the heavy metal cadmium (Cd). Despite the importance of understanding how L. chuanxiong responds to Cd stress, but little is currently known about the underlying mechanisms.
Methods: To address this gap, we conducted physiological and transcriptomic analyses on L. chuanxiong plants treated with different concentrations of Cd2+ (0 mg·L−1, 5 mg·L−1, 10 mg·L−1, 20 mg·L−1, and 40 mg·L−1).
Results: Our findings revealed that Cd stress inhibited biomass accumulation and root development while activating the antioxidant system in L. chuanxiong. Root tissues were the primary accumulation site for Cd in this plant species, with Cd being predominantly distributed in the soluble fraction and cell wall. Transcriptomic analysis demonstrated the downregulation of differential genes involved in photosynthetic pathways under Cd stress. Conversely, the plant hormone signaling pathway and the antioxidant system exhibited positive responses to Cd regulation. Additionally, the expression of differential genes related to cell wall modification was upregulated, indicating potential enhancements in the root cell wall’s ability to sequester Cd. Several differential genes associated with metal transport proteins were also affected by Cd stress, with ATPases, MSR2, and HAM3 playing significant roles in Cd passage from the apoplast to the cell membrane. Furthermore, ABC transport proteins were found to be key players in the intravesicular compartmentalization and efflux of Cd.
Discussion: In conclusion, our study provides preliminary insights into the mechanisms underlying Cd accumulation and tolerance in L. chuanxiong, leveraging both physiological and transcriptomic approaches. The decrease in photosynthetic capacity and the regulation of plant hormone levels appear to be major factors contributing to growth inhibition in response to Cd stress. Moreover, the upregulation of differential genes involved in cell wall modification suggests a potential mechanism for enhancing root cell wall capabilities in isolating and sequestering Cd. The involvement of specific metal transport proteins further highlights their importance in Cd movement within the plant.
1 Introduction
In recent years, the pollution of soil by Cd has become a growing concern due to various sources such as metal smelting, sewage treatment, phosphorus fertilizer application, and atmospheric sedimentation (Huang et al., 2019b; Du et al., 2020; Ren et al., 2022). Cd is a toxic heavy metal, and excessive accumulation of Cd in the human body can cause damage to bones, kidneys, and lungs, resulting in long-term negative effects on human health (Staessen et al., 1999; Godt et al., 2006; Nawrot et al., 2006; Johri et al., 2010). Ligusticum chuanxiong Hort. belongs to the genus Ligusticum in the family Umbelliferae and is a famous medicinal plant. The dried rhizome of L. chuanxiong is widely recognized for its medicinal properties, especially its effectiveness in the treatment of cardiovascular and cerebrovascular diseases (Chen et al., 2018). It has been observed in several studies that traditional growing areas of L. chuanxiong have high levels of Cd in the soil (Qi-qua, 2014; Li et al., 2017; Deng et al., 2019). The excessive presence of Cd in the soil inhibits the growth of L. chuanxiong and leads to excessive accumulation of Cd in L. chuanxiong. This limits the development of the L. chuanxiong industry and threatens human health (Zeng et al., 2020; Xiao et al., 2021). Considering the detrimental effects of Cd and its impact on L. chuanxiong, it is essential to conduct research on Cd stress to uncover the mechanisms of Cd accumulation and tolerance in L. chuanxiong for future implications.
Plant growth under Cd stress is commonly associated with various detrimental effects, such as stunted growth, reduced lateral root development, decreased root vigor, and impaired water uptake and nutrient utilization (Hasan et al., 2009). Moreover, Cd accumulation in the above-ground parts of plants also leads to plant dwarfing, chlorosis, reduced biomass, and inhibition of photosynthesis (Liu et al., 2010). To reduce the toxic effects of Cd, plants have evolved a series of detoxification strategies. The cell wall is an important line of defense for Cd entry into plant cells (Tenhaken, 2014). For instance, under Cd stress, pectin in the cell wall can form chelates with Cd, while the cell wall cellulose content increases to limit the diffusion of Cd. Additionally, the lignin content of the secondary wall increases, thereby strengthening the cell wall and enhancing the plant’s tolerance to Cd (Liu et al., 2007; Ma et al., 2014). The Cd that enters the apoplast can subsequently move into the cytoplasm through cation exchange across the plasma membrane. Notably, ATPases, COPT-type, ZIP, and other metal transport proteins play crucial roles in Cd transport during this stage (Li et al., 2018b; Li et al., 2021). Some amino acids and organic acids are also involved in Cd transport, such as nicotinamide, histidine, malic acid, citric acid, etc. (Pilon et al., 2009). Substances such as reduced glutathione (GSH), phytochelatins (PCs), and metallothioneins in the cytoplasm chelate with Cd2+ and further translocate it to the vesicles for sequestration via the ABC transporter protein family. In addition, Cd stress induces a large amount of reactive oxygen species (ROS) production in plants, causing oxidative damage. At this time, the activities of antioxidative enzymes such as catalase (CAT), peroxidase (POD), and superoxide dismutase (SOD) will increase in plants to scavenge excess ROS from cells and reduce oxidative damage (Rizwan et al., 2016; Jiang et al., 2020). Overall, under Cd stress, plants employ various mechanisms, including cell wall fortification, enhanced Cd transport, chelation, and antioxidant defense activation, to alleviate the toxic effects and ensure their survival and growth.
With the rapid development of RNA sequencing (RNA-Seq) technology, RNA-Seq has been extensively utilized to investigate the molecular mechanisms of medicinal plants in response to different stresses (Li et al., 2018a; Zhao et al., 2019a; Wang et al., 2020). Nowadays, Dendrobium officinale (Jiang et al., 2020), Verbena bonariensis (Wang et al., 2019b), Mulberry (Dai et al., 2020), Trifolium repens (Wu et al., 2022), Calotropis gigantea (Yang et al., 2022), and Hibiscus cannabinus L. (Chen et al., 2021b), among other medicinal plants, are being investigated for their molecular mechanisms in response to Cd stress (Fan et al., 2023). For example, Chen et al. (2021b) found that after 20 days of Cd treatment in mulberry seedlings, the expression of RNA regulation (MYB, bHLH, and WRKY transcription factors), abiotic stress (genes encoding heat-shock protein), and signaling pathways (calmodulin-binding transcriptional activators, the protein kinase family proteins, and G-proteins) was upregulated, while the expression of secondary metabolism (genes related to the metabolism of flavonoids, phenylpropanoids, and terpenoid), hormone regulation (genes for the biosynthesis of jasmonates, ethylene and auxin), and the cell wall pathway (genes related to the synthesis, modification, and degradation of the cell wall) were downregulated. Wu et al. (2022) analyzed Dendrobium’s differential genes associated with Cd uptake, transport, and detoxification metabolic pathways. Metal transporters, sulfate glutathione metabolism, cell wall metabolism, and phenylpropanoid metabolism have important effects on adaptation to Cd stress. However, the molecular mechanisms of Cd stress response in medicinal plants are complex, and different medicinal plants may exhibit diverse Cd stress response mechanisms. At present, there are few reports on the molecular mechanisms of Cd stress response in L. chuanxiong.
In this study, we compared the integrated changes in the leaves and roots of L. chuanxiong under Cd stress by analyzing morphological characteristics and physiological indices. The transcriptome analysis of the roots of L. chuanxiong seedlings under Cd stress was performed using the RNA-Seq technique. The objectives of this study were as follows: (1) to investigate the changes in morphological characteristics and physiological responses of L. chuanxiong under Cd stress and (2) to screen key differentially expressed genes that respond to Cd treatment and preliminarily reveal the molecular mechanism of Cd accumulation and Cd tolerance in L. chuanxiong.
2 Materials and methods
2.1 Preparation of plant material and Cd treatment
Seedlings of L. chuanxiong were cultivated at the Pharmacy College, Chengdu University of Traditional Chinese Medicine, located in Sichuan Province, China (30° 41′ N, 104° 48′ E). Cluster buds of L. chuanxiong derived from plant tissue culture were planted in glass bottles with a volume of approximately 470 cm3. Each bottle contained 50 mL of 1/2 MS medium (1 mg·L−1 IBA, 8 g·L−1 agar, and 30 g·L−1 sucrose; pH 5.85). The seedlings were rooted in a culture chamber with controlled conditions, including a 14/10-h light/dark cycle, 65% humidity, and a temperature of 23°C ± 1°C. After 20 days, uniform and rooted L. chuanxiong seedlings weighing approximately 0.5 g each were selected. They were then transferred to glass bottles containing various concentrations of CdCl2·2.5H2O and 1/2 MS medium (1/2 MS, 8 g·L−1 agar, and 30 g·L−1 sucrose; pH 5.85) for the Cd stress treatment. The Cd2+ concentrations used were 5 mg·L−1 (Cd 5), 10 mg·L−1 (Cd 10), 20 mg·L−1 (Cd 20), and 40 mg·L−1 (Cd 40), while a control group (CK) without any Cd2+ was included (Cd2+ concentration of 0 mg·L−1). The culture chamber conditions remained constant throughout the experiment. Each experimental group consisted of 25 L. chuanxiong seedlings, and after 20 days of Cd stress treatment, samples were collected and processed for further analysis.
2.2 Measurement of growth indicators and sample processing
After 20 days of Cd stress treatment, L. chuanxiong plant samples from each experimental group were harvested. The samples were thoroughly rinsed with running water, followed by several rinses with distilled water. The excess surface water was removed using sterile blotting paper. The height and fresh weight of the L. chuanxiong seedlings were then measured. Subsequently, the root system was scanned using a root scanner (ScanMaker i800 Plus, Shanghai Microtek Technology Co., Ltd., Shanghai, China). The resulting images were analyzed using the Plant Root Analysis System (LA-S, Hangzhou WSEEN Testing Technology Co., Hangzhou, China) to assess parameters such as total root length, root surface area, root volume, number of root tips, and average root diameter of L. chuanxiong. Samples of L. chuanxiong roots and leaves for physiological index and subcellular fraction determination of Cd and L. chuanxiong roots for gene expression analysis were immediately frozen in liquid nitrogen and stored at −80°C in the refrigerator. All samples for total Cd content determination were dried in an oven at 60°C for 2 days, crushed, passed through a 0.3- mm mesh sieve, and then stored in a dry place for further estimations. Three biological replicates were set up for all experiments.
2.3 Determination of physiological indicators
The chlorophyll content in the leaves of L. chuanxiong was quantified following the method described by Rajalakshmi and Banu (2015). The activities of SOD, POD, CAT, and glutathione (GSH), as well as the content of malondialdehyde (MDA), were determined using commercially available kits from Suzhou Cominbio Biotechnology Co. Ltd., Suzhou, China. Each assay was performed according to the manufacturer’s instructions, and the absorbance was measured using a microplate reader (SpectraMax ABS PLUS, Molecular Devices Shanghai Co., Ltd., Shanghai, China). Activity and content calculations were carried out based on the measured absorbance values. We have shown detailed calculations and units of measure in the Supplementary Documentation (Word S1). Statistical analysis of the experimental data was conducted using Excel 2019 and SPSS 17.0, and differences between groups were assessed using the least significant difference (LSD) method with a significance level set at p < 0.05.
2.4 Determination of total Cd and subcellular distribution of Cd
Dried powders of L. chuanxiong roots and leaves (0.5 g for each sample) were digested using an electric hot plate at 170°C in a mixture of 7 mL of HNO3 (CAS 7697-37-2, 99%) and 1 mL of HClO4 (CAS 7601-90-3, 99%) (Schützendübel et al., 2001). The digestion process was carried out until complete digestion was achieved. Subsequently, the digested sample was diluted to a final volume of 25 mL with deionized water. The content of Cd in the L. chuanxiong samples was determined using ICP-MS (NexION, 2000, PerkinElmer Inc., Massachusetts, USA). The enrichment factor (EF) was equal to the ratio of the Cd content in the roots to the Cd content in 1/2 MS medium. The translocation factor (TF) was the ratio of Cd content in leaves to Cd content in roots.
The method of Weigel and Jäger (1980) was used to determine the subcellular distribution of Cd in L. chuanxiong with a slight improvement. The fresh sample of 0.5 g was ground with prechilled homogenate [250 mmol·L−1 of sucrose, 50 mmol·L−1 of Tris-HCL (pH 7.5), and 1 mmol·L−1 of dithioerythritol (C4H10O2S2)] at 4°C according to the ratio of material to liquid 1:10. The homogenates were sequentially centrifuged at 1,000 × g for 10 min and 10,000 × g for 30 min. The precipitates obtained were cell wall and organelle fractions, and the supernatants were cell-soluble components (including vacuole and ribosomes and proteins, etc.). The precipitates obtained from each component were digested using the same digestion method described above to determine the total Cd content. The supernatant was diluted with deionized water and directly used for Cd content determination.
2.5 RNA extraction, library preparation, and sequencing
Total RNA was extracted from the tissue using TRIzol® Reagent (Invitrogen Corporation, California, USA), and genomic DNA was removed using DNaseI (Takara Bio Inc., Shiga, Japan) (Pandey et al., 2022). Then the integrity and purity of the total RNA quality were determined by the 2100 Bioanalyser (Agilent Technologies, Inc., California, USA) and quantified using the ND-2000 (NanoDrop Technologies, Inc., Delaware, USA).
RNA purification, reverse transcription, library construction, and sequencing were performed at Shanghai Majorbio Bio-pharm Biotechnology Co. Ltd. (Shanghai, China). The transcriptome library was prepared following the TruSeqTM RNA sample preparation Kit from Illumina (San Diego, CA, USA) using 1 μg of total RNA. Shortly, messenger RNA was isolated according to the polyA selection method by oligo (dT) beads and then fragmented by fragmentation buffer first. Secondly, double-stranded cDNA was synthesized using a SuperScript double-stranded cDNA synthesis kit (Invitrogen, CA) with random Hexamer primers (Illumina). The synthesized cDNA was then subjected to end-repair, phosphorylation, and “A” base addition according to Illumina’s library construction protocol. Libraries were size selected for cDNA target fragments of 300 bp on 2% Low Range Ultra Agarose, followed by PCR amplification using Phusion DNA polymerase (NEB) for 15 PCR cycles. After quantified by TBS380, the paired-end RNA-seq library was sequenced with the Illumina NovaSeq 6000 sequencer (2 bp × 150 bp read length).
2.6 De novo assembly and annotation
The raw paired-end reads were trimmed, and quality was controlled by FastP with default parameters. Subsequently, clean data from the samples were used to do de novo assembly with Trinity. The assembled transcripts were then assessed and optimized with BUSCO, TransRate, and CD-HIT. All the assembled transcripts were searched against the NCBI protein nonredundant (NR), Swiss-Prot, Pfam, Clusters of Orthologous Groups of proteins (COG), GO, and KEGG databases using BLASTX to identify the proteins that had the highest sequence similarity with the given transcripts to retrieve their function annotations, and a typical cutoff E-values of less than 1.0 × 10−5 was set.
2.7 Differential expression analysis and functional enrichment
To identify differentially expressed genes (DEGs) between two different groups, the expression level of each gene was calculated according to the transcripts per million reads (TPM) method. RSEM was used to quantify gene abundances (Smyth, 2005; Li and Dewey, 2011). Essentially, differential expression analysis was performed using the DEGseq2, DEGs with |log2 (foldchange)|≥1, and p-adjust ≤ 0.05 were considered to be significantly differentially expressed genes. In addition, functional-enrichment analysis including GO and KEGG was performed to identify which DEGs were significantly enriched in GO terms and metabolic pathways at p-adjust ≤ 0.05 compared with the whole-transcriptome background. GO functional enrichment and KEGG pathway analysis were carried out by Goatools (https://github.com/tanghaibao/Goatools) and KOBAS (http://kobas.cbi.pku.edu.cn/home.do) (Xie et al., 2011).
2.8 qRT-PCR validation
Roots of L. chuanxiong exposed to 20 mg·L−1 and 40 mg·L−1 Cd for 20 days were used in qRT-PCR. A total of 10 random DEGs were selected for qRT-PCR analysis, with each gene being analyzed with three replicates. The RNA used for qRT-PCR analysis was isolated from the same RNA sequencing samples mentioned earlier. qRT-PCR was performed on a Fluorescent Quantitative PCR Detection system (LineGene 4800, BIOER, Hangzhou, China) using SYBR-GREEN1 fluorescent reagents (Foregene, Chengdu, China). Thermocycler settings were as follows: 95°C for 3 min, 40 cycles of 95°C for 15 s, 60°C for 30 s, and 72°C for 10 s. In addition, melting curves from 65°C to 95°C were used to confirm reaction specificity. All primers are listed in Supplementary Table S1. Actin served as a normalization control. The 2−ΔΔCt method was used to calculate relative gene expression.
3 Results
3.1 Changes in morphological characteristics of L. chuanxiong under Cd treatment
Compared to the control group (CK), the L. chuanxiong plants treated with 5 mg·L−1 and 10 mg·L−1 showed no apparent signs of toxicity. However, the plants treated with 20 mg·L−1 and 40 mg·L−1 exhibited clear symptoms of toxicity, including yellowing and wilting (Figure 1A). Except for Cd 5, the height, total root length, and fresh weight of the aboveground and belowground parts of the L. chuanxiong plant were inhibited to different degrees with increasing Cd concentration (Figures 1B, C). Compared to the CK group, the height of plants exposed to 40 mg·L−1 Cd treatment was significantly reduced from 55.87 mm to 38.53 mm, total root length was significantly reduced from 100.52 mm to 19.98 mm, and aboveground mass and belowground biomass were significantly reduced by 61.15% and 32.54%, respectively. The root scan results showed that the trends of root surface area, root volume, and root tip number were consistent with the above results except for the mean root diameter (Figures 1D–G). Cd stress significantly affected the growth and development of L. chuanxiong, and the high concentration of Cd treatment inhibited the growth and reduced the biomass accumulation of L. chuanxiong.
Figure 1 Effect of different concentrations of Cd treatment on the growth of L. chuanxiong. (A) Growth of L. chuanxiong under different concentrations of Cd treatment; (B) changes in biomass of L. chuanxiong; (C) plant height and total root length; (D) root surface area; (E–G) root volume, root tip number and mean root diameter, in that order. Bars represent the mean ± standard deviation (n = 3). Values marked with different letters differ significantly (p < 0.05).
3.2 Effect of Cd treatment on the physiological indices of L. chuanxiong
The chlorophyll content (chlorophyll a, chlorophyll b, and total chlorophyll) in the leaves of L. chuanxiong decreased gradually as the Cd concentration increased, and it was significantly lower than the control at the 40-mg·L−1 treatment (Figure 2A). In the plant’s antioxidant system, the activities of antioxidant enzymes CAT, SOD, and POD and nonenzymatic antioxidant GSH, as well as the MDA content in leaves, were positively regulated by Cd stress overall and reached a maximum at 40 mg·L−1 Cd treatment. The trends of CAT, SOD, POD, and GSH in roots were basically the same as those in leaves, and the MDA content reached a maximum at 10 mg·L−1 and then gradually decreased (Figures 2B–F). The high concentration of Cd treatment reduced chlorophyll accumulation in the leaves of L. chuanxiong, resulting in a decrease in the photosynthetic capacity of L. chuanxiong seedlings. Cd stress activated the antioxidant system of L. chuanxiong, and the significant increase in the activity of antioxidant substances was beneficial to alleviating the toxic effects of Cd.
Figure 2 Changes in physiological indicators of chuanxiong after 20 days of treatment with different Cd concentrations. (A) Chlorophyll content; (B) CAT activity in leaves and roots; (C–F) indicate SOD, POD, GSH, and MDA, respectively. Bars represent the mean ± standard deviation (n = 3). Values marked with different letters are significantly different (p < 0.05). FW indicates fresh weight.
3.3 Total and subcellular distribution of Cd
The Cd content in the CK group was not measured due to being below the detection limit of the device, while the remaining groups showed a positive correlation with the Cd treatment concentration, both in roots and leaves. The Cd content in leaves and roots reached a maximum value of 843.22 mg·kg−1 and 2,389.83 mg·kg−1 at 40 mg·L−1 concentration, which was 757.98% and 730.06% higher compared to 5 mg·L−1, respectively (Figure 3A). The EF of L. chuanxiong gradually decreased with increasing Cd concentration and increased to 59.74 after reaching the lowest point of 39.62 at 20 mg·L−1 concentration treatment. The TF, on the other hand, showed no obvious trend and was stable between 0.32 and 0.36 (Table 1). The Cd content in the roots was significantly higher than that in the aboveground parts in all the different treatment groups, which showed that in L. chuanxiong, the Cd accumulated mainly in the roots. The enrichment and transfer coefficients of L. chuanxiong demonstrated its significant ability to enrich Cd in the growing environment and transfer it from the underground part to the aboveground part in a stable proportion.
Figure 3 Cd content and subcellular distribution of Cd in roots and leaves after 20 days of treatment with different Cd concentrations. (A) Cd content in leaves and roots; (B) subcellular distribution of cadmium in roots; (C) subcellular distribution of cadmium in leaves; (D, E) percentage of Cd content of each fraction in roots and leaves. Bars represent mean ± standard deviation (n = 3). Values marked with different letters are significantly different (p < 0.05). DW, dry weight; FW, fresh weight.
Table 1 Enrichment factor and translocation factor of L. chuanxiong under different Cd concentration treatments.
The results of the subcellular distribution of Cd in the roots and leaves of L. chuanxiong showed that the higher the concentration of Cd treatment, the greater the accumulation of Cd in the cell wall, organelle, and soluble fraction (Figures 3B, C). The shares of soluble fraction, cell wall fractions, and organelles in the roots were 60%–80%, 18%–33%, and 2%–5%, respectively. The shares in leaves were 62%–83%, 13%–31%, and 4%–7% (Figures 3D, E). Cd was mainly distributed in the cell-soluble fraction, followed by the cell wall, and the lowest percentage of Cd was found in the cell organelles. The first two contained up to 93%–98% of Cd in the whole cell, which was the main storage location of Cd in L. chuanxiong cells.
3.4 RNA-seq and annotation results
Combined with the changes in phenological characteristics and physiological analysis of L. chuanxiong under Cd treatment, the growth of L. chuanxiong was significantly affected by Cd treatment under Cd 20 and Cd 40 treatments. Therefore, we selected three groups, CK, Cd 20, and Cd 40, of Cd-treated roots of L. chuanxiong for 20 days for transcriptome sequencing analysis. A total of nine cDNA libraries were constructed and sequenced from three independent biological replicates. More than 41.36 million (6.04 Gb) clean reads were obtained for each sample (Supplementary Table S2). The quality fraction of Q30 levels ranged from 92.42% to 93.61%, and the average GC content ranged from 42.56% to 43.08% (Supplementary Table S2). The average length was 980 BP (Supplementary Table S2). In comparison to the six publicly available databases, 90,335 annotated unigenes were obtained. The number of successfully annotated unigenes in the six databases is shown in Figure 4A. Daucus carota L. provided 79.22% similarity to L. chuanxiong, followed by Apium graveolens L. (2.08%) (Figure 4B). This indicates high sequencing quality and sufficient sequencing depth for transcriptome analysis. The raw data obtained are stored in the National Center for Biotechnology Information Biology Project (NCBI) database (access number PRJNA960886, access link: https://www.ncbi.nlm.nih.gov/sra/PRJNA960886).
Figure 4 RNA-seq annotation results. (A) Number of unigenes successfully annotated to different databases; (B) species similarity.
3.5 Analysis of differentially expressed genes in GO term and KEGG metabolic pathways
To further understand the altered gene expression levels of L. chuanxiong in response to Cd stress, differential expression analysis was performed with DEGseq2. Compared with the control, 4,050 (2,125 upregulated, 1,925 downregulated) and 6,174 (2,861 upregulated, 3,313 downregulated) DEGs were obtained for Cd 20 and Cd 40, respectively. Both together had 1,173 upregulated and 1,332 downregulated differential genes (Figures 5A, B).
Figure 5 Annotation results of differentially expressed genes. (A, B) Venn analysis of differentially expressed genes: (A) upregulated; (B) downregulated. (C) Go term annotation results. (D) KEGG pathways annotation results.
Cd 20 and Cd 40 are annotated to 48 and 49 GO terms, respectively. Among them, 21, 15, 12, and 21 GO terms belong to biological processes (BP), 15 to molecular functions (MF), and 13 to cellular components (CC) (Supplementary Table S3). Among the GO terms shared by both, catalytic activity and binding represent the most abundant categories of MF; cell part and membrane part appear more in the category of CC; and cellular processes and metabolic processes account for the largest proportion of BP (Figure 5C). KEGG analysis revealed 118 (1,094 DEGs) and 127 (1,929 DEGs) pathways, respectively (Supplementary Table S4). Compared with the control group, the common pathways of Cd 20 and Cd 40 are mainly concentrated in energy metabolism, lipid metabolism, carbohydrate metabolism, signal transduction, and amino acid metabolism (Figure 5D). Among them, differential expression of metabolic pathways such as biosynthesis of unsaturated fatty acids, MAPK signaling pathway —plant, and glutathione metabolism genes were mostly upregulated. Photosynthesis, plant hormone signal transduction, photosynthesis-antenna proteins, phenylpropanoid biosynthesis, and carbon fixation in photosynthetic organisms were mostly downregulated (Supplementary Figure S5).
3.6 DEGs and metabolic pathways associated with Cd transport and detoxification
All 42 DEGs of the photosynthetic pathway were significantly downregulated (Supplementary Table S6), mainly associated with photosystem I (20), photosystem II (9), and ferredoxin (5). There were 19 DEGs of the photosynthetic antenna protein, all significantly downregulated, of which 14 DEGs were associated with chlorophyll a– b-binding proteins (Supplementary Table S6). The carbon fixation in photosynthetic organisms (four upregulated, 14 downregulated) was dominated by DEGs such as glyceraldehyde-3-phosphate dehydrogenase (two upregulated, three downregulated), fructose-bisphosphate aldolase (two downregulated), fructose-1,6-bisphosphatase (two downregulated), phosphoenolpyruvate carboxylase (two downregulated), and ribulose diphosphate carboxylase (two downregulated) (Supplementary Table S6).
There were 35 DEGs in the antioxidant system of L. chuanxiong that were affected by Cd stress (14 upregulated, 21 downregulated), of which 28 DEGs were associated with POD (11 upregulated, 17 downregulated). In the glutathione metabolic pathway, there were 18 genes whose expression was regulated by Cd (13 upregulated and five downregulated). These included glutathione S-transferase (four upregulated, two downregulated), glutathione transferase (three upregulated), glutathione peroxidase (two upregulated), glutathione (one upregulated), glutathione reductase (one upregulated), etc. (Supplementary Table S7).
Auxin is an important plant growth regulator that can induce root growth by enhancing cell division and elongation. In the plant hormone signal transduction pathway, 15 DEGs associated with the auxin synthesis pathway were significantly downregulated in most DEGs, except for SAUR family genes. Ethylene is a simple and very important gaseous phytohormone, and six ethylene signaling pathway-related genes (EBF1, ERF1, ETR) were upregulated in expression under Cd stress. In addition, abscisic acid signaling pathway (PYL, PP2C) genes, salicylic acid transcription factor (TGA2), and gibberellin receptor (GID1C) were also regulated by Cd (Supplementary Table S8).
In the plant cell wall modification pathway, the DEGs of pectin esterase (seven upregulated, four downregulated) and cellulose synthase (one upregulated) are regulated by Cd. The expression levels of DEGs related to the lignin biosynthesis pathway, such as laccase (3), caffeic acid O-methyltransferase, bergaptol O-methyltransferase, and cinnamaldehyde dehydrogenase, were upregulated (Supplementary Table S9).
Among the metal transporter family, we annotated 11 ABC transporter proteins (nine upregulated and two downregulated), four calcium transporter proteins (ATPase, three upregulated), magnesium transporter protein (MRS2), and Cd/Zn transporter ATPase (HAM3)-related genes differentially expressed under Cd stress (Supplementary Table S10).
3.7 qRT-PCR validation
To validate the RNA-seq results, 10 DEGs were randomly selected for qRT-PCR analysis to determine the mRNA levels after exposure to 20 mg·L−1 and 40 mg·L−1 Cd2+. The results showed that the relative gene expression trends of qRT-PCR were consistent with RNA-seq, demonstrating the reliability of the RNA-seq data (Figure 6).
Figure 6 Validation of RNA-Seq results using qRT-PCR. Comparing the log2FC values of RNA-seq and qRT-PCR (log2FC > 0 indicates upregulation, log2FC < 0 indicates downregulation), if the expression trends of RNA-seq and qRT-PCR are consistent, it indicates that the results of RNA-seq are available (n = 3).
4 Discussion
4.1 Effect of Cd treatment on the growth index of L. chuanxiong
We evaluated the effect of Cd treatment at different concentrations on the growth of L. chuanxiong. Low concentrations of Cd treatment stimulated the growth of L. chuanxiong, and similar findings were found for plants such as Oxalis Calotropis gigantea L. (Yang et al., 2022) and Lonicera japonica Thunb. (Jia et al., 2015). As the concentration of Cd increased, so did the growth of L. chuanxiong seedlings, which was inhibited with a decrease in plant height and biomass. The 40-mg·L−1 Cd treatment resulted in a significant reduction in total root length, root surface area, root volume, and root tip number in L. chuanxiong, and similar results were found in studies on Cd treatment of Solanum nigrum L. (Guo et al., 2009). Previous studies have demonstrated that plants exhibit the highest efficiency in cation absorption near the root tips (Piñeros et al., 1998; Huang et al., 2019a; Rivetta et al., 2023). The reduction in the number of root tips and root surface area significantly diminishes the effective area of the plant root system for Cd uptake, ultimately determining its uptake capacity (Hasan et al., 2009). The inhibition of root system growth is a protective mechanism that prevents L. chuanxiong from excessive contact with Cd, thus minimizing its uptake and accumulation. This approach is advantageous in mitigating Cd toxicity and enhancing L. chuanxiong’s tolerance toward Cd. However, the inhibition of root elongation also hampers L. chuanxiong’s ability to absorb and utilize water and nutrients, consequently limiting its overall growth (Hasan et al., 2009).
4.2 Cd accumulation and subcellular distribution of Cd in L. chuanxiong
Except for a few hyperaccumulating plants, most plants store Cd mainly in the roots (Zacchini et al., 2009). As with many plants, more than 63% of Cd was concentrated mainly in the roots of L. chuanxiong. The Cd content of all parts of L. chuanxiong showed a linear increase with increasing Cd treatment concentration. The Cd enrichment factor in its roots fluctuated between 39.62 and 59.74, with the translocation factor less affected by Cd concentration. The roots of L. chuanxiong have a strong Cd enrichment capacity and are highly susceptible to enrichment from the survival environment to higher concentrations of Cd, which are transferred to the leaves in a certain proportion. The high accumulation characteristics of Cd in L. chuanxiong, coupled with the increasing Cd content in the soil of traditional cultivation areas, have contributed to the increasing Cd pollution problem of L. chuanxiong in recent years (Ma et al., 2014; Qi-qua, 2014; Li et al., 2017). At present, implementing soil ecological remediation in traditional planting areas or selecting alternative planting sites with lower Cd content may be an effective measure to solve the problem of Cd contamination in L. chuanxiong.
The subcellular distribution of heavy metals is closely related to the tolerance and detoxification of heavy metals in plants. Our research findings reveal that approximately 18% –33% of Cd is stored in the cell wall. The cell wall contains many hydroxyl, thiol, and carboxyl functional groups, which can form Cd complexes with Cd and prevent the flow of Cd to the protoplasts (Harholt et al., 2010). The study by Feng et al. (2018) found that when plants are subjected to Cd stress, their cell walls thicken, thereby binding more Cd, which is consistent with our study. More than 62% of Cd is stored in the cell- soluble fraction. Vacuoles are the main soluble fraction of cells. When the binding capacity of Cd on the cell wall reached saturation, a large amount of Cd passing through the cell wall was transported to the vacuole for sequestration (Mendoza-Cózatl et al., 2010). These findings indicate that in L. chuanxiong, Cd mainly accumulates in the roots through vacuole compartmentalization and fixation within the cell wall. This strategy reduces the transfer of Cd to the aboveground parts of the plant, ensuring their growth and development, thereby serving as an important tolerance mechanism for L. chuanxiong.
4.3 Effect of Cd treatment on photosynthesis of L. chuanxiong seedlings
During photosynthesis in plants, chlorophyll plays a crucial role in the absorption and utilization of light energy. A high concentration of Cd significantly inhibited chlorophyll accumulation in the L. chuanxiong leaves. The reduction of chlorophyll content will affect the plant’s photosynthetic efficiency and lead to the failure of the plant to grow normally. Transcriptome sequencing results showed that most of the DEGs in photosynthesis, photosynthesis-antenna proteins, and carbon fixation in photosynthetic organism pathways were significantly downregulated under Cd stress treatment. Cd stress inhibited light energy uptake, electron transfer, and the Calvin cycle from disrupting photosynthesis in L. chuanxiong, resulting in a decrease in photosynthetic capacity and growth inhibition in L. chuanxiong. Similar findings were observed in a study on cotton (Farooq et al., 2016).
4.4 Cd treatment activates the antioxidant system of L. chuanxiong
Cd stress leads to the production of a large amount of reactive oxygen species (ROS) in plants and causes oxidative damage (Nagajyoti et al., 2010; Ehsan et al., 2014). Plants possess antioxidative enzymes (POD, SOD, CAT) that can convert toxic H2O2 and into harmless oxygen and water, effectively scavenging ROS (Haider et al., 2021). In our study, we observed that the activities of antioxidative enzymes (CAT, POD, SOD) in L. chuanxiong increased gradually with higher concentrations of Cd treatment. Additionally, we identified several differentially expressed genes related to the antioxidant system in L. chuanxiong that responded positively to Cd stress. This suggests that the antioxidant system of L. chuanxiong plays an important role in scavenging ROS and resisting oxidative damage caused by Cd stress. MDA is one of the breakdown products of polyunsaturated fatty acids in membranes and is commonly used to evaluate the degree of oxidative damage (Demiral and Türkan, 2005). In our study, we observed a gradual increase in MDA content in the leaves of L. chuanxiong, while it remained relatively stable in the roots. This suggests that the root system has a stronger capacity to resist oxidative damage compared to the leaves. GSH, as a nonenzymatic antioxidant, also plays an important role in the plant antioxidant system. On the one hand, it can directly reduce ROS and alleviate Cd-induced oxidative stress. On the other hand, Cd ions can also form complexes with GSH and then translocate to the vacuole or apoplast for detoxification (Yadav, 2010). The subcellular distribution of Cd in L. chuanxiong showed that Cd is mainly distributed in the soluble components of cells, which may be related to the chelation and transport of GSH. We annotated in our RNA-seq results that several DEGs (GST and GPX, etc.) of the glutathione pathway are upregulated. Some of the GST genes have been shown to have enhanced Cd tolerance and reduced Cd accumulation in plants (Zhao et al., 2009; Dixit et al., 2011; Wu et al., 2019).
4.5 Plant hormone signal transduction pathway positively responds to Cd treatment and may affect the growth and Cd accumulation of L. chuanxiong
Plant hormone signaling pathways (IAA, ethylene, SA, GA, JA, and ABA) were found to respond positively to Cd regulation in the present study. Under Cd stress, changes in plant hormone levels affect plant growth and development and mediate the accumulation of Cd in plants (Qin and Huang, 2018; Nguyen et al., 2021). Cd stress resulted in significant down regulation of gene expression levels of the GH3, AUX/IAA, and AFR families in the IAA biosynthetic pathway of L. chuanxiong. Cd severely disrupts the synthesis and transport of IAA in plants and inhibits the growth of their root meristematic tissues (Yuan and Huang, 2016; Fattorini et al., 2017; Zhang et al., 2018). Exogenous addition of IAA can alleviate Cd toxicity in plants by reducing Cd uptake and enhancing antioxidant activity (Agami and Mohamed, 2013; Singh and Prasad, 2015). The expression of five ethylene receptors (ETR) and the ethylene-responsive transcription factor (ERF1) was up regulated in L. chuanxiong under Cd stress, indicating the ethylene signaling pathway was active under Cd stress. High concentrations of Cd caused rapid ethylene synthesis in plants, inhibiting cell proliferation in the meristematic tissue zone and cell elongation in the elongation zone of plant roots to hinder root growth (Růzicka et al., 2007; Schellingen et al., 2015; Qin and Huang, 2018; Wakeel et al., 2018). Salicylic acid (SA) is a ubiquitous plant phenolic compound, and SA pretreatment can enhance the plant antioxidant system, improve photosynthetic capacity, and regulate the uptake and distribution of Cd in plants, thus reducing the inhibition of plant root growth caused by Cd toxicity (Guo et al., 2007; Liu et al., 2016; Emamverdian et al., 2020; Zheng et al., 2022). TGA is a key factor in the SA signaling pathway and determines the downstream responses. In the transcriptome sequencing results of the L. chuanxiong, five TGA family SA transcription factors were significantly down regulated, which may be the mechanism of the involvement of the SA signaling pathway in the regulation of growth and Cd accumulation in L. chuanxiong under Cd stress. Expression of the gibberellin (GA) receptor GID1 was up regulated under Cd stress. GA was able to affect the division and expansion of plant cells under Cd stress, increase the growth of radicle cells in meristematic tissues, and increase the germination rate and shoot and root length of plants affected by stress (Emamverdian et al., 2020). For example, GA application improved root morphology and vigor of lettuce under Cd treatment and attenuated the toxic effects of Cd on lettuce leaves by reducing root uptake of Cd and root crown translocation of Cd (Chen et al., 2021a). In addition, DEGs of the GA signaling pathway (GH3 and JAR1) and abscisic acid (ABA) signaling pathway (PYL and PP2C) were also regulated by Cd. Exogenous application of jasmonic acid (JA) helped enhance rice’s antioxidant capacity, promoted the chelation of Cd with soluble pectin in rice roots, increased cell wall compartmentalization of Cd, and inhibited the transfer of Cd to protoplasts, thus reducing Cd uptake in roots and aboveground parts (Singh and Shah, 2014; Li et al., 2022). Abscisic acid (ABA) has also been shown to reduce Cd uptake in strawberries effectively (Kocaman, 2023). In conclusion, various phytohormone signaling pathways play important roles in plant responses to Cd stress. These pathways regulate plant growth, Cd accumulation, and antioxidant activity, and their modulation (e.g., application of exogenous hormones) may alleviate the Cd toxicity of L. chuanxiong.
4.6 Cd treatment may enhance the cell wall structure of the root system of L. chuanxiong, in order to improve the plant’s tolerance to Cd
The cell wall is composed of cellulose, hemicellulose, pectin, and lignin (Cosgrove, 2005), and it plays a crucial role in plant defense against Cd toxicity. By limiting the transport of Cd across the cell membrane, the cell wall can reduce the Cd concentration in the protoplasts and protect plants from toxicity (Krzesłowska, 2011; Le Gall et al., 2015). The plant cell wall modification pathway showed upregulation of seven DEGs associated with pectin esterase synthesis. The increased activity of pectin esterase (PME) leads to higher levels of de-esterified pectin, resulting in increased negative charges available to bind Cd2+ ions (Xiong et al., 2009). For instance, a study on alfalfa seedlings demonstrated that aeration of the root system increased pectin content and PME activity, ultimately reducing Cd toxicity (Zhao et al., 2019b). In our study with L. chuanxiong seedlings, we observed upregulation of the cellulose synthase (CESA) gene involved in cellulose synthesis after Cd stress exposure. Similar results were found by Guo et al. (2016) in their investigation of Miscanthus sacchariflorus. CESA is a key enzyme involved in plant cell wall biosynthesis. Previous studies have shown that interference with CESA, such as in the rice CESA9 mutant created by Song et al. (2013), can reduce Cd accumulation in plants. Furthermore, lignin biosynthesis plays a role in Cd accumulation in plants (Finger-Teixeira et al., 2010). In our transcriptome sequencing results of L. chuanxiong, we observed upregulation of genes involved in lignin biosynthesis, including laccase, caffeic acid-O-methyltransferase, bergaptol O-methyltransferase, and cinnamaldehyde dehydrogenase cinnamyl alcohol dehydrogenase, in the roots of L. chuanxiong. Under Cd stress, lignin content in plant cell walls increased and cell wall structure was enhanced, resulting in reduced Cd penetration and enhanced plant tolerance to Cd (Rui et al., 2016; Zhou et al., 2016). Increased lignin content also enhances the lignification of the endodermis of plant roots, which can effectively prevent a higher proportion of Cd from entering the xylem through the apoplast, thus limiting Cd transport to the aboveground parts (Feng et al., 2018). To summarize, the cell wall components, including pectin esterase, cellulose synthase, and lignin, play crucial roles in plant defense against Cd toxicity. Enhancing their activities and expression levels can improve Cd tolerance. Further research in this area is warranted to deepen our understanding of the mechanisms underlying plant responses to Cd stress.
4.7 Metal transporter proteins play an important role in Cd transport and accumulation in L. chuanxiong
Heavy metal transporters are essential for the transport and accumulation of heavy metals in plants (Fan et al., 2023). Research has found that some metal transporters in the ZIP family, NRAMP family, and ABC transporter family in plants have the function of transporting Cd (Sterckeman and Thomine, 2020; Fan et al., 2023). For example, the NRAMP family can bind and transport Mn2+, Fe2+, and Zn2+ and participate in the absorption of Cd2+ by plant cells and the transport of Cd/Cd complexes (Oomen et al., 2009; Sasaki et al., 2012; Ehrnstorfer et al., 2014). The IRT1 of the ZIP family can transport various metal ions such as Cd2+, Zn2+, Mn2+, and Fe2+ (Guerinot, 2000; Long et al., 2018; Halimaa et al., 2019). Some Ca2+ transporter proteins in the cell membrane have been shown to have the ability to transport Cd2+ (Kabała and Kłobus, 2005; Demidchik et al., 2018; Rivetta et al., 2023). In the RNA-seq results, we annotated the upregulated expression of several Ca transporter protein-related genes (ATPase4 and ATPase13), and they may be involved in the endocytic flow of Cd in L. chuanxiong cells. The ABC superfamily (ABC), one of the largest and most diverse protein families in plants, actively participates in the active transport of metal chelates between biological membranes (Fan et al., 2023). The ABCB and ABCC protein subfamilies of the ABC transporter protein family, which are usually localized to the tonoplast, transport glutathione-Cd and phytochelatin-Cd complexes from the cytoplasm to the vacuole for sequestration, respectively, to enhance the tolerance of plants to Cd (Ortiz et al., 1992; Ortiz et al., 1995; Fu et al., 2017). In our results, several genes of ABCB1 and ABCC2 were significantly upregulated, and they played an important role in intravesicular transport and sequestration of Cd in L. chuanxiong, as well as a key factor in enhancing Cd tolerance and the accumulation capacity of L. chuanxiong. In addition to the ABC transporter family, a gene related to the Cd/Zn transporter (HMA3) has also been upregulated. HMA3 is located on the tonoplast and has been shown to have the function of transporting Cd, which may also be involved in the transfer of Cd from the cytoplasm of L. chuanxiong to the vacuoles (Ueno et al., 2011; Sasaki et al., 2014). ABCG protein is one of the few known Cd efflux proteins that can eliminate Cd from root epidermal cells and reduce Cd toxicity in plants. For example, overexpression of the PtoABCG36 gene reduced Cd accumulation in Arabidopsis (Wang et al., 2019a). In our results, two ABCG protein subfamily genes were differentially expressed, and they may play an important role in Cd efflux in L. chuanxiong. The above findings highlight the critical roles of multiple metal transporter proteins in Cd transport and accumulation in L. chuanxiong. These findings help to deepen our understanding of the mechanisms of Cd tolerance and accumulation in L. chuanxiong.
5 Conclusions
This study analyzed the growth characteristics, physiological indices, and transcriptomic data of L. chuanxiong under different concentrations of Cd treatment. We explored the effects of Cd stress on the growth and development of L. chuanxiong seedlings and proposed insights into Cd accumulation and tolerance mechanisms in L. chuanxiong. The decrease in photosynthetic capacity and the regulation of plant hormone levels appear to be major factors contributing to growth inhibition in response to Cd stress. Moreover, the upregulation of differential genes involved in cell wall modification suggests a potential mechanism for enhancing root cell wall capabilities in isolating and sequestering Cd. The involvement of specific metal transport proteins further highlights their importance in Cd movement within the plant. This study enhances our understanding of the impacts of Cd stress on L. chuanxiong and provides a basis for future research in this field.
Data availability statement
The data presented in the study are deposited in the National Center for Biotechnology Information Biology Project (NCBI) repository, accession number PRJNA960886.
Author contributions
ZZ: Conceptualization, Data curation, Formal Analysis, Methodology, Resources, Software, Validation, Writing – original draft. LZ: Conceptualization, Data curation, Funding acquisition, Project administration, Writing – original draft. WX: Formal Analysis, Investigation, Methodology, Software, Writing – review & editing. YD: Formal Analysis, Investigation, Software, Writing – review & editing. GH: Methodology, Resources, Supervision, Writing – review & editing. ZY: Methodology, Resources, Conceptualization, Writing – review & editing. DH: Conceptualization, Funding acquisition, Project administration, Supervision, Writing – review & editing. CZ: Funding acquisition, Project administration, Resources, Investigation, Writing – review & editing.
Funding
The authors declare financial support was received for the research, authorship, and/or publication of this article. This work is supported by the Natural Science Foundation of Sichuan Province (2022NSFSC0101, 2023NSFSC1994) and the National Natural Science Foundation of China (81673553, 81001610).
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.1263981/full#supplementary-material
References
Agami, R. A., Mohamed, G. F. (2013). Exogenous treatment with indole-3-acetic acid and salicylic acid alleviates cadmium toxicity in wheat seedlings. Ecotoxicol. Environ. Saf. 94, 164–171. doi: 10.1016/j.ecoenv.2013.04.013
Chen, P., Li, Z., Luo, D., Jia, R., Lu, H., Tang, M., et al. (2021b). 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., Yang, R. X., Zhang, X., Chen, Y. H., Xia, Y., Xu, X. M. (2021a). Foliar application of gibberellin inhibits the cadmium uptake and xylem transport in lettuce (Lactuca sativa L.). Scientia Hortic. 288, 110410. doi: 10.1016/j.scienta.2021.110410
Chen, Z., Zhang, C., Gao, F., Fu, Q., Fu, C., He, Y., et al. (2018). A systematic review on the rhizome of Ligusticum chuanxiong Hort. (Chuanxiong). Food Chem. Toxicol. 119, 309–325. doi: 10.1016/j.fct.2018.02.050
Cosgrove, D. J. (2005). Growth of the plant cell wall. Nat. Rev. Mol. Cell Biol. 6 (11), 850–861. doi: 10.1038/nrm1746
Dai, F., Luo, G., Li, Z., Wei, X., Wang, Z., Lin, S., et al. (2020). Physiological and transcriptomic analyses of mulberry (Morus atropurpurea) response to cadmium stress. Ecotoxicol. Environ. Saf. 205, 111298. doi: 10.1016/j.ecoenv.2020.111298
Demidchik, V., Shabala, S., Isayenkov, S., Cuin, T. A., Pottosin, I. (2018). Calcium transport across plant membranes: mechanisms and functions. New Phytol. 220 (1), 49–69. doi: 10.1111/nph.15266
Demiral, T., Türkan, S. (2005). Comparative lipid peroxidation, antioxidant defense systems and proline content in roots of two rice cultivars differing in salt tolerance. Environ. Exp. Bot. 53 (3), 247–257. doi: 10.1016/j.envexpbot.2004.03.017
Deng, S., Yu, J., Wang, Y., Xie, S., Ran, Z., Wei, W. (2019). Distribution, transfer, and time-dependent variation of Cd in soil-rice system: A case study in the Chengdu plain, Southwest China. Soil Tillage Res. 195, 104367. doi: 10.1016/j.still.2019.104367
Dixit, P., Mukherjee, P. K., Ramachandran, V., Eapen, S. (2011). Glutathione transferase from Trichoderma virens enhances cadmium tolerance without enhancing its accumulation in transgenic Nicotiana tabacum. PloS One 6 (1), e16360. doi: 10.1371/journal.pone.0016360
Du, J., Wang, Z., Liu, J., Zhong, S., Wei, C. (2020). Distribution Characteristics of Soil Heavy Metals, their Source Identification and their Changes Influenced by Anthropogenic Cultivation Activities in Purple Hilly Regions of Sichuan Basin, China. J. Soil Sci. Plant Nutr. 20 (3), 1080–1091. doi: 10.1007/s42729-020-00194-1
Ehrnstorfer, I. A., Geertsma, E. R., Pardon, E., Steyaert, J., Dutzler, R. (2014). Crystal structure of a SLC11 (NRAMP) transporter reveals the basis for transition-metal ion transport. Nat. Struct. Mol. Biol. 21 (11), 990–996. doi: 10.1038/nsmb.2904
Ehsan, S., Ali, S., Noureen, S., Mahmood, K., Farid, M., Ishaque, W., et al. (2014). Citric acid assisted phytoremediation of cadmium by Brassica napus L. Ecotoxicol. Environ. Saf. 106, 164–172. doi: 10.1016/j.ecoenv.2014.03.007
Emamverdian, A., Ding, Y., Mokhberdoran, F. (2020). The role of salicylic acid and gibberellin signaling in plant responses to abiotic stress with an emphasis on heavy metals. Plant Signal Behav. 15 (7), 1777372. doi: 10.1080/15592324.2020.1777372
Fan, P., Wu, L., Wang, Q., Wang, Y., Luo, H., Song, J., et al. (2023). Physiological and molecular mechanisms of medicinal plants in response to cadmium stress: Current status and future perspective. J. Hazard Mater 450, 131008. doi: 10.1016/j.jhazmat.2023.131008
Farooq, M. A., Ali, S., Hameed, A., Bharwana, S. A., Rizwan, M., Ishaque, W., et al. (2016). Cadmium stress in cotton seedlings: Physiological, photosynthesis and oxidative damages alleviated by glycinebetaine. South Afr. J. Bot. 104, 61–68. doi: 10.1016/j.sajb.2015.11.006
Fattorini, L., Ronzan, M., Piacentini, D., Della Rovere, F., De Virgilio, C., Sofo, A., et al. (2017). Cadmium and arsenic affect quiescent centre formation and maintenance in Arabidopsis thaliana post-embryonic roots disrupting auxin biosynthesis and transport. Environ. Exp. Bot. 144, 37–48. doi: 10.1016/j.envexpbot.2017.10.005
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 (2), 558–571. doi: 10.1111/pbi.12795
Finger-Teixeira, A., Ferrarese Mde, L., Soares, A. R., da Silva, D., Ferrarese-Filho, O. (2010). Cadmium-induced lignification restricts soybean root growth. Ecotoxicol. Environ. Saf. 73 (8), 1959–1964. doi: 10.1016/j.ecoenv.2010.08.021
Fu, X. Z., Tong, Y. H., Zhou, X., Ling, L. L., Chun, C. P., Cao, L., et al. (2017). Genome-wide identification of sweet orange (Citrus sinensis) metal tolerance proteins and analysis of their expression patterns under zinc, manganese, copper, and cadmium toxicity. Gene 629, 1–8. doi: 10.1016/j.gene.2017.07.072
Godt, J., Scheidig, F., Grosse-Siestrup, C., Esche, V., Brandenburg, P., Reich, A., et al. (2006). The toxicity of cadmium and resulting hazards for human health. J. Occup. Med. Toxicol. 1 (1), 22. doi: 10.1186/1745-6673-1-22
Guerinot, M. L. (2000). The ZIP family of metal transporters. Biochim. Biophys. Acta 1465 (1-2), 190–198. doi: 10.1016/s0005-2736(00)00138-3
Guo, H., Hong, C., Xiao, M., Chen, X., Chen, H., Zheng, B., et al. (2016). Real-time kinetics of cadmium transport and transcriptomic analysis in low cadmium accumulator Miscanthus sacchariflorus. Planta 244 (6), 1289–1302. doi: 10.1007/s00425-016-2578-3
Guo, B., Liang, Y. C., Zhu, Y. G., Zhao, F. J. (2007). Role of salicylic acid in alleviating oxidative damage in rice roots (Oryza sativa) subjected to cadmium stress. Environ. pollut. 147 (3), 743–749. doi: 10.1016/j.envpol.2006.09.007
Guo, Z., Yuan, H., Yansong, A. O. (2009). Effect of cadmium on photosynthesis and nutrient elements uptake of Solanum nigrum L. seedlings. Ecol. Environ. Sci. 18 (03), 824–829. doi: 10.16258/j.cnki.1674-5906.2009.03.004
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
Halimaa, P., Blande, D., Baltzi, E., Aarts, M. G. M., Granlund, L., Keinänen, M., et al. (2019). Transcriptional effects of cadmium on iron homeostasis differ in calamine accessions of Noccaea caerulescens. Plant J. 97 (2), 306–320. doi: 10.1111/tpj.14121
Harholt, J., Suttangkakul, A., Vibe Scheller, H. (2010). Biosynthesis of pectin. Plant Physiol. 153 (2), 384–395. doi: 10.1104/pp.110.156588
Hasan, S. A., Fariduddin, Q., Ali, B., Hayat, S., Ahmad, A. (2009). Cadmium: toxicity and tolerance in plants. J. Environ. Biol. 30 (2), 165–174.
Huang, L., Li, W. C., Tam, N. F. Y., Ye, Z. (2019a). Effects of root morphology and anatomy on cadmium uptake and translocation in rice (Oryza sativa L.). J. Environ. Sci. 75, 296–306. doi: 10.1016/j.jes.2018.04.005
Huang, Y., Wang, L., Wang, W., Li, T., He, Z., Yang, X. (2019b). Current status of agricultural soil pollution by heavy metals in China: A meta-analysis. Sci. Total Environ. 651 (Pt 2), 3034–3042. doi: 10.1016/j.scitotenv.2018.10.185
Jia, L., Liu, Z., Chen, W., Ye, Y., Yu, S., He, X. (2015). Hormesis effects induced by cadmium on growth and photosynthetic performance in a hyperaccumulator, lonicera japonica thunb. J. Plant Growth Regul. 34 (1), 13–21. doi: 10.1007/s00344-014-9433-1
Jiang, W., Wu, Z., Wang, T., Mantri, N., Huang, H., Li, H., et al. (2020). Physiological and transcriptomic analyses of cadmium stress response in Dendrobium officinale seedling. Plant Physiol. Biochem. 148, 152–165. doi: 10.1016/j.plaphy.2020.01.010
Johri, N., Jacquillet, G., Unwin, R. (2010). Heavy metal poisoning: the effects of cadmium on the kidney. BioMetals 23 (5), 783–792. doi: 10.1007/s10534-010-9328-y
Kabała, K., Kłobus, G.y. (2005). Plant ca2+-ATPases. Acta Physiol. Plant. 27 (4), 559–574. doi: 10.1007/s11738-005-0062-y
Kocaman, A. (2023). Effects of foliar application of abscisic acid on antioxidant content, phytohormones in strawberry shoots, and translocation of various heavy metals. Scientia Hortic. 314, 111943. doi: 10.1016/j.scienta.2023.111943
Krzesłowska, M. (2011). The cell wall in plant cell response to trace metals: polysaccharide remodeling and its role in defense strategy. Acta Physiol. Plant. 33 (1), 35–51. doi: 10.1007/s11738-010-0581-z
Le Gall, H., Philippe, F., Domon, J. M., Gillet, F., Pelloux, J., Rayon, C. (2015). Cell wall metabolism in response to abiotic stress. Plants (Basel) 4 (1), 112–166. doi: 10.3390/plants4010112
Li, B., Dewey, C. N. (2011). RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinf. 12, 323. doi: 10.1186/1471-2105-12-323
Li, J.-R., Liu, C.-C., Sun, C.-H., Chen, Y.-T. (2018a). Plant stress RNA-seq Nexus: a stress-specific transcriptome database in plant cells. BMC Genomics 19 (1), 966. doi: 10.1186/s12864-018-5367-5
Li, X., Mao, X., Xu, Y., Li, Y., Zhao, N., Yao, J., et al. (2021). Comparative transcriptomic analysis reveals the coordinated mechanisms of Populus x canadensis ‘Neva’ leaves in response to cadmium stress. Ecotoxicol. Environ. Saf. 216, 112179. doi: 10.1016/j.ecoenv.2021.112179
Li, Q., Wang, C., Dai, T., Shi, W., Zhang, X., Xiao, Y., et al. (2017). Prediction of soil cadmium distribution across a typical area of Chengdu Plain, China. Sci. Rep. 7 (1), 7115. doi: 10.1038/s41598-017-07690-y
Li, Y., Zhang, S., Bao, Q., Chu, Y., Sun, H., Huang, Y. (2022). Jasmonic acid alleviates cadmium toxicity through regulating the antioxidant response and enhancing the chelation of cadmium in rice (Oryza sativa L.). Environ. Pollut. 304, 119178. doi: 10.1016/j.envpol.2022.119178
Li, X., Zhang, X., Wu, Y., Li, B., Yang, Y. (2018b). Physiological and biochemical analysis of mechanisms underlying cadmium tolerance and accumulation in turnip. Plant Divers. 40 (1), 19–27. doi: 10.1016/j.pld.2017.12.005
Liu, Z., Ding, Y., Wang, F., Ye, Y., Zhu, C. (2016). Role of salicylic acid in resistance to cadmium stress in plants. Plant Cell Rep. 35 (4), 719–731. doi: 10.1007/s00299-015-1925-3
Liu, X., Peng, K., Wang, A., Lian, C., Shen, Z. (2010). Cadmium accumulation and distribution in populations of Phytolacca Americana L. and the role of transpiration. Chemosphere 78 (9), 1136–1141. doi: 10.1016/j.chemosphere.2009.12.030
Liu, C. P., Shen, Z. G., Li, X. D. (2007). Accumulation and detoxification of cadmium in Brassica pekinensis and B. chinensis. Biol. Plant. 51 (1), 116–120. doi: 10.1007/s10535-007-0023-y
Long, L., Persson, D. P., Duan, F., Jørgensen, K., Yuan, L., Schjoerring, J. K., et al. (2018). The iron-regulated transporter 1 plays an essential role in uptake, translocation and grain-loading of manganese, but not iron, in barley. New Phytol. 217 (4), 1640–1653. doi: 10.1111/nph.14930
Ma, Y., He, J., Ma, C., Luo, J., Li, H., Liu, T., et al. (2014). Ectomycorrhizas with Paxillus involutus enhance cadmium uptake and tolerance in Populus × canescens. Plant Cell Environ. 37 (3), 627–642. doi: 10.1111/pce.12183
Mendoza-Cózatl, D. G., Zhai, Z., Jobe, T. O., Akmakjian, G. Z., Song, W.-Y., Limbo, O., et al. (2010). Tonoplast-localized abc2 transporter mediates phytochelatin accumulation in vacuoles and confers cadmium tolerance*. J. Biol. Chem. 285, 40416–40426. doi: 10.1074/jbc.M110.155408
Nagajyoti, P. C., Lee, K. D., Sreekanth, T. V. M. (2010). Heavy metals, occurrence and toxicity for plants: a review. Environ. Chem. Lett. 8 (3), 199–216. doi: 10.1007/s10311-010-0297-8
Nawrot, T., Plusquin, M., Hogervorst, J., Roels, H. A., Celis, H., Thijs, L., et al. (2006). Environmental exposure to cadmium and risk of cancer: a prospective population-based study. Lancet Oncol. 7 (2), 119–126. doi: 10.1016/s1470-2045(06)70545-9
Nguyen, T. Q., Sesin, V., Kisiala, A., Emery, R. J. N. (2021). Phytohormonal roles in plant responses to heavy metal stress: implications for using macrophytes in phytoremediation of aquatic ecosystems. Environ. Toxicol. Chem. 40 (1), 7–22. doi: 10.1002/etc.4909
Oomen, R. J., Wu, J., Lelièvre, F., Blanchet, S., Richaud, P., Barbier-Brygoo, H., et al. (2009). Functional characterization of NRAMP3 and NRAMP4 from the metal hyperaccumulator Thlaspi caerulescens. New Phytol. 181 (3), 637–650. doi: 10.1111/j.1469-8137.2008.02694.x
Ortiz, D. F., Kreppel, L., Speiser, D. M., Scheel, G., McDonald, G., Ow, D. W. (1992). Heavy metal tolerance in the fission yeast requires an ATP-binding cassette-type vacuolar membrane transporter. EMBO J. 11 (10), 3491–3499. doi: 10.1002/j.1460-2075.1992.tb05431.x
Ortiz, D. F., Ruscitti, T., McCue, K. F., Ow, D. W. (1995). Transport of metal-binding peptides by HMT1, a fission yeast ABC-type vacuolar membrane protein. J. Biol. Chem. 270 (9), 4721–4728. doi: 10.1074/jbc.270.9.4721
Pandey, A., Khan, M. K., Hamurcu, M., Brestic, M., Topal, A., Gezgin, S. (2022). Insight into the Root Transcriptome of a Boron-Tolerant Triticum zhukovskyi Genotype Grown under Boron Toxicity. Agronomy 12 (10), 2421. doi: 10.3390/agronomy12102421
Pilon, M., Cohu, C. M., Ravet, K., Abdel-Ghany, S. E., Gaymard, F. (2009). Essential transition metal homeostasis in plants. Curr. Opin. Plant Biol. 12 (3), 347–357. doi: 10.1016/j.pbi.2009.04.011
Piñeros, M. A., Shaff, J. E., Kochian, L. V. (1998). Development, characterization, and application of a cadmium-selective microelectrode for the measurement of cadmium fluxes in roots of thlaspi species and wheat. Plant Physiol. 116 (4), 1393–1401. doi: 10.1104/pp.116.4.1393
Qin, H., Huang, R. (2018). Auxin controlled by ethylene steers root development. Int. J. Mol. Sci. 19 (11), 3656. doi: 10.3390/ijms19113656
Qi-qua, L. (2014). Contents and sources of cadmium in farmland soils of Chengdu plain, China. J. Agro-Environ. Sci. 33 (05), 898–906. doi: 10.11654/jaes.2014.05.011
Rajalakshmi, K., Banu, N. (2015). Extraction and estimation of chlorophyll from medicinal plants. Int. J. Sci. Res. 4 (11), 209–212. doi: 10.21275/v4i11.nov151021
Ren, S., Song, C., Ye, S., Cheng, C., Gao, P. (2022). The spatiotemporal variation in heavy metals in China’s farmland soil over the past 20 years: A meta-analysis. Sci. Total Environ. 806 (Pt 2), 150322. doi: 10.1016/j.scitotenv.2021.150322
Rivetta, A., Pesenti, M., Sacchi, G. A., Nocito, F. F., Cocucci, M. (2023). Cadmium transport in maize root segments using a classical physiological approach: evidence of influx largely exceeding efflux in subapical regions. Plants (Basel) 12 (5), 992. doi: 10.3390/plants12050992
Rizwan, M., Ali, S., Abbas, T., Zia-Ur-Rehman, M., Hannan, F., Keller, C., et al. (2016). Cadmium minimization in wheat: A critical review. Ecotoxicol. Environ. Saf. 130, 43–53. doi: 10.1016/j.ecoenv.2016.04.001
Rui, H., Chen, C., Zhang, X., Shen, Z., Zhang, F. (2016). Cd-induced oxidative stress and lignification in the roots of two Vicia sativa L. varieties with different Cd tolerances. J. Hazard Mater 301, 304–313. doi: 10.1016/j.jhazmat.2015.08.052
Růzicka, K., Ljung, K., Vanneste, S., Podhorská, R., Beeckman, T., Friml, J., et al. (2007). Ethylene regulates root growth through effects on auxin biosynthesis and transport-dependent auxin distribution. Plant Cell 19 (7), 2197–2212. doi: 10.1105/tpc.107.052126
Sasaki, A., Yamaji, N., Ma, J. F. (2014). Overexpression of OsHMA3 enhances Cd tolerance and expression of Zn transporter genes in rice. J. Exp. Bot. 65 (20), 6013–6021. doi: 10.1093/jxb/eru340
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 (5), 2155–2167. doi: 10.1105/tpc.112.096925
Schellingen, K., van der Straeten, D., Remans, T., Loix, C., Vangronsveld, J., Cuypers, A. (2015). Ethylene biosynthesis is involved in the early oxidative challenge induced by moderate Cd exposure in Arabidopsis thaliana. Environ. Exp. Bot. 117, 1–11. doi: 10.1016/j.envexpbot.2015.04.005
Schützendübel, A., Schwanz, P., Teichmann, T., Gross, K., Langenfeld-Heyser, R., Godbold, D. L., et al. (2001). Cadmium-induced changes in antioxidative systems, hydrogen peroxide content, and differentiation in Scots pine roots. Plant Physiol. 127 (3), 887–898. doi: 10.1104/pp.010318
Singh, S., Prasad, S. M. (2015). IAA alleviates Cd toxicity on growth, photosynthesis and oxidative damages in eggplant seedlings. Plant Growth Regul. 77 (1), 87–98. doi: 10.1007/s10725-015-0039-9
Singh, I., Shah, K. (2014). Exogenous application of methyl jasmonate lowers the effect of cadmium-induced oxidative injury in rice seedlings. Phytochemistry 108, 57–66. doi: 10.1016/j.phytochem.2014.09.007
Smyth, G. K. (2005). “limma: Linear Models for Microarray Data,” in Bioinformatics and Computational Biology Solutions Using R and Bioconductor. Eds. Gentleman, R., Carey, V. J., Huber, W., Irizarry, R. A., Dudoit, S. (New York, NY: Springer New York), 397–420.
Song, X. Q., Liu, L. F., Jiang, Y. J., Zhang, B. C., Gao, Y. P., Liu, X. L., et al. (2013). Disruption of secondary wall cellulose biosynthesis alters cadmium translocation and tolerance in rice plants. Mol. Plant 6 (3), 768–780. doi: 10.1093/mp/sst025
Staessen, J. A., Roels, H. A., Emelianov, D., Kuznetsova, T., Thijs, L., Vangronsveld, J., et al. (1999). Environmental exposure to cadmium, forearm bone density, and risk of fractures: prospective population study. Public Health and Environmental Exposure to Cadmium (PheeCad) Study Group. Lancet 353 (9159), 1140–1144. doi: 10.1016/s0140-6736(98)09356-8
Sterckeman, T., Thomine, S. (2020). Mechanisms of cadmium accumulation in plants. Crit. Rev. Plant Sci. 39 (4), 322–359. doi: 10.1080/07352689.2020.1792179
Tenhaken, R. (2014). Cell wall remodeling under abiotic stress. Front. Plant Sci. 5. doi: 10.3389/fpls.2014.00771
Ueno, D., Milner, M. J., Yamaji, N., Yokosho, K., Koyama, E., Clemencia Zambrano, M., et al. (2011). Elevated expression of TcHMA3 plays a key role in the extreme Cd tolerance in a Cd-hyperaccumulating ecotype of Thlaspi caerulescens. Plant J. 66 (5), 852–862. doi: 10.1111/j.1365-313X.2011.04548.x
Wakeel, A., Ali, I., Upreti, S., Azizullah, A., Liu, B., Khan, A. R., et al. (2018). Ethylene mediates dichromate-induced inhibition of primary root growth by altering AUX1 expression and auxin accumulation in Arabidopsis thaliana. Plant Cell Environ. 41 (6), 1453–1467. doi: 10.1111/pce.13174
Wang, M. Q., Bai, Z. Y., Xiao, Y. F., Li, Y., Liu, Q. L., Zhang, L., et al. (2019b). Transcriptomic analysis of Verbena bonariensis roots in response to cadmium stress. BMC Genomics 20 (1), 877. doi: 10.1186/s12864-019-6152-9
Wang, X., Li, N., Li, W., Gao, X., Cha, M., Qin, L., et al. (2020). Advances in transcriptomics in the response to stress in plants. Glob Med. Genet. 7 (2), 30–34. doi: 10.1055/s-0040-1714414
Wang, H., Liu, Y., Peng, Z., Li, J., Huang, W., Liu, Y., et al. (2019a). Ectopic expression of poplar ABC transporter ptoABCG36 confers cd tolerance in arabidopsis thaliana. Int. J. Mol. Sci. 20 (13), 3293. doi: 10.3390/ijms20133293
Weigel, H. J., Jäger, H. J. (1980). Subcellular distribution and chemical form of cadmium in bean plants. Plant Physiol. 65 (3), 480–482. doi: 10.1104/pp.65.3.480
Wu, X., Chen, J., Yue, X., Wei, X., Zou, J., Chen, Y., et al. (2019). The zinc-regulated protein (ZIP) family genes and glutathione s-transferase (GST) family genes play roles in Cd resistance and accumulation of pak choi (Brassica campestris ssp. chinensis). Ecotoxicol. Environ. Saf. 183, 109571. doi: 10.1016/j.ecoenv.2019.109571
Wu, F., Fan, J., Ye, X., Yang, L., Hu, R., Ma, J., et al. (2022). Unraveling Cadmium Toxicity in Trifolium repens L. Seedling: Insight into Regulatory Mechanisms Using Comparative Transcriptomics Combined with Physiological Analyses. Int. J. Mol. Sci. 23 (9), 4612. doi: 10.3390/ijms23094612
Xiao, X., Pei, M., Zhou, J., Sun, S., Li, C., Zhu, X., et al. (2021). Soil amendments inhibited the cadmium accumulation in Ligusticum striatum DC. and improved the plant growth. Environ. Sci. Pollut. Res. Int. 28 (47), 67788–67799. doi: 10.1007/s11356-021-15332-0
Xie, C., Mao, X., Huang, J., Ding, Y., Wu, J., Dong, S., et al. (2011). KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 39 (Web Server issue), W316–W322. doi: 10.1093/nar/gkr483
Xiong, J., An, L., Lu, H., Zhu, C. (2009). Exogenous nitric oxide enhances cadmium tolerance of rice by increasing pectin and hemicellulose contents in root cell wall. Planta 230 (4), 755–765. doi: 10.1007/s00425-009-0984-5
Yadav, S. K. (2010). Heavy metals toxicity in plants: An overview on the role of glutathione and phytochelatins in heavy metal stress tolerance of plants. South Afr. J. Bot. 76 (2), 167–179. doi: 10.1016/j.sajb.2009.10.007
Yang, J., Li, L., Zhang, X., Wu, S., Han, X., Li, X., et al. (2022). Comparative Transcriptomics Analysis of Roots and Leaves under Cd Stress in Calotropis gigantea L. Int. J. Mol. Sci. 23 (6), 3329. doi: 10.3390/ijms23063329
Yuan, H. M., Huang, X. (2016). Inhibition of root meristem growth by cadmium involves nitric oxide-mediated repression of auxin accumulation and signalling in Arabidopsis. Plant Cell Environ. 39 (1), 120–135. doi: 10.1111/pce.12597
Zacchini, M., Pietrini, F., Scarascia Mugnozza, G., Iori, V., Pietrosanti, L., Massacci, A. (2009). Metal tolerance, accumulation and translocation in poplar and willow clones treated with cadmium in hydroponics. Water Air Soil Pollut. 197 (1), 23–34. doi: 10.1007/s11270-008-9788-7
Zeng, J., Li, X., Wang, X., Zhang, K., Wang, Y., Kang, H., et al. (2020). Cadmium and lead mixtures are less toxic to the Chinese medicinal plant Ligusticum chuanxiong Hort. Than either metal alone. Ecotoxicol. Environ. Saf. 193, 110342. doi: 10.1016/j.ecoenv.2020.110342
Zhang, T., Li, R., Xing, J., Yan, L., Wang, R., Zhao, Y. (2018). The YUCCA-auxin-WOX11 module controls crown root development in rice. Front. Plant Sci. 9. doi: 10.3389/fpls.2018.00523
Zhao, Y. Y., He, C. X., Wu, Z. C., Liu, X. W., Cai, M. M., Jia, W., et al. (2019b). Selenium reduces cadmium accumulation in seed by increasing cadmium retention in root of oilseed rape (Brassica napus L.). Environ. Exp. Bot. 158, 161–170. doi: 10.1016/j.envexpbot.2018.11.017
Zhao, F. Y., Liu, W., Zhang, S. Y. (2009). Different responses of plant growth and antioxidant system to the combination of cadmium and heat stress in transgenic and non-transgenic rice. J. Integr. Plant Biol. 51 (10), 942–950. doi: 10.1111/j.1744-7909.2009.00865.x
Zhao, D., Zhang, X., Fang, Z., Wu, Y., Tao, J. (2019a). Physiological and transcriptomic analysis of tree peony (Paeonia section moutan DC.) in response to drought stress. Forests 10 (2), 135. doi: 10.3390/f10020135
Zheng, Y., Zhang, R., Zhu, Y., Ao, Q., Liu, H., Li, A., et al. (2022). Salicylic acid improves Nasturtium officinale phytoremediation capability for cadmium-contaminated paddy soils. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.1059175
Zhou, Q., Guo, J. J., He, C. T., Shen, C., Huang, Y. Y., Chen, J. X., et al. (2016). Comparative transcriptome analysis between low- and high-cadmium-accumulating genotypes of pakchoi (Brassica chinensis L.) in response to cadmium stress. Environ. Sci. Technol. 50 (12), 6485–6494. doi: 10.1021/acs.est.5b06326
Keywords: Ligusticum chuanxiong Hort., cadmium stress, subcellular localization, antioxidant system, transcriptomics
Citation: Zhang Z, Zhong L, Xiao W, Du Y, Han G, Yan Z, He D and Zheng C (2023) Transcriptomics combined with physiological analysis reveals the mechanism of cadmium uptake and tolerance in Ligusticum chuanxiong Hort. under cadmium treatment. Front. Plant Sci. 14:1263981. doi: 10.3389/fpls.2023.1263981
Received: 20 July 2023; Accepted: 30 August 2023;
Published: 22 September 2023.
Edited by:
Debasis Chakrabarty, CSIR-National Botanical Research Institute, IndiaReviewed by:
Nitika Kapoor, Guru Nanak Dev University, IndiaVaishali Yadav, Chaudhary Charan Singh University, India
Copyright © 2023 Zhang, Zhong, Xiao, Du, Han, Yan, He and Zheng. 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: Dongmei He, aGVkb25nbWVpQGNkdXRjbS5lZHUuY24=
†These authors have contributed equally to this work and share first authorship