Skip to main content

Physiological of biochar and α-Fe2O3 nanoparticles as amendments of Cd accumulation and toxicity toward muskmelon grown in pots

Abstract

Background

Due to the severe cadmium (Cd) pollution of farmland soil, effective measures need to be taken to reduce the Cd content in agricultural products. In this study, we added α-Fe2O3 nanoparticles (NPs) and biochar into Cd-contaminated soil to investigate physiological responses of muskmelon in the whole life cycle.

Results

The results showed that Cd caused adverse impacts on muskmelon (Cucumis melo) plants. For instance, the chlorophyll of muskmelon leaves in the Cd alone treatment was reduced by 8.07–32.34% in the four periods, relative to the control. The treatments with single amendment, α-Fe2O3 NPs or 1% biochar or 5% biochar, significantly reduced the soil available Cd content, but the co-exposure treatments (α-Fe2O3 NPs and biochar) had no impact on the soil available Cd content. All treatments could reduce the Cd content by 47.64–74.60% and increase the Fe content by 15.15–95.27% in fruits as compared to the Cd alone treatment. The KEGG enrichment results of different genes in different treatments indicated that single treatments could regulate genes related to anthocyanin biosynthesis, glutathione metabolism and MAPK signal transduction pathways to reduce the Cd toxicity.

Conclusions

Overall the combination of biochar and α-Fe2O3 NPs can alleviate Cd toxicity in muskmelon. The present study could provide new insights into Cd remediation in soil using α-Fe2O3 NPs and biochar as amendments.

Graphic Abstract

Synopsis

All the single or co-treatments could reduce the Cd content and increase the Fe content in fruit.

Introduction

The industry development and human activities have made tremendous contributions to improving living standards, but they have also caused environmental problems, such as the accumulation of heavy metals in soil and water resources [1]. Nowadays, Heavy metal pollution continues to be one of the most serious environmental problems, which has attracted great attention [2]. In China, about 1/5 of the soil has been contaminated with heavy metals [3]. Cadmium (Cd) is an injurious heavy metal that has no biological importance [4]. Cd is extremely toxic to plants and can subsequently cause crop reduction in agriculture [5, 6]. In the soil, plants absorb Cd via roots, and these Cds eventually enter the human body through the food chain, which pose a great threat to human health [7]. For example, the prevalence of itai-itai disease in Japan was caused by Cd pollution [8]. Effective methods are urgently needed to remediate Cd contamination.

In-situ remediation technology is an effective method to reduce Cd pollution because of low cost and high efficiency [9]. Many remediating agents can be used to improve soil Cd pollution, such as Cellulose-zeolitic imidazolate frameworks (CelloZIFs), polylactic acid (PLA)-hydroxyapatite (HAp) composite, biochar, etc [10,11,12]. Biochar is a carbon-enriched material synthesized by pyrolysis of plant materials, which can be obtained from agricultural by-products [13]. Biochar has the characteristics of porous structure and high pH and because of its harmlessness to the environment and ability to improve crop productivity, it has attracted attentions in agricultural applications [14]. For example, biochar has great potential in reducing the effectiveness of soil Cd migration [13]. In addition, it was also demonstrated that biochar increased the biomass of Cd treated spinach by reducing the Cd accumulation in spinach [15]. Similarly, biochar reduced the Cd accumulation in super japonica rice and lettuce (Lactuca sativa) grown in Cd contaminated soil [16,17,18]. Therefore, biochar has considerable potential in remediating Cd contamination in soil.

Iron is an essential trace element for plant growth, which is closely related to the formation of chlorophyll in photosynthesis [19]. Due to their large surface area and high catalytic activity, nanomaterials have been widely used in industry and agriculture [20]. Fe2O3 nanoparticles (NPs) can be used as potential iron fertilizers in agriculture [21]. Recent studies have shown that iron oxide NPs also have a certain effect on remediating heavy metal polluted environments. For example, each gram of the hollow flower-like α-Fe2O3 nanostructures can adsorb 39 mg of Pb2+ or 48 mg of Cu2+ [22]. Even in comparison with traditional mercury adsorption methods, Fe3O4 NPs functionalized with gold showed a high adsorption affinity for mercury, with an adsorption capacity of 79.59 mg/g [23].

Magnetite (γ-Fe2O3) and hematite (α-Fe2O3) are natural forms of iron oxide in nature [24]. Previous studies have shown that α-Fe2O3 NPs have less negative effects on the growth and physiology of melon plants than γ-Fe2O3 NPs because of the excessive adsorption of γ-Fe2O3 NPs onto soil colloids with subsequent low water extractable iron [25, 26]. If α-Fe2O3 NPs can also reduce Cd pollution in the soil, and thus it can be used as an iron fertilizer or a remediating agent in the soil. It is reported that combining biochar and inorganic remediating agents such as zero valent iron and sulfur-based amendments may be better than individual applications [18]. However, effects of combination of Fe2O3 NPs and biochar on alleviating Cd-induced toxicity to terrestrial plants is still unclear.

Muskelon (Cucumis melo) is an important economic crop because its edible portion and seeds have medicinal values [27, 28]. The present experiment was conducted under greenhouse conditions to investigate the role of co-existing biochar and α-Fe2O3 NPs in altering Cd accumulation and toxicity to muskmelon plants, and further to optimize the combining ratio of the two amendments for soil remediation. The experiment was lasted until the maturity of muskmelon. At harvest, physiological and biochemical endpoints were measured across all treatments. Overall, these findings could provide important information for the restoration of soil Cd pollution and the safe cultivation of muskmelon plants on Cd-contaminated substrates.

Materials and methods

Biochar synthesis and material characterization

The soil is ordinary bagged matrix soil, and its physical and chemical properties are listed in Additional file 1: Table S1. α-Fe2O3 NPs (60 nm) were purchased from Macklin, Shanghai, China. Muskelon (Cucumis melo) seeds were obtained from Hubei Academy of Agricultural Sciences, China.

