Next Article in Journal
Comparison of Policosanol Profiles of the Sprouts of Wheat Mutant Lines and the Effect of Differential LED Lights on Selected Lines
Next Article in Special Issue
Genome-Wide Identification and Expression Analysis of the TCP Gene Family Related to Developmental and Abiotic Stress in Ginger
Previous Article in Journal
Continuous Cropping Inhibits Photosynthesis of Polygonatum odoratum
Previous Article in Special Issue
Quantification of Glucose Metabolism and Nitrogen Utilization in Two Brassicaceae Species under Bicarbonate and Variable Ammonium Soil Conditions
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 12
Issue 19
10.3390/plants12193375
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

MeGLYI-13, a Glyoxalase I Gene in Cassava, Enhances the Tolerance of Yeast and Arabidopsis to Zinc and Copper Stresses

by
Ruimei Li
1,2,3,4,†,
Fenlian Tang
1,3,†,
Yannian Che
1,3,
Alisdair R. Fernie
4,
Qin Zhou
1,3,
Zhongping Ding
1,3,
Yuan Yao
1,2,
Jiao Liu
1,2,
Yajie Wang
1,2,
Xinwen Hu
3 and
Jianchun Guo
1,2,3,*
1
Key Laboratory of Biology and Genetic Resources of Tropical Crops, Institute of Tropical Bioscience and Biotechnology, Chinese Academy of Tropical Agricultural Sciences, Haikou 571101, China
2
Key Laboratory for Biology and Genetic Resources of Tropical Crops of Hainan Province, Hainan Institute for Tropical Agricultural Resources, Haikou 571101, China
3
College of Tropical Crops, Hainan University, Haikou 570228, China
4
Root Biology and Symbiosis, Max-Planck-Institute of Molecular Plant Physiology, Am Muhlenberg 1, 14476 Potsdam-Golm, Germany
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Plants 2023, 12(19), 3375; https://doi.org/10.3390/plants12193375
Submission received: 13 August 2023 / Revised: 19 September 2023 / Accepted: 22 September 2023 / Published: 25 September 2023
(This article belongs to the Special Issue Plant Stress Physiology and Molecular Biology)
Browse Figures
Versions Notes

Abstract

:
Although zinc and copper are the two essential nutrients necessary for plant growth, their excessive accumulation in soil not only causes environmental pollution but also seriously threatens human health and inhibits plant growth. The breeding of plants with novel zinc or copper toxicity tolerance capacities represents one strategy to address this problem. Glyoxalase I (GLYI) family genes have previously been suggested to be involved in the resistance to a wide range of abiotic stresses, including those invoked by heavy metals. Here, a MeGLYI-13 gene cloned from a cassava SC8 cultivar was characterized with regard to its potential ability in resistance to zinc or copper stresses. Sequence alignment indicated that MeGLYI-13 exhibits sequence differences between genotypes. Transient expression analysis revealed the nuclear localization of MeGLYI-13. A nuclear localization signal (NLS) was found in its C-terminal region. There are 12 Zn2+ binding sites and 14 Cu2+ binding sites predicted by the MIB tool, of which six binding sites were shared by Zn2+ and Cu2+. The overexpression of MeGLYI-13 enhanced both the zinc and copper toxicity tolerances of transformed yeast cells and Arabidopsis seedlings. Taken together, our study shows the ability of the MeGLYI-13 gene to resist zinc and copper toxicity, which provides genetic resources for the future breeding of plants resistant to zinc and copper and potentially other heavy metals.

1. Introduction

Given that heavy metal pollution has had, and will continue to have, a serious impact on the environment and plant, animal and microbial life, leading to the destruction of ecosystems, heavy metal pollution is attracting increasing worldwide attention [1]. There are many reasons for heavy metal pollution in farmland soil during agricultural production, including the use of wastewater for soil irrigation, the application of fertilizers, the spraying of pesticides, herbicides and so on [2]. Especially, the frequent and extensive use of various zinc (Zn)- and copper (Cu)-containing fertilizers and fungicides has resulted in the continuous accumulation of Zn and Cu in farmland soils [3]. Both Zn and Cu are essential micronutrients for plant growth and development. However, when their concentrations are too high, they result in toxicity to organisms [4,5].
Zinc is an essential nutrient for plants, because it is an important component of key protein structures and catalytic sites [6,7]. In plants, zinc has a significant effect on reactive oxygen metabolism and regulating defense reactions and signals [8,9]. Zinc deficiency reduces leaf biomass, increases ROS accumulation and affects plant hormone levels and TCA-related enzyme activities [9]. The leaves of Zn2+-deficient plants were notably yellow, and the levels of iron in their branches increased [6]. Zinc deficiency results in stunted growth, reduced tillers and low productivity of rice plants, with brown spots and stripes on the leaves [10]. Moreover, concerning zinc deficiency, some enzymes such as alkaline phosphatase (ALP), alcohol dehydrogenase (ADH) and polyphenol oxidase (PPO) activity decreased, while guaiacol peroxidase (GPX) activity increased [11]. Zinc toxicity also induces reactive oxygen species (ROS) accumulation [9,12]. Zinc toxicity additionally significantly reduces plant biomass, inhibits the accumulation of other metal ions, reduces the photosynthetic rate and renders plants more susceptible to disease [12].
Copper is an essential nutrient for the normal metabolism of aerobic life. As a cofactor of various redox enzymes, copper plays a key role in redox catalysis. Copper is also involved in the signal regulation, protein transport and metabolism of other metals in cells [13,14,15]. Copper deficiency in the growth environment additionally results in adverse effects on plant growth and development, such as inhibiting vegetative growth and reproductive development, causing leaf phenotypic distortion and color change, affecting pollen and embryo development, reducing crop yield and affecting cell wall formation and water transport [13]. Copper deficiency also affects photosynthesis and respiration in plants [16,17]. However, more copper is not always better, and more than the required levels can be toxic to the plant. Plant exposure to high concentrations of copper results in physiological disorder and dwarf phenotypes [18,19,20]. The excessive copper also inhibits root growth and alters root morphology [21]. The level of toxic ROS increases; the cell structure is damaged and the photochemical efficiency, chlorophyll content, stomatal conductance and transpiration rate are also affected in plants exposed to too-high levels of copper [18,19].
To cope with the damage of abiotic stress, plants have evolved a series of complex plant protection mechanisms to mitigate external damage. Some recent studies have found that the glyoxalase system is involved in plant protection against heavy metals by removing the excess methylglyoxal (MG) [22,23,24]. Glyoxalase I (GLYI) is the first enzyme in the MG detoxicant pathway [25]. The increased Cu tolerance in rice is accompanied by the activation of glyoxalase I and II and the decrease in the methylglyoxal content [24]. Similarly, it has been found that, in mung beans, glyoxalase was upregulated to reduce the Cu phytotoxicity [22]. The heterologous expression of the wheat TaGLYI gene in tobacco has improved the tolerance of transgenic tobacco to heavy metal Zn2+ stress [26]. Transgenic tobacco heterologously expressing the daylily HfGLYI-1 is tolerant to abiotic stress at high concentrations of ZnCl2 [27]. Glyoxalase is also involved in the defense of plants against other heavy metal stress, such as nickel (Ni), chromium (Cr), iron (Fe), arsenic (As), vanadium and cadmium (Cd) [23,26,28,29,30,31,32,33,34,35].
Cassava (Manihot esculenta) is a woody shrub that accumulates considerable quantities of starch in its tuberous roots. It is not only a food crop to ensure world food security but also as a high-quality raw material for the production of clean bioenergy. A previous study indicated that cassava can be grown in heavy metal (Cd, Cu, Pb and Zn)-contaminated fields without affecting its biomass production [36]. This indicated that cassava has a strong resistance to heavy metal stress. Therefore, mining heavy metal resistance genes in cassava is of great value to improve heavy metal resistance in cassava and other plants. We have previously identified 19 putative GLYI genes from cassava, in which MeGLYI-13 was found to play important roles in the resistance to Fe3+ stress [33]. However, whether the MeGLYI-13 gene is involved in other heavy metal stresses, particularly Zn and Cu toxicity, remains unclear. Further study of MeGLYI-13 sequences, structural characteristics, protein subcellular localization and tolerance to Zn and Cu stress is of great significance to broaden our understanding of the function of the MeGLYI-13 gene and its potential application in modern agriculture.

