# The Interactions among the Heavy Metals in Soils and in Weeds and Their Antioxidant Capacity under the Mining Activities in Thai Nguyen Province, Vietnam

Authors: Ngoc-Lien Nguyen; Cam-Tu Vu; Hien-Minh To; Hoang-Nam Pham; Hai-Dang Nguyen; Tien-Dat Nguyen; Kieu-Oanh Nguyen Thi
Journal: Journal of Chemistry
 (2020)

---

## Abstract

In this study, the relationship between heavy metals (HMs) concentrations in soils and several weeds includingCyclosorus parasiticus, Dicranopteris linearis, Pityrogramma calomelanos, and Pteris vittata in three mining sites (Cam Gia (Thai Nguyen city), Tan Long (Dong Hy district), and Ha Thuong (Dai Tu district)) in Thai Nguyen province, Vietnam, have been investigated. The levels of HMs varied among soil origins and showed the contaminations of As, Cu, and Pb in soil samples collected in Dong Hy and Dai Tu districts. In addition, the HM distribution and cocontamination phenomena in different soils significantly affected the HM residues and transportation abilities into different species as well as tissues. Moreover, based on the analysis of bioaccumulation factor (BF) and translocation factor (TF), C. parasiticus and D. linearis were found potentially for phytoextraction by roots, while P. calomelanos and P. vittata were suitable for hyperaccumulation in shoots and leaves. Consequently, the strongest antioxidant property by 1,1-diphenyl-2-picrylhydrazyl (DPPH) and superoxide anion (SRSA) radical scavenging assays were demonstrated in the methanol root extracts of C. parasiticus and P. vittata, respectively. In conclusion, the correlation among HM in soils and tissues with antioxidant property allows us to hypothesize that the presence of these elements can enhance the antioxidant activity of plant extracts, suggesting to apply the weeds in phytoremediation as well as in phytomedicine.


---

## Body

## 1. Introduction

Heavy-metal (HM) pollution is raising alarms since they can cause massive development in industrialization and mining, especially in developing countries. Although HM explorations bring many economic benefits, the serious pollution by exceeding HMs in the soil and water occurs as a matter of obsolete technology or the lack of systematic mining process and preliminary treatment. In the world, up to 80% of industrial polluted sites, which cause health-related threats and HM pollutants such as As, Pb, Cd, Zn, are listed in the top 12 red-flag contaminants [1]. In Northern Vietnam, several mining areas in Thai Nguyen province are famous for their rich mineral resources of Fe, Ti, Zn, Cu, Pb, ferrous, and nonferrous ores. The studies of soil analysis in these regions reportedly contained several HMs contaminated sites in the scenario of the foreign firm’s accelerating exploitation and outdated manual mining technology [2]. Analysis of soil samples in the Pb-Zn ore Dong Hy, Thai Nguyen showed concentrations of the elements Pb, Zn, Cd, and As that exceeded many times the permissible standards [3, 4]. In Dai Tu district, mining activities at local craft have created a significant amount of waste tailings and waste rock. Results of several waste rock samples have shown that the average As levels reached 5000 mg/kg. Other HM contents in the samples are also very high (Cu 1260 mg/kg; Pb 105 mg/kg; Cd 0.5 mg/kg; Se 17 mg/kg) which indicates that there is cocontamination of several HMs in a site [5]. Since As contaminated sites have been proved to simultaneously exist with high levels of Cu, Cr, Cd, Pb, and Zn [1, 6], the exceeding HMs distributed to agricultural soil (paddy, tea, and forest fields) surrounding tin and tungsten mines in Dai Tu district, present at up to 1000 mg/kg for As, 300 mg/kg for Cu, and abnormal high level of other HMs including Cd, Zn, Pb, Mn. Other soil analysis in the mining sites of Northern Vietnam also shows cocontamination of Pb, As, Zn, Cd, Cu with the factors of 314, 2990, 280, 28, 11 times, respectively, higher the level of Vietnamese standard for mining soil (QCVN 03-MT: 2015/BTNMT 2015) [7, 8].Excess of HMs in the soil can cause phytotoxic effects in plants in several ways including the excessive production of free radicals and reactive oxygen species (ROS). ROS being highly unstable could play a dual role: (i) damage cellular components and (ii) act as an important secondary messenger for inducing plant defense systems, which affect plant growth and biomass [9]. To counteract this damage, plants are equipped with enzymatic and nonenzymatic defense mechanisms. Several plants are known as hyperaccumulators since they are capable of growing in soils with very high concentrations of metals, absorbing these metals through their roots, and concentrating extremely high levels of metals in their tissues [10, 11]. These hyperaccumulators possess several antioxidative defense systems to scavenge toxic free radicals in order to protect themselves from the oxidative stress that is caused by HMs.Interestingly, many hyperaccumulators belong to pteridophytes that are important plants from a phylogenetic and evolutionary perspective [12]. Botanically, these temperate-thriving plants were wild and many of them were weeds a thousand years ago. As these plants have survived since Paleozoic times, pteridophytes have survived more disruptive changes in the environment than any other vascular plants. This can be explained as the reason why these ferns were considered as hyperaccumulators. C. parasiticus, D. linearis, P. calomelanos, and P. vittata are four ferns that dominate in mining sites in Thai Nguyen province, Vietnam, in which the two latter weeds, P. calomelanos, and P. vittata, have been reported as strong arsenic hyperaccumulators [6, 13–15].In order to explore the relationship between heavy-metal concentrations in the metalliferous soils on the heavy-metal uptake, capabilities of the corresponding local several dominant ferns includingC. parasiticus, D. linearis, P. calomelanos, and P. vittata collected in the three mining sites in Thai Nguyen province, Vietnam, were evaluated for their antioxidant activity via 1,1-diphenyl-2-picrylhydrazyl (DPPH) and superoxide anion (SRSA) radical scavenging assays. The screening results can be a good start for further studies of the antioxidants which can be applied in the health industry. Consequently, the investigation of HMs contaminated levels can be useful for applying some methods in order to reduce pollution and/or protect the environment. We also try to find out the relationship between antioxidant activity of plants and HMs levels in soil in order to incorporate the environment purpose and natural sources valorization.


## 2. Materials and Methods

### 2.1. Chemicals

Hydrochloric acid, HCl (36.5–38%), nitric acid, HNO3 (63%), sulfuric acid, H2SO4 (98%), HClO4, DPPH, methanol (MeOH, HPLC grade), ascorbic acid (>99%), catechin, and xanthine were purchased from Sigma Aldrich (Taufkirchen, Germany) in analytical grade. Ultrapure water (UPW) was prepared using a Millipore (Milli-Q) purification system (Millipore GmbH, Schwalbach, Germany, 18.2 MΩ·cm, TOC < 2μg/L).


### 2.2. Soil and Plant Sampling

Four dominant plant species were found at three mine sites that were described in Table S1.C. parasiticus (Thelypteridaceae) were collected at the Campus of Thai Nguyen Iron and Steel Corporation, Cam Gia Ward, Thai Nguyen City, Thai Nguyen province. Two other weeds, D. linearis (Gleicheniaceae) and P. calomelanos (Adiantaceae), were collected in Ha Thuong commune, Dai Tu district, Thai Nguyen province. P. vittata (Pteridaceae) was collected in Hich village, Tan Long commune, Dong Hy district, Thai Nguyen province, Vietnam. All plant samples were identified by Dr. Nguyen The Cuong, Institute of Ecology and Biological Resources, Vietnam Academy of Science and Technology. The voucher of specimens was deposited in the Department of Life Sciences, University of Science and Technology of Hanoi. The corresponding soil samples were dug in the rhizospheric area (up to 20 cm depth from the topsoil and 10 cm diameter from plant position) of the collected plants. All plant and soil samples were collected under dry weather conditions, March 9th, 2015, stored in polyethylene bags, labeled appropriately.For sample preparation, each plant was divided into 3 parts: leaves, shoots, and roots. All plant tissues were washed by tap water to remove all dust and soil residues and then rinsed by UPW. All parts were cut down into smaller parts and dried in the oven at 55°C for 48 h to achieve the stable weight of samples. Similarly, soil samples were also prepared in the same manner of drying in the oven at 105°C for 24 h until constant weight after removing all the moisture. Then, plant and soil samples were crushed by using the clean and dry mortar, followed by screening steps using 2 mm and 0.18 mm sieves, respectively. The prepared samples were finally stored at −80°C until analysis.


