Next Article in Journal
Novel SCN5A p.Val1667Asp Missense Variant Segregation and Characterization in a Family with Severe Brugada Syndrome and Multiple Sudden Deaths
Previous Article in Journal
Procaspase-Activating Compound-1 Synergizes with TRAIL to Induce Apoptosis in Established Granulosa Cell Tumor Cell Line (KGN) and Explanted Patient Granulosa Cell Tumor Cells In Vitro
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 22
Issue 9
10.3390/ijms22094698
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Structure and Expression Analysis of Sucrose Phosphate Synthase, Sucrose Synthase and Invertase Gene Families in Solanum lycopersicum

by
Yaoke Duan
1,
Lan Yang
1,
Haijia Zhu
1,
Jie Zhou
1,
Hao Sun
2,* and
Haijun Gong
1,*
1
Shaanxi Engineering Research Center for Vegetables/College of Horticulture, Northwest A & F University, Yangling 712100, China
2
Henan Key Laboratory of Ion-Beam Bioengineering, College of Agricultural Sciences, Zhengzhou University, Zhengzhou 450000, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(9), 4698; https://doi.org/10.3390/ijms22094698
Submission received: 12 April 2021 / Revised: 26 April 2021 / Accepted: 26 April 2021 / Published: 29 April 2021
(This article belongs to the Section Molecular Plant Sciences)
Browse Figures
Versions Notes

Abstract

:
Sucrose phosphate synthase (SPS), sucrose synthase (SUS) and invertase (INV) are all encoded by multigene families. In tomato (Solanum lycopersicum), a comprehensive analysis of structure characteristics of these family genes is still lacking, and the functions of individual isoforms of these families are mostly unclear under stress. Here, the structure characteristics of the three families in tomato were analyzed; moreover, as a first step toward understanding the functions of isoforms of these proteins under stress, the tissue expression pattern and stress response of these genes were also investigated. The results showed that four SPS genes, six SUS genes and nineteen INV genes were identified in tomato. The subfamily differentiation of SlSPS and SlSUS might have completed before the split of monocotyledons and dicotyledons. The conserved motifs were mostly consistent within each protein family/subfamily. These genes demonstrated differential expressions among family members and tissues, and in response to polyethylene glycerol, NaCl, H2O2, abscisic acid or salicylic acid treatment. Our results suggest that each isoform of these families may have different functions in different tissues and under environmental stimuli. SlSPS1, SlSPS3, SlSUS1, SlSUS3, SlSUS4, SlINVAN5 and SlINVAN7 demonstrated consistent expression responses and may be the major genes responding to exogenous stimuli.

1. Introduction

Plants are autotrophic organisms that can fix carbon dioxide and produce carbohydrates through photosynthesis. Sucrose is the main end product of photosynthesis in higher plants, and it is transported from the source leaves to sink organs [1]. Sucrose plays an important role in plant growth and development. On the one hand, it provides energy and structural components of plants; on the other hand, sucrose and its hydrolysis products may function as signaling molecules, which regulate the expression of genes involved in important physiological processes [2,3]. Sucrose also functions in osmotic adjustment [3] and thus plays an important role in stress adaption.
The key enzyme responsible for sucrose synthesis is sucrose phosphate synthase (SPS), which catalyzes the formation of 6-phosphate sucrose from UDP-glucose and 6-phosphate fructose [1]. The main enzymes involved in sucrose catabolism are invertase (INV) and sucrose synthase (SUS)—invertase catalyzes the degradation of sucrose into glucose and fructose, whereas SUS converts sucrose into UDP-glucose and fructose [4]. Invertases can be classified as acid invertase and alkaline/neutral (A/N) invertase based on their optimum pH [5]. Acid invertases are usually localized in the cell wall or vacuole, and alkaline/neutral invertases are localized in the cytosol, mitochondria or plastid [2]. SUS proteins are mainly localized in the cytosol or plasma membrane, and some are localized in the vacuole, cell wall or mitochondria [1].
SPS, SUS and INV are all encoded by multigene families and have been characterized in some plant species. SPS is encoded by a small family. For example, in the model plants Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa), there are four and five SPS members, respectively [6,7]. Maize (Zea mays) and apple (Malus domestica) both have six members of SPS genes [8,9], whereas Amborella trichopoda only has two SPS genes [10]. According to phylogenetic analysis, SPS genes are categorized into four families (A, B, C and D), and family D only exists in some monocots [6,8]. The member number of SUS genes varies greatly among plant species. In Amborella trichopoda, only two SUS genes have been identified [11]. There are six, eight, twelve and fourteen SUS genes in Arabidopsis, carrot (Daucus carota), soybean (Glycine max) and tobacco (Nicotiana tabacum), respectively [11,12]. In Chinese pear, it is reported that there were thirty SUS genes [13]. SUS genes can be classified into three subfamilies—SUS I, SUS II and SUS III [1]. The genes encoding invertase have been identified in some plant species. For instance, in Arabidopsis, eight acid invertase genes and nine A/N invertase genes have been identified [14]. There are eight acid invertase genes and six A/N invertase genes in sugarcane [15], and these corresponding numbers in soybean are nineteen and thirteen, respectively [16]. In tomato (Solanum lycopersicum), an important vegetable crop, four SPS genes [17], six SUS genes [11], eleven acid invertase genes [5] and eight alkaline/neutral invertase encoding genes [2] have also been identified. However, the structure characteristics of the genes have not been comprehensively analyzed, except those of one type of invertase—A/N invertase [2].
Up to date, studies on the functions of genes encoding SPS, SUS and INV have mostly focused on growth, development and fruit quality, especially the former two aspects. Bahaji et al. suggested that SPS is essential for Arabidopsis viability by using double and triple mutants of SPSA1, SPSA2, SPSB and SPSC [18]. In sugarcane, SoSPS1 overexpression increased the plant height and stalk number of some transgenic lines [19]. Heterologous expression of a spinach SPS gene in cotton improved the fiber quality [20], while the expression of the Arabidopsis family A SPS gene in tobacco increased the fiber length [21]. Park et al. reported that expression of Arabidopsis SPS gene in poplar not only increased the length of xylem fibers, but also delayed the onset of senescence and advanced the bud flush [22]. In addition, SPS genes may also be involved in the regulation of flowering time, flower number, pollen germination and fruit development [23,24,25]. SUS is considered as a biochemical marker of sink strength [11], and it plays roles in the development of vascular tissues and shoot apical meristem. SUS genes are highly expressed in plant vascular tissues and are involved in the synthesis of both cellulose and callose [1]. Barratt et al. observed that Arabidopsis double mutant sus5/sus6 had reduced callose levels in the sieve plates of phloem [26]. Wei et al. found that overexpression of poplar xylem SUS2 in tobacco increased the cellulose content and the thickness of the xylem cell wall [27]. In tomato, Goren et al. reported that suppression of SUS genes did not affect the carbohydrate levels in the fruits, but affected auxin signaling and the leaf morphology, indicating the possible roles of SUS in regulating early leaf development [28]. SUS may also play important roles in mutualism with symbiotic organisms and fruit ripening [29,30]. Invertase is also involved in the regulation of plant growth and development. Knockout mutation in A/N-InvC caused a reduction in shoot growth and oxygen consumption, suggesting that the invertase played an important role in the respiratory process in mitochondria [31]. In rice, the mutation of CYT-INV1 inhibited the root growth and delayed flowering [32]. Zanor et al. found that silencing the expression of LIN5, which encodes a cell wall invertase in tomato, changed the flower and fruit morphology, increased the petal and sepal numbers of flowers, and decreased pollen viability and seed size, further confirming the important role of invertase in reproductive growth of plants [33].
In plants, although the roles of some SPS, SUS and INV genes in growth and development have been investigated, responses of these genes to environmental stresses are largely unknown, and information on their roles under stress conditions is very limited. Takehara et al. reported that high expression of SUS3 in brown rice was associated with its heat tolerance [34]. Wang et al. found that inhibition of CsSUS3 expression decreased the tolerance of cucumber to hypoxic stress [35]. Qian et al. found that overexpression of tea CsINV5 enhanced the cold tolerance of Arabidopsis [36]. In tomato, Liu et al. and Xu et al. reported that manipulating the expression of cell wall invertase inhibitor gene could regulate fruit setting under heat stress and chilling tolerance [37,38]. Since individual members of SPS, SUS and INV family genes in plants have distinct tissue expression patterns and have been found to demonstrate different functions in growth and development [16,26], the functions of the individual isoforms of these enzymes under stress conditions may be different. In this study, comprehensive analyses on the structural characteristics of genes encoding SPS, SUS and INV in tomato were conducted. As a first step toward understanding the functions of isoforms of these enzymes under stress conditions, the tissue expression pattern and the expression response of these genes to different environmental stimuli were also investigated. The study may help us better understand the gene structure characteristics of the main sucrose metabolism related enzymes, and provide a basis for investigating the roles of individual isoforms of each enzyme in growth, development and stress tolerance of tomato plants.

2. Results

2.1. Identification of SPS, SUS and INV Genes in Tomato Genome

To identify the SPS genes in tomato, iterative protein BLAST analysis was performed using previously identified Arabidopsis SPS sequences as the queries. The candidate Solanum lycopersicum SPS (SlSPS) genes were further screened by BLASTp in Swiss-Prot database and Batch-CDD tests. As a result, four SlSPS genes were identified. The genes were renamed based on their chromosome numbering, and the basic characteristics of these genes and the corresponding proteins were analyzed (Table 1). These SlSPS genes had similar CDS lengths, and the corresponding proteins shared similar amino acid sequence lengths, molecular weights and isoelectric points. The online tool Kinasephos was used to predict the phosphorylation sites, and the results showed that the SPS proteins have 15–24 putative phosphorylation sites. All the SPS proteins were predicted to be localized at the plasma membrane.
Using the same methods, including two rounds of BLAST and screening by structure characteristics, six SlSUS genes were identified in the tomato genome (Table 1). The CDS of SlSUS genes had 2412–2676 bp, and the length of corresponding polypeptides was 803–891 Aa. The isoelectric point of these SUS proteins was 5.87–6.94, and the molecular weight range was 91.6–100.7 KDa. The six SUS proteins were all predicted to be localized at the plasma membrane. The SlSUS contained 9–20 phosphorylation sites, with SUS6 possessing the highest number.
Nineteen INV genes were identified in the tomato genome, including 11 acid invertase genes and eight A/N invertase genes (Table 1). According to the predicted subcellular localization of the corresponding proteins, the acid invertase genes can be further divided into cell wall invertase genes (INVCW1–9) and vacuolar invertase genes (INVVR1–2). We named these genes according to the property and subcellular localization of the corresponding proteins, and numbered them based on their localizations on chromosomes. The CDS length of the INVCW genes and the corresponding amino acid sequence and molecular weight of the proteins were similar, except those of INVCW2 and INVCW6, which were slightly shorter or lower. The isoelectric point and number of putative phosphorylation sites varied among the INVCWs, ranging from 4.96 to 9.23, and from 5 to 14, respectively. The two vacuolar invertase genes (INVVR1–2) had similar CDS lengths, and the isoelectric points and molecular weights of the corresponding proteins were also similar. The length of SlINVAN CDS and amino acid sequence varied among family members, with ranges of 1656–2019 bp and 551–672 Aa, respectively. Subcellular localization prediction of these INVAN proteins showed that two (INVAN1/2) were localized in the mitochondria, two (INVAN6/8) in the chloroplast, and others in the cytoplasm.