For biochar synthesis, the air-dried straw was ground into powder. The raw materials were pyrolyzed in a tube furnace under N2 atmosphere at 500 °C for 5 h. The obtained biochar was ground and passed through a 40-mesh screen prior to characterization. The scanning electron microscope (SEM; S-4800, Japan) images of biochar and α-Fe2O3 NPs are shown in Fig. 1. The XRD patterns of α-Fe2O3 NPs and biochar are shown in Additional file 1: Figure S1. Briefly, Crystallinity of the α-Fe2O3 NPs and biochar were tested by producing X-ray diffraction (XRD) patterns with an X-ray diffractometer (XRD; Bruker D8 Advance, Germany). The 2θ was within range of 10°–70°. Ni-filtered Cu Ka radiation was used as the X-ray source. Voltage and current were 40 kV and 40 mA, respectively. The hydrodynamic diameter and zeta potential of α-Fe2O3 NPs and the FTIR results of biochar were described in our previous study [24, 25].

Fig. 1
figure 1

SEM images of α-Fe2O3 NPs (A) and biochar (B)

Greenhouse experimental design

A greenhouse experiment was designed to investigate individual and combined effects of Fe2O3 NPs and biochar on Cd treated muskmelon. Our previous results have demonstrated that muskmelon is a Cd-tolerant plant; thus, the selected exposure concentrations of Cd and the amendments, Fe2O3 NPs and biochar, were increased [25]. Details of each treatment are listed in Table 1. Muskmelon seeds were germinated in deionized water. Then, seedlings with uniform size were transplanted into pots containing matrix soil amended with different analytes to make the designated treatment. Muskmelon leaves at the same position in each seedling were sampled to analyze physiological and biochemical indicators. In addition, plant height, hypocotyl thickness, and leaf length and leaf width of muskmelon plants were measured at week 3, 5, 9 and 10. The size and weight of the muskmelon fruit was measured at harvest, and soil samples across all the treatments were also collected for Cd and Fe content analysis.

Table 1 Details of nine treatments in the pot experiment

Physiological measurements

Plant height, hypocotyl thickness, leaf length and leaf width were measured at week 3, 5, 9 and 10. At each testing point, the leaves at node 8, 12, 16 and 20 of each muskmelon were sampled for physiological and biochemical analysis.

At harvest, the fresh weight of fruits was measured using a FA1004C electronic analytical balance (Shanghai Yueping Scientific Instrument Co., Ltd, China). The vertical diameter, horizontal diameter and peel thickness of the fruits were also measured across all the treatment. The sugar content of muskmelon fruits was measured using a PAL-1 sugar detector (ATAG0).

For protein and MDA content measurements, approximately 0.6 g leaves were ground in phosphate buffer (PBS, 0.05 mM, pH 7.8), and then the mixtures were centrifuged at 4000 rpm, 4 ºC for 10 min. The supernatant was diluted to 10 mL with PBS (pH 7.8). The soluble protein content was measured using Coomassie Brilliant Blue G250 at 595 nm [24]. The MDA content was determined by thiobarbituric acid at 450, 532 and 600 nm using a UV-752 N spectrophotometer (Shanghai Precision Scientific instrument Co Ltd, China) [25].

Pigment measurement

Chlorophyll content was measured using a SPAD-502 Plus chlorophyll meter (Konica Minolta). The leaves at three identical positions of each seedling were selected to measure the pigment content. For anthocyanin content measurements, approximately 0.1 g leaf tissues were ground in 1% (v:v) HCl (35%): ethanol (95%) solution. After the solution sits for 30 min, the mixtures were centrifuged at 4000 rpm for 10 min. The supernatant was diluted to 10 mL with 1% (v:v) HCl (35%): ethanol (95%) solution. The anthocyanin content was measured at 530 and 657 nm [29].

Antioxidant enzyme activity

The antioxidant enzyme extracts were the same as used in the MDA content measurement. Superoxide dismutase (SOD) activity was determined using riboflavin and nitrogen blue tetrazole at 560 nm. For the peroxidase (POD) activity, H2O2 was used to initiate the reaction between POD and guaiacol, the absorbance of generated brown substance was measured at 470 nm every 30 s over 4.5 min, and the slope of the increased absorbance readings was used to calculate the POD activity. The catalase (CAT) activity was measured by calculating the H2O2 decomposition at 240 nm every 10 s over 1.5 min, the slope of the decreased absorbance readings was used to calculate the CAT activity [30]. For the Ascorbic acid peroxidase (APX) activity, APX was used to catalyze the reaction between H2O2 and ascorbate (ASA), the absorbance of remaining ASA was measured at 290 nm every 10 s over 1.5 min, and the slope of the decreased absorbance readings was used to calculate the APX activity [31].

Cd and Fe content measurement

Plant tissues and soil samples were oven-dried and ground into fine powder. A weight of 0.15 mg sample was added into a tube containing 3 mL nitric acid and then heated at 100℃ for 1 h. After cooling to ambient temperature, 0.5 mL H2O2 were added for another half an hour at 100℃ to complete the digestion. The digesta were filtered and diluted with deionized water. For the available Cd in the soil, the soil samples were extracted with diethylenetriamine pentaacetic acid (DTPA), and then acidified with nitric acid for analysis [32]. The Cd and Fe concentration in digesta were determined by Avanta M atomic absorption spectrophotometer (GBC, Australia) [33].

SEM–EDS analysis of Cd and Fe element in muskmelon tissues

Scanning electron microscope-energy dispersive spectrometer (SEM–EDS) (SEM, S-4800, Japan; EDS, GENESIS, USA) was used to further verify the presence of Fe and Cd in leaves and fruits. The freeze-dried leaves and fruit tissues were ground into fine powder and then mounted onto a sample stage for SEM–EDS analysis.

RNA-Seq analysis of leaves upon exposure to Cd and two amendments