### 2.3. Atomic Absorption Spectrometry (AAS) Technique

To determine the Pb, Cu, Cd, and Zn concentrations in soil samples: 15 mL HCl 37% and 5 mL HNO 3 65% were added to 2 g of dried and ground soil samples which were contained in 50 mL baker. The mixtures were digested on a hot plate at 95°C until getting the transparent solutions and then let to room temperature and prepared for analysis. Especially for As analysis: 5 mL HNO 3 65% was added to 2 g of prepared soil samples, and then the mixture was treated in a microwave oven for 40 min. The sample was cooled down to room temperature. After the treatment, all five HMs in soil samples were formed in the dissolved solution and filtered through 0.45μm cellulose acetate membrane (Whatman, England) and ready for analysis.In order to determine the HM concentration in the plants, each 0.5 g of each tissue powder was digested by a mixture including 3 mL concentrated H2SO4 98%, 0.5 mL HClO4, and 10 mL of concentrated HNO3. The mixture was then heated in the oven at temperature between 150 and −200°C for at least 120 min. The operation was repeated until the digested solutions became clear and homogenous. Then, the obtained solution was normed up to 25 mL in the Erlenmeyer flask. All the final sample solutions were shaken well at 180 rpm on IKA® C-MAG HS 4 magnetic stirrer (Cologne, Germany) for 20 min, filtered by 0.45μm cellulose acetate membrane to remove all the particulars, and then diluted to 25 mL for measuring Cu, Pb, Zn, Cd concentration. For As concentration measurement, 15 mL of the solution on the surface was made up with 10 mL UPW to get finally 25 mL for measuring.All measurements were done in triplicate by the Atomic Absorption Spectrometry method, proposed by APHA (1995) based on the Varian 280FS and AA280Z (Australia). The results were reported by the average values with relative standard deviation (% RSD) lower than 10%.


### 2.4. DPPH Scavenging Assay

In order to evaluate the antioxidant activity of plant tissues, approximately 5.0 g of each tissue powder was extracted with MeOH solution (100 mL × 3 times) at room temperature in an ultrasonic bath Skymen JP-008 (Guangdong, China). The combined extracts were concentrated in Buchi Rotavapor (Flawil, Switzerland) to obtain crude extracts.The DPPH radical scavenging capacity was evaluated according to the previously described method [16]. Briefly, two concentrations of 100 μg/mL and 500 μg/mL of crude extracts dissolved in dimethyl sulfoxide (5 μl) were added to 195 μl DPPH solution (150 mM) in 96-well plates. The mixed solution was incubated at room temperature for 30 min. Then the absorbance of the reaction mixture was measured at 517 nm on a microplate reader iMark™ (BioRad, California, USA). Ascorbic acid was used as a positive control (SC 91.3% at 100 μg/mL). The differences in absorbance between the test sample and the blank (DPPH in MeOH) was calculated and expressed as (%) scavenging capacity:(1)SC%=Ab−AsAb×100,where As is the absorbance of the test sample and Ab is the absorbance of the blank. Samples with SC values above 50% were considered to be active.


### 2.5. Superoxide Anion Radical Scavenging Assay

The superoxide anion scavenging activity was measured based on the described method [17]. In this experiment, crude extracts were prepared at 2 concentrations as 100 μg/mL and 500 μg/mL in DMSO. The standard reaction mixture in a final volume of 1 mL contained 10 μL of sample solution, 490 μL of xanthine (1 mM; in potassium phosphate buffer 0.1 M, pH 7.8), and 490 μL of NBT (0.4 mM, in potassium phosphate buffer 0.1 M, pH 7.8). The reaction was initiated by the addition of 10 μL of xanthine oxidase. Incubation at room temperature for 20 minutes and the samples’ absorbance was read at 560 nm by using microplate reader iMark™ (BioRad, California, USA). Catechin was used as standard control (positive sample) and the blank well without sample was used as a negative sample. All experiments were performed in triplicate and the reported results are the averages of three independent measurements.Superoxide radical scavenging activity (SRSA) is calculated by the following equation:(2)SC%=Ac−AsAc−Ab×100.Ac is the absorbance of the positive sample, As is the absorbance of the test sample, and Ab is the absorbance of the negative sample. Samples with SC values above 50% were considered to be active.


### 2.6. Statistical Analysis

Statistics were performed in R version 3.5.3 [18]. Data was normalized by taking their logarithm base 10. In order to discriminate the accumulation ability of different species, the principal component analysis (PCA) was done by PCA and dudi.pca function in FactoMineR, factoextra 1.0.7, and ade4 R-package [19, 20]. To evaluate the effect of environmental factors (HM in soils) on the accumulation of HM in the plants, a function envfit in vegan 2.5–6 R-package was used [21]. The correlation between the HM concentrations in soil and plants with antioxidant activities was calculated and visualized by Performance Analytics R-package [22].


## 2.1. Chemicals

Hydrochloric acid, HCl (36.5–38%), nitric acid, HNO3 (63%), sulfuric acid, H2SO4 (98%), HClO4, DPPH, methanol (MeOH, HPLC grade), ascorbic acid (>99%), catechin, and xanthine were purchased from Sigma Aldrich (Taufkirchen, Germany) in analytical grade. Ultrapure water (UPW) was prepared using a Millipore (Milli-Q) purification system (Millipore GmbH, Schwalbach, Germany, 18.2 MΩ·cm, TOC < 2μg/L).


## 2.2. Soil and Plant Sampling

Four dominant plant species were found at three mine sites that were described in Table S1.C. parasiticus (Thelypteridaceae) were collected at the Campus of Thai Nguyen Iron and Steel Corporation, Cam Gia Ward, Thai Nguyen City, Thai Nguyen province. Two other weeds, D. linearis (Gleicheniaceae) and P. calomelanos (Adiantaceae), were collected in Ha Thuong commune, Dai Tu district, Thai Nguyen province. P. vittata (Pteridaceae) was collected in Hich village, Tan Long commune, Dong Hy district, Thai Nguyen province, Vietnam. All plant samples were identified by Dr. Nguyen The Cuong, Institute of Ecology and Biological Resources, Vietnam Academy of Science and Technology. The voucher of specimens was deposited in the Department of Life Sciences, University of Science and Technology of Hanoi. The corresponding soil samples were dug in the rhizospheric area (up to 20 cm depth from the topsoil and 10 cm diameter from plant position) of the collected plants. All plant and soil samples were collected under dry weather conditions, March 9th, 2015, stored in polyethylene bags, labeled appropriately.For sample preparation, each plant was divided into 3 parts: leaves, shoots, and roots. All plant tissues were washed by tap water to remove all dust and soil residues and then rinsed by UPW. All parts were cut down into smaller parts and dried in the oven at 55°C for 48 h to achieve the stable weight of samples. Similarly, soil samples were also prepared in the same manner of drying in the oven at 105°C for 24 h until constant weight after removing all the moisture. Then, plant and soil samples were crushed by using the clean and dry mortar, followed by screening steps using 2 mm and 0.18 mm sieves, respectively. The prepared samples were finally stored at −80°C until analysis.


## 2.3. Atomic Absorption Spectrometry (AAS) Technique

