Skip to main content

ORIGINAL RESEARCH article

Front. Genet., 22 February 2022
Sec. Evolutionary and Population Genetics
This article is part of the Research Topic Genomic Selection and Characterization in Cereals View all 20 articles

NHX-Type Na+/H+ Antiporter Gene Expression Under Different Salt Levels and Allelic Diversity of HvNHX in Wild and Cultivated Barleys

  • 1Department of Biosciences, COMSATS University Islamabad (CUI), Islamabad, Pakistan
  • 2Department of Agronomy, Institute of Crop Science, College of Agriculture and Biotechnology, Zhejiang University, Hangzhou, China

Salinity tolerance is a multifaceted trait attributed to various mechanisms. Wild barley is highly specialized to grow under severe environmental conditions of Tibet and is well-known for its diverse germplasm with high tolerance to abiotic stresses. The present study focused on determining the profile of the expression of isoforms of the HvNHX gene in 36 wild and two cultivated barley under salt stress. Our findings revealed that in leaves and roots, expression of HvNHX1 and HvNHX3 in XZ16 and CM72 was upregulated at all times as compared with sensitive ones. The HvNHX2 and HvNHX4 isoforms were also induced by salt stress, although not to the same extent as HvNHX1 and HvNHX3. Gene expression analysis revealed that HvNHX1 and HvNHX3 are the candidate genes that could have the function of regulators of ions by sequestration of Na+ in the vacuole. HvNHX1 and HvNHX3 showed a wide range of sequence variations in an amplicon, identified via single-nucleotide polymorphisms (SNPs). Evaluation of the sequencing data of 38 barley genotypes, including Tibetan wild and cultivated varieties, showed polymorphisms, including SNPs, and small insertion and deletion (INDEL) sites in the targeted genes HvNHX1 and HvNHX3. Comprehensive analysis of the results revealed that Tibetan wild barley has distinctive alleles of HvNHX1 and HvNHX3 which confer tolerance to salinity. Furthermore, less sodium accumulation was observed in the root of XZ16 than the other genotypes as visualized by CoroNa-Green, a sodium-specific fluorophore. XZ16 is the tolerant genotype, showing least reduction of root and leaf dry weight under moderate (150 mM) and severe (300 mM) NaCl stress. Evaluation of genetic variation and identification of salt tolerance mechanism in wild barley could be promoting approaches to unravel the novel alleles involved in salinity tolerance.

Introduction

Na+/H+, counter-transporters (NHXs) not only serve as essential membrane transporters but also help catalyze the neutral exchange of K+ or Na+ for H+ and have been found to play a significant role in pH and ion homeostasis, cell expansion, and salt tolerance. To develop crops with enhanced tolerance to abiotic stresses, the establishment of a better understanding of the underlying mechanisms is essential. With growing advances made within the last decade, researchers have revealed various mechanisms of adaptation and molecular details of responses triggered by salt stress. Although barley is tolerant to high salt stress compared to other cereals, its production is still hampered by salinity (Qui et al., 2011). Salinity is one of the most detrimental stress factors among the abiotic stresses (Stepien and Klobus, 2005) which affected nearly 800 and 32 million hectares of land around the world (Wani et al., 2020). Soil salinization causes severe reduction in barley production all around the world (Rengasamy et al., 2003; Mishra and Tanna 2017). In such a scenario, to go for soil amendments is near to impossible. To cope with soil salinity, the development of salt-tolerant cultivars could be the best possible solution. Hence, understanding salt-tolerant mechanisms is essential for the genetic improvement of crops. Having unique features among cereal crops, barley is widely used in physiological and genetic studies to unravel the mechanisms of salt tolerance (Munns and Mark, 2008). However, the process of domestication of the present cultivated barley (H. vulgare L.) resulted in loss of vital allelic variations (Russell et al., 2004), leading to limited genetic diversity in its gene pools in comparison to their wild ancestors. Therefore, a narrow genetic base could be a great obstacle in the development of cultivars more adapted to the environment (Ellis et al., 2000; Feuillet et al., 2008; Shavrukov et al., 2010). A wider range of genetic variations has been reported in barley populations under various stressful environments (Nevo et al., 1997). Moreover, Tibetan wild barley (H. spontaneum L.), which habited the Qinghai–Tibet Plateau of China years ago, is one of the ancestors of cultivated barley (Dai et al., 2012).

The involvement of both osmotic and ionic components has been elaborated in plant growth inhibition under salt stress (Munns et al., 2006), for example, Na+ specifies not only damage to leaf tissues but also causes many metabolic problems in plants due to its high accumulation in shoots (Munns et al., 2006). Low cytoplasmic Na+/K+ ratio is maintained by multiple mechanisms that include morphological and biochemical adaptations (Munns and Mark, 2008). Vacuolar Na+/H+ antiporters belonging to the NHX gene family, that is, AtNHX1, which serve to detoxify the cytoplasm by compartmentalization of Na+ into the vacuole, have emerged as a vital group of transporters that helped in Na+ tissue tolerance mechanisms (Pardo et al., 2006). According to previous studies, overexpression of the NHX gene and its homologs from other plant species could result in salt tolerance (Zhang and Blumwald 2001; He et al., 2005).

Based on subcellular localization and physiological roles, previous studies classified Na+ (K+)/H+ antiporters in plants into two families, namely, plasma membrane transporter (SOS) and intracellular transporter (NHX) (Brett et al., 2005). (Brett et al., 2005). The NHX family can be further divided into two distinct groups, class-I and class-II (Pardo et al., 2006). Class-I includes NHX isoforms of A. thaliana (AtNHX1-4) with strong vacuolar localization (Pardo et al., 2006; Barragan et al., 2012). However, NHX proteins belonging to class-II, namely, AtNHX5-6 exhibit endosomal localization in cells (Pardo et al., 2006; Rodríguez-Rosales et al., 2008). Cytosol is the main site of action for salts, affecting plant growth and development by disturbing important physiological and biochemical processes. Ion homeostasis, specific for each cellular compartment, is necessary for plant cells to ensure the availability of optimal conditions for gene expression, various enzymatic processes, and protein structure and function (Paroutis et al., 2004). Moreover, a range of physiological processes that include cell expansion, osmotic adjustment, ion regulation, pH homeostasis, and cellular stress responses are essential for plant cell homeostasis and are regulated by NHX-type (cation/H+) antiporters (Pardo et al., 2006; Bassil et al., 2012). These counter-transporters emerged in the early evolution of plants and have been found in entire sequences of plant genomes (Bassil et al., 2012; Chanroj et al., 2012).