The 5-week old Cd treated muskmelon leaves were selected for transcriptomic analysis as affected by biochar and α-Fe2O3 NPs. Three biological replicates in the control and each treatment were applied. Thus, total number of 27 samples from 9 treatments were analyzed for transcriptome sequencing at Wuhan Genoseq Technology Co., Ltd., China. The total RNA in the sample was extracted using RNAprep Pure Plant kit (Bio TeKe, Beijing, China). The purity and concentration of RNA samples were determined using NanoDrop 2000 spectrophotometer (Thermo Scientific, Wilmington, DE, USA). Then the same amount of RNA from each sample was combined to construct a cDNA library. After digestion, the Hiseq 4000 (Illumina, San Diego, CA, USA) platform was used to sequence the first strand of cDNA. After data quality control, the clean data were aligned to the reference genome sequences of the muskmelon genome DHL92-v3.6.1. The differentially expressed genes (DEG) were identified in the R software through deseq2, with an FDR (false discovery rate) 0.05 and |log2 (fold change)|> 1 between two samples. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of differentially expressed genes were carried out by phyper in the R software, and the 20 most significant enrichment pathways were selected to draw the KEGG enrichment scatter plot.

Statistical analysis

For each assay, the value is expressed as mean ± standard deviation (SD). A one-way analysis of variance followed by Duncan's multiple comparisons was used to test the significance at p < 0.05 (IBM SPSS version 22). Principal component analysis is a collection of the results of all the previous detection indicators, from which two main components are extracted to analyze the difference across all the treatments, and its visualization was completed by R 4.0.0. According to the results of transcriptome sequencing, the differentially expressed genes (DEG) were screened out according to the screening conditions of false discovery rate (FDR) < 0.05 and log2(fold change) 1, the DEG were enriched according to the different pathway, and then used R 4.0.0 to complete the visual drawing.

Results and discussion

Cd accumulation in plant tissues and soil

The Cd content in muskmelon fruits of all treatments except control was higher than the national standard concentration of Cd (0.05 mg/kg) because of the extremely high Cd dose used in the present study. The Cd content (Fig. 2) decreased following the order: soil > leaves > fruits. In the Cd-alone treatment, the fruit Cd content was 4.7% of the soil Cd content. In the FeNPs, BC and BFNC treatments, the fruit Cd content was only 1.19–2.46% of the soil Cd content, suggesting that the single or co-treatments could significantly lower the Cd accumulation in the edible portion.

Fig. 2
figure 2

Cd accumulation in soil and muskmelon tissues. Phenotypic images of muskmelon fruits across all the treatments (A), SEM scanning images of muskmelon fruits in Cd-alone treatment (B), EDS analysis of Cd in muskmelon leaves in Cd-alone treatment (C) and FeNPs treatment (D), Cd content in muskmelon fruits (E), leaves (F) and soil (G) across all the treatments. Each bar represents the mean ± SD of three replicates. Different lowercase letters indicate significant differences at p < 0.05

In addition, the biochar (1% or 5%) or α-Fe2O3 NPs treatment significantly reduced soil available Cd content, thereby reducing the Cd uptake by plants. The amendment of biochar potentially altered the hydrolysis balance of Cd and stabilized the soil Cd via complexation and precipitation [34]. Similarly, Fan et al. found that the main mechanism of the Cd adsorption process by biochar involved ion exchange, surface complexation, electrostatic interaction, surface co-precipitation, and other mechanisms [35]. After adding the nano iron oxide, Cd could bind directly to hydroxyl surface groups of Fe oxides [36], In addition, the iron plaque on root surfaces could act as a barrier to prevent contaminants from entering plant bodies [37, 38]. Thus, the effect of α-Fe2O3 NPs was well, too. However, the co-treatments with biochar and α-Fe2O3 NPs did not show an additive effect on the reduction of the soil available Cd content. Wu et al. (2019) reported that a portion of the micropores in the straw biochar was filled with γ-Al2O3 particles [39]. It is noted that the pore size of biochar (Fig. 1) was in micrometer, which was larger than the size of α-Fe2O3 NPs. Thus, it is reasonable that these large pores could be filled with nano-sized particles, and we also confirmed this view by SEM (Additional file 1: Figure S2). The presence of α-Fe2O3 NPs lowered the efficiency of biochar on the adsorption of the soil available Cd. In the BFNC treatments, the addition of 5% biochar reduced the Cd content in leaves to the similar level of the BC treatments because that the excess biochar might offset the negative effects caused by α-Fe2O3 NP adsorption by the biochar.

A decreasing pattern of the Cd content from leaves to fruits was evident (Fig. 2), indicating that more Cd was stored in the leaves. It is well-known that vacuole is the largest cell organelle and can store toxins to lower their toxicity to plants [40]. Previous study demonstrated that Cd complexed with chelatin was stored in the vacuoles, resulting the high concentration of Cd in the vacuole than in other cell organelles [41]. Energy-dispersive spectrum (EDS) analysis (Additional file 1: Figure S3 and Fig. 2C, D) further confirmed the Cd accumulation in the muskmelon fruits.

Fe accumulation in muskmelon plants and soil

The Cd addition inhibited the Fe accumulation in muskmelon fruits (Additional file 1: Figure S4). However, across all the treatments with single or co-treatments, the Fe level in the Cd treated fruits was significantly increased when comparing to the Cd-alone treatment. It could be ascribed as that the competition between Cd and Fe was determined by the regulation of metal-related transporters [42]. Cd-induced mineral nutrient displacement could further affect plant metabolisms [43]. All the single or co-treatments significantly increased the Fe content in the fruits as compared to the Cd-alone treatment.

Muskmelon growth and fruit traits

During the whole growth cycle, the plant height in all the treatments was significantly lower than the control at the early stage. With the growth of muskmelon plants, the plant height in all single or co-treatments was recovered to a certain extent compared to the control (Fig. 3A). The leaf width in all the treatments was significantly lower than that of the control during the whole growth period (Fig. 3C). Similarly, other physiological parameters, including fresh weight, sugar content, transverse diameter, and peel thickness, were all lower as compared to the control (Table 2), suggesting that the Cd-alone and the co-exposure with biochar and α-Fe2O3 NPs resulted in the certain levels of abiotic stresses to muskmelon plants. The size of the fruits across all the treatments was also consistent with the physiological results (Fig. 2A). In addition, no significant difference in the other growth indicators of all co-treatments as compared to the Cd-alone treatment, except for the fruit transverse diameter in the BFNC treatments, which was not only significantly smaller than that of the control, but also notably smaller than the ones in the Cd-alone treatment, indicating that the co-treatments with both biochar and α-Fe2O3 NPs might have an adverse effect on the Cd treated muskmelon to some extent.