To determine the Pb, Cu, Cd, and Zn concentrations in soil samples: 15 mL HCl 37% and 5 mL HNO 3 65% were added to 2 g of dried and ground soil samples which were contained in 50 mL baker. The mixtures were digested on a hot plate at 95°C until getting the transparent solutions and then let to room temperature and prepared for analysis. Especially for As analysis: 5 mL HNO 3 65% was added to 2 g of prepared soil samples, and then the mixture was treated in a microwave oven for 40 min. The sample was cooled down to room temperature. After the treatment, all five HMs in soil samples were formed in the dissolved solution and filtered through 0.45μm cellulose acetate membrane (Whatman, England) and ready for analysis.In order to determine the HM concentration in the plants, each 0.5 g of each tissue powder was digested by a mixture including 3 mL concentrated H2SO4 98%, 0.5 mL HClO4, and 10 mL of concentrated HNO3. The mixture was then heated in the oven at temperature between 150 and −200°C for at least 120 min. The operation was repeated until the digested solutions became clear and homogenous. Then, the obtained solution was normed up to 25 mL in the Erlenmeyer flask. All the final sample solutions were shaken well at 180 rpm on IKA® C-MAG HS 4 magnetic stirrer (Cologne, Germany) for 20 min, filtered by 0.45μm cellulose acetate membrane to remove all the particulars, and then diluted to 25 mL for measuring Cu, Pb, Zn, Cd concentration. For As concentration measurement, 15 mL of the solution on the surface was made up with 10 mL UPW to get finally 25 mL for measuring.All measurements were done in triplicate by the Atomic Absorption Spectrometry method, proposed by APHA (1995) based on the Varian 280FS and AA280Z (Australia). The results were reported by the average values with relative standard deviation (% RSD) lower than 10%.


## 2.4. DPPH Scavenging Assay

In order to evaluate the antioxidant activity of plant tissues, approximately 5.0 g of each tissue powder was extracted with MeOH solution (100 mL × 3 times) at room temperature in an ultrasonic bath Skymen JP-008 (Guangdong, China). The combined extracts were concentrated in Buchi Rotavapor (Flawil, Switzerland) to obtain crude extracts.The DPPH radical scavenging capacity was evaluated according to the previously described method [16]. Briefly, two concentrations of 100 μg/mL and 500 μg/mL of crude extracts dissolved in dimethyl sulfoxide (5 μl) were added to 195 μl DPPH solution (150 mM) in 96-well plates. The mixed solution was incubated at room temperature for 30 min. Then the absorbance of the reaction mixture was measured at 517 nm on a microplate reader iMark™ (BioRad, California, USA). Ascorbic acid was used as a positive control (SC 91.3% at 100 μg/mL). The differences in absorbance between the test sample and the blank (DPPH in MeOH) was calculated and expressed as (%) scavenging capacity:(1)SC%=Ab−AsAb×100,where As is the absorbance of the test sample and Ab is the absorbance of the blank. Samples with SC values above 50% were considered to be active.


## 2.5. Superoxide Anion Radical Scavenging Assay

The superoxide anion scavenging activity was measured based on the described method [17]. In this experiment, crude extracts were prepared at 2 concentrations as 100 μg/mL and 500 μg/mL in DMSO. The standard reaction mixture in a final volume of 1 mL contained 10 μL of sample solution, 490 μL of xanthine (1 mM; in potassium phosphate buffer 0.1 M, pH 7.8), and 490 μL of NBT (0.4 mM, in potassium phosphate buffer 0.1 M, pH 7.8). The reaction was initiated by the addition of 10 μL of xanthine oxidase. Incubation at room temperature for 20 minutes and the samples’ absorbance was read at 560 nm by using microplate reader iMark™ (BioRad, California, USA). Catechin was used as standard control (positive sample) and the blank well without sample was used as a negative sample. All experiments were performed in triplicate and the reported results are the averages of three independent measurements.Superoxide radical scavenging activity (SRSA) is calculated by the following equation:(2)SC%=Ac−AsAc−Ab×100.Ac is the absorbance of the positive sample, As is the absorbance of the test sample, and Ab is the absorbance of the negative sample. Samples with SC values above 50% were considered to be active.


## 2.6. Statistical Analysis

Statistics were performed in R version 3.5.3 [18]. Data was normalized by taking their logarithm base 10. In order to discriminate the accumulation ability of different species, the principal component analysis (PCA) was done by PCA and dudi.pca function in FactoMineR, factoextra 1.0.7, and ade4 R-package [19, 20]. To evaluate the effect of environmental factors (HM in soils) on the accumulation of HM in the plants, a function envfit in vegan 2.5–6 R-package was used [21]. The correlation between the HM concentrations in soil and plants with antioxidant activities was calculated and visualized by Performance Analytics R-package [22].


## 3. Results and Discussion

### 3.1. HMs Level in Soil

Five HMs As, Cd, Zn, Pb, and Cu were quantified in 25 soil samples collected at the same positions with 25 plant samples including 5 samples ofC. parasiticus (in Thai Nguyen city), 6 samples for each type of D. linearis and P. calomelanos plants (in Dai Tu district), and 8 samples of P. vittata (in Dong Hy district). According to each sampling station, the sample abbreviation was used interchangeably for the respective soil and plant samples.Table1 shows the highest As, Zn, Pb, and Cd concentrations in P. vittata correspondent soil while the most concentration of Cu was found in the soil growing C. parasiticus plant. According to the Vietnam National regulations on the allowable limits of HMs in the soils [7] where the maximum allowable limits of As, Zn, Pb, Cu, and Cd in industrial soils are 25; 300; 300; 300, and 10 mg/kg, respectively. It was investigated that soil samples collected in Dong Hy district (P. vittata corresponding soil) were contaminated by both As and Pb, whereas soils from Dai Tu was contaminated with Cu. Soils from Cam Gia, Thai Nguyen city, witnessed the Pb pollution. Significantly, there were no observations of Zn and Cd contaminations in all soil samples collected from different sampling sites.Table 1
The average of five HMs (As, Zn, Pb, Cd, and Cu) concentrations (in mg/kg) of soil samples corresponding to 4 plant species.
CorrespondenceAsZnPbCdCuC. parasiticus9.31 ± 4.21108 ± 31.9307 ± 1381.98 ± 0.70171 ± 10.0D. linearis2.54 ± 0.4584.1 ± 19.4206 ± 56.60.81 ± 0.82453 ± 313P. calomelanos2.92 ± 0.9864.3 ± 15.4239 ± 66.90.38 ± 0.16338 ± 102P. vittata53.9 ± 13.9265 ± 111450 ± 2482.34 ± 1.4611.2 ± 8.06It is found that the distribution of HMs in soils was following the different soil origins and plant bioaccumulation and/or absorption [23]. To understand more this phenomenon, the principal component analysis (PCA) revealed the difference among the HM levels of three locations where 25 soil samples were collected by distinguishing P. vittata, D. linearis, P. calomelanos, and C. parasiticus zones (Figure 1). The first two axes of the PCA explained 66.4% of the variance. The ordination arrow represents the dominance of As and Zn in the soil samples collected in the rhizosphere of P. vittata in Dong Hy district. More Pb and Cd than other HMs were distributed in soils sampled in Dong Hy district and Thai Nguyen city, while Cu significantly existed in soils collected in Dai Tu district. These distributions of HMs in soils have thus corresponded to their detected concentration in different soils (Table 1). In addition, the arrow direction of vectors, which represented the distribution tendency of HMs according to soil samples, were also indicated for the cocontamination relationship between HMs in the soils. Cu vectors tend to contradict the vectors of other metals, indicating that Cu has less ability to coaccumulate concurrently with other metals in the soil. This was identified more clearly based on the Spearman index (Figure S1) when the indices representing the correlation between Cu and other metals were always very low (from −0.76 to −0.36). This phenomenon can be explained by higher solubility and leachability of Cu at the ideal soil pH of 6.0–6.5, compared to other HMs [24]. In contrast, Cd appeared to have a stronger correlation in soil accumulation with other metals, except Cu, while the lowest concentrations of Cd were often found in all soil samples (Table 1). Consequently, the Spearman indices indicating the cocontamination between Cd and other HMs following the order Zn (0.76) > As (0.59) > Pb (0.37). This could be partly clarified for the highest concentrations of Zn, As, Pb, and Cd were found in Dong Hy soil as previous investigation (Table 1). Also, Pb had significant cocontamination with As and Zn in soil following Spearman indices of 0.39 and 0.51, respectively, confirming for high abilities of Zn, As, and Pb accumulations in soils corresponding to P. vittata plant. The analysis results, therefore, presented the significant effects of cocontaminations between soil origins and the HM residences in those soils.Figure 1
PCA of the concentration of HMs (As, Cd, Cu, Pb, and Zn) in 25 soil samples corresponding to 25 plants of 4 species (C. parasiticus, D. linearis, P. calomelanos, and P. vittata). Each point displays a sample; each vector displays the level of an HM. Soil samples are colored and shaped according to different species where they were collected.