2. Results

2.1. MeGLYI-13 Gene Sequence Analysis

The MeGLYI-13 gene was cloned from the genome of the SC8 cassava cultivar using the Manes.12G054100 sequence retrieved from the AM560-2 cultivar as a reference (Manihot esculenta v6.1 dataset). Its cDNA length is 495 bp, and it encodes a protein of 164 amino acids in length. It has 99.3% similarity with the Manes.12G054100 sequence from cassava variety AM560-2 in the cassava genome database. There is one base difference between MeGLYI-13 in SC8 and Manes.12G054100 in AM560-2, which is “205C” in SC8 and “205T” in AM560-2. At the protein level, this corresponds to ”69P” in SC8 but “69S” in AM560-2 (Figure S1). The MeGLYI-13 gene from SC8 was registered in NCBI GenBank (MW413417).

2.2. Prediction of Zn2+ and Cu2+ Binding Sites in MeGLYI-13 Protein Sequence

In order to understand the potential of MeGLYI-13 to bind zinc and copper ions, we predicted the Zn2+/Cu2+ ion binding sites in its protein sequence. The results showed that there were 12 Zn2+ binding sites and 14 Cu2+ binding sites in the MeGLYI-13 amino acid sequence. In addition, six sites: 57LYS, 66HIS, 83ASP, 106 ASP, 108GLU and 127GLU bound to both Zn2+ and Cu2+ (Figure 1).

2.3. MeGLYI-13 Protein Was Localized in the Nucleus

With the intention of assessing the localization of the MeGLYI-13 protein, pCambia1300-MeGLYI-13-GFP and pCambia1300-GFP plasmids were transformed into onion epidermal cells by Agrobacterium tumefaciens infection. The green fluorescent signal could be clearly observed in the cytoplasm and nucleus following infection with the control vector, whereas the fluorescent signal could only be detected in the nucleus after infection with the pCambia1300-MeGLYI-13-GFP vector. The nuclear localization was confirmed by DAPI staining (Figure 2).

2.4. The Ability of Yeast to Tolerate Cu2+ and Zn2+ Toxicity Is Enhanced by MeGLYI-13

The growth characteristics of INVSc1 yeast, containing the pYES2-MeGLYI-13 plasmid or pYES2 control vector, respectively, treated by Zn2+ or Cu2+, were compared. The pYES2-MeGLYI-13-transformed yeast showed better growth ability than the control yeast after the Zn2+ and Cu2+ treatments (Figure 3). No colonies appeared after 10−5 dilution of the yeast transformed by the control vector after the Zn2+ treatment. However, under the same conditions, the pYES2-MeGLYI-13-transformed yeast could still grow after 10−6 dilution. After treatment with Cu2+, the growth of the 10−3 dilution control yeast was completely suppressed, whereas the growth of the 10−6 dilution pYES2-MeGLYI-13 yeast was not completely suppressed (Figure 3). These results indicated that overexpression of the MeGLYI-13 gene raised the resistance of INVSc1 yeast to both Zn2+ and Cu2+ toxicity.

2.5. Expression of MeGLYI-13 in Arabidopsis

Arabidopsis was transformed with the pCambia1300-MeGLYI-13-GFP vector. A total of six transgenic lines were generated under hygromycin selection, and the presence of MeGLYI-13 was confirmed by PCR (Figure 4A). The six transgenic lines were named OE1, OE2, OE3, OE4, OE5 and OE6, respectively. To determine the MeGLYI-13 gene expression in the transgenic lines, qRT-PCR was performed. The results revealed that the expression level of the MeGLYI-13 gene varied in different transgenic lines. The transcriptional abundance of the MeGLYI-13 gene was highest in the OE1 and OE5 lines, followed by OE2. The transcriptional abundance of MeGLYI-13 was lowest in the OE4 transgenic line. There was no MeGLYI-13 gene expression in the WT Arabidopsis (Figure 4B). Because the vector had a GFP tag to further verify the MeGLYI-13 transgenic Arabidopsis, the seedlings of six transgenic Arabidopsis lines were exposed to a LUYOR-3415RG dual-wavelength fluorescent protein excitation light source. The results revealed that the leaves of the WT and OE4 with low MeGLYI-13 expression were red, while the leaves and roots of the other five transgenic lines showed a strong green fluorescence. This indicated that the MeGLYI-13 gene was successfully translated in lines OE1, OE2, OE3, OE5 and OE6 (Figure 4C).

2.6. MeGLYI-13-Overexpressing Arabidopsis Showed Better Tolerance to Zinc Shortage and Toxity

To test whether MeGLYI-13 is involved in the tolerance to Zn2+ or Cu2+ toxicity, three independent lines—namely, OE1, OE2 and OE5—of MeGLYI-13 transgenic Arabidopsis were selected for study. The wild type and MeGLYI-13-overexpressing in vitro-grown plants were separately treated with different concentrations of Zn2+ and Cu2+. This study revealed that, when grown on normal 1/2 MS medium (15 μmol/L Zn2+ and 0.05 μmol/L Cu2+), there was no significant difference in the phenotypes of the WT and MeGLYI-13-overexpressing lines (Figure 5A,B). However, following treatment with Zn2+, we found that, when there was no Zn2+ or only low levels of Zn2+ (4 μmol/L) in the medium, the growth of WT was inhibited compared to the MeGLYI-13-overexpressing lines. Indeed, the leaves were smaller, and the taproots and lateral roots were shorter than that of the transgenic lines. After exposure to high levels of zinc (60 μmol/L), the taproots of the WT were obviously shorter than those of the MeGLYI-13-overexpressing lines. The lateral roots of the WT were hard to cultivate, whereas those of the MeGLYI-13-overexpressing lines grew well. Meanwhile, the leaves of the WT were also smaller than those of the MeGLYI-13-overexpressing lines (Figure 5C,D). The above consequences indicated that the overexpression of MeGLYI-13 enhanced the survivability towards both the deficiency and toxicity of zinc.

2.7. MeGLYI-13-Overexpressing Arabidopsis Showed Better Tolerance to Copper Toxity

When grown on zero Cu2+ medium, the taproot length of the MeGLYI-13-overexpressing lines was distinctly shorter than that of the WT. When grown on the medium containing 10 μmol/L Cu2+, by contrast, the roots of the MeGLYI-13-overexpressing lines became longer than those of the WT, and the leaves were visibly larger. At high concentrations (60 μmol/L), however, the growth of both the WT and MeGLYI-13-overexpressing lines were significantly inhibited. However, the MeGLYI-13-overexpressing lines still grew better than the WT, with larger leaves, longer taproots and more root hair (Figure 5E,F). Therefore, it is speculated that the MeGLYI-13 gene is involved in Cu2+ tolerance.

3. Discussion