The NHX family found in barley consists of four isoforms, mainly localized in the vacuole, called Na+/H+ antiporter genes (HvNHX1 (AB089197.1), HvNHX2 (AY247791), HvNHX3 (DQ372061.1), and HvNHX4 (DQ314285) (Francisco et al., 2012). Expression of genes responsive to salt stress could be species-dependent, and it starts from significantly lower salt levels, for example, 50–100 mM NaCl, which is sufficient for most of the plant species. However, depending on the degree of salt sensitivity, even lower concentrations could be considered for higher salt-sensitive plants. On the contrary, halophyte species could sustain their growth at higher salt levels (Shavrukov., 2013). Plants have developed various mechanisms to get rid of higher concentrations of Na+ ions, which include its transport from the cytosol into the vacuole or out of the cell with the help of Na+/H+ exchanger machinery found in the vacuolar and plasma membranes, respectively (Apse and Blumwald 2007). Thus, enhanced efficiency of the vacuole to compartmentalize more Na+ is a promising strategy to overcome both Na+ toxicity and osmotic effect caused by high salinity (Vera-Estrella et al., 2005).

However, this mechanism could not be the sole remedy of the problem as vacuoles have a limitation for accommodating Na+. Therefore, a combination of more mechanisms needs to work as a unit of the salt tolerance strategy. For instance, limiting Na+ entry and increasing Na+ extrusion could effectively reduce Na+ accumulation in the cytosol. As salt tolerance, barely, is a complex trait controlled by many factors (Zhu, 2001), only one trait is not likely to result in any significant improvement. The salinity problem can only be solved successfully if several important traits are combined in a complementary manner.

Tibet provides a rich gene pool of wild barley with wider variation in adaptation to abiotic stresses, including drought and salinity (Wu et al., 2011; Dai et al., 2012). For the discovery of novel alleles that are potentially related to salt tolerance in wild barley, evaluation of the genetic diversity and identification of salt tolerance mechanisms has been elaborated as important approaches (Wu et al., 2011). Moreover, among various available DNA molecular markers, SNPs are the most abundant type and have been proven to be useful in genetic studies. Barley seeds have been used in genotyping studies since long and also act as important resources for the characterization of genetic variation, which could ultimately pave the way for developing cultivated varieties with enhanced tolerance to abiotic stresses and other production challenges (Nevo and Chen, 2010). SNP markers were initially developed for cultivated barley; however, its usage could also be expanded to wild resources (Rostoks et al., 2005). In the present research, we aim to 1) highlight the salt tolerance mechanism of four isoforms, (HvNHX1 to HvNHX4), located in the vacuole and differing in their level of K+ and Na+ accumulation and 2) determine allelic diversity of HvNHX isoforms in wild and cultivated barleys.

Materials and Methods

This experiment was conducted in two parts. First, seeds of two cultivated barleys, CM72 (tolerant) Gairdner (sensitive), and two wild barleys, XZ16 (tolerant) and XZ169 (sensitive), were used to evaluate the salinity tolerance. Second, 36 wild barley accessions were used from eight sub-populations of Tibetan wild barley accessions to investigate the allelic variation on the basis of single-nucleotide polymorphism (SNP). The seeds were surface-sterilized with 3% H2O2 for 15 min and washed five times with double distilled water. The seeds were transferred to germination boxes and incubated (22/18°C, day/night) for 10 days. Uniform seedlings at the two-leaf stage were transplanted onto 35 L rectangular pots. Nutrient solution was prepared according to the method given by Zahra et al., (2014). The pH of the hydroponic solution was maintained at 5.5–6. After 2 weeks post transplantation, barley plants were grown in three different levels of salt stress, (1) 0 (control); (2) 150 mM NaCl; and (3) 300 mM NaCl. For gene expression analysis, root and leaf tissues were collected at four different time points between 0 and 24 h under 150 and 300 mM salt stress from four genotypes. A complete randomized design with three replications was used in this experiment.

Determination of Sodium and Potassium Concentration in Plant Tissue

The plants were harvested after 3 weeks of salt treatment. Samples of roots and leaves were collected from plants of each treatment separately, and fresh weight was recorded. Similarly, dry weight was recorded by drying the plant samples at 70°C for 74 h respectively. For measurements of root and leaf sample Na+ and K+ concentrations, 0.1 g of dry root and leaf samples was dried to ash and then dissolved in 10 ml HNO3: H2O (1:1). An atomic absorption spectroscope (Shimadzu, Japan) was used for ion content measurements.

RNA to cDNA Synthesis for RT- PCR

RNA from roots and leaves of four genotypes was extracted using a kit method (Tiangen Technology Co., Ltd.DP432), according to the manufacturer’s protocol. Full-length strand cDNA was synthesized from 2 μg of RNA using the Takara Bio Inc., (RR037A). BIO-RAD Master Mix kit was used to carry out RT-PCR reactions according to the manufacturer’s protocol. Measurements were taken for two biological and three technical repeats. The relative gene expression levels were calculated by subtracting the threshold cycle (Ct) values for Gapdh from those of the target gene (to give ΔCt) and then calculating 2−△△Ct (Livak and Schmittgen, 2001). The primers used for RT-PCR were designed using Primer-BLAST (Supplementary Table S1).

Isolation and Sequence Analysis of HvNHX1 and HvNHX3 cDNA