### 3.2. HMs Levels in Plant Tissues

The statistical description of HM concentration (in mg/kg) in 75 plant tissue samples was summarized by the average concentration with standard deviation. For each HM, the bioaccumulation factor (BF) was calculated by the ratio between the HM concentration in tissue and those in corresponding soil which is described as the ability of plants for elemental accumulation from the substrate. The efficiency of phytoremediation can be quantified by calculating the translocation factor (TF). This index expresses the capacity of a plant to store the HMs in its upper part, which is defined as the ratio of metal concentration in the upper part to that in the roots [25].The total concentration of selected metals as well as BF and TF indices in plant roots, shoots, and leaves is shown in Table2. In plant roots, the highest As, Pb, and Zn concentrations were found in P. vittata, whereas the most Cu concentration was measured in P. calomelanos. The large amounts of As, Pb, Zn, and Cu found in plant roots of P. vittata and P. calomelanos could be explained by the high level of these metal concentrations detected in the soils where corresponding plants were cultivated (Table 1). On the other side, there was no significant difference in Cd concentration in the roots of 4 plants. This correlated to the low level of Cd found in all corresponding soils (Table 1). On average, the HM concentrations in the roots were much higher than those in the leaves and shoots, implying that fewer HMs were translocated into leaves and shoots for all four species.Table 2
The average concentration (±SD) of five HMs in plant tissues (root, shoot, and leave) from 4 species,C. parasiticus, D. linearis, P. calomelanos, and P. vittata, and corresponding bioaccumulation factor (BF) and translocation factor (TF).
Plant speciesPlant tissueAsCdCuPbZnConc.BFTFConc.BFTFConc.BFTFConc.BFTFConc.BFTFC. parasiticusRoot1.19 ± 1.770.14—0.58 ± 0.490.29—19.7 ± 5.920.12—4.34 ± 1.480.01—121 ± 73.61.12—Shoot1.19 ± 0.900.214.510.66 ± 0.300.3618.413.1 ± 11.10.080.882.79 ± 1.640.011.2828.0 ± 13.00.271.10Leave1.35 ± 1.580.184.430.77 ± 0.300.458.3715.2 ± 8.800.090.815.51 ± 1.950.020.61110 ± 50.71.100.28D. linearisRoot122 ± 32.250.0—0.19 ± 0.090.47—75.8 ± 53.00.47—37.3 ± 20.90.21—35.5 ± 10.80.45—Shoot18.3 ± 14.48.030.320.08 ± 0.050.247.6141.1 ± 22.20.240.388.85 ± 4.490.050.4722.1 ± 10.70.290.99Leave33.4 ± 30.514.90.570.32 ± 0.250.820.4428.4 ± 14.20.170.5915.4 ± 6.220.090.2543.7 ± 31.50.570.59P. calomelanosRoot469 ± 69.5169—0.32 ± 0.110.97—430 ± 75.01.37—160 ± 68.80.78—92.3 ± 19.01.52—Shoot227 ± 71.81020.610.08 ± 0.030.260.2950.1 ± 15.30.160.1213.6 ± 10.90.070.1024.0 ± 2.260.400.27Leave630 ± 27.52291.370.16 ± 0.050.480.5759.0 ± 79.00.180.1716.6 ± 4.360.080.1249.8 ± 9.260.830.55P. vittataRoot1732 ± 299.1403—0.27 ± 0.260.13—5.98 ± 1.301.06—244 ± 60.70.88—188 ± 30.50.35—Shoot331 ± 75.376.20.200.13 ± 0.060.080.845.09 ± 2.370.960.8465.4 ± 18.30.190.2854.2 ± 22.20.230.65Leave438 ± 17699.51.090.33 ± 0.280.181.496.13 ± 3.41.171.03141 ± 51.50.470.59267 ± 1251.133.18In plant shoots, the As concentration ranged from 1.19 ± 0.9 mg/kg (inC. parasiticus) to 331 ± 75.3 mg/kg (in P. vittata). The total concentration (mg/kg) of Cu, Pb, and Zn in plant shoots varied significantly, with the highest values being 59.1 ± 15.3 (in P. calomelanos), 65.4 ± 18.3 (in P. vittata), and 54.2 ± 22.2 (in P. vittata), respectively. Similarly, in the case of metals accumulation in leaves, higher metals concentrations of As, Cu, Pb, and Zn were found in plant leaves of P. vittata and P. calomelanos. On the other side, according to different plant types, a similar distribution tendency of measured Cd concentrations was found in plant roots and shoots, and the lower metals concentrations were often found in all tissues of C. parasiticus and D. linearis compared to other ones. It could be seen that the distribution of metal concentration in plant tissues varied greatly among plant species and had a significant correlation with HM content in soils from different sampling sites.Moreover, the bioaccumulation factor (BF) indicates the ability of metals stored in plant tissues compared to soil [26]. For all plant tissues, the BF values often lower than 1.0 for Cd, Cu, Pb and Zn. In contrast, the BF values of As which were extremely high in all tissues of almost all plants except C. parasiticus with BFAs lower than 1. This indicated the higher moving ability of As from the external environment into plant tissues through the soil exposure route compared to the other 3 metals. It is also found that P. vittata sampled at Dong Hy possessed a higher metal bioaccumulation factor (up to 402.5) than other plants indicating its hyperaccumulating ability. There are a number of reports indicating that P. vittata is an As hyperaccumulator [27, 28] capable of translocating As and storing it in the shoots [29, 30].Similarly, most shoot and root tissues of all almost plants had TF values lower than 1, except TFAs (inC. parasiticus) and TFCd (in C. parasiticus and D. linearis plant) values, which are much higher than 1. In more detail, the maximum values of TFCd and TFAs found in shoots of C. parasiticus plants were 18.4 and 4.51, respectively. Furthermore, the TFmetals indices in the shoots were often higher than that in leaves of C. parasiticus and D. linearis plants, whereas higher TF metals indices were found in leaves of P. calomelanos and P. vittata. High root to shoot translocation of these metals in C. parasiticus and D. linearis indicated vital characteristics to be used in phytoextraction of these metals in shoot [31, 32]. On the other side, the results in TF and BF values of P. calomelanos and P. vittata were considered suitable for phytoextraction, which generally requires translocation of metals easily harvestable leaves tissue to remove the metal contamination from the soil without destroying soil structure and fertility [33, 34]. In fact, plant species were judged as being suited for phytoremediation based on several criteria including wide distribution, high above-ground biomass, high bioaccumulation factors (hyperaccumulators), short life cycles, and high propagation rates [35]. Therefore, it could be concluded that P. vittata and P. calomelanos can be used in phytoremediation in mine sites of Thai Nguyen province as they were considered as native plants that meet all criteria for a hyperaccumulator.


### 3.3. Species-Tissue Discrimination of HMs Levels in Plants