Zinc and copper are two essential metal elements for plant growth and development. The lack of these two elements in the environment can inhibit the growth of crops, resulting in, amongst other consequences, reduced yields and quality. However, the presence of excessive amounts of these two elements can also be toxic to plants [4]. In recent years, the extensive use of zinc-containing and copper-containing pesticides and fertilizers, as well as the development of zinc and copper mining, means that the pollution of zinc and copper in soil has become increasingly serious [3,37]. In order to solve the damage caused by soil metal pollution to plants, people have tried many methods. One approach is to use soil amendments to limit the migration of metal ions so as not to affect plant growth [20]. Another approach is to increase a plant’s tolerance to metal ion toxicity by spraying it with chemicals [5]. In addition, another method is genetic engineering to improve the metal ion tolerance of plants. The improved plants can efficiently absorb metal ions in the soil and reduce the excessive accumulation of metal ions. Therefore, it is important to cultivate plants resistant to zinc or copper toxicity for the phytoremediation of soil contaminated with copper and zinc. In this respect, a two-phase strategy is required: firstly, the generation of plants that can cope with heavy metals and produce sustainable yields and, secondly, finding a mechanism by which heavy metals are stored in the part of the plants not used for human and/or animal consumption [38,39]. Cassava is a kind of woody shrub, with its tuberous roots full of starch, which is a good resource for producing bioenergy. A previous study indicated that cassava could grow well in Pb-, Cd-, Zn- and Cu-polluted soils. The soil mobile metal and bioavailable metal content decreased significantly after cassava planting [36]. Thus, cassava can be used as an important plant resource for phytoremediation—on the one hand, to improve the polluted soil environment, and on the other hand, to produce bioenergy raw materials. Furthermore, the key genes for heavy metal resistance in cassava can also be used to improve heavy metal resistance in plants.
To combat the stress of metal toxicity, plants have evolved various protective mechanisms. Since metal toxicity stress can induce cytotoxic MG production in plants, plants have evolved a glyoxalase system to detoxify intracellular MG. Some reports have shown that the physiological mechanism of chemical measures to improve heavy metal resistance in plants involves the mobilization of a glyoxalase system [22,23,29,32,34,35]. Glyoxalase I (GLYI) is in the first step of the glutathione-dependent glyoxalase pathway to scavenge excess MG into S-D-lactoylglutathione. GLYI family genes have been identified and characterized in many species, such as Arabidopsis [40], rice [19], sugarcane [41], Chinese cabbage [42], grapes [43] and date palms [44]. We have also identified 19 GLYI family gene members in cassava [33]. The function of some GLYI genes in abiotic stress has been continuously demonstrated by experiments. For example, a sugar beet M14 glyoxalase I gene could increase E. coli tolerance to MG and increase tobacco tolerance to MG, H2O2, salt and mannitol [45], whereas an OsGLYI gene from rice was proven to improve rice tolerance to NaCl, ZnCl2 and mannitol [46]. Similarly, a sugarcane GLYI gene (SoGloI) improved E. coli tolerance to NaCl, CuCl2, CdCl2 and ZnSO4 stress [47]. Mutants of nuclear-located AtGLYI-2 displayed an increase in the level of MG, a decrease in GLYI activity and were more susceptible to salt stress [19], while nuclear-located OsGLYI-8 could recue the sensitivity of the atglyi-2 mutant to MG and salt stress [19]. In view of previous reports on the role of GLYI genes in plant regulation of abiotic stress resistance, especially resistance to heavy metals, and our previous finding that the cassava MeGLYI-13 gene confers Fe3+ resistance [33], we speculate that MeGLYI-13 may also confer other heavy metal stress resistance, which is worthy of further study. Here, in our study, we found that the expression of the MeGLYI-13 gene was induced by Zn2+ stress, and the peak value appeared at 6 h in the leaves and at 24 h in the roots (Figure S2A). The expression of the MeGLYI-13 gene was also induced by Cu2+ stress, and the peak value appeared at 12 h in the roots, stem and leaves (Figure S2B). The function of MeGLYI-13 in high-concentration Zn2+ and Cu2+ stress tolerance was further confirmed by gene overexpression in yeast cells and Arabidopsis seedlings (Figure 5). These results indicate that the MeGLYI-13 gene has a potential role in broader heavy metal stress resistance. However, the specific mechanism by which MeGLYI-13 enhances Zn2+ and Cu2+ resistance is still unclear. As evolutionarily ubiquitous metalloenzymes, GLYIs require bivalent metal ions to bind to their active sites to activate the enzyme [30]. The activity of different GLYI enzymes depends on different metal ions. For example, in rice, the OsGLYI-11.2 enzyme was activated by Ni2+ and Co2+ but not by Zn2+ [30], which is similar to two GLYIs (GLXI;1 and GLXI;2) in Arabidopsis [48]. In Arabidopsis, another GLYI (GLXI;3) could be activated by Mn2+, Co2+ and Zn2+ [48]. Whether Cu2+ plays a role in promoting GLYI activity is unknown. In our current study, we found that there were more than 10 binding sites for Zn2+ or Cu2+ in the MeGLYI-13 amino acid sequence, and there were six common binding sites for both (Figure 1). These consequences forecast that MeGLYI-13 has the binding potential of Zn2+ and Cu2+. Moreover, Zn2+ and Cu2+ both might have a role in MeGLYI-13 activity, but the degree to which they activate the MeGLYI-13 enzyme needs to be verified experimentally.
In this study, we tested the heavy metal tolerance function of the MeGLYI-13 gene in model organisms of yeast and Arabidopsis, because these two organisms are commonly used in laboratory studies and are easy to manipulate. Of course, if other organisms have matured genetic transformation systems, they can also be used to produce transgenic lines that overexpress MeGLYI-13 to repair soil contaminated by heavy metals. In addition, virus-induced gene silencing (VIGS) is a powerful tool to trigger transient sequence-specific gene silencing in plants, which is very helpful for studying gene function. VIGS techniques have also been reported in cassava but mostly for cassava leaves and less for the roots [49,50,51,52]. It would be helpful to silence the MeGLYI-13 gene in SC8 cassava roots by VIGS to better understand the function of MeGLYI-13 in the tolerance of plants to Zn2+, Cu2+ and other heavy metals. Whilst the ultimate proof of function needs to be provided in future experiments in which loss-of-function of the MeGLYI-13 gene is assessed directly in cassava, the data obtained in Arabidopsis suggest that it is highly likely that this gene mediates the tolerance of plants to Zn2+ and Cu2+ toxicity.
In our study, the MeGLYI-13 gene cloned from the SC8 cassava variety has one nucleotide difference and one amino acid difference compared to the sequence from the AM560-2 cassava variety (Figure S1). This result suggests that the MeGLYI-13 gene may display genotypic differences. Therefore, it is necessary to conduct sequencing verification in the study of other genotypes of cassava. Interestingly, in rapeseed, the BnGLYI-3 gene from a heat-tolerant rapeseed cultivar was found to improve the heat and cold tolerance of yeast cells, while the BnGLYI-2 gene from a heat-sensitive rapeseed cultivar was not, although there are only two amino acid differences between the BnGLYI-2 and BnGLYI-3 sequences [53]. It can be seen that, sometimes, changes in a few amino acids can lead to large differences in gene function. Therefore, whether the difference in the sequence of cassava MeGLYI-13 will affect the function of the gene remains to be verified.
Protein subcellular localization is very important for the study of protein function. Onion epidermal cells are large, living, transparent monolayer cells, which are suitable for visual fluorescence fusion protein observation. The research presented in our study shows that MeGLYI-13 from SC8 cassava is located in the nucleus of onion epidermal cells (Figure 2). To our knowledge, MeGLYI-13 is the first glyoxalase I protein found to be localized in the nucleus in cassava. The nuclear location of MeGLYI-13 is suggestive of a regulatory role. A few reports also found nuclear localization of GLYIs, such as AtGLYI2-4 in Arabidopsis and OsGLYI-8 in rice [25]. Their nuclear localization signal (NLS) sequences were predicted and experimentally verified [25]. Meanwhile, MG has also been detected in the nucleus, which indicated that the MG detoxification mechanism presents in the nucleus. According to previous studies, we also used the same tool to predict a strong nuclear localization signal in the MeGLYI-13 amino acid sequence (Figure S3). These results confirmed the nuclear location of MeGLYI-13, suggesting that MeGLYI-13 may play a vital role in nuclear MG detoxification, protecting the nucleus from adverse stress.
It seems strange that MeGLYI-13 is localized in the nucleus and resistant to heavy metal ion toxicity. However, previous studies have confirmed that heavy metal ions can enter the nucleus after crossing the plasma membrane barrier, thus interfering with nucleic acids, chromatin and the internal components [54,55,56]. For example, a higher concentration of copper induced a change in the nucleus structure, higher chromatin density and lower DNA synthesis [56]. DNA damage was found in potatoes after cadmium treatment [55]. In order to reduce the damage of heavy metals to the nucleus, some genes in the nucleus actively respond; in addition to some common transcription factors, there are some other genes. For instance, cysteine synthase 1 (OsCS1) in rice was proven to be localized in the nucleus, and its overexpression transgenic plant obtained a higher tolerance to cadmium stress [57]. A rice heavy metal-associated isoprenylated plant protein gene (OsHIPP17) was localized in the nucleus, and overexpressing it in Arabidopsis enhanced the copper tolerance of transgenic plants [58]. Therefore, we speculate that MeGLYI-13 localization in the nucleus may be related to the clearance of excessive MG caused by heavy metals in the nucleus, and the binding sites for Zn2+ or Cu2+ may be used for binding Zn2+ or Cu2+ in order to activate MeGLYI-13 enzyme activity.

4. Materials and Methods

4.1. MeGLYI-13 Gene Analysis