Fig. 3
figure 3

The growth record of muskmelon plants in the whole life cycle. Dynamic changes of plant height (A), leaf length (B), leaf width (C) and hypocotyl thickness (D) of muskmelon plants. The data show the mean ± SD of three replicates. Different lowercase letters indicate significant differences at p < 0.05

Table 2 Size and quality of muskmelon fruits across all the treatments

Physiological responses of muskmelon plants

Chlorophyll content

During the whole life cycle, the chlorophyll content in the muskmelon leaves of the Cd alone treatment was significantly lower than that of the control (Fig. 4A). Consistently, R Mera et al. (2016) also demonstrated that the addition of Cd resulted in the chlorophyll reduction in the microalga Chlamydomonas moewusii [44]. Previous studies also demonstrated that Cd [2]+ ions can replace Mg2+ in chlorophyll (Chl) to form [Cd]-Chl [45] and subsequently compromised photosynthesis due to the Cd-induced reduction of the chlorophyll content, which was unstable [42]. The chlorosis of leaves in muskmelon plants treated with Cd alone was also evident (Additional file 1: Figure S5). The decrease in the chlorophyll content resulted in lowering the photosynthetic efficiency, and consequently inhibiting the muskmelon growth, which aligns with the results of physiological parameters throughout the whole life cycle. The chlorophyll content in the treatments with α-Fe2O3 NPs, 1%BC and 5%BC were no difference with the control at the eleventh week while the chlorophyll content in the Cd-alone treatment was still lower than these treatments, indicating that these three treatments had a positive effect on alleviating the decrease in chlorophyll content caused by cadmium stress.

Fig. 4
figure 4

Physiological responses of muskmelon as affected by Cd, biochar and α-Fe2O3 NPs during the whole life cycle. Dynamic changes of chlorophyll (A), anthocyanin (B), soluble protein (C) and MDA content (D) in muskmelon plants. The data show the mean ± SD of three replicates. Different lowercase letters indicate significant differences at p < 0.05

Soluble protein content

Soluble protein is an important osmotic adjustment substance and nutrient substance in plants and contains various enzymes involved in metabolisms [46, 47]. It is often used as a phyisological indicator of plant resistance as the soluble protein synthesis is related to the resistance mechanisms of plants under abiotic stresses [33]. At the fourth and eleventh week, the soluble protein content in muskmelon leaves in the Cd-alone treatment was significantly lower by 16.5 and 50.5% as compared to the control. At the fourth week, the soluble protein content in the BC treatments and FeNPs treatment recovered to a comparable level with the control. At the eleventh week, the soluble protein content of muskmelon leaves in the single or co-treatments were restored to the control level.

Antioxidant systems

Abiotic stresses in plants can generate excessive amounts of reactive oxygen species (ROS) accumulation, which can result in oxidative stresses, including lipid peroxidation, membrane leakage, enzyme inactivation, and DNA damage [48]. CAT, POD, SOD, APX and anthocyanins are primary antioxidant substances that are in charge of breaking down ROS [49, 50].

Excessive amounts of ROS in cells can lead to membrane lipid peroxidation and membrane leakage, resulting in a significant increase in MDA content [51]. However, in the present experiment, high dose of Cd did not significantly induce the MDA content (Fig. 4D).

At the 11th week, the SOD and APX activities in leaves treated with Cd-alone were significantly increased as compared to the control; while the POD activity at the sixth week and the CAT activity at the 4th and 11th weeks were significantly decreased (Fig. 5). Previous studies reported that a dose-dependent pattern of the antioxidant enzyme activity was evident; the lower the Cd doses were, the higher the activities of antioxidant enzymes were found in plants [52]. In a relatively high Cd dose scenario, the various responses of the activities of four antioxidant enzymes indicate inconsistent sensibility to the Cd-induced phytotoxicity [53]. The reduction of antioxidant enzyme activity will inevitably lead to an increase in intracellular reactive oxygen species, which can also be reflected by the increase in anthocyanin content (Fig. 4B).

Fig. 5
figure 5

Antioxidant enzyme activities of muskmelon plants as affected by Cd, biochar and α-Fe2O3 NPs during the whole life cycle. Dynamic changes of SOD (A), POD (B), CAT (C) and APX activities (D) in muskmelon plants. The data show the mean ± SD of three replicates. Different lowercase letters indicate significant differences at p < 0.05

At the fourth and eighth weeks, the anthocyanin content (Fig. 4B) in Cd-alone treated leaves was significantly increased by 30.03 and 39.36% as compared to the control, respectively. In the FeNPs treatment, the anthocyanin content was decreased to the control level. Similarly, 5% biochar treatment at week 4 and 1% biochar treatment at week 8 also produced the same effect. The SOD and APX (Fig. 5A and D) activities in the Cd-alone treated leaves increased significantly, but were only found at the eleventh week; while their activities in FeNPs and BC treatments were recovered to the control level. At the eleventh week, the SOD activity in leaves in the BFNC treatments was significantly lower than that in the Cd-alone treatment, but still significantly higher than the control. The APX activity in FeNPs and BC treatments were recovered to the control level at the eleventh week.

At the sixth week, the activity of POD in the Cd-alone treatment was significantly reduced by 50% compared to the control. With the 5% BFNC treatment as an exception, whose POD activity was significantly higher than the Cd-alone treatment (Fig. 5B), there was no significant difference in the other remediation treatment with Cd alone treatment. The POD activity in the Cd-alone treatment also decreased insignificantly compared to the control at the fourth and eleventh weeks. Similarly, at the fourth and eleventh weeks, the CAT activity in the Cd-alone treatment was significantly lower than that of the control. At the eleventh week, in the co-treatments with 1% and 3% biochar, the CAT activities were significantly elevated when comparing to the Cd-alone treatment.

In the experiment, the responses of the four antioxidant enzymes to Cd stress were not consistent, which may be due to differences in their sensitivity to Cd [54]. In addition, most of the oxidative stress caused by Cd stress occur in the seedling and mature stages of muskmelon, which may be attributed to the plant's own regulatory effects [55].

Principal component analysis

