Accumulation and fractionation of rare earth elements are conserved traits in the Phytolacca genus

Rare earth elements (REEs) are now considered emerging pollutants in the environment. Phytolacca americana, an REE hyperaccumulating plant, has been proposed for the remediation of REE-contaminated soils. However, there is no REE-related information for other Phytolacca species. Here, we examined five species (P. americana, P. acinosa, P. clavigera, P. bogotensis, and P. icosandra) for their response to REEs. REE accumulation and fractionation traits both occurred on the same order of magnitude among the five species. Heavy REEs were preferentially transferred to leaves relative to light REEs. Regardless of the species, lateral root length and chlorophyll content decreased under REE exposure, and lateral roots and foliar anthocyanins increased. However, plants did not experience or only slightly experienced oxidative stress. Finally, REE exposure strongly modulated the ionome of roots and, to a lesser extent, that of leaves, with a negative correlation between REE and Mn contents. In conclusion, our study provides new data on the response of several Phytolacca species to REEs. Moreover, we highlighted that the REE accumulation trait was conserved among Phytolacca species. Thus, we provide valuable information for the phytoremediation of REE-contaminated sites since the most appropriate Phytolacca species could be selected depending on the climatic/pedological area to be remediated.


Results and Discussion
Molecular confirmation of the Phytolacca species. Seeds from five species of Phytolacca were obtained from commercial or institutional origins as described in the materials and methods section and identified as P. americana, P. acinosa, P. bogotensis, P. clavigera, and P. icosandra. The identification of these species was based on morphological traits. However, due to the ambiguous taxonomy within the Phytolacca genus 28 , we performed a molecular investigation to confirm that the species tested were different. For this purpose, the ITS region was amplified, sequenced and compared with ITS sequences of Phytolacca species retrieved from GenBank. For most species, there was a single ITS sequence available in the database, and none were available for P. clavigera (Fig. 1). For P. bogotensis, P. icosandra and P. americana, our data matched the previously published ITS sequences for these species. However, it was less straightforward for P. acinosa and P. clavigera, where a mis-identification might have been made. Facing this problem, we analysed a second seed batch of P. acinosa obtained from a different supplier. Both plant morphology and ITS sequences were identical for the two batches, suggesting a correct identification of P. acinosa. Notably, different leaf morphologies were observed between P. acinosa and P. clavigera (Fig. 1), which shared the highest identity (97%). Leaves of P. acinosa were obovate, whereas the four other species had oval-shaped leaves. In conclusion, the molecular analyses confirmed that the five species tested were different. Sequence data were deposited in GenBank under the following accession numbers: MK602340 (P. americana), MK602343 (P. acinosa), MK602344 (P. bogotensis), MK602341 (P. clavigera), and MK602342 (P. icosandra).
the Ree accumulation trait was shared among Phytolacca species. We first investigated the potential of the different Phytolacca species to accumulate REEs in their tissues. To this end, Phytolacca seedlings were grown hydroponically and exposed to a mixture of 16 REEs supplied at equimolar concentrations. Diatloff et al. (1996) reported that in uncontaminated sites, REE concentrations from the soil solution can reach the micromolar range 29 . Colim et al. (2019) analyzed the REE concentrations in surface water samples from the Lavras do Sul (Brazil) mining region. The total REE concentrations measured ranged from 27 to 279 µM 30 . Therefore, we decided to use two different REE concentrations, a low and a high concentration (10 and 100 µM). These concentrations were also pre-determined from preliminary experiments (data not shown), for which biomass was either not affected or reduced. Indeed, at the concentration of 10 µM REEs in the nutrient solution (REE10), no negative effect was observed for root, shoot and total plant biomass for any species (Fig. 2). A significant positive effect on the shoot biomass of P. icosandra and the root biomass of P. americana was observed even at this lowest concentration. Conversely, at 100 µM REEs (REE100), biomass tended to decrease by approximately 50% for most species, despite not always being significant (Fig. 2). These results suggest a similar REE tolerance level of the five Phytolacca species.
Then, REE accumulation was determined by ICP-MS, and the results are given in Fig. 3. At REE10, the five species accumulated approximately 1500 mg/kg REEs in the roots, and no significant difference was observed between the species (Fig. 3a). However, at REE100, an overall higher accumulation was found, and differences between species could be noted (Fig. 3a). Phytolacca icosandra displayed the highest concentration of approximately 13,000 mg/kg REEs, while P. acinosa and P. bogotensis accumulated 2.4 and 2.3 times less REEs, respectively (Fig. 3a). P. clavigera and P. americana had intermediate concentrations in their roots.
In leaves, the patterns of REE accumulation differed from those of roots. At REE10, P. bogotensis and P. acinosa were the two species accumulating the highest concentrations of REEs, with 279 and 261 mg/kg, respectively. Phytolacca americana and Phytolacca clavigera accumulated slightly less, with 174 and 157 mg/kg, respectively, while the lowest concentration was recorded for P. icosandra, with less than 100 mg/kg (Fig. 3a). However, those variations were not observed at REE100. Despite a mean accumulation value ranging from 265 mg/kg (P. clavigera) to 837 mg/kg (P. bogotensis), these values were not significantly different (Fig. 3a). Our results with P. americana are consistent with those from a previous study, carried out with the same species also grown hydroponically, where a similar REE accumulation capacity was reported 21 . As a comparison, P. americana accumulated in the leaves at approximately 300 and 500 mg/kg REEs when exposed to 10 and 100 µM REEs, respectively 21 .
For each species, the translocation factor (TF) of REEs from roots to leaves was further calculated (Fig. 3b). Under both REE concentrations, P. acinosa and P. bogotensis had the highest TFs, whereas P. clavigera and P. icosandra had the lowest TFs. However, the TF values were in the same range, and relatively low differences were found between the species (Fig. 3b). Moreover, our data are in agreement with those described by Yuan et al. (2017) for P. americana, where a TF of approximately 0.1 was reported 21 .
Overall, our results suggest that the accumulation of REEs in the leaves and their translocation from roots to shoots are of the same order of magnitude in the five Phytolacca species.
HREEs were preferentially translocated to leaves in the five Phytolacca species. We further investigated whether various REE fractionation patterns (LREEs vs HREEs) could be mediated by the different Phytolacca species. Therefore, we analysed the fractionation pattern of the whole set of REEs in both roots and shoots of the five species. To compare the fractionation process for various REE species mediated by different plant species, we used our hydroponic system to ensure equimolar concentrations of all REE species in the nutritive medium. The data of REE fractionation are given in Fig. 4a. REEs are ordered by their decreasing ionic radii at their trivalent oxidation state. The non-lanthanide Y, generally included in the HREEs, thus ranged between Dy and Ho. Interestingly, the five Phytolacca species displayed very similar REE accumulation patterns. In roots, the REE concentration decreased from La to Lu (Fig. 4a). Conversely, in leaves, the REE concentration increased from La to Lu (Fig. 4a). These conclusions were true regardless of the Phytolacca species and the REE concentration (REE10, REE100). The only exceptions were praseodynium (Pr) and ytterbium (Yb), for which anomalies could be demonstrated. However, these anomalies were similar among the five Phytolacca species.   www.nature.com/scientificreports www.nature.com/scientificreports/ We further analysed the accumulation ratios between LREEs and HREEs ( Fig. 4b). At REE10, there was no particular enrichment of LREEs versus HREEs in roots. At REE100, and except for P. icosandra, a slight LREE enrichment in roots was found for the different Phytolacca species. Conversely, leaves were highly enriched with HREEs regardless of the Phytolacca species considered. At REE10, HREE concentration was indeed between 2 and 4 times higher than that of LREEs. This enrichment was more pronounced at REE100, where HREEs were from 5 (P. bogotensis) to 13 (P. clavigera) times more abundant than LREEs (Fig. 4b). Altogether, the data indicated a preferential root-to-leaf translocation of HREEs vs LREEs, a trait that was conserved in the five Phytolacca species.
The preferential translocation of HREEs to shoots had already been observed for P. americana grown either under hydroponic conditions or naturally in REE-mining areas 7,21 . In non-accumulating species, a higher transfer of HREEs vs LREEs has also been reported in several angiosperms, notably wheat 31,32 , soybean 33 and rice 34 . Conversely, LREE enrichment in shoots has been reported in fern species accumulating REEs. Likewise, non-accumulating ferns also had a higher content of LREEs than HREEs (Grosjean N., personal communication). Therefore, the different fractionation patterns described above are very unlikely to be related to the REE accumulation trait. More likely, angiosperms and pteridophytes, which are evolutionarily very distant, might not share the same mechanisms underlying REE translocation from roots to shoots. Such different fractionation processes could be explained by the production of specific compounds with distinct REE-chelating properties. This hypothesis was supported by previous works that reported that the REE-accumulating fern D. dichotoma produces a specific LREE-binding peptide 35 and that HREE enrichment in the shoots of P. americana was associated and total plants (c) exposed to 10 µM or 100 µM REEs (REE10 and REE100, respectively). The data are the means (±SE) of n = 3 (control, REE10) or n = 4 (REE100) plants. Significant differences from the control condition are indicated by asterisks (P < 0.05, t-test).

