Tianjie Sun1,2,3†
Nan Ma1,2,3†
Caiqing Wang1,2,3†Huifen Fan1,2,3Mengxuan Wang1,2,3Jie Zhang1,2,3*Jinfeng Cao4,5*Dongmei Wang1,2,3*- 1State Key Laboratory of North China Crop Improvement and Regulation, Baoding, China
- 2Key Laboratory of Hebei Province for Plant Physiology and Molecular Pathology, Baoding, China
- 3College of Life Sciences, Hebei Agricultural University, Baoding, China
- 4Hebei Key Laboratory of Crop Salt-Alkali Stress Tolerance Evaluation and Genetic Improvement, Cangzhou, China
- 5Academy of Agricultural and Forestry Sciences, Cangzhou, China
Salt stress caused by soil salinization, is one of the main factors that reduce soybean yield and quality. A large number of genes have been found to be involved in the regulation of salt tolerance. In this study, we characterized a soybean sodium/hydrogen exchanger gene GmNHX5 and revealed its functional mechanism involved in the salt tolerance process in soybean. GmNHX5 responded to salt stress at the transcription level in the salt stress-tolerant soybean plants, but not significantly changed in the salt-sensitive ones. GmNHX5 was located in the Golgi apparatus, and distributed in new leaves and vascular, and was induced by salt treatment. Overexpression of GmNHX5 improved the salt tolerance of hairy roots induced by soybean cotyledons, while the opposite was observed when GmNHX5 was knockout by CRISPR/Cas9. Soybean seedlings overexpressing GmNHX5 also showed an increased expression of GmSOS1, GmSKOR, and GmHKT1, higher K+/Na+ ratio, and higher viability when exposed to salt stress. Our findings provide an effective candidate gene for the cultivation of salt-tolerant germplasm resources and new clues for further understanding of the salt-tolerance mechanism in plants.
Introduction
Soybean is one of the most widely cultivated oil crops worldwide. The usage of soybean products includes human foods, animal foods, industrial products, ingredients, precursor materials, etc. (Gaonkar and Rosentrater, 2019). The production of soybean is restricted by a variety of adverse factors, among which salt stress caused by soil salinization is one of the most serious factors that harm soybean yield and quality (Hasanuzzaman et al., 2016). Salt stress impairs soybean growth in many ways. On the one hand, salt stress inhibits the formation of nodules to weaken the assimilation of nitrogen (Singleton and Bohlool, 1984); on the other hand, high osmotic pressure causes cells to lose water and suffer toxic effects (Massa and Melito, 2019). Salt stress disrupts the dynamic equilibrium established by the reactive oxygen species (ROS)-scavenging system (Moradi and Ismail, 2007; Esfandiari and Gohari, 2017), causing the increase of membrane permeability and malondialdehyde (MDA) content due to the oxidation of membrane lipids (Geilfus, 2019).
Plants have evolved multiple ways to enhance their adaptability under salt stress. Cells selectively absorb salt and compartmentalize salt into vacuoles to maintain a relatively stable internal environment, or transport excess salt to older tissues to avoid damage to the young tissues that are more sensitive to salt stress (Hines, 2008; Wang et al., 2019). The inorganic ions absorbed by plants mainly comprise Na+, K+, Ca2+, Cl–, among which the balance of intracellular K+/Na+ concentration is the key to ensure normal physiological metabolism (Maathuis and Amtmann, 1999; James et al., 2011). In addition to the regulation of ion levels, plants also accumulate a series of osmotic adjustment substances, such as soluble saccharide, proline, and betaine (Qi et al., 2014), to regulate the osmotic potential and prevent the damage caused by ROS (Bohnert and Jensen, 1996).
Salt tolerance is an important trait in soybean breeding, but there have not been many salt-tolerant genetic resources available for breeding, which limits the development of new salt-tolerant soybean varieties. Therefore, it is of great significance to explore the molecular basis of salt tolerance in plants and to develop genetic resources that are available for genetic breeding. Among the molecular pathways that are involved in salt stress response, including protein kinases such as CIPK (Weinl and Kudla, 2009) and CDPK (Boudsocq and Sheen, 2012); transcription factors such as NAC (Tran et al., 2004), WRKY (Jiang and Deyholos, 2008), bHLH (Jiang et al., 2009), SERF (Schmidt et al., 2013), bZIP (Yang et al., 2009), and MYB (Yoo et al., 2005); channel proteins such as HKT (Horie et al., 2007) and AKT (Wang P. et al., 2015); ion pumps such as HAK (Wang et al., 2018), and ion exchangers such as NHX (Aviv, 1993) and SOS (Qiu et al., 2002), have been found to respond to salt stress and participate in the regulation of plant salt tolerance. The functions of NHX (Na+/H+ exchanger) family members have attracted growing attention of physiology researchers.
NHXs are monovalent ion exchangers that localize in the membrane, which catalyze the movement of Na+ or K+ (Barragán et al., 2012) to the side with high ion concentration, and at the same time exchange with an H+ across the membrane to maintain local potential conservation (Rodríguez-Rosales et al., 2009). Such functions maintain the ion concentration of cells (or organelles) on the one hand, and regulate the stability of the pH environment of the chamber on the other hand (Gadsby, 2009). In Arabidopsis, the NHX family consists of eight members, which can be divided into three categories according to their subcellular distribution: plasma membrane NHXs (AtNHX7/SOS1 and AtNHX8), vacuolar NHXs (AtNHX1-4), and inner membrane NHXs (AtNHX5 and AtNHX6) (Qiu, 2016). Up to now, most of the studies on the function of NHX family genes and their involvement in signaling pathways have focused on vacuolar membrane NHXs and plasma membrane NHXs. At the histological level, the expression of SOS1 has been detected in epidermal cells at the root tip and in parenchyma cells at the xylem/symplast boundary of roots, stems, and leaves, controlling Na+ transport through the system of xylem system (Shi et al., 2002). At the cellular level, AtNHX1 transports Na+ (Pehlivan et al., 2016) or K+ (Barragán et al., 2012) to the vacuole in exchange for H+ to the cytoplasm, while SOS1 exports Na+ to Pehlivan et al. (2016), or import K+ from Wu et al. (1996), the extracellular space, thereby increase the K+/Na+ ratio in the cytoplasm. Our previous studies also demonstrated that heterologous expression of GmNHX1 (a soybean NHX that localizes in the vacuole membrane) in Arabidopsis increased the flow rate of Na+ in the roots by up-regulating the expression of SKOR, SOS1, and AKT1 genes under salt stress, thereby improve salt tolerance via increasing the K+/Na+ ratio (Sun et al., 2019).
To our knowledge, however, research on the NHXs distributed in the endomembrane system is rare. AtNHX5 and AtNHX6 are located in the trans-Golgi network and pre-vacuolar compartments, affecting the cargo sorting behavior of the trans-Golgi network by adjusting the pH in the vesicles (Reguera et al., 2015). AtNHX5 and AtNHX6 also maintain the steady-state of K+ and pH, allowing cells to accumulate higher K+ under acidic conditions (Wang L. et al., 2015). The nhx5 nhx6 double mutant showed intense sensitivity to salt stress (Bassil et al., 2011a), and the heterologous expression of AtNHX5 in Broussonetia papyrifera increases the tolerance to salt stress and drought (Li et al., 2011). It is worth noting that the molecular mechanism by which these NHXs related to the inner membrane system and participate in plant salt tolerance has not yet been fully revealed.
In this study, we characterize GmNHX5, a Golgi apparatus-localized NHX in soybean, and reveal the mechanism by which it participates in salt stress adaption in plants. The transcription level of GmNHX5 increases significantly upon salt treatment in salt-tolerant soybean cultivars, whereas remains unchanged in salt-sensitive ones. The activity of the GmNHX5 promoter is high in new leaves and vascular tissues and increases when exposed to salt stress. Hairy roots overexpressing GmNHX5 show stronger salt tolerance, while the GmNHX5 knocked-out hairy roots are more sensitive to salt stress than wild type. We further overexpress GmNHX5 in soybean and find that the transformed plants show higher osmotic adjustment substances accumulation, K+/Na+ ratio, stronger salt tolerance, and higher expression levels of GmSOS1, GmSKOR, and GmHKT1 under salt stress than wild type. Our findings provide a candidate gene for genetic breeding and offer clues to the molecular mechanism by which endometrial NHXs participate in the response to salt stress.
Materials and Methods
Plant Cultivation
Soybean cultivar Jidou-7 was obtained from the Institute of Grain and Oil Crops, Hebei Academy of Agricultural and Forestry Sciences, and cultivar Mustang was obtained from Plant Transformation Core Facility, University of Missouri. Soybean plants were cultivated in a greenhouse with a 14 h light/10 h dark cycle at a constant temperature of 25°C and 700 μmol photons m–2 s–1. Soybeans hairy roots were cultured in sterile petri dishes in a tissue culture room with the same light and temperature conditions as the greenhouse.
Isolation of Plant RNA and Detection of Gene Expression
Total RNA was isolated using an RNA extraction kit (UNlQ-10 Column Trizol Total RNA Isolation Kit, Sangon Biotech), then reverse transcribed into cDNA through a reverse transcription kit (PrimeScript RT reagent Kit with gDNA Eraser, TaKaRa). RT-PCR and RT-qPCR experiments were performed using primers listed in Supplementary Table 1. Actin was used as a reference gene in the RT-PCR and RT-qPCR experiments according to our previous report (Ma et al., 2013). Three biological replicates with their respective three technical replicates were conducted for each sample in the RT-qPCR experiments and calculated using the 2–ΔΔCT method.
Phylogenetic Analysis
Protein sequences were retrieved from UniProtKB through the accessions in Supplementary Table 2 before aligned using ClustalW (Thompson et al., 1994). The similarity between proteins was inferred by using the Maximum Likelihood method and JTT matrix-based model (Jones et al., 1992). The initial tree for the heuristic search was obtained automatically by applying Neighbor-Join and BioNJ algorithms to a matrix of pairwise distances estimated using a JTT model and then selecting the topology with superior log likelihood value. The tree was drawn to scale, with branch lengths measured in the number of substitutions per site. Evolutionary analyses were conducted in MEGA X (Kumar et al., 2018). The proteins in the phylogenetic tree were further colored based on their subcellular distribution analyzed based on their Gene Ontology (GO) terms (Gaudet et al., 2011).
Plasmid Construction
In order to construct a plasmid for histological analysis of GmNHX5, purified DNA from soybean variety Jidou-7 was used as a template, and primer 1381-GmNHX5p-F/R was used to amplify the 2000 bp sequence upstream of the GmNHX5 start codon. The amplified fragment was inserted between EcoRI and HindIII in plasmid pCAMBIA1381.
For the construction of the plasmid for GmNHX5 subcellular localization, primers GmNHX5GFP-F/R were used to amplify the CDS region of the GmNHX5 gene using the cDNA of soybean cultivar Jidou-7 as a template, and the amplified fragment was inserted into the plasmid pCAMBIA3301-GFP between NcoI and EcoO65I.
The plasmid used to silence GmNHX5 by VIGS was constructed by inserting a specific fragment of GmNHX5 amplified using the primers GmNHX5VIGS-F/R and inserted into the pTRV2 plasmid between EcoRI and XhoI.
For the construction of the GmNHX5-overexpressing plasmid, the primers GmNHX5OE-F/R were used to amplify the CDS region of the GmNHX5 gene using the cDNA of soybean cultivar Jidou-7 as a template, and the amplified fragment was inserted into the plasmid pCAMBIA3301 between NcoI and EcoO65I.
Targeted gene knockout was performed using the CRISPR/Cas9-based approach. The gene sequence of GmNHX5 with introns was aligned with the whole soybean genome (Wm82.a2.v1, Phytozome) (Schmutz et al., 2010), and the appropriate sgRNA was designed according to a previously published protocol (Xie et al., 2014). Two specific sgRNAs (sgRNA-1: 5′-CAGACACCGAGACTAATATC-3′; sgRNA-2: 5′-ACTCCTTTAGTACTCAGTCT-3′) were used for GmNHX5 knockout, according to the previously published method (Xing et al., 2014). To construct the GmNHX5 knockout plasmid, the primers GmNHX5-DT1-BsF/F0/R0/BsR were used to amplify the plasmid pCBC-DT1T2 and inserted the amplified fragment into the Cas9 expression vector pBSE401 according to a previously reported protocol (Xing et al., 2014). All primers are listed in Supplementary Table 1.
Detection of GmNHX5 Histological Distribution
The histological distribution of GmNHX5 promoter activity was detected by GUS staining. The T2 progeny of 5-day-old Arabidopsis seedlings transformed with GmNHX5p:GUS plasmid were transferred to MS medium containing 0, 100, or 200 mM NaCl, respectively, for 24 h. Histochemical GUS staining assays were performed using a GUS histochemical assay kit (Real-Times) following the manufacturer’s protocol. The seedlings were washed with 70% ethanol to dissolve chlorophyll and GUS staining was examined under a fluorescence microscope (BX53, Olympus). These experiments were repeated three times to ensure accuracy.
Detection of GmNHX5 Subcellular Localization
The subcellular localization of GmNHX5 was observed using laser scanning confocal microscopy. The T2 progeny of Arabidopsis seedlings transformed with plasmids carrying CaMV 35S:GmNHX5-GFP or CaMV 35S:GFP expression loci were observed under a laser scanning confocal microscope (FV1000, Olympus). Fluorochrome Golgi-Tracker Red (Beyotime) was applied to the roots of 5-day-old Arabidopsis seedlings to label the Golgi apparatus following the manufacturer’s protocol. Fluorescence of GFP was excited at 488 nm and detected at 520–540 nm. Fluorescence of Golgi-Tracker Red was excited at 589 nm and detected at 605–635 nm. Line profiling was performed with Fiji (ImageJ1), in accordance with a previous study (Li et al., 2012). These experiments were repeated three times to ensure accuracy.
Generation of GmNHX5-Silencing Plants
TRV-mediated VIGS was used to generate GmNHX5-silencing plants. The VIGS experiment was conducted according to a previously reported protocol (Muthappa and Mysore, 2014), with minor modifications. One-week-old soybean seedlings were used in the VIGS experiment. The empty plasmid pTRV2 was used as control (TRV: 00). Agrobacterium tumefaciens carrying pTRV1 and recombinant pTRV2 plasmids were resuspended in the infection buffer [50 mM MES, 2 mM Na3PO4, 28 mM D-glucose, 0.1 mM acetosyringone, 4.1 mM L-Cys, 0.02% (w/v) Silwet L-77] to OD600 = 0.5, mixed with 1:1 (v:v), and poured into the root of soybean seedlings with 5 ml per plant. The treatment was performed three times, with an interval of 5 days each time. RT-qPCR was used to examine the silencing effect in the first leaves in soybean seedlings 3 weeks after initial Agrobacterium tumefaciens treatment before subjected to salt stress treatment.
Hairy Roots Induction and Biomass Statistics
Hairy roots were prepared using Agrobacterium rhizogenes-mediated soybean cotyledon induction according to previously published protocols (Chen et al., 2014; An et al., 2017) with minor modification. After surface sterilization for 1 h using chlorine vapor, the soybean cultivar Mustang seeds were placed on germination medium (Gamborg’s B5 medium with 87.64 mM sucrose, 1.88 mM MES, and 8 g/L agar, pH 5.75) and cultivated at 25°C and 14/10 (L/D) photoperiod for 7 days. Agrobacterium rhizogenes strain K599 carrying the recombinant plasmid for GmNHX5-overexpression or knockout were used for soybean hairy root induction. Soybean seedlings were placed in a sterile petri dish, the cotyledons were cut off, and the cotyledons were soaked in a solution (87.64 mM sucrose, 1.88 mM MES, 0.2 mM acetosyringone, 1.66 mM L-cystine, and 1 mM DTT, pH 5.4) containing corresponding Agrobacterium rhizogenes strains at OD600 = 0.5 for 30 min, and cultured on sterile filter paper in the dark at 25°C for 5 days. The cotyledons were washed with sterile water and transferred to the induction medium (Gamborg’s B5 medium with 43.82 mM sucrose, 1.88 mM MES, 400 mg/L cefotaxime sodium, 100 mg/L carbenicillin disodium, 100 mg/L vancomycin hydrochloride, 100 mg/L timentin, and 8 g/L agar, pH 5.75) and incubated in the dark at 25°C for 15 days.
Hairy roots were cut off and used to detect the expression of GmNHX5 (for OE hairy roots) or gene editing (for KO hairy roots), and about 1 cm of hairy roots were reserved for subsequent experiments. The DNA of the transformed hairy roots was amplified by GmNHX5KD-F/R and sequenced to confirm the knockout effect. The partially cut hairy roots were used for RT-qPCR analysis or sequencing, and the remaining hairy roots (about 1 cm in length) with their cotyledons were continued to be cultured in the induction medium containing 0, 70, and 150 mM NaCl for another 15 days before analysis. For hairy root biomass statistics, the maximum root length was measured as the highest length of hairy roots produced on each cotyledon. Fresh weight was measured as the weight of all hairy roots for each cotyledon.
Generation of GmNHX5-Overexpressing Soybean Plants
Soybean plants overexpressing GmNHX5 were generated using Agrobacterium tumefaciens- mediated cotyledonary node transformation following a previously published protocol (Paz et al., 2004; Lee et al., 2013), and the soybean cultivar Mustang was used as the transformation recipient. Agrobacterium tumefaciens strain EHA105 carrying plasmid for GmNHX5-overexpression was used for soybean transformation. The glyphosate-screened soybean seedlings were further verified using PCR (with primers BarS/X), colloidal gold test, and Southern blot analysis to confirm the transformation event.
Colloidal Gold Test
The colloidal gold test was performed using colloidal gold test strips (QuickStix for PAT/bar, Envirologix) following the manufacturer’s protocol.
Southern Blot
DNA isolated from GmNHX5-overexpressing soybean seedling leaves were used for Southern blot analysis. The digoxigenin-labeled probe was synthesized using primers BarS/X with a labeling kit (PCR DIG probe synthesis kit, Roche) following the manufacturer’s protocol. The DNA samples digested with HindIII, together with a molecular marker (DIG-labeled DNA molecular weight marker II, Roche), were separated on a 1% agarose gel, transferred to a positively charged nylon membrane, and subjected to UV cross-linking and probe hybridization. The labeled DNA fragments were developed using a DIG detection kit (DIG-High prime DNA labeling and detection starter kit I, Roche).
Phenotyping and Physiological Analysis
Three-week-old WT and GmNHX5-overexpressing soybean seedlings were irrigated with Hoagland’s nutrient solution with 0 or 150 mM NaCl for 20 days, leaf dry weight and leaf fresh weight were measured as the total weight of fresh or kiln-dried leaves from single seedlings.
Soybean seedlings were exposed to 0 or 150 mM NaCl for 6 days before determination of free proline, betaine, and malondialdehyde (MDA). Free proline content was measured according to a previously published protocol (Bates et al., 1973). Betaine was extracted with ODS-C18 column and measured using an HPLC based method (Schmedes and Ga, 1989). The MDA content was determined according to a thiobarbituric acid (TBA) method (Schmedes and Ga, 1989).
Determination of Na+ and K+ Ion Content
Three-weeks-old seedlings of GmNHX5-overexpressing T3 soybean lines #3 and #4, and untransformed soybean cultivar Mustang was used for Na+ and K+ content determination. Seedlings were exposed to 0 or 150 mM NaCl for 6 h before analysis.
Ion content determination was performed as previously described (Sun et al., 2019), with minor modification. The surface of the leaves, roots, and stems of materials post-treatment were washed with sterile water and were baked in an oven at 120°C for 30 min, and then dried at 80°C to a constant weight. The dried soybean material was ground into a fine powder, adding 10 mL anhydrous sulfuric acid, boiled at 170°C for 20 min, 220°C for 40 min, and 320°C for 2 h until the solution became clear. During the boiling process, a small amount of 30% hydrogen peroxide was added intermittently until a large amount of smoke appeared. After cooling, the product was diluted to 50 ml, and 100 μL of the diluted product was further diluted to 6 ml with distilled water. The diluted sample was filtered through a 0.45 μm membrane, and the content of K+ and Na+ ions was measured with an atomic absorption spectrophotometer (ZA3000, Hitachi).
Accessions
Sequence information of genes involved in this study can be found in NCBI2 through following accessions. XM_006597582.2 (GmNHX5); NM_001258010.2 (GmSOS1); XM_003544313.4 (GmSKOR); XM_014764887.2 (GmHKT1); XM_003545450.3 (GmAKT1); XM_003553762.3 (GmHAK5); XM_026123971.1 (GmNHX1); XM_026124165.1 (Actin).
Statistical Analysis
Statistical calculations were done with Microsoft Excel. Detailed statistical analysis is described in legend of each figure.
Results
GmNHX5 Responds to Salt Stress
The Arabidopsis nhx5 nhx6 double mutant shows high sensitivity to salt stress, indicating that Golgi-localized NHXs may be involved in the tolerance of salt stress in plants (Bassil et al., 2011a). Phylogenetic analysis showed that GmNHX5 shared high sequence identity with AtNHX5 and AtNHX6 (Figure 1A), suggesting that GmNHX5 may also be localized in the Golgi apparatus and is relevant to salt tolerance. GmNHX5 possesses a typical sodium/hydrogen exchanger domain and has 10 transmembrane domains (Figure 1B), which indicates that GmNHX5 is also a membrane protein and may have similar functions to other NHXs, namely regulating the steady-state of H+ and Na+ on both sides of the membrane.
Figure 1. Characterization and expression of GmNHX5. (A) Phylogenetic analysis of GmNHX5 and its homologs from other plants. Protein sequences of NHXs from several plants were used for phylogenetic analysis using the maximum likelihood method and JTT matrix-based model. Species abbreviated as: Ta, Triticum aestivum; Os, Oryza sativa; Zm, Zea mays; Gm, Glycine max; Gs, Glycine soja; At, Arabidopsis thaliana. Red dot emphasizes GmNHX5 used in this study. Colored sections indicate potential subcellular localization of corresponding proteins retrieved from UniProtKB. The scale bar indicates the number of amino acid substitutions per site. (B) Structural domains of GmNHX5. The colored section indicates Na+/H+ exchanger functional domain, Braces indicate transmembrane (TM) domains. (C) Relative expression of GmNHX5 under salt stress. RT-qPCR was used to detect the expression of GmNHX5 in soybean Jidou-7 and Mustang at 200 mM NaCl for the corresponding time. Actin was used as a reference gene. Each value represents the mean ± SE of three biological replicates relative to 0 h. Asterisks above the bars indicate significant differences from the values measured in Mustang (P < 0.05 by Student’s t-test).
We further examined the association between the transcription level of GmNHX5 and salt tolerance. The abundance of GmNHX5 mRNA in the salt-tolerant soybean cultivar Jidou-7 increased upon salt treatment and reached its peak at 6 h after treatment, then slightly decrease afterward, while the expression in the salt-sensitive cultivar Mustang showed no significant changes (Figure 1C), indicating that GmNHX5 may be involved in salt tolerance in soybean.
GmNHX5 Actively Transcribes in New Leaves and Vascular Tissues Under Salt Stress
We used β-glucuronidase (GUS) as a reporter to assess GmNHX5 promoter activity in various tissues in Arabidopsis. The GUS activity driven by the GmNHX5 promoter was detected in new leaves at the stem apex, in the veins of mature leaves, and the vascular tissues of the stems and roots (Figure 2A). Upon salt treatment, the GUS activity in the new leaves increased, while decreased in the veins of the mature leaves and the vascular tissues of the stem, and the GUS activity in the vascular tissues of the root increased with the salt treatment. These results indicate that the transcription of GmNHX5 is regulated by salt stress, and the GmNHX5 gene is actively transcribed in new leaves and roots under salt stress. Given that the root system is the tissue where the plant directly contacts the soil solution, we speculate that the transcriptional up-regulation of GmNHX5 in root vascular tissues and new leaves may lead to the up-regulation of GmNHX5 protein abundance in these tissues, thereby increases the adaptability to salt stress in roots, and protects young tissues from osmotic stress, respectively.
Figure 2. Detection of histological distribution of subcellular localization of GmNHX5. (A) Histochemical analysis of GmNHX5 expression. GUS activity of GmNHX5p:GUS in response to salt stress in transformed Arabidopsis was detected by GUS staining. Five-day-old seedlings were treated with 0, 100, or 200 mM NaCl for 24 h before staining. (B) Fluorescence of GmNHX5-GFP or GFP in Arabidopsis roots. (C) Fluorescence of GmNHX5-GFP and its colocalization with fluorochrome Golgi-Tracker Red. Confocal imaging of the roots of 5-day-old Arabidopsis seedlings transformed with plasmids carrying CaMV 35S:GmNHX5-GFP or CaMV 35S:GFP expression loci. Representative images of three biological replicates are shown. Bars = 200 μm in panel (A), Bars = 10 μm in panels (B,C). (D) Line profile was used to illustrate colocalization between GmNHX5-GFP and Golgi-Tracker Red along the dotted line in panel (C). Green and red lines indicate GmNHX5-GFP and Golgi-Tracker Red fluorescence profiles, respectively.
GmNHX5 Localizes to the Golgi Apparatus
The subcellular localization of protein determines the possible way of its function. To reveal the subcellular localization of GmNHX5, we examined the roots of Arabidopsis expressing the GmNHX5-GFP fusion protein. Using the plants transformed with the plasmid that harbors GFP alone as a control, a granular-distributed fluorescence pattern was detected in cells expressing the GmNHX5-GFP fusion protein (Figure 2B). In order to confirm the possibility of GmNHX5 distribution in the Golgi apparatus, we introduced the fluorochrome Golgi-Tracker Red to label the cellular Golgi apparatus. Granular distribution of fluorescence was also observed in cells labeled with Golgi-Tracker Red, and could co-localize with most of the fluorescence of GmNHX5-GFP (Figure 2C), indicating that GmNHX5 was distributed in the Golgi apparatus. The co-localization efficiency was further demonstrated by line profile (Figure 2D). These results indicate that GmNHX5 was distributed in the Golgi apparatus intracellularly.
Silencing of GmNHX5 Reduces Salt Tolerance in Soybean
As a preliminary exploration of the potential function of GmNHX5 involved in plant salt tolerance, we employed TRV-mediated VIGS technology and observed the tolerance of silent plants to salt stress. TRV carrying a specific segment of GmNHX5 mRNA produced GmNHX5 gene silencing with an efficiency of over 50% in TRV-inoculated Jidou-7 plants (Figure 3A). The silencing of GmNHX5 did not significantly affect plant development (Figure 3B). After 200 mM salt treatment, GmNHX5-silencing soybean plants showed wilting at the whole plant level and obstructed apical development, simple leaves showed symptoms of chlorosis and wilting, ternately compound leaves showed obvious growth retardation, while the growth status of the control (TRV: 00) plants was significantly better, and less affected by salt stress (Figure 3B). These results indicate that the silencing of GmNHX5 reduces the adaptability of soybean to salt stress.
Figure 3. Response of GmNHX5-silencing soybean plants to salt stress. (A) Silencing efficiency of GmNHX5 in TRV-inoculated plants. The relative expression of GmNHX5 was detected using RT-qPCR, and normalized to the values measured in TRV: 00-#1. Actin was used as a reference gene. Each value represents the mean ± SE of three technical replicates. Different letters above the bars indicate significant differences (P < 0.05 by one-way ANOVA and Duncan’s multiple comparison). (B) Phenotype of GmNHX5-silencing plants under salt stress condition. TRV-inoculated soybean plants were subjected to 0- or 200 mM NaCl solution treatment for 5 days. The images on the right show the phenotype of simple (top) and ternately compound leaves (bottom) of indicated plants under salt stress.
GmNHX5 Positively Regulates Salt Tolerance in Soybean Hairy Roots
We next analyzed the function of GmNHX5 on salt tolerance in hairy roots induced by soybean cotyledons. The plasmid carrying the GmNHX5 open reading frame (ORF) driven by the CaMV 35S promoter was used for overexpression (OE), and the CRISPR/Cas9 technique was used to conduct GmNHX5 targeted knockout (KO). We specifically designed sgRNA-1 and sgRNA-2 that target two GmNHX5-specific sites to increase knockout efficiency (Figure 4A). After sequencing, we found that GmNHX5 in soybean hairy roots transformed with the plasmid expressing Cas9 and sgRNAs had frameshift mutations (Figure 4B), which will result in the malfunctional GmNHX5 to be expressed. Meanwhile, the relative expression of GmNHX5 in the OE hairy roots was significantly higher than that in the hairy roots transformed with the empty vector (EV) (Figure 4C).
Figure 4. Functional characterization of GmNHX5 in soybean hair roots. (A) Schematic diagram of the of sgRNA binding sites used for CRISPR/Cas9-mediated knockout of GmNHX5 on the genome. The start codon (ATG), Exons (blue box), stop codon (TGA), sequences of Exon 2 and Exon 3 and their flanking sequences are indicated, sgRNAs and PAMs are underlined. (B) Sequence alignment of the knockout (KO) hairy roots at the sgRNA1 and sgRNA2 positions. EV, hairy roots transformed with empty vector pBSE401. (C) Expression of GmNHX5 in GmNHX5-overexpressing (OE) hairy roots. The relative expression of GmNHX5 was detected using RT-qPCR, and normalized to the values measured in the hair roots transformed with empty vector pCAMBIA3301 (EV). Actin was used as a reference gene. Each value represents the mean ± SE of three technical replicates of one hairy root induced by a single cotyledon. Asterisks above the bars indicate significant differences from the values measured in EV (P < 0.05 by Student’s t-test). (D) Phenotype of GmNHX5-overexpressing and knockout hairy roots under salt stress. Hairy roots were cultured on 1/4 MS medium with 0, 70, and 150 mM NaCl for 15 days before observation. WT, hairy roots induced by Agrobacterium rhizogenes without plasmid. Representative images of three biological replicates are shown. Maximum length (E) and fresh weight (F) of soybean hairy roots in (D). The data in panels (E,F) shows the mean ± SD of 30 independent measurements. Different letters indicate significant difference (P < 0.05 by one-way ANOVA/Duncan).
The growth status and maximum root length of OE and KO hairy roots were similar on medium free of NaCl (Figures 4D,E), whereas the fresh weight of the OE and the wild-type (WT) were significantly higher than the KO (Figure 4F). The difference in fresh weight (Figure 4F) in 0 mM NaCl treatment between the WT-, OE-, and KO hairy roots suggest a possible K+ transport role of GmNHX5 in soybean, and the reduction in fresh weight in GmNHX5-KO hairy roots may be due to impaired K+-uptake and utilization where GmNHX5 was knocked out. On the medium containing 70 mM NaCl, the growth status of the OE hairy roots was superior to the others, the maximum root length and fresh weight of the OE were significantly higher than that of the KO and WT hairy roots. When the NaCl concentration increased to 150 mM, the maximum root length of the OE was not significantly different from the WT, but the fresh weight was significantly higher than that of the WT and the KO. These results indicate that GmNHX5 positively regulates salt tolerance in soybean hairy roots.
GmNHX5 Enhances the Adaptability to Salt Stress in Transformed Soybean
The fact that GmNHX5 positively regulates salt tolerance in hairy roots suggests that the gene may functions similarly at the organism level. To test this hypothesis, we generated GmNHX5-overexpressing plants in the salt-sensitive variety Mustang. Colloidal gold test on the phosphinothricin acetyl-transferase (PAT), a selective marker on the plasmid, showed that the protein was detected in #3 and #4, two transformed T3 soybean lines (Figure 5A). Southern blot analysis showed that the plasmid T-DNA carrying the CaMV 35S: GmNHX5 expression cassette was integrated into the soybean genome in a single copy (Figure 5B). The relative transcription level of GmNHX5 in #3 and #4 was significantly higher than that in the WT (Figure 5C). Under conventional conditions, the growth of GmNHX5-OE plants was not significantly different from that of WT. After treatment with 150 mM NaCl for 20 days, WT plants showed strong wilting, chlorosis, and growth obstacles, whereas the condition of GmNHX5-OE plants was significantly better, with a merely mild wilting phenotype (Figure 5D). Furthermore, biomass analysis showed that under conventional conditions, the leaf dry weight and leaf fresh weight of GmNHX5-OE plants were not significantly different from those of WT plants (Figures 5E,F); under salt stress, the leaf dry weight and leaf fresh weight of GmNHX5-OE plants were markedly higher than those of WT plants, which was consistent with the observed phenotype.
Figure 5. Overexpression of GmNHX5 enhances salt tolerance in soybean. (A) Identification of GmNHX5-overexpressing soybean plants. Leaf extract from independent GmNHX5-OE transformed T3 lines #3 and #4 were tested using colloidal gold strips. Phosphinothricin acetyl-transferase (PAT) refers to the herbicide resistance marker on the plasmid. –, untransformed soybean; +, positive control. (B) Southern blot analysis of the copy number of plasmids in the genome of GmNHX5-overexpressing plants. +, plasmids for soybean transformation. Numbers indicate molecular weight standards. (C) Expression of GmNHX5 in GmNHX5-overexpressing (OE) soybean plants. The relative expression of GmNHX5 was detected using RT-qPCR, and normalized to the values measured in the untransformed soybean (WT). Actin was used as a reference gene. Each value represents the mean ± SE of three biological replicates. Asterisks above the bars indicate significant differences from the values measured in WT (P < 0.05 by Student’s t-test). (D) Phenotypes of soybean plants overexpressing GmNHX5 under salt stress. Three-weeks-old WT and GmNHX5-overexpressing soybean seedlings were irrigated with Hoagland’s nutrient solution with 0 or 150 mM NaCl for 20 days. Arrows indicate partial close-ups of the plant. Leaf dry weight (E) and leaf fresh weight (F) were measured as the total leaf weight in independent GmNHX5-OE transformed T3 lines #3 and #4. The data in panels (E,F) shows the mean ± SD of triplicate measurements. Different letters indicate significant difference (P < 0.05 by one-way ANOVA/Duncan).
An important mechanism by which plants respond to and adapt to salt stress is the accumulation of osmotic adjustment substances, such as free proline and betaine. Under normal conditions, the contents of free proline (Figure 6A) and betaine (Figure 6B) in GmNHX5-OE plants were not significantly different from those in WT plants. However, when exposed to salt stress, the increase in the content of these two osmotic adjustment substances in GmNHX5-OE plants was much greater than that in WT plants. Meanwhile, GmNHX5-OE plants accumulated less MDA than WT plants under salt stress (Figure 6C), indicating that GmNHX5 attenuated lipid peroxidation caused by salt stress-induced ROS.
Figure 6. Physiological response of GmNHX5-overexpressing soybean to salt stress. Content of free proline (A), betaine (B), and MDA (C) were measured in indicated GmNHX5-overexpressing T3 lines 6-days after 0 or 150 mM NaCl treatment. The data shows the mean ± SD of triplicate measurements. Different letters indicate significant difference (P < 0.05 by one-way ANOVA/Duncan).
GmNHX5-Overexpressing Soybean Has Higher K+/Na+ Ratio Under Salt Stress
In addition to accumulating osmotic adjustment substances, the transportation and compartmentalization of K+ and Na+ at the cellular and histological levels are also one of the main manners for plants to detoxify salt stress (Qiu, 2016).
In order to inquire into the mechanism and characteristics of GmNHX5-mediated salt tolerance in plants, we examined the content of K+ and Na+ in #3 and #4, two independent GmNHX5-OE T3 lines. Under normal conditions, the content of Na+ in OE leaves was lower than that of WT leaves (Figure 7A); the Na+ content in OE roots was also lower than that in WT plants (Figure 7B), while the overexpression of GmNHX5 did not affect the Na+ content in stems (Figure 7C). In addition, under normal conditions, the K+ content of OE plants was higher than that of WT plants in all tissues tested (Figures 7D–F). When exposed to salt stress, Na+ content in all tissues detected in OE plants was significantly lower than that in WT plants, including leaves (Figure 7A), roots (Figure 7B), and stems (Figure 7C). The content of K+ in leaves (Figure 7D) and roots (Figure 7E) of OE plants were significantly higher than that of WT plants, and the content of K+ in stems of OE plants was lower than that of WT plants except (Figure 7F).
Figure 7. Na+, K+, and K+/Na+ ratio in soybean overexpressing GmNHX5. Content of Na+ (A–C) and K+ (D–F) was measured in #3 and #4, two independent GmNHX5-overexpressing soybean T3 lines 6 h after 0 or 150 mM NaCl treatment. Leaves, roots, and stems were examined to show the histologic distribution of Na+ and K+. K+/Na+ ratio was separately calculated in leaves (G), roots (H), and stems (I). The data shows the mean ± SD of triplicate measurements. Different letters indicate significant difference (P < 0.05 by one-way ANOVA/Duncan).
As a result, under normal conditions, the K+/Na+ ratio of GmNHX5-OE plants was higher than that of WT plants (Figures 7G–I). Under salt stress, the K+/Na+ ratio in leaves (Figure 7G), roots (Figure 7H), and stems (Figure 7I) of OE plants was also higher than that of WT plants. The higher K+/Na+ ratio of OE plants under salt stress may be the reason why they were more accommodating to salt stress.
GmNHX5 Regulates Salt-Responsive Genes
Our results suggested that GmNHX5 localized in the Golgi apparatus and regulated the accumulation of Na+ and K+, which implies that it may indirectly affect the uptake of these two ions. Transcriptional regulation of salt-related genes is a vital way for plants to strengthen their adaptability under salt stress. In order to reveal the regulatory effect of GmNHX5 on salt-responsive genes, we examined the transcriptional expression of several salt-responsive genes (Figures 8A–F), including K+ channel GmSKOR (Figure 8A), GmHKT1 (Figure 8C), and GmAKT1 (Figure 8D); K+ pump GmHAK5 (Figure 8E); and exchangers GmSOS1 (Figure 8B) and GmNHX1 (Figure 8F) in soybean overexpressing GmNHX5. It was found that under normal conditions, overexpression of GmNHX5 induced the expression of GmSKOR and GmHKT1; reduced the expression of GmNHX1; but did not affect the expression of GmSOS1, GmAKT1, and GmHAK5. Under salt stress, however, the overexpression of GmNHX5 increased the expression of GmSKOR, GmSOS1, and GmHKT1; decreased the expression of GmHAK5 and GmNHX1; and did not affect the expression of GmAKT1. These results indicate that the overexpression of GmNHX5 affects the expression level of salt-related genes, but the specific molecular mechanism by which GmNHX5 regulates the expression of these genes remains to be further studied.
Figure 8. Expression analysis of salt tolerance related genes. The relative transcription level of GmSKOR (A), GmSOS1 (B), GmHKT1 (C), GmAKT1 (D), GmHAK5 (E), and GmNHX1 (F) were detected in the GmNHX5-overexpressing T3 soybean seedlings 6 h after 0 or 150 mM NaCl treatment, and normalized to the values measured in WT exposed to 0 mM NaCl. Actin was used as a reference gene. Each value represents the mean ± SE of three biological replicates. Different letters indicate significant difference (P < 0.05 by one-way ANOVA/Duncan).
Discussion
Given that NHXs are ion/proton exchangers that directly regulate cellular Na+ or K+ uptake or compartmentalization, such genes have been shown to possess the capability to improve the salt tolerance in various plant species, including Medicago sativa (Zhang et al., 2015), Capsicum annuum (Bulle et al., 2016), Cucumis sativus (Wang et al., 2013), etc. The key to NHX-mediated increase in plant salt tolerance is to increase the ratio of cytoplasmic K+/Na+ to avoid excessive accumulation of these two ions, so as to ultimately reduce the damage to organelles. Although many studies have shown that members of the NHX family of plants have the ability to exchange Na+ (Shi et al., 2002), there is also evidence that they may function by transporting K+ (Bassil et al., 2011b; Wang L. et al., 2015). In view of the large number of ion channels and ion pumps in the membrane system of plant cells, there is currently no direct evidence of the exact ion molecule to be transported by NHXs. Although there are still in debate over the actual ions used by NHXs to exchange H+, it is certain that most of the NHXs has been discovered reduces the Na+ content in cytoplasm or maintain high K+/Na+ ratio.
According to the structural characteristics, NHXs are mainly classified into two categories, AtNHX1, AtNHX2, AtNHX3, AtNHX4, AtNHX7/SOS1, and AtNHX8 are the first type of NHX, which are located in the plasma membrane and vacuolar membrane. The main functional domain of these NHXs belongs to the “Na+-H+_Exchanger” superfamily and was considered to have an important role in maintaining the Na+ or K+ (Bassil et al., 2011b) levels, Na+/K+ homeostasis and pH stability (Dibrov and Fliegel, 1998; Rodríguez-Rosales et al., 2009). AtNHX5 and AtNHX6 belong to the second type of NHX, which is located on the vesicle membrane in the endoplasmic reticulum and the Golgi network, whose main domain belongs to the “MFS” superfamily, and was found to regulate pH and affects the sorting behavior of the Golgi network (Bassil et al., 2011a; Reguera et al., 2015; Date et al., 2016). According to phylogenetic analysis, GmNHX5 was highly similar to NHXs localized by the Golgi network/endosomal system represented by AtNHX5 and AtNHX6 (Figure 1A). The results of subcellular localization further indicate that GmNHX5 was localized on the Golgi network (Figure 2C), in accordance with the distribution pattern of AtNHX5 (Bassil et al., 2011a). The Golgi network is considered to be a center for cell component transportation, which connects to the endoplasmic reticulum, vacuole, and plasma membrane (Larkin et al., 2016). Under salt stress conditions, plants produce large amounts of secondary metabolites and discharge or isolate salt ions that destroy intracellular structures, which all depend on the role of the Golgi network (Ueda et al., 2016). Thereby, we speculate that the salt tolerance-related molecular mechanism of GmNHX5 may rely on the function of the Golgi network. In addition to the distribution within the cell, the histological distribution of the protein also plays a decisive role in its function. GUS activity driven by OsNHX1 promoter was localized to the guard cells and trichome, whereas OsNHX5 promoter-GUS activity was localized to the root tip and pollen grains (Fukuda et al., 2011). In this study, it was found that the GmNHX5 promoter activity was distributed in the young leaves and roots in the transformed Arabidopsis, and with the increase of the intensity of salt stress, the distribution of GUS signal enriched in the young leaves and the growing point, and the vascular bundles in the roots. Given that the Golgi apparatus connects vacuole and plasma membrane, the localization of GmNHX5 on the Golgi apparatus (Figure 2C) of the vascular (Figure 2A) indicating that it may participate in the transport of K+. Since the K+ content in the root system is higher than that in the stem and leaves under normal condition (Figures 7D–F), we speculate that GmNHX5 may participate in the transfer of K+ absorbed by the root to the aboveground organs to protect the young tissues, and at the same time regulate the utilization of K+.
TRV-VIGS-induced GmNHX5 silencing plants showed enhanced sensitivity to salt stress, which gave us a glimpse of the functions of GmNHX5 related to salt tolerance and prompted us to further explore the functions of GmNHX5 in diverse ways. The hairy roots of plants induced by Agrobacterium rhizogenes inherit the plasmid carried in the bacteria, and are considered to be an ideal material for genetic function verification at the histological level. Due to its fast and convenient characteristics, hairy roots have been used in the study of physiological mechanisms in various plant species, including Solanum lycopersicum (Ron et al., 2014), Saussurea involucrate (Li et al., 2006), and Glycine max (Kim et al., 2009). In this study, GmNHX5 was overexpressed or knocked out with CRISPR/Cas9 in soybean hairy roots, and it was found that hair roots overexpressing GmNHX5 had lower salt sensitivity, which was opposite to that of GmNHX5 was knocked out (Figure 4). In addition, in view of the significant increase of K+ content in the roots of GmNHX5-OE plants (Figure 7E), the decrease in the biomass of GmNHX5-KO hairy roots may also be related to the impaired of K+ uptake. These results are consistent with the function of AtNHX5 in Arabidopsis (Zhu et al., 2018) and prove that the expression of GmNHX5 was positively correlated with the salt tolerance of soybean. In addition to the hairy root experiments, this study also generated the T3 generation transformed lines overexpressing GmNHX5 on the soybean Mustang, and verified the transformation event from the protein level of the reporter gene PAT (bar), the relative transcription level of GmNHX5, and the presence of plasmid integration in the soybean genome (Figure 5). We also examined the influence of GmNHX5 on the physiological state of soybean plants from the contents of proline, betaine, and MDA (Figure 6). As well-known osmotic adjustment substances in plants, the content of free proline and betaine increases significantly when undergoing osmotic stress. Varieties with higher drought and salt tolerance tend to accumulate more free proline, it can be considered that the free proline content in plants reflects the osmotic tolerance to some extent (Trinchant et al., 2004; Nedjimi, 2013; Al-Farsi et al., 2020). In this study, we found that GmNHX5-overexpressing soybean plants accumulated more free proline and betaine when subjected to salt stress, but the contents of these two osmotic adjustment substances were not significantly different from WT under normal conditions, indicating that GmNHX5 does not directly increase the background levels of proline and betaine, but may regulate their accumulation in a salt stress-dependent manner. MDA is a product of oxidative damage to the cell membrane, and its accumulation level indicates the degree of cell damage (Wang, 2014). The MDA content of GmNHX5-overexpressing plants under salt stress was significantly lower than that in WT plants, indicating that GmNHX5 may protect cells from oxidative damage.
Na+, K+, and Cl– are the main substances that affect the osmotic pressure of plant cells. Whether in monocotyledonous or dicotyledonous plants, the presence of these ions constitutes 80–95% of the osmotic pressure (Hastings and Gutknecht, 1974; Adams et al., 1992). Our results suggest that GmNHX5 positively regulates soybean tolerance to salt stress, we speculate that this may be achieved by multiple ways: On the one hand, GmNHX5 may regulate the increase of K+/Na+ ratio (Figures 7G–I); On the other hand, it may also increase the content of osmotic adjustment substances. In view of the similarity of the electrochemical characteristics of Na+ and K+, we conjecture that the transport activity of GmNHX5 is suitable for both of these two ions. As a homolog of AtNHX6 which participates in the maturation process of the pre-vacuolar compartment and uptakes K+ into the Golgi network (Bassil et al., 2011a; Wang L. et al., 2015; Zhu et al., 2018), we speculate that GmNHX5 may function in the following ways: (1) Transporting Na+ to the trans-Golgi network and pre-vacuolar compartment to carry out the secretion of Na+ directly; (2) Transporting K+ in the cytoplasm to the Golgi apparatus, leading to acidification of the cytoplasm and activation of intracellular signals such as cytosolic pH. These signals may further act on transcription factors in the nucleus and regulate the expression of salt-responsive genes such as GmSKOR, GmSOS1, GmHKT1, and GmNHX1, and induce the contents of free proline and betaine, to combat ROS burst and organelle damage caused by salt stress. Since ROS induces the transcription of HAK5 (Huang et al., 2019), the reduced GmHAK5 expression in the GmNHX5-OE plants (Figure 8E) may be the mobilization of ROS-scavenging system. The heterologous expression of GmNHX1 can positively regulate the Na+ flow rate of Arabidopsis roots under salt stress, and enhance salt tolerance (Sun et al., 2019). Our results showed that the overexpression of GmNHX5 constitutively down-regulated the expression of GmNHX1 (Figure 8F). It is speculated that the regulation of cytoplasmic K+/Na+ balance by GmNHX5 reduces the requirement for the expression of the vacuole-localized GmNHX1, and GmNHX5 may play upstream of GmNHX1 in conferring salt tolerance. The induced expression of GmSOS1 (Figure 8B) in GmNHX5-OE plants may also be one of the reasons for the reduced intracellular Na+. The expression of AKT1 gene was decreased in the salt-tolerant cultivar of Fragaria variety under salt stress condition (Garriga et al., 2017), but we found that the expression of GmAKT1 was unchanged by salt stress, which may be due to the differences in species and salt tolerance mechanisms. The specific molecular mechanism of GmNHX5 in regulating the expression of these salt-responsive gene remains to be further studied. These results not only deepen our understanding of the role and molecular mechanism of NHXs located in the Golgi apparatus in response to salt stress, but also provide genetic resources for the cultivation of salt-tolerant plant materials.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author Contributions
DW, JC, and JZ supervised the project. DW and JZ conceived and designed the experiments. TS, NM, and CW performed the experiments and analyzed the data. HF and MW contributed reagents, materials, and analysis tools. TS and NM wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
The sources of funding for this study are provided by China Ministry of Agriculture, Hebei Provincial Department of Science and Technology, and Department of Human Resources and Social Security of Hebei Province through the following projects:
(1) Open project of Key Laboratory of crop salt tolerance evaluation and genetic improvement in Hebei Province (2020-006), JZ;
(2) Major scientific and technological projects for breeding new varieties of genetically modified organisms (2009ZX08004-001B), DW;
(3) Major scientific and technological projects for breeding new varieties of genetically modified organisms (2014ZX0800402B-001), DW;
(4) The project for the introducted overseas students funded by Hebei Province (CN201706), JZ;
(5) Hebei Natural Science Foundation (C2020301020), JZ.
The funding bodies had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Conflict of Interest
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
We thank Zhouliang Yu of the University of California, Berkeley, for his help in polishing the language of this article.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.638340/full#supplementary-material
Footnotes
References
Adams, P., Thomas, J. C., Vernon, D. M., Bohnert, H. J., and Jensen, R. G. (1992). Distinct cellular and organismic responses to salt stress. Plant Cell Physiol. 33, 1215–1223. doi: 10.1093/oxfordjournals.pcp.a078376
Al-Farsi, S., Nawaz, A., Anees ur, R., Nadaf, S., Al-Sadi, A., Siddique, K., et al. (2020). Effects, tolerance mechanisms and management of salt stress in lucerne (Medicago sativa). Crop Pasture Sci. 71, 411–428. doi: 10.1071/CP20033
An, J., Hu, Z., Che, B., Chen, H., Yu, B.-J., and Cai, W. (2017). Heterologous expression of Panax ginseng PgTIP1 confers enhanced salt tolerance of soybean cotyledon hairy roots, composite, and whole plants. Frontiers in Plant Science 8, 1232. doi: 10.3389/fpls.2017.01232
Aviv, A. (1993). “The Na+/H+ antiport: its role in cellular metabolism and in the pathophysiology of essential hypertension,” in Ionic Transport in Hypertension: New Perspectives, eds A. Coca and R. Garay (Boca Raton, FL: CRC Press), 119–132. Available online at: https://books.google.fr/books?id=zHfXqrFhr1AC&dq=The+Na%2B/H%2B+Antiport:+Its+Role+in+Cellular+Metabolism+and+in+The+Pathophysiology+of+Essential+Hypertension&lr=&hl=zh-CN&source=gbs_navlinks_s
Barragán, V., Leidi, E. O., Andrés, Z., Rubio, L., De Luca, A., Fernández, J. A., et al. (2012). Ion exchangers NHX1 and NHX2 mediate active potassium uptake into vacuoles to regulate cell turgor and stomatal function in Arabidopsis. Plant Cell 24, 1127–1142. doi: 10.1105/tpc.111.095273
Bassil, E., Ohto, M.-A., Esumi, T., Tajima, H., Zhu, Z., Cagnac, O., et al. (2011a). The Arabidopsis intracellular Na+/H+ antiporters NHX5 and NHX6 are endosome associated and necessary for plant growth and development. Plant Cell 23, 224–239. doi: 10.1105/tpc.110.079426
Bassil, E., Tajima, H., Liang, Y. C., Ohto, M. A., Ushijima, K., Nakano, R., et al. (2011b). The Arabidopsis Na+/H+ antiporters NHX1 and NHX2 control vacuolar pH and K+ homeostasis to regulate growth, flower development, and reproduction. Plant Cell 23, 3482–3497. doi: 10.1105/tpc.111.089581
Bates, L. S., Waldren, R. P., and Teare, I. D. (1973). Rapid determination of free proline for water-stress studies. Plant Soil 39, 205–207. doi: 10.1007/BF00018060
Bohnert, H., and Jensen, R. (1996). Strategies for engineering water-stress tolerance in plants. Trends Biotechnol. 14, 89–97. doi: 10.1016/0167-7799(96)80929-2
Boudsocq, M., and Sheen, J. (2012). CDPKs in immune and stress signaling. Trends Plant Sci. 18, 30–40. doi: 10.1016/j.tplants.2012.08.008
Bulle, M., Yarra, R., and Abbagani, S. (2016). Enhanced salinity stress tolerance in transgenic chilli pepper (Capsicum annuum L.) plants overexpressing the wheat antiporter (TaNHX2) gene. Mol. Breed. 36:36. doi: 10.1007/s11032-016-0451-5
Chen, L., Jiang, B., Wu, C., Sun, S., Hou, W., and Han, T. (2014). GmPRP2 promoter drives root-preferential expression in transgenic Arabidopsis and soybean hairy roots. BMC Plant Biol. 14:245. doi: 10.1186/s12870-014-0245-z
Date, S., Chen, C.-Y., Chen, Y., and Jansen, M. (2016). Experimentally optimized threading structures of the proton-coupled folate transporter. FEBS Open Bio 6, 216–230. doi: 10.1002/2211-5463.12041
Dibrov, P., and Fliegel, L. (1998). Comparative molecular analysis of Na+/H+ exchangers: a unified model for Na+/H+ antiport? FEBS Lett. 424, 1–5. doi: 10.1016/s0014-5793(98)00119-7
Esfandiari, E., and Gohari, G. (2017). Response of ROS-scavenging systems to salinity stress in two different wheat (Triticum aestivum L.) Cultivars. Not. Bot. Horti Agrobot. Cluj Napoca 45, 287–291. doi: 10.15835/nbha45110682
Fukuda, A., Nakamura, A., Hara, N., Toki, S., and Tanaka, Y. (2011). Molecular and functional analyses of rice NHX-type Na+/H+ antiporter genes. Planta 233, 175–188. doi: 10.1007/s00425-010-1289-4
Gadsby, D. C. (2009). Ion channels versus ion pumps: the principal difference, in principle. Nat. Rev. Mol. Cell Biol. 10, 344–352. doi: 10.1038/nrm2668
Gaonkar, V., and Rosentrater, K. A. (2019). “Soybean,” in Integrated Processing Technologies for Food and Agricultural by-Products, eds Z. Pan, R. Zhang, and S. Zicari (San Diego, CA: Academic Press), 73–104.
Garriga, M., Raddatz, N., Véry, A. A., Sentenac, H., Rubio-Meléndez, M. E., González, W., et al. (2017). Cloning and functional characterization of HKT1 and AKT1 genes of Fragaria spp.-Relationship to plant response to salt stress. J. Plant Physiol. 210, 9–17. doi: 10.1016/j.jplph.2016.12.007
Gaudet, P., Livstone, M. S., Lewis, S. E., and Thomas, P. D. (2011). Phylogenetic-based propagation of functional annotations within the Gene Ontology consortium. Brief. Bioinform. 12, 449–462. doi: 10.1093/bib/bbr042
Geilfus, C.-M. (2019). Controlled Environment Horticulture. Springer: Berlin, 69–80. Available online at: https://books.google.fr/books?id=9o-qDwAAQBAJ&dq=Controlled+environment+horticulture+Geilfus,+Christoph-Martin&lr=&hl=zh-CN&source=gbs_navlinks_s
Hasanuzzaman, M., Nahar, K., Rahman, A., Mahmud, J. A., Hossain, M. S., and Fujita, M. (2016). “Soybean production and environmental stresses,” in Environmental Stresses in Soybean Production, ed. M. Miransari (Cambridge, MA: Academic Press), 61–102. doi: 10.1016/b978-0-12-801535-3.00004-8
Hastings, D. F., and Gutknecht, J. (1974). “Turgor pressure regulation: modulation of active potassium transport by hydrostatic pressure gradients,” in Membrane Transport in Plants, eds U. Zimmermann and J. Dainty (Berlin: Springer), 79–83. doi: 10.1007/978-3-642-65986-7_10
Horie, T., Costa, A., Kim, T.-H., Han, M., Horie, R., Leung, H.-Y., et al. (2007). Rice OsHKT2;1 transporter mediates large Na+ influx component into K+-starved roots for growth. EMBO J. 26, 3003–3014. doi: 10.1038/sj.emboj.7601732
Huang, Y., Cao, H., Yang, L., Chen, C., Shabala, L., Xiong, M., et al. (2019). Tissue-specific respiratory burst oxidase homolog-dependent H2O2 signaling to the plasma membrane H+-ATPase confers potassium uptake and salinity tolerance in Cucurbitaceae. J. Exp. Bot. 70, 5879–5893. doi: 10.1093/jxb/erz328
James, R., Blake, C., Byrt, C., and Munns, R. (2011). Major genes for Na+ exclusion, Nax1 and Nax2 (wheat HKT1;4 and HKT1;5), decrease Na+ accumulation in bread wheat leaves under saline and waterlogged conditions. J. Exp. Bot. 62, 2939–2947. doi: 10.1093/jxb/err003
Jiang, Y.-Q., and Deyholos, M. (2008). Function characterization of Arabidopsis NaCl-inducible WRKY25 and WRKY33 transcription factors in abiotic stresses. Plant Mol. Biol. 69, 91–105. doi: 10.1007/s11103-008-9408-3
Jiang, Y.-Q., Yang, B., and Deyholos, M. (2009). Functional characterization of the Arabidopsis bHLH92 transcription factor in abiotic stress. Mol. Genet. Genomics 282, 503–516. doi: 10.1007/s00438-009-0481-3
Jones, D., Taylor, W., and Thornton, J. (1992). The Rapid generation of mutation data matrices from protein sequences. Comput. Appl. Biosci. 8, 275–282. doi: 10.1093/bioinformatics/8.3.275
Kim, S.-H., Lim, C.-W., Park, J.-Y., and Hwang, C.-H. (2009). Inhibition of SKTI synthesis in Agrobacterium rhizogenes-induced hairy root reduces the number of nodule in soybean. Korean J. Crop Sci. 54, 299–306. doi: 10.1097/00000542-199612000-00051
Kumar, S., Stecher, G., Li, M., Knyaz, C., and Tamura, K. (2018). MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35, 1547–1549. doi: 10.1093/molbev/msy096
Larkin, H., Costantino, S., Seaman, M., and Lavoie, C. (2016). Calnuc function in endosomal sorting of lysosomal receptors. Traffic 17, 416–432. doi: 10.1111/tra.12374
Lee, H., Park, S.-Y., and Zhang, Z. (2013). “An overview of genetic transformation of soybean,” in A Comprehensive Survey of International Soybean Research - Genetics, Physiology, Agronomy and Nitrogen Relationships, ed. J. E. Board (Rijeka: Intech press).
Li, F.-X., Jin, Z.-P., Zhao, D.-X., Cheng, L.-Q., Fu, C.-X., and Ma, F. (2006). Overexpression of the Saussurea medusa chalcone isomerase gene in S. involucrata hairy root cultures enhances their biosynthesis of apigenin. Phytochemistry 67, 553–560. doi: 10.1016/j.phytochem.2005.12.004
Li, M., Li, Y., Li, H.-Q., and Wu, G. (2011). Overexpression of AtNHX5 improves tolerance to both salt and drought stress in Broussonetia papyrifera (L.) Vent. Tree Physiol. 31, 349–357. doi: 10.1093/treephys/tpr003
Li, S., Lei, L., Somerville, C. R., and Gu, Y. (2012). Cellulose synthase interactive protein 1 (CSI1) links microtubules and cellulose synthase complexes. Proc. Natl. Acad. Sci. U.S.A. 109, 185–190. doi: 10.1073/pnas.1118560109
Ma, S., Niu, H., Liu, C., Zhang, J., Hou, C., and Wang, D. (2013). Expression stabilities of candidate reference genes for RT-qPCR under different stress conditions in soybean. PLoS One 8:e75271. doi: 10.1371/journal.pone.0075271
Maathuis, F., and Amtmann, A. (1999). K+ nutrition and Na+ toxicity: the basis of cellular K+/Na+ ratios. Ann. Bot. 84, 123–133. doi: 10.1006/anbo.1999.0912
Massa, D., and Melito, S. (2019). “Signaling molecules in ecophysiological response mechanisms of salt-stressed plants,” in Plant Signaling Molecules: Role and Regulation under Stressful Environments, eds M. I. R. Khan, P. S. Reddy, A. Ferrante and N. A. Khan (Sawston: Woodhead Publishing), 1–18. doi: 10.1016/b978-0-12-816451-8.00001-0
Moradi, F., and Ismail, A. (2007). Responses of photosynthesis, chlorophyll fluorescence and ROS-scavenging systems to salt stress during seedling and reproductive stages in rice. Ann. Bot. 99, 1161–1173. doi: 10.1093/aob/mcm052
Muthappa, S.-K., and Mysore, K. (2014). Tobacco rattle virus-based virus-induced gene silencing (TRV-VIGS) in Nicotiana benthamiana. Nat. Protoc. 9, 1549–1562. doi: 10.1038/nprot.2014.092
Nedjimi, B. (2013). “Involvement of proline in plant response to salt stress,” in Proline: Biosynthesis, Regulation and Health Benefits, ed. B. Nedjimi (Hauppauge, NY: Nova Science Publishers, Inc), 1–10. doi: 10.1007/978-1-4614-6108-1_1
Paz, M., Shou, H., Guo, Z., Zhang, Z., Banerjee, A., and Wang, K. (2004). Assessment of conditions affecting Agrobacterium-mediated soybean transformation using the cotyledonary node explant. Euphytica 136, 167–179. doi: 10.1023/B:EUPH.0000030669.75809.dc
Pehlivan, N., Sun, L., Jarrett, P., Yang, X., Mishra, N., Chen, L., et al. (2016). Co-overexpressing a plasma membrane sodium/proton antiporter and a vacuolar membrane sodium/proton antiporter significantly improves salt tolerance in transgenic Arabidopsis plants. Plant Cell Physiol. 57, 1069–1084. doi: 10.1093/pcp/pcw055
Qi, W., Zhang, L., Xu, H., Wang, L., and Jiao, Z. (2014). Physiological and molecular characterization of the enhanced salt tolerance induced by low-dose gamma irradiation in Arabidopsis seedlings. Biochem. Biophys. Res. Commun. 450, 1010–1015. doi: 10.1016/j.bbrc.2014.06.086
Qiu, Q. (2016). Plant endosomal NHX antiporters: activity and function. Plant Signal. Behav. 11:e1147643. doi: 10.1080/15592324.2016.1147643
Qiu, Q., Guo, Y., Dietrich, M., Schumaker, K., and Zhu, J.-K. (2002). Regulation of SOS1, a plasma membrane Na+/H+ exchanger in Arabidopsis thaliana, by SOS2 and SOS3. Proc. Natl. Acad. Sci. U.S.A. 99, 8436–8441. doi: 10.1073/pnas.122224699
Reguera, M., Bassil, E., Tajima, H., Wimmer, M., Chanoca, A., Otegui, M. S., et al. (2015). pH Regulation by NHX-type antiporters is required for receptor-mediated protein trafficking to the vacuole in Arabidopsis. Plant Cell 27, 1200–1217. doi: 10.1105/tpc.114.135699
Rodríguez-Rosales, M. P., Gálvez, F., Huertas, R., Aranda, N., Baghour, M., Cagnac, O., et al. (2009). Plant NHX cation/proton antiporters. Plant Signal. Behav. 4, 265–276. doi: 10.4161/psb.4.4.7919
Ron, M., Kajala, K., Pauluzzi, G., Wang, D., Reynoso, M., Zumstein, K., et al. (2014). Hairy root transformation using Agrobacterium rhizogenes as a tool for exploring cell type-specific gene expression and function using tomato as a model. Plant Physiol. 166, 455–469. doi: 10.1104/pp.114.239392
Schmedes, A., and Ga, H. (1989). A new thiobarbituric acid (TBA) method for determining free malondialdehyde (MDA) and hydroperoxides selectively as a measure of lipid peroxidation. J. Oil Fat Ind. 66, 813–817. doi: 10.1007/BF02653674
Schmidt, R., Mieulet, D., Hubberten, H.-M., Obata, T., Hoefgen, R., Fernie, A., et al. (2013). SALT-RESPONSIVE ERF1 regulates reactive oxygen species-dependent signaling during the initial response to salt stress in rice. Plant Cell 25, 2115–2131. doi: 10.1105/tpc.113.113068
Schmutz, J., Cannon, S. B., Schlueter, J., Ma, J., Mitros, T., Nelson, W., et al. (2010). Genome sequence of the palaeopolyploid soybean. Nature 463, 178–183. doi: 10.1038/nature08670
Shi, H., Quintero, F., Pardo, J., and Zhu, J.-K. (2002). The putative plasma membrane Na+/H+ antiporter SOS1 controls long-distance Na+ transport in plants. Plant Cell 14, 465–477. doi: 10.1105/tpc.010371.et
Singleton, P., and Bohlool, B. (1984). Effect of salinity on nodule formation by soybean. Plant Physiol. 74, 72–76. doi: 10.1104/pp.74.1.72
Sun, T.-J., Fan, L., Yang, J., Cao, R.-Z., Yang, C.-Y., Zhang, J., et al. (2019). A Glycine max sodium/hydrogen exchanger enhances salt tolerance through maintaining higher Na+ efflux rate and K+/Na+ ratio in Arabidopsis. BMC Plant Biol. 19:469. doi: 10.1186/s12870-019-2084-4
Thompson, J. D., Higgins, D. G., and Gibson, T. J. (1994). CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22, 4673–4680. doi: 10.1093/nar/22.22.4673
Tran, L.-S., Nakashima, K., Sakuma, Y., Simpson, S., Fujita, Y., Maruyama, K., et al. (2004). Isolation and functional analysis of Arabidopsis stress-inducible NAC transcription factors that bind to a drought-responsive cis-element in the Early responsive to dehydration stress 1 promoter. Plant Cell 16, 2481–2498. doi: 10.1105/tpc.104.022699
Trinchant, J.-C., Boscari, A., Spennato, G., Van de Sype, G., and Le Rudulier, D. (2004). Proline betaine accumulation and metabolism in alfalfa plants under sodium chloride stress. Exploring its compartmentalization in nodules. Plant Physiol. 135, 1583–1594. doi: 10.1104/pp.103.037556
Ueda, M., Tsutsumi, N., and Fujimoto, M. (2016). Salt stress induces internalization of plasma membrane aquaporin into the vacuole in Arabidopsis thaliana. Biochem. Biophys. Res. Commun. 474, 742–746. doi: 10.1016/j.bbrc.2016.05.028
Wang, J., Qiu, N., Wang, P., Zhang, W., Yang, X., Chen, M., et al. (2019). Na+ compartmentation strategy of Chinese cabbage in response to salt stress. Plant Physiol. Biochem. 140, 151–157. doi: 10.1016/j.plaphy.2019.05.001
Wang, L., Wu, X., Liu, Y., and Qiu, Q. (2015). AtNHX5 and AtNHX6 control cellular K+ and pH homeostasis in Arabidopsis: three conserved acidic residues are essential for K+ transport. PLoS One 10:e0144716. doi: 10.1371/journal.pone.0144716
Wang, L.-F. (2014). Physiological and molecular responses to drought stress in rubber tree (Hevea brasiliensis Muell. Arg.). Plant Physiol. Biochem. 83, 243–249. doi: 10.1016/j.plaphy.2014.08.012
Wang, P., Guo, Q., Wang, Q., Zhou, X.-R., and Wang, S. (2015). PtAKT1 maintains selective absorption capacity for K+ over Na+ in halophyte Puccinellia tenuiflora under salt stress. Acta Physiol. Plant. 37:100. doi: 10.1007/s11738-015-1846-3
Wang, S., Li, Z., Rui, R., Fan, G. S., and Lin, K. W. (2013). Cloning and characterization of a plasma membrane Na+/H+ antiporter gene from Cucumis sativus. Russ. J. Plant Physiol. 60, 330–336. doi: 10.1134/S102144371303014X
Wang, Y., Lü, J., Chen, D., Zhang, J., Qi, K., Cheng, R., et al. (2018). Genome-wide identification, evolution, and expression analysis of the KT/HAK/KUP family in pear. Genome 61, 755–765. doi: 10.1139/gen-2017-0254
Weinl, S., and Kudla, J. (2009). The CBL-CIPK Ca2+-decoding signaling network: function and perspectives. New Phytol. 184, 517–528. doi: 10.1111/j.1469-8137.2009.02938.x
Wu, S. J., Ding, L., and Zhu, J. K. (1996). SOS1, a genetic locus essential for salt tolerance and potassium acquisition. Plant Cell 8, 617–627. doi: 10.1105/tpc.8.4.617
Xie, K., Zhang, J., and Yang, Y. (2014). Genome-wide prediction of highly specific guide RNA spacers for CRISPR–Cas9-mediated genome editing in model plants and major crops. Mol. Plant 7, 923–926. doi: 10.1093/mp/ssu009
Xing, H.-L., Dong, L., Wang, Z., Zhang, H.-Y., Han, C.-Y., Liu, B., et al. (2014). A CRISPR/Cas9 toolkit for multiplex genome editing in plants. BMC Plant Biol. 14:327. doi: 10.1186/s12870-014-0327-y
Yang, O., Popova, O., Süthoff, U., Lüking, I., Dietz, K.-J., and Golldack, D. (2009). The Arabidopsis basic leucine zipper transcription factor AtbZIP24 regulates complex transcriptional networks involved in abiotic stress resistance. Gene 436, 45–55. doi: 10.1016/j.gene.2009.02.010
Yoo, J., Park, C. Y., Kim, J., Heo, W. D., Cheong, M., Park, H., et al. (2005). Direct interaction of a divergent CaM isoform and the transcription factor, MYB2, enhances salt tolerance in Arabidopsis. J. Biol. Chem. 280, 3697–3706. doi: 10.1074/jbc.M408237200
Zhang, Y. M., Zhang, H. M., Liu, Z. H., Li, H. C., Guo, X. L., and Li, G. L. (2015). The wheat NHX antiporter gene TaNHX2 confers salt tolerance in transgenic alfalfa by increasing the retention capacity of intracellular potassium. Plant Mol. Biol. 87, 317–327. doi: 10.1007/s11103-014-0278-6
Keywords: salt stress, VIGS, GmNHX5, CRISPR/Cas9, hairy roots
Citation: Sun TJ, Ma N, Wang CQ, Fan HF, Wang MX, Zhang J, Cao JF and Wang DM (2021) A Golgi-Localized Sodium/Hydrogen Exchanger Positively Regulates Salt Tolerance by Maintaining Higher K+/Na+ Ratio in Soybean. Front. Plant Sci. 12:638340. doi: 10.3389/fpls.2021.638340
Received: 06 December 2020; Accepted: 26 January 2021;
Published: 09 March 2021.
Edited by:
Honghong Wu, Huazhong Agricultural University, ChinaReviewed by:
Koushik Chakraborty, National Rice Research Institute (ICAR), IndiaJian Sun, Jiangsu Normal University, China
Copyright © 2021 Sun, Ma, Wang, Fan, Wang, Zhang, Cao and Wang. 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) and the copyright owner(s) 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: Jie Zhang, emhhbmdqaWV6aGFvbGhAMTYzLmNvbQ==; Jinfeng Cao, Y2FvamluZmVuZzIwMDNAMTYzLmNvbQ==; Dongmei Wang, ZG9uZ21laXdhbmc2M0AxMjYuY29t
†These authors have contributed equally to this work and share first authorship