After confirmation of the gene expression analysis, out of four isoforms, two (HvNHX1 and HvNHX3) were selected to further investigate allelic variation based on SNPs and small INDEL detection approach. To detect SNPs in the cDNA pool from 38 barley genotypes, four primer pairs on to HvNHX1 (Supplementary Figure S1; Supplementary Table S2) and three pairs of primers to the HvNHX3 gene were designed to amplify whole coding sequence (CDS) regions for sequencing (Supplementary Figure S2; Supplementary Table S2). PCR primers were designed using Primer-BLAST. Each 50 μL PCR was carried out using 25 μl (2x Easy Taq® PCR Super Mix), 2 μl (cDNA template), 1/1 μl (Primer forward/reverse 10 μM), and 21 μl (ddH2O). To confirm the primer amplification specificity, samples were tested for gel electrophoresis, and the required bandwidth obtained was then matched with a wide range of DNA markers. The samples were then sent to Shanghai Majorbio Co. Ltd., Shanghai, P.R. China for sequencing. Gene-specific primers were used for SNP identification. The sequences obtained from the company were aligned using ClustalX software to observe SNPs.

Visualization of Na+ Ions Through Fluorescence Dye

To observe tissue-specific Na ion accumulation, the roots were stained with 25 mM specific fluorescent probe (Coro Na-Green AM) and 0.02% pluronic acid (Invitrogen) for 3 h as described by in our previous study by Zahra et al., (2014). The roots were incubated in the dark for 3 h and then root tips were carefully washed with deionized water and observed using a confocal microscope (LSM 710 NLO Jena, Germany) at wavelengths of 492 and 516 nm.

Statistical Analysis

SPSS (17.0) software was used for data analyses. The sequences were assembled using DNAStar. TASSEL was used to identify single-nucleotide polymorphisms (SNPs) within the sequence of the HvNXH1 and HvNXH3 genes.

Results

The Difference in Tissue Dry Weight and Ion Content Among Four Genotypes Under Moderate and Severe Salt Stress

After 3 weeks of salt treatment, the treated plants showed a significant reduction in root and leaf dry weight; however, the level of reduction varied among the four genotypes (Table 1). XZ16 and CM72 are more tolerant genotypes, showing a minimal reduction of root and leaf dry weight under moderate (150 mM) and severe (300 mM) NaCl stress, whereas the two salt-sensitive genotypes, Gairdner and XZ169, exhibited a significanty reduction in root and leaf dry weight reduced relative to the control. However, the extent of reduction in dry weight of roots and leaves at both salinity levels was in order of XZ16 < CM72 < XZ169 < Gairdner (Table 1).

TABLE 1
www.frontiersin.org

TABLE 1. Comparison of root and leaf dry weight, Na+ and K+ content, and the Na+/K+ ratios among four barley genotypes under salt stress.

Na+ concentrations in both roots and leaves of four genotypes were noticeably enhanced under salt stress, while K+ concentrations in roots and leaves of the four genotypes showed a remarkable reduction with increasing NaCl levels. Hence, more increase of Na+ concentration and reduction of K+ concentration was observed in Gairdner and XZ169 than that in XZ16 and CM72. Therefore, Gairdner and XZ169 had a considerably higher Na+/K+ ratios than XZ16 and CM72 in both roots and leaves under moderate and severe salt stress (Table 1).

Correlation Between Relative dry Weight and Ionic Contents

Na+ concentrations in leaves and roots were significantly negatively correlated with relative root and leaf dry weight (Figures 1C,D). Moreover, K+ concentration in both the tissues had a positive correlation with relative leaf and root dry weight (Figures 1A,B). However, a significantly negative correlation was observed between shoot and root Na+/K+ ratio and relative leaf and root dry weight (Figures 1E,F).

FIGURE 1
www.frontiersin.org

FIGURE 1. Correlation between K+ and Na+ contents and Na+/K+ ratios and relative leaf (A,C,E) and root (B,D,F) dry weight is based in of four genotypes.

Expression of Na+/H+-Antiporters (HvNHX) Isoform in Response to Moderate and Severe Salt Stress in Roots

We studied the influence of salt stress on HvNHX1, HvNHX2, HvNHX3, and HvNHX4 gene expression in roots and leaves of four barley genotypes after 0, 6, 12, and 24 h of exposure to moderate and severe salt stress. Our results revealed that in the case of XZ16, the HvNHX1 gene showed a higher expression level in the root tissue at all time points under moderate and severe salt stress than the control. The expression level of the HvNHX1 gene in CM72 was also upregulated at all time points under moderate and severe salt stress; however, the expression level of HvNHX1 in XZ16 was markedly enhanced as compared with CM72 (Figures 2A,B), whereas in salt-sensitive genotypes XZ169 and Gairdner, downregulation was observed in the expression level of the HvNHX1 gene at all time points in XZ169, except at 6 h, whereas in Gairdner at both the salinity levels, downregulation was observed as compared with the control (Figures 2A,B).

FIGURE 2
www.frontiersin.org

FIGURE 2. Relative gene expression of HvNHX isoforms in leaves of four genotypes after moderate and severe salt stress. Data are expressed asmeans ± SD of at least three repeats.

The expression level of HvNHX2 was significantly upregulated in XZ16 at both salinity levels as compared with the control, but a higher expression level was observed at 6 and 24 h (Figures 2C,D) However, the expression level of the HvNHX2 gene in CM72 under moderate salt stress was slightly expressed as compared with the control, but at 12 and 24 h exposure of severe stress, the expression level in CM72 was more than that of XZ16 under 300 mM salt stress. In contrast, the expression level of HvNHX2 in Gairdner was downregulated under both salinity levels, while the expression level in XZ169 was not affected at the initial hour but then declined under moderate salt stress, while the reverse was true for severe salt stress (Figures 2C,D).

It was also found that the expression of the HvNHX3 gene in XZ16 was highly upregulated at both salinity levels as compared with other three genotypes, whereas the expression level in CM72 was upregulated in moderate salt stress, while at severe salt stress, slight increase was observed as compared with the control (Figures 2E,F). Moreover, the expression level in XZ169 and Gairdner was downregulated at moderate salt stress, and no change was observed in HvNHX3 gene expression in response to 300 mM NaCl (Figures 2E,F).