In order to see whether there are differences in the concentration of HMs level (i.e., variables) in samples (i.e., individuals), PCA biplots where both individuals and variables are graphed on the same plot were constructed. In the first two PCs, more than 60% of the variance in the data is explained for all PCA plots. In accordance with the PCA biplot (Figure2(a)), HM levels in 75 plant tissues were grouped in a different pattern based on species. The discrimination of these patterns indicates that four species show a different capacity in uptaking 5 mentioned HMs. The PCA biplots containing only data of HM levels in shoots or leaves showed better separation between four species (Figures 2(c) and 2(d)) while that of roots separated well 3 patterns corresponding to species D. linearis, P. calomelanos, and P. vittata (Figure 2(b)). This indicates that there is a contribution to the difference in accumulation ability among tissues on the discrimination of 75 samples based on species. This result would be useful for understanding the different translocation factors of 4 ferns calculated in Table 2. In addition, taking into account the variable vector in a PCA biplot, it can be seen that Cu was accumulated with a higher level in all tissue of P. calomelanos which confirms again the previous finding of the Cu accumulation ability of this fern [13, 36, 37]. In the current study, P. calomelanos does not show the arsenic accumulation capability, whereas many studies have reported that this fern can accumulate more than 2000 mg/kg in their tissues [12, 38], even better than P. vittata [39]. However, the arsenic accumulation behavior of these two ferns changed according to field conditions, i.e., subtropical or tropical, highly or lightly contaminated soils, drought and other weeds under field conditions [40]. Other studies also reported the ability of P. calomelanos in As-Pb coaccumulating [15, 41]. In this study, the poor performance of P. calomelanos in HM accumulation at Ha Thuong, Dai Tu may be due to low concentration of As, Pb, Cu in this site, compared to Tan Long, Dong Hy where P. vittata were collected. It can also be highlighted from Figure 3 that As, Pb, Zn were accumulated with the highest levels in P. vittata. This finding is consistent with the previous research conclusion that P. vittata is an arsenic hyperaccumulator [27], useful for the remediation of As polluted sites cocontaminated with Zn [11], Pb [15, 41], and Cd [6].Figure 2
Principal component analysis of HM levels in 75 samples (3a), in 25 roots (3b), in 25 shoots (3c), and in 25 leaves (3d) from four speciesC. parasiticus, D. linearis, P. calomelanos, and P. vittata. Samples are colored and shaped according to different species. The length of the HM vector shows how high the concentration of these HMs levels is in the samples where this vector is directed toward. (a) All tissues, (b) roots, (c) shoots, and (d) leaves.
(a)(b)(c)(d)Figure 3
Boxplots of the DPPH and SRSA scavenging inhibition of leaves, roots, and shoots from 4 studied ferns.
(a)(b)(c)(d)


### 3.4. Relationship between the HMs Levels in Soils and Plants

To evaluate whether the HMs levels in soils have impact on the difference in accumulation activity of plant tissues, environmental factors (i.e., HMs levels in soils) were fitted onto the ordination represented by the principal component analysis. As a result, the significant correlation between HMs levels in soils and samples similarities of the dataset for 75 samples and the subsets of roots, shoots, and leaves, withp values below 0.05, are reported in Table 3. For all 75 samples, the concentration of two elements Zn and Cu in soils showed a significant correlation (r2 = 0.42, p value < 0.001, and r2 = 0.40, p value < 0.001, respectively). Furthermore, concerning the impact of HMs in soils on the sample similarities in only one tissue, these two above-mentioned variables presented more important correlations (r2 = 0.57 and r2 = 0.65 in roots, r2 = 0.59 and r2 = 0.55 in shoots, r2 = 0.71 and r2 = 0.59 in leaves, respectively with all p value < 0.001). This suggested that the concentration of Zn and Cu in soils is most responsible for the discrimination based on species between all samples as well as roots, shoots, and leaves. The concentration of three remaining HMs (As, Cd, Pb) was less contributed to the variation in HMs accumulation of 4 weeds (r2 < 0.4 except the levels of Cd and Pb in shoots).Table 3
The correlation coefficient of 5 environmental factors on the HMs in plant tissues, roots, shoots, and leaves was calculated by envfit function in vegan R-package.
All tissuen = 75Rootn = 25Shootn = 25Leaven = 25r2pr2pr2pr2pCu0.40p<0.0010.65p<0.0010.55p<0.0010.59p<0.001Zn0.42p<0.0010.57p<0.0010.59p<0.0010.71p<0.001As0.180.0020.030.720.370.0080.340.009Cd0.30p<0.0010.390.0070.440.0050.350.009Pb0.100.0360.130.240.450.0030.120.244The correlation between the HMs levels in soils and in tissues was calculated by the Spearman method to reveal the relationship among these variables. As can be seen from Figure4, Pb level in tissues correlated very strongly with As concentration in tissues (r = 0.86, p value < 0.001) which indicates there is a coaccumulation of Pb and As of these ferns. To the best our knowledge, it is the first time these two elements were found to be proportionally uptaken in plant tissues. Moreover, Zn concentration in tissues correlated moderately with Zn level in soils (r = 0.43, p value < 0.001) indicating the Zn accumulation capacity of 4 ferns does not depend very much on the fern species.Figure 4
The significant pair to pair correlation of nine variables including As, Pb, Zn concentration in tissue, Zn concentration on soil, DPPH scavenging activity at 100μg/mL and 500 μg/mL, SRSA scavenging activity at 100 μg/mL and 500 μg/mL. The distribution of each variable is shown on the diagonal. On the bottom of the diagonal, the bivariate scatter plots with a fitted line are displayed. On the top of the diagonal, the value of the correlation plus the significance level is denoted as stars. Each significance level is associated with a symbol: ∗∗∗p value < 0.001, ∗∗p value < 0.01, and ∗p value <0.05. The Pearson method was used to compute the correlation coefficient.


### 3.5. Antioxidant Properties of Plant Extracts

The radical scavenging capacity of the crude plant extracts was determined from the reduction in the optical absorbance due to the scavenging of DPPH and SRSA radicals activity. Our study indicated that the roots possessed a high antioxidant activity, while shoots and leaves expressed the lower DPPH scavenging effects for all four species (Figure3). At 500 μg/mL of concentration, most of the samples were considered to have DPPH radical scavenging activity. At a concentration of 100 μg/mL, roots of C. parasiticus presented the highest DPPH scavenging activity with the inhibition of radicals generation value 88.1 ± 2.61% at the concentration 100 μg/mL (Table 4). To date, there is almost no research about the antioxidant of C. parasiticus except a study on the chemical composition of this fern that reveals the presence of some coumarin-chalcone hybrids in this weed [42]. Concerning SRSA scavenging activity, only the roots of P. vittata were considered to show an inhibitory effect on the O2-generation (SC% of 54.4 ± 5.81%). In this study, two different approaches were applied to measure the scavenging capacity against stable, nonbiological radicals and against specific ROS/RNS. Two levels of screening were applied to see whether the antioxidant properties of these weeds can be explained by different mechanisms related to enzymatic (SRSA assay) and nonenzymatic (DPPH assay) [43]. From the applied point of view, this property opens the way to benefit these ferns as phytosource of natural antioxidants that could be used in the chemical catalyst, food, nutraceuticals, cosmetic, and pharmaceutical industries, thereby restricting the use of synthetic antioxidants.Table 4
DPPH scavenging capacity (SC%) of the plant extracts.
Plant speciesPlant tissuesDPPH SC% at 100μg/mLDPPH SC% at 500μg/mLSRSA SC% at 100μg/mLSRSA SC% at 500μg/mLC. parasiticusRoots88.1 ± 2.6187.9 ± 1.969.65 ± 5.3814.6 ± 5.01Shoots29.2 ± 7.3480.4 ± 13.97.83 ± 1.0312.3 ± 9.53Leaves37.7 ± 10.177.9 ± 6.1011.6 ± 0.006.48 ± 11.1D. linearisRoots61.5 ± 26.181.1 ± 9.795.94 ± 3.2117.0 ± 7.46Shoots64.5 ± 28.976.6 ± 11.15.44 ± 3.0512.8 ± 5.05Leaves31.8 ± 16.066.0 ± 15.13.99 ± 1.7410.1 ± 0.90P. calomelanosRoots44.5 ± 24.269.2 ± 21.026.3 ± 2.5136.0 ± 4.82Shoots26.8 ± 12.965.0 ± 20.210.5 ± 6.4118.8 ± 11.8Leaves17.3 ± 5.0447.3 ± 8.2211.4 ± 5.4122.4 ± 8.39P. vittataRoots41.7 ± 20.083.5 ± 9.5329.9 ± 9.1354.4 ± 5.81Shoots31.0 ± 11.461.9 ± 16.310.8 ± 2.8130.3 ± 6.11Leaves46.2 ± 11.889.3 ± 0.6314.9 ± 6.0424.3 ± 5.67Ascorbic acid—91.3———Catechin———48.7—


### 3.6. Relationship between HM Level in Soil and in Plant Tissue and Their Antioxidant Activity