REE (ppm)
Leaves REE (ppm)  www.nature.com/scientificreports www.nature.com/scientificreports/ with the long-distance transport of HREE-organic acid complexes, such as HREE-citrate 21 . However, to test this hypothesis, further studies investigating the ligands associated with REEs in the xylem sap of a panel of angiosperms and pteridophytes are needed.

Roots
Rees impacted root architecture. It has already been reported that high concentrations of La can inhibit root elongation and induce lateral root development in Arabidopsis thaliana 36 . This mechanism was explained by the accumulation of reactive oxygen species (ROS) in the root tip, leading to the death of cells from the primary root tips and reorientation of auxin to lateral roots 36,37 . Thus, to verify whether root elongation and branching were also affected by a mixture of REEs, the root architecture of the five Phytolacca species was analysed after exposure to 10 and 100 µM REEs and compared to control plants (Fig. 5). The overall root architecture of the  Roots -REE100 30  and (c) lateral root length of the different Phytolacca species exposed to a mixture of 10 µM or 100 µM REEs (REE10 and REE100, respectively) or not exposed to REEs (control). Root systems of three (control, REE10) or four (REE100) plants per species and per treatment were analysed. Within a given species, significant differences between treatments are indicated by different letters (P < 0.05, Kruskal-Wallis www.nature.com/scientificreports www.nature.com/scientificreports/ different species was relatively similar under the control condition (Fig. 5a). The addition of REEs triggered modifications in the root architecture that were characterized by shorter lateral roots compared to the control condition ( Fig. 5a). At REE10, morphological differences were noticed for P. clavigera, P. icosandra, and P. acinosa when compared to non-exposed plants. These differences were much less obvious in P. bogotensis and P. americana. However, at REE100, the root architecture of all species was more impacted, even if P. americana was the less affected species (Fig. 5a). The quantification of both the density and length of lateral roots confirmed the visual observations (Fig. 5b,c). The number of lateral roots indeed increased under REE exposure. However, for both P. americana and P. bogotensis, the values were only significantly different when comparing the control and REE100 conditions (Fig. 5b). Similarly, with increasing REE concentrations, the lateral root length decreased, with the exception of P. americana, which was not affected at REE10 when compared to the control (Fig. 5c). In conclusion, analysis of the root architecture of the different species allowed us to show that, despite a similar REE accumulation rate for the species tested, P. americana was the least affected by REEs.
Phytolacca species did not experience or only slightly experienced oxidative stress under REE exposure. We further selected two markers related to oxidative stress, namely, the malondialdehyde (MDA) concentration and the total antioxidant capacity (TAC), to investigate the effect of REE exposure in the five species (Table 1).
TAC represents the non-enzymatic antioxidant capacity and is indicative of the ability to counteract oxidative stress-induced damage in cells. In both roots and leaves of the different species, only a few differences were found among the different treatments (Table 1). Indeed, REE exposure only impacted the TAC (1.3-fold compared to the control) in the leaves of P. clavigera. Similarly, in the roots of P. acinosa, the TAC content increased by 1.6-fold at REE10 when compared to the control ( Table 1).
The quantification of MDA, one of the end products of the chain reaction of lipid peroxidation, allows us to estimate potential oxidative damage. No significant change was observed in roots under REE exposure, regardless of the species considered. In the leaves of P. acinosa and P. bogotensis, the MDA concentration decreased by 2.5-fold at REE10 and to a lesser extent (1.3-fold) in P. acinosa at REE100. Conversely, it slightly increased (1.6-fold) in the leaves of P. icosandra at REE10 (Table 1).
Combined, MDA and TAC analyses suggest that Phytolacca species were not affected or were only slightly affected by oxidative stress when exposed to REEs. Conversely, several studies have demonstrated the generation of REE-induced oxidative stress in plant species that are not REE accumulators (e.g., Nymphoides peltata and Hydrilla verticillata), with, for example, a progressive increase in the concentration of MDA with increasing REE concentrations [38][39][40] . Therefore, our data suggest that the five Phytolacca species are most likely highly tolerant to REEs. However, additional studies will be needed to compare the REE tolerance levels of different Phytolacca species.
Rees impacted the pigment contents and nitrogen balance index. Given the high REE accumulation in leaves, we investigated its impact on the content of different pigments (chlorophyll, flavonoids, and anthocyanins) and on the nitrogen balance index (NBI). Differences between the treatments but also between the species were recorded (Fig. 6). Regardless of the species, the chlorophyll index moderately decreased at REE100 when compared to the control. Indeed, the chlorophyll content was reduced from 1.1 to 1.6 times in P. bogotensis and P. icosandra, respectively (Fig. 6a). The lower chlorophyll content under REE exposure could either be due to a reduction of Mg supply 40 or to chlorophyll degradation caused by lipid peroxidation 39 . Since the MDA concentration did not increase at REE100 in the different Phytolacca species, this second hypothesis is very unlikely. However, it is noteworthy that the chlorophyll content increased for P. icosandra and P. americana at REE10 www.nature.com/scientificreports www.nature.com/scientificreports/ (Fig. 6a). Such a result could suggest a stimulation of photosynthesis. Several studies reported growth stimulation in rice 41 , tobacco 42 and soybean 43 at low REE concentrations.
Compared to chlorophylls, the flavonoid content was poorly affected by REE exposure (Fig. 6b). When compared to the control, it slightly increased at REE10 by 1.1 and 1.2 times in P. icosandra and P. acinosa, respectively. Conversely, it decreased in P. americana and P. bogotensis at either REE10 or REE100. The anthocyanin index, usually indicative of plant stress 44 , was also analysed and is reported in Fig. 6c. A similar pattern was found for the different species, where all species showed a significant increase in anthocyanins at REE100 when compared to the control (Fig. 6c). These data suggest that this high dose of REEs, associated with a subsequent high REE accumulation in leaves (Fig. 3a), could lead to a modification of anthocyanin biosynthesis. In addition, the increase in anthocyanin content could be caused by an increased limitation of phosphate in the plant due to the high reactivity of REEs with phosphates. Surprisingly, at REE10, P. americana was the only species for which there was a decrease of anthocyanins (Fig. 6c), suggesting that at this low REE exposure, the leaf stress was less important to P. americana than to its close relatives. The different pigment indexes measured can provide indications of the health or stress status of plants under tested conditions. It has been demonstrated in Arabidopsis thaliana that an increase of anthocyanins combined with a decrease of the chlorophyll content was triggered in response to metallic stress 45 . Indeed, anthocyanins can play a major protective role against metal stresses by acting as antioxidants 45 . However, since no difference was obtained from the MDA and TCA measurements, strong oxidative stress is unlikely, and a higher anthocyanin index of all five species measured at the highest REE concentration might reflect phosphorus deficiency, as suggested above and as previously reported 46 .
We also analysed the NBI in plants exposed to REEs. NBI is a marker that corresponds to the ratio of chlorophyll content to that of flavonoids and gives an indication of the nitrogen status of plants 47 . While no difference was noticed for P. bogotensis, the NBI slightly increased at REE10 for P. americana, P. clavigera, and P. icosandra (Fig. 6d). Except for P. bogotensis, the NBI decreased at REE100 in the four other species.
The relatively lower root branching, the increase of the chlorophyll content, along with the increased NBI for P. americana at REE10, could suggest that P. americana is slightly more tolerant to low REE exposure compared to the other species tested (Figs. 5, 6). However, since the growth of other Phytolacca species was not more affected by REEs than that of P. americana for both REE10 and REE100 (Fig. 2), REE tolerance levels of the different species should be similar. Further dose-response studies would, however, be needed to shed more light on the tolerance levels of the different Phytolacca species.
Rees modulate the ionome of Phytolacca species. Finally, the ionome of the different species under the three different exposure conditions was established (Supplementary Tables S1, S2). The two different compartments, leaves (Supplementary Tables S1a, S2a) and roots (Supplementary Tables S1b, S2b), were treated separately. Regarding the accumulation and fractionation of REEs, no differences were observed between the five species tested. Consequently, a global analysis was carried out using the ionomes of the different Phytolacca species to identify putative correlations between REE accumulation and essential element composition (Fig. 7). Two principal component analyses (PCA) (Fig. 7a,b) and correlation matrices were generated for both leaves and roots (Fig. 7c,d).
In the PCA for the leaves, the two first principal components accounted for 78.4% of the observed variations, with the first principal component accounting for 54.6% (Fig. 7a). The different conditions tested (control, REE10, and REE100) were well distinguished in the PCA. As expected, all the REEs were positively correlated together (Fig. 7a). Interestingly, in the second principal component analysis, we observed that LREEs and HREEs were grouped separately (Fig. 7a). This is in agreement with the HREE enrichment relative to the LREEs observed in www.nature.com/scientificreports www.nature.com/scientificreports/ the leaves (Fig. 4a). The first axis of the PCA was mostly structured by the REEs, contributing more than 82% each but also by P with 57% and Mn for 40%. In the second axis, S and Mg were the most structuring variables, with 91% of explained inertia, followed by Cu (85%), Ca (84%), B (80%), K (70%), Fe and Mn (65%). The segregation of REEs (first axis) from the other elements (second axis) on the PCA (Fig. 7a) is supported by the Pearson correlation matrix (Fig. 7c). Indeed, relatively few correlations were found between REE accumulation and the concentration of essential elements. Whereas no significant correlation was found between any of the REEs and Fe or Na, negative correlations were found between the heaviest REEs (Tm to Lu) and Cu, Ca, Mg, and Zn (Fig. 7c). A similar negative correlation was observed for Mn with HREEs (Gd excepted) (Fig. 7c). Conversely, positive correlations were found between P and all REEs except La. Similar positive correlations were also found between some LREEs and B and S (Fig. 7c). Interestingly, the Pearson's r coefficient increased (P only) or decreased (the other elements), with the atomic number of REEs (Fig. 7c). Again, Sc did not follow the global rule observed for the other REEs, supporting its difference compared to the lanthanides and Y.
Correlations were more abundant and stronger between macro-or micronutrient elements and REEs in the root compartment (Fig. 7b-d). In the PCA, the two first principal components accounted for 80.5% of the observed variations, with the first principal component accounting for 71.2% (Fig. 7b). The first axis of the PCA was mostly structured by REEs, K, and P, contributing more than 84% each, and Mg contributed 68%. In the second axis, Ca was the most structuring variable, with 74% of explained inertia, followed by Na (65%), Mg (59%), B (56%), and S (47%). REEs covariated together and were also positively correlated with P in both roots and leaves. However, since (i) the nutritive solution lacked P during the exposure phase of REEs and (ii) growth at REE100 was reduced compared to the control, it was expected that P concentration in the biomass would be subjected to a dilution effect in the control treatment compared to the REE-exposed treatments. The observations reported in Supplementary Table S1 and Fig. 7c,d are in agreement with this hypothesis. While no correlation was found in the leaves between Fe or Na and REEs, positive correlations were observed in the roots (Fig. 7b-d). The other elements, namely, Cu, Mn, S, Mg, Ca, Zn, K, and B, were all negatively correlated with REEs. These correlations were highly influenced by the highest exposure condition (REE100) (Fig. 7b-d, Supplementary Tables S1, S2).
REEs have a similar ionic radius to that of Ca. This characteristic has resulted in the use of REEs as Ca-channel blockers 48,49 , and later on, it was suggested that La and Eu could enter into plant cells through Ca-channels 50 . Therefore, the negative correlation found between Ca and REEs is not unexpected. Such antagonistic effects between Ca and La have been previously reported in plants 41,51 . Yuan et al. (2017) also demonstrated that increasing concentrations of Ca inhibited REE accumulation in P. americana 21 . Along with Ca, Mg also had a negative correlation with REE contents. Such deprivation had already been reported for P. americana exposed to REEs 21 but also in Brassica napus exposed to Ce 52 . Mg is necessary for chlorophyll synthesis and functioning; thus, the   decrease in chlorophyll content could be related to the Mg reduction due to REE exposure. A negative correlation between REE and Mn was also found in leaves and roots (Fig. 7). Phytolacca americana [17][18][19] and P. acinosa 25,26 are two known Mn hyperaccumulators. It is interesting to note the reduced Mn content that occurred concomitantly with the accumulation of REEs. This might suggest a competitive uptake between these elements in Phytolacca species as in P. acinosa between Mn and Cd 53 .
conclusions Phytolacca species are fast growing and high-biomass producing plants. P. americana is the first species of the genus that was demonstrated to be an REE accumulator. By testing five Phytolacca species under two different REE concentrations, this study therefore brings new information regarding the effect of REEs on the physiology and ionome of the REE-accumulating Phytolacca species. We used here hydroponic conditions to finely compare the potential of the different Phytolacca species to accumulate the whole set of REE species. Further research should be performed on the growth and REE accumulation potential of the different Phytolacca species with various REE-contaminated soils. These two parameters need to be evaluated in conditions of contamination to determine in situ the phytoextraction potential of the different Phytolacca species. It is interesting to note that the REE accumulating trait is conserved within the Phytolacca genus, as reported for the Carya genus 16 . Such information is very useful for the phytoremediation of REE-contaminated sites. Indeed, the different Phytolacca species tested here originate from different locations. For example, P. bogotensis is native to Colombia, P. icosandra to Mexico, and P. acinosa to Tibet 28 . Therefore, the conserved REE accumulation trait among these different species can be used to choose the most appropriate species according to the climatic and/or pedological context of the area to be remediated. Finally, to better understand the molecular and genetic mechanisms that underlie this particular trait, further work is needed.  54 . Two-litre black containers were filled with nutritive solution that was continuously aerated. Six-week-old seedlings were subsequently transferred to the same nutritive solution containing REEs or not (control medium) and lacking KH 2 PO 4 to avoid REE precipitation with phosphates 55 . Promethium was omitted since it is only available under its radioactive form. Two different total concentrations (10 and 100 µM) of a mixture of 16 REEs were used. All REEs were therefore applied at equimolar concentrations (625 nM or 6.25 µM). REE tri-chloride salts were used (Sc, Y, La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb and Lu) and were purchased from Sigma-Aldrich (Saint-Quentin-Fallavier, France). Three or four plants (control and 10 µM REE treatments, n = 3; 100 µM REE treatments, n = 4) per species and per treatment were grown individually in hydroponic containers in a growth chamber under controlled conditions: 18/23 °C, 12/12 h dark/light cycles (light intensity 260 µmol photons/m 2 /s), respectively.

Materials and Methods
After 21 days of exposure, plants were harvested. Roots were first rinsed in water, followed by a 15 min incubation in 5 mM ice-cold CaCl 2 , washed twice in deionized water and finally blotted dry carefully with tissue paper. Root architecture and leaf pigments were analysed (see below). Subsequently, plants were separated into roots, stems and leaves and weighed. Finally, the samples were ground into thin powder using liquid nitrogen and stored at −80 °C until use. Aliquots of powder samples were also dried at 70 °C for 2 days, weighed and further used for elemental analysis (see below).
The solutions were analyzed (see below) at the beginning and at the end of the experiment. The nominal and measured concentrations were very similar at the beginning of the experiment. Indeed, when compared to the nominal concentrations of 10 or 100 µM, the measured total REE concentrations varied only by a mean of 5.3 and 5.5%, respectively. At the end of the experiment, the REE concentrations were similar to the initial concentrations. As an example, after three weeks of incubation, the concentrations varied from 6.5% (100 µM treatment) to 11.4% (10 µM treatment) from the initial concentrations (data not shown). Even though plants took up a fraction of the REEs, this part did not contribute to deplete the medium and was counterweighed by a slight reduction of the volume of the growth medium due to evapotranspiration. As a matter of fact, the plants were exposed to similar concentrations all over the incubation period.
Species identification. Cellular extracts from the different Phytolacca species were obtained using the REDextract-N-Amp ™ Plant PCR kit (Sigma-Aldrich) according to the manufacturer's protocol. The internal transcribed spacer (ITS) region was amplified using the primers ITS1 and ITS4 56 . Twenty microliters of a mixture containing 0.2 µM of each primer, 4 µl of plant extract, 5.6 µl of H 2 O and 10 µl of REDExtract-N-Amp PCR ready mix were used for PCR amplification. The following PCR programme was used: 3 min at 94 °C, followed by 40 cycles of 94 °C for 1 min, 55 °C for 1 min and 72 °C for 1 min, and 10 min at 72 °C. PCR products were cloned into the pGEM-T ® Easy vector system I (Promega, Charbonnières-les-Bains, France), and the corresponding plasmids were sequenced. Related Phytolacca sequences were retrieved from the GenBank database using the NCBI BLAST program, and ClustalW was used to align sequences. The maximum likelihood tree was drawn with the kimura 2-parameter model using MEGA 7 57 . Bootstrap analysis was performed with 1,000 replications. (2019) 9:18458 | https://doi.org/10.1038/s41598-019-54238-3 www.nature.com/scientificreports www.nature.com/scientificreports/ Analysis of biomarkers. A series of biomarkers were used to assess the potential toxic effects of REEs on Phytolacca species. First, leaf pigments (chlorophylls, flavonoids and anthocyanins) and the nitrogen balance index (NBI) were quantified using a DUALEX ® SCIENTIFIC + apparatus 58 (Force-A, Orsay, France) before harvesting. Forty technical replicates were carried out per plant. Second, root architecture was analysed. Root systems were photographed at harvest. The analyses of lateral root length and lateral root density (number of lateral roots per cm of primary root) were carried out using the RootNav software 59 . Finally, the leaf and root contents of malondialdehyde (MDA) and total antioxidant compounds (TAC) were assayed from 50 mg of frozen powder using commercial kits (MAK085 and MAK187, Sigma-Aldrich, France) according to the manufacturer's protocols. Two technical replicates were performed per sample.
Elemental analysis. The quantification of elements was performed using 250 mg of dried plant powder. Mineralization was carried out using 3.5 mL nitric acid (HNO 3 , 65%, analytical grade, Optima) and 1 mL hydrogen peroxide, and samples were incubated in a heating block digestion system (DigiPREP, SCP Sciences, Courtaboeuf, France). A gradual heating mode was used to achieve a final temperature of 100 °C (total run of 265 min). Then, ultrapure water (18.2 MΩ/cm 2 , Millipore Milli-Q Integral 3 system, Molsheim, France) was added to a final volume of 25 mL, and samples were filtrated to 1 μm. Elemental concentrations were determined by ICP-AES (Radial ICAP 6500 Model, Thermo Fischer Scientific, Courtaboeuf, France) or ICP-MS (X Series II Model, Thermo Fischer Scientific, Courtaboeuf, France) for REEs. Oriental basma tobacco leaves (INCT-OBTL-5, LGC Promochem, Molsheim, France) were used as certified reference material.
Statistical analyses. One-way ANOVA and Tukey's HSD post hoc test were used when normality of data (Shapiro-Wilks test) and homoscedasticity (Levene's test) conditions were respected. When at least one condition was rejected, data were subsequently analysed using one-way non-parametric Kruskal-Wallis test followed by Wilcoxon post hoc test. The rejection level was set at α = 0.05 in all analyses. Percentage values were arcsin-transformed before statistical analyses. All analyses were performed using R software (version 3.4.1).