The expression level of HvNHX4 was upregulated in XZ16 and CM72 under both levels of salt stress, but the expression level was higher in XZ16 than CM72 under moderate salt stress, whereas the reverse was true for CM72 under severe salt stress. In the roots of XZ169 and Gairdner, the expression of HvNHX4 was downregulated at moderate salt stress, whereas it was slightly upregulated under severe salt stress (Figures 2G,H).

Expression of Na+/H+-Antiporter (HvNHX) Isoform in Response to Moderate and Severe Salt Stress in Leaves

The expression level of HvNHX1 in leaves of XZ16 and CM72 was upregulated at all time points, but markedly increased in XZ16 compared with CM72 under moderate and severe salt stress (Figures 3A,B). In contrast, the expression level of HvNHX1 in Gairdner was downregulated at 24 h of moderate salt stress, while no change was observed at severe salt stress compared with the control. However, the expression level in XZ169 was slightly upregulated under both salinity levels (Figures 3A,B Figures 3A,B). Moreover, the expression level of HvNHX2 was significantly upregulated in XZ16 and CM72 at both salinity levels. However, this gene was more highly expressed in CM72 than XZ16 under both salinity levels, except at 24 h of severe salt stress (Figures 3C,D). In contrast, HvNHX2 in Gairdner was downregulated under both salinity levels, whereas the expression level in XZ169 was slightly increased at 24 h of moderate salt stress as compared with the control (Figures 3C,D). As shown in Figures 3E,F, the expression of the HvNHX3 gene in XZ16 was highly upregulated at both salinity levels as compared with the other three genotypes, except at 24 h of severe salt stress, where expression of HvNHX3 was slightly downregulated as compared with CM72. In contrast, the expression of HvNHX3 in Gairdner was downregulated under both salinity levels. Moreover, the expression level in XZ169 remained unchanged under moderate salt stress but slightly induced at 6 and 12 h and then declined at 24 h of severe salt stress as compared with the control (Figures 3E,F).

FIGURE 3
www.frontiersin.org

FIGURE 3. Relative gene expression of HvNHX isoforms in roots of four genotypes after moderate and severe salt stress. Data are expressed as means ± SD of at least three repeats.

The expression of HvNHX4 was upregulated in XZ16 and CM72 under moderate and severe salt stress as compared with the control. The expression was more pronounced in CM72 than XZ16, except at 24 h of severe salt stress in ZX16. In leaves, the overall expression level of HvNHX4 in XZ169 remained unchanged but slightly induced at 6 h under moderate salt stress, whereas under severe salt stress, the expression of HvNHX4 in XZ169 was slightly upregulated at all time points (Figures 3G,H).

In general, both moderate and severe salt treatment caused fluctuations in the expression of all HvNHX isoforms in roots and leaves of both the salt-tolerant genotypes as compared with sensitive ones. The strongest expression level in XZ16 roots and leaves was ranked as HvNHX1>HvNHX3>HvNHX4>HvNHX2 and in CM72 as HvNHX2>HvNHX4> HvNHX1>HvNHX3. HvNHX isoforms were also induced in XZ169, the salt-sensitive genotype, but to a lesser extent. However, all the HvNHX isoforms were down regulated in Gairdner.

It was also observed that less Na+ was accumulated in the roots of XZ16 than the other genotypes visualized by CoroNa-Green, a sodium-specific fluorophore (Figure 4). Bright images showed that more Na+ content was accumulated in plant roots.

FIGURE 4
www.frontiersin.org

FIGURE 4. Fluorescence detection of Na+ accumulation in roots of four genotypes after moderate and severe salt stress. Scale bar10 µm.

SNP Detection and Validation

In the present study, 36 wild barley accessions (from eight sub-population of Tibetan wild barley accession) and two cultivated barley genotypes were used to investigate the allelic variation in HvNHX1 and HvNHX3 based on SNPs. The sequence of genes HvNHX1 (AB089197.1) and HvNHX3 (DQ372061.1) was downloaded from NCBI. Complete coding regions of mRNA of HvNHX1 (2564 bp) and HvNHX3 (1794 bp) were cloned. Four primer pairs to HvNHX1 and three pairs of primers to the HvNHX3 gene were designed to amplify the complete CDS regions of candidate genes for sequencing. In order to confirm the primer amplification specificity, few of the representative samples were tested for gel electrophoresis, and the required bandwidth obtained was then matched with a wide range of DNA markers (Figures 5, 6). SNPs were detected using the sequencing and alignment method. We successfully amplified and sequenced HvNHX1 and HvNHX3 genes. Details of the nucleotide variations among the barley accessions and/or genotypes for HvNHX1 genes are presented in Table 2 and HvNHX3 genes are presented in Table 3.

FIGURE 5
www.frontiersin.org

FIGURE 5. Pre-sequencing sample amplification trial to verify the gene and amplicon size for the selected primer pairs of HvNHX1 primer 1 (A), primer 2 (B), primer 3(C), and primer 4(D). M = wide range DNA marker (100–2000 bp).

FIGURE 6
www.frontiersin.org

FIGURE 6. Pre-sequencing sample amplification trial to verify the gene and amplicon size for the selected primer pairs of HvNHX3 primer 1 (A), primer 2 (B), and primer 3(C). M = wide range DNA marker (100–2000 bp).

TABLE 2
www.frontiersin.org

TABLE 2. Sequencing results showing single-nucleotide polymorphisms (SNPs) and insertion–deletions (INDEL) pattern in the HvNHX1 gene locus having Na+/H+ activity in barley; del: deletion. Dots indicate the same nucleotide with the reference sequence; the letters in each sample represent nucleotide substitution sites. SPS = SNP per sample. The letters in each sample represent nucleotide substitution sites.

TABLE 3
www.frontiersin.org

TABLE 3. Sequencing results showing single-nucleotide polymorphisms (SNPs) and insertion–deletion (INDEL) pattern in the HvNHX3 gene locus having Na+/H+ activity in barley. H: haplotype; del: deletion. Dots indicate the same nucleotide with the reference sequence; the letters in each sample represent nucleotide substitution sites. SPS = SNP per sample. The letters in each sample represent nucleotide substitution sites.

