Photosynthesis Performance and Antioxidative Enzymes Response of Melia azedarach and Ligustrum lucidum Plants Under Pb–Zn Mine Tailing Conditions

Lead–zinc (Pb–Zn) mine tailings pose a great risk to the natural environment and human health because of their high toxicity. In this study, the responses of photosynthesis, chlorophyll fluorescence, and antioxidative enzyme of Melia azedarach and Ligustrum lucidum in the soil contaminated by Pb–Zn mine tailings were investigated. Results showed that Pb–Zn mine tailings significantly reduced net photosynthetic rates and leaf photosynthetic pigment content of both trees, and the reduction of net photosynthetic rates was mainly caused by their biochemical limitation (BL). The chlorophyll fluorescence parameters from Pb–Zn tailing stressed leaves indicated that Pb–Zn tailings affected PSII activity which was evident from the change values of energy fluxes per reaction center (RC): probability that an electron moves further than QA− (ETO/TRO), maximum quantum yield for primary photochemistry (TRO/ABS), the density of PSII RC per excited cross-section (RC/CSO), the absorption of antenna chlorophylls per PSII RC (ABS/RC), and the turnover number of QA reduction events (N). Pb–Zn mine tailings also affected the oxidation and reduction of PSI, which resulted in a great increase of reactive oxygen species (ROS) contents and then stimulated the rate of lipid peroxidation. Both trees exhibited certain antioxidative defense mechanisms as elevated superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT) activities, then declined under high level of Pb–Zn tailing treatment. Comparatively, L. lucidum showed less extent effect on photosynthesis and higher antioxidative enzyme activities than M. azedarach; thus L. lucidum was more tolerant than M. azedarach at least under the described Pb–Zn tailing treatment. These results indicate that the effect of Pb–Zn mine tailings on photosynthesis performance mainly related to imbalance of the PSII activity and PSI redox state in both trees. We propose that M. azedarach and L. lucidum could relieve the oxidative stress for phytoremediation under the appropriate Pb–Zn mine tailing content.


INTRODUCTION
Plants require certain heavy metals (HMs) as essential elements for their growth, development, and yield production, but excess amount of these metals can become phytotoxic and cause adverse effects on plant biomass production, crop yield, and food safety (Pierattini et al., 2017). The major sources of HMs originated from anthropogenic activities; mining activities seem to be the largest contributor of HM pollutions in many places, and it is a particular case in China (Ha et al., 2011;Han et al., 2013;Zhuang et al., 2014;Chen et al., 2015). It was recently reported that the mining activities have resulted in about 40,000 ha mining wastelands in this country, and the wastelands have been continuously expanded at a rate of 3,300 ha per year (Li, 2006;Luo et al., 2015;Yu et al., 2019). Mining activities generated a large amount of mine tailings in the mining sites where high concentrations of Pb, Zn, and other HMs were detected in local environments, which caused a wide range of environmental problems (Han et al., 2013). Therefore, it is urgent and necessary to reestablish vegetation in mining wasteland for ecological restoration (Yu et al., 2019).
Various methods have been identified and employed to the restoration of HM-polluted soils (Teng et al., 2015). Phytoremediation technology has been widely considered as an efficient, inexpensive, and environmental-friendly approach to clean up the contaminated environments by HMs (Han et al., 2016). Trees have the great potential use in remediation of HMpolluted soil in view of large biomass and massive root systems (Pulford and Watson, 2003). However, Tang et al. (2019) demonstrated that the Pb-Zn mine tailing significantly reduced plant growth and chlorophyll contents in three woody plants (Castanopsis fissa, Daphniphyllum calycinum and Pinus massoniana). The limitation in tree growth might thus be associated with photosynthesis which was sensitive to HM stress (Çicȩk et al., 2017;Zhong et al., 2018;Huang et al., 2019;Liang et al., 2019). Evidence showed that the reduction of leaf photosynthesis under HM toxicity might be attributed to the limitation of stomatal opening (Gs) and CO 2 diffusion in mesophyll (Gm), the suppression of photochemistry and biochemistry, or synthetic combinations of these factors (Kosobrukhov et al., 2004;Kola and Wilkinson, 2005;Krantev et al., 2008;Deng et al., 2014). Sagardoy et al. (2010) reported the main limitation to photosynthesis rate in sugar beet was likely due to the reduction of Gs under Zn stress. Lin and Jin (2018) found similar results in an experiment with three vegetables (B. chinensis, C. coronarium, and B. alboglabra) under Cu stress. However, Velikova et al. (2011) studied the impact of Ni on the photosynthesis of Populus nigra and found the limitation to photosynthesis rates of this species resulted from the restriction in Gm, not in Gs. Gm influenced photosynthetic capacity and determined the available CO 2 concentration for photosynthesis in the chloroplasts (Centritto et al., 2009). Some studies showed that the electron transfer was inhibited by Pb pollutant at the photosystem I (PSI) donor side because Pb affected activity of PSI (Belatik et al., 2013;Bernardini et al., 2016), but some studies suggested the activity of PSI in Microcystis aeruginosa and Chlorella pyrenoidosa had no inhibition with Cd treatment (Zhou et al., 2006;Wang et al., 2013). And some studies indicated that Cd affected the whole electron transport in photosystem II (PSII): on the donor side, it inhibited the OEC, while on the acceptor side, it inhibited electron transport between Q A − and Q B − . (Mallick and Mohn, 2003;Sigfridsson et al., 2004;Chu et al., 2018). Additionally, a number of studies also reported that HMs exhibited less effect of photosynthetic rate, electron transport, conversion of light energy, and photochemical efficiency in tolerant plant species than those in sensitive ones (Guo et al., 2018;Sorrentino et al., 2018;Huang et al., 2019).
Reactive oxygen species (ROS) can highly increase in the chloroplasts when photosynthetic electron transport chain was blocked by HM toxicity . ROS like superoxide (O 2 − ) and hydrogen peroxide (H 2 O 2 ) could lead to lipid peroxidation, protein oxidation, membrane and nucleic acid damage, and inactivation of enzymes (Bi et al., 2016;Bezerril et al., 2017;Lu et al., 2017). To prevent an oxidative damage, plants activated enzymatic ROS scavenging mechanisms, such as superoxide dismutase (SOD), peroxidase (POD) and catalase (CAT), ascorbic acid and glutathione, to keep ROS at a basal non-toxic level (Israr et al., 2011;Sťefanićet al., 2018;Du et al., 2020). In the past two decades, many studies have estimated the direct link between the oxidative stress in metal toxicity and metal-tolerant plants (Boominathan and Doran, 2003;Chiang et al., 2006). A large proportion of studies have indicated that metal tolerant plants were linked to superior constitutive antioxidative defense (Srivastava and Doran, 2005;Chiang et al., 2006;Nie et al., 2016). Although great progress has been made in supporting the HM toxicity to plant photosynthesis and its redox balance, the effects of HM stress on photosynthetic performance varied depending on the plant species, metal ion and concentration. Meanwhile, these available data provide a restricted view on single metal contamination or herbaceous plants (Mobin and Khan, 2007;Wali et al., 2016). Therefore, the photosynthesis and redox responses of trees to HM tailing stress still need to be studied further.
Melia azedarach is a fast-growing, deciduous broad-leaved tree species and is widely distributed in the southern regions of China. It is mainly planted for reforestation as a useful timber production species or as an ornamental plant. Ligustrum lucidum is a commonly-seen evergreen broad-leaved tree species in South China. This species is often planted as an ornamental tree in the urban. Moreover, both M. azedarach and L. lucidum are native plants in the southern regions of China, Freŕot et al. (2006) pointed out that native plants could be a useful option to phytoremediation because native plants are better adapted to local climate conditions than plants introduced from other environments, and they were previously found to have high tolerance in Cd or Mn contaminated soils (Trikshiqi and Rexha, 2015;Su et al., 2017;Liang et al., 2019). However, the metal toxicity manifestations and mechanisms behind tolerance need further investigation, so the specific goals of this study were: 1) to investigate the contributions of Gs, Gm, and biochemical component respectively to the photosynthesis reduction in two tested plants by photosynthesis limitation analysis, 2) to identify the impact of Pb-Zn stress on the whole photosynthetic electron flow chain from PSII to PSI and the state of PSI, 3) to explore the tolerant mechanisms of both trees to Pb-Zn tailing stress in terms of the redox responses of photosystems and antioxidative enzymes. Based on the above studies, the differences of tolerance between M. azedarach and L. lucidum were discussed. Results could provide a theoretical basis for selecting and breeding the resistant trees species grown in Pb-Zn polluted environments.