2.2. Phylogenetic Analysis of the SPS, SUS and INV Families

To analyze the phylogenetic relationship of SPSs, SUSs or INVs from tomato with those from Arabidopsis and rice, which are, respectively, dicotyledonous and monocotyledonous model plants, an unrooted maximum likelihood tree was constructed using the MEGA software. As shown in Figure 1A, the SlSPS had a closer relationship with AtSPS than OsSPS. SPS from these three plant species could be divided into four distinct families: A, B, C and D. SlSPS1 and SlSPS2 were grouped into family A, SlSPS3 was grouped into family B, and SlSPS4 belonged to family C. No SPS from tomato or Arabidopsis was grouped into family D (Figure 1A).
The evolutionary relationship of SUS proteins in tomato, Arabidopsis and rice was analyzed using the amino acid sequences of 19 SUSs from these species (Figure 1B). These SUS proteins can be divided into three groups: SUS I, SUS II and SUS III. SlSUS1, SlSUS3 and SlSUS5 belonged to the SUS I group and were clustered together, and they were close to AtSUS1 and AtSUS4 in Arabidopsis. SlSUS4 belonged to SUS II, and it was closer to a SUS (LOC_Os03g22120) in rice, rather than to that in Arabidopsis. SlSUS6 and SlSUS7 belonged to SUS III, and they were close to AtSUS5 and AtSUS6.
To explore the phylogenetic relationship of acid invertase protein in tomato, Arabidopsis and rice, we constructed a phylogenetic tree based on the amino acid sequences of 11 SlINVs, 8 AtINVs and 10 OsINVs (Figure 1C). The acid invertases were grouped into two clades: 23 INVCWs were predicted to be localized in the cell wall and six INVVRs in the vacuole. SlINVCW3–9 and SlINVCW1–2 were clustered together, respectively. All the acid INVs in tomato were closer to those in Arabidopsis, rather than to those in rice. The phylogenetic relationship of A/N invertases from tomato, Arabidopsis and rice was also analyzed (Figure 1D). According to the phylogenetic tree, these invertases can be classified into two groups—α and β. SlINVAN3, 4, 5, and 7 belonged to the β group, and the remaining four belonged to the α group. The α group could be further divided into α1 and α2 subgroups according to their subcellular localization. In tomato, SlINVAN6 and 8 belonged to the α1 group, whereas SlINVAN1 and 2 belonged to the α2 group.

2.3. Gene Structure and Conserved Protein Motif Analysis of SPS, SUS and INV Gene Families