Overall, evaluation of the sequencing data of 38 barley genotypes showed mutation, including SNPs and small insertion–deletion (INDEL) sites, for targeted gene HvNHX1. Thirty-nine SNPs were observed in total 38 barley genotypes/accessions as shown in Table 2. Moreover, the range of polymorphic sites was 1–3, with an average of one SNP site per barley genotype (Table 2). The HvNHX3 gene was amplified by three primers. The range of polymorphic sites was 1–5, with an average of two SNP sites per barley genotype. Eighty-four SNPs were detected for the HvNHX3 gene in 38 barley genotypes/accessions as shown in Table 3.

Moreover, SNP analysis revealed that in CM72 (salt-tolerant genotype), A is deleted at 300bp and replaced by G at 1978bp while in XZ16 (salt-tolerant genotype), C was replaced by T at 1843bp. For salt-sensitive Gairdner genotype, G, A, and C were deleted at 267, 1,069, and 1,816 bp, respectively. In wild XZ169 salt-sensitive genotype, C was deleted at 1,816 bp and G was replaced by A at 1,820 bp as compared with the reference HvNHX1 gene. Moreover, in the HvNHX3 gene in CM72, A and T was replaced by T and G at 373 and 1221bp, respectively, while in XZ16, G was deleted at 11bp and replaced by T at 45bp while A, T, and G were replaced by T, G, and A at 373, 1,221, and 1,661 bp, respectively (Table 3). For salt-sensitive Gairdner genotype, A was replaced by T at 373bp. In XZ169, G was replaced by T at 11 bp while A was replaced by T at 373 bp and at 1190 bp. T was deleted as compared with the reference HvNHX3 gene. It may be concluded that variations of SNP in salt-tolerant wild barley might offer elite alleles for the development of salt-tolerant barley.

Discussion

Salinity tolerance in plants is a complex multigenic trait, including physiological and molecular aspects (Flowers, 2004). In salt stress conditions, Na+ first enters the cytosol of a plant and disturbs important physiological, biochemical, and molecular processes; consequently, it restricts plant growth, disruption of ion homeostasis (Mane et al., 2010; Basu et al., 2020), and development, therefore posing a serious threat to crop production (Zhu, 2007). One of the key mechanisms in plants to cope with salinity stress is their ability to reduce sodium ion (Na+) transport at both the tissue and cellular level, either by emitting Na+ into tissues (Tester and Davenport 2003), or by maintaining ion homeostasis within the cell, including Na+ compartmentalization in the vacuole (Tester and Davenport 2003; Flowers and Colmer 2008). Barley is a well-known salt-tolerant crop (Steppuhn et al., 2005; Jabeen et al., 2021); however, in cultivated barley, the increase in sensitivity to salinity stress is observed due to the increasingly narrow genetic diversity (Zhu, 2001). However, the Tibet plateau provides a rich pool of wild barleys with a high degree of contrast in salt tolerance that could be attributed to high genetic diversity compared with cultivated barley (Shavrukov et al., 2010).

In this study, the salt tolerance mechanism of four isoforms (HvNHX1 to HvNHX4) and their allelic diversity in wild and cultivated barleys was evaluated after application of 150 and 300 mM salt stress in hydroponic conditions. Salt stress causes several physiological, morphological, and biochemical changes in plants (Uçarlı and Gürel 2020). The current results indicate that XZ16 and CM72 are more tolerant genotypes with minimal reduction of root and leaf dry weight under moderate (150 mM) and severe (300 mM) NaCl stress. However, within tolerant genotypes, wild barley XZ16 showed comparatively less reduction than cultivated barley CM72. This suggested the involvement of a rather different mechanism of salt tolerance in XZ16 compared with CM72. The result was consistent with the findings observed in our previous study (Wu et al., 2011; Zahra et al., 2014). A high expression level of Na+/H+ antiporters in tolerant plants could lead to enhanced Na+ compartmentalization into vacuoles and ultimately improve plant growth by defending the cytoplasm from harmful effects of Na+ (Blumwald et al., 2000; Galvez et al., 2012). Moreover, in the present study, the increase in Na+ and decrease in K+ concentrations in Gairdner and XZ169 was distinctly higher than those in XZ16 and CM72. A similar trend was observed in the studies conducted by Qui et al., 2011; Wu et al., 2011, describing that the reduction in growth was caused by enhanced Na+ and reduced K+ tissue content, which caused ion toxicity and damaged plant metabolism and growth (Table 1).The tissue-specific Na ion accumulation was also confirmed through fluorescence dye, which was directly proportional to the Na+ ion accumulation in roots (Figure 4). Genotypes with the lowest Na+ tissue accumulation produced more biomass and vice versa (Munns and Mark, 2008). In the present study, the expression level of HvNHX isoforms follows a complex pattern, but the gene expression was more induced in salt tolerant genotypes, indicating the important role of these genes in the salt tolerance mechanism. In leaves and roots, the expression of HvNHX1 and HvNHX4 in XZ16 and CM72 was upregulated at all time points as compared with sensitive ones, NHX1 (Quintero et al., 2009; Fukuda et al., 2004; Saqib et al., 2005; Brini et al., 2007) and NHX4 (Gálvez et al., 2012) in wheat and Arabidopsis, and a high expression pattern of NHX1 (Quintero et al., 2009; Fukuda et al., 2004; Saqib et al., 2005; Brini et al., 2007) and NHX4 genes in tomato (Galvez et al., 2012) was reported to be involved it in better plant growth because NHX1 and NHX4 are the key molecular players in maintaining plant cell homeostasis by regulating several physiological processes, such as cell expansion, osmotic adjustment, cell volume, pH, and ion regulation (Pardo et al., 2006).

