Concentrations of heavy metals continue to increase in soil environments as a result of both anthropogenic activities and natural processes. Cadmium (Cd) is one of the most toxic heavy metals and poses health risks to both humans and the ecosystem. Herein, we explore the impacts of Cd on a soil-plant system composed of oilseed rape (Brassica napus and Brassica juncea) and bacteria. The two species of oilseed rape displayed a similar variation trend under Cd treatment. Cd accumulation within plant tissues enhanced with increasing concentrations of Cd in soils, and Cd treatment decreased chlorophyll content and suppressed plant growth. Meanwhile, Cd stress induced the changes of antioxidative enzymes activities including elevating SOD and POD activities and reducing CAT activity. The impact of Cd on the bacterial communities of soils was greater than bacterial communities of plants (phyllosphere and endophyte). The α-diversity of bacterial community in soils declined significantly under higher Cd concentration (30 mg/kg). In addition, soil bacterial communities composition and structure were altered in the presence of higher Cd concentration. Meanwhile, the bacterial community of bulk soil was significantly correlated with Cd, while the variation of rhizosphere soil bacterial community was markedly correlated with Cd and other environmental factors of both soils and plants. These results suggested that Cd could affect both the growth of plants and the indigenous bacterial community in soil-plant system, which might further change ecosystem functions in soils.
Research Article
The Growth of Plants and Indigenous Bacterial Community Were Significantly Affected by Cadmium Contamination in Soil-Plant System
https://doi.org/10.21203/rs.3.rs-181936/v1
This work is licensed under a CC BY 4.0 License
Journal Publication
published 10 Jul, 2021
You are reading this latest preprint version
Cd
Bacterial community
Oilseed rape
Soil
Phyllosphere
Endophyte
Heavy metals (HMs) in soils have become a serious environmental issue due to their poisonousness and bioaccumulation within the food chain. They can deteriorate soil quality, reduce food safety, and threaten human health (Li et al., 2014). Moreover, these metals are not degraded by chemical or biological methods and are persistent in soil (Cheraghi-Aliakbari et al., 2020). Cadmium (Cd), one of the most toxic heavy metals, and nonessential for humans and animals, is widely present in soil due to anthropogenic activities such as smelting, mining and battery disposal (DalCorso et al., 2019; Mitra et al., 2018). Cd accumulation in food could cause serious health problems in humans and animals (Khan et al., 2017). Phytoremediation is an efficient and environmental-friendly method to remove Cd from soil and could be used to remediate Cd pollution soil (Liu et al., 2020)..
Under Cd stress, plants are damaged by photosynthesis inhibition and disruption of nutrition absorption (Li et al., 2018), affecting the growth of sensitive plants (Wang et al., 2020a). However, Cd hyperaccumulators have higher tolerance to Cd, these include many types of Brassica species that can be grown under heavy metal stress in soil, with the whole plant being used for biofuel production (Rizwan et al., 2018). Brassica napus and Brassica juncea are widely grown around the world and have been investigated extensively for the remediation of Cd (Goswami and Das, 2015; Rossi, 2002).
Soil is the main medium for terrestrial ecosystems, supporting productive activities, regulating nutrient flow, and maintaining ecosystem health with microbes playing a significant role in these processes (Xia et al., 2018). However, many studies have found that microbial diversity, abundance, and composition in soils are strongly affected by Cd contamination (Hou et al., 2018; Wood et al., 2016); As an essential part of ecosystems, microbial communities play a important role in global biogeochemical cycle (Beattie et al., 2018). Moreover, microorganisms are essential in biogeochemical cycling of HMs (Jing and Kjellerup, 2018), with their responses to pollution having profound ecological effects and can serve as biological indicators of heavy metals such as Cd toxicity. Microorganisms are an important component of phytoremediation technology (Zhang et al., 2012).
Phyllosphere (Jia et al., 2018) and endophytic bacteria (Wang et al., 2020b) are also beneficial in phytoremediation. However, less attention has been paid to the effect of Cd on plant microbial community (phyllospheric and endophytic bacterial communities) in the soil-plant ecosystem.
In this study, we chose two species of oilseed rapes (B. napus and B. juncea) to investigate the effect of Cd on plants and compared their Cd-tolerance. A 16S rRNA gene amplicon Illumina Miseq approach was performed to examine effects of Cd contamination on bacterial communities in the soil-plant ecosystem. Our results provide a detailed understanding of the effects of Cd on plants and soil-plant system microbial communities and help improve phytoremediation systems.
Greenhouse experiments
The test soil was agricultural topsoil (0-20 cm) taken from a suburb of Hunan province (27°54′15″N, 112°55′06′'E ) and the main properties of the soil prior to any treatment were as follows: pH 5.56, total organic carbon (TOC) 1.62%, total nitrogen (TN) 1803.74 mg/kg, total phosphorus (TP) 921.04 mg/kg, available phosphorus (AP) 78.14 mg/kg, available potassium (AK) 135.41 mg/kg, ammoniacal nitrogen (NH3_N ) 17.65 mg/kg and nitrate nitrogen (NO3_N) 19.44 mg/kg and a background Cd concentration of 0 mg/kg.
The seeds of B. napus and B. juncea are Zhong-shuang 11 and Purple-leaf B. juncea, respectively. The plants were grown in a greenhouse located at the Hunan University of Science and Technology (27°54′15″N, 112°55′06′'E, Hunan, China) between October 11, 2018 and Noveber 30, 2018. Soils were disposed by air-dried, ground and sieved by a 5mm mesh. Then Cd aqueous solution (CdCl2·2.5H2O) was carefully added into the soil to obtain three Cd concentrations (Control (CK), 0 mg/kg Cd; concentration 1 (C10), 10 mg/kg Cd; and concentration 2 (C30), 30 mg/kg) were applied. After mixing and being allowed to stabilize for three weeks, the soil (7 kg) was transferred into plastic pots (14 cm tall, 20 cm width, 49 cm height). Seeds were sown directly into pots and ten plants were maintained in each pot finally, they were planted in the greenhouse at 20°C ±10°C with suitable humidity. Every treatment was set to six replicates.
Sample collection
Soil samples were divided into two parts, one for measuring soil physico-chemical properties while the other for conducting molecular tests. Plant samples were obtained 50 days after planting, the aboveground and belowground parts were harvested dividually by sterile scissors. Microorganisms in the soil-plant ecosystems (including the rhizosphere, bulk soil, phyllospheric, and endophytic bacterial communities) were collected on the basis of Kong et al. (2018).
Measure of Cd content in soil and plant tissues
The samples were oven-dried, then ground and finally digested using HNO3 in a microwave (PyNN 140899, Peian, Beijing, China), and the contents of Cd were measured using flame atomic absorption spectrophotometer (Agilent 200 AA, Agilent Technology Co. LTD).
The Bioaccumulation Factor (BAF) and Translocation Factor (TF) of B.napus and B. juncea were calculated as follows:
Physiological index of plant tissues
Leaf samples were used for the determination of biochemical indices, including soluble sugar (Anthrone colorimetry method), superoxide dismutase ( SOD, hydroxylamine method), peroxidase (POD, colorimetric method ), and catalase (CAT, visible light method), which were used commercial kits (Nanjing Jiancheng Bioengineering Institute, China).The content of soluble protein and chlorophyll in leaf samples was measured by the Coomassie Brilliant Blue G250 staining method and acetone extraction method, respectively.
Physicochemical property of soils
Soil pH and the content of TN, TP, TOC, NH3_N, NO3_N, AP and AK were measured according to Kong, et al. (2018).
DNA extraction, PCR amplification, and sequencing
A total of 180 samples (divided equally into 5 parts including phyllosphere, leaf endophyte, root endophyte, rhizosphere soil and bulk soil samples) were sequenced by the following steps. DNA extraction, amplification and sequencing were carried by according to Kong, et al. (2018).The samples were sequenced using Illumina Miseq platform at BeiJing Fixgene Co., LTD.
Quantitative PCR (qPCR)
For determining abundances of bacteria, qPCR was performed by using the primer pair 799F (5'-AACMGGATTAGATACCCKG-3') and 1115R (5'-AGGGTTGCGCTCGTTG-3') together with a CFX Connect TM Real-Time PCR Detection System (BioRad). The qPCR reaction mixture which used MonAmpTM SYBR® qPCR Mix (Monad Biotech Co.Ltd) was done in a volume of 20µl consisting of 10 μL of qPCR Mix, 1 μL of each primer, 1 μL of template DNA and 7 μL of nuclease-free water. qPCR procedure was set as follows: 95°C for 30 s, 40 cycles of 95°C for 5 s, 57°C for 10 s, 72°C for 30 s. Every sample was amplified in triplicate.
Sequencing data process
Processing of the raw sequencing data were performed on a publically available Galaxy pipeline (http://mem.rcees.ac.cn:8080) at Research Center for Eco-Environmental Science, China Academy of Sciences. The procedures were as follows: raw reads were assigned to different samples using barcodes, followed by removal of primers sequences, forward and reverse sequences were combined by FLASH (Magoč and Salzberg, 2011), and sequences shorter than 200 bp were removed using Btrim (Kong, 2011). Extraction of FASTA data from FASTQ data, checking for and removal of chimeras, and assignment of sequences with 97% identity to the same operational taxonomic unit (OTU) were performed using the UPARSE algorithm (Edgar, 2013). The RDP Classifier database was used to Taxonomic assignment (Wang et al., 2007). In order to avoid the effects caused by different sequencing depth, the data were resampled randomly with the minimum number of sequences (30000). The following analysis used the resampled OTU table.
The raw reads generated in this study have been deposited in the NCBI sequence Read Archive (accession no. SRP283176).
Statistical analysis
The top 10 phyla and top 50 genera were selected in line with the results of phyllospheric, leaf and root endophytic, rhizosphere and bulk soil bacterial communities sample species annotations respectively. One-way ANOVA method was conducted to analyze the significance difference (P<0.05) between treatments. The relationship between environmental factors and α-diversities was investigated by pearson's correlation analysis. ANOVA method and Pearson's correlation test were performed using SPSS 21 software. The dissimilarity test was performed to evaluate the significance of clustering. Weighted principal coordinate analysis (PCoA) on the UniFrac matrix was applied to compare the different samples' bacterial community structure. Mantel test, canonical correspondence analysis (CCA) was used to show which environmental factors (EFs) significantly impact microbial community structure and CCA-based variation partitioning analysis (VPA) was used to determine the contributions of EFs to bacterial community.
Effect of Cd on physiological properties of oilseed rapes and Cd accumulation
Compared to the control (CK), Cd treatment suppressed plant growth (Fig. 1), specifically plant height, fresh weight, and total leaf area were significantly reduced with increasing Cd concentration (P<0.05). At CK and lower Cd concentration (10 mg/kg), the biomass of B. napus was significantly higher than B. juncea, but at higher Cd concentration (30 mg/kg) those tendencies were reversed. Pearson correlation analysis showed that Cd content in tissues was negatively correlated with plant height, weight, and leaf area in B. napus and B. juncea (P<0.01) (Additional file 1: Table S1 and Table S2).
Chlorophyll was significantly decreased under higher Cd concentration, with B. juncea being significantly less affected than B. napus (P<0.05) (Fig. 1D). Soluble sugar content was markedly decreased under Cd stress in both of B. napus and B. juncea, while the content of soluble protein was markedly increased under Cd treatment when compared with CK (P<0.05) (Fig. 1E , F). SOD and POD activity values showed similar trends, decreasing at first and then increasing. The SOD activity at 30 mg/kg was significantly higher than in other treatments (P<0.05). However, CAT activity was notably inhibited at the higher Cd concentration (P<0.05). SOD and POD activity were higher in B. napus than B. juncea, while CAT activity values were the reverse (Fig. 1G, H, I). Pearson correlation test showed that plant Cd content was negatively related with chlorophyll, soluble sugar, and CAT (P < 0.01) and positively related with soluble protein and SOD (P < 0.05) in B. napus and B. juncea (Additional file 1: Table S1 and Table S2).
The concentration of Cd in oilseed rape leaves and roots were significantly higher with increasing Cd levels, and Cd content in B. napus tissues was significantly higher than B. juncea (P<0.05) (Additional file 2: Figure S1). TF in the two oilseed rapes species significantly decreased with increasing Cd and was higher in B. napus than B. juncea (P<0.05) (Additional file 1: Table S3). BAF in leaves and roots of both B. napus and B. juncea was bigger than 1 and higher in 30 mg/kg compared to 10 mg/kg, especially in roots (Additional file 1: Table S3).
Effect of Cd on physicochemical properties of soils
In B. napus, the pH, TN, TP, and NO3_N were significantly decreased under the higher Cd concentration (P<0.05) and TOC first increased and then decreased in both rhizosphere and bulk soils (Fig. 2), meanwhile TOC, TN, NH3_N, and AP were higher, and pH lower, in rhizosphere than bulk soil. pH and TN were markedly reduced in the 30 mg/kg Cd treatment (P<0.05) and NH3_N was first increased and then decreased in both rhizosphere and bulk soils, while AK was significantly increased in bulk soil (P<0.05) in B. juncea under the higher Cd concentration (Fig. 2). Furthermore, TOC, TN, TP, NH3_N, NO3_N, and AK were higher, and pH was lower, in rhizosphere than bulk soil. Compared with B. napus, most soil nutrients were higher in B. juncea samples under the higher Cd treatment.
Pearson correlation analysis showed that for B. napus, the Cd content was negatively correlated with pH, TN, TP, and NO3_N content in both rhizosphere and bulk soils and with TOC of the rhizosphere (Additional file 1: Table S4 and Table S6). Cd was negatively correlated with pH, TN, and NH3_N content in rhizosphere and bulk soils and negatively correlated with AP of the rhizosphere, but positively correlated with AK of bulk soil in B. juncea (P < 0.05) (Additional file 1: Table S5 and Table S7).
Effect of Cd on bacterial numbers in soils
The results revealed that total bacterial numbers slightly increased with increasing Cd levels (Additional file 2: Fig. S2) and the Pearson analysis showed that Cd was positively correlation with bacterial numbers in bulk soil of B. napus (Additional file 1: Table S8). Meanwhile, no significant difference was observed between B. napus and B. juncea (P > 0.05).
Effect of Cd on the α-diversity of bacterial community
After removing low quality reads and chimaeras in 108 plant samples and 72 soil samples, a total of 13,352,813 high-quality 16s rRNA gene reads were obtained, which were clustered into 14359 phylotypes (OTUs) by grouping at a 97% identity threshold. The sequencing depths of all samples were appropriate for downstream analyses (Additional file 2: Fig. S3).
In the plant samples, Cd mainly affected root endophytic community of B. napus and phyllosphere community of B. juncea. Shannon index and richness of B. napus’s root endophytes and the richness and Chao1 of B. juncea’s phyllosphere decreased significantly at the higher Cd concentration (P<0.05) (Additional file 2: Fig. S4). Pearson correlation analysis demonstrated that plant physiological factors mainly correlated with α-diversity indexes of root endophytes in B. napus (Additional file 1: Table S10) and that the majority of plant physiological properties were cardinally correlated with OTU numbers (richness and Chao1) of phyllosphere in B. juncea (Additional file 1: Table S11). pH and TOC were significantly positive correlated with and Cd was significantly negative correlated with the α-diversity of root endophytic bacterial communities in B. napus (Additional file 1: Table S10).
However, high level of Cd significantly depressed the α-diversities in the soil bacterial communities. Inverse Simpson index and richness of rhizosphere and inverse Simpson indexes of bulk soil were markedly reduced in the 30 mg/kg Cd treatment in B. napus (P<0.05) (Fig. 3). Shannon and inverse Simpson indexes of rhizosphere and Shannon, inverse Simpson, richness and Chao1 of bulk soils in B. juncea were significantly decreased under the higher Cd treatment (P<0.05) (Fig. 3). Most α-diversity indexes between the two species of oilseed rapes had no significant differences (Fig. 3 and Additional file 2: Fig. S4).
Pearson tests showed that plant’s physiological factors mainly influenced α-diversity indexes of the rhizosphere in B. napus and bacterial diversity (Shannon and inverse Simpson indexes) of the rhizosphere in B. juncea (Additional file 1: Table S12 and Table S13). pH, TOC, TN, TP, and NO3_N were positively correlated with the α-diversity in rhizosphere of B. napus. Meanwhile, TN and NO3_N were positively correlated with bacterial diversity (Shannon and inverse Simpson indexes) of bulk soil in B. napus samples (Additional file 1: Table S12). In B. juncea, pH was positively, and AK negatively, correlated with α-diversity in bulk soil (Additional file 1:Table S13). However, Cd concentration showed a significant negative correlation with α-diversity of soil bacteria communities for both species of oilseed rapes (P<0.05) (Additional file 1: Table S12 and Table S13).
Effect of Cd on bacterial community composition and structure
Cd could affect the composition of bacterial communities in soil-plant ecosystem, particularly under higher levels of Cd stress (Additional file 1: Table S14, Table S15, Table S16, Table S17 and Table S18). Under 30 mg/kg Cd treatment, the bacterial relative abundance on the phylum level of Gemmatimonadetes, and Chloroflexi were significantly decreased in the B. juncea phyllosphere (P<0.05). Actinobacteria was significantly decreased in B. juncea’s root endophytic community and the B. napus rhizosphere (P<0.05). Gemmatimonadetes and Verrucomicrobia were significantly reduced in both rhizosphere and bulk bacterial community (P<0.05), but Proteobacteria and Bacteroidetes were significantly increased in the B. napus rhizosphere and Firmicutes was significantly increased in bulk soil of both of two rapeseed species.
At the genus level (Fig. 4), the relative abundance of some genera was altered under Cd treatment. In the B. napus phyllosphere samples, the relative abundances of Massililia sp., Rhodanobacter sp., and Rickettsia sp. were increased, and Buchera sp., Achromobacter sp., and Acinetobacter sp. were decreased under Cd treatment. While in B. juncea phyllosphere samples, Lysobacter sp., Stenotrophomonas sp., and Gibbsiella sp. were increased, while Gaiell sp., Telluria sp., and Herbaspirillum sp. were decreased under Cd treatment. In leaf endophyte samples, Brochothrix sp. and Acinetobacter sp. were increased in B. napus but decreased in B. juncea under Cd treatment. In root endophyte samples, Chryseobacterium sp. and Pantoea sp. were increased and Caulobacter Ideonella sp. and Herbaspirillum sp. were decreased in B. napus. Sphingomonas sp., Ralstonia sp., and Methylobacterium sp. were increased, and Rhizobium sp., Rhodanobacter sp., and Duganella sp. were decreased in B. juncea under Cd treatment.
In the rhizosphere, Niastella sp., Methylotenera sp., and Lystobacter sp. were increased and Arthrobacter sp., Gemmatimanas sp., and Haliangium sp. were decreased in B. napus under 30mg/kg Cd treatment. Massilia sp., Ralstonia sp., and Streptomyces sp. were increased and GP2 sp., Terriglobus sp., and Candidatus Solibacter sp. were decreased in B. juncea under 30mg/kg Cd treatment. In bulk soil, Sphingomonas sp., Rhodanobacter sp., and Roseateles sp. were increased and Arthrobacter sp., Gemmatimonas sp., and Terriglobus sp. were decreased in B. napus. Streptomyces sp., Pseudomocardia sp., and Blastococcus sp. were increased and Haliangium sp., Phenylobacterium ap., and Gemmatimonas sp. were decreased in B. juncea under Cd treatment.
The principal co-ordinates analysis (PCoA) (Fig. 5) and dissimilarity analysis (Additional file 1: Table S19 and Table S20) indicated that the bacterial community structures of both rhizosphere and bulk in both B. napus and B. juncea were significantly changed under higher Cd concentration compared to control (P < 0.05), but not significantly affected plant bacterial community structures.
Relationship between microbial community structure and environmental factors
The result of Mantel test showed that there are no significant association between most environment factors and phyllosphere or leaf endophyte bacterial communities (Additional file 1: Table S21 and Table S22). Biomass (height, weight, and leaf area), TOC and root Cd had significant association with root endophyte bacterial community in B. napus (Additional file 1: Table S21). Biomass, pH, NO3_N and Soil_Cd had significant association with rhizosphere soil and TN and Soil_Cd had significant association with bulk soil in B. napus (Table 1). Biomass, TN and Soil_Cd were significant correlated with rhizosphere soil bacterial communityin B. juncea and pH and Soil_Cd were significant correlated with bulk soils bacterial community in B. juncea (Additional file 1: Table S23).
The CCA model of root endophyte, rhizosphere and bulk soil bacterial community were significant (P<0.05, Additional file 2: Fig. S5A, Fig. S6A and Fig. 6A). The results of VPA indicated that biomass, pH, soil nutrients, and Root_Cd explained 14.3%, 5.4%, 37.4%, and 7.2% of variation in B. napus (Additional file 2: Fig. S5B) and 15.5%, 6.9%, 40.2%, and 4.2% of variation in B. juncea (Additional file 2: Fig. S5C) in root endophyte bacterial community, respectively.
CCA-based VPA indicated that biomass, pH, soil nutrients, and Cd concentration of rhizosphere soil bacterial community explained 13.8%, 4.7%, 32.6%, and 3.6% of variation in B. napus (Additional file 2: Fig. S6B), and 14.7%, 3.4%, 35.7%, and 5.3% of variation in B. juncea (Additional file 2: Fig. S6C), respectively. For bulk soils bacterial community, VPA indicated that pH, soil nutrients, and Cd explained 4.9%, 39.4%, and 4.4% variation in B. napus (Fig. 6B) and 5.1%, 42.3%, and 6.1% variation in B. juncea (Fig. 6C), respectively.
In this research, we comprehensively studied the effect of Cd contamination on soil-plant system by comparing the variance of bacterial community under different Cd treatments. The results reflected that Cd could inhibit the growth of both species of oilseed rapes, and change the bacterial community structure in soil-plant ecosystem, especially in soil.
Cd could be accumulated in plants from soil via root absorption (Khan, et al., 2017). In the current study, Cd content in B. napus and B. juncea tissues both increased with increasing Cd levels (Additional file 2: Fig. S1). Hyperaccumulators generally have BAF values higher than 1 (Gascó et al., 2019). The results revealed that BAF in both B. napus and B. juncea, leaves and roots, were higher than 1 indicating that they are accumulator plants (shoots and roots) for Cd (Additional file 1: Table S3). TF in two oilseed rapes species was significantly decreased with elevated levels of Cd, which indicated the translocation capability of rapeseeds decreased under higher Cd treatment (Additional file 1: Table S3). Zeng et al. (2020) also found that TF values were decreased when the Cd concentration was elevated. It should be noted that Cd accumulation in tissues and TF values in B. napus were significantly higher than B. juncea. It may be due to different adsorption coefficients of Cd by different plant (Rattan et al., 2005).
Additionally, the accumulation of Cd in plants could also cause plant growth inhibition (Fig. 1). The biomass was declined with increasing Cd levels. This inhibition might be due to Cd inhibiting photosynthesis, retarding plant growth, since the content of chlorophyll decreased with increasing Cd concentration (Fig. 1D). One recent report demonstrated that total chlorophyll concentrations decreased with increasing Cd concentrations, which was accompanied by the reduction of biomass (Liu et al., 2011). However, in comparison between the two species of oilseed rapes, biomass of B. juncea was markedly higher than B. napus under the higher Cd concentration (30 mg/kg) while B. napus was significantly higher than B. juncea under CK treatment. The result indicated that B. juncea might have stronger Cd resistance than B.napus.
When plants experience Cd stress, antioxidative enzymes are induced (Belimov et al., 2007). In the current study, SOD and POD activities were elevated under the higher Cd concentration compared with CK, especially for SOD (Fig. 1G), while CAT was opposite (Fig. 1I). SOD is the first line to defend against ROS, so its activity increasing means that ROS production was enhanced. Increases in SOD activity could enhance H2O2 production, which is accompanying with an inhibition of CAT (Giansoldati et al., 2012). In this study, SOD and POD activities were higher and CAT activity was lower in B. napus than B. juncea, which indicates Cd may cause more damage to cells in B. napus.
The above results revealed that B. juncea may have relatively higher Cd-tolerance than B. napus. Molnár et al. (2020) found that higher levels of zinc oxide nanoparticles (ZnO NPs) were toxic to oilseed rapes and B. juncea has higher tolerance to ZnO NPs than B. napus. This may be due to the higher chlorophyll content of B. juncea under all treatments (Fig. 1D).
Cd, which is extensively found in agricultural soils (Guo et al., 2017), might interfere with the metabolic activity of local resident microorganisms, resulting in modified soil conditions (Lu et al., 2013). In this study, a portion of soil physiochemical properties was altered under 30 mg/kg Cd contamination (Fig. 2). The result showed that pH was lower with the increase of Cd levels and it might be closely correlated with microbial activity (Deng et al., 2018). Nitrogen is an essential nutrient for all life forms and could limit the primary productivity in many ecosystems (Frey et al., 2004). Several types of nitrogen content (NH3-N, NO3-N, and TN) were significantly decreased under higher Cd concentration (30mg/kg) (Fig. 2), indicating that Cd may detrimentally influence soil quality, possibly be due to a detrimental effect of the heavy metal on the nitrogen cycle. Work by Sarria Carabalí et al. (2020) showed that Cd contamination has serious effects on N mobilization.
The soil microorganism plays a role in soil ecosystem functions (Schulz et al., 2013). The α-diversity of soil bacterial communities decreased under higher Cd concentration in soils (Fig. 3). The relative abundance of Actinobacteria was significantly decreased and Proteobacteria and Bacteroidetes were significantly increased in the B. napus rhizosphere. This is in agreement with Wu et al. (2018) who found that Actinobacteria abundance was significantly decreased under Cd stress, with a increase of Proteobacteria and Firmicutes.
Cd pollution could change the relative abundance on genus level (Fig. 4B) including Sphingomonas sp., Lysobacter sp., and Ralstonia sp. which increased with the increase of Cd levels and these genera could resist Cd pollution (Guo, et al., 2017; Nies, 2000; Tipayno et al., 2018). Additionally, the relative abundance of Gemmatimonas sp., Arthrobacter sp., Haliangium sp., and Terriglobus sp., decreased with increasing Cd concentration. The increase of resistant groups and the reduction of sensitive groups might cause the change (Singh et al., 2014).
Results from dissimilarity tests and PCoA indicated that the soil bacterial community structures were obviously altered under the higher Cd concentration (30 mg/kg) (Additional file 1: Table S19 and Table S20; Fig. 5B). Hou, et al. (2018) found that α-diversity was significantly decreased and bacterial community structure was altered in rice when Cd levels were elevated. Wood, et al. (2016) also found that Cd significantly changed bacterial community structure with no reduction in bacterial number in the Cd-accumulating plant. Results of Mantel test, CCA, and CCA-based VPA (Table 1 and Additional file 1: Table S23, Additional file 2: Fig. S6 and Fig. 6) suggested the bacterial community in soil samples is most significantly and negatively correlated with Cd stress, which was in accordance with previous findings (Beattie, et al., 2018).
There were no significant affects to the bacterial α-diversity or community structure under 10 mg/kg Cd stress in soils, possibly due to the presence of the oilseed rapes. Wang et al. (2018) found that Cd stress reduced microbial α-diversity, while planting with oilseed rape increased the microbial α-diversity, alleviating the toxic effects of Cd in soils. The result revealed that both of B. napus and B. junca could reduce the toxicity of Cd to bacteria.
Although many studies have concerned the influence of Cd on soil microbial communities, less research has paid attention on plant microbial community (phyllospheric and endophytic) affected by contaminant. The phyllosphere and endosphere bacteria are important in remediation of heavy metal pollutions (Jia, et al., 2018; Wang, et al., 2020b) and some of them could promote the plant growth (Bulgarelli et al., 2013). In this study, there are no significant variation for phyllospheric and leaf endophytic bacterial diversities, but root endophytic bacterial diversities under Cd stress was altered (Additional file 2: Fig. S4). Meanwhile, Cd contamination could change the relative abundance of some genera (Fig. 4A). Some of these bacteria had reported that are beneficial for soil-plant ecosystems. For example, Ralstonia sp. and Methylobacterium sp. were Cd-resistant bacteria (Lodewyckx et al., 2002; Nies, 2000). The variation of the above bacteria under Cd imposition may lead to the change of the ecological function in the soil-plant ecosystem.
This study demonstrated the effect of Cd on the soil-plant ecosystem. Cd suppressed the growth of oilseed rapes (B.napus and B. juncea) and the physiological traits of plants were changed under Cd stress. Meanwhile, Cd primarily affected pH and nitrogen in soil physicochemical properties. Accumulation of Cd in plants increased and TF was decreased with increasing Cd level in soils, and B. napus and B. juncea were accumulators for Cd because of BAF higher than 1. The bacterial community in soil were significantly altered under higher Cd treatments while no significant difference was found for plant bacterial community (except for root endophytic bacteria). Additionally, B. juncea may have stronger Cd tolerance than B. napus under higher Cd concentration. This research offers a new perspective for the impact of contaminants on soil-plant system and might be help improve phytoremediation efficiency.
Authors` Contribution
D.J. and M.Y. conceived the experiments; Y.D., D.Zhang and D.Zhou designed the experiments; Y.D., D.Zhang and J.W. performed the experiments; Y.D., D.Zhou, J.W., L.L. and H.C. analyzed the data; Y.D., wrote the paper and D.J., D.Zhang with M.Y. revised and perfected the paper. All authors have read and agreed to the published version of the manuscript.
Acknowledgements
The authors are thankful to School of Life Sciences, Hunan University of Science and Technology and Key Laboratory of Environmental Biotechnology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences for their support in this experiment.
Competing interests
The authors declear that they have no competing interests.
Avaiability of data and material
All data obtained have been included into the manuscript and its additional files.
Consent for publication
Not applicable.
Ethics approval and consent participate
Not applicable.
Funding
This work was supported by National Key Research and Development Program of China (grant number 2018YFD1000904), the Scientific Research Fund of Hunan Provincial Education Department (grant number 17K035), the opening funding of Hunan Key Laboratory of Economic Crops Genetic Improvement and Integrated Utilization (grant number E22011), Guangxi Key Laboratory of Biology for Crop Diseases and Insect Pests (grant number 2019-KF-03) and the opening funding of Jiangsu Key Laboratory for Food Quality and Safety-State Key Laboratory Cultivation Base, Ministry of Science and Technology (grant number 028074911709).
Publisher´s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations
Open Access
This article is distributed under the terms of the Creative Commons Attribution 4.0 International License(http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were mad.
- Beattie RE, Henke W, Campa MF, Hazen TC, McAliley LR, Campbell JH (2018) Variation in microbial community structure correlates with heavy-metal contamination in soils decades after mining ceased. Soil Biol Biochem 126: 57-63.
- Belimov AA, Safronova VI, Demchinskaya SV, Dzyuba OO (2007) Intraspecific variability of cadmium tolerance in hydroponically grown Indian mustard (Brassica juncea (L.) Czern.) seedlings. Acta Physiol Plant 29(5): 473-478.
- Bulgarelli D, Schlaeppi K, Spaepen S, Ver Loren van Themaat E, Schulze-Lefert P (2013) Structure and functions of the bacterial microbiota of plants. Annu Rev Plant Biol 64: 807-38.
- Cheraghi-Aliakbari S, Beheshti-Alagha A, Ranjbar F, Nosratti I (2020) Comparison of Myagrum perfoliatum and Sophora alopecuroides in phytoremediation of Cd- and Pb-contaminated soils: A chemical and biological investigation. Chemosphere 259: 127450.
- DalCorso G, Fasani E, Manara A, Visioli G, Furini A (2019) Heavy Metal Pollutions: State of the Art and Innovation in Phytoremediation. Int J Mol Sci 20(14):
- Deng S, Ke T, Li L, Cai S, Zhou Y, Liu Y, Guo L, Chen L, Zhang D (2018) Impacts of environmental factors on the whole microbial communities in the rhizosphere of a metal-tolerant plant: Elsholtzia haichowensis Sun. Environ Pollut 237: 1088-1097.
- Edgar RC (2013) UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods 10(10): 996-8.
- Frey SD, Knorr M, Parrent JL, Simpson RT (2004) Chronic nitrogen enrichment affects the structure and function of the soil microbial community in temperate hardwood and pine forests. For Ecol Manag 196(1): 159-171.
- Gascó G, Álvarez ML, Paz-Ferreiro J, Méndez A (2019) Combining phytoextraction by Brassica napus and biochar amendment for the remediation of a mining soil in Riotinto (Spain). Chemosphere 231: 562-570.
- Giansoldati V, Tassi E, Morelli E, Gabellieri E, Pedron F, Barbafieri M (2012) Nitrogen fertilizer improves boron phytoextraction by Brassica juncea grown in contaminated sediments and alleviates plant stress. Chemosphere 87(10): 1119-25.
- Goswami S, Das S (2015) A Study on Cadmium Phytoremediation Potential of Indian Mustard, Brassica juncea. Int J Phytoremediation 17(1-6): 583-8.
- Guo H, Nasir M, Lv J, Dai Y, Gao J (2017) Understanding the variation of microbial community in heavy metals contaminated soil using high throughput sequencing. Ecotoxicol Environ Saf 144: 300-306.
- Hou D, Wang R, Gao X, Wang K, Lin Z, Ge J, Liu T, Wei S, Chen W, Xie R, Yang X, Lu L, Tian S (2018) Cultivar-specific response of bacterial community to cadmium contamination in the rhizosphere of rice (Oryza sativa L.). Environ Pollut 241: 63-73.
- Jia T, Guo T, Cao M, Chai B (2018) Effects of Heavy Metals on Phyllosphere and Rhizosphere Microbial Community of Bothriochloa ischaemum. Appl Sci 8(9):
- Jing R, Kjellerup BV (2018) Biogeochemical cycling of metals impacting by microbial mobilization and immobilization. J Environ Sci 66: 146-154.
- Khan MA, Khan S, Khan A, Alam M (2017) Soil contamination with cadmium, consequences and remediation using organic amendments. Sci Total Environ 601-602: 1591-1605.
- Kong X, Jin D, Jin S, Wang Z, Yin H, Xu M, Deng Y (2018) Responses of bacterial community to dibutyl phthalate pollution in a soil-vegetable ecosystem. J Hazard Mater 353: 142-150.
- Kong Y (2011) Btrim: a fast, lightweight adapter and quality trimming program for next-generation sequencing technologies. Genomics 98(2): 152-3.
- Li Y, Wang S, Zhang Q, Zang F, Nan Z, Sun H, Huang W, Bao L (2018) Accumulation, interaction and fractionation of fluoride and cadmium in sierozem and oilseed rape (Brassica napus L.) in northwest China. Plant Physiol Biochem 127: 457-468.
- Li Z, Ma Z, van der Kuijp TJ, Yuan Z, Huang L (2014) A review of soil heavy metal pollution from mines in China: pollution and health risk assessment. Sci Total Environ 468-469: 843-53.
- Liu H, Xie Y, Li J, Zeng G, Li H, Xu F, Feng S, Xu H (2020) Effect of Serratia sp. K3 combined with organic materials on cadmium migration in soil-vetiveria zizanioides L. system and bacterial community in contaminated soil. Chemosphere 242: 125164.
- Liu YT, Chen ZS, Hong CY (2011) Cadmium-induced physiological response and antioxidant enzyme changes in the novel cadmium accumulator, Tagetes patula. J Hazard Mater 189(3): 724-31.
- Lodewyckx C, Mergeay M, Vangronsveld J, Clijsters H, Van der Lelie D (2002) Isolation, characterization, and identification of bacteria associated with the zinc hyperaccumulator Thlaspi caerulescens subsp. calaminaria. Int J Phytoremediation 4(2): 101-15.
- Lu M, Xu K, Chen J (2013) Effect of pyrene and cadmium on microbial activity and community structure in soil. Chemosphere 91(4): 491-7.
- Magoč T, Salzberg SL (2011) FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27(21): 2957-63.
- Mitra S, Pramanik K, Ghosh PK, Soren T, Sarkar A, Dey RS, Pandey S, Maiti TK (2018) Characterization of Cd-resistant Klebsiella michiganensis MCC3089 and its potential for rice seedling growth promotion under Cd stress. Microbiol Res 210: 12-25.
- Molnár Á, Papp M, Zoltán Kovács D, Bélteky P, Oláh D, Feigl G, Szőllősi R, Rázga Z, Ördög A, Erdei L, Rónavári A, Kónya Z, Kolbert Z (2020) Nitro-oxidative signalling induced by chemically synthetized zinc oxide nanoparticles (ZnO NPs) in Brassica species. Chemosphere 251: 126419.
- Nies DH (2000) Heavy metal-resistant bacteria as extremophiles: molecular physiology and biotechnological use of Ralstonia sp. CH34. Extremophiles 4: 77-82.
- Rattan RK, Datta SP, Chhonkar PK, Suribabu K, Singh AK (2005) Long-term impact of irrigation with sewage effluents on heavy metal content in soils, crops and groundwater-a case study. Agric Ecosyst Environ 109(3-4): 310-322.
- Rizwan M, Ali S, Zia Ur Rehman M, Rinklebe J, Tsang DCW, Bashir A, Maqbool A, Tack FMG, Ok YS (2018) Cadmium phytoremediation potential of Brassica crop species: A review. Sci Total Environ 631-632: 1175-1191.
- Rossi G, Figliolia, A., Socciarelli, S., Pennelli, B. (2002) Capability of Brassica napus to Accumulate Cadmium, Zinc and Copper from Soil. Acta Biotechnol 22(1-2): 133-440.
- Sarria Carabalí MM, García-Oliva F, Cortés Páez LE, López-Lozano NE (2020) Effect of cadmium contamination on the rhizosphere bacterial diversity of Echinocactus platyacanthus. Rhizosphere 13:
- Schulz S, Brankatschk R, Dümig A, Kögel-Knabner I, Schloter M, Zeyer J (2013) The role of microorganisms at different stages of ecosystem development for soil formation. Biogeosciences 10(6): 3983-3996.
- Singh BK, Quince C, Macdonald CA, Khachane A, Thomas N, Al-Soud WA, Sorensen SJ, He Z, White D, Sinclair A, Crooks B, Zhou J, Campbell CD (2014) Loss of microbial diversity in soils is coincident with reductions in some specialized functions. Environ Microbiol 16(8): 2408-20.
- Tipayno SC, Truu J, Samaddar S, Truu M, Preem JK, Oopkaup K, Espenberg M, Chatterjee P, Kang Y, Kim K, Sa T (2018) The bacterial community structure and functional profile in the heavy metal contaminated paddy soils, surrounding a nonferrous smelter in South Korea. Ecol Evol 8(12): 6157-6168.
- Wang P, Peng H, Liu J, Zhu Z, Bi X, Yu Q, Zhang J (2020a) Effects of exogenous dissolved organic matter on the adsorption-desorption behaviors and bioavailabilities of Cd and Hg in a plant-soil system. Sci Total Environ 728: 138252.
- Wang Q, Garrity GM, Tiedje JM, Cole JR (2007) Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol 73(16): 5261-7.
- Wang Q, Ge C, Xu S, Wu Y, Sahito ZA, Ma L, Pan F, Zhou Q, Huang L, Feng Y, Yang X (2020b) The endophytic bacterium Sphingomonas SaMR12 alleviates Cd stress in oilseed rape through regulation of the GSH-AsA cycle and antioxidative enzymes. BMC Plant Biol 20(1): 63.
- Wang X, Bai J, Wei T, Feng G, Zhao H, Wei W, Wang M, Zhao Y (2018) Oilseed rape cultivation increases the microbial richness and diversity in soils contaminated with cadmium. J Soils Sediments 18(7): 2451-2462.
- Wood JL, Zhang C, Mathews ER, Tang C, Franks AE (2016) Microbial community dynamics in the rhizosphere of a cadmium hyper-accumulator. Sci Rep 6: 36067.
- Wu B, Hou S, Peng D, Wang Y, Wang C, Xu F, Xu H (2018) Response of soil micro-ecology to different levels of cadmium in alkaline soil. Ecotoxicol Environ Saf 166: 116-122.
- Xia C, Luo Y, Yao J, Liu W, Wang F, Wu X (2018) Joint effects of Cd and thioglycollic acid on soil microbial activity. Int Biodeterior Biodegradation 128: 164-170.
- Zeng X, Zou D, Wang A, Zhou Y, Liu Y, Li Z, Liu F, Wang H, Zeng Q, Xiao Z (2020) Remediation of cadmium-contaminated soils using Brassica napus: Effect of nitrogen fertilizers. J Environ Manage 255: 109885.
- Zhang WH, Huang Z, He LY, Sheng XF (2012) Assessment of bacterial communities and characterization of lead-resistant bacteria in the rhizosphere soils of metal-tolerant Chenopodium ambrosioides grown on lead-zinc mine tailings. Chemosphere 87(10): 1171-8.
Table 1 Mantel analysis of the relationship between the soil bacterial community structure in B. napus and environmental factors based on Bray Curtis (BC) and Jaccard (JC) method.
|
Samples |
Environmental factors |
r.BC |
p.BC |
r.JC |
p.JC |
|
Rhizosphere soil |
height |
0.6394 |
0.001 |
0.3089 |
0.001 |
|
weight |
0.5874 |
0.001 |
0.2704 |
0.001 |
|
|
leaf_area |
0.5709 |
0.001 |
0.2964 |
0.001 |
|
|
pH |
0.2294 |
0.031 |
0.0347 |
0.348 |
|
|
TOC |
0.1250 |
0.192 |
0.0323 |
0.336 |
|
|
TN |
0.1993 |
0.056 |
0.0037 |
0.421 |
|
|
TP |
0.0518 |
0.293 |
-0.0074 |
0.475 |
|
|
NH3_N |
0.1902 |
0.075 |
-0.0443 |
0.567 |
|
|
NO3_N |
0.2159 |
0.028 |
-0.0811 |
0.711 |
|
|
AP |
-0.1334 |
0.899 |
0.2172 |
0.081 |
|
|
AK |
-0.0977 |
0.777 |
-0.1445 |
0.804 |
|
|
Soil_Cd |
0.6319 |
0.001 |
0.3198 |
0.003 |
|
|
Bulk soil
|
pH |
0.1710 |
0.092 |
0.0154 |
0.412 |
|
TOC |
0.0624 |
0.361 |
0.0554 |
0.326 |
|
|
TN |
0.2064 |
0.035 |
0.1129 |
0.158 |
|
|
TP |
-0.0074 |
0.507 |
-0.0625 |
0.631 |
|
|
NH3_N |
-0.1643 |
0.886 |
-0.0580 |
0.608 |
|
|
NO3_N |
0.0963 |
0.143 |
0.0007 |
0.482 |
|
|
AP |
0.1708 |
0.14 |
-0.0058 |
0.481 |
|
|
AK |
0.1677 |
0.108 |
-0.0081 |
0.481 |
|
|
Soil_Cd |
0.2343 |
0.018 |
0.0764 |
0.243 |