In order to investigate whether HMs in soils and HMs in tissues can affect the antioxidant activity of plants, we analyzed the relationship between the concentrations of five metals in soil and plant samples and the percentages of DPPH and SRSA scavenging of 25 plants growing in these soils. Pair to pair correlations among 14 variables were calculated by the Spearman method. The correlation coefficient was calculated as the statistical measurement of how strong a linear relationship of the paired data is. Taking into consideration the significant correlations (p value < 0.01) and important correlations that were contributed by antioxidant activities, the simplified pairwise correlations between 5 variables (Pb in tissues, Cd in soils, Zn in tissues, Zn in soils, and As in tissues) and 4 antioxidant variables (SC% DPPH at 100 and 500 μg/mL SC% SRSA at 100 and 500 μg/mL) were shown in Figure 4. Noticeably, the concentration of As in tissues was proportional to both SRSA scavenging activity at 100 and 500 μg/mL (r = 0.62 and r = 0.67, respectively, p<0.001) but not correlated significantly with SC% DPPH. This finding has been consistent with the previous study of Cao et al., 2004, which showed that at low levels of As exposure (0–20 mg/kg), enzymatic antioxidants (which may contribute to SRSA activity) play a more important role in As detoxification than nonenzymatic antioxidants (which may contribute to DPPH activity) in P. vittata [44]. By other research, P. calomelanos also presented a remarkable change in amino acid levels in As-treated plants, including those related to the biosynthesis of enzymatic antioxidants of this weed [9]. Similarly, it was also observed that Pb concentration in tissues correlated significantly and strongly to SC% of SRSA at two studied concentrations (r = 0.66, p value < 0.001 and r = 0.70, p value < 0.001, respectively) while did not correlate to SC% of DPPH. This allows hypothesizing that the excess of Pb accumulation significantly enhanced the synthesis of some enzymatic antioxidants in plants, thus induced the antioxidant activity of plants.


## 3.1. HMs Level in Soil

Five HMs As, Cd, Zn, Pb, and Cu were quantified in 25 soil samples collected at the same positions with 25 plant samples including 5 samples ofC. parasiticus (in Thai Nguyen city), 6 samples for each type of D. linearis and P. calomelanos plants (in Dai Tu district), and 8 samples of P. vittata (in Dong Hy district). According to each sampling station, the sample abbreviation was used interchangeably for the respective soil and plant samples.Table1 shows the highest As, Zn, Pb, and Cd concentrations in P. vittata correspondent soil while the most concentration of Cu was found in the soil growing C. parasiticus plant. According to the Vietnam National regulations on the allowable limits of HMs in the soils [7] where the maximum allowable limits of As, Zn, Pb, Cu, and Cd in industrial soils are 25; 300; 300; 300, and 10 mg/kg, respectively. It was investigated that soil samples collected in Dong Hy district (P. vittata corresponding soil) were contaminated by both As and Pb, whereas soils from Dai Tu was contaminated with Cu. Soils from Cam Gia, Thai Nguyen city, witnessed the Pb pollution. Significantly, there were no observations of Zn and Cd contaminations in all soil samples collected from different sampling sites.Table 1
The average of five HMs (As, Zn, Pb, Cd, and Cu) concentrations (in mg/kg) of soil samples corresponding to 4 plant species.
CorrespondenceAsZnPbCdCuC. parasiticus9.31 ± 4.21108 ± 31.9307 ± 1381.98 ± 0.70171 ± 10.0D. linearis2.54 ± 0.4584.1 ± 19.4206 ± 56.60.81 ± 0.82453 ± 313P. calomelanos2.92 ± 0.9864.3 ± 15.4239 ± 66.90.38 ± 0.16338 ± 102P. vittata53.9 ± 13.9265 ± 111450 ± 2482.34 ± 1.4611.2 ± 8.06It is found that the distribution of HMs in soils was following the different soil origins and plant bioaccumulation and/or absorption [23]. To understand more this phenomenon, the principal component analysis (PCA) revealed the difference among the HM levels of three locations where 25 soil samples were collected by distinguishing P. vittata, D. linearis, P. calomelanos, and C. parasiticus zones (Figure 1). The first two axes of the PCA explained 66.4% of the variance. The ordination arrow represents the dominance of As and Zn in the soil samples collected in the rhizosphere of P. vittata in Dong Hy district. More Pb and Cd than other HMs were distributed in soils sampled in Dong Hy district and Thai Nguyen city, while Cu significantly existed in soils collected in Dai Tu district. These distributions of HMs in soils have thus corresponded to their detected concentration in different soils (Table 1). In addition, the arrow direction of vectors, which represented the distribution tendency of HMs according to soil samples, were also indicated for the cocontamination relationship between HMs in the soils. Cu vectors tend to contradict the vectors of other metals, indicating that Cu has less ability to coaccumulate concurrently with other metals in the soil. This was identified more clearly based on the Spearman index (Figure S1) when the indices representing the correlation between Cu and other metals were always very low (from −0.76 to −0.36). This phenomenon can be explained by higher solubility and leachability of Cu at the ideal soil pH of 6.0–6.5, compared to other HMs [24]. In contrast, Cd appeared to have a stronger correlation in soil accumulation with other metals, except Cu, while the lowest concentrations of Cd were often found in all soil samples (Table 1). Consequently, the Spearman indices indicating the cocontamination between Cd and other HMs following the order Zn (0.76) > As (0.59) > Pb (0.37). This could be partly clarified for the highest concentrations of Zn, As, Pb, and Cd were found in Dong Hy soil as previous investigation (Table 1). Also, Pb had significant cocontamination with As and Zn in soil following Spearman indices of 0.39 and 0.51, respectively, confirming for high abilities of Zn, As, and Pb accumulations in soils corresponding to P. vittata plant. The analysis results, therefore, presented the significant effects of cocontaminations between soil origins and the HM residences in those soils.Figure 1
PCA of the concentration of HMs (As, Cd, Cu, Pb, and Zn) in 25 soil samples corresponding to 25 plants of 4 species (C. parasiticus, D. linearis, P. calomelanos, and P. vittata). Each point displays a sample; each vector displays the level of an HM. Soil samples are colored and shaped according to different species where they were collected.


## 3.2. HMs Levels in Plant Tissues

