Xiaokang Fu1†Zhengying Lu2†Hengling Wei1Jingjing Zhang1Xu Yang1
Aimin Wu1Liang Ma1Meng Kang1Jianhua Lu1
Hantao Wang1*
Shuxun Yu1*- 1State Key Laboratory of Cotton Biology, Institute of Cotton Research, Chinese Academy of Agricultural Sciences, Anyang, China
- 2Handan Academy of Agricultural Sciences, Handan, China
The sodium/hydrogen antiporter (NHX) gene family with the Na+/H+ exchange protein domain is a transporter of sodium and hydrogen ions and plays an important role in the response of plants to salt stress. Studying the response of cotton to salt stress through comprehensive identification and analysis of NHX genes in several species and their roles in salt tolerance mechanisms is of great significance. In this study, 23, 24, 12, and 12 NHX genes were identified from Gossypium hirsutum (Gh), G. barbadense, G. arboreum and G. raimondii, respectively. Phylogenetic analysis showed that these genes were mainly divided into three clades with significant subcellular localization, namely, endosome (Endo-class), plasma membrane (PM-class) and vacuole (Vac-class). By analyzing the structure of NHX genes and proteins, each branch of the NHX gene family was found to be structurally conserved, and collinearity analysis showed that NHX genes were mainly expressed through whole genome and segmental duplication. The non-synonymous (Ka)/synonymous (Ks) values showed that the NHX gene family experienced strong purifying selection during long-term evolution. Cis-acting element analysis showed that the NHX gene family may be related to the regulation of abscisic acid (ABA) and methyl jasmonate (MeJA) hormones. Additionally, transcriptomic data analysis and qRT-PCR showed that GhNHXs exhibited different expression patterns in each tissue and under different salinities. These results provide an important reference for us to further understand and analyze the molecular regulation mechanism of cotton NHX genes.
Introduction
Soil salinization is one of the limiting factors of global crop yield, and approximately 800 million serious salinization events occur worldwide (Rengasamy, 2006; Munns and Tester, 2008). Increasingly serious soil salinization is an urgent problem for China’s agriculture (Flowers, 1999; Tomoaki and Julian, 2004; Heng et al., 2012; Deinlein et al., 2014). Salt stress can seriously affect plant growth and development and other normal life activities, resulting in the loss of crop yield and quality (Sun et al., 2019). Therefore, mining salt-tolerant genes, studying their functions through molecular analysis, and improving crop salt tolerance through molecular breeding are of great significance.
Plants activate a series of coping mechanisms when under salt stress, such as the use of hormones to regulate plant growth and metabolism, ion homeostasis and osmotic regulation (Wu G. Q. et al., 2019). Previous studies have found that the sodium/hydrogen antiporter (NHX) protein is key to the regulation mechanism of ion homeostasis (Pehlivan et al., 2016; Cao et al., 2020). The NHX genes mainly use two proton pumps formed by the H+-ATP enzyme and H+-PPase to produce H+ electrochemical gradients to transport Na+ from the cytoplasm to vacuoles or outside the cell, thereby maintaining the stability of Na+ ions and avoiding the toxic effect of Na+ accumulation in cells (Munns, 2002; Zorb et al., 2005; Sze and Chanroj, 2018). NHX proteins are ubiquitous transmembrane proteins belonging to the monovalent cation/proton antiporter 1 (CPA1) superfamily (Brett et al., 2005; Chanroj et al., 2012), and most NHX proteins contain 10–12 transmembrane helix domains (TMs) (Yamaguchi et al., 2003; Wu G. Q. et al., 2019). The NHX proteins are mainly localized in vacuoles, endosomes and plasma membranes (Aharon et al., 2003; Pehlivan et al., 2016). The first plant NHX gene was identified by Roberto et al. (1999) in Arabidopsis thaliana (At) (Roberto et al., 1999), and a total of 8 NHXs have been identified in this species to date. Of these NHXs, 2 genes (AtNHX7 and AtNHX8) are located in the PM-class (plasma membranes), 2 genes (AtNHX5 and AtNHX6) are located in the Endo-class (endosomes), and 4 genes (AtNHX1-4) are located in the Vac-class (vacuoles) (Shi et al., 2000; Aharon et al., 2003; Brett et al., 2005; Bassil et al., 2011).
NHX genes are involved in the salt stress response, cell expansion, membrane vesicle trafficking, and pH homeostasis (Dragwidge et al., 2019; Jegadeeson et al., 2019). Developmental disorders and cell division abnormalities in embryo and root tissues were observed after double knockout of AtNHX5 and AtNHX6 in A. thaliana (Dragwidge et al., 2018). In rice, salt stress, hyperosmotic stress and abscisic acid (ABA) stress induced expression of OsNHX1, OsNHX2, OsNHX3 and OsNHX5, and the differential expression of NHX genes in different tissues may be an important factor determining the salt tolerance of rice (Fukuda et al., 2011). Salt-tolerant perennial ryegrass was obtained by transforming the rice vacuolar membrane OsNHX1 gene, and the plant was found to be able to survive under high salt conditions (350 mM) for 10 weeks (Wu et al., 2005). Transgenic Brassica napus plants overexpressing a vacuolar NHX gene were able to grow, flower, and produce seeds in the presence of 200 mM NaCl (Zhang et al., 2001). Vac-class membrane gene AgNHX1 of Atriplex gmelinii was transferred into rice, and the transgenic plants could survive under conditions of 300 mM NaCl for 3 days, while the wild-type rice plants could not (Ohta et al., 2002). In addition, transfer of the AtNHX1 gene from A. thaliana into tomato (Eduardo et al., 2009) and transfer of the Vac-class membrane gene SsNHX1 (Salsola collina Pall.) into Medicago sativa (Li et al., 2011) revealed that all plants with high expression of receptors significantly improved their salt tolerance.
Cotton is an important economic crop, and with the gradual transfer of cotton agriculture to Xinjiang (China) and saline-alkali beaches, improving the salt tolerance and stress resistance of cotton is very important (Gao et al., 2018; Dong et al., 2019). Some studies have shown that under a high concentration of NaCl (200 mM), NHX genes may play an important role in the salt tolerance of cotton (Kumar et al., 2017; Elsawy et al., 2018; Mushke et al., 2019). For example, cotton overexpressing AtNHX1 can produce more biomass and fiber (He et al., 2005), and the expression level of GhSOS1 (a plasma membrane NHX gene) in cotton roots is significantly increased (Chen et al., 2017). However, at present, research on the number and salt tolerance of NHX genes in cotton is limited. Therefore, we selected four cotton species [Gossypium arboretum (Ga), G. raimondii (Gr), G. hirsutum (Gh) and G. barbadense (Gb)] as the research objects, and a total of 71 NHX genes were identified. Through systematic analysis of the gene distribution, gene structure, phylogenetic relationships and cis-acting elements of the NHX family, the purpose of this paper is to provide a theoretical reference for the study of salt tolerance and the improvement of new cotton varieties in the future.
Results
Chromosomal Distribution of NHXs in Gossypium spp.
A total of 23, 24, 12, and 12 NHX family genes were identified in G. hirsutum, G. barbadense, G. raimondii, and G. arboreum, respectively. The NHX genes were named according to their homology with A. thaliana genes from high to low (Supplementary Table S1). The position of each gene on the chromosome is shown in Figure 1. Through identification and analysis of the NHX genes of the four cotton species, their locations and distributions in chromosomes and the numbers of genes were found to be conserved in the evolution from diploid to tetraploid cotton species.
Figure 1. NHXs are distributed in the chromosomes of four cotton species (G. arboreum, G. raimondii, G. barbadense, and G. hirsutum), and the chromosome numbers are shown on the left side of each vertical bar.
Identification of the Physicochemical Properties of NHX Genes
The allotetraploid species G. hirsutum and G. barbadense have 11 and 12 NHX genes in the At subgenome, respectively, and 12 NHX genes in the Dt subgenome. Diploid species G. arboreum and G. raimondii each have 12 NHX genes (Figure 1). The transcripts of 71 genes are 963–7055 bp, and the protein sequences are 320–1166 aa. The isoelectric point (pI) of the protein is between 5.237 and 10.181, and the molecular weight is between 35.548 kDa and 129.861 kDa. Among 71 proteins, the number of TMs varies from 5 to 12, and that of 53 proteins ranges from 10 to 12 (Supplementary Table S2).
Through analysis of the number of phosphorylation sites of 75 protein sequences in G. hirsutum, G. barbadense, G. arboreum, and G. raimondii, 17 (GhNHX_21) – 76 (GhNHX_1 and GbNHX_1) serine (S) and 4 (GhNHX_20) – 24 (GhNHX_2 and GbNHX_2) threonine (Th) sites were found; GhNHX_20 and GbNHX_24 did not contain tyrosine (Ty), and GhNHX_1, GhNHX2, GbNHX_1, GbNHX_2, GaNHX_1, and GrNHX_1 each contained 9 tyrosine sites, which was the highest number of sites. Among the three kinds of proteins, the content of serine was higher than those of threonine and tyrosine, indicating that serine is a common phosphorylation site in the upland cotton NHX family. In addition, the phosphorylation degrees of protein kinase C (PKC) and UNSP in NHX proteins were the highest, and the phosphorylation degrees of PKB, SRC, GSK3, and ATM were lower (Supplementary Table S3).
Phylogenetic Analysis of NHX Sequences
To study the evolutionary relationships of the NHX gene family, a total of 106 NHX protein sequences from six species [G. hirsutum, G. raimondii, G. barbadense, G. arboreum, Arabidopsis and Populus trichocarpa (Pt)] were used to construct a phylogenetic tree (Figure 2). Based on the results of the phylogenetic tree combined with those of previous studies in Arabidopsis, 106 genes were mainly divided into three categories, namely, class I (Endo-class), class II (PM-class) and class III (Vac-class). Class III (Vac-class) was the largest group among the 6 species, which included 71 proteins (17 NHXs in G. hirsutum, 18 NHXs in G. barbadense, 9 NHXs in G. arboreum, 9 NHXs in G. raimondii, 4 NHXs in Arabidopsis and 14 NHXs in P. trichocarpa). Seventeen and 18 NHXs were distributed in class I (Endo-class) and class II (PM-class), respectively.
Figure 2. Phylogenetic tree of the NHX proteins in six species. Three categories of NHXs are indicated using different colors by iTOL.
Evolutionary Analysis of NHX Genes in Gossypium spp.
The number of NHX genes in the two diploid cotton species, G. raimondii (12) and G. arboreum (12), was approximately half that in the tetraploid cotton G. hirsutum (23) and G. barbadense (24) (Figure 2). The NHX gene family was further verified to have emerged due to chromosome doubling during cotton evolution from diploid to tetraploid. Therefore, we speculate that during evolution, diploids may have become tetraploids through whole-genome duplication (WGD). To verify this hypothesis, we used MCScanX software to analyze the collinearity of the At and Dt subgenomes of G. hirsutum and the corresponding ancestral A and D diploid genomes. The results showed that the homologous gene pairs of GhNHX genes were in the collinearity analysis section between different genomes (Figure 3). These results provide strong evidence for the hypothesis that NHX genes were obtained by whole-genome duplication (WGD) or segmental duplication, and we analyzed the duplication types of the NHX genes in G. hirsutum and G. barbadense, as shown in Supplementary Table S4. The results showed that 46 of 47 NHXs were obtained by WGD or segmental duplication during evolution from diploid to tetraploid cotton species.
Figure 3. Collinearity analysis between G. arboreum (AA), G. hirsutum (At and Dt subgenomes) and G. raimondii (DD) by MCScanX. Red and gray lines depict NHX homologous gene pairs and WGD pairs, respectively.
To study the selective pressure of NHX genes, homologous gene pairs were determined from the results of the gene sequence comparison program BLASTP and collinearity analysis. In this study, the homologous genes from diploid species to the heterotetraploid subgenome were mainly analyzed. A total of 44 pairs of homologous genes were identified by collinearity analysis between the allotetraploid (G. barbadense and G. hirsutum) genomes and the corresponding diploid genomes (G. raimondii and G. arboreum) (Figure 3 and Supplementary Table S5). Eleven homologous gene pairs each were identified in the G. barbadense and G. hirsutum subgenomes and the corresponding G. arboreum (Ga) genomes, except for GaNHX_7. The other 11 NHX genes all had homologous gene pairs in tetraploid cotton species, and 22 homologous gene pairs were found in G. raimondii and the corresponding G. barbadense and G. hirsutum Dt subgenomes. Among the 44 homologous gene pairs identified, the selection pressure of these gene pairs was evaluated by calculating non-synonymous (Ka)/synonymous (Ks) values, and 42 pairs had a Ka/Ks ratio < 1 (Supplementary Table S5). The results showed that during evolution, most NHX genes were strongly purified and selected, and only a few genes differentiated and produced new biological functions.
Conserved Sequence and Gene Structure Analysis
To understand the structural diversity and structural characteristics of NHX genes, motif analysis of 71 NHX amino acid sequences of four cotton genomes (G. hirsutum, G. barbadense, G. arboreum, and G. raimondii) was carried out by the MEME program (Figure 4). A total of 6 motifs were identified in the NHXs, and only motif 3 existed in all genes. Motifs 2, 3, 4, and 6 appeared in PM-class NHX genes. Most of the Endo-class genes contained motifs 1, 3, 4, 5, and 6, while GhNHX_23 contained only motifs 1, 3, and 5. Of the genes in the Vac-class, 39 genes, 10 genes, 3 genes and 1 gene contained 6, 5, 4, and 3 motifs, respectively. Figure 4 shows that the numbers and types of motifs in the NHXs were significantly different between groups and were almost similar within groups, proving that the amino acid sequences and gene structures of NHX genes were highly conserved. These results serve as a basis for predicting the function of the NHX gene family.
Figure 4. Gene structures and conserved protein motifs of NHXs. (I) NJ phylogenetic tree analysis of G. hirsutum. (II) Distributions of the predicted motifs in the NHX genes. (III) The numbers, lengths, and positions of exons and introns within NHX genes.
The exons and introns of G. hirsutum, G. barbadense, G. arboreum and G. raimondii were analyzed and compared by the online program GSDS1 (Figure 4). The numbers of exons of GhNHX_1, GhNHX_2, GaNHX_1, GbNHX_1, GbNHX_2, and GrNHX_1 located on the cell membrane were the highest (23–24). A total of 12 genes in the four species were located in the Endo-class, 11 of which had between 18 and 23 exons, and only GrNHX_8 had 24 exons, while the number of exons of the other genes in the Vac-class varied between 8 and 16. The numbers of exons in the different classes were obviously different, ranking from lowest to highest as PM-class < Endo-class < Vac-class. The numbers of exons and introns between species were relatively conserved.
Analysis of Cis-Acting Elements in NHX Promoter Regions
To better study gene expression and transcriptional regulation, 71 NHX promoter sequences of the four cotton species were identified by PlantCARE software (Figure 5 and Supplementary Table S6), and mainly stress- and hormone-related cis-acting elements were analyzed. Among the identified cis-acting elements, 545 were related to hormones, including ABA, salicylic acid (SA), auxin (IAA), gibberellin (GA) and jasmonate. Of these, ABA (163) and methyl jasmonate (MeJA) (246) were the largest types, accounting for 29.9 and 45% of all identified cis-acting elements of hormones, respectively. A total of 690 cis-acting elements were related to stress, such as low temperatures, drought, trauma, and anaerobic conditions. Among these elements, the contents of antioxidant- and defense-related reaction elements were higher. In addition, among the four cotton species of G. hirsutum, G. barbadense, G. arboreum, and G. raimondii, W BOX (a salt-responsive element) predicted 9, 10, 4, and 6 elements, respectively (Figure 5B and Supplementary Table S6). The results of these cis-acting elements suggest that NHXs may play an important role in hormone regulation and stress responses in cotton.
Figure 5. Distribution characteristics of phytohormone-responsive and stress-related cis-acting elements in NHX promoter fragments of four cotton species (G. hirsutum, G. barbadense, G. arboreum and G. raimondii). (A) The sum of five phytohormone-responsive cis-acting elements [methyl jasmonate (MeJA), abscisic acid (ABA), auxin (IAA), gibberellin (GA), and salicylic acid (SA)] in the NHX promoter fragments of four cotton species. (B) The number of cis-acting elements related to stress in the four cotton species.
Expression Characterization of GhNHXs in Different Tissues
To study the biological function of GhNHXs in different tissues, we analyzed spatio-temporal expression patterns in root, stem, leaf, anther, filament, pistil, petal, and torus tissues from transcriptome data (Figure 6A). The Figure 6A showed that 2 genes (GhNHX_12 and GhNHX_13) had low or no expression in the 8 tissues. GhNHX_20 and GhNHX_21 were predominantly expressed in pistil; GhNHX_17 and GhNHX_22 were predominantly expressed in petal; GhNHX_2 and GhNHX_7 were predominantly expressed in torus. GhNHX_1, GhNHX_4, GhNHX_8, GhNHX_9, GhNHX_15, and GhNHX_23 were highly expressed in root and stem; GhNHX_3, GhNHX_10, GhNHX_11, and GhNHX_18 were highly expressed in anther, filament or petal.
Figure 6. Expression profiles of GhNHXs in different tissues (A) and responses to different stresses (B). The tissues or treatments are shown at the bottom, genes are shown on the right, and the phylogenetic relationships are shown on the left.
Clustering analysis of the 23 GhNHX genes treated with polyethylene glycol (PEG) or NaCl (Figure 6B) showed that the expression patterns in transcriptome data from the two treatment conditions were different. In PEG treatment, the expression of GhNHX_1, GhNHX_2, GhNHX_6, GhNHX_11, and GhNHX_19 reached the highest level at 12 h after treatment. In salt treatment, the expression of GhNHX_16 and GhNHX_17 showed an overall upward trend from 0 h to 24 h. And GhNHX_3, GhNHX_9, GhNHX_10, GhNHX_18 and GhNHX_22 responded quickly to salt stress at 1 h. The expression of the GhNHX_1, GhNHX_2, GhNHX_8, GhNHX_11, GhNHX_15, GhNHX_16, GhNHX_17, GhNHX_19, GhNHX_20, GhNHX_21, and GhNHX_23 were high at 6 or 12 h. At the same time, the Figure 6B also suggested that 4 genes (GhNHX_1, GhNHX_2, GhNHX_6, and GhNHX_19) showed obvious responses under PEG and NaCl treatments, indicating that these genes may be regulated by drought and salt stresses.
According to the transcriptome data, six genes were selected for qRT-PCR analysis (Figure 7). GhNHX_1 and GhNHX_2 belong to the PM-class, GhNHX_20 and GhNHX_23 belong to the Endo-class, and GhNHX_3 and GhNHX_5 belong to the Vac-class. The results of qRT-PCR in 4 tissues showed that GhNHX_2 and GhNHX_3 were highly expressed in petal; GhNHX_1 and GhNHX_20 were more expressed in stem; GhNHX_5 was more expressed in leaf and petal; GhNHX_23 was more expressed in root and stem. The results showed that GhNHXs were expressed in different tissues.
Figure 7. Expression analysis of GhNHXs in different tissues. Error bars show the standard deviation of three biological replicates.
Expression Analysis of GhNHX Genes Under Salt Stress
By analyzing the cis-acting elements of the promoter, we speculated that GhNHXs may be related to the stress response. According to the analysis of transcriptome data, some genes have relatively high expression under different stresses (Figure 6B). Most genes responded to NaCl stress, and different genes had different responses to different treatment times. For example, the expression of GhNHX_3 and GhNHX_18 was relatively high under NaCl stress for 1 h and decreased gradually with increasing treatment time, and the expression of GhNHX_16 and GhNHX_17 gradually increased with increasing treatment time. In addition to the response to salt treatment, a response to PEG treatment was also observed. These results showed that GhNHXs can improve resistance to abiotic stress through different expression patterns in response to stress and provide candidate genes for further investigation of salt stress.
To further analyze the response of GhNHXs to salt stress, GhNHX_11 and GhNHX_20 (Endo-class), GhNHX_1 and GhNHX_2 (PM-class), and GhNHX_5 and GhNHX_17 (Vac-class) were selected for verifying by qRT-PCR. As shown in Figure 8, GhNHX genes were induced during salt treatment. Under the same salt concentration, the expression trends in roots and leaves were different. Under low-salt treatment (50 mM NaCl), the expression of GhNHX_2 and GhNHX_5 increased at first and then decreased in leaves, and the expressions of GhNHX_1, GhNHX_2, GhNHX_5, GhNHX_17, and GhNHX_20 were the highest at 48 h in roots. Under high-salt stress (200 mM NaCl), the expressions of GhNHX_1, GhNHX_2, GhNHX_11, GhNHX_17, and GhNHX_20 reached the maximum at the later stage of treatment (12–48 h), and the expressions of GhNHX_11 and GhNHX_17 increased gradually with the extension of treatment time in leaves. In roots, the expressions of GhNHX_1 and GhNHX_2 were higher at 3 h and lower at 24 h, and the expressions of GhNHX_5 and GhNHX_17 were higher at 9 and 48 h. The expression of GhNHX_11 maintained a high level from 6 to 48 h. These results suggested that these genes might play an important role in response to salt stress.
Figure 8. Relative expression levels of GhNHX genes in the leaves and roots of G. hirsutum TM-1 seedlings subjected to 50 and 200 mM NaCl for 0, 3, 6, 9, 12, 24, and 48 h. The 2–ΔΔCt method was used to calculate the expression levels of GhNHX genes at different times. Error bars show the standard deviation of three biological replicates.
Responses of GhNHX Genes to MeJA and ABA Treatment
Among the predicted promoters, the contents of MeJA and ABA hormone response elements were found to be the highest in the GhNHX gene family (Figure 5A). To further explore the effect of hormones on the NHX family, 6 GhNHX genes containing these two hormone response elements were selected for qRT-PCR analysis. The results showed that after MeJA hormone treatment, compared with the control at the same time point, all 6 genes reached maximum levels at the time point after MeJA treatment, 3 of which (GhNHX_2, GhNHX_3, and GhNHX_21) showed a trend of continuous upregulation (Figure 9), reaching maximum levels at 12 h, while the other three genes (GhNHX_4, GhNHX_8, and GhNHX_17) reached maximum levels at 9 h. The results from ABA treatment showed that GhNHX_2 and GhNHX_4 reached maximum levels at 6 h after treatment, GhNHX_3 and GhNHX_21 decreased slightly, and GhNHX_8 showed an upward trend after treatment (Figure 9). These results suggest that the transcriptional levels of GhNHXs may be closely related to the regulation of MeJA and ABA.
Figure 9. Relative expression levels of GhNHX genes under ABA and MeJA treatments. Two true leaves were sprayed with 200 mM ABA (red) and 100 mM MeJA (green), and water (blue) was used as a blank control. Error bars show the standard deviation of three biological replicates.
Discussion
The NHX family plays an important role in plant cells and is indispensable for normal life activities. NHX proteins can help absorb K+ ions under normal conditions and can be stored in vacuoles for cellular storage, turgor generation and pH regulation. Under salt or osmotic stress, NHX proteins fulfill a protective function through vacuolar compartmentalization of K+ as well as Na+ in some cases, thereby preventing toxic Na+-K+ ratios in the cytosol while accruing solutes to balance both (Jiang et al., 2010).
Gene Expansion and Duplication of NHXs in the Evolution of Cotton
The production of polyploids is an expression of plant adaptability to the environment and an important mechanism for the formation of new species (Ramsey and Schemske, 1998). The replication of a single gene, chromosome or whole genome is the main force of plant genome evolution (Paterson et al., 2012). Data indicate an origin for Gossypium 5–15 million years ago (mya) and rapid early diversification of the major genome groups (Wendel and Cronn, 2003). Allotetraploid cotton (G. hirsutum) is the result of genomic hybridization and doubling approximately 1–1.5 mya (Ks peaks at 0.005 and 0.008, respectively) from diploid A (G. arboreum) originating in Africa and diploid D (G. raimondii) originating in Mexico (Wendel and Cronn, 2003; Li et al., 2015; Zhang et al., 2015), and gene loss is most likely an ongoing process in allotetraploid cotton (Zhang et al., 2015), which may have resulted in G. hirsutum (23) lacking an NHX gene compared with G. barbadense (24) and provides strong support for the study of cotton polyploidy. In this research, 23, 24, 12, and 12 NHX genes were identified in G. hirsutum, G. barbadense, G. arboreum and G. raimondii, respectively. The results showed that the number of NHXs in two tetraploid cotton species was almost the sum of the NHXs in two diploid cotton species, which resulted in more gene duplication in the tetraploid NHX gene family than in diploid cotton species, and the tetraploid genome was further verified to have been formed by doubling of the chromosomes of the two diploid genomes (Paterson et al., 2012). Although more NHXs were present in tetraploid cotton species than in diploid cotton species, genomic arrangements and chromosome doubling may have also led to some gene loss in the process of doubling from diploid to tetraploid cotton species (Paterson et al., 2012). The duplication types of 47 genes in G. hirsutum and G. barbadense showed that 46 NHX genes were formed by WGD or segmental duplication based on collinearity analysis (Supplementary Table S4). Additionally, cotton MADS-Box, the GT47 family and soybean WRKY showed enlargement as the result of WGD and segmental duplication (Yin et al., 2013; Ren et al., 2017; Wu A. et al., 2019), indicating that gene duplication plays an important role in the evolution of gene families (Cannon et al., 2004). The collinearity analysis of NHXs showed that the A genome and D genome of NHXs were similar to the At subgenome and Dt subgenome, respectively (Figure 3), indicating that they share the same ancestor (Qanmber et al., 2019). In conclusion, WGD or segmental duplication should be the main driving force for the NHX family from the diploid genome to the tetraploid genome, and these findings also provide a theoretical reference for the study of gene duplication, genetic information transmission between species and species evolution.
Evolutionary Conservation and Divergence of NHX Genes
In this study, 8, 27, 23, 24, 12, and 12 NHXs were identified in the genomes of Arabidopsis, P. trichocarpa, G. hirsutum, G. barbadense, G. arboreum, and G. raimondii, respectively. The differences in the number of NHX genes in plants may be attributed to gene duplication and loss specific to different subfamilies of the NHX gene over the course of evolution (Wu G. Q. et al., 2019). In Arabidopsis, 8 NHXs were divided into three groups according to subcellular location, including 2 PM-class genes, 2 Endo-class genes and 4 Vac-class genes (Brett et al., 2005). In our research, 106 NHX genes from 6 species were also classified into 3 categories through phylogenetic tree analysis (Figure 2). Among them, the numbers of GhNHXs in these categories were 17, 2, and 4, respectively, and the number of genes in the Vac-class was the largest, indicating that NHXs in the Vac-class may play an important role in cells. Conservation of the gene sequence structure is the basis of the conservation of biological function. In this study, to explore the characteristics of NHX gene sequences in cotton, 71 NHX genes in 4 cotton species were analyzed for exon/intron structures and motifs (Figure 4), and the results showed that the structures and motifs of the NHX genes were similar. In G. hirsutum, the numbers of exons of GhNHX_1 and GhNHX_2 in class II (PM-class) were the largest, the numbers of exons of GhNHX_23, GhNHX_20, GhNHX_19, and GhNHX_11 were the second largest, and the numbers of exons of GhNHX_18, GhNHX_21 and GhNHX_23 were the lowest, all of which were less than or equal to 16. The numbers of exons of the three NHX types (PM-class, Endo-class and Vac-class) in poplar were 23, 22, and 14, respectively (Tian et al., 2017). The largest number of exons in sugar beet was found in BvNHX4 (Endo-class), with 26 exons, followed by BvNHX5 (PM-class) with 23 exons (Wu G. Q. et al., 2019). In summary, structural similarities and quantitative diversity exist in exons and introns among different species, not only indicating that the NHX gene family conserved its gene function but also reflecting the diversity of selection in the process of evolution.
Previous studies have shown that most of the N-terminal domains of NHX proteins contain 9–12 TMs, which are highly conserved (Yamaguchi et al., 2003; Wu G. Q. et al., 2019). The domain composed of the amino acid sequence FFIYLLPPI was called the amiloride-binding site, which was sensitive to the Vac-class membrane NHX in the presence of the drug amiloride or its derivatives and was responsible for transport (Wu et al., 2011). In this study, 12 GhNHXs had this binding site in G. hirsutum, almost all of which were in the N-terminal TM3 of the Vac-class and 9–12 transmembrane fragments. The same results were also observed in AtNHX1-4 of Arabidopsis (Aharon et al., 2003) and ZxNHX1 of Zanthoxylum (Wu et al., 2011). These results indicated that the NHX reverse transporters of the Vac-class were highly conserved during evolution. In this study, through analysis of selective stress in the evolution of the NHX gene family, the results showed that the Ka/Ks ratio was less than 1 (Supplementary Table S5), indicating that the NHX gene family was strongly purified and selected during long-term evolution and is functionally conserved.
Cis-Acting Elements of NHX in Upland Cotton
In eukaryotes, transcriptional regulation is the main mechanism of plant gene expression regulation, and cis-acting elements are involved in the transcriptional regulation of genes and play an important role in various biological processes, such as hormonal responses, abiotic stress responses and development (Ding et al., 2018; Verma et al., 2019; Wu A. et al., 2019). Different cis-acting elements have particularities and consistency in the promoter. The cis-acting element containing the AuxRE and DR5 motifs is usually an auxin-induced promoter (Ulmasov et al., 1997), light-induced promoters usually contain G-box, AT-rich, GT1-motif and I-box elements (Gilmartin et al., 1992; Foster et al., 1994), and the promoter containing CATGTG and CACG cis-acting elements is drought-induced (Tran et al., 2004). In plants, hormones such as ABA, ethylene, SA and IAA play important roles in plant growth and development and responses to stress (Díaz et al., 2002; Mishra et al., 2014; Liu et al., 2015; Shu et al., 2019; Wang and Huang, 2019). In this study, 71 NHX promoter fragments from four cotton species were extracted, and cis-acting elements were identified. We mainly statistically analyzed stress- and hormone-related cis-acting elements, as shown in Figure 5. The results showed that the cis-acting elements related to hormone responses were IAA, GA, SA, ABA, and MeJA; these five hormones existed in all four cotton species, and the contents of ABA and MeJA were the highest (Supplementary Table S6). qRT-PCR analysis revealed that exogenous application of MeJA and ABA had different effects on the induction of GhNHX gene transcription, especially exogenous MeJA treatment, and six genes were upregulated to different degrees and reached maximum levels at different time points (Figure 9). NHXs have been suggested to be regulated by many hormones, and ABA and MeJA hormones play more extensive roles in hormone regulation. In addition, cis-acting elements related to stress responses were identified, mainly including responses to low temperature, drought, wounding, MYB binding and anoxia. The proportion of each cis-acting element in different cotton species was similar, as shown in Supplementary Table S6, indicating that the NHX family may have the same regulation mode in cotton (Wu A. et al., 2019). In addition, W BOX, a cis-acting element associated with development processes and the salt response (Xu et al., 2018; Wu G. Q. et al., 2019), was predicted in all four cotton species (G. hirsutum, G. barbadense, G. arboreum, and G. raimondii) (Figure 5B), which can be recognized by the family of WRKY transcription factors. The W BOX cis-element may have important reference value in salt stress in the NHX family. In summary, the NHX family may be involved in the regulation of a variety of hormones and stress responses. These cis-acting elements enable plants to rapidly respond to external biological and abiotic stresses and improve the resistance of plants to the outside world, facilitating further understanding of the related mechanism of cis-acting regulatory elements in the stress response.
Expression Analysis of GhNHX Genes Under Multiple Stresses
One important method to alleviate salt stress in plants is to isolate Na+ into vacuoles to maintain Na+ homeostasis and reduce the toxicity of Na+ to the cytoplasm (Li et al., 2017; Wu, 2018). Under PEG treatment, the expression of GhNHX_1, GhNHX_2, GhNHX_6, GhNHX_11, and GhNHX_19 reached the highest level at 12 h, and GhNHX_1, GhNHX_2, and GhNHX_19 also responded to salt stress in varying degrees, demonstrating these genes may be associated with salt stress and drought stress. Some studies verified that ABA and MeJA were related to salt tolerance (Zhang et al., 2019; Zhao et al., 2019). Depending on the qRT-PCR results of exogenous MeJA and ABA treatment, it was showed that GhNHX_3, GhNHX_8, GhNHX_17, and GhNHX_21 were highly expressed at one specific time after MeJA treatment; GhNHX_4 was more expressed at 6 h, which was 7–8 times higher than that at 12 h after ABA treatment (Figure 9), illustrating these genes may respond to MeJA and ABA hormone regulation. In addition, we identified the relative expression of 6 GhNHXs under salt stress by qRT-PCR. Under high-salt conditions (200 mM NaCl), there were differences between roots and leaves in the expression trends of GhNHX genes. The expression of GhNHX_1 and GhNHX_2 (PM-class) increased gradually within 12 h of salt treatment in leaves and remained at a high level from 12 to 48 h, while the expression in roots reached the highest level at 3 h, decreased gradually, and then increased again at 48 h (Figure 8). Similar results were observed in the BvNHX genes of sugar beet (Wu G. Q. et al., 2019). Similarly, the expression of the wheat gene TaNHX3 increased within 24 h in roots or leaves and decreased gradually within 48 h, and transgene analysis confirmed that TaNHX3 has important functions in regulating the plant tolerance to salt stress (Lu et al., 2014). The expression of the Porteresia coarctata gene PcNHX1 in roots increased gradually within 24 h of salt treatment and then decreased gradually from 36 to 48 h and transgenic tobacco PcNHX1 seedlings show growth advantage in increasing NaCl (Jegadeeson et al., 2019). Interestingly, during root treatment (Figure 8), the expression of GhNHX_5, GhNHX_17, and GhNHX_20 decreased rapidly from 0 to 3 h under high-salt conditions (200 mM NaCl) and increased with increasing salt treatment time in contrast to the relatively low expression at 0 h under low-salt conditions, which may be related to regulation of the salt tolerance of upland cotton TM-1 with different concentrations of salt treatment. Under low-salt conditions (50 mM NaCl), The expression level of GhNHX_1, GhNHX_2, GhNHX_5, GhNHX_17, and GhNHX_20 was the highest at 48h. The expression of BvNHX1 and BvNHX3 in sugar beet were significantly upregulated in NaCl solutions of different concentrations within 48 h (Wu G. Q. et al., 2019). According to the above analysis results, it was found that these genes (GhNHX_1, GhNHX_2, GhNHX_5, GhNHX_17, and GhNHX_20) were closely related to the response to salt stress, but their regulatory mechanism needs to be further studied.
Materials and Methods
Identification and Characterization of the NHX Genes
Eight Arabidopsis ATNHX sequences were obtained from the Arabidopsis thaliana database TAIR (Supplementary Table S7) (Lamesch et al., 2012). The protein and genome sequences of G. arboreum (CRI, version 1.0), G. raimondii (JGI, version 2.1), G. hirsutum (HAU, version 1) and G. barbadense (HAU, version 1) were downloaded from the cotton genome database CottonFGD (Zhu et al., 2017). The protein and genome sequences of P. trichocarpa were downloaded from the NCBI database.2 The hidden Markov model (HMM) profile was constructed using eight ATNHX genes, and the NHX genes in four cotton species and P. trichocarpa were identified by HMMER 3.0 software with an E value threshold (E) < 10–40 (Potter et al., 2018). Furthermore, the conserved domains of all the candidate NHX protein sequences were identified using the online Simple Modular Architecture Research Tool (SMART) (Letunic et al., 2015).
Transmembrane domain analysis of NHX sequences was performed using TMHMM 2.03 (Sonnhammer et al., 1998; Krogh et al., 2001). Phosphorylation sites were predicted using netphos 3.1 (Blom et al., 2004). The physicochemical properties of NHX proteins were analyzed by the ExPASy Proteomics Server (Bjellqvist et al., 1993).
Phylogenetic Analysis
To study phylogenetic relationships among different species, multiple sequence alignments of NHX family proteins from several species, including A. thaliana, P. trichocarpa, G. hirsutum, G. barbadense, G. arboreum, and G. raimondii, were carried out using the ClustalW function of MEGA 7.0 software, and the NJ method was used to compare the results. The bootstrap method setting value for the phylogenetic tree was 1000, and the default values were used for other parameters (Kumar et al., 2016). Visualization of developmental trees was performed by iTOL4 (Letunic and Bork, 2016).
Chromosome Distribution of NHX Genes and Analysis of the Ka/Ks Ratio
According to analysis of the cotton genome database, the physical locations of NHX genes in four cotton species [G. hirsutum (Gh), G. barbadense (Gb), G. arboreum (Ga), and G. raimondii (Gr)] were determined. Collinearity analysis of protein sequences of the four cotton species was performed using BLASTP and MCScanX software (Wang et al., 2012). TBtools software was used to visualize the chromosome location and collinearity results, and the NG function was used to estimate the Ks and Ka replacement rates (Yang and Nielsen, 2000; Chen et al., 2018).
Analysis of Conserved Motifs, Gene Structures, and Cis-Acting Elements
The conserved domains of NHX genes were analyzed using the online software MEME5.1.0,5 the maximum number of motifs was set to 6, and other parameters were set to default values (Bailey et al., 2009).
GSDS was used to analyze the members of the gene family by comparing the CDS sequence of NHX genes with the corresponding genome sequence (Hu et al., 2015). CDS sequence and genome sequence data packages were downloaded from CottonFGD,6 and the result of the TBtools software analysis was used to visualize the image (Chen et al., 2018). The 2000-bp upstream regions from the initiation codons (ATG) of NHXs were analyzed by PlantCARE software, and the cis-acting elements in the promoter were determined (Lescot et al., 2002).
Plant Materials and Treatments
The Upland Cotton Cultivar “Texas Marker-1” (TM-1) was planted in Anyang City, Henan Province, China. RNA was extracted from roots, stems, leaves and petals. After soaking TM-1 seeds in distilled water overnight, uniform seeds were transferred to plastic containers containing Hoagland nutrient solution and cultured in a greenhouse (light and dark cycle: 28°C for 16 h/22°C for 8 h) to study the reaction to NaCl treatment. When two true leaves were cultured, Hoagland nutrient solution containing 200 mM NaCl was used to treat them. Roots and leaves treated for 0, 3, 6, 9, 12, 24, and 48 h were used for RNA extraction.
Similarly, TM-1 seeds were planted in a greenhouse, and when they grew to two flat true leaves, they were sprayed with 100 mM MeJA or 200 mM ABA for treatment, and water was used as a control. Three leaves of seedlings were collected for RNA extraction at 3, 6, 9, and 12 h of treatment. All samples were immediately frozen in liquid nitrogen and stored at −80°C.
Transcriptome Data Analysis and Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR)
Raw RNA-seq data (SRA: PRJNA490626) were downloaded from the NCBI Sequence Read Archive. RNA-seq expression data were analyzed by TopHat2 (Kim et al., 2013) and cufflinks (Trapnell et al., 2010), and gene expression was measured in fragments per kilobase million (FPKM). A heat map was created by TBtools software. RNAprep Pure Plant Kit (TIANGEN, Beijing, China) was used to extract total RNA from collected samples. A PrimeScript RT Reagent kit (Takara, Japan) was used to synthesize cDNA. The sequences of qRT-PCR primers are shown in Supplementary Table S8. GhACTIN was used as a constituent expression control in qRT-PCR experiments. An ABI 7500 real-time PCR system (Applied Biosystems, United States) was used to perform qRT-PCR (Promega, Madison, WI, United States) with three biological replicates. The relative expression level of GhNHXs was calculated by the 2–ΔΔCT method.
Conclusion
The NHX gene family plays an important role in the study of plant salt tolerance. In this study, 71 members of the NHX family from four cotton species were identified, and their gene sequences, phylogeny and biological characteristics were analyzed by comprehensive bioinformatics. The purpose of this study was to explore the conservative structures, evolutionary relationships and regulatory elements of the NHX family and finally reveal this family’s important role in the regulation of salt stress in cotton. Through qRT-PCR analysis, we found that the expression of some NHXs showed an upward trend in cotton under salt stress, and the expression patterns of some genes were different under low- and high-salt conditions. These results indicate that the NHX family plays an important role in the response of cotton to salt stress, and the systematic bioinformatics analysis also serves as a solid foundation for the exploration of cotton NHX genes.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found here: SRA: PRJNA490626.
Author Contributions
XF and HWa conceived and designed the study and prepared the manuscript. XF, ZL, JZ, XY, HWe, and MK performed the experiments. HWa, AW, LM, and JL assisted with the analysis and interpretation of the data. SY participated in the design of the experiments and provided a critical review. All authors contributed to the article and approved the submitted version.
Funding
This research was financed by the National Key Research and Development Program of China (2017YFD0101603) and the China’s Agricultural Research System (CARS-15-06).
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2020.00964/full#supplementary-material
TABLE S1 | Nomenclature and sequencing of NHX genes in G. hirsutum, G. barbadense, G. raimondii G. arboreum, and P. trichocarpa.
TABLE S2 | Detailed physicochemical characteristics of NHX proteins of G. hirsutum.
TABLE S3 | Phosphorylation sites of GhNHXs.
TABLE S4 | Duplicate type of NHXs in G. hirsutum and G. barbadense.
TABLE S5 | The Ka/Ks ratios of duplicated gene pairs.
TABLE S6 | Analysis of cis-acting elements in NHX gene promoter regions.
TABLE S7 | Eight Arabidopsis ATNHX sequences.
TABLE S8 | Primer pairs used in the qRT-PCR analysis.
Footnotes
- ^ http://gsds.gao-lab.org/
- ^ https://www.ncbi.nlm.nih.gov/genome/?term=Populus+trichocarpa
- ^ http://www.cbs.dtu.dk/services/TMHMM/
- ^ https://itol.embl.de/index.shtml
- ^ http://meme-suite.org/tools/meme
- ^ https://cottonfgd.org/
References
Aharon, G. S., Apse, M. P., Duan, S., Hua, X., and Blumwald, E. (2003). Characterization of a family of vacuolar Na+/H+ antiporters in Arabidopsis thaliana. Plant Soil 253, 245–256. doi: 10.1023/a:1024577205697
Bailey, T. L., Boden, M., Buske, F. A., Frith, M., Grant, C. E., Clementi, L., et al. (2009). MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 37, W202–W208.
Bassil, E., Tajima, H., Liang, Y., Ohto, M., Ushijima, K., Nakano, R., et al. (2011). 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
Bjellqvist, B., Hughes, G. J., Pasquali, C., Paquet, N., Ravier, F., Sanchez, J., et al. (1993). The focusing positions of polypeptides in immobilized pH gradients can be predicted from their amino acid sequences. Electrophoresis 14, 1023–1031. doi: 10.1002/elps.11501401163
Blom, N., Sicheritzponten, T., Gupta, R., Gammeltoft, S., and Brunak, S. (2004). Prediction of post-translational glycosylation and phosphorylation of proteins from the amino acid sequence. Proteomics 4, 1633–1649. doi: 10.1002/pmic.200300771
Brett, C. L., Donowitz, M., and Rao, R. (2005). Evolutionary origins of eukaryotic sodium/proton exchangers. Am. J. Physiol. Cell Physiol. 288, C223–C239.
Cannon, S. B., Mitra, A., Baumgarten, A., Young, N. D., and May, G. (2004). The roles of segmental and tandem gene duplication in the evolution of large gene families inArabidopsis thaliana. BMC Plant Biol. 4:10. doi: 10.1186/1471-2229-4-10
Cao, B., Xia, Z., Liu, C., Fan, W., Zhang, S., Liu, Q., et al. (2020). New insights into the structure-function relationship of the endosomal-Type Na(+), K(+)/H(+) antiporter NHX6 from mulberry (Morus notabilis). Int. J. Mol. Sci. 21:428. doi: 10.3390/ijms21020428
Chanroj, S., Wang, G., Venema, K., Zhang, M. W., Delwiche, C. F., and Sze, H. (2012). Conserved and diversified gene families of monovalent cation/h(+) antiporters from algae to flowering plants. Front. Plant Sci. 3:25. doi: 10.3389/fpls.2012.00025
Chen, C., Xia, R., Chen, H., and He, Y. (2018). TBtools, a Toolkit for Biologists integrating various HTS-data handling tools with a user-friendly interface. bioRxiv [Preprint], doi: 10.1101/289660
Chen, X., Lu, X., Shu, N., Wang, D., Wang, S., Wang, J., et al. (2017). GhSOS1, a plasma membrane Na+/H+ antiporter gene from upland cotton, enhances salt tolerance in transgenic Arabidopsis thaliana. PLoS One 12:e0181450. doi: 10.1371/journal.pone.0181450
Deinlein, U., Stephan, A. B., Horie, T., Luo, W., Xu, G., and Schroeder, J. I. (2014). Plant salt-tolerance mechanisms. Trends Plant Ence 19, 371–379.
Díaz, J., Ten Have, A., and Van Kan, J. (2002). The role of ethylene and wound signaling in resistance of tomato to Botrytis cinerea. Plant Physiol. 129, 1341–1351. doi: 10.1104/pp.001453
Ding, X., Li, J., Pan, Y., Zhang, Y., Ni, L., Wang, Y., et al. (2018). Genome-wide identification and expression analysis of the UGlcAE gene family in tomato. Intern. J. Mol. Sci. 19:1583. doi: 10.3390/ijms19061583
Dong, Y., Hu, G., Yu, J., Thu, S. W., Grover, C. E., Zhu, S., et al. (2019). Salt-tolerance diversity in diploid and polyploid cotton (Gossypium) species. Plant J. 101, 1135–1151. doi: 10.1111/tpj.14580
Dragwidge, J., Scholl, S., Schumacher, K., and Gendall, A. (2019). NHX-type Na(K)/H antiporters are required for TGN/EE trafficking and endosomal ion homeostasis in Arabidopsis thaliana. J. Cell. Sci. 132:jcs226472. doi: 10.1242/jcs.226472
Dragwidge, J. M., Ford, B. A., Ashnest, J. R., Das, P., and Gendall, A. R. (2018). Two endosomal NHX-Type Na+/H+ antiporters are involved in auxin-mediated development in Arabidopsis thaliana. Plant Cell Physiol. 59, 1660–1669. doi: 10.1093/pcp/pcy090
Elsawy, H., Mekawy, A., Elhity, M., Abdel-Dayem, S., Abdelaziz, M., Assaha, D., et al. (2018). Differential responses of two Egyptian barley (Hordeum vulgare L.) cultivars to salt stress. Plant Physiol. Biochem. 127, 425–435. doi: 10.1016/j.plaphy.2018.04.012
Foster, R., Izawa, T., and Chua, N. H. (1994). Plant bZIP proteins gather at ACGT elements. FASEB J. 8, 192–200. doi: 10.1096/fasebj.8.2.8119490
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
Gao, W., Xu, F. C., Guo, D. D., Zhao, J. R., Liu, J., Guo, Y. W., et al. (2018). Calcium-dependent protein kinases in cotton: insights into early plant responses to salt stress. BMC Plant Biol. 18:15. doi: 10.1186/1471-2229-4-15
Gilmartin, P. M., Memelink, J., Hiratsuka, K., Kay, S. A., and Chua, N. H. (1992). Characterization of a gene encoding a DNA binding protein with specificity for a light-responsive element. Plant Cell 4, 839–849. doi: 10.1105/tpc.4.7.839
He, C., Yan, J., Shen, G., Fu, L., Holaday, A. S., Auld, D., et al. (2005). Expression of an Arabidopsis vacuolar sodium/proton antiporter gene in cotton improves photosynthetic performance under salt conditions and increases fiber yield in the field. Plant Cell Physiol. 46, 1848–1854. doi: 10.1093/pcp/pci201
Heng, Z., Bing, H., Tai, W., Sixue, C., Haiying, L., Yuhong, Z., et al. (2012). Mechanisms of plant salt response: insights from proteomics. J. Proteome Res. 11, 49–67. doi: 10.1021/pr200861w
Hu, B., Jin, J., Guo, A., Zhang, H., Luo, J., and Gao, G. (2015). GSDS 2.0: an upgraded gene feature visualization server. Bioinformatics 31, 1296–1297. doi: 10.1093/bioinformatics/btu817
Jegadeeson, V., Kumari, K., Pulipati, S., Parida, A., and Venkataraman, G. (2019). Expression of wild rice Porteresia coarctata PcNHX1 antiporter gene (PcNHX1) in tobacco controlled by PcNHX1 promoter (PcNHX1p) confers Na+-specific hypocotyl elongation and stem-specific Na+ accumulation in transgenic tobacco. Plant Physiol. Biochem. 139, 161–170. doi: 10.1016/j.plaphy.2019.03.014
Jiang, X., Leidi, E. O., and Pardo, J. M. (2010). How do vacuolar NHX exchangers function in plant salt tolerance? Plant Signal. Behav. 5, 792–795. doi: 10.4161/psb.5.7.11767
Kim, D., Pertea, G., Trapnell, C., Pimentel, H., Kelley, R., and Salzberg, S. L. (2013). TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 14, 1–13.
Krogh, A., Larsson, B., Von Heijne, G., and Sonnhammer, E. L. L. (2001). Predicting transmembrane protein topology with a hidden markov model: application to complete genomes. J. Mol. Biol. 305, 567–580. doi: 10.1006/jmbi.2000.4315
Kumar, S., Kalita, A., Srivastava, R., and Sahoo, L. (2017). NHX1Co-expression of Arabidopsis and improves the tolerance to salinity, oxidative stress, and herbicide in transgenic mungbean. Front. Plant Sci. 8:1896. doi: 10.3389/fpls.2012.01896
Kumar, S., Stecher, G., and Tamura, K. (2016). MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 33, 1870–1874. doi: 10.1093/molbev/msw054
Lamesch, P., Berardini, T. Z., Li, D., Swarbreck, D., Wilks, C., Sasidharan, R., et al. (2012). The Arabidopsis information resource (TAIR): improved gene annotation and new tools. Nucleic Acids Res. 40, 1202–1210.
Lescot, M., Dehais, P., Thijs, G., Marchal, K., Moreau, Y., De Peer, Y. V., et al. (2002). PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 30, 325–327. doi: 10.1093/nar/30.1.325
Letunic, I., and Bork, P. (2016). Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 44, W242–W245.
Letunic, I., Doerks, T., and Bork, P. (2015). SMART: recent updates, new developments and status in 2015. Nucleic Acids Res. 43, D257–D260.
Li, F., Fan, G., Lu, C., Xiao, G., Zou, C., Kohel, R., et al. (2015). Genome sequence of cultivated Upland cotton (Gossypium hirsutum TM-1) provides insights into genome evolution. Nat. Biotechnol. 33, 524–530.
Li, N., Wang, X., Ma, B., Du, C., Zheng, L., and Wang, Y. (2017). Expression of a Na + /H + antiporter RtNHX1 from a recretohalophyte Reaumuria trigyna improved salt tolerance of transgenic Arabidopsis thaliana. J. Plant Physiol. 218, 109–120. doi: 10.1016/j.jplph.2017.07.015
Li, W., Wang, D., Jin, T., Chang, Q., Yin, D., Xu, S., et al. (2011). The Vacuolar Na + /H + antiporter gene SsNHX1 from the halophyte salsola soda confers salt tolerance in transgenic alfalfa (Medicago sativa L.). Plant Mol. Biol. Rep. 29, 278–290. doi: 10.1007/s11105-010-0224-y
Liu, X., Rockett, K., Kørner, C., and Pajerowska-Mukhtar, K. (2015). Salicylic acid signalling: new insights and prospects at a quarter-century milestone. Essays Biochem. 58, 101–113. doi: 10.1042/bse0580101
Lu, W., Guo, C., Li, X., Duan, W., Ma, C., Zhao, M., et al. (2014). Overexpression of TaNHX3, a vacuolar Na(+)/H(+) antiporter gene in wheat, enhances salt stress tolerance in tobacco by improving related physiological processes. Plant Physiol. Biochem. 76, 17–28. doi: 10.1016/j.plaphy.2013.12.013
Mishra, S., Shukla, A., Upadhyay, S., Sanchita, Sharma, P., Singh, S., et al. (2014). Identification, occurrence, and validation of DRE and ABRE Cis-regulatory motifs in the promoter regions of genes of Arabidopsis thaliana. J. Integrat. Plant Biol. 56, 388–399. doi: 10.1111/jipb.12149
Munns, R. (2002). Comparative physiology of salt and water stress. Plant Cell Environ. 25:239. doi: 10.1046/j.0016-8025.2001.00808.x
Munns, R., and Tester, M. (2008). Mechanisms salinity tolerance. Annu. Rev. Plant Biol. 59, 651–681.
Mushke, R., Yarra, R., and Kirti, P. (2019). Improved salinity tolerance and growth performance in transgenic sunflower plants via ectopic expression of a wheat antiporter gene (TaNHX2). Mol. Biol. Rep. 46, 5941–5953. doi: 10.1007/s11033-019-05028-7
Eduardo, O. L., Verónica, B., Lourdes, R., Abdelaziz, E.-H. M., Teresa, R., Beatriz, C., et al. (2009). The AtNHX1 exchanger mediates potassium compartmentation in vacuoles of transgenic tomato. Plant J. 61, 495–506. doi: 10.1111/j.1365-313x.2009.04073.x
Ohta, M., Hayashi, Y., Nakashima, A., Hamada, A., Tanaka, A., Nakamura, T., et al. (2002). Introduction of a Na + /H + antiporter gene from Atriplex gmelini confers salt tolerance to rice. FEBS Lett. 532:282.
Paterson, A. H., Wendel, J. F., Gundlach, H., Guo, H., Jenkins, J., Jin, D., et al. (2012). Repeated polyploidization of Gossypium genomes and the evolution of spinnable cotton fibres. Nature 492, 423–427.
Pehlivan, N., Sun, L., Jarrett, P., Yang, X., Mishra, N., Chen, L., et al. (2016). Co-overexpressing a plasma membrane 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
Potter, S. C., Luciani, A., Eddy, S. R., Park, Y., Lopez, R., and Finn, R. D. (2018). HMMER web server: 2018 update. Nucleic Acids Res. 46, W200–W204.
Qanmber, G., Lu, L., Liu, Z., Yu, D., Zhou, K., Huo, P., et al. (2019). Genome-wide identification of GhAAI genes reveals that GhAAI66 triggers a phase transition to induce early flowering. J. Exp. Bot. 70, 4721–4736. doi: 10.1093/jxb/erz239
Ramsey, J., and Schemske, D. W. (1998). Pathways, mechanisms, and rates of polyploid formation in flowering plants. Annu. Rev. Ecol. Evol. Syst. 29, 467–501. doi: 10.1146/annurev.ecolsys.29.1.467
Ren, Z., Yu, D., Yang, Z., Li, C., Qanmber, G., Li, Y., et al. (2017). Gossypium hirsutum genome-wide identification of the MIKC-Type MADS-box gene family in L. Unravels their roles in flowering. Front. Plant Sci. 8:384. doi: 10.3389/fpls.2012.00384
Rengasamy, P. (2006). World salinization with emphasis on Australia. J. Exper. Bot. 57, 1017–1023. doi: 10.1093/jxb/erj108
Roberto, A., Gaxiola, R. R., Amir, S., Paula, G., and Seth, L. A. (1999). The Arabidopsis thaliana proton transporters, AtNhx1 and Avp1, can function in cation detoxification in yeast. Proc. Natl. Acad. Sci. U.S.A. 1999, 1480–1485. doi: 10.1073/pnas.96.4.1480
Shi, H., Ishitani, M., Kim, C., and Zhu, J. (2000). The Arabidopsis thaliana salt tolerance gene SOS1 encodes a putative Na+/H+ antiporter. Proc. Natl. Acad. Sci. U.S.A. 97, 6896–6901. doi: 10.1073/pnas.120170197
Shu, W., Zhou, H., Jiang, C., Zhao, S., Wang, L., Li, Q., et al. (2019). The auxin receptor TIR1 homolog (PagFBL 1) regulates adventitious rooting through interactions with Aux/IAA28 in Populus. Plant Biotechnol. J. 17, 338–349. doi: 10.1111/pbi.12980
Sonnhammer, E. L. L., Von Heijne, G., and Krogh, A. (1998). A hidden markov model for predicting transmembrane helices in protein sequences. Proc. Int. Conf. Intell. Syst. Mol. Biol. 6, 175–182.
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/1471-2229-4-469
Sze, H., and Chanroj, S. (2018). Plant endomembrane dynamics: studies of K(+)/H(+) Antiporters provide insights on the effects of pH and Ion homeostasis. Plant Physiol. 177, 875–895. doi: 10.1104/pp.18.00142
Tian, F., Chang, E., Li, Y., Sun, P., Hu, J., and Zhang, J. (2017). Expression and integrated network analyses revealed functional divergence of NHX-type Na + /H + exchanger genes in poplar. Sci. Rep. 7, 2607–2607.
Tomoaki, H., and Julian, I. S. (2004). Sodium transporters in plants, diverse genes and physiological functions. Plant Physiol. 136, 2457–2462. doi: 10.1104/pp.104.046664
Tran, L. S., Nakashima, K., Sakuma, Y., Simpson, S. D., 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
Trapnell, C., Williams, B. A., Pertea, G., Mortazavi, A., Kwan, G., Van Baren, M. J., et al. (2010). Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 28, 511–515. doi: 10.1038/nbt.1621
Ulmasov, T., Murfett, J., Hagen, G., and Guilfoyle, T. J. (1997). Aux/IAA proteins repress expression of reporter genes containing natural and highly active synthetic auxin response elements. Plant Cell 9, 1963–1971. doi: 10.1105/tpc.9.11.1963
Verma, D., Lakhanpal, N., and Singh, K. (2019). Genome-wide identification and characterization of abiotic-stress responsive SOD (superoxide dismutase) gene family in Brassica juncea and B. rapa. BMC Genomics 20:227. doi: 10.1186/1471-2229-4-227
Wang, J., and Huang, R. (2019). Modulation of ethylene and ascorbic acid on reactive oxygen species scavenging in plant salt response. Front. Plant Sci. 10:319. doi: 10.3389/fpls.2012.00319
Wang, Y., Tang, H., Debarry, J., Tan, X., Li, J., Wang, X., et al. (2012). MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 40:e49. doi: 10.1093/nar/gkr1293
Wendel, J. F., and Cronn, R. C. (2003). Polyploidy and the evolutionary history of cotton. Adv. Agron. 78, 139–186. doi: 10.1016/s0065-2113(02)78004-8
Wu, A., Hao, P., Wei, H., Sun, H., Cheng, S., Chen, P., et al. (2019). Genome-Wide identification and characterization of glycosyltransferase family 47 in cotton. Front. Genet. 10:824. doi: 10.3389/fpls.2012.00824
Wu, G. Q., Wang, J. L., and Li, S. J. (2019). Genome-wide identification of Na(+)/H(+) antiporter (NHX) genes in sugar beet (Beta vulgaris L.) and their regulated expression under salt stress. Genes 10:401. doi: 10.3390/genes10050401
Wu, G., Wang, Q., Bao, A., and Wang, S. (2011). Amiloride reduces sodium transport and accumulation in the succulent Xerophyte Zygophyllum xanthoxylum under salt conditions. Biol. Trace Elem. Res. 139, 356–367. doi: 10.1007/s12011-010-8662-9
Wu, H. (2018). Plant salt tolerance and Na + sensing and transport. Crop J. 6, 215–225. doi: 10.1016/j.cj.2018.01.003
Wu, Y. Y., Chen, Q. J., Chen, M., Chen, J., and Wang, X. C. (2005). Salt-tolerant transgenic perennial ryegrass (Lolium perenne L.) obtained by Agrobacterium tumefaciens -mediated transformation of the vacuolar Na + /H + antiporter gene. Plant Sci. 169:73.
Xu, Z., Raza, Q., Xu, L., He, X., Huang, Y., Yi, J., et al. (2018). ArabidopsisGmWRKY49, a salt-responsive nuclear protein, improved root length and governed better salinity tolerance in transgenic. Front. Plant Sci. 9:809. doi: 10.3389/fpls.2012.00809
Yamaguchi, T., Apse, M. P., Shi, H., and Blumwald, E. (2003). Topological analysis of a plant vacuolar Na+/H+ antiporter reveals a luminal C terminus that regulates antiporter cation selectivity. Proc. Natl. Acad. Sci. U.S.A. 100, 12510–12515. doi: 10.1073/pnas.2034966100
Yang, Z., and Nielsen, R. (2000). Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models. Mol. Biol. Evolut. 17, 32–43. doi: 10.1093/oxfordjournals.molbev.a026236
Yin, G., Xu, H., Xiao, S., Qin, Y., Li, Y., Yan, Y., et al. (2013). The large soybean (Glycine max) WRKY TF family expanded by segmental duplication events and subsequent divergent selection among subgroups. BMC Plant Biol. 13:148. doi: 10.1186/1471-2229-13-148
Zhang, H. X., Hodson, J. N., Williams, J. P., and Blumwald, E. (2001). Engineering salt-tolerant Brassica plants: characterization of yield and seed oil quality in transgenic plants with increased vacuolar sodium accumulation. Proc. Natl. Acad. Sci. U.S.A. 98, 12832–12836. doi: 10.1073/pnas.231476498
Zhang, T., Hu, Y., Jiang, W., Fang, L., Guan, X., Chen, J., et al. (2015). Sequencing of allotetraploid cotton (Gossypium hirsutum L. acc. TM-1) provides a resource for fiber improvement. Nat. Biotechnol. 33, 531–537.
Zhang, W., Liu, S., Li, C., Zhang, P., and Zhang, P. (2019). Transcriptome sequencing of antarctic moss under salt stress emphasizes the important roles of the ROS-scavenging system. Gene 696, 122–134. doi: 10.1016/j.gene.2019.02.037
Zhao, M., Yin, L., Ma, J., Zheng, J., Wang, Y., Lan, J., et al. (2019). GmERF135 the roles of in improving salt tolerance and decreasing ABA sensitivity in soybean. Front. Plant Sci. 10:940. doi: 10.3389/fpls.2012.00940
Zhu, T., Liang, C., Meng, Z., Sun, G., Meng, Z., Guo, S., et al. (2017). CottonFGD: an integrated functional genomics database for cotton. BMC Plant Biol. 17:101. doi: 10.1186/s12870-017-1039-x
Keywords: Gossypium, Na+/H+ antiporter protein, gene expansion, salt tolerance, stress response
Citation: Fu X, Lu Z, Wei H, Zhang J, Yang X, Wu A, Ma L, Kang M, Lu J, Wang H and Yu S (2020) Genome-Wide Identification and Expression Analysis of the NHX (Sodium/Hydrogen Antiporter) Gene Family in Cotton. Front. Genet. 11:964. doi: 10.3389/fgene.2020.00964
Received: 03 April 2020; Accepted: 31 July 2020;
Published: 18 August 2020.
Edited by:
Maximiller Dal-Bianco, Universidade Federal de Viçosa, BrazilReviewed by:
Shang Haihong, Zhengzhou University, ChinaHai Zhou, South China Agricultural University, China
Copyright © 2020 Fu, Lu, Wei, Zhang, Yang, Wu, Ma, Kang, Lu, Wang and Yu. 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: Hantao Wang, w.wanghantao@163.com; Shuxun Yu, ysx195311@163.com
†These authors have contributed equally to this work