The total RNA from 40-day tissue culture seedlings of SC8 cassava (Manihot esculenta Crantz) was isolated following the instructions of the Plant Total RNA Isolation Kit Plus (FOREGENE). The RTIII Super Mix (with dsDNase) kit (MonScript, Monad, Suzhou, China) was used to clear away possible residual DNA fragments in the total RNA solution, and after, the first strand of cDNA was synthesized. According to the sequence of Manes.12G054100 in the AM560-2 cassava genome data, Primer Premier 6.0 was used to design a pair of specific primers with a predicted length of 556 bp (Table S1). The 50 μL PCR system included PrimeSTAR HS (premix) 25 μL, template cDNA (1 μg/μL) 1μL, 10 μmol/L upstream and downstream primers 2 μL each of ddH2O 20 μL. The PCR procedure included 94 °C pre-degeneration for 5 min, 98 °C degeneration for 10 s, 55 °C annealing for 15 s, 72 °C extension for 1 min for 32 cycles and 72 °C extension for 10 min. PCR production was then detected using 1% agarose gel electrophoresis. The target fragment was purified by the use of a Gel Extraction Kit (Omega) and then sequenced by Sangon Biotech. The sequence alignment was performed using Jellyfish software version 4.1.1.0. The Zn2+ and Cu2+ binding sites were predicted by the online tool MIB2 (http://bioinfo.cmu.edu.tw/MIB2/, accessed on 21 November 2022). The nuclear localization signal on the MeGLYI-13 protein sequence was predicted by using the online tool cNLS Mapper (http://nls-mapper.iab.keio.ac.jp/cgi-bin/NLS_Mapper_form.cgi, accessed on 6 January 2023).

4.2. Subcellular Localization of MeGLYI-13

SnapGene software version 4.2.4.0 was used to design homologous recombination primers. SpeI restriction sites were introduced upstream and KpnI restriction sites downstream of ORF of the MeGLYI-13 gene. By using the Seamless Assembly Cloning Kit (Clone Smarter, USA), the fusion expression vector pCAMBIA1300-35S-MeGLYI-13: GFP was obtained. The vector was then transferred into DH-5α E. Coli for sequencing. Both the pCAMBIA1300-35S-MeGLYI-13: GFP vector and the control pCAMBIA1300-35S-GFP vector were transformed into living onion epidermal cells by an Agrobacterium Tumefaciens LBA4404-mediated protocol. After incubation at 28 °C on solid MS medium (pH5.8) for 16–24 h, the onion cells were studied using a Laser Scanning Confocal Microscope FV1000 (Olympus, Tokyo, Japan).

4.3. Expression Analysis of MeGLYI-13 in Response to Zinc and Copper Stress

To explore the response of the MeGLYI-13 gene in zinc or copper stress, 30-day-old cassava seedlings of the SC8 variety were treated with 60 µmol/L Zn(CH3COOH)2 or 100 µmol/L CuSO4 for 0, 2, 6, 12 and 24 h, respectively. Three plants and three biological replicates were set for each treatment. The leaves, roots and stems were harvested and used to extract RNA and perform a qRT-PCR analysis according to the method in our previous report [33].

4.4. Yeast Transformation and Stress Treatments

In order to construct a yeast expression vector, a restriction site (KpnI) was introduced upstream of the MeGLYI-13 ORF, and an additional restriction site (EcoRI) was introduced downstream of the gene. Afterward, a fragment of MeGLYI-13 gene was inserted into the yeast expression vector pYES2 by the double restriction method to generate the pYES2-MeGLYI-13 yeast expression vector. The sequences of the primers pYES2-MeGLYI-13 forward (F) and pYES2-MeGLYI-13 reverse (R) are given in Supplementary Table S1. The empty plasmid pYES2 and recombinant plasmid pYES2-MeGLYI-13 were transformed into E. coli DH5α using the chemical transformation method, and the presence of the transgene was detected by PCR and subsequent sequencing. The pYES2 and pYES2-MeGLYI-13 plasmids in E. coli DH5α with correct sequencing were extracted and transferred into INVSc1 yeast strains using the PEG/LiAc transfer method. In order to investigate the response of MeGLYI-13 transgenic yeast to Zn2+ and Cu2+ ion stress, transgenic yeast and control yeast were cultured in SC-Ura liquid medium (2% galactose, 25 mg/L Amp) at 28 °C and 180 r/min for 36 h. The bacteria were then collected by centrifugation and then resuspended to OD600 = 0.2 by SC-Ura liquid medium containing different Zn2+ (0 mmol/L, 6 mmol/L, 12 mmol/L, 18 mmol/L, 24 mmol/L and 30 mmol/L) and Cu2+ (0 mmol/L, 0.6 mmol/L, 1.2 mmol/L, 1.8 mmol/L, 2.4 mmol/L and 3.0 mmol/L). Then, 2 mL of each bacterial solution was inserted into a new tube and incubated at 28 °C, 180 r/min for 30 h, and the OD600 value was measured. Three biological replicates were performed per treatment. At the same time, the survival difference between the recombinant yeast and the control yeast after 30 mmol/L Zn2+ and 3 mmol/L Cu2+ stress was assessed. Then, 5 mL of bacterial solution cultured for 36 h was centrifuged to collect the bacteria, respectively, and suspended in the treatment solution and cultured for 6 h at 28 °C, 180 r/min. After that, the treated bacterial suspension was attenuated into the six 10−1, 10−2, 10−3, 10−4, 10−5 and 10−6 gradient dilutions. Then, 5 μL of diluted bacterial suspension was plated on SC-Ura solid medium. The bacteria were cultured at 28 °C, and two days later, the survival difference of the yeast was observed and photographed.

4.5. Transformation of Arabidopsis with MeGLYI-13

The Agrobacterium LBA4404 containing the pCAMBIA1300-MeGLYI-13:GFP expression vector was shock cultured until the OD600 value of its suspension reached 0.9~1.0 at 28 °C and with 250 r/min shaking with YEP liquid medium containing 50 mg/L kanamycin and 50 mg/L rifampicin. The suspension was then centrifuged to collect the bacteria at 4000× g and 4 °C. After washing one time in 1/2 MS liquid medium (containing 0.03% Silwet-L77), the bacteria were resuspended to an OD600 value of 0.6~0.8 and then kept away from light for 3~12 h. Inflorescences containing unopened flowers of Arabidopsis were dipped into the above bacteria suspension for 3 min, then cultured in the dark for 16~24 h. After that, the Arabidopsis were cultured under 16 h light/8 h dark at 22 °C conditions for three days. The mature seeds were collected as the T0 generation. The seeds were disinfected with 75% ethanol (containing 0.01% Tween 20) for 5 min and washed by absolute ethanol twice. Next, the seeds were evenly spread on 1/2 MS solid medium (containing 80 mg/L hygromycin) and vernalized at 4 °C under darkness for two to four days and then cultured under 12 h light/12 h dark at 22 °C. Ten days later, the Arabidopsis seedlings that grew well on the resistance medium were transplanted into soil for pot culturing.

4.6. Detection of Transgenic Arabidopsis

To detect transgenic Arabidopsis, we performed DNA-level PCR. Then, 2 mm2 leaves of each Arabidopsis seedling were collected separately and treated with 20 μL lysate (M5 mix Kit, MF848, Mei5bio, Beijing, China) until the liquid turned green. The liquid was kept at 95 °C for 5 min and shortly centrifuged for 1 min at 10,000× g, and then, the supernatant was used as a template for the PCR. The sequences of primer pairs 1300F and 1300R used to detect the transgenic plants are given in Supplementary Table S1. The 20 μL PCR reaction liquid contained 10 μL 2 × M5 HiPer plus Taq HiFi, 0.5 μL 1300F, 0.5 μL 1300R, 2 μL template DNA and 7 μL ddH2O. The PCR procedure included 95 °C for 3 min, 94 °C for 25 s, 53 °C for 25 s, 72 °C for 45 s for 34 cycles and a 72 °C extension for 5 min. The PCR product was detected by 1% agarose electrophoresis.

4.7. Expression Level Detection of MeGLYI-13 in Transgenic Arabidopsis

To detect the expression level of MeGLYI-13 in PCR-positive lines by qRT-PCR, the total RNA from the PCR-positive Arabidopsis lines was extracted, purified and reverse-transcribed by successively using the RNA Isolation Kit (Foregene, Chengdu, China) and RTIII Super Mix with dsDNase kit (Monad, Suzhou, China). To perform the reaction for qRT-PCR, the reaction system was prepared according to the manual of SYBR® Premix Ex TaqTM II (Takara, Otsu, Japan). For detecting the MeGLYI-13 gene, the primer pairs of QMeGLYI-13F and QMeGLYI-13R were used. For reference, the Actin gene primer pairs ActinF and ActinR were used (see Supplementary Table S1). QRT-PCR was performed and recorded using a 7900 HT Fast Real-Time PCR System (ABI, Foster City, CA, USA). The relative expression levels of MeGLYI-13 in different Arabidopsis lines were calculated using the 2−∆∆CT method. The wild type Arabidopsis (WT) was used as the control. Three biological replicates were performed in this experiment.

4.8. Fluorescence Detection of Transgenic Arabidopsis

The T1 generation seeds were germinated on 1/2 MS medium after vernalization, as described in Section 2.4. Then, a dual-wavelength fluorescent protein excitation light source (Luyor-3415RG, Luyor Corporation, Shanghai, China) was used to irradiate the germinated seedlings. Positive transgenic Arabidopsis seedlings were selected, as they showed green fluorescence under light, while the wild type and transgenic lines with too-low gene expression showed a red color [59].

4.9. Znic and Copper Stress Treatment of Transgenic Arabidopsis

After identification by PCR, qRT-PCR and fluorescence, some transgenic lines were selected for further experiments. The WT seeds and T4 generation seeds of three selected transgenic Arabidopsis lines were sterilized and seeded on 1/2 MS medium at 4 °C for 2 days of vernalization and then germinated at 22 °C for 12 h light/12 h dark for 8 days of germination. Then, the Arabidopsis seedlings were transferred to 1/2 MS medium containing diverse concentrations of Zn(CH3COOH)2 (0, 4 and 60 μmol/L) and CuCl2 (0, 10 and 60 μmol/L) for stress treatment under 22 °C for 12 h light/12 h dark cultural conditions for 12 days. The normal 1/2 MS medium (Zn2+ 15 μmol/L and Cu2+ 0.05 μmol/L) was set as the control. After treatment, any phenotypic changes were observed and photographed. The length of the taproot was measured. For each treatment, three plates were set up, and three plants were used per line in each plate.

4.10. Statistical Analysis

All data analyses and graphic output were performed using GraphPad Prism 8.0 software. The one-way analysis of variance was performed to analyze the data, and Dunnett’s multiple comparisons test was used to assess the significant differences.

5. Conclusions

In summary, a MeGLYI-13 gene was cloned from SC8 cassava with a cDNA length of 495 bp and protein length of 164 amino acids. The MeGLYI-13 from SC8 cassava germplasm was different to that from AM560 cassava germplasm with one nucleobase difference in the cDNA sequence and one amino acid in the protein sequence. Subcellular localization of the MeGLYI-13 protein was demonstrated in the nucleus. There, six binding sites were predicted to bind to both Zn2+ and Cu2+ in the MeGLYI-13 protein sequence. The overexpression of MeGLYI-13 enhanced the Zn2+ and Cu2+ toxicity in both yeast cells and Arabidopsis seedlings. The results of this study highlight the function of one MeGLYI-13 gene from cassava in resistance to multiple heavy metals stress, possibly by stabilizing the MG balance in the nucleus. It can thus likely be used as a target for genetic modification to improve crop resistance to heavy metals. Therefore, we will continue to explore the mechanisms by which MeGLYI-13 regulates Zn and Cu toxicity tolerance in the future.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/plants12193375/s1: Figure S1: Sequence alignments of MeGLYI-13 CDS (A) and amino acids (B) between SC8 and AM560-2 cassava genome data. Figure S2: Expression profile of MeGLYI-13 in cassava seedlings after 60 μmol/L Zn2+ and 100 μmol/L Cu2+ treatments. Figure S3: Nuclear localization signal (NLS) prediction of MeGLYI-13 in cassava SC8. The red characters indicate the NLS sequence. Table S1: Sequences of the primers.

Author Contributions

Methodology, R.L. and F.T.; Investigation, R.L., F.T., Y.C., Q.Z. and Z.D.; Writing—original draft, R.L. and F.T.; Writing—review and editing, A.R.F., Y.Y., J.L., Y.W., X.H. and J.G. and Supervision, J.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Provincial Natural Fund of Hainan Province: 323RC537, the Earmarked Fund for China Agriculture Research System: CARS-11-HNGJC, the Hainan Yazhou Bay Seed Laboratory project: B21HJ0303, the National Natural Science Foundation of China: 31601359, the State Scholarship Fund of China Scholarship Council: CSC202003260020, the Major Science and Technology plan of Hainan Province: ZDKJ2021012 and the Central Public-Interest Scientific Institution Basal Research Fund for Chinese Academy of Tropical Agricultural Sciences: 1630052022008.

Data Availability Statement

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Acknowledgments

We are grateful for the patience and comments of the editors and the reviewers.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ruffing, A.M.; Anthony, S.M.; Strickland, L.M.; Lubkin, I.; Dietz, C.R. Identification of Metal Stresses in Arabidopsis thaliana Using Hyperspectral Reflectance Imaging. Front. Plant Sci. 2021, 12, 624656. [Google Scholar] [CrossRef]
  2. Alsafran, M.; Saleem, M.H.; Rizwan, M.; Al Jabri, H.; Usman, K.; Fahad, S. An overview of heavy metals toxicity in plants, tolerance mechanism, and alleviation through lysine-chelation with micro-nutrients—A novel approach. Plant Growth Regul. 2023, 100, 337–354. [Google Scholar] [CrossRef]
  3. Rehman, M.; Liu, L.J.; Wang, Q.; Saleem, M.H.; Bashir, S.; Ullah, S.; Peng, D.X. Copper environmental toxicology, recent advances, and future outlook: A review. Environ. Sci. Pollut. Res. 2019, 26, 18003–18016. [Google Scholar] [CrossRef]
  4. Behtash, F.; Abedini, F.; Ahmadi, H.; Mosavi, S.B.; Aghaee, A.; Morshedloo, M.R.; Lorenzo, J.M. Zinc Application Mitigates Copper Toxicity by Regulating Cu Uptake, Activity of Antioxidant Enzymes, and Improving Physiological Characteristics in Summer Squash. Antioxidants 2022, 11, 1688. [Google Scholar] [CrossRef]
  5. Saleem, M.H.; Fahad, S.; Adnan, M.; Ali, M.; Rana, M.S.; Kamran, M.; Ali, Q.; Hashem, I.A.; Bhantana, P.; Ali, M.; et al. Foliar application of gibberellic acid endorsed phytoextraction of copper and alleviates oxidative stress in jute (Corchorus capsularis L.) plant grown in highly copper-contaminated soil of China. Environ. Sci. Pollut. Res. 2020, 27, 37121–37133. [Google Scholar] [CrossRef]
  6. Liao, F.; Lilay, G.H.; Castro, P.H.; Azevedo, H.; Assuncao, A.G.L. Regulation of the Zinc Deficiency Response in the Legume Model Medicago truncatula. Front. Plant Sci. 2022, 13, 916168. [Google Scholar] [CrossRef]
  7. Stanton, C.; Sanders, D.; Kramer, U.; Podar, D. Zinc in plants: Integrating homeostasis and biofortification. Mol. Plant 2022, 15, 65–85. [Google Scholar] [CrossRef]
  8. Cabot, C.; Martos, S.; Llugany, M.; Gallego, B.; Tolra, R.; Poschenrieder, C. A Role for Zinc in Plant Defense against Pathogens and Herbivores. Front. Plant Sci. 2019, 10, 1171. [Google Scholar] [CrossRef]
  9. Paradisone, V.; Navarro-Leon, E.; Albacete, A.; Ruiz, J.M.M.; Esposito, S.; Blasco, B. Improvement of the physiological response of barley plants to both Zinc deficiency and toxicity by the application of calcium silicate. Plant Sci. 2022, 319, 111259. [Google Scholar] [CrossRef]
  10. Lu, X.; Liu, S.; Zhi, S.; Chen, J.; Ye, G. Comparative transcriptome profile analysis of rice varieties with different tolerance to zinc deficiency. Plant Biol. 2021, 23, 375–390. [Google Scholar] [CrossRef]
  11. Shakouri, M.; Mandoulakani, B.A.; Norouzi, M. Effect of Soil Zinc Deficiency on Morphological Traits, Yield, and Yield Components and Activity of Zinc-Containing Enzymes in Zn-Efficient and -Inefficient Bread Wheat Cultivars. Commun. Soil. Sci. Plant Anal. 2022, 53, 2529–2542. [Google Scholar] [CrossRef]
  12. Li, X.; Zhang, L.; Ren, H.; Wang, X.; Mi, F. Zinc toxicity response in Ceratoides arborescens and identification of CaMTP, a novel zinc transporter. Front. Plant Sci. 2022, 13, 976311. [Google Scholar] [CrossRef]
  13. Chen, G.; Li, J.; Han, H.; Du, R.; Wang, X. Physiological and Molecular Mechanisms of Plant Responses to Copper Stress. Int. J. Mol. Sci. 2022, 23, 12950. [Google Scholar] [CrossRef] [PubMed]
  14. Nazir, F.; Hussain, A.; Fariduddin, Q. Hydrogen peroxide modulate photosynthesis and antioxidant systems in tomato (Solanum lycopersicum L.) plants under copper stress. Chemosphere 2019, 230, 544–558. [Google Scholar] [CrossRef] [PubMed]
  15. Leng, X.; Mu, Q.; Wang, X.; Li, X.; Zhu, X.; Shangguan, L.; Fang, J. Transporters, chaperones, and P-type ATPases controlling grapevine copper homeostasis. Funct. Integr. Genom. 2015, 15, 673–684. [Google Scholar] [CrossRef]
  16. Hunter, C.; Ware, M.A.; Gleason, S.M.; Pilon-Smits, E.; Pilon, M. Recovery after deficiency: Systemic copper prioritization and partitioning in the leaves and stems of hybrid poplar. Tree Physiol. 2022, 42, 1776–1785. [Google Scholar] [CrossRef]
  17. Schulten, A.; Pietzenuk, B.; Quintana, J.; Scholle, M.; Feil, R.; Krause, M.; Romera-Branchat, M.; Wahl, V.; Severing, E.; Coupland, G.; et al. Energy status-promoted growth and development of Arabidopsis require copper deficiency response transcriptional regulator SPL7. Plant Cell 2022, 34, 3873–3898. [Google Scholar] [CrossRef]
  18. Zhang, T.; Wang, Y.; Ma, X.; Ouyang, Z.; Deng, L.; Shen, S.; Dong, X.; Du, N.; Dong, H.; Guo, Z.; et al. Melatonin Alleviates Copper Toxicity via Improving ROS Metabolism and Antioxidant Defense Response in Tomato Seedlings. Antioxidants 2022, 11, 758. [Google Scholar] [CrossRef]
  19. Raza Khan, A.; Fan, X.; Salam, A.; Azhar, W.; Ulhassan, Z.; Qi, J.; Liaquat, F.; Yang, S.; Gan, Y. Melatonin-mediated resistance to copper oxide nanoparticles-induced toxicity by regulating the photosynthetic apparatus, cellular damages and antioxidant defense system in maize seedlings. Environ. Pollut. 2022, 316 Pt 2, 120639. [Google Scholar] [CrossRef] [PubMed]
  20. Rehman, M.; Saleem, M.H.; Fahad, S.; Bashir, S.; Peng, D.X.; Deng, G.; Alamri, S.; Siddiqui, M.H.; Khan, S.M.; Shah, R.A.; et al. Effects of rice straw biochar and nitrogen fertilizer on ramie (Boehmeria nivea L.) morpho-physiological traits, copper uptake and post-harvest soil characteristics, grown in an aged-copper contaminated soil. J. Plant Nutr. 2021, 45, 11–24. [Google Scholar] [CrossRef]
  21. Yemets, A.; Horiunova, I.; Blume, Y. Cadmium, nickel, copper, and zinc influence on microfilament organization in Arabidopsis root cells. Cell Biol. Int. 2021, 45, 211–226. [Google Scholar] [CrossRef] [PubMed]
  22. Abdulmajeed, A.M.; Alnusairi, G.S.H.; Alharbi, M.H.; Almushhin, A.; Hasan, M.M.; Soliman, M.H. Alleviation of copper phytotoxicity by acetylsalicylic acid and nitric oxide application in mung bean involves the up-regulation of antioxidants, osmolytes and glyoxalase system. J. Plant Interact. 2021, 16, 201–212. [Google Scholar] [CrossRef]
  23. Zaid, A.; Mohammad, F.; Wani, S.H.; Siddique, K.M.H. Salicylic acid enhances nickel stress tolerance by up-regulating antioxidant defense and glyoxalase systems in mustard plants. Ecotoxicol. Environ. Saf. 2019, 180, 575–587. [Google Scholar] [CrossRef]
  24. Mostofa, M.G.; Hossain, M.A.; Fujita, M.; Tran, L.S. Physiological and biochemical mechanisms associated with trehalose-induced copper-stress tolerance in rice. Sci. Rep. 2015, 5, 11433. [Google Scholar] [CrossRef] [PubMed]
  25. Kaur, C.; Tripathi, A.K.; Nutan, K.K.; Sharma, S.; Ghosh, A.; Tripathi, J.K.; Pareek, A.; Singla-Pareek, S.L.; Sopory, S.K. A nuclear-localized rice glyoxalase I enzyme, OsGLYI-8, functions in the detoxification of methylglyoxal in the nucleus. Plant J. 2017, 89, 565–576. [Google Scholar] [CrossRef]
  26. Lin, F.; Xu, J.; Shi, J.; Li, H.; Li, B. Molecular cloning and characterization of a novel glyoxalase I gene TaGly I in wheat (Triticum aestivum L.). Mol. Biol. Rep. 2010, 37, 729–735. [Google Scholar] [CrossRef]
  27. Tan, Y.X.; Liang, J.; Zhang, Z.; Dong, S.Q.; Cao, J.F.; Qin, Q.P.; Zhang, Z.G.; Ni, D.A. Cloning and molecular characterisation of a putative glyoxalase I Gene (HfGlX I-1) of Daylily (Hemerocallis spp.). J. Plant Biochem. Biotechnol. 2022, 32, 265–273. [Google Scholar] [CrossRef]
  28. Li, Z.G. Methylglyoxal and Glyoxalase System in Plants: Old Players, New Concepts. Bot. Rev. 2016, 82, 183–203. [Google Scholar] [CrossRef]
  29. Kharbech, O.; Ben Massoud, M.; Sakouhi, L.; Djebali, W.; Mur, L.A.J.; Chaoui, A. Exogenous application of hydrogen sulfide reduces chromium toxicity in maize seedlings by suppressing NADPH oxidase activities and methylglyoxal accumulation. Plant Physiol. Biochem. 2020, 154, 646–656. [Google Scholar] [CrossRef]
  30. Mustafiz, A.; Ghosh, A.; Tripathi, A.K.; Kaur, C.; Ganguly, A.K.; Bhavesh, N.S.; Tripathi, J.K.; Pareek, A.; Sopory, S.K.; Singla-Pareek, S.L. A unique Ni2+-dependent and methylglyoxal-inducible rice glyoxalaseI possesses a single active site and functions in abiotic stress response. Plant J. 2014, 78, 951–963. [Google Scholar] [CrossRef]
  31. Singla-Pareek, S.L.; Yadav, S.K.; Pareek, A.; Reddy, M.K.; Sopory, S.K. Transgenic tobacco overexpressing glyoxalase pathway enzymes grow and set viable seeds in zinc-spiked soils. Plant Physiol. 2006, 140, 613–623. [Google Scholar] [CrossRef] [PubMed]
  32. Mohsin, S.M.; Hasanuzzaman, M.; Parvin, K.; Hossain, M.S.; Fujita, M. Protective role of tebuconazole and trifloxystrobin in wheat (Triticum aestivum L.) under cadmium stress via enhancement of antioxidant defense and glyoxalase systems. Physiol. Mol. Biol. Plants 2021, 27, 1043–1057. [Google Scholar] [CrossRef]
  33. Tang, F.; Li, R.; Zhou, Y.; Wang, S.; Zhou, Q.; Ding, Z.; Yao, Y.; Liu, J.; Wang, Y.; Hu, X.; et al. Genome-Wide Identification of Cassava Glyoxalase I Genes and the Potential Function of MeGLYⅠ-13 in Iron Toxicity Tolerance. Int. J. Mol. Sci. 2022, 23, 5212. [Google Scholar] [CrossRef] [PubMed]
  34. Kaya, C.; Ashraf, M.; Alyemeni, M.N.; Corpas, F.J.; Ahmad, P. Salicylic acid-induced nitric oxide enhances arsenic toxicity tolerance in maize plants by upregulating the ascorbate-glutathione cycle and glyoxalase system. J. Hazard. Mater. 2020, 399, 123020. [Google Scholar] [CrossRef] [PubMed]
  35. Altaf, M.M.; Diao, X.-P.; Altaf, M.A.; Ur Rehman, A.; Shakoor, A.; Khan, L.U.; Jan, B.L.; Ahmad, P. Silicon-mediated metabolic upregulation of ascorbate glutathione (AsA-GSH) and glyoxalase reduces the toxic effects of vanadium in rice. J. Hazard. Mater. 2022, 436, 129145. [Google Scholar] [CrossRef]
  36. Shen, S.L.; Chen, J.Q.; Chang, J.J.; Xia, B.C. Using bioenergy crop cassava (Manihot esculenta) for reclamation of heavily metal-contaminated land. Int. J. Phytoremediat. 2020, 22, 1313–1320. [Google Scholar] [CrossRef]
  37. Gasparatos, D. Soil Contamination by Heavy Metals and Metalloids. Environments 2022, 9, 32. [Google Scholar] [CrossRef]
  38. Chen, Y.; Chao, Z.-F.; Jin, M.; Wang, Y.-L.; Li, Y.; Wu, J.-C.; Xiao, Y.; Peng, Y.; Lv, Q.-Y.; Gui, S.; et al. A heavy metal transporter gene ZmHMA3a promises safe agricultural production on cadmium-polluted arable land. J. Genet. Genom. 2022, 50, 130–134. [Google Scholar] [CrossRef]
  39. Liu, Y.; Nour-Eldin, H.H.; Zhang, L.; Li, Z.; Fernie, A.R.; Ren, M. Biotechnological detoxification: An unchanging source-sink balance strategy for crop improvement. Trends Plant Sci. 2022, 28, 135–138. [Google Scholar] [CrossRef]
  40. Balparda, M.; Schmitz, J.; Duemmel, M.; Wuthenow, I.C.; Schmidt, M.; Alseekh, S.; Fernie, A.R.; Lercher, M.J.; Maurino, V.G. Two plant glyoxalase systems with different evolutionary origins detoxify diverse reactive carbonyl species. Plant Physiol. 2022, 191, 1214–1233. [Google Scholar] [CrossRef]
  41. Manoj, V.M.; Anunanthini, P.; Swathik, P.C.; Dharshini, S.; Narayan, J.A.; Manickavasagam, M.; Sathishkumar, R.; Suresha, G.S.; Hemaprabha, G.; Ram, B.; et al. Comparative analysis of glyoxalase pathway genes in Erianthus arundinaceus and commercial sugarcane hybrid under salinity and drought conditions. BMC Genom. 2019, 19, 986. [Google Scholar] [CrossRef] [PubMed]
  42. Yan, G.; Xiao, X.; Wang, N.; Zhang, F.; Gao, G.; Xu, K.; Chen, B.; Qiao, J.; Wu, X. Genome-wide analysis and expression profiles of glyoxalase gene families in Chinese cabbage (Brassica rapa L). PLoS ONE 2018, 13, e0191159. [Google Scholar] [CrossRef]
  43. Li, T.; Cheng, X.; Wang, Y.; Yin, X.; Li, Z.; Liu, R.; Liu, G.; Wang, Y.; Xu, Y. Genome-wide analysis of glyoxalase-like gene families in grape (Vitis vinifera L.) and their expression profiling in response to downy mildew infection. BMC Genom. 2019, 20, 362. [Google Scholar] [CrossRef] [PubMed]
  44. Jana, G.A.; Yaish, M.W. Functional characterization of the Glyoxalase-I (PdGLX1) gene family in date palm under abiotic stresses. Plant Signal Behav. 2020, 15, 1811527. [Google Scholar] [CrossRef] [PubMed]
  45. Wu, C.; Ma, C.; Pan, Y.; Gong, S.; Zhao, C.; Chen, S.; Li, H. Sugar beet M14 glyoxalase I gene can enhance plant tolerance to abiotic stresses. J. Plant Res. 2013, 126, 415–425. [Google Scholar] [CrossRef]
  46. Zeng, Z.; Xiong, F.; Yu, X.; Gong, X.; Luo, J.; Jiang, Y.; Kuang, H.; Gao, B.; Niu, X.; Liu, Y. Overexpression of a glyoxalase gene, OsGly I, improves abiotic stress tolerance and grain yield in rice (Oryza sativa L.). Plant Physiol. Biochem. 2016, 109, 62–71. [Google Scholar] [CrossRef]
  47. Wu, Q.; Gao, S.; Pan, Y.B.; Su, Y.; Grisham, M.P.; Guo, J.; Xu, L.; Que, Y. Heterologous expression of a Glyoxalase I gene from sugarcane confers tolerance to several environmental stresses in bacteria. PeerJ 2018, 6, e5873. [Google Scholar] [CrossRef]
  48. Schmitz, J.; Dittmar, I.C.; Brockmann, J.D.; Schmidt, M.; Hüdig, M.; Rossoni, A.W.; Maurino, V.G. Defense against Reactive Carbonyl Species Involves at Least Three Subcellular Compartments Where Individual Components of the System Respond to Cellular Sugar Status. Plant Cell 2017, 29, 3234–3254. [Google Scholar] [CrossRef]
  49. Lentz, E.M.; Kuon, J.E.; Alder, A.; Mangel, N.; Zainuddin, I.M.; McCallum, E.J.; Anjanappa, R.B.; Gruissem, W.; Vanderschuren, H. Cassava geminivirus agroclones for virus-induced gene silencing in cassava leaves and roots. Plant Methods 2018, 14, 73. [Google Scholar] [CrossRef]
  50. Luo, X.Q.; An, F.F.; Xue, J.J.; Zhu, W.L.; Wei, Z.W.; Ou, W.J.; Li, K.M.; Chen, S.B.; Cai, J. Integrative analysis of metabolome and transcriptome reveals the mechanism of color formation in cassava (Manihot esculenta Crantz) leaves. Front. Plant Sci. 2023, 14, 1181257. [Google Scholar] [CrossRef]
  51. Tuo, D.C.; Zhou, P.; Yan, P.; Cui, H.G.; Liu, Y.; Wang, H.; Yang, X.K.; Liao, W.B.; Sun, D.; Li, X.Y.; et al. A cassava common mosaic virus vector for virus-induced gene silencing in cassava. Plant Methods 2021, 17, 74. [Google Scholar] [CrossRef]
  52. Zeng, H.Q.; Xie, Y.W.; Liu, G.Y.; Wei, Y.X.; Hu, W.; Shi, H.T. Agrobacterium-Mediated Gene Transient Overexpression and Tobacco Rattle Virus (TRV)-Based Gene Silencing in Cassava. Int. J. Mol. Sci. 2019, 20, 3976. [Google Scholar] [CrossRef] [PubMed]
  53. Yan, G.; Lv, X.; Gao, G.; Li, F.; Li, J.; Qiao, J.; Xu, K.; Chen, B.; Wang, L.; Xiao, X.; et al. Identification and Characterization of a Glyoxalase I Gene in a Rapeseed Cultivar with Seed Thermotolerance. Front. Plant Sci. 2016, 7, 150. [Google Scholar] [CrossRef]
  54. Sas-Nowosielska, H.; Pawlas, N. Heavy metals in the cell nucleus—Role in pathogenesis. Acta Biochim. Pol. 2015, 62, 7–13. [Google Scholar] [CrossRef] [PubMed]
  55. Gichner, T.; Patkova, Z.; Szakova, J.; Znidar, I.; Mukherjee, A. DNA damage in potato plants induced by cadmium, ethyl methanesulphonate and gamma-rays. Environ. Exp. Bot. 2008, 62, 113–119. [Google Scholar] [CrossRef]
  56. Doncheva, S.; Nikolov, B.; Ogneva, V. Effect of copper excess on the morphology of the nucleus in maize root meristem cells. Physiol. Plant. 1996, 96, 118–122. [Google Scholar] [CrossRef]
  57. Shen, J.; Su, Y.; Zhou, C.; Zhang, F.; Zhou, H.; Liu, X.; Wu, D.L.; Yin, X.C.; Xie, Y.J.; Yuan, X.X. A putative rice L-cysteine desulfhydrase encodes a true L-cysteine synthase that regulates plant cadmium tolerance. Plant Growth Regul. 2019, 89, 217–226. [Google Scholar] [CrossRef]
  58. Shi, Y.; Jiang, N.; Wang, M.T.; Du, Z.Y.; Chen, J.; Huang, Y.Y.; Li, M.Y.; Jin, Y.F.; Li, J.H.; Wan, J.; et al. OsHIPP17 is involved in regulating the tolerance of rice to copper stress. Front. Plant Sci. 2023, 14, 1183445. [Google Scholar] [CrossRef]
  59. Lin, X.; Li, R.; Zhou, Y.; Tang, F.; Wang, Y.; Lu, X.; Wang, S.; Yao, Y.; Liu, J.; Hu, X.; et al. Overexpression of Cassava MeAnn2 Enhances the Salt and IAA Tolerance of Transgenic Arabidopsis. Plants 2021, 10, 941. [Google Scholar] [CrossRef]
Plants 12 03375 g001
Figure 1. The Zn2+ and Cu2+ ion binding site analysis in the MeGLYI-13 amino acid sequence. (A) The Zn2+ binding sites in MeGLYI-13. (B) The Cu2+ binding sites in MeGLYI-13.
Figure 1. The Zn2+ and Cu2+ ion binding site analysis in the MeGLYI-13 amino acid sequence. (A) The Zn2+ binding sites in MeGLYI-13. (B) The Cu2+ binding sites in MeGLYI-13.
Plants 12 03375 g001
Plants 12 03375 g002
Figure 2. Subcellular localization of MeGLYI-13 in onion epidermal cells.
Figure 2. Subcellular localization of MeGLYI-13 in onion epidermal cells.
Plants 12 03375 g002
Plants 12 03375 g003
Figure 3. Overexpression of MeGLYI-13 enhanced the yeast resistance to Zn2+ and Cu2+ toxicity. (A) Growth activity of yeast after different concentrations of Zn2+ treatments. (B) The growth of MeGLYI-13 transgenic yeast and empty yeast cells after the 30 mmol/L Zn2+ treatment. (C) Growth activity of yeast after different concentrations of Cu2+ treatments. (D) The growth of MeGLYI-13 transgenic yeast and empty yeast cells after the 3 mmol/L Cu2+ treatment. The error lines represent ±SD, n = 3; * signals a significant difference of p ≤ 0.05 and ** p ≤ 0.01.
Figure 3. Overexpression of MeGLYI-13 enhanced the yeast resistance to Zn2+ and Cu2+ toxicity. (A) Growth activity of yeast after different concentrations of Zn2+ treatments. (B) The growth of MeGLYI-13 transgenic yeast and empty yeast cells after the 30 mmol/L Zn2+ treatment. (C) Growth activity of yeast after different concentrations of Cu2+ treatments. (D) The growth of MeGLYI-13 transgenic yeast and empty yeast cells after the 3 mmol/L Cu2+ treatment. The error lines represent ±SD, n = 3; * signals a significant difference of p ≤ 0.05 and ** p ≤ 0.01.
Plants 12 03375 g003
Plants 12 03375 g004
Figure 4. Generation of MeGLYI-13 overexpressing Arabidopsis. (A) PCR analysis of transgenic Arabidopsis. M, DL2000 DNA Marker; 1~6, the transgenic Arabidopsis lines; +, the positive control (pCambia1300-MeGLYI-13-GFP vector plasmid); −, the negative control (wild type Arabidopsis). (B) Expression analysis of the different transgenic lines. WT, wild type Arabidopsis; OE1~OE6, the transgenic Arabidopsis lines. The error lines represent ±SD, n = 3; the character on the top of each bar represents a significant difference. (C) Fluorescence detection of the different transgenic lines. WT, wild type Arabidopsis; OE1~OE6, the transgenic Arabidopsis lines.
Figure 4. Generation of MeGLYI-13 overexpressing Arabidopsis. (A) PCR analysis of transgenic Arabidopsis. M, DL2000 DNA Marker; 1~6, the transgenic Arabidopsis lines; +, the positive control (pCambia1300-MeGLYI-13-GFP vector plasmid); −, the negative control (wild type Arabidopsis). (B) Expression analysis of the different transgenic lines. WT, wild type Arabidopsis; OE1~OE6, the transgenic Arabidopsis lines. The error lines represent ±SD, n = 3; the character on the top of each bar represents a significant difference. (C) Fluorescence detection of the different transgenic lines. WT, wild type Arabidopsis; OE1~OE6, the transgenic Arabidopsis lines.
Plants 12 03375 g004
Plants 12 03375 g005
Figure 5. Zinc and copper tolerance analysis of MeGLYI-13 transgenic Arabidopsis. (A) Phenotypes of Arabidopsis grown on normal 1/2 MS medium. (B) Statistics of the taproot length of Arabidopsis grown on normal 1/2 MS medium. (C) Phenotypes of Arabidopsis under different Zn2+ concentrations. (D) Statistics of the taproot length of Arabidopsis under different Zn2+ concentration treatments. (E) Phenotypes of Arabidopsis under different Cu2+ concentration stresses. (F) Statistics of the taproot length of Arabidopsis under different Cu2+ concentration treatments. The error lines represent ±SD, n = 3; ** indicates significant differences of p ≤ 0.01, respectively.
Figure 5. Zinc and copper tolerance analysis of MeGLYI-13 transgenic Arabidopsis. (A) Phenotypes of Arabidopsis grown on normal 1/2 MS medium. (B) Statistics of the taproot length of Arabidopsis grown on normal 1/2 MS medium. (C) Phenotypes of Arabidopsis under different Zn2+ concentrations. (D) Statistics of the taproot length of Arabidopsis under different Zn2+ concentration treatments. (E) Phenotypes of Arabidopsis under different Cu2+ concentration stresses. (F) Statistics of the taproot length of Arabidopsis under different Cu2+ concentration treatments. The error lines represent ±SD, n = 3; ** indicates significant differences of p ≤ 0.01, respectively.
Plants 12 03375 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, R.; Tang, F.; Che, Y.; Fernie, A.R.; Zhou, Q.; Ding, Z.; Yao, Y.; Liu, J.; Wang, Y.; Hu, X.; et al. MeGLYI-13, a Glyoxalase I Gene in Cassava, Enhances the Tolerance of Yeast and Arabidopsis to Zinc and Copper Stresses. Plants 2023, 12, 3375. https://doi.org/10.3390/plants12193375

AMA Style

Li R, Tang F, Che Y, Fernie AR, Zhou Q, Ding Z, Yao Y, Liu J, Wang Y, Hu X, et al. MeGLYI-13, a Glyoxalase I Gene in Cassava, Enhances the Tolerance of Yeast and Arabidopsis to Zinc and Copper Stresses. Plants. 2023; 12(19):3375. https://doi.org/10.3390/plants12193375

Chicago/Turabian Style

Li, Ruimei, Fenlian Tang, Yannian Che, Alisdair R. Fernie, Qin Zhou, Zhongping Ding, Yuan Yao, Jiao Liu, Yajie Wang, Xinwen Hu, and et al. 2023. "MeGLYI-13, a Glyoxalase I Gene in Cassava, Enhances the Tolerance of Yeast and Arabidopsis to Zinc and Copper Stresses" Plants 12, no. 19: 3375. https://doi.org/10.3390/plants12193375

APA Style

Li, R., Tang, F., Che, Y., Fernie, A. R., Zhou, Q., Ding, Z., Yao, Y., Liu, J., Wang, Y., Hu, X., & Guo, J. (2023). MeGLYI-13, a Glyoxalase I Gene in Cassava, Enhances the Tolerance of Yeast and Arabidopsis to Zinc and Copper Stresses. Plants, 12(19), 3375. https://doi.org/10.3390/plants12193375

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop