Proteomic Profiling of the Interactions of Cd/Zn in the Roots of Dwarf Polish Wheat (Triticum polonicum L.)

Cd and Zn have been shown to interact antagonistically or synergistically in various plants. In the present study of dwarf polish wheat (DPW)roots, Cd uptake was inhibited by Zn, and Zn uptake was inhibited by Cd, suggesting that Cd and Zn interact antagonistically in this plant. A study of proteomic changes showed that Cd, Zn, and Cd+Zn stresses altered the expression of 206, 303, and 190 proteins respectively. Among these, 53 proteins were altered significantly in response to all these stresses (Cd, Zn, and Cd+Zn), whereas 58, 131, and 47 proteins were altered in response to individual stresses (Cd, Zn, and Cd+Zn, respectively). Sixty-one differentially expressed proteins (DEPs) were induced in response to both Cd and Zn stresses; 33 proteins were induced in response to both Cd and Cd+Zn stresses; and 57 proteins were induced in response to both Zn and Cd+Zn stresses. These results indicate that Cd and Zn induce differential molecular responses, which result in differing interactions of Cd/Zn. A number of proteins that mainly participate in oxidation-reduction and GSH, SAM, and sucrose metabolisms were induced in response to Cd stress, but not Cd+Zn stress. This result indicates that these proteins participate in Zn inhibition of Cd uptake and ultimately cause Zn detoxification of Cd. Meanwhile, a number of proteins that mainly participate in sucrose and organic acid metabolisms and oxidation-reduction were induced in response to Zn stress but not Cd+Zn stress. This result indicates that these proteins participate in Cd inhibition of Zn uptake and ultimately cause the Cd detoxification of Zn. Other proteins induced in response to Cd, Zn, or Cd+Zn stress, participate in ribosome biogenesis, DNA metabolism, and protein folding/modification and may also participate in the differential defense mechanisms.