The statistical description of HM concentration (in mg/kg) in 75 plant tissue samples was summarized by the average concentration with standard deviation. For each HM, the bioaccumulation factor (BF) was calculated by the ratio between the HM concentration in tissue and those in corresponding soil which is described as the ability of plants for elemental accumulation from the substrate. The efficiency of phytoremediation can be quantified by calculating the translocation factor (TF). This index expresses the capacity of a plant to store the HMs in its upper part, which is defined as the ratio of metal concentration in the upper part to that in the roots [25].The total concentration of selected metals as well as BF and TF indices in plant roots, shoots, and leaves is shown in Table2. In plant roots, the highest As, Pb, and Zn concentrations were found in P. vittata, whereas the most Cu concentration was measured in P. calomelanos. The large amounts of As, Pb, Zn, and Cu found in plant roots of P. vittata and P. calomelanos could be explained by the high level of these metal concentrations detected in the soils where corresponding plants were cultivated (Table 1). On the other side, there was no significant difference in Cd concentration in the roots of 4 plants. This correlated to the low level of Cd found in all corresponding soils (Table 1). On average, the HM concentrations in the roots were much higher than those in the leaves and shoots, implying that fewer HMs were translocated into leaves and shoots for all four species.Table 2
The average concentration (±SD) of five HMs in plant tissues (root, shoot, and leave) from 4 species,C. parasiticus, D. linearis, P. calomelanos, and P. vittata, and corresponding bioaccumulation factor (BF) and translocation factor (TF).
Plant speciesPlant tissueAsCdCuPbZnConc.BFTFConc.BFTFConc.BFTFConc.BFTFConc.BFTFC. parasiticusRoot1.19 ± 1.770.14—0.58 ± 0.490.29—19.7 ± 5.920.12—4.34 ± 1.480.01—121 ± 73.61.12—Shoot1.19 ± 0.900.214.510.66 ± 0.300.3618.413.1 ± 11.10.080.882.79 ± 1.640.011.2828.0 ± 13.00.271.10Leave1.35 ± 1.580.184.430.77 ± 0.300.458.3715.2 ± 8.800.090.815.51 ± 1.950.020.61110 ± 50.71.100.28D. linearisRoot122 ± 32.250.0—0.19 ± 0.090.47—75.8 ± 53.00.47—37.3 ± 20.90.21—35.5 ± 10.80.45—Shoot18.3 ± 14.48.030.320.08 ± 0.050.247.6141.1 ± 22.20.240.388.85 ± 4.490.050.4722.1 ± 10.70.290.99Leave33.4 ± 30.514.90.570.32 ± 0.250.820.4428.4 ± 14.20.170.5915.4 ± 6.220.090.2543.7 ± 31.50.570.59P. calomelanosRoot469 ± 69.5169—0.32 ± 0.110.97—430 ± 75.01.37—160 ± 68.80.78—92.3 ± 19.01.52—Shoot227 ± 71.81020.610.08 ± 0.030.260.2950.1 ± 15.30.160.1213.6 ± 10.90.070.1024.0 ± 2.260.400.27Leave630 ± 27.52291.370.16 ± 0.050.480.5759.0 ± 79.00.180.1716.6 ± 4.360.080.1249.8 ± 9.260.830.55P. vittataRoot1732 ± 299.1403—0.27 ± 0.260.13—5.98 ± 1.301.06—244 ± 60.70.88—188 ± 30.50.35—Shoot331 ± 75.376.20.200.13 ± 0.060.080.845.09 ± 2.370.960.8465.4 ± 18.30.190.2854.2 ± 22.20.230.65Leave438 ± 17699.51.090.33 ± 0.280.181.496.13 ± 3.41.171.03141 ± 51.50.470.59267 ± 1251.133.18In plant shoots, the As concentration ranged from 1.19 ± 0.9 mg/kg (inC. parasiticus) to 331 ± 75.3 mg/kg (in P. vittata). The total concentration (mg/kg) of Cu, Pb, and Zn in plant shoots varied significantly, with the highest values being 59.1 ± 15.3 (in P. calomelanos), 65.4 ± 18.3 (in P. vittata), and 54.2 ± 22.2 (in P. vittata), respectively. Similarly, in the case of metals accumulation in leaves, higher metals concentrations of As, Cu, Pb, and Zn were found in plant leaves of P. vittata and P. calomelanos. On the other side, according to different plant types, a similar distribution tendency of measured Cd concentrations was found in plant roots and shoots, and the lower metals concentrations were often found in all tissues of C. parasiticus and D. linearis compared to other ones. It could be seen that the distribution of metal concentration in plant tissues varied greatly among plant species and had a significant correlation with HM content in soils from different sampling sites.Moreover, the bioaccumulation factor (BF) indicates the ability of metals stored in plant tissues compared to soil [26]. For all plant tissues, the BF values often lower than 1.0 for Cd, Cu, Pb and Zn. In contrast, the BF values of As which were extremely high in all tissues of almost all plants except C. parasiticus with BFAs lower than 1. This indicated the higher moving ability of As from the external environment into plant tissues through the soil exposure route compared to the other 3 metals. It is also found that P. vittata sampled at Dong Hy possessed a higher metal bioaccumulation factor (up to 402.5) than other plants indicating its hyperaccumulating ability. There are a number of reports indicating that P. vittata is an As hyperaccumulator [27, 28] capable of translocating As and storing it in the shoots [29, 30].Similarly, most shoot and root tissues of all almost plants had TF values lower than 1, except TFAs (inC. parasiticus) and TFCd (in C. parasiticus and D. linearis plant) values, which are much higher than 1. In more detail, the maximum values of TFCd and TFAs found in shoots of C. parasiticus plants were 18.4 and 4.51, respectively. Furthermore, the TFmetals indices in the shoots were often higher than that in leaves of C. parasiticus and D. linearis plants, whereas higher TF metals indices were found in leaves of P. calomelanos and P. vittata. High root to shoot translocation of these metals in C. parasiticus and D. linearis indicated vital characteristics to be used in phytoextraction of these metals in shoot [31, 32]. On the other side, the results in TF and BF values of P. calomelanos and P. vittata were considered suitable for phytoextraction, which generally requires translocation of metals easily harvestable leaves tissue to remove the metal contamination from the soil without destroying soil structure and fertility [33, 34]. In fact, plant species were judged as being suited for phytoremediation based on several criteria including wide distribution, high above-ground biomass, high bioaccumulation factors (hyperaccumulators), short life cycles, and high propagation rates [35]. Therefore, it could be concluded that P. vittata and P. calomelanos can be used in phytoremediation in mine sites of Thai Nguyen province as they were considered as native plants that meet all criteria for a hyperaccumulator.


## 3.3. Species-Tissue Discrimination of HMs Levels in Plants

In order to see whether there are differences in the concentration of HMs level (i.e., variables) in samples (i.e., individuals), PCA biplots where both individuals and variables are graphed on the same plot were constructed. In the first two PCs, more than 60% of the variance in the data is explained for all PCA plots. In accordance with the PCA biplot (Figure2(a)), HM levels in 75 plant tissues were grouped in a different pattern based on species. The discrimination of these patterns indicates that four species show a different capacity in uptaking 5 mentioned HMs. The PCA biplots containing only data of HM levels in shoots or leaves showed better separation between four species (Figures 2(c) and 2(d)) while that of roots separated well 3 patterns corresponding to species D. linearis, P. calomelanos, and P. vittata (Figure 2(b)). This indicates that there is a contribution to the difference in accumulation ability among tissues on the discrimination of 75 samples based on species. This result would be useful for understanding the different translocation factors of 4 ferns calculated in Table 2. In addition, taking into account the variable vector in a PCA biplot, it can be seen that Cu was accumulated with a higher level in all tissue of P. calomelanos which confirms again the previous finding of the Cu accumulation ability of this fern [13, 36, 37]. In the current study, P. calomelanos does not show the arsenic accumulation capability, whereas many studies have reported that this fern can accumulate more than 2000 mg/kg in their tissues [12, 38], even better than P. vittata [39]. However, the arsenic accumulation behavior of these two ferns changed according to field conditions, i.e., subtropical or tropical, highly or lightly contaminated soils, drought and other weeds under field conditions [40]. Other studies also reported the ability of P. calomelanos in As-Pb coaccumulating [15, 41]. In this study, the poor performance of P. calomelanos in HM accumulation at Ha Thuong, Dai Tu may be due to low concentration of As, Pb, Cu in this site, compared to Tan Long, Dong Hy where P. vittata were collected. It can also be highlighted from Figure 3 that As, Pb, Zn were accumulated with the highest levels in P. vittata. This finding is consistent with the previous research conclusion that P. vittata is an arsenic hyperaccumulator [27], useful for the remediation of As polluted sites cocontaminated with Zn [11], Pb [15, 41], and Cd [6].Figure 2
Principal component analysis of HM levels in 75 samples (3a), in 25 roots (3b), in 25 shoots (3c), and in 25 leaves (3d) from four speciesC. parasiticus, D. linearis, P. calomelanos, and P. vittata. Samples are colored and shaped according to different species. The length of the HM vector shows how high the concentration of these HMs levels is in the samples where this vector is directed toward. (a) All tissues, (b) roots, (c) shoots, and (d) leaves.
(a)(b)(c)(d)Figure 3
Boxplots of the DPPH and SRSA scavenging inhibition of leaves, roots, and shoots from 4 studied ferns.
(a)(b)(c)(d)


## 3.4. Relationship between the HMs Levels in Soils and Plants