Na+/H+ antiporter protein has different membrane positions in the cell, and its function may be affected by ion accumulation. Previously, the topological studies elucidated that the position of N-terminal of AtNHX1 is facing toward the cytosol, and its C-terminal, hydrophilic region, residing in the vacuolar lumen could protect the cytoplasm from deleterious effects of Na+ (Yamaguchi and Blumwald, 2005, Flowers and Colmer, 2008. The function of Na+/H+ antiporters may not only be related to the regulation of gene expression but could also be involved in transcriptional modification of the proteins. The activity of antiporters could be regulated by phosphorylation through the interaction of various kinases with other cellular proteins. So, the differential response of these binding factors in a species-dependent manner could alter the activity of the Na+/H+ antiporters. For instance, the binding of AtCaM15, a calmodulin-like protein 15 (localized in plant vacuolar compartment), to the C-terminal domain of AtNHX1 (a tonoplast transporter) changed the Na+/K+ selectivity of the antiporter in Arabidopsis (Yamaguchi and Blumwald, 2005; Munns and Mark, 2008). Moreover, similar to the overexpression of HvNHX1 in tolerant genotypes, HvNHX3 was also expressed in wild barley XZ16 and CM72. These results suggest that HvNHX1 and HvNHX3 may encode a putative vacuolar NHX that could play an important role in salt tolerance by mediating K+/H+ exchange in plants (Liu et al., 2010).

In plants, H+V-PPase and H+-ATPase are the two different vacuolar pumps that help the Na+/H+ antiporter in the vacuolar lumen to transport Na+ from the cytoplasm to vacuoles by generating the electrochemical gradient force. Generally, under salinity stress, salt-tolerant plants maintain higher K+/Na+ or lower Na+/K+ ratios in the cytoplasm and regulate the osmotic balance of the cells by sequestering Na+ in the vacuoles (Maeshima, 2000). Our results are consistent with those of a previous study that a lower Na+/K+ ratio was observed in tolerant genotypes (Table 1) (Wu et al., 2011; Zahra et al., 2014). The HvNHX2 and HvNHX4 isoforms were also induced by salt stress, although not to the same extent as HvNHX1 and HvNHX3. Expression of HvNHX2 and HvNHX4 in roots under moderate salt stress was more upregulated in XZ16 than in CM72, whereas high expression levels of HvNHX2 and HvNHX4 under severe salt stress in both roots and leaves were observed in CM72. The results from previous studies showed that HvNHX2 and HvNHX4 play an important role in the maintenance of K+ concentration in plant tissues (Venema et al., 2003; Rodriguez Rosales et al., 2008).

To study the polymorphism in terms of SNP detection in the genes of interest is a powerful technique to investigate the gene(s) function and get desirable mutations for crop breeding. In our previous studies, more allelic variation in Tibetan wild barley accessions was observed for HvCBF3, HvCBF4, and HvHKT genes responsible for salt tolerance (Qui et al., 2011; Wu et al., 2011). These findings are further confirmed by the present study that Tibetan wild barley could provide rich source of allelic variation for the salt-responsive gene(s) as compared with salt-tolerant cultivated barley CM72. So, evaluation of genetic variation and identification of salt tolerance mechanism in wild barley are important steps to unravel the novel alleles involved in salinity tolerance. In conclusion, physiological and gene expression analysis revealed that HvNHX1 and HvNHX3 are the candidate genes that function as regulating ions by sequestration of Na+ in the vacuole. Moreover, Tibetan wild barley could be used as a rich source of genetic variation to explore the dynamics of abiotic stress tolerance in barley and other cereal crops.

Data Availability Statement

The authors acknowledge that the data presented in this study must be deposited and made publicly available in an acceptable repository, prior to publication. Frontiers cannot accept a article that does not adhere to our open data policies.

Author Contributions

ZJ: Contributed to conception and design of study and manuscript writing; FI: performed the data analysis; NH: contributed to sampling, data analysis and manuscript writing; YH: performed the data analysis; GZ: contributed to experimental design, manuscript writing, and editing.

Funding

This research was funded by National Natural Science Foundation of China (31330055, 31771687).

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.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.809988/full#supplementary-material

References

Apse, M. P., and Blumwald, E. (2007). Na+ Transport in Plants. FEBS Lett. 581, 2247–2254. doi:10.1016/j.febslet.2007.04.014

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Bassil, E., Coku, A., and Blumwald, E. (2012). Cellular Ion Homeostasis: Emerging Roles of Intracellular NHX Na+/H+ Antiporters in Plant Growth and Development. J. Exp. Bot. 63, 5727–5740. doi:10.1093/jxb/ers250

CrossRef Full Text | Google Scholar

Basu, S., Kumar, A., Benazir, I., and Kumar, G. (2020). Reassessing the Role of Ion Homeostasis for Improving Salinity Tolerance in Crop Plants. Physiologia Plantarum 171, 502–519. doi:10.1111/ppl.13112

PubMed Abstract | CrossRef Full Text | Google Scholar

Brett, C. L., Donowitz, M., and Rao, R. (20052004). Evolutionary Origins of Eukaryotic Sodium/proton Exchangers. Am. J. Physiol. Cel Physiol. 288, C223–C239. doi:10.1152/ajpcell.0036010.1152/ajpcell.00360.2004

CrossRef Full Text | Google Scholar

Brini, F., Hanin, M., Mezghani, I., Berkowitz, G. A., and Masmoudi, K. (2007). Overexpression of Wheat Na+/H+ Antiporter NHX1 and H- Pyrophosphatase TVP1 Improve Salt- and Drought-Stress Tolerance in Arabidopsis thaliana Plants. J. Exp. Bot. 58, 301–308.

Google Scholar

Dai, F., Nevo, E., Wu, D., Comadran, J., Zhou, M., Qiu, L., et al. (2012). Tibet Is One of the Centers of Domestication of Cultivated Barley. Proc. Natl. Acad. Sci. 109, 16969–16973. doi:10.1073/pnas.1215265109

PubMed Abstract | CrossRef Full Text | Google Scholar

Ellis, R. P., Forster, B. P., Robinson, D., Handley, L. L., Gordon, D. C., Russell, J. R., et al. (2000). Wild Barley: a Source of Genes for Crop Improvement in the 21st century. J. Exp. Bot. 51342, 9–17. doi:10.1093/jexbot/51.342.9

CrossRef Full Text | Google Scholar

Feuillet, C., Langridge, P., and Waugh, R. (2008). Cereal Breeding Takes a Walk on the Wild Side. Trends Genet. 24, 24–32. doi:10.1016/j.tig.2007.11.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Flowers, T. J., and Colmer, T. D. (2008). Salinity Tolerance in Halophytes*. New Phytol. 179, 945–963. doi:10.1111/j.1469-8137.2008.02531.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Flowers, T. J. (2004). Improving Crop Salt Tolerance. J. Exp. Bot. 55, 307–319. doi:10.1093/jxb/erh003

CrossRef Full Text | Google Scholar

Fukuda, A., Nakamura, A., Tagiri, A., Tanaka, H., Miyao, A., Hirochika, H., et al. (2004). Function, Intracellular Localization and the Importance in Salt Tolerance of a Vacuolar Na+/H+ Antiporter from Rice. Plant Cel Physiol 45, 146–159. doi:10.1093/pcp/pch014

CrossRef Full Text | Google Scholar

Gálvez, F. J., Baghour, M., Hao, G., Cagnac, O., Rodríguez-Rosales, M. P., and Venema, K. (2012). Expression of LeNHX Isoforms in Response to Salt Stress in Salt Sensitive and Salt Tolerant Tomato Species. Plant Physiol. Biochem. 51, 109–115. doi:10.1016/j.plaphy.2011.10.012

PubMed Abstract | CrossRef Full Text | Google Scholar

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 Fieldfiber Yield in the Field. Plant Cel Physiol 46, 1848–1854. doi:10.1093/pcp/pci201

CrossRef Full Text | Google Scholar

Jabeen, Z., Fayyaz, H. A., Irshad, F., Hussain, N., Hassan, M. N., Li, J., et al. (2021). Sodium Nitroprusside Application Improves Morphological and Physiological Attributes of Soybean (Glycine max L.) Under Salinity Stress. PLoS ONE 16 (4), e0248207. doi:10.1371/journal.pone.0248207

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, H., Tang, R., Zhang, Y., Wang, C., Lv, Q., Gao, X., et al. (2010). AtNHX3 Is a Vacuolar K+/H+ Antiporter Required for Low-Potassium Tolerance in Arabidopsis thaliana. Plant Cel Environ 33, 1989–1999. doi:10.1111/j.1365-3040.2010.02200.x

CrossRef Full Text | Google Scholar

Livak, K. J., and Schmittgen, T. D. (2001). Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2−ΔΔCT Method. Methods 25 (4), 402–408. doi:10.1006/meth.2001.1262

PubMed Abstract | CrossRef Full Text | Google Scholar

Maeshima, M. (2000). Vacuolar H+-pyrophosphatase. Biochim. Biophys. Acta (Bba) - Biomembranes 1465, 37–51. doi:10.1016/s0005-2736(00)00130-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Mane, A. V., Karadge, B. A., and Samant, J. S. (2010). Salinity Induced Changes in Photosynthetic Pigments and Polyphenols of Cymbopogon Nardus (L.) Rendle. J. Chem. Pharm. Res. 2, 338–347.

Google Scholar

Mishra, A., and Tanna, B. (2017). Halophytes: Potential Resources for Salt Stress Tolerance Genes and Promoters. Front. Plant Sci. 8, 829. doi:10.3389/fpls.2017.00829

PubMed Abstract | CrossRef Full Text | Google Scholar

Munns, R., James, R. A., and Läuchli, A. (2006). Approaches to Increasing the Salt Tolerance of Wheat and Other Cereals. J. Exp. Bot. 57, 1025–1043. doi:10.1093/jxb/erj100

CrossRef Full Text | Google Scholar

Munns, R., and Tester, M. (2008). Mechanisms of Salinity Tolerance. Annu. Rev. Plant Biol. 59, 651–681. doi:10.1146/annurev.arplant.59.032607.092911

PubMed Abstract | CrossRef Full Text | Google Scholar

Nevo, E., Apelbaum-Elkaher, I., Garty, J., and Beiles, A. (1997). Natural Selection Causes Microscale Allozyme Diversity in Wild Barley and a Lichen at 'Evolution Canyon', Mt. Carmel, Israel. Heredity 78, 373–382. doi:10.1038/hdy.1997.60

CrossRef Full Text | Google Scholar

Nevo, E., and Chen, G. (2010). Drought and Salt Tolerances in Wild Relatives for Wheat and Barley Improvement. Plant Cel. Envir. 33, 670–685. doi:10.1111/j.1365-3040.2009.02107.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Pardo, J. M., Cubero, B., Leidi, E. O., and Quintero, F. J. (2006). Alkali Cation Exchangers: Roles in Cellular Homeostasis and Stress Tolerance. J. Exp. Bot. 57, 1181–1199. doi:10.1093/jxb/erj114

CrossRef Full Text | Google Scholar

Paroutis, P., Touret, N., and Grinstein, S. (2004). The pH of the Secretory Pathway: Measurement, Determinants, and Regulation. Physiology 19, 207–215. doi:10.1152/physiol.00005.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiu, L., Wu, D., Ali, S., Cai, S., Dai, F., Jin, X., et al. (2011). Evaluation of Salinity Tolerance and Analysis of Allelic Function of HvHKT1 and HvHKT2 in Tibetan Wild Barley. Theor. Appl. Genet. 122, 695–703. doi:10.1007/s00122-010-1479-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodriguez-Rosales, M. P., Jiang, X., Galvez, F. J., Aranda, M. N., Cubero, B., and Venema, K. (2008). Overexpression of the Tomato K+/H+ Antiporter LeNHX2 Confers Salt Tolerance by Improving Potassium Compartmentalization. New Phytol. 179, 366e377. doi:10.1111/j.1469-8137.2008.02461.x

CrossRef Full Text | Google Scholar

Rostoks, N., Mudie, S., Cardle, L., Russell, J., Ramsay, L., Booth, A., et al. (2005). Genome-wide SNP Discovery and Linkage Analysis in Barley Based on Genes Responsive to Abiotic Stress. Mol. Genet. Genomics. 274, 515–527. doi:10.1007/s00438-005-0046-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Russell, J., Booth, A., Fuller, J., Harrower, B., Hedley, P., Machray, G., et al. (2004). A Comparison of Sequence-Based Polymorphism and Haplotype Content in Transcribed and Anonymous Regions of the Barley Genome. Genome 47, 389–398. doi:10.1139/g03-125

PubMed Abstract | CrossRef Full Text | Google Scholar

Saqib, M., Zörb, C., Rengel, Z., and Schubert, S. (2005). The Expression of the Endogenous Vacuolar Na+/H+ Antiporters in Roots and Shoots Correlates Positively with the Salt Resistance of Wheat (Triticum aestivum L.). Plant Sci. 169, 959–965. doi:10.1016/j.plantsci.2005.07.001

CrossRef Full Text | Google Scholar

Shavrukov, Y., Gupta, N. K., Miyazaki, J., Baho, M. N., Chalmers, K. J., Tester, M., et al. (2010). HvNax3-a Locus Controlling Shoot Sodium Exclusion Derived from Wild Barley (Hordeum Vulgare Ssp. Spontaneum). Funct. Integr. Genomics 10, 277–291. doi:10.1007/s10142-009-0153-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Shavrukov, Y. (2013). Salt Stress or Salt Shock: Which Genes Are We Studying. Exbotj 64, 119–127. doi:10.1093/jxb/ers316

PubMed Abstract | CrossRef Full Text | Google Scholar

Stepien, P., and Klobus, G. (2005). Antioxidant Defense in the Leaves of C3 and C4 Plants under Salinity Stress. Physiol. Plant 125, 31–40. doi:10.1111/j.1399-3054.2005.00534.x

CrossRef Full Text | Google Scholar

Steppuhn, H. M., Genuchten, T., and Grieve, C. M. (2005). Root-zone Salinity Selecting a Product-Yield index and Response Functions for Crop Tolerance. Crop Sci. 45, 209–220. doi:10.2135/cropsci2005.0209

CrossRef Full Text | Google Scholar

Tester, M., and Davenport, R. (2003). Na+ Tolerance and Na+ Transport in Higher Plants. Ann. Bot. 91, 503–527. doi:10.1093/aob/mcg058

PubMed Abstract | CrossRef Full Text | Google Scholar

Uçarlı, C., and Gürel, F. (2020). Differential Physiological and Molecular Responses of Three-Leaf Stage Barley (Hordeum Vulgare L.) under Salt Stress within Hours. Plant Biotechnol. Rep. 14, 89

Google Scholar

Venema, K., Belver, A., Marin-Manzano, M. C., Rodriguez-Rosales, M. P., and Donaire, J. P. (2003). A Novel Intracellular K+/H+ Antiporter Related to Na+/H+ Antiporters Is Important for K Ion Homeostasis in Plants. J. Biol. Chem. 278, 22453e22459. doi:10.1074/jbc.m210794200

CrossRef Full Text | Google Scholar

Vera-Estrella, R., Barkla, B. J., García-Ramírez, L., and Pantoja, O. (2005). Salt Stress in Thellungiella Halophila Activates Na+ Transport Mechanisms Required for Salinity Tolerance. Plant Physiol. 139, 1507–1517. doi:10.1104/pp.105.067850

PubMed Abstract | CrossRef Full Text | Google Scholar

Wani, S. H., Kumar, V., Khare, T., Guddimalli, R., Parveda, M., Solymosi, K., et al. (2020). Engineering Salinity Tolerance in Plants: Progress and Prospects. Planta 251, 76–29. doi:10.1007/s00425-020-03366-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, D., Qiu, L., Xu, L., Ye, L., ChenSun, M. D., Sun, D., et al. (2011). Genetic Variation of HvCBF Genes and Their Association with Salinity Tolerance in Tibetan Annual Wild Barley. PLoS ONE 6, e22938. doi:10.1371/journal.pone.0022938

PubMed Abstract | CrossRef Full Text | Google Scholar

Yamaguchi, T., and Blumwald, E. (2005). Developing Salt-Tolerant Crop Plants: Challenges and Opportunities. Trends Plant Sci. 10, 615–620. doi:10.1016/j.tplants.2005.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Zahra, J., Nazim, H., Cai, S., Han, Y., Wu, D., Zhang, B., et al. (2014). The Influence of Salinity on Cell Ultrastructures and Photosynthetic Apparatus of Barley Genotypes Differing in Salt Stress Tolerance. Acta Physiol. Plant 36, 1261–1269. doi:10.1007/s11738-014-1506-z

CrossRef Full Text | Google Scholar

Zhang, H.-X., and Blumwald, E. (2001). Transgenic Salt-Tolerant Tomato Plants Accumulate Salt in Foliage but Not in Fruit. Nat. Biotechnol. 19, 765–768. doi:10.1038/90824

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, J.-K. (2001). Plant Salt Tolerance. Trends Plant Sci. 6, 66–71. doi:10.1016/s1360-1385(00)01838-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, J. K. (2007). “Plant Salt Stress,” in Encyclopedia of Life Sciences. Editor J. Wiley. doi:10.1002/9780470015902.a0001300.pub2

CrossRef Full Text | Google Scholar

Keywords: salinity stress, wild barley, gene expression, polymorphism, SNP, vacuole

Citation: Jabeen Z, Irshad F, Hussain N, Han Y and Zhang G (2022) NHX-Type Na+/H+ Antiporter Gene Expression Under Different Salt Levels and Allelic Diversity of HvNHX in Wild and Cultivated Barleys. Front. Genet. 12:809988. doi: 10.3389/fgene.2021.809988

Received: 05 November 2021; Accepted: 28 December 2021;
Published: 22 February 2022.

Edited by:

Muhammad Abdul Rehman Rashid, Government College University, Faisalabad, Pakistan

Reviewed by:

Meixue Zhou, University of Tasmania, Australia
Essa Ali, University of Poonch, Pakistan
Mehnaz Ghaffar, KPK, Pakistan

Copyright © 2022 Jabeen, Irshad, Hussain, Han and Zhang. 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: Zahra Jabeen, zahra.jabeen@comsats.edu.pk; Guoping Zhang, zhanggp@zju.edu.cn

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.