The Physicochemical Properties of Pb-Zn Tailings and Experimental Soil
The study was conducted in the Central South University of Forestry and Technology (CSUFT), Changsha City, Hunan Province, China (28°8′12″N, 112°59′36″E). The Pb-Zn mine tailing samples were collected from a Pb-Zn mining site in Suxian district, Chenzhou City, Hunan Province, China (25°30′ 38″-25°00′19″N, 112°16′41″-112°53′23″E). The soil samples were collected from the top soil (5-20 cm) in a garden field of CSUFT campus. Both Pb-Zn mine tailing samples and top soil samples were air-dried at room temperature. The large debris, stones, and pebbles were manually removed before being applied to the pot experiment. The soil pH value, soil total carbon (TC), total nitrogen (TN), total phosphorus (TP), and soil heavy metal content were measured in laboratory according to Bao (2000), the basic physicochemical properties of soil were as follows: pH 4.95, total C 5.82 g kg −1 , total N 0.33 g kg −1 , total P 0.16 g kg −1 , Pb 0.002 g kg −1 , Zn 0.003 g kg −1 . The properties of Pb-Zn mine tailing samples were determined: pH 3.89, total C 13.89 g kg −1 , total N 1.25 g kg −1 , P 0.82 g kg −1 , Pb 8.92 g kg −1 , Zn 14.41 g kg −1 .

Plant Seedlings
Two-year-old young plants of M. azedarach (mean tree height: 73.5 cm, mean stem base diameter:~0.6 cm) and L. lucidum (mean tree height:~130.0 cm, mean stem base diameter:~0.9 cm) were purchased from a local nursery.

Experimental Design
Pb-Zn mine tailings represent a poor spoil substrate characterized by high contents of Pb and Zn, low levels of organic nutrients, and poor physical structure and water retention capacity (Meeinkuirt et al., 2012); few plants can survive in such harsh environmental condition. Mixing mine tailings with garden soil helps decrease the bioavailability of metals and improves mine tailing structure and ultimately upgrades the physical properties and nutrient status of mine tailings. Moreover, this method has been commonly used in many studies of phytoremediation in mine tailing Han et al., 2016;Tang et al., 2019). In the present study, four treatments were set up in the pot experiment with different weighted proportions of Pb-Zn mine tailings and garden soils. Each pot contained 10 kg mixed soils. The four treatments were: (1) 90% garden soils + 10% Pb-Zn mine tailings (designed as L1, the 10 kg mixed soils were 9 kg garden soil and 1 kg mine tailings, the same as below); (2) 75% garden soils + 25% Pb-Zn mine tailings (L2); (3) 50% garden soils + 50% Pb-Zn mine tailings (L3), and (4) 100% garden soils + 0% Pb-Zn mine tailings as the control (C). The Pb-Zn mine tailings and garden soils were completely mixed and then the M. azedarach and L. lucidum young plants were transplanted into the pots with one plant per pot. Each treatment was replicated six times. The temperatures of the greenhouse were set at 30/25°C (day temperature for 10 h and night temperature for 14 h), and relative humidity was 65/ 85%. During the study period, each pot was supplied with equal quantity of pure water every 1 to 2 days until the young plants were harvested. Both trees were measured for all physiological parameters after growing for 425 days (from June 15 th , 2017 to September 13 th , 2018) in Pb-Zn tailing treatments.