To evaluate whether the HMs levels in soils have impact on the difference in accumulation activity of plant tissues, environmental factors (i.e., HMs levels in soils) were fitted onto the ordination represented by the principal component analysis. As a result, the significant correlation between HMs levels in soils and samples similarities of the dataset for 75 samples and the subsets of roots, shoots, and leaves, withp values below 0.05, are reported in Table 3. For all 75 samples, the concentration of two elements Zn and Cu in soils showed a significant correlation (r2 = 0.42, p value < 0.001, and r2 = 0.40, p value < 0.001, respectively). Furthermore, concerning the impact of HMs in soils on the sample similarities in only one tissue, these two above-mentioned variables presented more important correlations (r2 = 0.57 and r2 = 0.65 in roots, r2 = 0.59 and r2 = 0.55 in shoots, r2 = 0.71 and r2 = 0.59 in leaves, respectively with all p value < 0.001). This suggested that the concentration of Zn and Cu in soils is most responsible for the discrimination based on species between all samples as well as roots, shoots, and leaves. The concentration of three remaining HMs (As, Cd, Pb) was less contributed to the variation in HMs accumulation of 4 weeds (r2 < 0.4 except the levels of Cd and Pb in shoots).Table 3
The correlation coefficient of 5 environmental factors on the HMs in plant tissues, roots, shoots, and leaves was calculated by envfit function in vegan R-package.
All tissuen = 75Rootn = 25Shootn = 25Leaven = 25r2pr2pr2pr2pCu0.40p<0.0010.65p<0.0010.55p<0.0010.59p<0.001Zn0.42p<0.0010.57p<0.0010.59p<0.0010.71p<0.001As0.180.0020.030.720.370.0080.340.009Cd0.30p<0.0010.390.0070.440.0050.350.009Pb0.100.0360.130.240.450.0030.120.244The correlation between the HMs levels in soils and in tissues was calculated by the Spearman method to reveal the relationship among these variables. As can be seen from Figure4, Pb level in tissues correlated very strongly with As concentration in tissues (r = 0.86, p value < 0.001) which indicates there is a coaccumulation of Pb and As of these ferns. To the best our knowledge, it is the first time these two elements were found to be proportionally uptaken in plant tissues. Moreover, Zn concentration in tissues correlated moderately with Zn level in soils (r = 0.43, p value < 0.001) indicating the Zn accumulation capacity of 4 ferns does not depend very much on the fern species.Figure 4
The significant pair to pair correlation of nine variables including As, Pb, Zn concentration in tissue, Zn concentration on soil, DPPH scavenging activity at 100μg/mL and 500 μg/mL, SRSA scavenging activity at 100 μg/mL and 500 μg/mL. The distribution of each variable is shown on the diagonal. On the bottom of the diagonal, the bivariate scatter plots with a fitted line are displayed. On the top of the diagonal, the value of the correlation plus the significance level is denoted as stars. Each significance level is associated with a symbol: ∗∗∗p value < 0.001, ∗∗p value < 0.01, and ∗p value <0.05. The Pearson method was used to compute the correlation coefficient.


## 3.5. Antioxidant Properties of Plant Extracts

The radical scavenging capacity of the crude plant extracts was determined from the reduction in the optical absorbance due to the scavenging of DPPH and SRSA radicals activity. Our study indicated that the roots possessed a high antioxidant activity, while shoots and leaves expressed the lower DPPH scavenging effects for all four species (Figure3). At 500 μg/mL of concentration, most of the samples were considered to have DPPH radical scavenging activity. At a concentration of 100 μg/mL, roots of C. parasiticus presented the highest DPPH scavenging activity with the inhibition of radicals generation value 88.1 ± 2.61% at the concentration 100 μg/mL (Table 4). To date, there is almost no research about the antioxidant of C. parasiticus except a study on the chemical composition of this fern that reveals the presence of some coumarin-chalcone hybrids in this weed [42]. Concerning SRSA scavenging activity, only the roots of P. vittata were considered to show an inhibitory effect on the O2-generation (SC% of 54.4 ± 5.81%). In this study, two different approaches were applied to measure the scavenging capacity against stable, nonbiological radicals and against specific ROS/RNS. Two levels of screening were applied to see whether the antioxidant properties of these weeds can be explained by different mechanisms related to enzymatic (SRSA assay) and nonenzymatic (DPPH assay) [43]. From the applied point of view, this property opens the way to benefit these ferns as phytosource of natural antioxidants that could be used in the chemical catalyst, food, nutraceuticals, cosmetic, and pharmaceutical industries, thereby restricting the use of synthetic antioxidants.Table 4
DPPH scavenging capacity (SC%) of the plant extracts.
Plant speciesPlant tissuesDPPH SC% at 100μg/mLDPPH SC% at 500μg/mLSRSA SC% at 100μg/mLSRSA SC% at 500μg/mLC. parasiticusRoots88.1 ± 2.6187.9 ± 1.969.65 ± 5.3814.6 ± 5.01Shoots29.2 ± 7.3480.4 ± 13.97.83 ± 1.0312.3 ± 9.53Leaves37.7 ± 10.177.9 ± 6.1011.6 ± 0.006.48 ± 11.1D. linearisRoots61.5 ± 26.181.1 ± 9.795.94 ± 3.2117.0 ± 7.46Shoots64.5 ± 28.976.6 ± 11.15.44 ± 3.0512.8 ± 5.05Leaves31.8 ± 16.066.0 ± 15.13.99 ± 1.7410.1 ± 0.90P. calomelanosRoots44.5 ± 24.269.2 ± 21.026.3 ± 2.5136.0 ± 4.82Shoots26.8 ± 12.965.0 ± 20.210.5 ± 6.4118.8 ± 11.8Leaves17.3 ± 5.0447.3 ± 8.2211.4 ± 5.4122.4 ± 8.39P. vittataRoots41.7 ± 20.083.5 ± 9.5329.9 ± 9.1354.4 ± 5.81Shoots31.0 ± 11.461.9 ± 16.310.8 ± 2.8130.3 ± 6.11Leaves46.2 ± 11.889.3 ± 0.6314.9 ± 6.0424.3 ± 5.67Ascorbic acid—91.3———Catechin———48.7—


## 3.6. Relationship between HM Level in Soil and in Plant Tissue and Their Antioxidant Activity

In order to investigate whether HMs in soils and HMs in tissues can affect the antioxidant activity of plants, we analyzed the relationship between the concentrations of five metals in soil and plant samples and the percentages of DPPH and SRSA scavenging of 25 plants growing in these soils. Pair to pair correlations among 14 variables were calculated by the Spearman method. The correlation coefficient was calculated as the statistical measurement of how strong a linear relationship of the paired data is. Taking into consideration the significant correlations (p value < 0.01) and important correlations that were contributed by antioxidant activities, the simplified pairwise correlations between 5 variables (Pb in tissues, Cd in soils, Zn in tissues, Zn in soils, and As in tissues) and 4 antioxidant variables (SC% DPPH at 100 and 500 μg/mL SC% SRSA at 100 and 500 μg/mL) were shown in Figure 4. Noticeably, the concentration of As in tissues was proportional to both SRSA scavenging activity at 100 and 500 μg/mL (r = 0.62 and r = 0.67, respectively, p<0.001) but not correlated significantly with SC% DPPH. This finding has been consistent with the previous study of Cao et al., 2004, which showed that at low levels of As exposure (0–20 mg/kg), enzymatic antioxidants (which may contribute to SRSA activity) play a more important role in As detoxification than nonenzymatic antioxidants (which may contribute to DPPH activity) in P. vittata [44]. By other research, P. calomelanos also presented a remarkable change in amino acid levels in As-treated plants, including those related to the biosynthesis of enzymatic antioxidants of this weed [9]. Similarly, it was also observed that Pb concentration in tissues correlated significantly and strongly to SC% of SRSA at two studied concentrations (r = 0.66, p value < 0.001 and r = 0.70, p value < 0.001, respectively) while did not correlate to SC% of DPPH. This allows hypothesizing that the excess of Pb accumulation significantly enhanced the synthesis of some enzymatic antioxidants in plants, thus induced the antioxidant activity of plants.


## 4. Conclusion

This study firstly showed the investigation of HM concentrations in soils that had a significant variation according to mining soils and corresponding plant types. As accumulated at high levels in soils and plant tissues (shoot, roots, leaves). Cd had a strong cocontamination relationship with other metals As, Zn, Pb according to high Spearman values (from 0.37 to 0.76). Cu contaminated in Dai Tu district and had a low ability to cocontaminated with four other HMs. For all weeds, with BF and TF values, almost lower than 1, Cu, Zn, Pb, and Cd had limited ability to uptake in plant tissues. In addition, the distribution or discrimination of HM followed the different species, tissues, and soil sampling regions. Indeed, it is found that the roots ofP. vittata of C. parasiticus possess the antioxidant activity via DPPH and SRSA radical scavenging assay, respectively. Therefore, the correlation between antioxidant activity of plants and HMs content has opened a prospect in the identification of hyperaccumulator species with high potential in pharmaceutical as well as phytoremediation research.



---

*Source: 8010376-2020-10-31.xml*