The gene structure display server (GSDS) online program was used to analyze the intron-exon structure. Analysis on the conserved domain structure was based on the online program of NCBI BATCH CD-search tool (https://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi, accessed on 12 April 2020). The MEME search tool was employed to predict the conserved protein motifs. The results showed that the SlSPS genes had 12–14 exons and 11–13 introns (Figure S1). Glucosyl transferase Glycos-transf-1 domain (pfam00534) and Glyco-trans-4-4 domain (pfam13579) belong to the Glucosyl transferase GTB-type-superfamily (cl10013). All the SlSPS members had Glycos-transf-1 and S6PP domains (cl37722) (Figure S2), and the Glyco-trans-4-4 domain only existed in SlSPS3 (Figure 2A). The Glycos-transf-1 domain is related to the transfer of glucosyl, and S6PP domains (SPP-like domain) may be the site of SPS binding SPP [46,47]. Ten conserved motifs were predicted in each SPS protein, and most of the motifs were repeated once, except motif 6, which was repeated twice only in SlSPS2 (Figure 2A and Figure S3A). Taken the analysis results of the motif and domain together, it was found that motifs 3, 1, 10 and 2 belonged to the Glycos-transf-1 domain, while motif 9 belonged to the S6PP domain. In SPS1 and SPS2, motifs 7, 4, 8 and 6 belonged to the GTB-type-superfamily; in SPS4, motifs 5, 7, 4 and 8 belonged to the GTB-type-superfamily; and in SPS3, motifs 7, 4, 8 and 6 belonged to the Glyco-trans-4-4 domain.
The six SlSUS genes contained 11–15 exons (Figure S1). These genes shared two typical domains: sucrose synthase domain (pfam00862) and glucosyl transferase domain (pfam00534) (Figure 2B and Figure S4). Ten conserved motifs existed in these SUS proteins and each motif was repeated once (Figure 2B and Figure S3B). The motifs 9, 10, 3, 4, 5 and 1 belonged to the sucrose synthase domain, whereas motifs 8, 6, 2 and 7 belonged to glucosyl-transferase domain (Figure 2B).
To explore the structure characteristics of acid invertase genes in tomato, the intron-exon structure and conserved motifs were analyzed. Except INVCW2, which had eight exons, the other INVCW genes all had six exons (Figure S1). Both INVVR1 and INVVR2 contained seven exons (Figure S1). It is interesting to note that the second exon of all acid invertase genes contained only nine nucleotide acids (Figure S1). The N-terminal of acid invertase is glycosylated and it belongs to glycohydrolase family 32 (Glycoside Hydrolase Family 32, GH32) [48]. The cysteine catalytic MWECP/V and β-furosidase motif NDPNG/A were conservative in all the tomato acid invertase sequences (Figure S5). Both INVCWs and INVVRs in tomato had the Glyco-hydro-32N and Glyco-hydro-32C domains, and the DUF3357 domain was only present in INVVRs (Figure 2C). Motif analysis showed that SlINVCWs contained motifs 1–10 except INVCW1, where motif 9 was lacking (Figure 2C and Figure S3C). Another interesting finding is that some motifs were only conserved in SlINVCWs. For instance, motif 9 and motif 10 were present in SlINVCWs but missing in SlINVVR proteins. Motifs 1, 5, 8, 6, 4 and 7 belonged to the Glyco-hydro-32N domain, while motifs 7, 9, 2, 10 and 3 belonged to Glyco-hydro-32C domain (Figure 2C).
There were eight A/N invertase genes in tomato (Table 1). Among these, SlINVAN1, 2, 6 and 8, which encoded α group invertases, contained 6–7 exons, whereas SlINVAN3, 4, 5 and 7, which encoded β group invertases, had four exons (Figure S1). A/N invertases are nonglycosylated proteins, belonging to the glycosyl hydrolase family 100 (GH100) [48]. The Glyco-hydro-100 domain was observed in all the SlINVAN genes (Figure 2D). The eight SlINVAN proteins shared 10 conserved motifs, which had similar arrangement (Figure 2D and Figure S3D). It is noted that the amino acid sequence before the first motif (motif 8) was longer in the α group invertases than that in the β group invertases (Figure 2D).

2.4. Chromosome Distribution, Synteny Analysis of SPS, SUS and INV Genes

Chromosome distribution of SPS, SUS and INV gene families was analyzed base on the physical location of the GCF_000188115.4_SL3.0_genomic database from the NCBI website. Gene duplication plays an important role in the amplification of gene family and subsequent evolution, with segmental and tandem duplications being the main causes of gene family expansion in plants [49]. Duplication events of these genes were analyzed using the MCScanX software. The syntenic relationship of SPS, SUS and INV genes in different plant species was analyzed using the Dual Synteny Plotter software.
Our results showed that the SlSPS genes were distributed on four chromosomes—SPS1, SPS2, SPS3 and SPS4 on chromosomes 7, 8, 9 and 11, respectively—and no tandem duplication or segmental duplication was detected (Figure 3A). Synteny analysis between SPSs in tomato and those in both Arabidopsis and rice indicated that AtSPS1F-SlSPS1, SlSPS1-LOC_Os06g43630 and AtSPS3F-SlSPS3 were, respectively, identified as syntenic gene pairs (Figure 3B).
The six tomato SUS genes were distributed on five chromosomes (chromosomes 2, 3, 7, 9 and 12) (Figure 4A). Both SlSUS3 and SlSUS5 were located on chromosome 7, and they were identified as tandem duplication genes (Figure 4A). SlSUS1 and SlSUS5 were identified as segmental duplication gene pairs. Synteny analysis showed that there were three syntenic gene pairs: AtSUS1-SlSUS1, AtSUS1-SlSUS5 and SlSUS1- LOC_Os07g42490.
There were 19 INV genes in the tomato genome, including 11 acid INV and eight A/N INV genes; these genes were unevenly distributed on eight chromosomes (Figure 5A). For instance, on chromosome 10, there were five INV genes, whereas there was only one INV gene on chromosomes 4 and 8, respectively (Figure 5A). Gene duplication analysis showed that there were three tandem duplication pairs (INVCW3/INVCW4, INVCW5/INVCW6 and INVCW8/INVCW9) and two segmental duplication gene pairs (INVCW3/INVCW5, INVAN6/INVAN8). To further determine the origin and evolution dynamics of tomato INV genes, we investigated the syntenic relationship between INV genes in tomato and those in Arabidopsis thaliana and rice. The result showed that 10 syntenic gene pairs between tomato and Arabidopsis and nine pairs between tomato and rice were identified, and they were listed as follows: AtAN_INVC-SlINVAN1-LOC_Os03g20020, AtAN_INVB-SlINVAN3, AtAN_INVD-SlINVAN5, AtAN_INVF-SlINVAN5, AtAN_INVG-SlINVAN5, AtCWINV2-SlINVCW3-LOC_Os04g33720, AtCWINV4-SlINVCW3-LOC_Os04g33720, AtCWINV2 -SlINVCW5-LOC_Os04g33720, AtCWINV2-SlINVCW5-LOC_Os02g33110, AtCWINV4-SlINVCW5-LOC_Os04g33720, AtCWINV4-SlINVCW5- /LOC_Os02g33110, AtAN_INVE-SlINVAN8-LOC_Os02g32730, AtAN_INVE-SlINVAN8-LOC_Os04g33490, SlINVAN6-LOC_Os02g32730. SlINVAN6-LOC_Os04g33490. The INV genes in tomato which showed syntenic relationships with those in Arabidopsis or rice were mainly located at chromosomes 1, 6, 9, 10 and 11.

2.5. Tissue Expression Profile of SUS, INV and SPS in Tomato

Both bioinformatic and experimental approaches were employed to determine the tissue expression profiles of SPS, SUS and INV genes in tomato. The expression data collected from public RNA-seq repositories demonstrated that each family gene showed differential tissue expression patterns (Figure 6). For instance, the expression of SPS1 was higher than that of any other SPS in the fruit, columella, locular tissue, pericarp, placenta and septum (Figure 6A). In the SUS family, compared with other members, SUS1 and SUS3 were both highly expressed in the fruit, columella, outer pericarp, ovary wall, ovule and pericarp (Figure 6B). The expression of SUS5 was much higher in roots than any other tissues (Figure 6B). Among the INV family genes, by and largely, INVVR1 and INVAN4 (especially the former) were highly expressed in most tissues (Figure 6C). In addition, INVCW4 appeared to be specifically expressed in the flower and pollen (Figure 6C).
A quantitative RT-PCR technique was also introduced to analyze the expression of SUS, INV and SPS genes in the leaf and root of tomato. The results showed that all the SPS genes in tomato were mainly expressed in the leaf, and both SPS1 and SPS3 (especially the latter) demonstrated much higher expression than the other two SPS genes (Figure 7A). Except SUS4, most SUS genes showed higher expression in the root than the leaf, with SUS3 and SUS5 being the top two highly expressed SUS genes (Figure 7B). We also analyzed the tissue expression of all mined SlINVAN and SlINVVR genes, as well as SlINVCW5 and SlINVCW6—two highly expressed SlINVCW genes observed in the collected public data (Figure 6C and Figure 7C). Among the SlINVANs, SlINVAN4 was highly expressed in both leaves and roots (Figure 7C). SlINVAN6 and SlINVAN7 also had relatively high expression in the leaf and root, respectively (Figure 7C). SlINVVR1 demonstrated high expression in the leaf. In brief, by and large, these results are consistent with the public RNA-seq data.

2.6. Expression of SPS, SUS and INV Genes under Simulated Stress Conditions in Tomato

PEG and NaCl are usually used to simulate drought and salinity stress, respectively. H2O2, one of the main reactive oxygen species and signaling molecules, usually accumulates under stress conditions. ABA and SA are also signaling molecules and respond to environmental stresses. In this study, PEG, NaCl, H2O2, ABA, and SA treatments were applied to simulate various stress conditions. The results showed that the expression of SPS1 was largely induced by all the treatments, especially in the root (Figure 8A,B). SPS2 expression was increased in the leaf but decreased in the root by ABA treatment (Figure 8C,D). PEG and NaCl treatments for 72 h decreased SPS2 expression in the leaf and root (Figure 8C,D). The expression of SPS3 was decreased in the leaf by all the treatment at various time points, while in the root, the expression was decreased by PEG and NaCl treatments at 12 h, but it was increased later (Figure 8E,F). ABA treatment also stimulated SPS3 expression in the root at later stress periods (Figure 8F). SPS4 expression in the leaf was slightly decreased at 12 and 24 h by all the five treatments, but it was differentially increased at 72 h, while in the root, the expression was increased by H2O2 and SA treatments at 12 and 24 h, but decreased by PEG and ABA treatments at a later period.
The SUS family members also demonstrated differential expression responses to the chemical treatments. The expression of SUS1, SUS3 and SUS4 in the leaf and root was largely increased in response to the treatment of five chemicals, except at 12 h, when SUS1 and SUS4 expression was decreased by H2O2 and SA treatments, respectively (Figure 9A–F). Under PEG stress, the expression of SUS5 was significantly decreased in the leaf and root (Figure 9G,H). In the leaf, SUS5 expression was increased by H2O2 and SA treatments, whereas the expression was decreased by NaCl and ABA treatments in the root (Figure 9G,H). The expression of SUS6 and SUS7 in the leaf and root was stimulated by NaCl treatment at the early stress period, but they were decreased at the later period (Figure 9I,L). The expression of these two genes was inhibited at a late period of PEG treatment (Figure 9I–L). Neither the expression of SUS6 nor that of SUS7 was affected by SA treatment (Figure 9I–L).
The changes in the expression of all tomato A/N invertase and vacuolar acid invertase genes as well as two major cell wall invertase genes were investigated under the different treatment conditions (Figure 10A–X). The results showed that the expression of INVAN1 in the leaf was significantly decreased by NaCl, H2O2, ABA and SA treatments at early stages, and it was recovered to the control level or even increased after 72 h of treatment (Figure 10A). The expression of INVAN2 and INVAN3 was increased by PEG treatment in the leaf and root (Figure 10C–F). In the leaf, NaCl, H2O2, ABA and SA treatments all decreased the INVAN3 expression at 12 and 24 h, but the expression was unaltered or increased at 72 h (Figure 10E), while in the root, the INVAN3 expression was stimulated by salt stress (Figure 10F). Salt stress also promoted the expression of INVAN4 in the leaf (Figure 10G). The expression of INVAN5 was inhibited in both leaves and roots under all treatment conditions (Figure 10I,J). The expression of INVAN6 and INVAN8 was stimulated in the root by PEG and NaCl treatments at 12 h and they were also increased in the leaf by ABA treatment at 72 h (Figure 10K,L,O,P). The INVAN7 expression was decreased by different chemical treatments, especially in the root (Figure 10M,N). In the leaf, the expression of INVVR1 was promoted only by PEG, H2O2, ABA at 72 h, while in the root, the expression was stimulated by PEG, NaCl and ABA treatments. All the treatments inhibited the expression of INVVR2 and INVCW6 in the leaf and root as well as the INVCW5 expression in the leaf at early stages (Figure 10S,T).

3. Discussion

SPS, SUS and INV are the main enzymes that regulate sugar metabolism. With the development of sequencing technology and bioinformatics, the family genes encoding these enzymes have been identified and the structure characteristics have also been analyzed in some plants, such as Arabidopsis and rice [6,7,11,14]. In tomato, although some genes of SPS, SUS and INV families have been identified [2,5,11,17], a comprehensive analysis of the structure characteristics and evolution of these genes is still lacking. Moreover, the expression responses of these genes to environmental stimuli are largely unknown. In this study, comprehensive analyses on the gene characteristics and evolution of SPS, SUS and INV in tomato were conducted, and the expression responses of these genes to different stress stimuli were also investigated.

3.1. Evolution and Structure of SPS Genes

SPS is the key enzyme responsible for sucrose synthesis [1]. The SPS gene family in plants is relatively small: there are 3–5 SPS genes in most species [50]. In this study, four SPS genes were identified in tomato (Table 1), which is consistent with most plants. The relatively small SPS family in plants may be partly related to the lack of duplication events of this gene, at least in tomato, as observed here (Figure 3A). In this study, synteny analysis demonstrated that AtSPS1F-SlSPS1, SlSPS1-LOC_Os06g43630 and AtSPS3F-SlSPS3 were syntenic gene pairs, respectively (Figure 3B), implying the interspecific conservation of these genes.
The identified four SlSPS members together with the SPS in Arabidopsis and rice were categorized into four families (Figure 1A). SlSPS1 and SlSPS2 were grouped into family A, while SlSPS3 and SlSPS4 belonged to families B and C, respectively (Figure 1A). The presence of family D SPS in rice but absence in tomato and Arabidopsis (Figure 1A) is consistent with the finding of Castleden et al. [8], who suggested that families A, B and C of SPS widely exist in monocotyledons and dicotyledons, whereas family D SPS only exist in the Poaceae (monocotyledons). In addition, the existence of families B and C SPS in both monocotyledon (rice) and dicotyledons (Arabidopsis and tomato) implies that the differentiation of SPS families was completed before the differentiation of monocotyledons and dicotyledons. It has been suggested that the SPS family differentiation event occurred about 200 million years ago [51].
Analysis of intron-exon structure is helpful to study the possible origin and relationship of genes [11]. Here, in tomato, the four SlSPS genes had similar exon-intron structures, and the numbers of both exons and introns were the same in SlSPS1 and SlSPS2 (Table 1; Figure S1). Compared with SlSPS1 and SlSPS2, SlSPS3 and SlSPS4 experienced intron loss or gain events. SlSPS3 had a longer first exon due to the loss of the equivalence of first intron, while the equivalence of the fifth exon of SlSPS4 was inserted by the 238 bp intron, resulting in an additional exon compared with SlSPS1 and SlSPS2 (Figure S1A). Similar intron loss and gain events were also observed in litchi [50], reflecting the evolutionary relationship of these individual SPS genes plants.
All the four SPS members had glycosyl transferase and S6PP domains (Figure 2A), which are unique to SPS family [47]. Ten conserved motifs were predicted in each SPS protein, and most of the motifs were repeated once (Figure 2A and Figure S3A). It is interesting to note that motif 6, which belongs to the GTB-type-superfamily, was repeated at the N-terminal of SlSPS2 (Figure 2A). The functional significance of this additional repeat remains to be investigated.

3.2. Evolution and Structure of SUS Genes

In plants, the number of SUS gene members varies greatly among species. For instance, only two SUS genes were identified in Amborella trichopoda [11], whereas in Chinese pear, thirty SUS genes were reported [13]. The difference in family member numbers may be related to the difference in family gene expansion [52]. In this study, six SUS genes were identified in tomato (Table 1). SlSUS3 and SlSUS5 were identified as tandem duplication genes, while SlSUS1 and SlSUS5 were identified as segmental duplication gene pairs (Figure 4A). Therefore, the SUS duplication events have contributed to the expansion of this gene family.
In plants, SUS can be classified into three subfamilies—SUS I, SUS II and SUS III [1]. Phylogenetic analysis of SlSUS together with the SUS in Arabidopsis and rice demonstrated that in each subfamily, there were SUS from both monocotyledonous (rice) and dicotyledonous (Arabidopsis and tomato) plants (Figure 1B). This suggests that the subfamily differentiation of SUS happened before the split of monocotyledons and dicotyledons.
A previous study has revealed that SUS genes were highly conservative in many dicotyledons and monocotyledons [11]. Consistent with this, our study shows that the exon-intron structure and arrangement of the SUS family genes in tomato were also generally conservative (Figure S1B and Figure 2B). However, there were still some differences in the exon-intron structure. For instance, the equivalences of the 5th and 12th introns in SlSUS4 were lost in all the SUS I genes (SlSUS1, SlSUS3 and SlSUS5), while two SUS III genes (SlSUS6 and SlSUS7) lost the equivalence of the 12th but not the 5th intron (Figure S1B). Intriguingly, compared with other subfamily members, two SUS III subfamily genes, SlSUS6 and SlSUS7, respectively, had one and two more exons at the 3′ end (Figure S1B), and the corresponding amino acid sequences were longer (Table 1). In addition, it is noted that SlSUS1 and SlSUS3, respectively, had a longer UTR in the 5′ end. The significance of the differential intron loss, 3′ end extension and different UTR length still remains to be explored in the future.

3.3. Evolution, Structure and Classification of INV Genes

The number of INV family genes also varies greatly among species. There were 14 INV genes in sugarcane [15], while this number in soybean is 32 [16]. In the present study, nineteen INV genes were identified in tomato, including eleven acid invertase genes and eight A/N invertase genes (Table 1). The interspecific variation in the number of INV genes should be associated with the difference in gene duplication events as found in this study—there were three tandem duplication pairs (INVCW3/INVCW4, INVCW5/INVCW6 and INVCW8/INVCW9) and two segmental duplication gene pairs (INVCW3/INVCW5, INVAN6/INVAN8) in tomato (Figure 5A). In addition, ten syntenic gene pairs between tomato and Arabidopsis and nine pairs between tomato and rice were identified (Figure 5B), suggesting the evolutionary conservation of INV genes in plants.
In this study, eleven acid invertase genes including nine cell wall invertase genes and two vacuolar invertase genes were identified (Table 1). Most SlINVCW genes had six exons, except SlINVCW2, which had eight exons, while the SlINVVR genes contained seven exons (Table 1; Figure S1C). Compared with most SlINVCW genes, SlINVCW2 contained two introns while the SlINVVR genes had one intron in the equivalence of the third exon (Figure S1C). Phylogenetic analysis indicated that SlINVCW1/2 and SlINVVR genes evolved earlier than other SlINVCW genes (Figure 1C and Figure 2C). This suggests that the acid invertase genes might have experienced intron losses during evolution. The second exon of all the 11 identified INV genes was a 9 bp long mini exon, which encoded three amino acids (Figure 2C and Figure S1C). Such a mini exon is a typical structure characteristic of acid invertase genes in plants [53], and it is also the smallest exon that has even been found in plants [54]. It is noted that the first exon was longer in SlINVVR genes than SlINVCW genes (Figure S1C). The DUF3357 domain was encoded by the first exon, and it was present in SlINVVRs but absent in SlINVCWs, resulting in longer N-terminal of the former (Figure 2C).
Eight A/N invertase genes were identified in tomato (Table 1). According to the sequence homology, these invertases could be divided into two group: the α group (SlINVAN1, 2, 6 and 8) and β group (SlINVAN3, 4, 5 and 7) (Figure 1D). The two groups had different exon numbers: 6–7 exons in the α group and four in the β group (Figure S1). Our results are generally consistent with previous findings that the α group and β group genes typically have six and four exons, respectively [55,56]. However, there was an exception in SlINVAN1, which contained seven exons (Figure S1D). It should be pointed out that the exon numbers of SlINVAN1 and SlINVAN6 found in this study were inconsistent with the study of Pan et al. [2]. The differences may be related to the differences in the genomic annotations, which remain to be further explored. The two groups of tomato invertases also differed in the subcellular localization and the length of amino acid sequence. In this study, the α group invertases were localized either in mitochondria (SlINVAN1 and 2) or cytoplasm (excluding cytosol, SlINVAN6 and 8), whereas the β group (SlINVAN3, 4, 5 and 7) were localized in the cytosol (Table 1). The α group proteins had 653–672 amino acids, whereas the β group contained 551–570 amino acids (Table 1). The longer protein sequence in the α group was due to the longer N-terminal region (Figure 2D). Whether there were more signal peptide sequences in the long N-terminal region and, if there were, whether they were related to the multilocalization of the α group invertases (Table 1) remains to be investigated.

3.4. Tissue Expression Pattern of SPS, SUS and INV Genes

Comprehensive expression analysis of all gene family members may help understand their functions. In tomato, the expression profiles of SPS, SUS and INV genes have not been comprehensively analyzed in different tissues, except those of SUS and part of INV genes, but on in a few tissues [2,11]. In this study, in silico analysis on the expression of SPS, SUS and INV genes was conducted in 19 tissues of tomato, and the results showed that there were expression differences among family members and tissues (Figure 6A–C).
Among the four SlSPS genes, SlSPS1 demonstrated relatively higher expression in different tissues except pollen (Figure 6A). In the leaves, the expression of SlSPS1 and SlSPS3 was higher than that of SlSPS2 and SlSPS4, while in the root, the expression of SlSPS1 was the highest among the family members (Figure 6A). These results were confirmed by quantitative PCR in ‘Alisa Craig’ (Figure 7A). Bahaji et al. observed that, compared with the wild type, Arabidopsis mutants spsa1/spsc and spsa1/spsa2/spsc had smaller rosettes [18], flowers and siliques; moreover, the seeds of spsa1/spsb/spsc and spsa1/spsa2/spsb/spsc mutants germinated poorer and the plants were sterile. Therefore, the highly expressed SlSPS, especially SlSPS1 in different tissues, may play key roles in tomato growth and development, such as seed germination, vegetative growth, flower and fruit development.
In this study, it is noted that SlSUS1 demonstrated the highest expression in fruits, and the expression of SlSUS3 was also relatively high (Figure 6B). In the fruit, both SlSUS1 and SlSUS3 had the highest expression in the outer pericarp (Figure 6B). D’Aoust et al. reported that antisense inhibition of TOMSSF (SlSUS1 in this study; Table 1) reduced sucrose import into the tomato fruit and decreased fruit setting [40]. Zhao et al. found that downregulating the expression of FaSS1, a SUS gene in strawberry, delayed fruit ripening [30]. These studies suggest that SlSUS1 and SlSUS3 may play important roles in tomato fruit development. Among the six SlSUS genes, both SlSUS3 and SlSUS5 had the highest expression in the root (Figure 6B and Figure 7B), indicating that these two genes mainly contributed to the strength of this sink tissue. In addition, in the vascular tissue and leaf, SlSUS3 had a higher expression (Figure 6B). The highly expressed SlSUS3 gene in vascular tissues may be involved in the synthesis of both cellulose and callose, as suggested previously in Arabidopsis [26] and poplar [27].
Among the cell wall-localized acid invertase genes, SlINVCW3 demonstrated relatively high expression in different tissues, especially vascular tissue, flower and fruit (Figure 6C). These results are consistent with previous studies [43,57]. SlINVCW3 (named LIN5 in previous studies, Table 1) has been suggested to function in ovary-to-fruit transition in the sieve elements by generating a glucose signal, and thus regulating cell division [57]. Silencing this gene not only affected the fertility and fruit development, but also influenced the fruit hormone level [33]. These studies suggest that SlINVCW3 plays an important role in the reproductive growth of tomato plants. In this study, it was observed that SlINVCW4 was highly expressed in pollen (Figure 6C). The high expression of this gene has been confirmed to be required for pollen development [58]. It is interesting to note that SlINVVR1 was highly expressed in different tissues, especially in parenchyma and fruit (Figure 6C). Qin et al. reported that the inhibitor of VI (SlINVVR1 in this study) could mediate sucrose metabolism and affect fruit ripening, suggesting the role of SlINVVR1 in these processes [17]. The function of the highly expressed SlINVVR1 in parenchyma remains to be investigated. Among the eight SlINVAN genes, SlINVAN4 demonstrated relatively higher expression in different tissues (Figure 6C). Recently, Leskow et al. found that silencing this gene impaired the growth phenotype, delayed flowering and reduced fruit setting [45]. This suggests that the high expression of this gene is required for the normal growth, flowering and fruit development. SlINVAN6 had the highest expression in pistil (Figure 6C), and its role remains to be investigated.
In a word, the highly expressed genes, including SlSPS1, SlSUS1, SlSUS3, SlINVCW3, SlINVVR1 and SlINVAN4 in most parts of the fruit (Figure 6A–C), may be involved in the so called ‘futile cycles’ of sucrose [59], which regulate sugar accumulation and/or sugar signaling, and thus fruit development. Despite the high expression in fruits, these genes also demonstrated differential expressions in different parts of the fruit (Figure 6A–C). Therefore, the exact regulatory functions of these genes in sugar accumulation in fruits remain to be further investigated. Although the functions of some members of these gene families have been clarified in flowers and fruits, those of other members, especially the SlINVAN genes, are mostly unclear. Moreover, the genes of these three families may play important roles in other physiological processes except flower and fruit development. For instance, it has been shown that SlSUS participates in the regulation of early leaf morphology development [28]. However, relevant information in these aspects is still limited, and more work is needed in the future.

3.5. Expressions of SPS, SUS and INV Genes under Stresses

Although the roles of some (not all) SPS, SUS and INV genes in flower and fruit development of tomato have been investigated, the responses of these genes to environmental stimuli and their roles in stress tolerance still remain largely unknown. In this study, different chemicals were applied to tomato seedlings to simulate various stress conditions, and the expression responses of the SPS, SUS and INV family genes were investigated. In general, we observed that within each gene family, the members demonstrated differential expression responses to environmental stimuli and in different tissues (Figure 8, Figure 9 and Figure 10). It is interesting to see that some genes demonstrated consistent expression responses to different treatments in the leaves and/or roots. For example, the SlSPS3 expression in leaves was downregulated and the SlSPS1 expression in roots was upregulated under all treatment conditions (Figure 8B,E). Under PEG, NaCl, H2O2 and ABA treatments, the expression of SlSUS1, SlSUS3 and SlSUS4 was mostly increased in both leaves and roots. In addition, both SlINVAN5 and SlINVAN7 demonstrated downregulated expression under the five treatment conditions (Figure 10I,J,M,N). These results suggest that SlSPS1, SlSPS3, SlSUS1, SlSUS3, SlSUS4, SlINVAN5 and SlINVAN7 may be the major genes responding to the different exogenous stimuli in tomato. It is also noted that the expression of some genes was not altered or changed much by the treatments, such as SlSPS2 and SlINVAN6 in leaves under NaCl treatment, SlINVAN4 in roots under NaCl stress, SlSUS4 and SlINVAN6 in roots, and SlSUS6, SlSUS7 and SlINVAN2 in both leaves and roots under SA treatment. The differential expression responses of these genes to environmental stresses have been reported previously. For instance, Hu et al. observed that in perennial ryegrass, the expression of SPS was slightly increased in the roots but decreased in the stem under salt stress, while in the leaves, the expression was increased in the salt-sensitive accession but remained unchanged in the tolerant accession [60]. Solís-Guzmán et al. reported that the expression of AtSPS2F and AtSPS4F were up-regulated under osmotic stress, whereas those of AtSPS1F and AtSPS3F were down-regulated [61]. Yang et al. found that in apple leaves, drought stress increased the expression of SPS2, SPS3, SPS4, SUSY1, SUSY3, SUSY5, NINV1 and NINV2, but decreased the expression of SPS1, SUSY2, CWINV1, CWINV2, NINV3 and AINV1. They also observed that the expression changes in SPS5, SPS6, SUSY4 and AINV3 expression were dependent on stress duration—increased in the early stress stage but decreased in the later stage [3]. Dahro et al. reported that PtrA/NINV expression in Poncirus trifoliata was upregulated by low temperature, salt, dehydration, sucrose and ABA, but downregulated by glucose. These studies suggest that each isoform of these gene families may have distinct functions under different environmental stimuli and in different tissues [62].
It is noticed that under the exogenous treatments, the changes in the expression of SlSUS genes were different from those of SlINV genes. In most cases, the expression of SlSUS genes was upregulated by different treatments (Figure 9), and the upregulation magnitudes were higher than those of SlINV genes (Figure S6). This seems to suggest that the SlSUS genes may play more important roles than the SlINV genes in the response to environmental stimuli. In future, more work is needed to investigate the roles of individual genes of the three families in stress tolerance in tomato plants.

4. Materials and Methods

4.1. Identification of SPS, SUS and INV Genes in Tomato

To identify the SPS genes in tomato (Solanum lycopersicum), the annotated tomato genome was downloaded from the website of National Centre for Biotechnology Information (http://www.ncbi.nlm.nih.gov, accessed on 26 February 2020). The sequences of all known AtSPS proteins in Arabidopsis (At5g20280 (AtSPS1F), At5g11110 (AtSPS2F), At1g04920 (AtSPS3F) and At4g10120 (AtSPS4F)) were used as queries to search the tomato genome by BLASTp. Subsequently, the candidate sequences were subjected to the second round BLASTp against SWISS-PROT database to check whether the corresponding sequences were close to the identified SPS family in other plants. Batch-CD was then employed to screen whether the candidate proteins had the glucosyl transferase glycos-transf-1 domain and S6PP domain that are unique to sucrose phosphate synthase family, and partial and defective sequences were eliminated during manual verification. When there were multiple transcripts for a gene, the longest variant was selected. The sequence length, molecular weight, and isoelectric point information of the identified SPS proteins were obtained using the online tool—ExPASy (https://web.expasy.org/protparam, accessed on 10 March 2020). The online tool PLANT-PLOC (http://www.csbio.sjtu.edu.cn/bioinf/plant, accessed on 12 March 2020) was employed to predict the protein subcellular localizations [63], and Kinasephos (http://kinasephos.mbc.nctu.edu.tw, accessed on 15 March 2020) was used to predict the phosphorylation sites.
Identification of SUS and INV genes in tomato was as described for the identification of SPS genes. The query sequences for SUS protein identification were from Arabidopsis (At5g20830 (AtSUS1), At5g49410 (AtSUS2), At4g02280 (AtSUS3), At3g43190 (AtSUS4), At5g37180 (AtSUS5) and At1g73370 (AtSUS6)). The accessions of query proteins for INV identification were At3g52600 (AtCWINV2), At2g36190 (AtCWINV4), At3g13784 (AtCWINV5), At5g11920 (AtCWINV6), At1g55120 (AtFRUCT5), At3g13790 (AtBFRUCT1), At1g62660, At1g12240 (AtATV12), At1g56560 (AtAN-INVA), At4g34860 (AtAN-INVB), At3g06500 (AtAN-INVC), At1g35580 (AtAN-INVG), At4g09510 (AtAN-INV1), At1g22650 (AtAN-INVD), At1g72000 (AtAN-INVF), At5g22510 (AtAN-INVE), At3g05820 (AtAN-INVH). For SUS identification, the candidate sequences were screened with Batch-CDD online tool to see if they share the sucrose synthase and glucosyl-transferase domains that are unique to sucrose synthase [64]. For INV identification, the sequences that contain a glyco-hydro-32N, glyco-hydro-32C domain, a cysteine catalytic domain MWECP/V and a β-furosidase motif NDPNG/A were identified as acid invertase, whereas if the candidate sequences only contain a glyco-hydro-100 domain, they were identified as alkaline/neutral invertases [54].

4.2. Phylogenetic Relationship, Gene Structure and Protein Motif Analysis

A phylogenetic tree was constructed by MEGA-X using the maximum likelihood method with Poisson model and 1000 bootstrap replications [65]. SPS, SUS and INV protein sequences of Arabidopsis and rice were downloaded from UniProt website (https://sparql.uniprot.org, accessed on 2 March 2020). The gene structures of SlSPS, SlSUS and SlINV were analyzed by gene structure display server (GSDS) online program (http://gsds.cbi.pku.edu.cn, accessed on 7 April 2020). The NCBI BATCH CD-search tool (https://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi, accessed on 12 April 2020) was used to analyze the conserved domain structure based on the corresponding protein sequence, and Pfam database was chosen for searching. The conserved motifs were predicted by the MEME Suite tools (http://meme-suite.org, accessed on 12 April 2020) [66], and the number of motif parameters was limited to less than 10 manually.

4.3. Chromosome Distribution, Gene Duplication and Syntenic Relationship Analysis

The physical location information of SlSPS, SlSUS and SlINV genes was obtained from the GCF_000188115.4_SL3.0_genomic database in NCBI, and all the identified genes were mapped to the tomato chromosomes by TBtools V1.089 [67]. Multiple collinear scanning toolkits (MCScanX) with default parameters were used to analyze the gene duplication events [68]. The syntenic relationship of genes in tomato, Arabidopsis and rice was analyzed using Dual Synteny Plotter software [67].

4.4. Plant Preparation and Chemical Treatment

Tomato (‘Alisa Craig’) seeds were sterilized in warm water at 55 °C for 15 min, and then washed with sterile water at room temperature. The sterilized seeds were placed on two layers of filter paper in a petri dish and cultured in darkness at 30 °C for 2 d. The germinated seeds were sown in plugs and seedlings at two-leaf stage were transplanted into dark plastic boxes containing 1/4 Hoagland nutrient solution (pH 5.8) [69]. Seven days later, when the seedlings had three fully expanded leaves, the culture medium was replaced with 1/2 Hoagland solution (pH 5.8), in the absence or presence of 20% (w/w) polyethylene glycerol 6000 (PEG-6000), 150 mM NaCl, 1.5 mM H2O2, 50 μM abscisic acid (ABA) or 200 μM salicylic acid (SA). Seedlings cultured in 1/2 Hoagland solution were employed as control. After 12, 24 and 72 h of treatment, all the roots and leaves were separately collected, immediately frozen in liquid nitrogen and then stored at −80 °C for analysis of gene expression.

4.5. RNA Extraction and Quantitative RT-PCR Analysis

Total RNA of leaves or roots was extracted with an E.Z.N.A. Plant RNA Kit (Omega, GA, USA). Genomic DNA in the RNA samples was erased by both adsorption column in the E.Z.N.A plant RNA Kit (Omega) and DNase in HiScript® II Q RT SuperMix for qPCR (+gDNA wiper, Vazyme, Nanjing, China). The first strand of cDNA was synthesized with a HiScript® II Q RT SuperMix for qPCR. Quantitative RT-PCR was performed on a Real-Time PCR system (QuantStudioTM 5, Applied Biosystems, Waltham, MA, USA) using AceQ® qPCR SYBR Green Master Mix Kit (Vazyme). Actin, TIP41 and 17sRNA were employed as internal controls [42,70]. All the primers used were listed in Table S1. The relative expression levels of SlSUS, SlINV and SlSPS genes were calculated using the 2−∆∆Ct method [71].

4.6. Tissue Expression Characteristics of SPSs, SUSs and INVs in Tomato

The tissue expression patterns of SPS, SUS and INV genes in tomato were determined bioinformatically and experimentally. The expression values in different tissues were retrieved from large scales of RNA-Seq with Genevestigator (https://genevestigator.com, accessed on 11 May 2020), a software that provides the average signal intensity values of a gene from a high diversity of experiments covering different tissues, ages, and treatments [72,73]. The IDs of all SPS, SUS or INV genes were given as inputs in the Anatomy tool of Genevestigator and searched against the RNA-seq database of Solanum lycopersicum involving seed, roots, shoot apex, leaf, vascular tissue, parenchyma, floral meristem, flower, pollen, pistil, ovule, ovary wall, fruit, outer pericarp, pericarp, inner pericarp (columella), locular tissue, placenta and septum. The source of RNA-seq data are provided in Table S2. The expression SlSUSs, SlINVs and SlSPSs in leaves and roots was also analyzed by quantitative RT-PCR.

4.7. Statistical Analysis

The data were subjected to one-way analysis of variance with Duncan’s multiple range test at p < 0.05 using SPSS 18.0 (IBM, Armonk, NY, USA).

5. Conclusions

Four SPS genes, six SUS genes and nineteen INV genes were identified in tomato. Tandem and segmental duplications contributed to the expansion of SlSUS and SlINV families, and no gene duplication was found in the SlSPS family. SlSPS could be classified into three families (A, B and C), and SlSUS could be classified into three subfamilies (I, II and III). The SlINV genes included eleven acid invertase genes and eight A/N invertase genes. The differentiation of SlSPS and SlSUS might have completed the differentiation of monocotyledons and dicotyledons. The conserved motifs were mostly consistent within each protein family or subfamily. SlSPS, SlSUS and SlINV genes showed differential expressions among family members and tissues, and in response to different environmental stimuli. SlSPS1, SlSPS3, SlSUS1, SlSUS3, SlSUS4, SlINVAN5 and SlINVAN7 showed consistent expression responses to different treatments in the leaves and/or roots, and thus may be the major genes responding to exogenous stimuli. The results suggest that the individual isoforms of these gene families may demonstrate differential functions under environmental stimuli and in different tissues. The study not only provides the structure characteristics and evolution information of genes encoding sucrose metabolism related enzymes, but also lays a foundation for investigating the functions of individual isoforms of each enzyme in plant growth, development and stress tolerance.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijms22094698/s1.

Author Contributions

Conceptualization, H.G. and Y.D.; methodology, Y.D. and H.S.; software, Y.D. and L.Y.; validation, H.Z. and J.Z.; formal analysis, Y.D. and H.G.; investigation, Y.D.; writing—original draft preparation, Y.D. and H.S.; writing—review and editing, H.G. and Y.D.; visualization, Y.D. and H.G.; project administration, H.G.; funding acquisition, H.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Key Research and Development Program of China (2018YFD1000800) and National Natural Science Foundation of China (32072561, 31772290).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Stein, O.; Granot, D. An overview of sucrose synthases in plants. Front. Plant. Sci. 2019, 10, 95. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Pan, L.; Guo, Q.; Chai, S.; Cheng, Y.; Ruan, M.; Ye, Q.; Wang, R.; Yao, Z.; Zhou, G.; Li, Z.; et al. Evolutionary conservation and expression patterns of neutral/alkaline invertases in solanum. Biomolecules 2019, 9, 763. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Yang, J.; Zhang, J.; Li, C.; Zhang, Z.; Ma, F.; Li, M. Response of sugar metabolism in apple leaves subjected to short-term drought stress. Plant Physiol. Biochem. 2019, 141, 164–171. [Google Scholar] [CrossRef]
  4. Ruan, Y.-L. Sucrose metabolism: Gateway to diverse carbon use and sugar signaling. Ann. Rev. Plant. Biol. 2014, 65, 33–67. [Google Scholar] [CrossRef] [PubMed]
  5. Wan, H.; Wu, L.; Yang, Y.; Zhou, G.; Ruan, Y.-L. Evolution of sucrose metabolism: The dichotomy of invertases and beyond. Trends Plant. Sci. 2017, 23, 163–177. [Google Scholar] [CrossRef]
  6. Langenkämper, G.; Fung, R.W.M.; Newcomb, R.D.; Atkinson, R.G.; Gardner, R.C.; MacRae, E.A. Sucrose phosphate synthase genes in plants belong to three different families. J. Mol. Evol. 2002, 54, 322–332. [Google Scholar] [CrossRef]
  7. Okamura, M.; Aoki, N.; Hirose, T.; Yonekura, M.; Ohto, C.; Ohsugi, R. Tissue specificity and diurnal change in gene expression of the sucrose phosphate synthase gene family in rice. Plant Sci. 2011, 181, 159–166. [Google Scholar] [CrossRef]
  8. Castleden, C.K.; Aoki, N.; Gillespie, V.J.; Macrae, E.A.; Quick, W.P.; Buchner, P.H.; Foyer, C.H.; Furbank, R.T.; Lunn, J.E. Evolution and function of the sucrose-phosphate synthase gene families in wheat and other grasses. Plant Physiol. 2004, 135, 1753–1764. [Google Scholar] [CrossRef] [Green Version]
  9. Li, M.; Feng, F.; Cheng, L. Expression patterns of genes involved in sugar metabolism and accumulation during apple fruit development. PLoS ONE 2012, 7, 3. [Google Scholar]
  10. Ma, P.; Zhang, X.; Chen, L.; Zhao, Q.; Zhang, Q.; Hua, X.; Wang, Z.; Tang, H.; Yu, Q.; Zhang, M.; et al. Comparative analysis of sucrose phosphate synthase (SPS) gene family between Saccharum officinarum and Saccharum spontaneum. BMC Plant. Biol. 2020, 20, 1–15. [Google Scholar] [CrossRef]
  11. Xu, X.; Yang, Y.; Liu, C.; Sun, Y.; Zhang, T.; Hou, M.; Huang, S.; Yuan, H. The evolutionary history of the sucrose synthase gene family in higher plants. BMC Plant. Biol. 2019, 19, 1–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Wang, Z.; Wei, P.; Wu, M.; Xu, Y.; Li, F.; Luo, Z.; Zhang, J.; Chen, A.; Xie, X.; Cao, P.; et al. Analysis of the sucrose synthase gene family in tobacco: Structure, phylogeny, and expression patterns. Planta 2015, 242, 153–166. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Abdullah, M.; Cao, Y.; Cheng, X.; Meng, D.; Chen, Y.; Shakoor, A.; Gao, J.; Cai, Y. The sucrose synthase gene family in chinese pear (Pyrus bretschneideri rehd.): Structure, expression, and evolution. Molecules 2018, 23, 1144. [Google Scholar] [CrossRef] [Green Version]
  14. Ji, X.; Ende, W.V.; Laere, A.V.; Cheng, S.; Bennett, J. Structure, evolution, and expression of the two invertase gene families of rice. J. Mol. Evol. 2005, 60, 615–634. [Google Scholar] [CrossRef] [PubMed]
  15. Wang, L.; Zheng, Y.; Ding, S.; Zhang, Q.; Chen, Y.; Zhang, J. Molecular cloning, structure, phylogeny and expression analysis of the invertase gene family in sugarcane. BMC Plant. Biol. 2017, 17, 1. [Google Scholar] [CrossRef] [Green Version]
  16. Su, T.; Han, M.; Min, J.; Chen, P.; Mao, Y.; Huang, Q.; Tong, Q.; Liu, Q.; Fang, Y. Genome-wide survey of invertase encoding genes and functional characterization of an extracellular fungal pathogen-responsive invertase in Glycine max. Int. J. Mol. Sci. 2018, 19, 2395. [Google Scholar] [CrossRef] [Green Version]
  17. Qin, G.; Zhu, Z.; Wang, W.; Cai, J.; Chen, Y.; Li, L.; Tian, S. A tomato vacuolar invertase inhibitor mediates sucrose metabolism and influences fruit ripening. Plant Physiol. 2016, 172, 1596–1611. [Google Scholar] [CrossRef] [Green Version]
  18. Bahaji, A.; Baroja-Fernández, E.; Ricarte-Bermejo, A.; Sánchez-López, Á.M.; Muñoz, F.J.; Romero, J.M.; Ruiz, M.T.; Baslam, M.; Almagro, G.; Sesma, M.T.; et al. Characterization of multiple sps knockout mutants reveals redundant functions of the four arabidopsis sucrose phosphate synthase isoforms in plant viability, and strongly indicates that enhanced respiration and accelerated starch turnover can alleviate the blockage of sucrose biosynthesis. Plant. Sci. 2015, 238, 135–147. [Google Scholar]
  19. Anur, R.M.; Mufithah, N.; Sawitri, W.D.; Sakakibara, H.; Sugiharto, B. Overexpression of sucrose phosphate synthase enhanced sucrose content and biomass production in transgenic sugarcane. Plants (Basel) 2020, 9, 200. [Google Scholar] [CrossRef] [Green Version]
  20. Haigler, C.H.; Singh, B.; Zhang, D.; Hwang, S.; Wu, C.; Cai, W.X.; Hozain, M.; Kang, W.; Kiedaisch, B.; Strauss, R.E.; et al. Transgenic cotton over-producing spinach sucrose phosphate synthase showed enhanced leaf sucrose synthesis and improved fiber quality under controlled environmental conditions. Plant. Mol. Biol. 2007, 63, 815–832. [Google Scholar] [CrossRef] [PubMed]
  21. Park, J.-Y.; Canam, T.; Kang, K.-Y.; Ellis, D.D.; Mansfield, S.D. Over-expression of an Arabidopsis family A sucrose phosphate synthase (SPS) gene alters plant growth and fibre development. Transgenic Res. 2008, 17, 181–192. [Google Scholar] [CrossRef]
  22. Park, J.-Y.; Canam, T.; Kang, K.-Y.; Unda, F.; Mansfield, S.D. Sucrose phosphate synthase expression influences poplar phenology. Tree Physiol. 2009, 29, 937–946. [Google Scholar] [CrossRef] [Green Version]
  23. Baxter, C.J.; Foyer, C.H.; Turner, J.; Rolfe, S.A.; Quick, W.P. Elevated sucrose-phosphate synthase activity in transgenic tobacco sustains photosynthesis in older leaves and alters development. J. Exp. Bot. 2003, 54, 1813–1820. [Google Scholar] [CrossRef] [Green Version]
  24. Hirose, T.; Hashida, Y.; Aoki, N.; Okamura, M.; Yonekura, M.; Ohto, C.; Terao, T.; Ohsugi, R. Analysis of gene-disruption mutants of a sucrose phosphate synthase gene in rice, OsSPS1, shows the importance of sucrose synthesis in pollen germination. Plant. Sci. 2014, 225, 102–106. [Google Scholar] [CrossRef] [PubMed]
  25. Tian, H.; Ma, L.; Zhao, C.; Hao, H.; Gong, B.; Yu, X.; Wang, X. Antisense repression of sucrose phosphate synthase in transgenic muskmelon alters plant growth and fruit development. Biochem. Biophys. Res. Commun. 2010, 393, 365–370. [Google Scholar] [CrossRef] [PubMed]
  26. Barratt, D.H.P.; Derbyshire, P.; Findlay, K.; Pike, M.; Wellner, N.; Lunn, J.; Feil, R.; Simpson, C.; Maule, A.J.; Smith, A.M. Normal growth of Arabidopsis requires cytosolic invertase but not sucrose synthase. Proc. Natl. Acad. Sci. USA 2009, 106, 13124–13129. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Wei, Z.; Qu, Z.; Zhang, L.; Zhao, S.; Bi, Z.; Ji, X.; Wang, X.; Wei, H. Overexpression of poplar xylem sucrose synthase in tobacco leads to a thickened cell wall and increased height. PLoS ONE 2015, 10, e0120669. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Goren, S.; Lugassi, N.; Stein, O.; Yeselson, Y.; Schaffer, A.A.; David-Schwartz, R.; Granot, D. Suppression of sucrose synthase affects auxin signaling and leaf morphology in tomato. PLoS ONE 2017, 12, e0182334. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Gordon, A.J.; Minchin, F.R.; James, C.L.; Komina, O. Sucrose synthase in legume nodules is essential for nitrogen fixation. Plant Physiol. 1999, 120, 867–878. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Zhao, C.; Hua, L.-N.; Liu, X.-F.; Li, Y.-Z.; Shen, Y.-Y.; Guo, J.-X. Sucrose synthase FaSS1 plays an important role in the regulation of strawberry fruit ripening. Plant. Growth Regul. 2017, 81, 175–181. [Google Scholar] [CrossRef]
  31. Martín, M.L.; Lechner, L.; Zabaleta, E.J.; Salerno, G.L. A mitochondrial alkaline/neutral invertase isoform (A/N-InvC) functions in developmental energy-demanding processes in Arabidopsis. Planta 2013, 237, 813–822. [Google Scholar] [CrossRef] [PubMed]
  32. Jia, L.; Zhang, B.; Mao, C.; Li, J.; Wu, Y.; Wu, P.; Wu, Z. OsCYT-INV1 for alkaline/neutral invertase is involved in root cell development and reproductivity in rice (Oryza sativa L.). Planta 2008, 228, 51–59. [Google Scholar] [CrossRef] [PubMed]
  33. Zanor, M.I.; Osorio, S.; Nunes-Nesi, A.; Carrari, F.; Lohse, M.; Usadel, B.; Kühn, C.; Bleiss, W.; Giavalisco, P.; Willmitzer, L.; et al. RNA interference of lin5 in tomato confirms its role in controlling brix content, uncovers the influence of sugars on the levels of fruit hormones, and demonstrates the importance of sucrose cleavage for normal fruit development and fertility. Plant Physiol. 2009, 150, 1204–1218. [Google Scholar] [CrossRef] [Green Version]
  34. Takehara, K.; Murata, K.; Yamaguchi, T.; Yamaguchi, K.; Chaya, G.; Kido, S.; Iwasaki, Y.; Ogiwara, H.; Ebitani, T.; Miura, K. Thermo-responsive allele of sucrose synthase 3 (Sus3) provides high-temperature tolerance during the ripening stage in rice (Oryza sativa L.). Breed. Sci. 2018, 68, 336–342. [Google Scholar] [CrossRef] [Green Version]
  35. Wang, H.; Sui, X.; Guo, J.; Wang, Z.; Cheng, J.; Ma, S.; Li, X.; Zhang, Z. Antisense suppression of cucumber (Cucumis sativus L.) sucrose synthase 3 (CsSUS3) reduces hypoxic stress tolerance. Plant. Cell Environ. 2014, 37, 795–810. [Google Scholar] [CrossRef]
  36. Qian, W.; Xiao, B.; Wang, L.; Hao, X.; Yue, C.; Cao, H.; Wang, Y.; Li, N.; Yu, Y.; Zeng, J.; et al. CsINV5, a tea vacuolar invertase gene enhances cold tolerance in transgenic Arabidopsis. BMC Plant. Biol. 2018, 18, 228. [Google Scholar] [CrossRef]
  37. Liu, Y.-H.; Offler, C.E.; Ruan, Y.-L. Cell wall invertase promotes fruit set under heat stress by suppressing ros-independent cell death. Plant Physiol. 2016, 172, 163–180. [Google Scholar] [CrossRef] [Green Version]
  38. Xu, X.-x.; Hu, Q.; Yang, W.-n.; Jin, Y. The roles of cell wall invertase inhibitor in regulating chilling tolerance in tomato. BMC Plant. Biol. 2017, 17, 195. [Google Scholar] [CrossRef]
  39. Wang, F.; Smith, A.G.; Brenner, M.L. Isolation and sequencing of tomato fruit sucrose synthase cDNA. Plant Physiol. 1993, 103, 1463–1464. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. D’Aoust, M.-A.; Yelle, S.; Nguyen-Quoc, B. Antisense inhibition of tomato fruit sucrose synthase decreases fruit setting and the sucrose unloading capacity of young fruit. Plant. Cell 1999, 11, 2407–2418. [Google Scholar] [CrossRef] [Green Version]
  41. Chengappa, S.; Loader, N.; Shields, R. Cloning, expression, and mapping of a second tomato sucrose synthase gene, Sus3. Plant Physiol. 1998, 118, 1533. [Google Scholar]
  42. Goren, S.; Huber, S.C.; Granot, D. Comparison of a novel tomato sucrose synthase, SlSUS4, with previously described SlSUS isoforms reveals distinct sequence features and differential expression patterns in association with stem maturation. Planta 2011, 233, 1011–1023. [Google Scholar] [CrossRef]
  43. Godt, D.E.; Roitsch, T. Regulation and tissue-specific distribution of mRNAs for three extracellular invertase isoenzymes of tomato suggests an important function in establishing and maintaining sink metabolism. Plant Physiol. 1997, 115, 273–282. [Google Scholar] [CrossRef] [Green Version]
  44. Klann, E.; Yelle, S.; Bennett, A.B. Tomato fruit acid invertase complementary dna nucleotide and deduced amino acid sequences. Plant Physiol. 1992, 99, 351–353. [Google Scholar] [CrossRef]
  45. Leskow, C.C.; Conte, M.; del Pozo, T.; Bermúdez, L.; Lira, B.S.; Gramegna, G.; Baroli, I.; Burgos, E.; Zavallo, D.; Kamenetzky, L.; et al. The cytosolic invertase NI6 affects vegetative growth, flowering, fruit set and yield in tomato. J. Exp. Bot. 2020, 72, 2525–2543. [Google Scholar] [CrossRef] [PubMed]
  46. Lunn, J.E.; Ashton, A.R.; Hatch, M.D.; Heldt, H.W. Purification, molecular cloning, and sequence analysis of sucrose-6F-phosphate phosphohydrolase from plants. Proc. Natl. Acad. Sci. USA 2000, 97, 12914–12919. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Lunn, J.E.; MacRae, E. New complexities in the synthesis of sucrose. Curr. Opin. Plant. Biol. 2003, 6, 208–214. [Google Scholar] [CrossRef]
  48. Lammens, W.; Roy, K.L.; Schroeven, L.; Laere, A.V.; Rabijns, A.; Ende, W.V.d.E. Structural insights into glycoside hydrolase family 32 and 68 enzymes: Functional implications. J. Exp. Bot. 2009, 60, 727–740. [Google Scholar] [CrossRef] [Green Version]
  49. Cannon, S.B.; Mitra, A.; Baumgarten, A.; Young, N.D.; May, G. The roles of segmental and tandem gene duplication in the evolution of large gene families in Arabidopsis thaliana. BMC Plant. Biol. 2004, 4, 10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Wang, D.; Zhao, J.; Hu, B.; Li, J.; Qin, Y.; Chen, L.; Qin, Y.; Hu, G. Identification and expression profile analysis of the sucrose phosphate synthase gene family in Litchi chinensis Sonn. PeerJ 2018, 6, e4379. [Google Scholar] [CrossRef] [Green Version]
  51. Mitchell-Olds, T.; Clauss, M.J. Plant evolutionary genomics. Curr. Opin.Plant. Biol. 2002, 5, 74–79. [Google Scholar] [CrossRef]
  52. Jiang, S.-Y.; Chi, Y.-H.; Wang, J.-Z.; Zhou, J.-X.; Cheng, Y.-S.; Zhang, B.-L.; Ma, A.; Vanitha, J.; Ramachandran, S. Sucrose metabolism gene families and their biological functions. Sci. Rep. 2015, 5, 17583. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Bocock, P.N.; Morse, A.M.; Dervinis, C.; Davis, J.M. Evolution and diversity of invertase genes in Populus trichocarpa. Planta 2008, 227, 565–576. [Google Scholar] [CrossRef] [PubMed]
  54. Sturm, A. Invertases. Primary structures, functions, and roles in plant development and sucrose partitioning. Plant Physiol. 1999, 121, 1–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Shen, L.-B.; Yao, Y.; He, H.; Qin, Y.-L.; Liu, Z.-J.; Liu, W.-X.; Qi, Z.-Q.; Yang, L.-J.; Cao, Z.-M.; Yang, Y. Genome-wide identification, expression, and functional analysis of the alkaline/neutral invertase gene family in pepper. Int. J. Mol. Sci. 2018, 19, 224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Yao, S.G.; Kodama, R.; Wang, H.; Ichii, M.; Taketa, S.; Yoshida, H. Analysis of the rice SHORT-ROOT5 gene revealed functional diversification of plant neutral/alkaline invertase family. Plant. Sci. 2009, 176, 627–634. [Google Scholar] [CrossRef]
  57. Palmer, W.M.; Ru, L.; Jin, Y.; Patrick, J.W.; Ruan, Y.-L. Tomato ovary-to-fruit transition is characterized by a spatial shift of mRNAs for cell wall invertase and its inhibitor with the encoded proteins localized to sieve elements. Mol. Plant. 2015, 8, 315–328. [Google Scholar] [CrossRef] [Green Version]
  58. Proels, R.K.; González, M.-C.; Roitsch, T. Gibberellin-dependent induction of tomato extracellular invertase Lin7 is required for pollen development. Funct. Plant. Biol. 2006, 33, 547–554. [Google Scholar] [CrossRef]
  59. Nguyen-Quoc, B.; Foyer, C.H. A role for ‘futile cycles’ involving invertase and sucrose synthase in sucrose metabolism of tomato fruit. J. Exp. Bot. 2001, 52, 881–889. [Google Scholar] [CrossRef] [Green Version]
  60. Hu, T.; Hu, L.; Zhang, X.; Zhang, P.; Zhao, Z.; Fu, J. Differential responses of CO2 assimilation, carbohydrate allocation and gene expression to nacl stress in perennial ryegrass with different salt tolerance. PLoS ONE 2013, 8, e66090. [Google Scholar]
  61. Solís-Guzmán, M.G.; Argüello-Astorga, G.; López-Bucio, J.; Ruiz-Herrera, L.F.; López-Meza, J.E.; Sánchez-Calderón, L.; Carreón-Abud, Y.; Martínez-Trujillo, M. Arabidopsis thaliana sucrose phosphate synthase (SPS) genes are expressed differentially in organs and tissues, and their transcription is regulated by osmotic stress. Gene Expr. Patterns 2017, 25, 92–101. [Google Scholar] [CrossRef]
  62. Dahro, B.; Wang, F.; Peng, T.; Liu, J.-H. PtrA/NINV, an alkaline/neutral invertase gene of Poncirus trifoliata, confers enhanced tolerance to multiple abiotic stresses by modulating ROS levels and maintaining photosynthetic efficiency. BMC Plant. Biol. 2016, 16, 76. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Chou, K.-C.; Shen, H.-B. Large-scale plant protein subcellular location prediction. J. Cell. Biochem. 2007, 100, 665–678. [Google Scholar] [CrossRef] [PubMed]
  64. Punta, M.; Coggill, P.C.; Eberhardt, R.Y.; Mistry, J.; Tate, J.G.; Boursnell, C.; Pang, N.; Forslund, K.; Ceric, G.; Clements, J.; et al. The Pfam protein families database. Nucleic Acids Res. 2000, 30, 276–280. [Google Scholar] [CrossRef] [PubMed]
  65. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef] [PubMed]
  66. Bailey, T.L.; Bodén, M.; Buske, F.A.; Frith, M.C.; Grant, C.E.; Clementi, L.; Ren, J.; Li, W.W.; Noble, W.S. MEME Suite: Tools for motif discovery and searching. Nucleic Acids Res. 2009, 37, 202–208. [Google Scholar] [CrossRef] [PubMed]
  67. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An integrative toolkit developed for interactive analyses of big biological data. Mol. Plant. 2020, 13, 1194–1202. [Google Scholar] [CrossRef]
  68. Wang, Y.; Tang, H.; DeBarry, J.D.; Tan, X.-f.; Li, J.; Wang, X.; Lee, T.-H.; Jin, H.; Marler, B.; Guo, H.; et al. MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012, 40, e49. [Google Scholar] [CrossRef] [Green Version]
  69. Hoagland, D.R.; Arnon, D.I. The water culture method for growing plants without soil. Calif. Agric. Exp. Stn. Circ. 1950, 347, 32. [Google Scholar]
  70. Expósito-Rodríguez, M.; Borges, A.A.; Borges-Pérez, A.; Pérez, J.A. Selection of internal control genes for quantitative real-time RT-PCR studies during tomato development process. BMC Plant. Biol. 2008, 8, 131. [Google Scholar] [CrossRef] [Green Version]
  71. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative pcr and the 2(-delta delta C(T)) method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  72. Hruz, T.; Laule, O.; Szabo, G.; Wessendorp, F.; Bleuler, S.; Oertle, L.; Widmayer, P.; Gruissem, W.; Zimmermann, P. Genevestigator v3: A reference expression database for the meta-analysis of transcriptomes. Adv. Bioinform. 2008, 2008, 420747. [Google Scholar] [CrossRef] [PubMed]
  73. Zimmermann, P.; Hirsch-Hoffmann, M.; Hennig, L.; Gruissem, W. GENEVESTIGATOR. Arabidopsis microarray database and analysis toolbox. Plant Physiol. 2004, 136, 2621–2632. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Ijms 22 04698 g001 550
Figure 1. Phylogenetic analysis of SPS (A), SUS (B), acid INV proteins (C) and alkaline/neutral INV proteins (D). Trees were built using amino acid sequences from Arabidopsis thaliana, tomato (Solanum lycopersicum) and rice (Oryza sativa). The phylogenetic tree was constructed using the maximum likelihood method with MEGA-X.
Figure 1. Phylogenetic analysis of SPS (A), SUS (B), acid INV proteins (C) and alkaline/neutral INV proteins (D). Trees were built using amino acid sequences from Arabidopsis thaliana, tomato (Solanum lycopersicum) and rice (Oryza sativa). The phylogenetic tree was constructed using the maximum likelihood method with MEGA-X.
Ijms 22 04698 g001
Ijms 22 04698 g002 550
Figure 2. Gene structure and conserved motif analysis of SPS (A), SUS (B), acid INV (C) and alkaline/neutral INV (D) in tomato. In the gene structure, the filled boxes and lines represent exons and introns, respectively. The green and yellow boxes, respectively, represent UTR sequence and CDS. In the motif map, different conserved motifs were shown in different colors.
Figure 2. Gene structure and conserved motif analysis of SPS (A), SUS (B), acid INV (C) and alkaline/neutral INV (D) in tomato. In the gene structure, the filled boxes and lines represent exons and introns, respectively. The green and yellow boxes, respectively, represent UTR sequence and CDS. In the motif map, different conserved motifs were shown in different colors.
Ijms 22 04698 g002
Ijms 22 04698 g003 550
Figure 3. Chromosome localization of tomato SPS genes (A) and synteny analysis of SPS genes from tomato, Arabidopsis and rice (B). Chr, chromosome.
Figure 3. Chromosome localization of tomato SPS genes (A) and synteny analysis of SPS genes from tomato, Arabidopsis and rice (B). Chr, chromosome.
Ijms 22 04698 g003
Ijms 22 04698 g004 550
Figure 4. Chromosome localization of tomato SUS genes (A) and synteny analysis of SUS genes from tomato (B), Arabidopsis and rice. The red arc line indicates tandem duplicated gene pair, and the red straight line indicates segmental duplicated gene pair. Chr, chromosome.
Figure 4. Chromosome localization of tomato SUS genes (A) and synteny analysis of SUS genes from tomato (B), Arabidopsis and rice. The red arc line indicates tandem duplicated gene pair, and the red straight line indicates segmental duplicated gene pair. Chr, chromosome.
Ijms 22 04698 g004
Ijms 22 04698 g005 550
Figure 5. Chromosome localization of tomato INV genes (A) and synteny analysis of INV genes from tomato, Arabidopsis and rice (B). The red arc lines indicate tandem duplicated gene pairs, and the red straight line indicates segmental duplicated gene pair. Chr, chromosome.
Figure 5. Chromosome localization of tomato INV genes (A) and synteny analysis of INV genes from tomato, Arabidopsis and rice (B). The red arc lines indicate tandem duplicated gene pairs, and the red straight line indicates segmental duplicated gene pair. Chr, chromosome.
Ijms 22 04698 g005
Ijms 22 04698 g006 550
Figure 6. Tissue expression pattern of SPS, SUS and INV genes in tomato. The expression values were retrieved from large scale RNA-seq with Genevestigator. (A) SPS expression; (B) SUS expression; (C) INV expression.
Figure 6. Tissue expression pattern of SPS, SUS and INV genes in tomato. The expression values were retrieved from large scale RNA-seq with Genevestigator. (A) SPS expression; (B) SUS expression; (C) INV expression.
Ijms 22 04698 g006
Ijms 22 04698 g007 550
Figure 7. Expression of SPS, SUS and INV genes in leaves and roots of tomato. Expression was monitored using quantitative RT-PCR at three-leaf stage. (A) SPS expression; (B) SUS expression; (C) INV expression.
Figure 7. Expression of SPS, SUS and INV genes in leaves and roots of tomato. Expression was monitored using quantitative RT-PCR at three-leaf stage. (A) SPS expression; (B) SUS expression; (C) INV expression.
Ijms 22 04698 g007
Ijms 22 04698 g008 550
Figure 8. Changes in the expression of SPS family genes under different exogenous stimuli conditions. The three-leaf age seedlings were treated with 20% (w/w) polyethylene glycerol-6000 (PEG-6000), 150 mM NaCl, 1.5 mM H2O2, 50 μM abscisic acid (ABA) or 200 μM salicylic acid (SA) for 12 to 72 h, after which the leaves and roots were collected for gene expression analysis using quantitative RT-PCR. (A) SPS1 expression in leaves; (B) SPS1 expression in roots; (C) SPS2 expression in leaves; (D) SPS2 expression in roots; (E) SPS3 expression in leaves; (F) SPS3 expression in roots; (G) SPS4 expression in leaves; (H) SPS4 expression in roots. Actin, TIP41 and 17sRNA were used as internal controls. Data are means ± SD (n = 3). Different letters above bars indicate a significant difference at p < 0.05. CT, control.
Figure 8. Changes in the expression of SPS family genes under different exogenous stimuli conditions. The three-leaf age seedlings were treated with 20% (w/w) polyethylene glycerol-6000 (PEG-6000), 150 mM NaCl, 1.5 mM H2O2, 50 μM abscisic acid (ABA) or 200 μM salicylic acid (SA) for 12 to 72 h, after which the leaves and roots were collected for gene expression analysis using quantitative RT-PCR. (A) SPS1 expression in leaves; (B) SPS1 expression in roots; (C) SPS2 expression in leaves; (D) SPS2 expression in roots; (E) SPS3 expression in leaves; (F) SPS3 expression in roots; (G) SPS4 expression in leaves; (H) SPS4 expression in roots. Actin, TIP41 and 17sRNA were used as internal controls. Data are means ± SD (n = 3). Different letters above bars indicate a significant difference at p < 0.05. CT, control.
Ijms 22 04698 g008
Ijms 22 04698 g009 550
Figure 9. Changes in the expression of SUS family genes under different exogenous stimuli conditions. The three-leaf age seedlings were treated with 20% (w/w) polyethylene glycerol-6000 (PEG-6000), 150 mM NaCl, 1.5 mM H2O2, 50 μM abscisic acid (ABA) or 200 μM salicylic acid (SA) for 12 to 72 h, after which the leaves and roots were collected for gene expression analysis using quantitative RT-PCR. (A) SUS1 expression in leaves; (B) SUS1 expression in roots; (C) SUS3 expression in leaves; (D) SUS3 expression in roots; (E) SUS4 expression in leaves; (F) SUS4 expression in roots; (G) SUS5 expression in leaves; (H) SUS5 expression in roots; (I) SUS6 expression in leaves; (J) SUS6 expression in roots; (K) SUS7 expression in leaves; (L) SUS7 expression in roots. Actin, TIP41 and 17sRNA were used as internal controls. Data are means ± SD (n = 3). Different letters above bars indicate a significant difference at p < 0.05. CT, control.
Figure 9. Changes in the expression of SUS family genes under different exogenous stimuli conditions. The three-leaf age seedlings were treated with 20% (w/w) polyethylene glycerol-6000 (PEG-6000), 150 mM NaCl, 1.5 mM H2O2, 50 μM abscisic acid (ABA) or 200 μM salicylic acid (SA) for 12 to 72 h, after which the leaves and roots were collected for gene expression analysis using quantitative RT-PCR. (A) SUS1 expression in leaves; (B) SUS1 expression in roots; (C) SUS3 expression in leaves; (D) SUS3 expression in roots; (E) SUS4 expression in leaves; (F) SUS4 expression in roots; (G) SUS5 expression in leaves; (H) SUS5 expression in roots; (I) SUS6 expression in leaves; (J) SUS6 expression in roots; (K) SUS7 expression in leaves; (L) SUS7 expression in roots. Actin, TIP41 and 17sRNA were used as internal controls. Data are means ± SD (n = 3). Different letters above bars indicate a significant difference at p < 0.05. CT, control.
Ijms 22 04698 g009
Ijms 22 04698 g010 550
Figure 10. Changes in the expression of INV family genes under different exogenous stimuli conditions. The plants were treated with 20%(w/w) polyethylene glycerol-6000 (PEG-6000), 150 mM NaCl, 1.5 mM H2O2, 50 μM abscisic acid (ABA) or 200 μM salicylic acid (SA) for 12 to 72 h, after which the leaves and roots were collected for gene expression analysis using quantitative RT-PCR. (A) INVAN1 expression in leaves; (B) INVAN1 expression in roots; (C) INVAN2 expression in leaves; (D) INVAN2 expression in roots; (E) INVAN3 expression in leaves; (F) INVAN3 expression in roots; (G) INVAN4 expression in leaves; (H) INVAN4 expression in roots; (I) INVAN5 expression in leaves; (J) INVAN5 expression in roots; (K) INVAN6 expression in leaves; (L) INVAN6 expression in roots; (M) INVAN7 expression in leaves; (N) INVAN7 expression in roots; (O) INVAN8 expression in leaves; (P) INVAN8 expression in roots; (Q) INVVR1 expression in leaves; (R) INVVR1 expression in roots; (S) INVVR2 expression in leaves; (T) INVVR2 expression in roots; (U) INVCW5 expression in leaves; (V) INVCW5 expression in roots; (W) INVCW6 expression in leaves; (X) INVCW6 expression in roots. Actin, TIP41 and 17sRNA were used as internal controls. Data are means ± SD (n = 3). Different letters above bars indicate a significant difference at p < 0.05. CT, control. UD, undetected. *** p < 0.001.
Figure 10. Changes in the expression of INV family genes under different exogenous stimuli conditions. The plants were treated with 20%(w/w) polyethylene glycerol-6000 (PEG-6000), 150 mM NaCl, 1.5 mM H2O2, 50 μM abscisic acid (ABA) or 200 μM salicylic acid (SA) for 12 to 72 h, after which the leaves and roots were collected for gene expression analysis using quantitative RT-PCR. (A) INVAN1 expression in leaves; (B) INVAN1 expression in roots; (C) INVAN2 expression in leaves; (D) INVAN2 expression in roots; (E) INVAN3 expression in leaves; (F) INVAN3 expression in roots; (G) INVAN4 expression in leaves; (H) INVAN4 expression in roots; (I) INVAN5 expression in leaves; (J) INVAN5 expression in roots; (K) INVAN6 expression in leaves; (L) INVAN6 expression in roots; (M) INVAN7 expression in leaves; (N) INVAN7 expression in roots; (O) INVAN8 expression in leaves; (P) INVAN8 expression in roots; (Q) INVVR1 expression in leaves; (R) INVVR1 expression in roots; (S) INVVR2 expression in leaves; (T) INVVR2 expression in roots; (U) INVCW5 expression in leaves; (V) INVCW5 expression in roots; (W) INVCW6 expression in leaves; (X) INVCW6 expression in roots. Actin, TIP41 and 17sRNA were used as internal controls. Data are means ± SD (n = 3). Different letters above bars indicate a significant difference at p < 0.05. CT, control. UD, undetected. *** p < 0.001.
Ijms 22 04698 g010
Table
Table 1. Gene structure and protein characteristics of SUSs, INVs and SPSs in tomato.
Table 1. Gene structure and protein characteristics of SUSs, INVs and SPSs in tomato.
GeneGene StructureProtein CharacteristicsGene Name in the Literature
NameLocusChrSequence Length (bp)CDS Length (bp)Exon No.Intron No.Length (Aa)IPMW (KDa)Subcellular LocalizationPhosphorylation Site No.(References)
SPS1Solyc07g007790785113165131210546.05118.5PM24SPSA1 [17]
SPS2Solyc08g042000872833147131210486.26117.9PM21SPSA2 [17]
SPS3Solyc09g092130982703195121110646.13119.6PM15SPSB [17]
SPS4Solyc11g0451101175343090141310296.59116.2PM21SPSC [17]
SUS1Solyc12g009300125834241813128055.9492.5PM9SUS1 [39]; TOMMSSF [40]
SUS3Solyc07g04255075663241813128055.9692.5PM13SUS3 [41]
SUS4Solyc09g09859096464243915148125.9192.9PM18SUS4 [42]
SUS5Solyc07g04252073624241211108035.9791.6PM13SS5 [17]
SUS6Solyc03g09829034379267615148915.87100.7PM20SS6 [17]
SUS7Solyc02g08130024757255014138496.9496.3PM14SS7 [17]
INVCW1Solyc03g121680327161716655716.864.2CW13
INVCW2Solyc06g064620628341596875314.9659.7CW8
INVCW3Solyc09g010080938631755655849.267.2CW *8LIN5 [43]
INVCW4Solyc09g010090935301752655836.9366.2CW14LIN7 [43]
INVCW5Solyc10g0832901055131749655829.2365.9CW9LIN6 [43]
INVCW6Solyc10g0833001045281560655198.6958.8CW8LIN8 [43]
INVCW7Solyc10g0853601044591737655788.765.9CW5
INVCW8Solyc10g0856401034981731655767.2665.7CW5
INVCW9Solyc10g0856501036591737655786.3966.6CW10
INVVR1Solyc03g083910341871911766365.5470.1Vac8VI [44]
INVVR2Solyc08g079080841941959766526.2172.8Vac9LIN9 [17]
INVAN1Solyc01g058010156942010766695.8476.0Mito10NI1 [17]
INVAN2Solyc01g100810158521962656538.1874.5Mito13NI2 [4,17]
INVAN3Solyc01g111100151891707435686.4564.9Cyt13CIN2 [4]
INVAN4Solyc04g081440447131713435705.9765.2Cytosol *13CIN3 [4]; NI6 [45]
INVAN5Solyc06g065210657881656435516.1662.8Cyt17NI3 [17]; CIN4 [4]
INVAN6Solyc11g0072701160871968656555.8473.5Chl12CIN5 [4]
INVAN7Solyc11g0206101161751659435526.0662.6Cyt16NI4 [17]; CIN6 [4]
INVAN8Solyc11g0670501142912019656726.8976.2Chl8NI5 [17]; CIN7 [4]
* Verified experimentally. Abbreviations: Aa, amino acid; CDS, coding DNA sequence; Chl, chloroplast; Chr, chromosome; CW, cell wall; Cyt, cytoplasmic; IP, isoelectric point; Mito, mitochondrial; MW, molecular weight; PM, plasma membrane; Vac, vacuole.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Duan, Y.; Yang, L.; Zhu, H.; Zhou, J.; Sun, H.; Gong, H. Structure and Expression Analysis of Sucrose Phosphate Synthase, Sucrose Synthase and Invertase Gene Families in Solanum lycopersicum. Int. J. Mol. Sci. 2021, 22, 4698. https://doi.org/10.3390/ijms22094698

AMA Style

Duan Y, Yang L, Zhu H, Zhou J, Sun H, Gong H. Structure and Expression Analysis of Sucrose Phosphate Synthase, Sucrose Synthase and Invertase Gene Families in Solanum lycopersicum. International Journal of Molecular Sciences. 2021; 22(9):4698. https://doi.org/10.3390/ijms22094698

Chicago/Turabian Style

Duan, Yaoke, Lan Yang, Haijia Zhu, Jie Zhou, Hao Sun, and Haijun Gong. 2021. "Structure and Expression Analysis of Sucrose Phosphate Synthase, Sucrose Synthase and Invertase Gene Families in Solanum lycopersicum" International Journal of Molecular Sciences 22, no. 9: 4698. https://doi.org/10.3390/ijms22094698

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

Article Metrics

Back to TopTop