Relative Growth Rate (RGR)
Plant biomass was measured using a harvesting method. The plant was divided into leaves, stem, and root components; the fresh weights of each component were measured by using an electronic balance and then dried at 70°C until constant weight was reached. The dry weight (DW) of each component was determined by using an electronic balance. Six individual plants were selected for biomass measurement for each of M. azedarach and L. lucidum species at the beginning of the treatment, respectively, as initial dry weight (DW). The biomass of each component was determined as initial dry weight (DW). After 425 days of treatment (from June 15 th , 2017 to September 13 th , 2018), all examined plants were harvested and the biomass was measured as final DW. The determination of relative growth rate (RGR) is based on the method of Environment Canada (2007) guidelines as follows: where X i and X j represent the values of final DW and initial DW, respectively; t j and t i represent initial time and final time, respectively.

Gas Exchange Measurements
Three new, similar size and fully expanded leaves per plant were chosen for the measurements of photosynthesis for each treatment. The leaf gas exchange and Chl fluorescence were measured by using an open gas exchange system (LI-COR 6400XT, Lincoln, USA) at the same time with an integrated Chl fluorescence chamber head in the morning (8:00-10:00 am). For all measurements, the following conditions were set up: leaf temperature of 25-32°C, PAR of 1,000 mmol (photon) m −2 s −1 , and vapor pressure deficit (VPD) of 2.0 ± 0.2 kPa.
Mesophyll conductance to CO 2 (Gm) was calculated by using 'variable J method' (Harley et al., 1992;Sagardoy et al., 2010) as: where, P n and Ci were obtained from gas exchange measurements under saturating light, G * was the CO 2 concentration at the compensation point in the absence of mitochondrial respiration and was obtained according to Bernacchi et al. (2002), R d was calculated based on the A-Ci curve on the same leaf according to Pinelli and Loreto (2003).
J flu represented the rate of electron transport and was calculated as: where, a was total leaf absorbance in the visible light range (taken as 0.85, Grassi and Magnani, 2005), 0.5 was a factor to account for the distribution of light between the two photosystems (Laisk and Loreto, 1996). The actual chloroplastic CO 2 concentration (Cc) was calculated from the g m value as Cc = Ci − (P n /Gm) (Harley et al., 1992).

Photosynthesis Limitation Analysis
The inhibition of Pb-Zn stress on plant photosynthesis was further assessed based on the contribution made by various functional components to the photosynthetic limitations (Sagardoy et al., 2010). The functional components were stomatal (SL), mesophyll conductance (MCL), and leaf biochemical characteristics (BL). The relative contributions of these functional components to the photosynthesis limitations were evaluated when compared with a reference status in which the photosynthesis limitations were ignorable, and the Gs, Gm, and V cmax were at their maximum. In this study, the corresponding values of Gs, Gm, and V cmax taken from C treatments were as reference values; thus the photosynthesis limitations were set to 0.

