Mingying Liu1,2
Wenming Qiu1,2Xuelian He1,2,3Liu Zheng1,2,3Xixi Song1,2Xiaojiao Han1,2Jing Jiang1,2Guirong Qiao1,2Jian Sang1,2Mingqing Liu1,4
Renying Zhuo1,2*- 1State Key Laboratory of Tree Genetics and Breeding, Chinese Academy of Forestry, Beijing, China
- 2Key Laboratory of Tree Breeding of Zhejiang Province, The Research Institute of Subtropical Forestry, Chinese Academy of Forestry, Hangzhou, China
- 3Biotechnology Research Center of China Three Gorges University, Yichang, China
- 4Vocational Secondary Specialized School of Hedong District, Linyi, China
Cadmium is a major toxic heavy-metal pollutant considering their bioaccumulation potential and persistence in the environment. The hyperaccumulating ecotype of Sedum alfredii Hance is a Zn/Cd co-hyperaccumulator inhabiting in a region of China with soils rich in Pb/Zn. Investigations into the underlying molecular regulatory mechanisms of Cd tolerance are of substantial interest. Here, library screening for genes related to cadmium tolerance identified a gene resembling the rubber elongation factor gene designated as SaREFl. The heterologous expression of SaREFl rescued the growth of a transformed Cd-sensitive strain (ycf1). Furthermore, SaREFl-expressing Arabidopsis plants were more tolerant to cadmium stress compared with wild type by measuring parameters of root length, fresh weight and physiological indexes. When under four different heavy metal treatments, we found that SaREFl responded most strongly to Cd and the root was the plant organ most sensitive to this heavy metal. Yeast two-hybrid screening of SaREFl as a bait led to the identification of five possible interacting targets in Sedum alfredii Hance. Among them, a gene annotated as prenylated Rab acceptor 1 (PRA1) domain protein was detected with a high frequency. Moreover, subcellular localization of SaREF1-GFP fusion protein revealed some patchy spots in cytosol suggesting potential association with organelles for its cellular functions. Our findings would further enrich the connotation of REF-like genes and provide theoretical assistance for the application in breeding heavy metal-tolerant plants.
Introduction
Nowadays, as we are enjoying the conveniences brought by industrial developments, we are also being challenged with the threats of environmental contaminations caused by hazardous environmental pollutants (Shahid et al., 2014). Among these pollutants, the heavy metals, a loosely defined group of elements exhibiting metallic properties (Rascio and Navari-Izzo, 2011), are exerting persistent potential harmful effects on the environment and living beings as they might be bioaccumulated in food webs because of difficulties in detoxification and long residence times in soils (Ahmad and Ashraf, 2011). Those elements that resemble essential minerals, such as mercury (Hg), lead (Pb), and cadmium (Cd), are deemed to be more threatening as they are more probable to enter the cell through the existing mineral uptake machinery (Lin and Aarts, 2012). As a Class B metal, cadmium represents a major toxic heavy metal and cadmium pollution has aroused worldwide concerns because excessive Cd exposure would cause emphysema and osteoporosis, leading to irreversible damages to lungs, kidneys, and bones in humans (Straif et al., 2009). When plants encounter with an excess level of a heavy metal, they will display several biological and physiological changes including reduced biomass, leaf chlorosis, inhibited root growth, morphological alterations, and more seriously plant death (Yadav, 2010; Su et al., 2014).
Based on the previous studies, the mechanisms utilized by plants to resist the toxic effect of a heavy metal can be summarized into two major strategies (Verbruggen et al., 2009a; Anjum et al., 2012). The first one is the excluder strategy, in which the plants attempt to prevent heavy metals from entering the roots, for instance by limiting metal bioavailability or by regulating expression of proteins involved in assimilating and transporting metals (Lin and Aarts, 2012). The other strategy is the tolerance tactics, which rely on the coordination and maintenance of homeostasis, for instance by controlling of metal ions through compartmentalization and detoxification and elevating expression of stress-related genes (Zhang et al., 2001; Lin and Aarts, 2012). These ecophysiological adaptations to metalliferous environment produce heavy-metal hypertolerant plants which can withstand exposure to high metal and its uptake to the shoot by restricting their accumulation to the root. Heavy metal hyperaccumulators that accumulate metals preferentially transport them to the above-ground parts (Jaffré et al., 1976; Pollard et al., 2002).
Until now, more than 450 plant species have been characterized from a number of different families such as the Asteraceace, Brassicaceae, Caryophyllaceae, Poaceae, Violaceae, and Fabaceae and most of them are nickel (Ni) tolerant (Verbruggen et al., 2009a; Krämer, 2010). Among these hyperaccumulators, only three plant species are classified to be Cd hyperaccumulators: Arabidopsis halleri (Bert et al., 2002), Noccaea caerulescens (Brooks et al., 1977; Baker et al., 1994) (formerly Thlaspi caerulescens) and Sedum alfredii Hance (Yang et al., 2004). At the forefront of studies concerning the hyperaccumulation process T. caerulescens and N. caerulescens have been used (Milner and Kochian, 2008). Due to the relatively low biomass and slow growth, they primarily served as model systems to investigate and identify the underlying molecular and physiological mechanisms of heavy metal hyperaccumulation. Ultimately the research aims is to apply the knowledge gained about these mechanisms to improve higher biomass of crop plants (Verbruggen et al., 2009a; Krämer, 2010). Their potential for use in decontamination of large-scale Cd-contaminated soils is not practically applicable. Moreover, in recent years, numerous studies have been focused on determining the feasibility of applying the mechanisms that metal-hyperaccumulator plant species employ in the breeding of plants suited for the remediation or phytoremediation of metal-contaminated soils. Therefore, it is essential to investigate hyperaccumulation mechanisms of other Cd hyperaccumulators, particularly non-Brassica species with potential practical applications.
The hyperaccumulating ecotype (HE) of Sedum alfredii Hance is a Zn/Cd co-hyperaccumulator originally found in a region of China with soils rich in PB/Zn, possessing remarkable traits of accumulating up to 9,000 μg g-1 Cd and 29,000 μg g-1 Zn in its shoots without any toxicity symptoms (Yang et al., 2002, 2004). It is a species native to China and the only one metal hyperaccumulator plant from the Crassulaceae family. The potential of S. alfredii Hance for practical decontamination is related to its large biomass, rapid growth, and asexual propagation. Also it could serve as an ideal candidate for obtaining a functional understanding of the intra-family hyperaccumulation mechanisms (Sun et al., 2009). As a species not been studied extensively and profoundly, the underlying mechanisms for hyper-tolerance in S. alfredii Hance (HE) is poorly elucidated, due in part to the lack of genes identified to be involved in heavy metal resistance or accumulation.
Based on our previous studies, we screened for functional genes related to Cd tolerance through expressing cDNA library in yeasts and obtained a gene resembling the rubber elongation factor (REF) gene. REF, previously identified in rubber biosynthesis from Hevea brasiliensis and termed by Dennis and Light (1989), is an abundant protein associated with the lipid membrane surrounding the rubber particles. Previous studies have found that this protein is an essential element required for rubber molecule elongation which interacts with prenyl transferase to alter the stereochemistry of isopentenyl pyrophosphate (IPP) addition from the normal trans addition to cis and overrides the normal termination after two trans additions to affect the formation of cis-polyisoprene (Dennis et al., 1989). Although proteins belonging to a larger stress-related protein (SRP) family from plants, and several SRPP non-latex producing homologous genes and genes containing REF domains have been reported to be involved in stress response (Qi et al., 2009), direct evidences confirming REF or SRPP as SRPs were still deficient (Berthelot et al., 2014). Most findings until present are embracing the context of their roles in rubber biosynthesis while functions in other processes such as stress response are seldom reported (Kush et al., 1990; Han et al., 2000; Aoki et al., 2014).
In our previous study, a gene containing the REF domain in S. alfredii Hance (HE) was detected to improve yeast tolerance to cadmium. Herein, we performed further studies to validate and elucidate functions of this gene in cadmium resistance in S. alfredii Hance (HE). The enhanced expression of this gene may contribute to Cd hyperaccumulation and hypertolerance in S. alfredii Hance.
Materials and Methods
Plants Materials and Growth Condition
Hyperaccumulating ecotype of S. alfredii Hance (HE), first identified by Yang et al. (2002), was found to be flourishing in the deserted Pb/Zn mining district located in Quzhou city, Zhejiang Province, PR China. S. alfredii Hance (HE) seedlings were cultured in buckets filled with tap water and exposed to 25°C and long days (16 h light/8 h dark each day) in an artificial climate growth chamber. The S. alfredii seedlings for the stress treatment were vegetatively propagated to ensure homogeneity and grown in half-strength Hoagland–Arnon solution (Hoagland and Arnon, 1950) for about 2 weeks until the relatively vigorous roots grew out.
Arabidopsis thaliana (ecotype Col-0) plants were grown on potted soil in a climate-controlled greenhouse under controlled environments (16 h light/8 h dark cycles, white light intensity approximately 125 μmol m-2 s-1 at 22°C and 50% relative humidity) till plant completed its life cycle. In experiments that included homozygous transgenic plants (T3), seeds were surface sterilized and germinated on 1/2 MS (Murashige and Skoog, 1962) agar plates containing 25 mg L-1 hygromycin. Three-week-old wild-type, homozygous transgenic seedlings were used in Cd resistance and other experiments.
Cloning and Sequence Analysis of SaREFl
In order to identify S. alfredii Hance (HE) genes that confer Cd tolerance, we constructed a cDNA library of S. alfredii Hance (HE) with a SMART cDNA Library Construction Kit (CLONTECH, Palo Alto, CA, USA) according to the recommended protocol and transformed the library into the INVSc1 yeast strain. This strain does not grow on half-strength synthetic galactose agar plates lacking uracil and containing 80 μM CdCl2. We selected transformed colonies that grew in the presence of 80 μM CdCl2 and then isolated and sequenced the respective plasmids.
Among these colonies, there was a sequence of interest encoding a 248 amino acid polypeptide with homology to REF which had rarely been reported to be related to heavy-metal stress. The sequence was analyzed for coding probability with the DNATools program (Rehm, 2001). Comparison against the GenBank protein database was performed using the BLAST network server at the National Center for Biotechnology Information (Altschul et al., 1997), based on which the sequence was designated as S. alfredii Hance (HE) REF like gene (SaREFl). Multiple protein sequences representing the homologous genes from other species with the Genbank accession numbers were listed in Table 1. Homology of SaREFl with the REF family proteins was aligned using the MegAlign program by the CLUSTALW method in the DNASTAR software package (Burland, 2000). A phylogenetic tree was constructed using the neighbor-joining method within the package PHYLIP 3.5c software package (Felsenstein, 1993) using 1000 bootstrap replicates.
TABLE 1. Names and accession numbers of the REF protein family.
Moreover, the genomic sequence of SaREFl was amplified through PCR using genomic DNA as templates with primers SaREFO-F and SaREFO-R listed in Table 2. In order to compare the genomic structure, we searched the genomic sequences of Populus trichocarpa REF (PtREF, POPTR_0013s01800g), Seteria italic REF (SiREF, NW_004675963.1), Arabidopsis thaliana REF (AtREF, AT3G05500F22F7.5) and Ricinus communis REF (RcREF, RCOM_1433280) from GenBank and analyzed the intron-exon structure.
TABLE 2. Primers designed for the experiments in the study.
Assays for Effects of Different Heavy-Metal Ions on Expression of SaREFl
In order to validate the expression profiles of SaREFl under heavy-metal stress, S. alfredii Hance (HE) seedlings were subjected to different heavy-metal ions treatment and the expression profile of SaREFl was examined by means of qRT-PCR. The S. alfredii Hance (HE) seedlings were obtained by asexual propagation. We cut the shoot tops of S. alfredii Hance (HE), grew them further in a greenhouse for 2 months, and then vigorous and uniform seedlings were selected for use. The selected seedlings were transferred to half-strength Hoagland–Arnon solution (pH 6.0) which was aerated continuously and renewed every 2 days. After 2 weeks of growth, the seedlings with relatively flourishing roots were randomly divided into experimental and control groups in a study on the effects of heavy-metal stress. The concentrations of the heavy metals were set according to a previous study characterizing the expression of a reliable internal gene for S. alfredii Hance (Sang et al., 2013). Briefly, the experimental groups were treated with 400 μM CdCl2, 400 μM Pb(NO3)2, 500 μM CuSO4 and 500 μM ZnSO4, respectively, while the control ones were cultured in water.
Three different plant organs (roots, stems, and leaves) were collected at different time points including 0, 0.5, 6, 12, 24, 48, and 72 h. The respective organs of the seedlings in the control group were also sampled. All experiments were performed in triplicates. All plant materials collected were promptly frozen in liquid nitrogen and kept at -80°C until total RNA extraction.
Total RNA was prepared according to the manual of the Total RNA Purification Kit (Norgan Biotek Corp., Thorold, ON, Canada) and then treated with RNase-free DNase I (NEB BioLabs, Ipswich, MA, USA) to remove potential genomic DNA contamination. The quality of total RNA was first assessed by agarose gel electrophoresis and then further quantified by a NanoDrop2000 spectrophotometer (Thermo Scientific, Wilmington, DE, USA). The first strand cDNA was generated from 3 μg of RNA sample with the SuperScript III First-Strand Synthesis System (Invitrogen, Carlsbad, CA, USA) using oligo d(T) primer followed by the RNase H step (Invitrogen, Carlsbad, CA, USA). Following quantification of the cDNA templates and examination of primers (SaREF-F/SaREF-R, Table 2), qRT-PCR was performed on ABI 7300 real-time PCR system (Applied Biosystems, USA) as described by Liu et al. (2014). The housekeeping TUB gene was used as the reference for internal standardization validated by Sang et al. (2013) and the detailed primer sequences (TUB-F/TUB-R) were listed in Table 2.
Cloning and Construction of Expression Vectors
To assess Cd tolerance in yeast cells, purified full-length SaREFl products without stop codon were cloned into the pYES2.1/V5-His-TOPO vector (Invitrogen, Carlsbad, CA, USA) and positive clones were sequenced by general primers GAL1 Forward and V5 C-term Reverse sequencing primers provided by the kit (Table 2).
For the gene overexpressing and subcellular location in Arabidopsis, the cDNA of SaREFl was also amplified by PCR using High Fidelity KOD-Plus DNA Polymerase (Toyobo, Japan) with primers pENTR-F and pENTR-R listed in Table 2. The purified PCR products were then cloned into the Gateway entry vector pENTR and positive clones were further sequenced to verify the direction and sequence accuracy. The verified plasmid was then recombined in pH2GW7 and pK7WGF2.0 to generate pH2GW7-SaREFl for overexpression and pK7WGF2.0-SaREFl for determination of their subcellular locations, respectively.
All the constructs above were verified by sequencing prior to use to confirm error- free amplification and cloning. Primers used during the construction and their uses for preparation of specific constructs were listed in Table 2.
Heterologous Expression of SaREFl and Cd Tolerance test in Saccharomyces cerevisiae
To assess Cd tolerance in yeast cells, the correct pYES2.1/V5-His-TOPO construct was transformed into Cd-sensitive mutant yeast strain DTY167 (MATa ura3-52 his6 leu2-3,-112 his3-D200 trp1-901 lys2-801 suc2-D, ycf1::hisG) by the lithium acetate method (Szczypka et al., 1994). Transformed yeast cells were selected on half-strength synthetic dextrose agar plates lacking uracil.
For assessment of growth inhibition in response to Cd, SaREFl expressing yeast lines were grown up to 2 OD600 nm units followed by serial dilution and spotting on half-strength synthetic galactose-uracil (SG-ura) agar plates in the absence or presence of 30 μM CdCl2 at 30°C for 3 days before photographs were taken. To examine the relative growth inhibition, yeast cultures on logarithmic phase were firstly diluted to OD600nm of approximately 0.4 unit with half-strength SG-ura supplemented with 30 μM CdCl2, and then the yeast cultures were kept on growing at 30°C for 72 h. Relative growth was determined through measuring OD600 nm every 12 h. Empty vector (EV) pYES2.1/V5-His-TOPO transformed yeast cells were considered as control in each experiment.
To determine whether expression of SaREFl conferred tolerance to other heavy metals, the constructed pYES2.1/V5-His-TOPO vector containing SaREFl was introduced into yeast strain INVSc1 using a lithium acetate method (Szczypka et al., 1994). Two lines of transformed yeast cells containing the vectors pYES2.1/V5-His-TOPO (control) or pYES2.1-SaREFl were cultivated on medium supplemented with 5 mM Cu2+, 1.5 mM Pb2+ and 15 mM Zn2+ at 30°C for 3 days before being photographed. The concentrations of the above three heavy metals were set according to a previous study characterizing the functions of a heavy metal-associated protein (AcHMA1) from the Halophyte, Atriplex canescens (Pursh) Nutt. (Sun et al., 2014).
Overexpression of SaREFl in Arabidopsis and Cadmium Stress Treatment
Arabidopsis thaliana Columbia was transformed with the plasmid pH2GW7-SaREFl with the Agrobacterium-mediated floral dip method as described previously (Hsiao et al., 2007). Several independent lines (T1) were selected on 1/2 MS agar plates containing 25 mg L-1 hygromycin and transgenic lines were further validated through PCR. The PCR validation was based on genomic DNA template primed with general primers of antibiotic (Hyg-F/Hyg-R) and specific gene primers (SaREFO-F/SaREFO-R) listed in Table 2. The expression of SaREFl in transgenic A. thaliana was asserted by RT-PCR using specific gene primers (SaREFO-F/SaREFO-R). Total RNA was extracted from young leaves of transgenic and WT plants and first strand cDNA was synthesized as described above. RT-PCR was carried out with parameters set as below: denaturation at 94°C for 5 min, 30 cycles of amplifications (94°C for 30 s, 55°C for 30 s, and 72°C for 60 s), followed by a final cycle of 72°C for 7 min. Three T3 homozygous lines (OE-1, OE-2, and OE-3) with single copy insertions were used to analyze the growth characteristics under normal and heavy-metal stress conditions in the following tolerance study.
For root length measurements, seeds of WT and three transgenic A. thaliana homozygous lines expressing SaREFl were germinated and grown vertically on half-strength MS (½ MS) media in the presence or absence of 100 μM cadmium for 10 d after stratification (4°C for 48 h in dark) (Thomine et al., 2003). After the indicated days of growth, A. thaliana seedlings were spread on a filter paper and individual root length was determined using a vernier caliper. The fresh weight (FW) was subsequently weighed and recorded. The experiments were replicated three times. The variables were expressed as mean standard deviation (SD) of data collected for 15–20 independent plants for each line. Statistical differences between WT and three transgenic lines were calculated by student t-test (P < 0.05).
For adult plant experiments, four-leaf stage wild-type and SaREFl-expressing Arabidopsis seedlings grown on MS agar were replanted into potted soils and grown under a 16/8-h photoperiod at 22°C for 3 weeks. Then the plants were watered with 500 μM CdCl2 and sampled 10 days later. Healthy young leaves of stress-treated and control plants were taken for further physio-biochemical analyses. Lipid peroxidation was determined based on the concentration of malondialdehyde (MDA) produced by thiobarbituric acid (TBA) reaction (Draper and Hadley, 1990). In brief, leaf samples were extracted in 2 ml 0.1% trichloroacetic acid (TCA) and 0.5 ml leaf extract was reacted with 2.0 ml of TBA reagent followed by boiling at 95°C for 30 min. The reaction mixtures were cooled on ice, centrifuged at 10,000 g for 5 min and absorbance of the supernatants was measured at 440, 532, and 600 nm.
To determine electrolyte leakage (EL), leaves were firstly washed thoroughly with deionized water to remove surface-adhered electrolytes, and then they were kept in closed vials containing 10 ml of deionised water and incubated at 25°C on a rotary shaker for 24 h. Subsequently, the electrical conductivity (EC) of the solution (Lt) was determined using a conductivity meter (SevenEasy, Mettler Toledo AG 8603, Switzerland). Leaf extracts were autoclaved at 120°C for 20 min, followed by cooling to 25°C and then EC (L0) was determined. The EL was determined: EL (%) = (Lt/L0) X 100.
Cadmium content measurements were also performed. Stress-treated plants were washed in distilled water three times and desorbed for 10 min with 2 mM CaSO4 and 10 mM ethylenediaminetetraacetic acid to remove surface and apoplastic heavy metals, then the shoots and leaves were excised from the plants, followed by the drying at 105°C for 30 min and 80°C for 48 h (Cailliatte et al., 2010). After weighing, the dry tissues were ground into powder and then digested in HNO3/HClO4 (4:1, v/v). The acid digests were diluted with deionized water and Cd2+ concentration was determined using atomic absorption spectrophotometry (Perkin-Elmer model 3300). The experiment was replicated three times, with each replicate consisting of at least 12 plants.
Yeast Two-Hybrid Experiments
In order to identify potential genes interacting with SaREFl, the yeast two-hybrid experiments using yeast mating were carried out following the instructions supplied by Clontech Matchmaker Gold Yeast Two-Hybrid System. Firstly, for the preparation of cDNA library, total RNA was extracted from fresh tissues including root, stem, and leaf using an RNeasy Plant Minikit (Qiagen) and treated with DNase I (Ambion). Matchmaker Library Construction and Screening kit (Clontech, Mountain View, CA, USA) was applied to generate an activation domain (AD)-tagged cDNA library, which was subsequently cloned into the plasmid pGADT7-sfiI (~2.5 × 106 clones). The mixed plasmids of positive inserts were further amplified in yeast strain Y187. A GAL4 DNA-BD fusion was generated by cloning SaREFl in frame with the GAL4 DNA binding domain of the bait plasmid pGBKT7 and transformed into cells of the yeast strain Y2HGold. It is imperative to confirm that the bait SaREFl does not autonomously activate the reporter genes in Y2HGold, in the absence of a prey protein. The absence of self-activation was verified by transformation of the bait alone, selected on minimal medium lacking tryptophan and histidine.
For library screening, the Clontech Matchmaker Gold Yeast Two-Hybrid System was applied to perform the yeast two-hybrid experiments using yeast mating. A concentrated overnight culture of the bait strain (Y2HGold [pGBKT7 + SaREFl]) was combined with the library strain and incubated at 30°C for approximately 20 h with gentle shaking until zygotes could be observed under a phase contrast microscope (40X). The positive colonies that activated at least two reporter genes (His 3, Ade2 or LacZ) were plated onto new plates and further identified by sequencing. In order to confirm that the positive interactions are genuine, the library plasmid responsible for activation of reporters was rescued and then co-transformation of pGBKT7/SaREFl and positive AD/library plasmids into Y2HGold was conducted. Those independent colonies that could activate two reporter genes (His 3, Ade2) were preliminarily verified as being genuine and the prey insert could be identified by sequencing. Y2HGold [pGBKT7-53] mating with Y187 [pGADT7-T] was applied as positive control. Negative control was performed using pGBKT7-Lam and pGADT7-T. All control constructs were provided by the manufacturer (Clontech).
Transient Expression and Subcellular Localization of SaREFl-GFP
The verified construct of pK7WGF2.0-SaREFl for subcellular localization study was used to transform onion epidermal cell layers using Helium Biolistic gene transformation system (Du Pont). Inner epidermal layers were peeled and placed inside up on Petri dishes containing MS basal medium containing 2.5 mg L-1 amphotericin B (Sigma, St. Louis, MO, USA) and 6% (w/v) agar. DNA samples containing pK7WGF2.0-SaREFl were prepared as described by the manufacturer (Du Pont) which were coated with 1 μM gold particles and followed by the introduction into onion epidermal cells using particle bombardment (PDS-100/He particle delivery system, Bio-rad). Then, Petri dishes were sealed with parafilm and incubated for 18–24 h at 28°C in the dark. After 20 h, fluorescent cells were imagined. The green fluorescence was observed through a laser scanning confocal microscope (LSM510; Karl Zeiss, Jena, Germany) excitated at 488 nm for GFP.
Statistical Analysis
Each experiment was performed under a completely randomized design with independent experiments with at least three replications. Level of significance was calculated using student t-test (P < 0.05).
Results
Sequence and Phylogenetic Analysis of SaREFl
To identify S. alfredii Hance (HE) genes that confer Cd tolerance, we transformed the yeast strain INVSc1 with a S. alfredii cDNA library. We selected colonies that grew in the presence of 80 μM CdCl2 and sequenced the respective plasmids. Among these colonies, a transcript of interest displayed a high homology to REF. The BLASTp searching at NCBI showed that the protein encoded by the cDNA that has a REF domain at residues 11–239. Therefore, the cDNA encoded an S. alfredii Hance (HE) REF-like gene, and we designated it as SaREFl. REF was a protein named in the context of their roles in rubber biosynthesis and had received the most attention because they were the most abundant rubber membrane particles. Interestingly, it had rarely been reported to be related to heavy-metal stress. Multiple analysis tools were applied to uncover the preliminary bioinformatics features of SaREFl.
The full-length cDNA of SaREFl is 968 bp, and its longest open reading frame codes for a protein of 248 amino acids with a predicted molecular mass of approximately 26.97 kDa. The 5′-untranslated region (UTR) is 78 bp, and the 3′-UTR is 137 bp with a polyadenylation tail (Figure 1A). The initiation codon (ATG) was assigned on the basis that there is no ATG in the 5′-UTR 78 nucleotides and the DNA surrounding the initiation codon ATG has a purine at positions of both -3 and +4, which is in accordance with the Kozak consensus sequence (A/GXXATGG) and is suggested to optimize translational efficiency (Kozak, 1987). SignalP software analysis by the Signal IP 3.0 server (Nielsen et al., 1997) revealed that SaREFl had no putative N-terminal signal peptide. The SaREFl sequence was submitted to the GenBank database and assigned with an accession number KP771801.
FIGURE 1. Bioinformatics analysis of SaREFl. (A) Nucleotide and deduced amino acid sequences of Sedum alfredii REFl cDNA (accession number in GenBank: KP771801). The in-frame stop codon in 5′UTR is marked by Δ. The start codon is boxed. The asterisk represents the stop codon. The polyadenylation signal is underlined by a single line. (B) Phylogenetic trees of SaREFl with other representative members of REF protein family. The phylogenetic trees of SaREFl as well as their counterparts of other REFs and REF protein family members constructed by the neighbor-joining method within the package PHYLIP 3.5c. Bootstrap majority consensus values on 1000 replicates are indicated at each branch point in percent. Accession numbers for sequences used are listed in Table 1. (C) Diagram of the genomic structures of REF genes from Sedum alfredii, Populus trichocarpa, Seteria italic, Arabidopsis thaliana, and Ricinus communis. Solid boxes represent the region coding for the structural sequences of the protein. Thin lines indicate the sequences of introns. The number of nucleotides above the thin lines represent the size of the introns, and the number below the solid boxes the size of the exons. The genomic sequences used in the analysis were listed as follows: Populus trichocarpa REF (PtREF, POPTR_0013s01800g), Seteria italic REF (SiREF, NW_004675963.1), Arabidopsis thaliana REF (AtREF, AT3G05500F22F7.5), and Ricinus communis REF (RcREF, RCOM_1433280).
Rubber elongation factor belongs to the REF family which is classified into three highly related categories including REF, small rubber particle protein (SRPP) and SRP. Like other REFs, SaREFl was similar in terms of structure and characterized by the presence of REF domain covering approximately 225 amino acids. A phylogenetic tree was constructed using the amino acid sequence of SaREFl as well as the counterparts of other representative REFs and REF family members to investigate the relationship among them. Basically, the homologous proteins analyzed were classified into three major groups (REF, SRP and SRPP) although several exceptions existed. SaREFl was clustered together with other REFs being obviously distinct from SRP and SRPP classes and located separately at the base of the subfamily containing most REF proteins (Figure 1B).
Analysis of the genomic structure exhibited that SaREFl gene consisted of three exons and two introns which was in accordance with findings of Hevea brasiliensis REF (Priya et al., 2006). The three exons of 45, 255, and 473 bp, respectively, were interspaced by two introns of 477 and 104 bp, which all began with GT and ended with an AG dinucleotide (Figure 1C). These sequences have been thought to be necessary for correct RNA splicing of various other eukaryotic genes (Breathnach et al., 1978). It was notable that the Populus trichocarpa, Seteria italic, Arabidopsis thaliana, and Ricinus communis REF genes all had three exons and two introns, indicating that SaREFl gene had an exon-intron organization similar to that of representative species of woody and herbaceous plants.
Further comparison with representative members of REF family demonstrated that SaREFl was much more closely identical to REFs comparing with that of SRPs and SRPPs. Among REFs, SaREFl shared with 61.8, 61.3, 61.2, and 61% identity to REFs of R. communis, Populus trichocarpa, Theobroma cacao and Morus alba, respectively, while it had less than 40% identity to the REFs of Brassica campestris and Medicago truncatula (Figure 2).
FIGURE 2. Sequences of SaREFl and proteins belonging to REF family were aligned using the MegAlign program by the CLUSTAL method in the DNASTAR software package. Identical residues are shaded in black and gaps were introduced into sequences to optimize alignment. GenBank accession numbers of the 18 amino acid sequences were listed in Table 1.
Transcript Analysis of SaREFl under Different Heavy-Metal Stress
When S. alfredii Hance (HE) seedlings were challenged with heavy-metal ions, how would SaREFl react? To answer this question, the expression profiles of SaREFl were investigated by a time-course treatment containing eight time points with four heavy metals. Among the four heavy-metal ions (Cd, Zn, Cu, Pb), we found that SaREFl reacted more dramatically to Cd treatment while no obvious up-regulation was observed under Zn treatment (Figure 3). When comparing the transcript levels in root, stem, and leaf under Cd excess, SaREFl showed a higher sensitivity in roots than other tissues, displaying a large-scale elevation of expression from the very beginning, maintaining the trend and reaching the peak at 96 h (Figure 3A). A similar trend was detected in the stems, with a high level of the transcript accumulation from 48 h and keeping rising until the treatment ended (Figure 3A). In the leaves, a relatively steady expression profile was found (Figure 3A).
FIGURE 3. Expression profiles of SaREFl upon treatment of Sedum alfredii seedlings with four heavy-metal ions (Cd, Pb, Cu, Zn) for different amounts of time. (A) Real-time qPCR analysis of SaREFl expression level under cadmium stress. (B) The effect of Pb treatment on SaREFl expression. (C) Expression profile of SaREFl under Cu treatment. (D) Changes of SaREFl transcript level under Zn stress. The time of exposure is indicated on the x-axis; the y-axis indicates the mRNA expression relative to TUB which was selected as an internal reference gene. The normalized mRNA levels without treatment (y-axis “Relative mRNA expression”) were set arbitrarily to 1. Vertical bars represent means ± SD (n = 3). Student’s t-tests was applied to evaluate statistical significance and significant differences (P < 0.05 and P < 0.01) are indicated by * and **respectively.
When the plants were treated with Pb, the expression of SaREFl was up-regulated instantly in the roots and the expression summit occurred at 24 h, then returned back to the level of 0.5 h at 48h and 96 h after treatment (Figure 3B). No dramatic response was detected in the tissues of stem except at 96 h (Figure 3B) and the induced level in leaf showed slight elevation just at 0.5, 6 and 96 h (Figure 3B). The profile of SaREFl responding to Cu stress was basically similar to that under Pb stress with divergences mainly occurring in the responding time points and amplitude of variation (Figure 3C). Interestingly, although S. alfredii was characterized with Zn/Cd hyper-accumulating traits, SaREFl demonstrated a steady-state expression to Zn stress except at three time points in the leaf (Figure 3D).
Based on the above results, we found that the expression of SaREFl was induced by heavy metals and divergent patterns of responses were observed among different metal irons.
SaREFl Enhances Cd Tolerance in the Yeast Assays
To verify the capability of SaREFl to rescue Δycf1 phenotype, yeast expression vector pYES2.1 containing SaREFl was transformed into Δycf1, which does not grow in medium containing 60–80 μM CdCl2 due to the defect of Cd sequestration into the vacuole (Shim et al., 2009). In addition, the EV pYES2.1 transformed yeast cells were considered as control in each experiment. Both Δycf1 strains carrying the SaREFl expression vector and EV pYES2.1 grew well under non-stressed conditions. In contrast, on half-strength SG agar medium supplemented with 30 μM CdCl2, the Δycf1 mutant transformed with pYES2.1-SaREFl construct was able to survive at different dilutions whereas growth of control was arrested (Figure 4A). Similarly, in SG-ura liquid medium supplemented with 30 μM CdCl2, Δycf1 cells expressing SaREFl exhibited dramatically enhanced growth when compared with Δycf1 cells transformed with the EV (Figure 4B). Moreover, transformed INVSc1 yeast cells expressing SaREFl were applied to test the tolerance to Zn2+, Cu2+ and Pb2+. No significant differences were observed comparing the growth of yeast strain carrying empty vector with that of yeast cells transformed with pYES2.1-SaREFl (Supplementary Figure S1). These observations suggest that SaREFl was able to rescue inability of Δycf1 to grow in medium containing a certain amount of cadmium which shows its relevance for cadmium tolerance.
FIGURE 4. Verification of cadmium tolerance of SaREFl in yeast. (A) Exhibition of SaREFl-conferred cadmium resistance in the Δycf1 yeast mutant. Expression of SaREFl brought about enhanced growth of SaREFl-transformed yeast on half-strength SG-ura agar plates supplemented with 30 μM CdCl2 compared with yeast strains transformed with empty vectors. OD600, optical density of the yeast suspension at 600 nm. (B) Time-dependent growth of yeast strains in SG-ura liquid medium supplemented with 30 μM CdCl2. The data are means and SD from three independent experiments.
Transgenic SaREFl Expressing Arabidopsis Lines Displayed Cd Tolerance and Accumulation
To determine whether cadmium tolerance could be enhanced by overexpressing SaREFl in planta, we transformed wild-type Arabidopsis with a construct in which SaREFl was driven by the 35S promoter. Seven independent T3 homozygous transgenic lines were obtained (OE-1 to OE-7) and the expression of SaREFl in different transgenic plants was assessed using RT-PCR. The results showed that the expression of S. alfredii SaREFl gene could be detected in all lines (Figure 5A), and three independent transgenic lines (OE-1, OE-2, and OE-3) were chosen for further tolerance study. Seeds of WT and SaREFl expressing Arabidopsis lines were germinated and grown on ½ MS media with or without cadmium. Seedlings were sampled for the measurement of root length and FW, which were considered as tolerance/sensitivity parameters for WT and transgenic lines. Under normal growing conditions, no significant difference was observed between the WT and the SaREFl expressing seedlings (Figure 5B). When grown on ½ MS media containing 100 μM CdCl2, both WT and SaREFl expressing seedlings exhibited similar reduced growth. However, transgenic lines showed greater tolerance to Cd stress (Figure 5B). Quantitative analysis showed that root length and FW of the SaREFl-expressing seedlings were significantly (P < 0.05) higher than those of WT (Figure 5C). Moreover, the weight of whole seedlings was also measured and more biomass was observed in the transgenic lines in comparison to WT plants grown on medium containing 100 μM CdCl2 (Figure 5D). The above findings suggest that SaREFl-expressing seedlings were more tolerant to cadmium stress comparing to wild type plants.
FIGURE 5. Expression of SaREFl in Arabidopsis confers tolerance for Cd. (A) RT-PCR detected the expression of SaREFl in transgenic plants. (B) Phenotypic changes of WT and transgenic lines (OE-1, OE-2, and OE-3) grown vertically for 10 days on ½ MS media with or without cadmium. (C,D) Root length and fresh weight of plants described in (B). Bars represent the mean of root length of 15–20 plants of each genotype and error bars represent ± SD. Statistical significance was determined by Student’s t-tests and significant differences (P < 0.05 and P < 0.01) are indicated by * and ** respectively. Experiments are performed more than three times and result is the representative of one set.
For the purpose of verifying whether SaREFl could confer cadmium tolerance under natural conditions, we grew WT and SaREFl-expressing Arabidopsis lines on soil for 3 weeks and then irrigated the plants with 500 μM CdCl2 for an additional 10 days. Two physiological indexes of plant indicating cell membrane stability and reactive oxygen species damage were assessed to monitor the growing conditions. The EL of both transgenic lines and WT showed an increase caused by the stress treatment. However, transgenic lines overexpressing SaREFl displayed significantly reduced EL values compared to WT plants (Figure 6A) implying the possible role of SaREFl gene in cell membrane protection. It was further supported by lipid peroxidation analysis, in which high MDA contents were found in WT plants compared to SaREFl-expressing transgenic lines under stress condition (Figure 6B). Furthermore, quantification of Cd contents in the aboveground plant parts was performed to study whether tolerance is accompanied with accumulation of heavy metals. The results showed that the three SaREFl-expressing Arabidopsis lines accumulated more Cd compared to WT (Figure 6C). Taken these results together, these observations supported that expression of SaREFl in Arabidopsis is related to Cd tolerance.
FIGURE 6. Analysis of SaREFl transgenic Arabidopsis lines under cadmium stress condition. (A) Measurement of electrolyte leakage of WT and transgenic lines (OE-1, OE-2, and OE-3); (B) Estimation of MDA in WT and transgenic lines (OE-1, OE-2, and OE-3); (C) Assessment of Cd accumulation in WT and transgenic lines (OE-1, OE-2, and OE-3) under stress conditions. Values are the average of three independent experiments error bars represent _SD (*P, 0.05; **P, 0.01).
Characterization of SaREFl-Interacting Proteins
In order to identify potential proteins interacting with SaREFl, we screened a S. alfredii Hance (HE) cDNA library using yeast two-hybrid technology with SaREFl as bait. Based on the results of self-activation assay, SaREFl did not self-activate the expression of His and Ade reporter genes. An initial screen was performed under stringent selective conditions for HIS3 activity and those colonies that could grow on SD-Leu-Trp-His medium were further tested for the activation of the Ade reporter gene. The screening among cDNA library produced prey inserts of 30 strongly growing colonies (Figure 7B). These His+ and Ade+ positive colonies were sequenced, among which nine colonies contained sequences encoding the prenylated Rab acceptor (PRA1) family proteins.
FIGURE 7. Screening of probable proteins interacting with SaREFl was performed through yeast two-hybrid technique. (A) The diagram displayed continuous numbering representing for different colonies screened by yeast mating. “a” represented for positive control (pGBKT7-53 interacting with pGADT7-T); “b” represented for negative control (pGBKT7-Lam interacting with pGADT7-T); “c” represented for self-activation control (pGBKT7-REF interacting with pGADT7). (B) Screened colonies grew on medium SD-Leu-Trp; (C) Screened colonies grew on medium SD-Leu-Trp-His; (D) Screened colonies grew on medium SD-Leu-Trp-His-Ade.
To determine the positive interactions, the bait (pGBKT7 + SaREFl) and the rescued prey plasmids were re-transformed into yeast Y2HGold. Subsequently, 12 independent colonies that could activate both the His reporter gene and the Ade reporter genes were distinguished from false positive interactions (Figures 7C,D). Among them, eight colonies were sequenced to be the prenylated Rab acceptor (PRA1) family proteins numbered as 3, 8, 9, 12, 14, 15, 17, 18, respectively as depicted in Figure 7A. The rest were COP9 signalosome complex subunit 5 (XP_008445090.1, numbered as 23), REF/SRPP-like protein (XP_008389863.1, numbered as 24), NAD(P)-binding Rossmann-fold superfamily protein isoform 1 (XP_007040707.1, numbered as 26), chromo domain protein (XP_009776768.1, numbered as 29) and binding partner of ACD11 1 (XP_004491782.1, numbered as 30).
Subcellular Localization of SaREFl
The localization of SaREFl protein was predicted using the TargetP program1 which inferred that it was mainly orientated in the cytosol. We fused the coding region of SaREFl at the N-terminus of a green fluorescent protein (GFP) fragment under the control of the CaMV 35S promoter. Subcellular localization of SaREF1-GFP fusion protein was examined by the biolistic assay in onion epidermal cell. Visualized fluorescence indicated that the fluorescent signals of SaREFl-GFP signal were dispersed throughout the cell similar to free GFP (Figures 8D–F), however, some patchy spots in cytosol were detected suggesting potential association with organelles for its cellular functions (Figures 8A–C).
FIGURE 8. Subcellular localization of SaREFl in onion epidermal cells. (A–C) SaREFl fused to GFP was transiently expressed in onion epidermal cells. (D–F) GFP alone was transiently expressed in epidermal cells. (A,D) Bright filed images; (B,E) GFP fluorescence images; (C,F) Artificially merged images. Bars, 100 μm.
Discussion
Heavy-metal environmental pollution is gradually becoming a serious and threating issue confronting society, scientists, and regulators worldwide (Nwaichi et al., 2014). While humans and animals still have an opportunity to avoid heavy metal-contaminated areas as they can move freely, plants cannot. Hence they need to evolve ways to combat the heavy metals they encounter in their living environment. Various defense mechanisms have been developed to prevent heavy-metal-induced injuries, including controlling influx and efflux of metals, binding metals to the cell wall, prohibiting toxic metal accumulation in sensitive tissues, chelation by organic molecules and sequestration in vacuoles (Verbruggen et al., 2009b; Tang et al., 2010; Álvarez et al., 2012). These mechanisms, which help plants to sustain their cellular homeostasis, are the outcome of long-term evolution and natural selection. Evolving to adaption is a path of survival living beings adopted to earn a place in the nature. The hyper-accumulating ecotype of Sedum alfredii Hance was not an exception either. As an organism inhabiting in the deserted Pb/Zn mine district, the traits of heavy-metal hyper-accumulating and resistance were precious gifts endowed by nature. Although no phylogenetic analysis was performed to explore gene evolution in this species, it is hypothesized that there would be numerous genes evolving and obtaining novel functions to cope with the severe environment.
Rubber elongation factor initially received most attention as the major membrane protein of rubber particles and its intervention in the coagulation of rubber particles was regarded as via mechanisms of agglutination or aggregation more than via an enzymatic activity (Berthelot et al., 2012). But there is also no convincing evidence to exclude the possibility of REF participating in the recruitment of rubber enzymes to the surface of latex globules (Berthelot et al., 2012). REF was found to be associated with latex allergy with several IgE epitopes, although its role in the immune response leading to latex hypersensitivity has not yet been determined (Banerjee et al., 2000). Although REF had been found to exhibit significant matches with SRP, the direct role of REF in combating stress has rarely been reported.
In our previous study, a transcript showing high homology to REF was detected to enhance yeast tolerance to cadmium. Preliminary bioinformatics analysis including phylogenetic clustering indicated that it was more identical to REFs in other species rather than SRPs. Therefore, we designated it as SaREFl and focused on the elucidation of the role it might play in cadmium tolerance in this study.
There have been numerous reports demonstrating that expression of genes related to stress could be induced by stress treatment such as transporters, heat-shock proteins, metallothioneins and sulfate-metabolizing proteins (Zhang et al., 2001; Norton et al., 2008; Song et al., 2010; Kumar et al., 2013, 2015; Li et al., 2014). In our study, SaREFl expression profiles were investigated in S. alfredii plants exposed to different heavy metals treatments. It was found that SaREFl exhibited different modes of transcriptional regulation upon heavy metal excess and the most sensitive tissue was the root. The transcript abundance of SaREFl demonstrated strong elevation in response to Cd, moderate up-regulation to Pb and Cu, but no significant change in response to Zn. Based on the above results, it could be observed that although SaREFl was not endowed with a tissue-specific expression pattern (data not shown), however, it did exist a tissue-priority or preference in the response to the stress. This may be partly due to the facts that roots are the primary site of perception and injury when the plants were confronted with stresses. Moreover, SaREFl exhibited a differential regulation to different heavy metal irons. It was somewhat surprising that it showed no dramatic response to Zn stress although S. alfredii Hance (HE) is a hyperaccumulator of Zn and Cd. These findings provide us potential research directions, more significantly, to a large degree, it may be an exhibition of heavy-metal resistance characteristics of SaREFl.
To investigate the cellular function of SaREFl, the gene was expressed in S. cerevisiae strain DTY167 (ycf1), a cadmium-sensitive yeast stain resulted from the deletion of the yeast cadmium factor gene (YCF1) (Szczypka et al., 1994; Li et al., 1996). Using this Cd-sensitive strain (ycf1), it was found that the growth of Δycf1 mutant transformed with the empty vector was inhibited by Cd whereas cells harboring the SaREFl expression vector exhibited better vitality. The results suggested that SaREFl was able to complement Cd sensitivity and rescued or partially rescued the Cd -sensitivity phenotype in the mutant yeast strain, which further supported the role of SaREFl in Cd tolerance. Moreover, metal tolerance assays using the yeast strain INVSc1 showed that expression of the SaREFl gene in yeast cells exhibiting a certain specificity to iron provided no obvious tolerance to Cu2+, Pb2+ or Zn2+. Although it was hard to assess the precise function of SaREFl into specific categories such as uptake or translocation, we inferred that the SaREFl-induced Cd tolerance may be due to the exertion of direct effects to several major pathways or taking parts in assembling other targets that contribute to cadmium tolerance. To address this question, the analysis of differential gene expression patterns present in SaREFl-expressing and wild-type yeast at the whole genome level needed to be investigated, which would also enable the characterization of novel targets of SaREFl and may uncover new aspects of Cd tolerance.
The trait of enhancing cadmium tolerance was further investigated through overexpression of SaREFl in Arabidopsis. Whereas retarded growth occurred in both the wild type and the transgenic lines in Cd-containing media, quantification of the FW and root length of seedlings confirmed that a significant difference could be observed between them, suggesting that the transgenic lines displayed a higher level of tolerance than wild types. Under stress, mature transgenic plant lines overexpressing SaREFl performed much better than wild types as far as the measurements of cell membrane stability indicated by electrolyte leakage and lipid peroxidation (MDA content). Moreover, in the transgenic lines, significantly higher accumulation of Cd was found in comparison to wild type. It is well known that the metabolic homestasis of a cell would be disturbed by abiotic stress resulting in a cascade of events including the enhanced production of ROS and the damage of cell membranes (Verbruggen et al., 2009b; Zwiewka et al., 2015). The findings here seemed to support a role of SaREFl in cadmium tolerance which is possibly achieved through maintaining the metabolic homeostasis during exposure to heavy metal stress.
Although the function of SaREFl in Cd tolerance has been assessed by heterologous expression, how SaREFl acts in cadmium tolerance is still unclear. The confocal microscopic analysis of the SaREFl-GFP protein confirmed that SaREFl was localized to the cytosol, suggesting that it might play roles in intercellular processes and alleviate the cytological effects of CdCl2 stress. An important intracellular metal detoxification is accomplished through metal chelation by peptide ligands such as metallothioneins (MT) and phytochelatins (PCs) (Rascio and Navari-Izzo, 2011). Metallothioneins (MTs) coordinate heavy-metal ions owing to the metal-thiolate bonds contributed by their abundant cysteine residues (Gil-Moreno et al., 2015) while phytochelatins (PCs) are enzymatically synthesized peptides consisting of only three amino acids, glutamine (Glu), cysteine (Cys), and glycine (Gly) with the Glu, and Cys residues linked through a γ-carboxylamide bond (Pal and Rai, 2010). The aminoacid sequence of SaREFl contained only two cysteine residues and lacked metal binding motifs. Hence, it would be hard to classify SaREFl as chelators. In order to characterize possible target genes interacting with SaREFl, we performed a screening among S. alfredii cDNA library using the yeast two-hybrid technology with SaREFl as a bait. We found that a gene named Prenylated Rab acceptor 1 (PRA1) domain protein was detected with high frequency. Numerous studies have shown that prenylated rab acceptor 1 (PRA1) domain proteins may functionally link various vesicle trafficking proteins, regulate vesicle trafficking and participate in protein transport through the Golgi apparatus (Bucci et al., 1999; Lin et al., 2001; Sivars et al., 2003; Kamei et al., 2008; Heo et al., 2010, Jung et al., 2011). Vesicular trafficking is an important step involved in the secretory process (Kohorn, 2000, Willats et al., 2001; Lee et al., 2008; Vassilieva and Nusrat, 2008). This dynamic homeostasis could be influenced by a variety of stimuli, for instance, various stress factors (Örvar et al., 2000; Halperin et al., 2003; Fan et al., 2011). The interaction of SaREFl with PRA suggested that SaREFl might possibly get involved in regulation or recovery of vesicle trafficking to combat the imposed stress. Although it still needs further validation by other protein-interaction techniques, it does provide some insight into the precise mechanism of SaREFl-related to Cd tolerance.
In summary, the present study assessed the characteristics of SaREFl associated with Cd tolerance. Several findings of interest include the different modes of transcriptional regulation of SaREFl upon heavy metal excess. Although there are still many puzzles around SaREFl, it seems that SaREFl is endowed with novel functions alongside the evolution of S. alfredii Hance (HE) being a Cd/Zn co-hyperaccumulator. The findings here could contribute toward enriching the connotation of REF-like genes and provide theoretical guidance for the application in breeding heavy metal-tolerant crop plants.
Author Contributions
RZ and MyL planned and designed the research. MyL, WQ, XS and XH performed experiments. LZ, XH, JJ, GQ, JS and MqL contributed analytical tools. MyL wrote the paper. All authors read and approved the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
This work was supported by National Natural Science Foundation of China (No. 31200465), National Nonprofit Institute Research Grant of CAF (No. TGB2013008, No. RISF2014010), National High Technology Research and Development Program of China (No. 2013AA102701-3).
Supplementary Material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2016.00965
FIGURE S1 | Assessment of metal tolerance of SaREFl to Zn2+, Cu2+ and Pb2+. Two lines of yeast cells expressing SaREFl or empty vector were cultured on medium supplemented with 5 mM Cu2+, 1.5 mM Pb2+ and 15 mM Zn2+.
Footnotes
References
Ahmad, M. S. A., and Ashraf, M. (2011). Essential roles and hazardous effects of nickel in plants. Rev. Environ. Contam. Toxicol. 214, 125–167. doi: 10.1007/978-1-4614-0668-6_6
Altschul, S. F., Madden, T. L., Schäffer, A. A., Zhang, J., Zhang, Z., Miller, W., et al. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389–3402. doi: 10.1093/nar/25.17.3389
Álvarez, R., del Hoyo, A., García-Breijo, F., Reig-Armiñana, J., del Campo, E. M., Guéra, A., et al. (2012). Different strategies to achieve Pb-tolerance by the two Trebouxia algae coexisting in the lichen Ramalina farinacea. J. Plant Physiol. 169, 1797–1806. doi: 10.1016/j.jplph.2012.07.005
Anjum, N. A., Ahmad, I., Mohmood, I., Pacheco, M., Duarte, A. C., Pereira, E., et al. (2012). Modulation of glutathione and its related enzymes in plants’ responses to toxic metals and metalloids-a review. Environ. Exp. Bot. 75, 307–324.
Aoki, Y., Takahashi, S., Takayama, D., Ogata, Y., Sakurai, N., Suzuki, H., et al. (2014). Identification of laticifer-specific genes and their promoter regions from a natural rubber producing plant Hevea brasiliensis. Plant Sci. 225, 1–8. doi: 10.1016/j.plantsci.2014.05.003
Baker, A. J. M., Reeves, R. D., and Hajar, A. S. M. (1994). Heavy-metal accumulation and tolerance in British populations of the metallophyte Thlaspi caerulescens J & C Presl (Brassicaceae). New Phytol. 127, 61–68. doi: 10.1111/j.1469-8137.1994.tb04259.x
Banerjee, B., Kanitpong, K., Fink, J. N., Zussman, M., Sussman, G. L., Kelly, K. J., et al. (2000). Unique and shared IgE epitopes of Hev b 1 and Hev b 3 in latex allergy. Mol. Immunol. 37, 789–798. doi: 10.1016/S0161-5890(00)00095-X
Bert, V., Bonnin, I., Saumitou-Laprade, P., de Laguerie, P., and Petit, D. (2002). Do Arabidopsis halleri from nonmetallicolous populations accumulate zinc and cadmium more effectively than those from metallicolous populations? New Phytol. 155, 47–57. doi: 10.1046/j.1469-8137.2002.00432.x
Berthelot, K., Lecomte, S., Estevez, Y., Coulary-Salin, B., Bentaleb, A., Cullin, C., et al. (2012). Rubber elongation factor (REF), a major allergen component in Hevea brasiliensis latex has amyloid properties. PLoS ONE 7:e48065. doi: 10.1371/journal.pone.0048065
Berthelot, K., Lecomte, S., Estevez, Y., and Peruch, F. (2014). Hevea brasiliensis REF (Hev b 1) and SRPP (Hev b 3): an overview on rubber particle proteins. Biochimie 106, 1–9. doi: 10.1016/j.biochi.2014.07.002
Breathnach, R., Benoist, C. K. O. H., O’hare, K., Gannon, F., and Chambon, P. (1978). Ovalbumin gene: evidence for a leader sequence in mRNA and DNA sequences at the exon–intron boundaries. Proc. Natl. Acad. Sci. U.S.A. 75, 4853–4857. doi: 10.1073/pnas.75.10.4853
Brooks, R. R., Lee, J., Reeves, R. D., and Jaffre, T. (1977). Detection of nickeliferous rocks by analysis of herbarium specimens of indicator plants. J. Geochem. Explor. 7, 49–57. doi: 10.1016/0375-6742(77)90074-7
Bucci, C., Chiariello, M., Lattero, D., Maiorano, M., and Bruni, C. B. (1999). Interaction cloning and characterization of the cDNA encoding the human prenylated rab acceptor (PRA1). Biochem. Biophys. Res. Commun. 258, 657–662. doi: 10.1006/bbrc.1999.0651
Burland, T. G. (2000). DNASTAR’s Lasergene sequence analysis software. Methods Mol. Biol. 132, 71–91.
Cailliatte, R., Schikora, A., Briat, J. F., Mari, S., and Curie, C. (2010). High-affinity manganese uptake by the metal transporter NRAMP1 is essential for Arabidopsis growth in low manganese conditions. Plant Cell 22, 904–917. doi: 10.1105/tpc.109.073023
Dennis, M. S., Henzel, W. J., Bell, J., Kohr, W., and Light, D. R. (1989). Light Amino acid sequence of rubber elongation factor protein associated with rubber particles in Hevea latex. J. Biol. Chem. 264, 18618–18626.
Dennis, M. S., and Light, D. R. (1989). Rubber elongation factor from Hevea brasiliensis. Identification, characterization, and role in rubber biosynthesis. J. Biol. Chem. 264, 18608–18617.
Draper, H. H., and Hadley, M. (1990). Malondialdehyde determination as index of lipid peroxidation. Methods Enzymol. 186, 421–431. doi: 10.1016/0076-6879(90)86135-I
Fan, J. L., Wei, X. Z., Wan, L. C., Zhang, L. Y., Zhao, X. Q., Liu, W. Z., et al. (2011). Disarrangement of actin filaments and Ca2+ gradient by CdCl2 alters cell wall construction in Arabidopsis thaliana root hairs by inhibiting vesicular trafficking. J. Plant Physiol. 168, 1157–1167. doi: 10.1016/j.jplph.2011.01.031
Felsenstein, J. (1993). PHYLIP (Phylogeny Interference Package), version 3.5c, Department of Genetics. Seattle, WA: University of Washington.
Gil-Moreno, S., Jiménez-Martí, E., Palacios,Ò., Zerbe, O., Dallinger, R., Capdevila, M., et al. (2015). Does variation of the inter-domain linker sequence modulate the metal binding behaviour of Helix pomatia Cd-metallothionein? Int. J. Mol. Sci. 17:E6. doi: 10.3390/ijms17010006
Halperin, S. J., Gilroy, S., and Lynch, J. P. (2003). Sodium chloride reduces growth and cytosolic calcium, but does not affect cytosolic pH, in root hairs of Arabidopsis thaliana L. J. Exp. Bot. 54, 1269–1280. doi: 10.1093/jxb/erg134
Han, K. H., Shin, D. H., Yang, J., Kim, I. J., Oh, S. K., and Chow, K. S. (2000). Chow genes expressed in the latex of Hevea brasiliensis. Tree Physiol. 20, 503–510. doi: 10.1093/treephys/20.8.503
Heo, J. B., Bang, W. Y., Kim, S. W., Hwang, S. M., Son, Y. S., Im, C. H., et al. (2010). OsPRA1 plays a significant role in targeting of OsRab7 into the tonoplast via the prevacuolar compartment during vacuolar trafficking in plant cells. Planta 232, 861–871. doi: 10.1007/s00425-010-1226-6
Hoagland, D. R., and Arnon, D. I. (1950). The Water-Culture Method for growing Plants Without Soil, Vol. 347, 2nd Edn. Berkeley, CA: College of Agriculture, 32.
Hsiao, P., Su, R. C., Da Silva, J. A. T., and Chan, M. T. (2007). Plant native tryptophan synthase beta 1 gene is a non-antibiotic selection marker for plant transformation. Planta 225, 897–906. doi: 10.1007/s00425-006-0405-y
Jaffré, T., Brooks, R. R., Lee, J., and Reeves, R. D. (1976). Sebertia acuminata: a hyperaccumulator of nickel from New Caledonia. Science 193, 579–580. doi: 10.1126/science.193.4253.579
Jung, C. J., Hui Lee, M., Ki Min, M., and Hwang, I. (2011). Localization and trafficking of an isoform of the AtPRA1 family to the Golgi apparatus depend on both N- and C-terminal sequence motifs. Traffic 12, 185–200. doi: 10.1111/j.1600-0854.2010.01140.x
Kamei, C. L. A., Boruc, J., Vandepoele, K., Van den Daele, H., Maes, S., Russinova, E., et al. (2008). The PRA1 gene family in Arabidopsis. Plant Physiol. 147, 1735–1749. doi: 10.1104/pp.108.122226
Kohorn, B. D. (2000). Plasma membrane-cell wall contacts. Plant Physiol. 124, 31–38. doi: 10.1104/pp.124.1.31
Kozak, M. (1987). An analysis of 5′-noncoding sequences from 699 vertebrate messenger RNAs. Nucleic Acids Res. 15, 8125–8148. doi: 10.1093/nar/15.20.8125
Krämer, U. (2010). Metal hyperaccumulation in plants. Annu. Rev. Plant Biol. 61, 517–534. doi: 10.1146/annurev-arplant-042809-112156
Kumar, S., Asif, M. H., Chakrabarty, D., Tripathi, R. D., Dubey, R. S., and Trivedi, P. K. (2013). Differential expression of rice Lambda class GST gene family members during plant growth, development, and in response to stress conditions. Plant Mol. Biol. Rep. 31, 569–580. doi: 10.1007/s11105-012-0524-5
Kumar, S., Asif, M. H., Chakrabarty, D., Tripathi, R. D., Dubey, R. S., and Trivedi, P. K. (2015). Comprehensive analysis of regulatory elements of the promoters of rice sulfate transporter gene family and functional characterization of OsSul1; 1 promoter under different metal stress. Plant Signal. Behav. 10:e990843. doi: 10.4161/15592324.2014.990843
Kush, A., Goyvaerts, E., Chye, M. L., and Chua, N. H. (1990). Laticifer specific gene expression in Hevea brasiliensis (rubber tree). Proc. Natl. Acad. Sci. U.S.A. 87, 1787–1790. doi: 10.1073/pnas.87.5.1787
Lee, Y. J., Szumlanski, A., Nielsen, E., and Yang, Z. (2008). Rho-GTPase-dependent filamentous actindynamics coordinate vesicle targeting and exocytosis during tip growth. J. Cell Biol. 181, 1155–1168. doi: 10.1083/jcb.200801086
Li, J., Liu, B., Cheng, F., Wang, X., Aarts, M. G., and Wu, J. (2014). Expression profiling reveals functionally redundant multiple-copy genes related to zinc, iron and cadmium responses in Brassica rapa. New Phytol. 203, 182–194. doi: 10.1111/nph.12803
Li, Z. S., Szczypka, M., Lu, Y. P., Thiele, D. J., and Rea, P. A. (1996). The yeast cadmium factor protein (YCF1) is a vacuolar glutathione S-conjugate pump. J. Biol. Chem. 271, 6509–6517. doi: 10.1074/jbc.271.11.6509
Lin, J., Liang, Z., Zhang, Z., and Li, G. (2001). Membrane topography and topogenesis of prenylated Rab acceptor (PRA1). J. Biol. Chem. 276, 41733–41741. doi: 10.1074/jbc.M103475200
Lin, Y. F., and Aarts, M. G. (2012). The molecular mechanism of zinc and cadmium stress response in plants. Cell. Mol. Life Sci. 69, 3187–3206. doi: 10.1007/s00018-012-1089-z
Liu, M., Qiao, G., Jiang, J., Han, X., Sang, J., and Zhuo, R. (2014). Identification and expression analysis of salt-responsive genes using a comparative microarray approach in Salix matsudana. Mol. Biol. Rep. 41, 6555–6568. doi: 10.1007/s11033-014-3539-1
Milner, M. J., and Kochian, L. V. (2008). Investigating heavy-metal hyperaccumulation using Thlaspi caerulescens as a model system. Ann. Bot. 102, 3–13. doi: 10.1093/aob/mcn063
Murashige, T., and Skoog, F. (1962). A revised medium for rapid growth and bio assays with tobacco tissue cultures. Physiol. Plant. 15, 473–497. doi: 10.1111/j.1399-3054.1962.tb08052.x
Nielsen, H., Engelbrecht, J., Brunak, S., and von Heijne, G. (1997). Identification of prokaryotic and eukaryotic signal peptides and prediction of their cleavage sites. Protein Eng. 10, 1–6. doi: 10.1093/protein/10.1.1
Norton, G. J., Lou-Hing, D. E., Meharg, A. A., and Price, A. H. (2008). Rice-arsenate interactions in hydroponics: whole genome transcriptional analysis. J. Exp. Bot. 59, 2267–2276. doi: 10.1093/jxb/ern097
Nwaichi, E. O., Wegwu, M. O., and Nwosu, U. L. (2014). Distribution of selected carcinogenic hydrocarbon and heavy metals in an oil-polluted agriculture zone. Environ. Monit. Assess. 186, 8697–8706. doi: 10.1007/s10661-014-4037-6
Örvar, B. L., Sangwan, V., Omann, F., and Dhindsa, R. S. (2000). Early steps in cold sensing by plant cells: the role of actin cytoskeleton and membrane fluidity. Plant J. 23, 785–794. doi: 10.1046/j.1365-313x.2000.00845.x
Pal, R., and Rai, J. P. (2010). Phytochelatins: peptides involved in heavy metal detoxification. Appl. Biochem. Biotechnol. 160, 945–963. doi: 10.1007/s12010-009-8565-4
Pollard, A. J., Powell, K. D., Harper, F. A., and Smith, J. A. C. (2002). The genetic basis of metal hyperaccumulation in plants. Crit. Rev. Plant Sci. 21, 539–566. doi: 10.1080/0735-260291044359
Priya, P., Venkatachalam, P., and Thulaseedharan, A. (2006). Molecular cloning and characterization of the rubber elongation factor gene and its promoter sequence from rubber tree (Hevea brasiliensis): a gene involved in rubber biosynthesis. Plant Sci. 171, 470–480. doi: 10.1016/j.plantsci.2006.05.009
Qi, X., Cui, Q., Luo, Y., Guo, C., and Chai, T. (2009). Zn stress-induced inhibition of bean PvSR2-GUS fusion gene splicing is gene-specific in transgenic tobacco. J. Plant Physiol. 166, 1223–1227. doi: 10.1016/j.jplph.2009.01.010
Rascio, N., and Navari-Izzo, F. (2011). Heavy metal hyperaccumulating plants: how and why do they do it? And what makes them so interesting? Plant Sci. 180, 169–181. doi: 10.1016/j.plantsci.2010.08.016
Rehm, B. (2001). Bioinformatic tools for DNA/protein sequence analysis, functional assignment of genes and protein classification. Appl. Microbiol. Biotechnol. 57, 579–592. doi: 10.1007/s00253-001-0844-0
Sang, J., Han, X., Liu, M., Qiao, G., Jiang, J., and Zhuo, R. (2013). Selection and validation of reference genes for real-time quantitative PCR in hyperaccumulating ecotype of Sedum alfredii under different heavy metals stresses. PLoS ONE 8:e82927. doi: 10.1371/journal.pone.0082927
Shahid, M., Austruy, A., Echevarria, G., Arshad, M., Sanaullah, M., Aslam, M., et al. (2014). EDTA-enhanced phytoremediation of heavy metals: a review. Soil Sediment Contam. Int. J. 23, 389–416. doi: 10.1080/15320383.2014.831029
Shim, D., Hwang, J. U., Lee, J., Lee, S., Choi, Y., An, G., et al. (2009). Orthologs of the class A4 heat shock transcription factor HsfA4a confer cadmium tolerance in wheat and rice. Plant Cell 21, 4031–4043. doi: 10.1105/tpc.109.066902
Sivars, U., Aivazian, D., and Pfeffer, S. R. (2003). Yip3 catalyses the dissociation of endosomal Rab-GDI complexes. Nature 425, 856–859. doi: 10.1038/nature02057
Song, W. Y., Park, J., Mendoza-Cózatl, D. G., Suter-Grotemeyer, M., Shim, D., Hörtensteiner, S., et al. (2010). Arsenic tolerance in Arabidopsis is mediated by two ABCC-type phytochelatin transporters. Proc. Natl. Acad. Sci. U.S.A. 107, 21187–21192.
Straif, K., Benbrahim-Tallaa, L., Baan, R., Grosse, Y., Secretan, B., El Ghissassi, F., et al. (2009). A review of human carcinogens-part C: metals, arsenic, dusts, and fibres. Lancet Oncol. 10, 453–454. doi: 10.1016/S1470-2045(09)70134-2
Su, C., Jiang, L., and Zhang, W. (2014). A review on heavy metal contamination in the soil worldwide: situation, impact and remediation techniques. Environment. Skeptics Critics 3, 24–38.
Sun, X. H., Yu, G., Li, J. T., Jia, P., Zhang, J. C., Jia, C. G., et al. (2014). A heavy metal-associated protein (AcHMA1) from the halophyte, Atriplex canescens (Pursh) Nutt., confers tolerance to iron and other abiotic stresses when expressed in Saccharomyces cerevisiae. Int. J. Mol. Sci. 15, 14891–14906. doi: 10.3390/ijms150814891
Sun, Y. B., Zhou, Q. X., Liu, W. T., and Liu, R. (2009). Chelator-enhanced phytoextraction of heavy metals from contaminated soil irrigated by industrial wastewater with the hyperaccumulator plant (Sedum alfredii Hance). Geoderma 150, 106–112. doi: 10.1016/j.geoderma.2009.01.016
Szczypka, M. S., Wemmie, J. A., Moye-Rowley, W. S., and Thiele, D. J. (1994). A yeast metal resistance protein similar to human cystic fibrosis transmembrane conductance regulator (CFTR) and multidrug resistance-associated protein. J. Biol. Chem. 269, 22853–22857.
Tang, K., Zhan, J. C., Yang, H. R., and Huang, W. D. (2010). Changes of resveratrol and antioxidant enzymes during UV-induced plant defense response in peanut seedlings. J. Plant Physiol. 167, 95–102. doi: 10.1016/j.jplph.2009.07.011
Thomine, S., Lelièvre, F., Debarbieux, E., Schroeder, J. I., and Barbier-Brygoo, H. (2003). AtNRAMP3, a multispecific vacuolar metal transporter involved in plant responses to iron deficiency. Plant J. 34, 685–695. doi: 10.1046/j.1365-313X.2003.01760.x
Vassilieva, E. V., and Nusrat, A. (2008). Vesicular trafficking: molecular tools and targets. Methods Mol. Biol. 440, 3–14. doi: 10.1007/978-1-59745-178-9_1
Verbruggen, N., Hermans, C., and Schat, H. (2009a). Mechanisms to cope with arsenic or cadmium excess in plants. Curr. Opin. Plant Biol. 12, 364–372. doi: 10.1016/j.pbi.2009.05.001
Verbruggen, N., Hermans, C., and Schat, H. (2009b). Molecular mechanisms of metal hyperaccumulation in plants. New Phytol. 181, 759–776. doi: 10.1111/j.1469-8137.2008.02748.x
Willats, W. G., McCartney, L., and Knox, J. P. (2001). In situ analysis of pectic polysaccharides in seed mucilage and at the root surface of Arabidopsis thaliana. Planta 213, 37–44. doi: 10.1007/s004250000481
Yadav, S. K. (2010). Heavy metals toxicity in plants: an overview on the role of glutathione and phytochelatins in heavy metal stress tolerance of plants. South Afr. J. Bot. 76, 167–179. doi: 10.1016/j.sajb.2009.10.007
Yang, X. E., Long, X., Ni, W., and Fu, C. (2002). Sedum alfredii H: a new Zn hyperaccumulating plant first found in China. Chin. Sci. Bull. 47, 1634–1637. doi: 10.1111/nph.12168
Yang, X. E., Long, X. X., Ye, H. B., He, Z. L., Calvert, D. V., and Stoffella, P. J. (2004). Cadmium tolerance and hyperaccumulation in a new Zn-hyperaccumulating plant species (Sedum alfredii Hance). Plant Soil 259, 181–189. doi: 10.1023/B:PLSO.0000020956.24027.f2
Zhang, Y., Chai, T. Y., Dong, J., Zhao, W., An, C. C., Chen, Z. L., et al. (2001). Cloning and expression analysis of the heavy-metal responsive gene PvSR2 from bean. Plant Sci. 161, 783–790. doi: 10.1016/S0168-9452(01)00470-8
Keywords: cadmium, Cd tolerance, hyperaccumulator, rubber elongation factor (REF), Sedum alfredii Hance
Citation: Liu M, Qiu W, He X, Zheng L, Song X, Han X, Jiang J, Qiao G, Sang J, Liu M and Zhuo R (2016) Functional Characterization of a Gene in Sedum alfredii Hance Resembling Rubber Elongation Factor Endowed with Functions Associated with Cadmium Tolerance. Front. Plant Sci. 7:965. doi: 10.3389/fpls.2016.00965
Received: 09 December 2015; Accepted: 16 June 2016;
Published: 29 June 2016.
Edited by:
David W. M. Leung, University of Canterbury, New ZealandReviewed by:
Woei-Jiun Guo, National Cheng Kung University, TaiwanUener Kolukisaoglu, University of Tübingen, Germany
Copyright © 2016 Liu, Qiu, He, Zheng, Song, Han, Jiang, Qiao, Sang, Liu and Zhuo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Renying Zhuo, zhuory@gmail.com