Principal component analysis (PCA) is a method of dimensionality reduction analysis, referring to extracting one or two main components from a large number of indicators [56]. Through principal component analysis, we can better understand the differences between samples and the correlation between indicators.

In the experiment, we extracted two main components, covering more than 50% of all the changes across all the treatments and parameters. In the fourth and sixth weeks, Cd content was negatively correlated with leaf length, leaf width, plant height and hypocotyl thickness, indicating that Cd stress inhibited the growth of muskmelon plants, which was consistent with our previous results. Cd content and chlorophyll content were negatively correlated throughout the growth cycle, indicating that Cd inhibited the synthesis of chlorophyll in muskmelon leaves. Regarding the antioxidant systems, in general, Cd stress was positively correlated with the activities of SOD and APX during most of the growth period of muskmelon plants, but negatively correlated with the activities of CAT and POD.

From the distribution of the treatment groups in the Fig. 6, Cd-alone treatment had a significant impact on the growth of plants, and all Cd-containing treatments had a greater impact on the growth of muskmelon. From the week 4 to the week 11, the groups with single or co-treatments gradually moved away from the Cd treatment group, indicating that with the growth of muskmelon plants, both single or co-treatments alleviated the Cd stress to a certain extent. In most of the time (Fig. 5A, B and C), single and co-treatments were in different quadrants and were far apart, indicating that in contrast, single treatments had a better repair effect on Cd stress, which was consistent with our previous analysis.

Fig. 6
figure 6

Principal component analysis of all physiological parameters in muskmelon plants. Principal component analysis was conducted on the physiological indicators of muskmelon plants in the fourth (A), sixth (B), eighth (C), and eleventh (D) weeks. The horizontal axis is the extracted first and second principal components, respectively

Transcriptome analysis

RNA-Seq technique was used to reveal the transcriptomic response of muskmelon leaves as a function of Cd, biochar and α-Fe2O3. In comparison with the control, 394 genes were significantly different in the Cd-alone treatment. Among them, 301 genes were up-regulated by more than twofold. The results of KEGG enrichment (Fig. 7) show that differentially expressed genes are related to multiple signal transduction pathways and biochemical metabolic processes. Among them, in the Cd-alone treatment, genes associated with photosynthesis were significantly down-regulated; while genes related to anthocyanin synthesis and glutathione metabolisms, both of which can scavenge ROS in plants [33], were significantly up-regulated. Genes related to sulfur metabolism were significantly down-regulated. Cysteine, an important hub of sulfur metabolism, is a substrate of glutathione and is involved in maintaining the level of antioxidant glutathione [57, 58]. In addition, genes related to the mitogen-activated protein kinase (MAPK) signaling pathway were significantly up-regulated. Previous study also demonstrated that the MAPK pathway could respond to the Cd-induced abiotic stress [59]. Heavy metal stresses can alter the patterns of gene regulation in plants [60]. MAPK is one of the main pathways that plants transduce extracellular stimuli and other signals to cause intracellular responses [61]. Given the evidence on the significant up-regulation of the gluthation metabolism associated genes, it is reasonable to speculate that exposure to Cd induces strong oxidative stress in the muskmelon at the early stage. Through MAPK or other signal transduction pathways, the expression of photosynthesis-related genes was down-regulated, and the synthesis and metabolism of antioxidant-related substances were triggred/activated to break down ROS in plants. However, the expression of genes related to sulfur metabolism was significantly down-regulated, and thus, the gluthatione biosynthesis as the down-stream metabolism pathways of sulfur assimilation might be affected as well.

Fig. 7
figure 7

Enrichment of differentially expressed genes in muskmelon plants. The enrichment results of different genes in Cd (A), FeNPs (B), 1%BC (C), 3%BC (D), 5%BC (E), 1%BFNC (F), 3%BFNC (G) and 5%BFNC (H) treated muskmelon plants as compared to the control

In the remediation treatments, up-regulation of genes related to anthocyanin synthesis was found as compared to the control, indicating that even if biochar or Fe is added to the Cd pollution, thereis still a certain amount of stress on the muskmelon plants compared with the control. In the 1%BC treatment (Fig. 7C), the ABC transporter family related genes were found to be significantly up-regulated as compared to the control. AtPDR8 (member of the pleiotropic drug resistance) in the ABC transporter family can act as an efflux pump for Cd2+ ions or Cd conjugates in Arabidopsis [62], implying that these two remediation treatments might indirectly alter the Cd distribution by effluxing Cd out of the plant body and subsequently alleviate the Cd-induced phytotoxicity. Similar KEGG enrichments in anthocyanin biosynthesis, glutathione metabolism and MAPK signal transduction pathways were also found in the treatments with α-Fe2O3 NPs and 3%BC. Therefore, it can be speculated that after the addition of the remediation agent, in addition to the aforementioned regulation of the muskmelon plant itself, it may also activate metal ion transporters such as the ABC transporter family, thereby expelling Cd ions to reduce the toxicity of Cd to cells. The differential expression of a large number of genes appeared in the 1%BFNC and 5%BFNC treatments, which was significantly more than other treatments. Through enrichment analysis of differential genes, we found that there were significant differential expressions of genes related to ABC transporter, anthocyanin biosynthesis, glutathione metabolism and MAPK signal transduction pathways, which could alleviate the Cd stress. It should be noted that the 3% BFNC treatment has fewer differential genes as compared to the control, implying 50 mg/kg α-Fe2O3 NPs and 3% biochar was a relatively good composite ratio. But combined with the previous growth and physiological indicators, single treatments had a better repair effect than co-treatments on Cd stress.

Conclusions