INTRODUCTION
Environmental toxicity from non-essential heavy metals such as cadmium (Cd), which is released from human activities and other environmental causes, is rapidly increasing (Ahsan et al., 2009). In humans, Cd causes diseases such as osteoporosis and emphysema by damaging the lungs, kidneys, and bones (Kazantzis, 2004;Straif et al., 2009). In plants, Cd damages the photosynthetic apparatus, interrupts respiratory and nitrogen metabolism, and unbalances water and nutrient uptake (Herbette et al., 2006;Balen et al., 2011), ultimately reducing biomass, causing leaf chlorosis, inhibiting root growth, and even leading to plant death (Lin et al., 2007;Yadav, 2010;Lin and Arats, 2012). Additionally, plants can accumulate high Cd contents in their edible parts, which poses a potentially major hazard to human health (Satarug et al., 2003).
Zinc (Zn), an essential metal and a cofactor of numerous plant proteins and enzymes, plays several crucial roles in protein binding, enzyme activity, transcriptional and translational regulation, and signal transduction (Broadley et al., 2007;Lin and Arats, 2012). However, excess Zn can also cause toxicity, as it can damage DNA replication and disrupt enzyme activities and protein folding and function, ultimately inducing chlorosis and inhibiting plant growth and development (Broadley et al., 2007;Lin and Arats, 2012;Schneider et al., 2013).
Cd and Zn coexist naturally in the soil. Due to their physical and chemical similarities (Chesworth, 1991), their uptake and transport in plants use similar pathways (Grant et al., 1998). Many metal transporters that transport both Cd and Zn have been identified, including AtNRAMP3 and AtNRAMP4 (Thomine et al., 2000;Lanquar et al., 2010). In response to Cd and Zn stresses, plants have developed strategies to prevent Cd-induced damage and maintain Zn homeostasis. Therefore, researchers have investigated the various synergistic and/or antagonistic interactions of Cd/Zn and found these interactions to depend on species, external bioavailable metal concentration, tissue type, and developmental stage. For example, some durum and bread wheat show antagonistic interactions of Cd/Zn in which Cd uptake is inhibited by Zn and Zn uptake is inhibited by Cd in roots, stems, and leaves (Hart et al., 2002(Hart et al., , 2005Sun et al., 2005). Conversely, some wheat under field conditions has shown synergistic interactions in which Cd and Zn uptake are promoted by each other (Nan et al., 2002).
However, previous studies on Cd/Zn interactions mainly focused on their transport and biochemical responses (Hart et al., 2002(Hart et al., , 2005Nan et al., 2002;Hassan et al., 2005;Sun et al., 2005). Although proteomic changes in response to Cd or Zn have been successfully investigated using a proteomics approach (Kieffer et al., 2008(Kieffer et al., , 2009Ahsan et al., 2009;Fukao et al., 2011;Lin and Arats, 2012;Schneider et al., 2013), the molecular mechanisms of Cd/Zn interactions are unknown, which limits our understanding of the interactions of Cd/Zn. Polish wheat (2n = 4x = 28, AABB, Triticum polonicum L.), which has low genetic similarity with T. aestivum (Wang et al., 2013;Michalcová et al., 2014), accumulates high concentrations of Zn, Fe, and Cu and therefore has attracted the interest of producers and breeders (Wiwart et al., 2013). Meanwhile, dwarf polish wheat (DPW), collected from Tulufan, Xingjiang, China, shows high tolerance to Cd and Zn because its growth is not affected by the accumulation of high concentrations of these metals in seedlings (Wang X. et al., in press). However, molecular responses to Cd and Zn remain unknown. Since DPW accumulates high concentrations of Cd and Zn in seedlings, it is a useful system for studying Cd/Zn interactions.
The purposes of this study are to understand molecular responses to Cd and Zn stresses, to investigate Cd/Zn interactions in DPW seedlings and to understand the molecular mechanisms of Cd/Zn interactions in DPW roots using isobaric tags for relative and absolute quantification (iTRAQ). iTRAQ is a highthroughput proteomic technology (Karp et al., 2010) that has been successfully used to reveal plant responses to heavy metals (Ahsan et al., 2009;Fukao et al., 2011).

Plant Material and Growth Conditions
DPW seeds were sterilized with 1% NaOCl and germinated in the dark for 5 days. The seedlings were cultured in full Hoagland nutrient solution in a growth chamber at 25 • C with a 16 hlight/8 h-dark cycle. At the two leaf stage, the seedlings were treated with null (CK), 40 µM CdSO 4 (Cd), 800 µM ZnCl 2 (Zn), or 40 µM CdSO 4 + 800 µM ZnCl 2 (Cd+Zn). Two days after treatments, the roots (three biological replications, each replication including 15 plants) were washed with 0.1 µM EDTA and ddH 2 O, snap frozen in liquid nitrogen and stored at −80 • C for iTRAQ analysis. Other roots and leaves were dried for 2 days at 80 • C for measuring metal concentrations.

Measurement of Cd and Zn Concentrations
Cd and Zn concentrations were measured as described by Wang et al. (2014). Briefly, the dried roots and leaves were ground to particle powders. Then, 0.2 g of powder was digested using concentrated sulfuric acid and hydrogen peroxide at 320 • C and then diluted to 50 ml. Metal concentrations were then determined using an atomic absorption spectrometer, Analyst 400 (PerkinElmer, CT, USA). Standard solutions of Cd and Zn were purchased from Fisher Scientific Ltd. (China). All data and figures were analyzed (t-test was conducted for the statistical analysis) and drawn using Sigmaplot 12.0.

Total Protein Extraction
Roots (two randomly selected biological replications) with 0.1 mg of polyvinylpyrrolidone (PVPP) were ground into powders using liquid nitrogen and then homogenized in Tris-phenol (pH 8.0) and protein extraction buffer (0.7 M sucrose, 0.1 M KCl, 50 mM EDTA, 0.5 M Tris, pH 7.5, 2% β-mercaptoethanol, and 1 mM PMSF). After centrifuging for 20 min at 6000 rpm, the supernatants were collected and re-purified using protein extraction buffer. Proteins were precipitated using ammonium acetate methanol and then washed with methanol and acetone. Finally, protein samples were diluted using RIPA reagent, and protein concentrations were measured using a BCA Assay Kit (Biotech).
iTRAQ Labeling and LC-MS Analysis iTRAQ labeling was performed according to Wu et al. (2013) with modifications. Briefly, 200 µg of protein from each sample (two biological replications) was reduced, alkylated and then subjected to tryptic hydrolysis. iTRAQ labeling was performed using an iTRAQ R reagents-8plex Kit (Applied Biosystems). Peptides of CK, Cd, Zn, and Cd+Zn samples were labeled singly with the iTRAQ reporters 113, 114, 115, and 116, respectively. LC-MS (TripleTOF5600, Applied Biosystems) analysis was performed as described by Wu et al. (2013).

Protein Identification and Quantification
Protein identification and relative quantification were also performed according to Wu et al. (2013). Protein Pilot software v. 4.0 (Applied Biosystems) was used to convert the raw data (.wiff) into peak lists (.mgf). Each MS/MS spectrum was searched against the protein database Uniprot-147389. The search parameters were as follows: Paragon method: iTRAQ-8plex, Cys alkylation: MMTS, Digestion: Trypsin, Instrument: TripleTOF 5600, ID focus: Biological modifications and Amino acids substitutions, Detected Protein Threshold [Unused ProtScore (Confidence)]: ≥ 1.3, Competitor Error Margin (ProtScore): 2.0, and No. Distinct Peptides (Confidence): ≥ 95%. The tolerances were specified as ± 0.05 Da for peptides and ± 0.05 Da for MS/MS fragments. The relative abundance (fold-change ratios of differential abundance between labeled samples), P-value, error factor, lower confidence interval and upper confidence interval were calculated using the ProteinPilot software. Proteins containing at least two distinct peptides and fold change ratios ≥ 1.5 or ≤ 0.67 were considered as more abundant or less abundant proteins, respectively.

Metal Concentrations
No Cd was detected in CK (control) samples ( Figure 1A). Two days after treatment, the Cd concentration in roots treated with Cd (752.55 ± 6.51 mg/Kg) was significantly higher (P < 0.01) than that in roots treated with Cd+Zn (76.75 ± 3.312 mg/Kg; Figure 1A). Meanwhile, the Cd concentration in leaves under Cd stress (40.87 ± 3.69 mg/Kg) was also significantly higher (P < 0.01) than that in leaves under Cd+Zn stress (9.20 ± 1.24 mg/Kg; Figure 1A). These results indicate that Zn inhibits Cd uptake in roots as well as its transport from roots to shoots.
Zn concentrations in the roots were always higher than those in the leaves ( Figure 1B). Zn concentrations in leaves were similar between Zn (139.26 ± 32.12 mg/Kg) and Cd+Zn (147.00 ± 20.15 mg/Kg) stresses ( Figure 1B). In roots treated with Zn (675.36 ± 41.67 mg/Kg), the Zn concentration was significantly higher (P < 0.01) than that in roots treated with Cd+Zn (557.63 ± 26.30 mg/Kg; Figure 1B). These results suggest that Cd inhibits Zn uptake in roots but does not affect its transport from root to shoot.

A Total of 432 Proteins Were Altered by Cd, Zn, or Cd+Zn Stresses
A total of 960 proteins with one or more distinct peptides and an Unused ProtScore ≥ 1.3 (with a peptide confidence ≥ 95%) (Data Sheet 1) were identified from the protein database Uniprot-147389. Compared with null, the expression levels of 206 (Data Sheet 2), 303 (Data Sheet 3), and 190 (Data Sheet 4) proteins were altered by Cd, Zn, and Cd+Zn stresses, respectively. Further analysis indicated that these proteins could be grouped into seven sub-groups (Table 1, Figure 2).  The relative abundances of 53 proteins were altered significantly by all three stresses (Cd, Zn, and Cd+Zn; Figure 2, Data Sheet 5). Among these, 13 noteworthy proteins participated in either sucrose metabolism (5 proteins), glutathione (GSH) metabolism (5 proteins), or the oxidation-reduction process (3 proteins; Table 1). However, the relative abundances of other proteins were differentially altered by Cd, Zn, and combined Cd+Zn stresses (Table 1). For example, the relative abundance of glucose-6-phosphate isomerase (protein 554) was increased by Cd stress but was decreased by both Zn and Cd+Zn stresses. Contrary results were observed for lactoyglutathione lyase (protein 424), as its relative abundance was decreased by Cd stress but was increased by both Zn and Cd+Zn stresses. Further, the relative abundance of cytochrome c oxidase subunit 6B (protein 616) was increased by both Cd and Cd+Zn stresses but was decreased by Zn stress (Table 1). Thus, our analysis revealed differential molecular responses to Cd, Zn, and Cd+Zn stresses.   FIGURE 2 | Numbers of altered proteins which were classified into differentially interactive groups in response to Cd, Zn, and Cd+Zn stresses.

Proteins Were Induced only in Response to Cd Stress
We identified 58 proteins whose relative abundances were induced only in response to Cd stress (Figure 2, Data Sheet 6).
Of these, the relative abundances of 23 proteins were increased, and those of 35 proteins were decreased (Figure 2, Data Sheet 6). These proteins were not induced by either Zn or combined Cd+Zn stress, which suggests that they might participate in Zn inhibition of Cd uptake and transport. Among the 58 proteins we identified were 18 noteworthy proteins that participated in the oxidation-reduction process (4 down and 4 up), GSH metabolism (2 up and 1 down), sucrose metabolism (2 up), calcium (Ca) metabolism (3 up), or S-adenosyl-L-methionine (SAM) metabolism (1 down and 1 up; Table 1).

Proteins Were Induced only in Response to Zn Stress
The relative abundances of 131 proteins were induced only in response to Zn stress (46 up and 85 down; Figure 2, Data Sheet 7). That these proteins were not induced in response to either Cd or combined Cd+Zn stress suggests that they might participate in Cd inhibition of Zn uptake. Among the 131 proteins we identified, we classified 26 noteworthy DEPs into four functional groups (Table 1): Sucrose metabolism (6 down and 2 up), organic acid metabolism (4 down and 2 up), the oxidationreduction process (4 up and 7 down), and cation transport (1 down).

Proteins Were Induced in Response to both Cd and Zn Stresses
We observed 61 DEPs whose relative abundances were altered in response to both Cd and Zn stresses (Figure 2, Data Sheet 8).
However, under Cd stress, the relative abundances of 24 proteins were increased and 37 were decreased, whereas under Zn stress, 29 were increased and 32 were decreased. We also observed 12 proteins whose relative abundances were altered inversely in response to Cd and Zn stress (Data Sheet 8, marked by yellow). These proteins were not induced in response to combined Cd+Zn stress, indicating that they might be involved in the mutual inhibition of Cd/Zn. Among the 61 proteins we identified were 5 noteworthy proteins that we divided into three functional pathways ( Table 1): The oxidation-reduction process (2), sucrose metabolism (1) and SAM metabolism (2).

Proteins Were Induced in Response to both Cd and Cd+Zn Stresses
We identified 33 proteins whose relative abundances were induced in response to both Cd and combined Cd+Zn stresses (Figure 2, Data Sheet 9). Under Cd stress, 13 proteins were more abundant, and 20 were less abundant, whereas under Cd+Zn stress, 18 were more abundant, and 15 were less abundant. These results indicate that Cd and Cd+Zn stresses induce differential molecular responses. We also identified 11 proteins whose relative abundances were altered inversely in response to Cd and Cd+Zn, including nicotianamine synthase 2 (NAS; Data Sheet 9, marked by yellow). That these proteins were not induced in response to Zn stress suggests that they might participate in Cd detoxification of Zn. Among the 33 proteins we identified, 2 proteins were key enzymes in SAM metabolism, 1 was involved in Ca metabolism, and 1 participated in sucrose metabolism ( Table 1).

Proteins Were Induced in Response to both Zn and Cd+Zn Stresses
We identified 57 proteins whose relative abundances were induced in response to both Zn and combined Cd+Zn stresses (Figure 2, Data Sheet 10). Under Zn stress, 28 proteins were more abundant, and 29 were less abundant, whereas under Cd+Zn stress, 39 were more abundant, and only 18 were less abundant (Data Sheet 10). These results indicate that Zn and Cd+Zn stresses induce differential molecular responses. We also identified 28 proteins whose relative abundances were altered inversely in response to Zn and Cd+Zn, including glutaredoxin-C8 and glucan 1,3-beta-glucosidase, (Data Sheet 10, marked by yellow). These proteins were not induced in response to Cd stress, which suggests that they might participate in Zn detoxification of Cd. Among the 57 proteins we identified, 4 proteins functioned in GSH metabolism, 3 were peroxidases involved in the oxidationreduction process, and 1 was involved in sucrose metabolism ( Table 1).

Proteins Were Induced only in Response to Cd+Zn Stress
We identified 47 proteins whose relative abundances were induced only in response to Cd+Zn stress (Figure 2,  Data Sheet 11). These proteins were not induced in response to either Cd or Zn stress alone, which suggests that the molecular response induced by Cd+Zn stress differs from that induced by Cd and Zn stresses individually. Of the 47 proteins we identified, 24 proteins were more abundant, and 23 proteins were less abundant (Data Sheet 11). Among these, we identified 11 noteworthy proteins that functioned in Ca metabolism, the oxidation-reduction process, organic acid metabolism, and sucrose metabolism (Table 1).

DISCUSSION
Interactions between Cd and Zn have previously been shown to be antagonistic and/or synergistic in various plants (Hart et al., 2002(Hart et al., , 2005Sun et al., 2005;Tkalec et al., 2014). In the present study, Cd uptake was inhibited by Zn and Zn uptake was inhibited by Cd in DPW roots (Figure 1). Cd transport from root to shoot was inhibited by Zn ( Figure 1A) but was promoted by Zn after 5 days after treatment (unpublished data). Meanwhile, Zn transport from root to shoot was not affected by Cd ( Figure 1B). These results indicate that Cd and Zinc interact antagonistically in DPW seedlings, as previously reported in bread and durum wheat (Hart et al., 2002(Hart et al., , 2005Sun et al., 2005) and unlike the synergistic interactions that have been reported in other wheat under field conditions (Nan et al., 2002). Proteomic changes in the roots implicated several proteins in the antagonistic interactions of Cd/Zn (Data Sheets 1-11). Two days after treatment, the relative abundances of 206 (Data Sheet 2), 303 (Data Sheet 3), and 190 (Data Sheet 4) proteins were induced in response to Cd, Zn, and Cd+Zn stresses, respectively (Figure 2). Among these, 53 proteins were induced in response to all three treatments, and 58, 131, and 47 proteins were induced in response to only Cd, Zn, or Cd+Zn stresses, respectively (Figure 2). We grouped these proteins into different interactions of Cd/Zn (Figure 2). Our results indicate that although Cd and Zn have similar physical and chemical properties (Chesworth, 1991) and pathways for uptake (Grant et al., 1998), they induce differential molecular responses (Lin and Arats, 2012), which result in the antagonistic interactions of Cd/Zn in DPW roots (Figure 1) and the high tolerance of DPW to Cd and Zn toxicity (Wang X. et al., in press). Some proteins identified in this study that are involved in noteworthy processes are discussed below.
To overcome oxidative toxicity caused by heavy metal stresses (Ranieri et al., 2005;Lin et al., 2007;Kieffer et al., 2008;Di Baccio et al., 2011;Zeng et al., 2011), plants utilize an effective antioxidant system that protects their cells against oxidative damage (Kieffer et al., 2008;Di Baccio et al., 2011) by inducing the expression of oxidation-reduction-related proteins (Lin et al., 2007;Kieffer et al., 2008Kieffer et al., , 2009Di Baccio et al., 2011;Zeng et al., 2011;Schneider et al., 2013). In this study, 31 oxidationreduction-related proteins were observed ( Table 1). Of these, 8, 11, and 4 proteins were altered in response to Cd, Zn, and Cd+Zn stresses, respectively ( Table 1). These results suggest that Cd, Zn, and Cd+Zn stresses cause differential oxidative threats which are detoxified through the induction of different oxidation-reduction-related proteins. Conversely, 8 Cd-induced proteins, 11 Zn-induced proteins, and 2 proteins induced by both Cd and Zn stresses were not induced in response to combined Cd+Zn stress, which suggests that the oxidative threats caused by Cd and Zn stresses are not the same as those caused by Cd+Zn stress. These results indicate that Cd and Zn detoxify each other in combined Cd+Zn stress, resulting in their uptakes being inhibited by each other. As described in previous studies (Kieffer et al., 2008(Kieffer et al., , 2009Schneider et al., 2013), Cd and Zn induced a greater abundance of some oxidative stress-related proteins but also induced a lower abundance of other oxidative stressrelated proteins ( Table 1). Among these, 3 proteins were induced by all 3 stresses (Cd, Zn, and Cd+Zn) ( Table 1), including ascorbate peroxidase (protein 285), L-ascorbate peroxidase 2 (protein 573), and peroxisome type ascorbate peroxidase (protein 511), which are key peroxide detoxification enzymes (Raven et al., 2004). These results suggest that ascorbate mediates Cd-and Zn-induced oxidative stress in plants (Kieffer et al., 2008).
As a precursor of GSH, SAM plays important roles in protecting against Cd stress-induced reactive oxygen species (ROS) (Noriega et al., 2007). In the present study, protein levels of serine hydroxymethyltransferase (SHMT) and spermidine synthase 1, key enzymes in SAM metabolism, were altered in response to both Cd and Zn stresses, suggesting that SAM plays important roles in protecting against these stresses. Sadenosylmethionine synthase (SAMS) synthesizes SAM, which is a precursor of nicotianamine (NA) (Schneider et al., 2013). The protein level of nicotianamine synthase 2 (NAS), which synthesizes nicotianamine (NA) from SAM, increased in response to Cd stress. NA is an essential compound for cell-to-cell transport of Zn, Fe, and Cu (Takahashi et al., 2003;Klatte et al., 2009). However, the abundances of both SAMS and NAS decreased in response to combined Cd+Zn stress, whereas a previous report in N. caerulescens showed increased SAMS and NAS levels in response to Zn stress (Schneider et al., 2013). Our results partially illustrate Cd inhibition of Zn uptake. Cd stress also causes the lignification of roots (Finger-Teixeira et al., 2010). SAM provides the methyl donor to caffeic acid 3-Omethyltransferase (COMT) in lignin biosynthesis (Wang Y. et al., 2016). COMT levels were increased only in response to Cd stress ( Table 1), suggesting that Cd also causes root lignification.
Cellulose and pectic polysaccharides are major components of the plant cell wall (Cosgrove, 2005), which can be modified by some heavy metals. For example, Cd enhances the contents of glucose and polysaccharides in cell walls . Further, exogenous glucose alleviates Cd toxicity by increasing Cd fixation in root cell walls (Shi et al., 2015). In this study, several sucrose metabolism-related proteins were induced by Cd, Zn, or Cd+Zn stress ( Table 1), suggesting that glucose and/or polysaccharides participate in Cd and Zn fixation, exclusion or sequestration in root cell walls Shi et al., 2015). However, 8 sucrose metabolism-related proteins were induced in response to Zn stress but not combined Cd+Zn stress (Table 1), which suggests that Cd detoxifies excess Zn by inhibiting its uptake, resulting in Cd-induced inhibition of Zn modification of sucrose metabolism.
Zn stress also affects the expression of P-type ATPases and several other metal transporters (Schneider et al., 2013). P-type ATPases, such as AtHMA4 from Arabidopsis, GmHMA3 from soybean and AhHMA3 from A. halleri, have Zn uptake activity (Becher et al., 2004;Hussain et al., 2004;Wang et al., 2012). Further, AtHMA2, 3, and 4 have been shown to transport Zn from root to shoot (Williams and Mills, 2005). In this study, the abundance of a P-type proton pump ATPase (protein 679) decreased in response to Zn stress but not combined Cd+Zn stress (Table 1), a result which contradicts previous work in N. caerulescens showing increased abundance of two P-type ATPases in response to Zn stress (Schneider et al., 2013). However, we found that other Zn-induced metal transporters were not observed 2 days after treatment, but their transcripts were regulated by Cd, Zn, and Cd+Zn 5 days after treatment (unpublished data).
Finally, as reported by previous studies (Di Baccio et al., 2011;Zeng et al., 2011), proteins that were similarly or differentially induced in response to Cd, Zn, and/or Cd+Zn stresses also participated in other processes, including ribosome biogenesis, DNA metabolism, protein folding/modification (all SFiles), suggesting that these proteins might contribute to differential defense mechanisms against these stresses (Zeng et al., 2011).

CONCLUSION
Taken together, our results indicate that Cd and Zn interact antagonistically in DPW seedlings. Although 206, 303, and 190 proteins were induced in response to Cd, Zn, and Cd+Zn stresses, respectively, only 53 proteins were induced in response to all three stresses. 58, 131, and 47 proteins were induced only in response to Cd, Zn, and Cd+Zn stresses, respectively (Figure 2). These proteins could be divided into groups that resulted in different Cd/Zn interactions. Our results suggest that Zn and Cd stresses cause differential molecular responses in DPW. Under these stresses, oxidative stress-related proteins, metal chelators, metabolism-related proteins, sucrose metabolismrelated proteins, and metal transporters are differentially induced to participate in metal detoxification, which ultimately causes antagonistic interactions and enhanced tolerance of Cd and Zn.

AUTHOR CONTRIBUTIONS
YW, XW, XX, and YZ conceived and designed research, and wrote the manuscript. YW, XW, XX, CW, FP, and RW conducted experiments. YW, XW, JZ, HK, XF, LS, and HZ analyzed data. All authors read and approved the manuscript.

ACKNOWLEDGMENTS
The authors thank the National Natural Science Foundation of China (No. 31301349, 31270243, and 31470305), Bureau of Science and Technology and Bureau of Education of Sichuan Province, China. We would like to thank Lu Gao (Guangzhou Fitgene Biotechnology Co., Ltd) for useful advice and discussion.

SUPPLEMENTARY MATERIAL
The