The Kinetics of Prompt Fluorescence and Modulated 820 nm Reflection
Fast chlorophyll a fluorescence was measured by using M-PEA (Multifunctional Plant Efficiency Analyser, Hansatech Instrument, UK). After 1 h dark adaptation using dark adaptation clips, leaves were exposed to a pulse of saturating red light (5,000 mmol m −2 s −1 , peak at 625 nm, duration from 50 ms to 2 s, records of 128 points), and the measurements were carried out during a period of 8:30-11:00 am. The OJIP transient was analyzed based on the JIP test (Strasser et al., 2004). The values of fluorescence intensity from the original measurements were used in this study: fluorescence intensity at 20 ms (at the O step, considering as the minimum fluorescence, F O ); 300 ms (F 300ms ) used for calculation of the initial slope of the relative variable fluorescence kinetics; 2 ms (at the J step, F J ); 30 ms (at the I step, F I ), and the P step (considering as the maximum fluorescence, F m ). The description and calculation of standardization formula of OJIP transients and formula of parameters were listed in Table 1.
Modulated 820 nm reflection was also expressed as MR/MR O . The value of MR and MR O were determined by the method of Zhou et al. (2019). The PSI redox was denoted by V PSI (maximum PSI oxidation rate) and V PSII-PSI (maximum PSI reduction rate, respectively, and were obtained according to Gao et al. (2014).

Evaluation of Chlorophyll and Carotenoid Contents
Leaf chlorophylls and carotenoid content were extracted using 95% ethanol and then placed in the dark for at least 72 h at 4°C. The extracts were measured at wavelengths of 665, 649, and 470 nm spectrophotometrically and was then calculated as milligram per gram of fresh weight (Ji et al., 2018).
The O 2 − production rate was measured using a reagent kit (Beijing Solarbio Science and Technology, China). H 2 O 2 content was measured using the method of Okuda et al. (1991). The MDA level was assayed as a thiobarbituric acid reactive substance according to the method of Dhindsa et al. (1981). A 0.25 g fresh leaf sample was homogenized; supernatant was collected and measured at the wave length of 532, 600, and 450 nm using a UV/visible spectrophotometer.

Estimation of Proline and Antioxidative Enzymes
Proline content was determined using Pure Pro as a standard (Gao, 2006). A 0.1 g fresh leaf sample was homogenized with 5 ml of 3% aqueous sulfosalicylic acid and was then extracted in a boiling bath.
For determination of SOD, POD, and CAT, a 0.3 g fresh leaf sample was grounded in 3 ml extraction buffer containing 25 mM Hepes, 2% polyvinyl-pyrrolidone (PVP), 0.2 mM EDTA, and 2 mM ascorbate (pH 7.8) and was centrifuged at 12,000 g for 20 min at 4°C. The supernatants were collected for enzyme analysis.
SOD activity was determined according to Giannopolitis and Ries (1977). The 3 ml reaction mixture contained 0.05 M phosphate buffer solution (pH 7.8), 0.06 M Riboflavin, 0.195 M Met, 0.003 M EDTA, 1.125 mM NBT, and 0.2 ml enzyme supernatants. The measurements were performed at 25°C. The tested samples were incubated for 10-20 min under 10,000 lx irradiance; inhibition rate of nitro blue tetrazolium (NBT) reached to 50% represented one unit of SOD activity as by spectrophotometer at 560 nm.
POD activity was assessed as in Beffa et al. (1990). The reaction mixture contained 0.1 M sodium-acetic buffer (pH 5.0), 0.25% (w/v) guaiacol, 0.75% H 2 O 2 , and 0.05 ml enzyme supernatants. One unit of POD activity was represented as an increase of 0.01 DOD value a minute at 470 nm.
CAT activity was determined according to Aebi (1984). The reaction mixture contained 0.05 M phosphate buffer (pH 7.0) and 0.45 M H 2 O 2 and 0.2 ml enzyme supernatants. One unit of CAT activity was represented as decrease of DOD value a minute at 240 nm. The measurement for each antioxidative enzyme was repeated three times.

Assay of Leaf Total Rubisco Activity
Ribulose-1, 5-bisphosphate carboxylase/oxygenase (Rubisco) of leaves from the tested plants was assayed according to the method of Chen et al. (2005). After centrifugation at 13,000 r·min −1 40 s at 2°C, the supernatant was used immediately for assays of Rubisco activity (Jiang et al., 2008).

Total Metal Content
For the determination of the metal contents in plant, the different dried plant components were powdered and passed through an 80-mesh sieve. About 0.1 g of plant material was digested with 1 ml of HNO 3 and H 2 O 2 (8:2, v/v). Then, the tube with plant material was on the block at 200°C for 0.75-1 h. The residue was taken up in 10 ml of demineralized water. The concentrations of Pb and Zn were determined by atomic absorption spectrophotometry (AA-6800, Shimadzu, Kyoto, Japan) . The measurement for each sample was repeated six times, and the mean value was calculated.
The metal translocation factor (TF) in the leaves was represented as the metal concentration ratio of plant leaves and stems to roots (Han et al., 2016).

Statistical Analyses
Two-way analysis of variance (ANOVA) was performed on the data using SPSS (20.0). Differences among the eight treatment combinations (two species × four Pb-Zn mine tailing treatments) were analyzed by two-way analysis of variance; eight means were separated by Duncan's new multiple range test at P < 0.05 level. Data were presented as means ± SD (n = 6). The principal component analysis (PCA) used CANOCO version 5.0.

RGR, Pb/Zn Contents in Roots, Stems, and Leaves
As Pb-Zn tailing portions increased, RGR of M. azedarach and L. lucidum was decreased progressively compared with the control group (P < 0.05) ( Table 2), the RGR values were reduced by 10-90% and 6-70% in M. azedarach and L. lucidum, respectively, in Pb-Zn treatments when compared with the C plants ( Table 2). The concentrations of Pb and Zn in the leaves, stems, and roots increased with the increase of the proportion of Pb-Zn tailings in both tested plants compared to the C ( Table 2). The Pb and Zn concentrations were significantly higher in the roots than in the leaves and stems. The Pb and Zn concentrations showed as Zn > Pb in the leaves, stems, and roots in all Pb-Zn tailing treatments. TFs of Pb and Zn were higher in M. azedarach than those in L. lucidum. Maximal recorded fluorescence intensity, at the peak P of OJIP N the turnover number of QA reduction events Probability that an electron moves Quantum yield for reduction of the end electron acceptors on the PSI acceptor side (RE) Probability that an electron is transported from the reduced intersystem electron acceptors to the final electron acceptors of PSI (RE)

Pigment Contents
The photosynthetic pigment contents of the two tested plants were found to decline with increasing Pb-Zn tailing portions ( Figures 1A-D). A significant reduction in chlorophylls a, b, and total chlorophyll in the two tested plants was observed under Pb-Zn tailing stress as compared to the C, respectively. Carotenoid content was notably affected by Pb-Zn tailing stress when the Pb-Zn treatment exceeded a portion of 10%. Compared with L. lucidum, M. azedarach showed a faster increase and a greater extent of chlorophylls a, b, total chlorophyll, and carotenoid content. When compared to the C groups, Chl a decreased by a range of 18 to 50%, Chl b from 9 to 35%, total chlorophyll from 16 to 46%, and the carotenoids from 5 to 37% in L. lucidum in L1-L3 treatments. The corresponding values decreased in M. azedarach from 4to 42% in Chl a, from 5 to 21% in Chl b, from 5 to 37% in total chlorophyll, and from 3 to 34% in carotenoids.

Gas Exchange
The Pb-Zn tailing treatments had an inhibitory effect on Pn, Gs, Gm, V cmax , J max , and Rubisco activity, but had promotion effect on Ci and Cc in both examined plants ( Figure 2). Pn, V cmax , J max significantly decreased after Pb-Zn mine tailing treatment. The Gm of the leaves from the two plants gradually decreased with the increase of Pb-Zn tailing portions. For M. azedarach, Gs showed a significant difference between C and L1 treatments and then gradually decreased when the plants were grown in L2 treatments, and eventually reached the minimum value at L3 treatments. In L. lucidum, Gs was decreased slightly as Pb-Zn tailing portions increased. Both plants exhibited the similar tendency in Ci and Cc as Pb-Zn tailing portions increased. Ci and Cc decreased slightly when the plants were in L1 treatments and then significantly increased in L2 and L3 treatments (P < 0.05), but the change of Ci was higher in L. lucidum (increased from −5 to 58%) than in M. azedarach (increased from −4 to 5%). The decrease of P n , Gm, V cmax , and J max in M. azedarach showed a greater extent than that in L. lucidum. When compared to the C groups, the Pn decreased by a range of 12 to 56%, the Gm from 7 to 18%, the V cmax from 9 to 43%, the J max from 19 to 64%, and the Rubisco activity from 10 to 52% in L. lucidum in L1-L3 treatments.
The corresponding values decreased in M. azedarach from 31 to 67% in Pn, from 10 to 28% in Gm, from 16 to 59% in V cmax , from 21 to 72% in J max , and from 31 to 62% in the Rubisco activity. Under Pb-Zn tailing treatments, biochemical limitations (BLs) increased in the two tested plants. The BL values of M. azedarach increased 12, 34, and 53% in the L1, L2, and L3 treatments compared to the C, respectively ( Table 3). The corresponding values increased 6, 23, and 36% in L1, L2, and L3 treatments in L. lucidum, respectively. The stomatal conductance limitation (SL) and the mesophyll conductance (MCL) of leaves in both trees increased slightly under Pb-Zn treatments. The SL values increased 4-6% and the MCL values increased 0.4-0.8% for M. azedarach and L. lucidum when compared to C.

Prompt Fluorescence OJIP Transient Analysis
The Pb-Zn stress had a considerable effect on fluorescent OJIP transients in two tested plants (Figure 3). The F o values were 0.53 ± 0.11c TF = the metal concentration ratio of plant leaves and stems to roots. Data are means ± SE, n = 6. Values with different letters within the same column indicate significant differences at the p < 0.05 level between Pb-Zn treatments. C: 100% garden soil, L1: 90% garden soil + 10% Pb-Zn tailings, L2: 75% garden soil + 25% Pb-Zn tailings, L3: 50% garden soil + 50% Pb-Zn tailings.
significantly reduced under Pb-Zn treatments in M. azedarach compared to the C ( Figure 3A, inset). However, the F o values of L. lucidum had no significant increase under Pb-Zn tailing treatments compared to the C ( Figure 3B, inset). Three Pb-Zn tailing treatments did not affect W k in L. lucidum compared to the C ( Figure 4A). In M. azedarach, W k was similar for both C and L1 treatments and significantly increased in L2 and L3 treatments ( Figure 4A). Pb-Zn stress treatments significantly decreased performance index (PIabs) in M. azedarach and L. lucidum ( Figure 4B) but distinctly increased the turnover number of Q A reduction events (N) and absorption flux (ABS/ RC) ( Figures 4F, G)

MR/MR O Transient Analysis
After Pb-Zn tailing treatments, the shape of the MR/MR O kinetics was obviously changed in the two tested plants. The lowest points represented a turning point of PSI oxidation state ( Figure 5). When compared to the C, the lowest points on the reflection curve of both tested plants increased with rising Pb-Zn tailing portions. The V PSI and V PSII-PSI of both tested plants declined significantly with the rising Pb-Zn tailing portions ( Table 4). V PSI and V PSII-PSI of both tested plants significantly decreased in the treated groups than in the C groups. Compared with L. lucidum, M. azedarach showed a greater decrease of V PSI and V PSII-PSI .
PI abs , DV IP in Relation to V cmax A significant positive relationship of PI abs , DV IP and V cmax was observed in both tested plants under Pb-Zn mine tailing treatments ( Figure 6). As V cmax decreased, PI abs , DV IP decreased linearly in M. azedarach and L. lucidum, respectively.

O 2 − Production Rate, H 2 O 2 , MDA and Proline Content in Leaves
With an increase of Pb-Zn mine tailing portions, the O 2 − production rate, H 2 O 2 , and MDA contents of two tested plant leaves increased notably ( Figures 7A, C, D). Compared with L. lucidum, M. azedarach showed a faster increase and a greater extent O 2 − production rate, H 2 O 2 , and MDA contents. The proline content in M. azedarach and L. lucidum reached the peak at L2 treatments ( Figure 7B), where the proline content was 1.96 and 1.36 times higher than that in the C treatments, proportions of Pb-Zn mine tailings. C (control), L1, L2, and L3 represent 100% garden soil, 90% garden soil + 10% Pb-Zn tailings, 75% garden soil + 25% Pb-Zn tailings, 50% garden soil + 50% Pb-Zn tailings, respectively. Different letters above the bars indicate a significant difference at P < 0.05. Values are means of n = 6; bar indicates standard error. (E) V cmax , maximum carboxylation rate allowed by Rubisco; (F) J max , rate of photosynthetic electron transport based on NADPH requirement; (G) Cc, chloroplastic CO2 concentration; (H) Rubisco activity. C (control), L1, L2, and L3 represent 100% garden soil, 90% garden soil + 10% Pb-Zn tailings, 75% garden soil + 25% Pb-Zn tailings, 50% garden soil + 50% Pb-Zn tailings, respectively. Different letters above the bars indicate a significant difference at P < 0.05. Values are means of n = 6; bar indicates standard error. respectively. Under L3 treatments, the proline content in M. azedarach and L. lucidum decreased by 16 and 13% of that in the C, respectively.

Antioxidant Enzyme Activity in Leaves
The SOD, POD, and CAT displayed similar trends in activity in both tested trees with the increase of Pb-Zn mine tailing portions ( Figures 8A-C). Compared with the C, the activities of SOD and POD in both tested trees significantly increased under L1 and L2 treatments and decreased under L3 treatment. The CAT activity in both tested trees under all Pb-Zn treatments was significantly higher than that in the C plants. Notably, the activities of SOD, POD, and CAT in L. lucidum were significantly higher than that in M. azedarach under all Pb-Zn mine tailing treatments.

Principal Component Analysis
PCA was used to understand Pb-Zn tailings affecting the photosynthetic characteristics including RGR, PSII performance, PSI content, net photosynthesis rates (Pn), and the antioxidative enzymes of leaves in both trees ( Figure 9). The first two components comprised 90.0% (50.3% for PC1 and 39.7% for PC2), 91.4% (62.9% for PC1 and 28.5% for PC2), and 90.9% (36.7% for PC1 and 54.2% for PC2) of the total variations in L1, L2, and L3 treatments, respectively. Under L1 and L2 treatments, PIabs, DV IP , and Pn were the most influential in the PC1; and SOD, POD and CAT in the PC2. Under L3 treatments, SOD, POD, and CAT were the most influential in the PC1, and PIabs and DV IP in the PC2. Evidently, two tested plans were more closely to photosynthetic parameters (PIabs, DV IP , and Pn) than to antioxidative enzymes (SOD, POD and CAT) under L1 and L2 treatments, but were more closely to antioxidative enzymes than to photosynthetic parameters under L3 treatment. In addition, RGR was positively related to the photosynthetic parameters in the tested plants.

The Reduction of Plant Growth Is Attributed to the Depression of Photosynthesis
When plants were grown in Pb-Zn contaminated soils, plant growth and development were retarded and eventually the biomass production decreased (Ha et al., 2011;Han et al., 2013). In this study, the RGR values for the two tested trees decreased under Pb-Zn treatments ( Table 2). It was well known that plant biomass productivity was dependent on the photosynthetic assimilating accumulation which provided energy and carbon sources . A positive correlation was found between Pn and RGR in this study (Figure 9), suggesting the photosynthetic processes were inhibited by Pb-Zn stress in the tested tree species.

The Photosynthesis Limitation Is Mainly Attributed to the Biochemical Limitation
A significant decrease of gas exchange parameters was observed in the Pb-Zn treated plants ( Figure 2). Photosynthesis limitations might be derived from different physiological  represent 100% garden soil, 90% garden soil + 10% Pb-Zn tailings, 75% garden soil + 25% Pb-Zn tailings, 50% garden soil + 50% Pb-Zn tailings, respectively. In the insets of the two panels, the F O of the fluorescence rise is compared between measurements on control and Pb-Zn treated leaves.
processes, such as stomatal opening, mesophyll conductance to CO 2 , and the carboxylation capacity (Flexas et al., 2008;Centritto et al., 2009). In this study, although Gs and Gm significantly decreased ( Figures 2B, D), they probably played a minor role in limiting photosynthesis because Ci and Cc increased simultaneously (Table 3, Figures 2C, G). The biochemical limitation was likely the major factor affecting photosynthesis in M. azedarach and L. lucidum (Table 3). Change for the results were in line with the findings of other studies (Seregin and Kozhevnikova, 2006;Velikova et al., 2011). The biochemical limitations were related to the considerable declines in V cmax and J max with increasing Pb-Zn stresses. The M. azedarach suffered more metabolic damage than that in L. lucidum ( Figures 2E, F), indicating that the supply of energy source and carbon skeleton (H) j Ro , quantum yield for reduction of the end electron acceptors on the PSI acceptor side (RE); (I) d Ro , probability that an electron is transported from the reduced intersystem electron acceptors to the final electron acceptors of PSI (RE); (J) DV IP , amplitude of the I to P phase of the OJIP fluorescence transient (associated with PSI reaction center content). C (control); L1, L2, and L3 represent 100% garden soil, 90% garden soil + 10% Pb-Zn tailings, 75% garden soil + 25% Pb-Zn tailings, 50% garden soil + 50% Pb-Zn tailings, respectively. Different letters above the bars indicate a significant difference at P < 0.05. Values are means of n = 6, bar indicates standard error.
for plant growth and development was lower in M. azedarach than in L. lucidum (Imsande and Touraine, 1994). The decrease of V cmax values might be ascribed to a reduction of Rubisco activity in the photosynthetic process ( Figure 2H). The influence of Rubisco activity due to adverse environmental conditions has been reported by other studies (Bah et al., 2010). The J max inhibition in the leaves under Pb-Zn tailing treatments observed from the changes in A/Cc curves indicated the photochemical limitation occurred in our study. The results were in line with the finding which Velikova et al. (2011) reported that alterations in the electron transport rate from PSII to PSI under heavy metals stresses result in a J max limitation.
In addition, there was a decrease in chlorophyll and carotenoid content with the rising Pb-Zn stress in M. azedarach and L. lucidum, which illustrated that Pb-Zn stress had harmed the photosynthetic apparatus for these two tree species. The decrease of chlorophyll content might be attributed to HMs interfering Fe metabolism, inhibiting chlorophyll synthetase activity and enhancing chlorophyll enzyme activity (Rana, 2015). As a result, chlorophyll degradation ultimately inactivated photosynthesis (Hajihashemi and Ehsanpour, 2013). We found that the decrease degree of chlorophyll content was higher than that of carotenoid content by HMs. The results were in line with the previous finding by Amir et al. (2020) who reported that chlorophylls were more susceptible to the negative impact of HMs as compared to carotenoids.  , and (C) L3 represent 100% garden soil, 90% garden soil + 10% Pb-Zn tailings, 75% garden soil + 25% Pb-Zn tailings, 50% garden soil + 50% Pb-Zn tailings, respectively.

The Photosynthetic Electron Transport Is Altered by the Pb-Zn Tailing Stress
Changes in environmental factors could affect photosynthetic performance (Ji et al., 2018). The impacts of HM stresses on the electron transfer and energy balance in photosynthetic processes are quantified by OJIP fluorescence (Li and Zhang, 2015;Guo et al., 2018). In this study, the chlorophyll fluorescence transient curves of the two tested plants were all modified by Pb-Zn mine tailings (Figure 3). With the increase of proportions of Pb-Zn mine tailings, F O values were significantly reduced in M. azedarach, but increased in L. lucidum compared to the control ( Figure 3); the value of F P decreased in the two tested plants. The F O decrease in M. azedarach was likely attributed to the photoinhibitory damage on PSII acceptor (Setlik et al., 1990), and the increase of F O in L. lucidum might be caused by increasing amount of free chlorophylls to the PSII RC (Gilmore et al., 1996). The decrease in F p was observed in other heat and salt stress studies due to the increased fraction of inactive RCs (Jedmowski and Brüggemann, 2015;Oukarroum et al., 2015). W K was the donor site parameter of PSII, which has been widely used to analyze damage to OEC by HM stresses (Li and Zhang, 2015). In this study, an increase of W k value in M. azedarach ( Figure 4A) suggested that Pb-Zn stress not only impaired the OEC; the process of electron transfer was also impaired at the P680 donor site (Jiang et al., 2008;Dąbrowski et al., 2016). However, the OEC and the donor side of PSII were , and L3 represent 100% garden soil, 90% garden soil + 10% Pb-Zn tailings, 75% garden soil + 25% Pb-Zn tailings, 50% garden soil + 50% Pb-Zn tailings, respectively. Different letters above the bars indicate a significant difference at P < 0.05. Values are means of n = 6; bar indicates standard error.
not impaired by Pb-Zn stress in L. lucidum since W k values did not show a significant change.
Here, several important JIP-test parameters based on the OJIP transient were used to detect and quantify the changes of the photosystem status under Pb-Zn mine tailings stresses. The reduction of the maximum quantum yield of primary photochemistry (TRo/ABS) ( Figure 4D) and the electron numbers moved further than the Q A − (ET O /TR O ) ( Figure 4E) indicating that Pb-Zn stress resulted in a fast accumulation of Q A − in the PSII reaction centers. The results suggested the photosynthetic electron transport was impaired in the acceptor side of PSII, and the electron flow beyond Q A − was blocked (Bernardini et al., 2016). It was reported that photoinhibition process can be more accurately assessed by the changes between flux ratios and specific fluxes per reaction center (RC) (Baker and Rosenqvist, 2004;Mlinarićet al., 2017). It was worth noting that the RC/CS O was diminished and ABS/RC was risen in the two tested plants ( Figures 4C, F), suggesting that although the numbers of active PSII RCs were decreased under Pb-Zn stress, the photochemical reaction efficiency of the remaining RCs was improved by maintaining the regular photosynthesis ability (Paunov et al., 2018). Due to the inactivation of a number of reaction centers, the remaining RCs had to increase their turnover in order to completely reduce the PQ pool (Zhou et al., 2019) as shown by a rise in parameter N ( Figure 4G). M. azedarach had lower values of PIabs, RC/CS O , TR O /ABS, and ET O /TR O and higher values of ABS/RC than that in L. lucidum, indicating that Pb-Zn stress may affect, to a greater extent, the PSII activity in M. azedarach species.
DV IP has been used to describe the changes in content of PSI reaction center under detrimental conditions (Oukarroum et al., 2009;Ceppi et al., 2012). The drought and Zn stress conditions declined DV IP in Hordeum spontaneum (Oukarroum et al., 2009) and Phragmites australis (Bernardini et al., 2016). In this study, DV IP was significantly lower for the two plants grown in Pb-Zn contaminated soil when compared to C ( Figure 4J), indicating that Pb-Zn stress decreased the PSI reaction center in the two tested tree species. The significant reductions of jR o and dR o (Figures 4H, I) indicated that the effect of Pb-Zn stress on electron flow occurred in the acceptor side of PSI because of the reduction of the PSI content (Hu et al., 2018). In addition, the values of V PSI of the two tested plants decreased in the Pb-Zn treated groups when compared to the C ( Figure 5, Table 4), indicating that the oxidization reactions were inhibited due to the Pb-Zn stress. The decline in the V PSII-PSI suggested that the reduction activity of PSI activity was injured by the Pb-Zn stress. Therefore, the decreases in maximum oxidization and reduction activity of PSI reaction center were observed since Pb-Zn stress affected PSI content.
In this study, a positive linear relation was found between PI abs and V cmax in this study suggesting that PSII performance was positively related to CO 2 assimilation ( Figure 6). The result was in accordance with the suggestions provided by Ghotbi-Ravandi et al. (2014); they reported that the reduction in PSII performance played an important role in the decrease of CO 2 assimilation rate in Morocco under mild and severe drought stresses. Meanwhile, DV IP also decreased linearly with the decrease of V cmax ; the result suggested that PSI content might limit the capacity of CO 2 assimilation ( Figure 6). Because the decrease of PSI content could disturb the electron flow from PSII to PSI and limited the synthesis of ATP and NADPH (Brestic et al., 2015), the inhibition of CO 2 assimilation rate can be suggested as a result of the excessive excitation energy that damage the photosystems and especially impeded the photochemical activity of PSI (Zivcak et al., 2013). Specifically, M. azedarach suffered more serious effects on the photosynthetic electron transport chain than L. lucidum under Pb-Zn stress. This finding was supported by the fact that photosynthetic parameters (PIabs, DV IP , and Pn) were more affected in M. azedarach than that in L. lucidum in PCA analysis (Figure 9).

L. lucidum May Be More Tolerant to Pb-Zn Stress Than M. azedarach
It is previously demonstrated that adverse conditions impaired PSII electron transport; the excess electron resulted in increasing the levels of ROS and consequently caused oxidative stress (Pospıśǐl, 2009;Guo et al., 2018). The malondialdehyde (MDA) content represents level of lipid peroxidation, which was able to indicate A B C FIGURE 9 | Principal component analysis (PCA) of antioxidant enzymes, PSII performance, PSI content, net photosynthesis rates (Pn), and the relative growth rate (RGR) of M. azedarach and L. lucidum exposed to different proportions of Pb-Zn mine tailings. (A) L1, (B) L2, and (C) L3 represent 100% garden soil, 90% garden soil + 10% Pb-Zn tailings, 75% garden soil + 25% Pb-Zn tailings, 50% garden soil + 50% Pb-Zn tailings, respectively.
oxidative damage of membrane lipids under stress conditions (Sharma and Dubey, 2005). In this study, with the increase of proportion of Pb-Zn mine tailings, the O 2 − production rate and H 2 O 2 content in both tested tree leaves increased significantly, contributing to a significant increase of MDA content. This result was in line with the previous reports where MDA content increased in Peganum harmala, Kandelia obovate, and Alternanthera bettzickiana under HMs stress (Lu et al., 2010;Cheng et al., 2017). In this study, proline contents increased in M. azedarach and L. lucidum leaves in L1 and L2 treatments but decreased at L3 treatments when compared to the C ( Figure 7B). The phenomena suggested the increase of proline contents at low and medium Pb-Zn treatments could maintain the integrity of cellular membranes, protect proton pump, and eliminate ROS (Dhir et al., 2012). The decrease of proline contents in L3 treatments might be explained by the fact that the physiological functions and metabolisms of plants were seriously damaged at high concentration of Pb-Zn toxicity (Chen et al., 2003). When the equilibrium between ROS generation and detoxification was disrupted by abiotic stresses, the induction of antioxidative enzyme defense activities played a crucial role in HM tolerance in plants (Ashraf, 2009). SOD was the first defense line against oxidative stresses, and superoxide radicals (O 2 − ) were scavenged by SOD. CAT was recognized as the most important enzyme for scavenging H 2 O 2 produced in plant cells and was primarily associated with the maintenance steady of cellular. POD enzyme can detoxify H 2 O 2 and thereby was conducive to maintaining the integrity of cellular membranes . In this study, SOD, POD, and CAT activity in M. azedarach and L. lucidum initially increased and then declined under Pb-Zn stress ( Figure 8). The same patterns were found in Iris halophila exposed to Pb mine tailing treatments (Han et al., 2016). The SOD, POD, and CAT activities were higher in the L1 and L2 treatments than those in the C, indicating that the antioxidative system can effectively mitigate oxidative damage (Nie et al., 2016). The reduction of antioxidative enzyme activity under the highest Pb-Zn tailing treatments might be ascribed to the fact that the gene expression of SOD, POD and CAT enzymes was impacted by strong Pb-Zn stress (Hu et al., 2015). Another account for the decrease in SOD, POD, and CAT activity was that these variables were exhausted to alleviate the detrimental effects of ROS under higher Pb-Zn stress (Nie et al., 2016).
Under Pb-Zn stress, L. lucidum had a lower level of lipid peroxidation and higher activities of the antioxidative enzyme when compared to M. azedarach. As the tolerance to HMs was linked with the low level of lipid peroxidation and high activities of antioxidative enzymes (Ma et al., 2001;Sharma and Dietz, 2009), L. lucidum likely possessed a higher tolerability to HMs than M. azedarach. Additionally, less influence occurred on PSII and PSI in L. lucidum than in M. azedarach (Figure 4), which was probably attributed to the lower MDA content and higher antioxidative enzyme activities in L. lucidum (Iqbal et al., 2019). It should be noted that the TFs for Pb and Zn in stems and leaves were lower in L. lucidum than in M. azedarach, and more Pb and Zn were retained in L. lucidum roots, indicating that L. lucidum has more excellent metal exclusion strategy under Pb-Zn mine tailing treatments.

CONCLUSION
The current study showed that Pb-Zn mine tailing had a crucial negative influence on two tested plants' development and biomass production by inhibiting the chlorophyll synthesis and photosynthetic metabolism. The reduction of net photosynthetic rates in M. azedarach and L. lucidum due to HM stress was mainly caused by their biochemical limitation, including decreases of V cmax and J max . Pb-Zn stress damaged multiple locations along the photosynthetic electron transport chain. Specifically, it impaired the OEC (only in M. azedarach) and blocked electron flow of acceptor side of PSII and disturbed the PSI oxidation and reduction in both tested trees. Moreover, the increase of ROS content in both plant species was directly related to the obstruction of the electron transfer. Meanwhile, M. azedarach and L. lucidum could maintain high levels of SOD, POD, CAT, and proline contents to effectively relieve oxidative stress. The more tolerance of L. lucidum might be attributed to these facts: (1) higher RGR; (2) more accumulation of Pb-Zn in roots; (3) a less extent effect occurred on PSII and PSI activity; and (4) lower ROS and MDA content and higher antioxidative enzymes activities.

DATA AVAILABILITY STATEMENT
All datasets presented in this study are included in the article.

AUTHOR CONTRIBUTIONS
Idea and study designed: FZ. Performed the experiments: ML and HX. Wrote the paper: XH. Helped revise original paper: XC, ZH, and GW.

ACKNOWLEDGMENTS
We thank Dr. DaYong Fan for stimulating discussions and critical readings of the manuscript.