Our research found that Cd pollution reduced the absorption of iron elements by muskmelon (Cucumis melo) plants (Additional file 1: Fig. S4), and also caused the growth of muskmelon plants to be inhibited (Fig. 3) and the size and quality of fruits to be reduced (Table 2). In addition, the high concentration of Cd pollution will also lead to a series of physiological and biochemical changes in muskmelon plants, including a decrease in chlorophyll and soluble protein content, an increase in anthocyanin content, and a decrease in POD and CAT activity. The addition of α-Fe2O3 NPs or biochar alone significantly reduced the absorption of Cd by plants, but the effect of co-treatments was not as good as that of the single treatments. The results of gene differential expression analysis showed that the expression of Cd transport related proteins in the repair treatments did increase significantly, which confirmed that the repair treatments had a certain positive effect on alleviating the Cd poisoning of plants. In addition, all treatments could significantly reduce the Cd content in fruit, and single treatments could alleviate the Cd toxicity according to various indicators. It should be noted that the compound treatment in our research does not seem to produce better results, but 50 mg/kg α-Fe2O3 NPs and 3% biochar was a relatively good composite ratio. Our research shows that biochar and α-Fe2O3 NPs can significantly reduce the effective Cd content in the soil and the Cd content in the fruit, which provides some ideas for the safe cultivation of muskmelon crops in Cd -contaminated soil. Our research provides some considerations for agricultural production and environmental protection under Cd pollution.

Availability of data and materials

All data associated with this study are present in the paper or additional materials.

References

  1. Ali A, Guo D, Jeyasundar P, Li Y, Zhang Z. Application of wood biochar in polluted soils stabilized the toxic metals and enhanced wheat (Triticum aestivum) growth and soil enzymatic activity. Ecotoxicology and Environmental Safety. 2019;184:109635.

    Article  CAS  PubMed  Google Scholar 

  2. Wu T, Liu C, Kong B, Sun J, Gong Y, Liu K, Xie J, Pei A, Cui Y. Amidoxime-functionalized macroporous carbon self-refreshed electrode materials for rapid and high-capacity removal of heavy metal from water. ACS Cent Sci. 2019;5(4):719–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Pa ND, Liu C, Yu H, Li F. A paddy field study of arsenic and cadmium pollution control by using iron-modified biochar and silica sol together. Environ Sci Pollut Res. 2019;26(24):24979–87.

    Article  Google Scholar 

  4. Qayyum MF, Rehman RA, Liaqat S, Ikram M, Hussain Q. Cadmium immobilization in the soil and accumulation by spinach ( Spinacia oleracea ) depend on biochar types under controlled and field conditions. Arab J Geoences. 2019;12(15):493.

    Article  Google Scholar 

  5. Ali S, Rizwan M, Noureen S, Anwar S, Ali B, et al. Combined use of biochar and zinc oxide nanoparticle foliar spray improved the plant growth and decreased the cadmium accumulation in rice (Oryza sativa L.) plant. Environm Sci Pollut Res. 2019;26:11288–99.

    Article  CAS  Google Scholar 

  6. Khan MA, Islam MR, Panaullah GM, Duxbury JM, Loeppert J. Accumulation of arsenic in soil and rice under wetland condition in Bangladesh. Plant Soil. 2010;333:263–74.

    Article  CAS  Google Scholar 

  7. Rizwan M, Ali S, Rehman M, Adrees M, Arshad M, et al. Alleviation of cadmium accumulation in maize (Zea mays L.) by foliar spray of zinc oxide nanoparticles and biochar to contaminated soil. Environ Pollut. 2019;248:358–67.

    Article  CAS  PubMed  Google Scholar 

  8. Ogawa T, Kobayashi E, Okubo Y, Suwazono Y, Kido T, Nogawa K. Relationship among prevalence of patients with Itai-itai disease, prevalence of abnormal urinary findings, and cadmium concentrations in rice of individual hamlets in the Jinzu River basin, Toyama prefecture of Japan. Int J Environ Health Res. 2004;14(4):243–52.

    Article  CAS  PubMed  Google Scholar 

  9. Kamran M, Malik Z, Parveen A, Zong Y, Ali M. Biochar alleviates Cd phytotoxicity by minimizing bioavailability and oxidative stress in pak choi (Brassica chinensis L.) cultivated in Cd-polluted soil. J Environ Manag. 2019;250:109500.

    Article  CAS  Google Scholar 

  10. Liu X, Zeng Z, Tie B, Chen Q, Chen Z, Zheng R. Effects of biochar and lime application on soluble Cd, Pb, As release and non-point loads of rice agroecosystem by in situ field experiment, Central Hunan Province Mining Area. Hamdard Medicus. 2013;32(1):28–35.

    Google Scholar 

  11. Abdelhamid HN, Mathew AP. Cellulose-zeolitic imidazolate frameworks (CelloZIFs) for multifunctional environmental remediation: Adsorption and catalytic degradation. Chem Eng J. 2021;426:131733.

    Article  Google Scholar 

  12. Fijoł N, Abdelhamid HN, Pillai B, et al. 3D-printed monolithic biofifilters based on a polylactic acid (PLA)-hydroxyapatite (HAp) composite for heavy metal removal from an aqueous medium[J]. RSC Adv. 2021;11(51):32408–18.

    Article  Google Scholar 

  13. Nzediegwu C, Prasher S, Elsayed E, Dhiman J, Mawof A, Patel R. Effect of biochar on heavy metal accumulation in potatoes from wastewater irrigation. J Environ Manag. 2018;232(15):153–64.

    Google Scholar 

  14. Jin HP, Choppala GK, Bolan NS, Chung JW, Chuasavathi T. Biochar reduces the bioavailability and phytotoxicity of heavy metals. Plant Soil. 2011;348(1):439–51.

    Google Scholar 

  15. Younis U, Malik SA, Rizwan M, Qayyum MF, Ok YS, Shah M, Rehman RA, Ahmad N. Biochar enhances the cadmium tolerance in spinach (Spinacia oleracea) through modification of Cd uptake and physiological and biochemical attributes. Environ Sci Pollut Res. 2016;23(21):1–10.

    Article  Google Scholar 

  16. Zhang ZY, Meng J, Dang S, Chen WF. Effect of biochar on relieving cadmium stress and reducing accumulation in Super japonica rice. J Integr Agric. 2014;13(3):547–53.

    Article  CAS  Google Scholar 

  17. Woldetsadik D, Drechsel P, Keraita B, Marschner B, Itanna F, Gebrekidan H. Effects of biochar and alkaline amendments on cadmium immobilization, selected nutrient and cadmium concentrations of lettuce (Lactuca sativa) in two contrasting soils. Springerplus. 2016;5:397.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Rajendran M, Shi L, Wu C, Li W, An W, Liu Z, Xue S. Effect of sulfur and sulfur-iron modified biochar on cadmium availability and transfer in the soil–rice system. Chemosphere. 2019;222:314–22.

    Article  CAS  PubMed  Google Scholar 

  19. Liu XM, Zhang FD, Feng ZB, Zhang SQ, Xu-Sheng HE, Wang RF, Wang YJ. Effects of Nano-ferric oxide on the growth and nutrients absorption of peanut. Plant Nutr Fertil. 2005;11(4):551–5.

    Google Scholar 

  20. Zhang R, Chen W. Graphene-non-noble metal hybrid nanomaterials as advanced electrocatalysts. 2014.

  21. Palmqvist N, Seisenbaeva GA, Svedlindh P, Kessler VG. Maghemite nanoparticles acts as nanozymes, improving growth and abiotic stress tolerance in Brassica napus. Nanoscale Res Lett. 2017;12(1):631.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Zhu D, Jian Z, Song J, Wang H, Zheng Y, Shen Y, Xie A. Efficient one-pot synthesis of hierarchical flower-like α-Fe2O3 hollow spheres with excellent adsorption performance for water treatment. Appl Surf Sci. 2013;284(1):855–61.

    Article  CAS  Google Scholar 

  23. Maia L, Santos MS, Andrade TG, Hott R, Rodrigues JL. Removal of mercury(II) from contaminated water by gold-functionalised Fe3O4 magnetic nanoparticles. Environ Technol. 2018;4:1–40.

    Google Scholar 

  24. Li J, Wan F, Guo W, Huang J, Wang Y. Influence of α- and γ-Fe2O3 Nanoparticles on Watermelon (Citrullus lanatus) Physiology and Fruit Quality. Water Air Soil Pollut. 2020;231(4):143.

    Article  CAS  Google Scholar 

  25. Zou Z, Wang Y, Huang J, Lei Z, Wang F, Dai Z, Yi L, Li J. A study on the mixture repairing effect of biochar and nano iron oxide on toxicity of Cd toward muskmelon. Environ Pollut. 2020;266:1.

    Article  Google Scholar 

  26. Alidoust D, Isoda A. Phytotoxicity assessment of γ-Fe2O3 nanoparticles on root elongation and growth of rice plant. Environ Earth Sci. 2014;71(12):5173–82.

    Article  CAS  Google Scholar 

  27. Fang G, Grumet R. Transformation in Muskmelon (Cucumis Melo L.). Plant Protoplasts Genetic Eng IV. 1993;23:209–14.

    Article  Google Scholar 

  28. Geng X, Li K, Wang F, Pan T, Diao Y, Lin Y. Extraction, purification and pharmacological activity of melon pedicle active ingredients. Asian J Chem. 2014;26(22):7569–74.

    Article  Google Scholar 

  29. Lopez-Martinez LX, Oliart-Ros RM, Valerio-Alfaro G, Lee CH, Parkin KL, Garcia HS. Antioxidant activity, phenolic compounds and anthocyanins content of eighteen strains of Mexican maize. Food Sci Technol. 2009;42(6):1187–92.

    CAS  Google Scholar 

  30. Wang Y, Wang S, Xu M, Xiao L, Dai Z, Li J. The impacts of γ-Fe2O3 and Fe3O4 nanoparticles on the physiology and fruit quality of muskmelon (Cucumis melo) plants. Environ Pollut. 2019;249:1011–8.

    Article  CAS  PubMed  Google Scholar 

  31. Ma D. Studies on salt tolerance of transgenic sweetpotato which harbors two genes expressing CuZn superoxide dismutase and ascorbate peroxidase with the stress-inducible SWPA2 promoter. Plant Gene Trait. 2011;3:6–12.

    CAS  Google Scholar 

  32. Haq AU, Bates TE, Soon YK. Comparison of extractants for plant-available zinc, cadmium, nickel, and copper in contaminated soils. Soil Sci Soc Am J. 1980;44(4):772–7.

    Article  CAS  Google Scholar 

  33. Xiao L, Guo H, Wang S, Li J, Wang Y, Xing B. Carbon dots alleviate the toxicity of cadmium ions (Cd2+) toward wheat seedlings. Environ Nano. 2019;6:1493.

    Article  CAS  Google Scholar 

  34. Ding W, Zhu Q, Zeng X, Dan WU. Biochars from different pyrolytic temperature amending lead and cadmium contaminated soil. Sci Technol Rev. 2011;29:22–5.

    CAS  Google Scholar 

  35. Fan S, Li H, Wang Y, et al. Cadmium removal from aqueous solution by biochar obtained by co-pyrolysis of sewage sludge with tea waste. Res Chem Intermed. 2018;44(1):135–54.

    Article  CAS  Google Scholar 

  36. Muehe EM, Obst M, Hitchcock AP, et al. Fate of Cd during microbial Fe(III) mineral reduction by a novel and Cd-tolerant Geobacter species. Environ Sci Technol. 2013;47(24):14099–109.

    Article  CAS  PubMed  Google Scholar 

  37. Colmer TD. Long-distance transport of gases in plants: a perspective on internal aeration and radial oxygen loss from roots. Plant, Cell Environ. 2010;26(1):17–36.

    Article  Google Scholar 

  38. Chen RF, Shen RF, Gu P, Dong XY, Du CW, Ma JF. Response of rice (Oryza sativa) with root surface iron plaque under aluminium stress. Ann Bot. 2006;98(2):389–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Wu P, Cui P, Marcelo E, et al. Interactive effects of rice straw biochar and γ-Al2O3 on immobilization of Zn. J Hazard Mater. 2019;373:250–7.

    Article  CAS  PubMed  Google Scholar 

  40. Krotz RM, Evangelou BP, Wagner GJ. Relationships between cadmium, zinc, Cd-peptide, and organic acid in tobacco suspension cells. Plant Physiol. 1989;91:780–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Salt DE, Rauser WE. MgATP-dependent transport of phytochelatins across the tonoplast of oat roots. Plant Physiol. 1995;107:1293–301.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Andresen E, Küpper H. Cadmium toxicity in plants. Met Ions Life Sci. 2013;11:395–413.

    Article  CAS  PubMed  Google Scholar 

  43. Nazar R, Iqbal N, Masood A, Khan M, Khan NA. Cadmium toxicity in plants and role of mineral nutrients in its alleviation. Am J Plant Sci. 2012;3(10):1476–89.

    Article  Google Scholar 

  44. Mera R, Torres E, Abalde J. Influence of sulphate on the reduction of cadmium toxicity in the microalga Chlamydomonas moewusii. Ecotoxicol Environ Saf. 2016;128:236–45.

    Article  CAS  PubMed  Google Scholar 

  45. Hendrik K, Frithjof K, Martin S. Environmental relevance of heavy metal-substituted chlorophylls using the example of water plants. J Exp Bot. 1996;2:259–66.

    Google Scholar 

  46. Wang K, Guo JJ, Wang DQ, Liu S, Ni PP. Responses of roots and needles of Pinus sylvestris var. mongolica and Pinus tabuliformis to spring drought stress. Chin J Ecol. 2015;8:1477–81.

    Google Scholar 

  47. Zhao MZ, Liu ZL, Chen W, Cai SY. Responses of Trifolium repens cv. rivendel seeds to cadmium stress in terms of electrolyte leakage and soluble protein content changes. Adv Mater Res. 2014;1:955–9.

    Article  Google Scholar 

  48. Halliwell B, Gutteridge J. Free radicals in biology and medicine. Biochem Soc Trans. 1985;13(5):976.

    Article  Google Scholar 

  49. Anjum SA, Tanveer M, Hussain S, Shahzad B, Ashraf U, Fahad S, Hassan W, Jan S, Khan I, Saleem MF. Osmoregulation and antioxidant production in maize under combined cadmium and arsenic stress. Environ Sci Pollut Res. 2016;23(12):11864–75.

    Article  CAS  Google Scholar 

  50. Xu Z, Rothstein SJ. ROS-Induced anthocyanin production provides feedback protection by scavenging ROS and maintaining photosynthetic capacity in Arabidopsis. Plant Signal Behav. 2018;13(3):e1451708–e1451708.

    Article  PubMed  PubMed Central  Google Scholar 

  51. García A. Vermicompost humic acids modulate the accumulation and metabolism of ROS in rice plants. J Plant Physiol. 2016;192(15):56–63.

    Article  PubMed  Google Scholar 

  52. Sandalio LM, Dalurzo HC, Gómez M, et al. Cadmium-induced changes in the growth and oxidative metabolism of pea plants. J Exp Bot. 2001;364:2115–26.

    Article  Google Scholar 

  53. Xue H, Zhang H, Chang H, Jiao Y, Li H. Effect of cadmium on antioxidase system in process of maize germination. Asian J Chem. 2015;27(11):4022–8.

    Article  CAS  Google Scholar 

  54. Batool M, Abdullah S, Abbas K. Antioxidant enzymes activity during acute toxicity of chromium and cadmium to Channa marulius and Wallago attu. Pak J Agric Sci. 2014;51(4):1017–23.

    Google Scholar 

  55. Wang Y, Hu J, Dai Z, Li J, Huang J. In vitro assessment of physiological changes of watermelon (Citrullus lanatus) upon iron oxide nanoparticles exposure. Plant Physiol Biochem. 2016;108:353–60.

    Article  CAS  PubMed  Google Scholar 

  56. Amini A, Wainwright M. High-dimensional analysis of semidefinite relaxations for sparse principal components. Ann Statist. 2008;37(5B):2877–921.

    Google Scholar 

  57. Yanmei C, Yuanqing C, Yaying L, Qingqi L, Jun B. Survival strategies of the plant-associated bacterium Enterobacter sp. strain EG16 under cadmium stress. Appl Environ Microbiol. 2016;82(6):1734–44.

    Article  Google Scholar 

  58. Capaldi FR, Gratao PL, Reis AR, Lima LW, Azevedo RA. Sulfur metabolism and stress defense responses in plants. Trop Plant Biol. 2015;8(3–4):60–73.

    Article  CAS  Google Scholar 

  59. Jonak C, Nakagami H, Hirt H. Heavy metal stress. Activation of distinct mitogen-activated protein kinase pathways by copper and cadmium. Plant Physiol. 2004;136(2):3276.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Lin CW, Chang HB, Huang HJ. Zinc induces mitogen-activated protein kinase activation mediated by reactive oxygen species in rice roots. Plant Physiol Biochem. 2005;43(10–11):963–8.

    Article  CAS  PubMed  Google Scholar 

  61. Zhao FY, Han MM, Zhang SY, Ren J, Hu F, Wang X. MAPKs as a cross point in H2O2 and auxin signaling under combined cadmium and zinc stress in rice roots. Russ J Plant Physiol. 2014;61:608–18.

    Article  CAS  Google Scholar 

  62. Kim DY, Bovet L, Maeshima M, Martinoia E, Lee Y. The ABC transporter AtPDR8 is a cadmium extrusion pump conferring heavy metal resistance. Plant J. 2010;50(2):207–18.

    Article  Google Scholar 

Download references

Acknowledgements

We thank the Core Facility and Technical Support at the Hubei Academy of Agricultural Science for technical support.

Funding

This work was supported by China Agriculture Research System of MOF and MARA [Grant No. CARS-25]; Key Research and Development Program of Hubei Province, China [Grant No. 2020BBA037]; the Fundamental Research Funds for the Central Universities [Grant No. 2019IB005, 2021IA005].

Author information

Authors and Affiliations

Authors

Contributions

YW, ZZ and JL designed the experiments, performed the experiments, curated the data, analyzed the data, and wrote the manuscript; XS, FW, YZ, ZL perform some of the experiments and edited the manuscript; LY and ZD managed the project and edited the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Junli Li.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

: Additional figures and tables.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, Y., Zou, Z., Su, X. et al. Physiological of biochar and α-Fe2O3 nanoparticles as amendments of Cd accumulation and toxicity toward muskmelon grown in pots. J Nanobiotechnol 19, 442 (2021). https://doi.org/10.1186/s12951-021-01187-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12951-021-01187-7

Keywords