Characterization of FBA genes in potato (Solanum tuberosum L.) and expression patterns in response to light spectrum and abiotic stress

Fructose-1, 6-bisphosphate aldolase (FBA) plays vital roles in plant growth, development, and response to abiotic stress. However, genome-wide identification and structural characterization of the potato (Solanum tuberosum L.) FBA gene family has not been systematically analyzed. In this study, we identified nine StFBA gene members in potato, with six StFBA genes localized in the chloroplast and three in the cytoplasm. The analysis of gene structures, protein structures, and phylogenetic relationships indicated that StFBA genes were divided into Class I and II, which exhibited significant differences in structure and function. Synteny analysis revealed that segmental duplication events promoted the expansion of the StFBA gene family. Promoter analysis showed that most StFBA genes contained cis-regulatory elements associated with light and stress responses. Expression analysis showed that StFBA3, StFBA8, and StFBA9 showing significantly higher expression levels in leaf, stolon, and tuber under blue light, indicating that these genes may improve photosynthesis and play an important function in regulating the induction and expansion of microtubers. Expression levels of the StFBA genes were influenced by drought and salt stress, indicating that they played important roles in abiotic stress. This work offers a theoretical foundation for in-depth understanding of the evolution and function of StFBA genes, as well as providing the basis for the genetic improvement of potatoes.

1 Introduction

Plants utilize various enzymes within their bodies to absorb solar energy and convert CO₂ into organic carbon in nature (Wildman, 2002). Ribulose-1, 5-bisphosphate carboxylase/oxygenase (Rubisco) plays a crucial role in carbon fixation and regulates biomass synthesis during photosynthesis (Ashida et al., 2019). However, the low catalytic efficiency of Rubisco remains a major limiting factor in photosynthesis (Bracher et al., 2017). With the rise in global temperatures and CO₂ concentrations, Rubisco is susceptible to oxidation, resulting in substantial loss of photosynthetic products. Some plants conserve energy by reducing the synthesis of Rubisco proteins when grown in high CO₂ environments (Jordan and Ogren, 1984; Long et al., 2004). Previous studies have shown that increasing the levels of enzymes associated with the regeneration phase of ribulose-1, 5-bisphosphate (RuBP) in the Calvin cycle, such as sedoheptulose-1, 7-bisphosphatase (SBP), transketolase, and aldolase, can enhance plant carbon assimilation efficiency (Raines, 2003; Lefebvre et al., 2005). Among them, aldolase reaches its highest activity level under CO₂ saturation (Haake et al., 1999). Overexpression of a gene-encoding plastid aldolase in transgenic tobacco enhanced the activity of the enzyme, accelerated the regeneration of RuBP, and the biomass significantly increased with increasing CO₂ levels (Uematsu et al., 2012).

Fructose-1, 6-bisphosphate aldolase (FBA; EC 4.1.2.13) is an important enzyme involved in plant photosynthesis and energy metabolism (Carrera et al., 2021). It catalyzes the condensation of glyceraldehyde-3-phosphate (GAP) and dihydroxyacetone phosphate (DHAP) to fructose-1, 6-bisphosphate (FBP) during glycolysis and its reversible reaction, which is ultimately used for the regeneration of RuBP (Zgiby et al., 2000; Capodagli et al., 2014). FBA can be divided into Class I and Class II. Class I aldolase is primarily found in eukaryotes, while Class II aldolase is widely distributed in bacteria, saccharomycetes, and fungi. There are significant differences in the structure and function of these classes (Marsh and Lebherz, 1992; Izard and Sygusch, 2004; Galkin et al., 2007). In higher plants, Class I aldolase is composed of homotetramers (Lebherz et al., 1984; Perham, 1990) and its protein tertiary structure contains a typical conserved structural domain, called the TIM barrel (Hester et al., 1991; Blom and Sygusch, 1997). The TIM barrel is a conserved protein fold pattern, which consists of eight repetitive α-helices and eight β-folds, which are arranged alternately, forming a characteristic α/β-barrel structure present in the protein’s main peptide segment (Lorentzen et al., 2003; Frenkel and Trifonov, 2005). In addition, there are two subtypes of Class I aldolase: chloroplast aldolase, and cytoplasmic aldolase (Pelzer-Reith et al., 1993; Tsutsumi et al., 1994). Chloroplast aldolase is a bifunctional enzyme in the Calvin cycle, involved in the synthesis of FBP and sedoheptulose-1, 7-bisphosphate (SuBP). It catalyzes the generation of metabolic products in the process of starch synthesis (Sonnewald et al., 1994; Flechner et al., 1999). Furthermore, cytoplasmic aldolase is an important metabolic enzyme in the pathways of glycolysis and gluconeogenesis, where it catalyzes the production of FBP and participates in the synthesis of sucrose (Gross et al., 1999; Fan et al., 2009). Mutations in Arabidopsis cytoplasmic aldolase AtFBA3 have been reported to affect the plastid glycolytic metabolism in the roots of plants, thereby limiting the synthesis of essential compounds such as amino acids (Carrera et al., 2021).

Previous studies reported that a slight reduction in FBA activity significantly inhibited the photosynthetic and carbon assimilations of antisense transgenic plants, which affected the synthesis of starch and severely inhibited the growth of plants (Jacquot et al., 1995; Haake et al., 1998; Bryant, 2000; Henkes et al., 2001), while reduced activity of FBA in tomato plants led to smaller seed size and slower plant growth (Cai et al., 2018). Suppressing the expression of the chloroplast FBA gene affects the development of tomatoes (Obiadalla-Ali et al., 2004). In rice, the OsALD-Y (fructose-1, 6-bisphosphate aldolase) encodes an FBA protein, which affects photosynthesis and sugar metabolism (Zhang et al., 2016). Furthermore, the FBA gene is also involved in a variety of physiological and biochemical processes, such as the regulation of secondary metabolism (Zeng et al., 2014), and regulation of plant growth and development (Haake et al., 1998; Obiadalla-Ali et al., 2004; Mutuku and Nose, 2012), responding to various biotic and abiotic stress (Lu et al., 2012; Cai et al., 2018), to regulate plant responses to hormones such as abscisic acid (Osakabe et al., 2005) and gibberellin (Konishi et al., 2004), as well as responses to environmental light signals (Oelze et al., 2014). It is therefore possible that FBA genes play a crucial role in enhancing the photosynthetic efficiency of plants, to improve crop yield and quality.

Potato is one of the four largest food crops worldwide (Zhang et al., 2017), and a significant crop for both food and as a vegetable in our country. The available genome sequences of potatoes (Pham et al., 2020; Zhou et al., 2020; Zhang et al., 2021) facilitate their use in potato growth and development. Members of the FBA gene family have been identified in different species, with eight identified in Arabidopsis thaliana (Lu et al., 2012), eight in Solanum lycopersicum (Cai et al., 2016), 22 in Brassica napus (Zhao et al., 2019), 16 in Nicotiana tabacum (Zhao et al., 2021), and 17 in cotton (Li et al., 2021). When expression of potato chloroplast aldolase (Registry Number: Y10380) was inhibited using antisense repression of expression technology, the photosynthesis of plants was decreased to significantly affect the sugar and starch contents of potatoes (Haake et al., 1999). However, there is still a lack of systematic genome-wide identification and functional investigation of potato FBA gene family members. In this study, we identified nine StFBA genes using bioinformatics tools, to systematically analyze gene structures, protein structures, chromosome distributions, phylogenetic relationships, synteny, collinearities, and cis-acting elements in promoter regions. Furthermore, we determined the expression patterns of StFBA in different tissues and in response to abiotic stress. The results provide a theoretical basis for the in-depth use of StFBA functions to improve photosynthetic efficiency and tuber development in potatoes.

2 Materials and methods

2.1 Plant materials

Potato (cv. Shepody) was used as experimental material, provided by the Crop Biotechnology Research Center of Anhui University of Science and Technology. Under sterile conditions, potato plantlets were cut into 1–1.5 cm segments with leaves and then placed into medium in tissue culture bottles (inner diameter, 63 mm; height, 85 mm). Each bottle contained four stem segments and was placed in a tissue culture room with relative humidity of 65% ± 5%, daytime temperature of 22°C ± 2°C, night temperature of 18°C ± 2°C, photoperiod uniformly set to 8 h light/16 h dark, and photosynthetic photon flux density of 65 μmol m⁻² s⁻¹. The medium used for propagation of plantlets in vitro was solid Murashige and Skoog (MS) medium (4%, w/v, sucrose, 0.9%, w/v, agar). The medium used for induction of microtubers was solid, higher sucrose MS medium (8%, w/v, sucrose, 0.9%, w/v, agar). Potato plantlets grown for 30 days were acclimated and transplanted to pots containing vegetative soil and vermiculite (1:1) and placed in a plant growth room where the relative humidity was 65% ± 5%, the temperature was 22°C ± 2°C, and the photoperiod was uniformly set at 12 h light/12 h dark. Light intensity was 200 μmol m⁻² s⁻¹. The flower tissue, petal, stamen, and pistil at 45 days of growth, as well as the leaf, petiole, stem, stolon, root, and tuber at 60 days of growth were collected for tissue expression analysis of the StFBA gene.

2.2 Identification of StFBA members in potato

The whole genome sequence, protein and encode sequences of potato v6.1 were downloaded from the Phytozome v13 website (https://phytozome.jgi.doe.gov/pz/portal.html). The protein sequences of the Arabidopsis FBA gene family were sourced from the TAIR database (http://www.arabidopsis.org). Genomic data for rice, tomato, eggplant, tobacco, and wheat were acquired from the Rice Genome Annotation Project (RGAP, http://rice.plantbiology.msu.edu/) and Ensemble database (http://plants.ensembl.org/index.html), reserved for subsequent analysis. Algorithm based on BLASTP was performed using the amino acid sequences of A. thaliana FBA proteins as queries in the protein databases of S. tuberosum, with an E-value ≤ 1e-5 and other parameters as default values, and redundant sequences were manually deleted. The Hidden Markov Model (HMM) of Glycolytic (PF00274) and fructose-bisphosphate aldolase class II domain (PF01116) were downloaded from InterPro database (https://www.ebi.ac.uk/interpro/search/sequence/) to identify all FBA proteins in the potato genome. To further verify the reliability of the candidate sequences, the structural domains of the candidate FBA protein sequences were identified in the Pfam database (http://pfam.xfam.org/) and SMART database (http://smart.embl-heidelberg.de). The website Cell-Ploc 2.0 (http://www.csbo.sjtu.edu.cn/bioinf/Cell-PLoc-2) was used to predict the subcellular localization of StFBA proteins. The online software ExPASy (https://web.expasy.org/protparam/) was used to analyze the physicochemical properties of the StFBA members.

2.3 Phylogenetic, gene structure and conserved motifs analyses of StFBA

The nine StFBA protein sequences were subjected to multiple sequence alignment using ClustaW software (Thompson et al., 2002) with default parameters. Neighbor-Joining (NJ) phylogenetic tree was constructed using MEGA11 software with 1,000 bootstrap replicates. The same method was used for constructing and analyzing the phylogenetic tree of potato FBA gene family and the FBA genes of eight in Arabidopsis (Lu et al., 2012), eight in tomato (Cai et al., 2016), seven in eggplant, 16 in tobacco (Zhao et al., 2021), seven in rice (Zhang et al., 2016) and 21 in wheat (Lv et al., 2017). The exon-intron structures were visualized using the Gene Structure Display Server (GSDS, http://gsds.cbi.pku.edu.cn) (Hu et al., 2015). The motif prediction was performed on the MEME program (http://meme-suite.org/) (Bailey et al., 2009), with the motif length ranging from 10 to 60 amino acid residues and maximum motif number set to 20, with all other parameters being default.

2.4 Synteny analysis and chromosome localization

TBtools software (Chen et al., 2020) was used to obtain information such as chromosome length and gene density in the potato genome and to draw the chromosome location maps. The synteny relationships of the FBA genes within and among different species were calculated using MCScanX, with a threshold of E < 1 × 10⁻⁵ (Wang et al., 2012).

2.5 Structure prediction of StFBA proteins

SOPMA software (http://npsa-pbil.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopma.html) was used to predict the secondary structure of StFBA protein. Tertiary structure prediction of StFBA proteins was performed using the latest artificial intelligence algorithm, AlphaFold2 (Jumper et al., 2021), which utilizes large-scale structural databases and deep-learning models that are able to provide more accurate protein structure prediction results. The tertiary structure visualization of StFBA proteins was performed with PyMOL2.5 software (https://pymol.org/2/).

2.6 Analysis of cis-acting elements in the promoter of StFBA genes

The sequence spanning 2000 bp upstream of the initiation codon of each StFBA gene was submitted to PlantCARE (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/) (Lescot et al., 2002) for predicting cis-elements. TBtools software was used to draw the cis-acting elements in the promoter region.

2.7 Abiotic stress treatments of potato plants

Two batches of potato plantlets with 30 days of healthy and uniform growth were selected. The middle parts of one batch of plantlets were cut into 1–1.5 cm stem segments with leaves and inoculated into glass tubes. They were first placed in dark conditions for cultivation and after 2 days, they were exposed to blue light at a wavelength of 460 nm, red light at a wavelength of 620 nm, and white light (control) by LEDs (Opt-run Biotechnology Co., Nanjing, China). Leaves and stolon of 35 days potato plantlets and 80 days tuber under different light conditions were then collected. Another batch of potato plantlets was transferred to hydroponic conditions for 1 week and then placed in 150 mM NaCl solution and 15% polyethylene glycol (PEG6000) solution for salt and drought stress treatments at 0, 3, 6, 12, and 24 h. Then, the leaves from the top of the potato plants were collected. All tissue samples collected were immediately placed in liquid nitrogen and stored at −80°C for later use, with each treatment group containing three biological replicates.

2.8 Quantitative real-time PCR analysis

RNA extraction and reverse transcription were referred to Huang et al. (2017), and qRT-PCR was referred to Jin et al. (2019) method. PCR was performed using an ABIViiA7 real-time PCR instrument (Life Technologies, United States). Specific amplification primers for the StFBA gene were designed using NCBI Primer-BLAST (https://www.ncbi.nlm.nih.gov/) with potato elongation factor-1alpha (EF1α) as an internal control (Supplementary Table S4) (Tang et al., 2017). The qRT-PCR for three independent biological replicate experiments, with three replicates for each PCR reaction. The relative expression of StFBA genes in different tissues and different treatments was calculated using the 2^−ΔΔCT and 2^−ΔΔCT method, respectively (Livak and Schmittgen, 2001). The one-way ANOVA of variance was used to conduct the statistical analyses of data results by SPSS software, and the significance of differences was analyzed using Duncan’s Multiple Range Test (DMRT) as a post hoc test (p < 0.05). GraphPad software was used for statistical analysis and drawing of data.

3 Results

3.1 Genome-wide identification of StFBA gene family members

A total of nine StFBA genes were identified from the potato genome using genome-wide analysis (Table 1). The StFBA genes were randomly and unevenly distributed across seven chromosomes; they were named StFBA1∼StFBA9 according to their order in the potato chromosome (Supplementary Figure S1). Analysis of the amino acid characteristics of StFBA family proteins showed that the number of amino acids (aa) of StFBA proteins ranged from 311 (StFBA3) to 1,379 (StFBA4), where the sequence length of the amino acid sequences of class II was significantly longer than that of class I. This difference may be due to the variations in their structure and function. The protein molecular weights ranged from 33,653.82 (StFBA3) to 148,195.01 (StFBA4) Da, while the theoretical isoelectric point varied between 5.94 (StFBA6) to 8.76 (StFBA5). Prediction of subcellular localizations of proteins showed that StFBA1∼StFBA5 and StFBA8 were localized in chloroplasts, and StFBA6, StFBA7 and StFBA9 were localized in the cytoplasm.

Table 1

Table 1. Profiles of the identified FBA gene family in S. tuberosum.

3.2 Phylogeny of StFBA and other plant FBA genes

To determine evolutionary relationships of the FBA gene family, 76 FBA protein amino acid sequences from seven species, including Arabidopsis, tomato, tobacco, rice, eggplant, wheat, and potato, were used to construct a phylogenetic tree (Supplementary Table S1). The results indicated that the StFBA gene, like other species, was classified into class I and II subclades (Figure 1). Class I could be further divided into four subclasses. Among the nine StFBA proteins, one was in Class I-a, two in Class I-b, one in Class I-c, four in Class I-d, and one in Class II. Of these, the identified wheat TaFBA19-TaFBA21 belonged to Class II, which was consistent with the clustering characteristics of StFBA4. The phylogenetic tree showed that FBA proteins were more closely related and similar between potato and other Solanaceae plants. In terms of evolutionary relationships, the FBA genes of potato and tomato were more closely related than those of other species, suggesting that they might have similar biological functions.

Figure 1

Figure 1. Phylogenetic tree of FBA proteins from potato and other plant species. Abbreviations represent the following species: S. tuberosum (St), A. thaliana (At), S. lycopersicum (Sl), S. melongena (Sm), N. tabacum (Nt), Oryza sativa (Os) and Triticum aestivum (Ta). Different colored squares represented different species. Different colored arcs indicate different classes or subclasses.

3.3 Gene structure and protein motif analyses of the StFBA gene family

Phylogenetic analysis indicated that FBA was divided into two classes in potatoes. Eight StFBA genes were identified in Class I and only StFBA4 was identified in Class II (Figure 2A). Analyses of gene structures revealed significant variations in exon count and sequence length among different subgroups of StFBA genes (Figure 2B). Class I had 3–6 exons, with the fewest exons found in StFBA7 and StFBA9, and the most in StFBA2 and StFBA5. Class II contained 42 exons, which is significantly more than Class I members, suggesting that StFBA4 might have a more complex structure, potentially including functional domains or regulatory elements. The motif distributions in FBA proteins were analyzed, and 20 conserved motifs were predicted (Supplementary Figure S2). Motif 1, motif 5, motif 6, and motif 10 were exclusive to Class I. StFBA1 and StFBA2 had 12 common motifs. StFBA6, StFBA7, and StFBA9 had 11 common motifs, with motif3 is a conserved motif specific to them. Moreover, some genes are missing more motifs, for example, StFBA8 contains only 6 motifs, suggesting that its function may be altered. Except for motif 9, the remaining eight motifs were specific to StFBA4 in Class II (Figure 2C). Overall, the type and distribution of conservative base sequences of Class I StFBA showed a high degree of uniformity, indicating that the structures of StFBA proteins were highly conserved in specific subclasses. In contrast, the gene structure and function of members of the potato FBA gene family might be varied between subclasses.

Figure 2

Figure 2. Phylogenetic relationships, gene structure, and conserved motifs of the StFBA gene family: (A) Phylogenetic relationships of the nine StFBA genes. (B) Gene structures of StFBA. UTR, untranslated region; CDS, coding sequence. (C) Distributions of the conserved motifs within StFBA proteins.

3.4 Tertiary structure of the StFBA proteins

Analyses of tertiary structures revealed that the StFBA proteins mainly consisted of α-helices, β-turns, irregular coils, and extended chains. The percentage of α-helix was the highest (42.77%–53.22%), followed by irregular coils (27.17%–33.76%), while extended chains and β-folds accounted for lower percentages (Supplementary Table S2). Comparison of the predicted tertiary structures of StFBA proteins revealed significant differences in the protein structures of different subclasses. The three-dimensional conformations of eight Class I StFBA gene family members were comparable (Figure 3). StFBA1, StFBA2, and StFBA5 were predicted to have more consistent protein structures. Compared to other members, the StFBA3 protein had a simpler structure, with fewer α-helices and β-turns, implying possible changes in its function. The remaining members had subtle differences in protein structures. StFBA4 had a larger number of amino acid sequences and more complex protein tertiary structures, suggesting that it might be involved in more intricate biological functions.

Figure 3

Figure 3. Predicted tertiary structures of FBA proteins in S. tuberosum.

3.5 Synteny analysis of FBA gene family members

To study the contraction and expansion of StFBA gene family members during evolution, the FBA genes of potato, tomato, and Arabidopsis were compared using collinearity analyses (Figure 4). The results showed that 11 StFBA members were found to have synteny with AtFBA, but StFBA2, StFBA3, and StFBA6 had no homologous genes in Arabidopsis. There were 14 pairs of syntenic FBA genes between potatoes and tomatoes, suggesting that the FBA genes were more homologous and conserved in Solanaceae than in dicotyledons. It is noteworthy that we have found that StFBA4 is homologous to AT1G18270 in Arabidopsis and Solyc03g118640.3.1 in tomato. Although AT1G18270 (NCBI Reference Sequence: NP_001117303.1) is annotated as ketose-bisphosphate aldolase class-II family protein, and Solyc03g118640.3.1 (NCBI Reference Sequence: XP_004235744.1) is annotated as uncharacterized protein LOC101261901. They all contain fructose-bisphosphate aldolase class-II (PF01116) domain and is further assigned as fructose-bisphosphate aldolase. Thus, AT1G18270 and Solyc03g118640.3.1 are actually both fructose-bisphosphate aldolase class II family proteins. Additionally, there were four pairs of duplicate segments, StFBA1-StFBA2, StFBA2-StFBA3, StFBA1-StFBA3, and StFBA1-StFBA8, in the potato genome region (Supplementary Figure S3). Expansion of potato FBA family members was mainly associated with segmental duplication events (Supplementary Table S3).

Figure 4

Figure 4. Collinear relationships of FBA genes in S. tuberosum, A. thaliana, and S. lycopersicum. The gray lines in the background indicate collinearity between potato genes and those of other species, and the green and yellow lines indicate collinear FBA gene pairs.

3.6 Expression patterns of StFBA genes in different tissues

To further understand the potential function of StFBA genes, the expression levels of nine StFBA genes in a variety of tissues were analyzed (Figure 5). It was found that some class I StFBA genes showed tissue specificities. For example, StFBA1 was highly expressed in leaves, StFBA2 was mainly expressed in green organs compared to other tissues, while StFBA5 was mainly expressed in floral organs. StFBA6 and StFBA9 were specifically expressed in petiole and had similar expression characteristics. StFBA3 was expressed at a lower level and was not expressed in many tissues. Most of the FBA genes were expressed in all tissues, indicating that they might be involved in the growth and development of potatoes. Specifically, StFBA7 was highly expressed in floral organs, tubers, and stolon, suggesting it might serve a pivotal function during flowering and tuber formation in potatoes. The class II StFBA4 was expressed in most tissues, but at a lower expression level.

Figure 5

Figure 5. Expression analyses of StFBA genes in different tissues. Significant differences among the groups were compared based on Duncan’s test (p < 0.05). The data points represent the mean ± SD.

3.7 Expression patterns of StFBA genes in response to different light

Multiple cis-regulatory elements were detected in the promoter region of StFBA genes (Figure 6). Based on functions of the cis-acting elements, the selected 20 elements were categorized into four groups: light response, stress response, hormone response and plant growth and development. Hormonal response elements included responses to auxin, abscisic acid, gibberellin, and salicylic acid; and stress included responses to drought and cold stress elements. We found that multiple light-responsive elements existed in the promoter regions of each StFBA gene, indicating that the StFBA genes were an important component of light responses in potatoes. Thus, the expression levels of StFBA in leaves, stolon, and microtubers were analyzed under different light treatments. The expression patterns of FBA genes exhibited discrepancies under different spectra. Due to monochromatic red light exposure, the leaves development of the microtubers was severely inhibited, hence there was no leaf sample under red light (Figure 7A). Under blue light, StFBA7 showed lower expression in leaves, but significantly increased expression in stolon and microtubers. Under red light, StFBA3, StFBA7, StFBA8, and StFBA9 showed higher expressions, when compared to under white light in stolon. In contrast, the expression levels of StFBA1, StFBA2, StFBA3, and StFBA9 were significantly higher than that of the other treatments in tubers. Class II StFBA4 gene showed lower expression in all tissues under red and blue light, suggesting that the gene might be spectrally insensitive. Moreover, the expression levels of StFBA1 and StFBA2 in all tissues were highly similar under various spectral treatments, suggesting that these genes might have similar regulatory mechanisms in response to different lighting conditions (Figure 7B).

Figure 6

Figure 6. Distribution of cis-elements in the promoter regions of StFBA genes. The heat map represents classifications and statistics of cis-acting components, the cis-elements were divided into four broad categories. Different colors represent different types of elements, and the ruler at the bottom indicates the direction and length of the sequence.

Figure 7

Figure 7. Growth of potato plantlet under different light spectrum and expressions of nine StFBA genes in different tissues under red and blue light treatments: (A) Growth of potato plantlet under different light spectrum. Blue, 460 nm blue light; red, 620 nm red light; control, white light. (B) The qRT-PCR was performed to analyze the relative expression levels of StFBA genes. Error bars represent standard deviations of the means from three independent experiments. Significant differences among the groups were compared based on Duncan’s test (p < 0.05). The data points represented the mean ± SD.

3.8 Expression profiles of the StFBA genes in response to abiotic stresses

The expression patterns of the StFBA gene ins response to abiotic stress such as drought and salt were determined, indicating that the expression levels of the nine StFBA genes varied to different degrees under different stress treatments. Under the 15% PEG6000 treatment, StFBA3 mainly showed a decreasing trend, while expression levels of StFBA5 and StFBA7 were significantly upregulated, so they might be important candidate genes involved in drought stress. StFBA4 expression was minimally affected. The expression levels of StFBA5 and StFBA7 were significantly upregulated (p < 0.01). The expression patterns of the remaining genes were more complicated. For example, StFBA1, StFBA2, and StFBA6 were highly induced at 3 h, then decreased at 6 h, then finally increased again at 12 h treatment, StFBA8 and StFBA9 expression levels showed a tendency to decrease and then increase, and all of these genes gradually recovered to their untreated expression levels at 24 h (Figure 8). Under salt conditions (150 mM NaCl), expression of StFBA3 significantly decreased. The expressions of StFBA1, StFBA2, StFBA4, and StFBA6 peaked at 3 h, while the expressions of StFBA5, StFBA7, StFBA8, and StFBA9 peaked at 6 h, then all of them showed a decreasing trend (Figure 9). Together, the results further confirmed the possible role of the StFBA gene family in potato responses to various abiotic stresses.

Figure 8

Figure 8. Expression profiles of the StFBA genes after 15% PEG6000 simulated drought stress. PEG, polyethylene glycol. Significant differences among the groups were compared based on Duncan’s test (p < 0.05). The data points represented the mean ± SD.

Figure 9

Figure 9. Expression profiles of the StFBA genes under salt stress treatment using 150 mM NaCl. Significant differences among the groups were compared based on Duncan’s test (p < 0.05). The data points represent the mean ± SD.

4 Discussion

Organic compounds produced by photosynthesis can promote the formation and enlargement of potatoes (Fernie and Willmitzer, 2001). As a key enzyme in the Calvin cycle and glycolysis, FBA participates in the metabolic process of starch biosynthesis. Decreasing the activity of this enzyme can result in stunted potato growth, reduced photosynthetic efficiency, and an impact on tuber starch content, ultimately affecting potato yield and quality (Kossmann et al., 1995). Hence, FBA plays a vital regulatory role in the growth and tuber yield of potatoes. Our study identified nine FBA genes in the complete genome sequence of potatoes, which were divided into Class I and Class II members, based on their two structural domains (Table 1). Potatoes, like Arabidopsis and tomatoes, possess eight Class I members, indicating the conservation of Class I StFBA genes during their evolution. Class I and Class II FBA have no similarities in gene sequence, protein structure, and function, and are often considered to have evolved from different origins (Schnarrenberger et al., 1990; Zgiby et al., 2000). We identified one Class II FBA gene in S. tuberosum, which has 1,379 amino acids and 42 exons, a significantly higher number than Class I. This finding aligns with observations in other plant species, three have been identified in Triticum aestivum (Lv et al., 2017), two in B. napus (Zhao et al., 2019), and four in cotton (Li et al., 2021). The amino acid numbers of these genes ranged from 1,352 to 1,383, with exon numbers between 41 and 42. This indicates that Class II FBAs are conserved among different species and might have similar functions. According to the conserved motif analysis results of StFBA, we found that apart from StFBA1 and StFBA2, which contain all 10 motifs, other members exhibited a certain degree of motif loss, with the most loss in StFBA8, suggesting possible functional changes in this gene. Class II FBA shared fewer motifs in common with Class I (Figure 2). Overall, members within the same subclass showed similarities in sequence, gene structure, and conserved motif composition, indicating that StFBA family members may have undergone similar genetic evolution processes.

Based on the phylogenetic tree, eight Class I StFBA members were further divided into four subclasses (Figure 1), the Class I-d subclass had the most clustered StFBA genes. Closely related gene members in the phylogenetic tree may have similar structures and functions. Because both S. tuberosum and S. lycopersicum belong to the Solanaceae plant family, they were closer in the evolutionary tree. It is worth noting that SlFBA7 and StFBA7 clustered together, had the same number of amino acids, and were located in the cytoplasm. Overexpression of SlFBA7 can significantly increase the FBA activity and the expression level of other enzymes in the Calvin cycle of tomato seedling leaves, along with an increase in regeneration of RuBP (Cai et al., 2016). Therefore, we speculate that StFBA7 may play a pivotal role in enhancing the photosynthetic efficiency of potatoes. At present, the functions of some FBA genes have been thoroughly studied in model plants such as A. thaliana and O. sativa (Lu et al., 2012; Simkin et al., 2017). For example, AtFBA2 has the highest activity and expression level in leaves, and overexpression of AtFBA2 boosts the photosynthetic capacity and biomass of the plants. However, fba2 mutant plants are dwarfed and the starch level is significantly reduced. AtFBA3 is mainly expressed in the vascular tissues of leaves and roots, so the absence of AtFBA3 limits the biosynthesis of essential compounds such as amino acids. Similarly, the fba3 mutants grow more slowly and show a decrease in photosynthesis (Carrera et al., 2021).

Collinearity analysis between S. tuberosum and A. thaliana revealed that AtFBA2 was homologous with StFBA1 and StFBA8, while AtFBA3 shared homology with StFBA5 (Figure 4). We speculate that StFBA1, StFBA5, and StFBA8 might contribute to enhancing tuber yield. However, their specific functions and mechanisms of action still require further confirmation. In plants, gene duplication events primarily consist of tandem duplication and segmental duplication, which play key roles in biological evolution, gene family formation, and the transmission of genetic information (Kong et al., 2007; Paterson et al., 2010). The members of the StFBA gene family were closely associated with segmental duplication events. The multiple occurrences of StFBA1 in repetitive events suggested that this gene may have played a significant role in driving the evolution of the FBA gene family, leading to gene amplifications and functional changes. The tertiary structure of proteins determines their function in cells and how they interact with other molecules (Kuhlman and Bradley, 2019). Based on the physicochemical properties and predicted protein structure of StFBA, we showed that the structures of StFBA3 and StFBA8 were simpler than other members (Figure 3). These two proteins exhibited the lowest number of amino acids and the smallest percentage of α-helices (Supplementary Table S2). It is therefore possible that during the process of segmental duplication, these two genes underwent loss or alterations in some of their conserved domains and functions. StFBA6, StFBA7, and StFBA9 were located in the cytoplasm, had very similar amino acid numbers, and showed a high degree of similarity in their secondary and tertiary structures, indicating that they may be homologous genes.

Through the analysis of transcriptome data and qRT-PCR validation, we obtained the expression levels of StFBA genes in different potato tissues. StFBA7 exhibited higher expressions in tubers, stolon, and flower organs, indicating that it may be a key gene regulating the processes of microtuber formation and expansion (Figure 5). The expression patterns of StFBA family members in different tissues demonstrated that FBA genes had specific roles during potato growth and development. The cis-acting element analysis revealed an abundance of light-responsive elements within the StFBA gene family (Figure 6). Throughout the growth and development of potatoes, light influenced tuber induction and expansion via the processes of the photoperiod, light intensity, and spectrum. Among them, light spectrum is one of the key signals regulating tuberization in potato plants. It has been reported that the transcription level of StBEL5, a primary signal for potato tuberization, was mostly affected by the light spectrum, with the red-blue spectrum inducing the initiation of StBEL5 synthesis (Chatterjee et al., 2007). Further studies found that blue light at 460 nm wavelength aided tuber expansion, while 620 nm red light was more favorable for tuber induction (Li et al., 2020). The microtubers under blue light were significantly larger than other spectra in our study (Supplementary Figure S4). The qRT-PCR results indicated that StFBA genes were highly responsive to spectral changes, with StFBA3, StFBA8, and StFBA9 showing significantly higher expression levels in leaf, stolon, and tuber, when compared to the white light control. These results implied that they may improve photosynthesis and play an important function in regulating the induction and expansion of microtubers (Figure 7).

Potato is highly sensitive to drought and salt stress. Drought affects the aboveground plant height and leaf area index, leading to a reduction in photosynthetic efficiency and adversely affecting potato yield and quality through shortened growth cycles and decreased underground tuber quantity and size (Deblonde and Jean-François, 2001). High salt stress disrupts the source-sink balance in potatoes, impairing the production of antioxidants and antioxidant enzyme activity, and accelerating the occurrence of plant diseases (Albacete et al., 2014; Chourasia et al., 2021). Studies have shown that FBA genes played a role in responses to various abiotic stresses, including drought, high salt (Yamada et al., 2000), and cold stress (Liu et al., 2019). We found that under drought stress, the expression of StFBA5 significantly increased, and it maintained a high expression level even after 24 h (Figure 8). When exposed to high salt stress, the expression levels of most StFBA genes showed an initial increase followed by a sharp decline (Figure 9). These results indicated that StFBA members played a crucial role in responding to different environmental stress. In addition, FBA has been found to have a wide range of functions in other species. For example, AtFBA6 participates in the memory of high temperature stress in the apical meristem of Arabidopsis (Olas et al., 2021). AtFBA8 responds to environmental humidity changes and interacts with actin, involving guard cell opening and closing (Van der Linde et al., 2011). The soybean gmfbacgmfbac2 double mutant exhibits slow growth, narrowed leaves, dwarfed seedlings, and severe imbalance in carbon-nitrogen metabolism (Qiu et al., 2023). It is hypothesized that the StFBA gene may have functional diversity, but the specific functional mechanism needs to be verified through future experiments.

5 Conclusion

In this study, eight class I and one class II StFBA gene members were identified using genome-wide analyses. We found that the FBA of potatoes and tomatoes were closely related in the evolutionary tree and in synteny relationships. StFBA1 was a key gene in the fragment duplication event within the gene family. StFBA6, StFBA7, and StFBA9 were involved in growth and development of potatoes. StFBA3, StFBA8, and StFBA9 may influence photosynthesis and play a critical function in regulating the induction and expansion of microtubers, while StFBA5 and StFBA9 played essential roles in responses to drought and salt stress. The results of this study provided a basis for understanding the potential functions of the StFBA gene family, and presented ideas for the improvement of photosynthetic efficiency and response to adversity stress of potatoes.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author contributions

TL: Data curation, Visualization, Writing–original draft. XnH: Visualization, Writing–original draft, Data curation. ZS: Data curation, Writing–original draft. BM: Data curation, Writing–original draft. XW: Data curation, Writing–original draft. TF: Supervision, Writing–review and editing. HA: Supervision, Writing–review and editing. XaH: Conceptualization, Project administration, Supervision, Writing–review and editing. RL: Conceptualization, Funding acquisition, Project administration, Supervision, Writing–review and editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research was funded by the Natural Science Foundation of Anhui (No. 2108085QC125), the Key Project of Anhui Provincial Education Department (Nos 2022AH051627 and 2023AH051886), the Talent Introduction Start-up Fund Project of Anhui Science and Technology University (No. NXYJ202102), the National Innovation and Entrepreneurship Training Program for College Students (No. 202310879007), and the Provincial Graduate Full English/Bilingual Demonstration Course (No. 2022qyw/sysfkc040). The funding institution was not involved in the design of the study, collection, analysis, and interpretation of data, and in writing the manuscript.

Acknowledgments

We thank International Science Editing (http://www.internationalscienceediting.com) for editing this manuscript.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

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.2024.1364944/full#supplementary-material

References

Albacete, A. A., Martínez-Andújar, C., and Pérez-Alfocea, F. (2014). Hormonal and metabolic regulation of source-sink relations under salinity and drought: from plant survival to crop yield stability. Biotechnol. Adv. 32, 12–30. doi:10.1016/j.biotechadv.2013.10.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Ashida, H., Mizohata, E., and Yokota, A. (2019). Learning RuBisCO's birth and subsequent environmental adaptation. Biochem. Soc. Trans. 47, 179–185. doi:10.1042/BST20180449

PubMed Abstract | CrossRef Full Text | Google Scholar

Bailey, T. L., Boden, M., Buske, F. A., Frith, M., Grant, C. E., Clementi, L., et al. (2009). MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 37, W202–W208. doi:10.1093/nar/gkp335

PubMed Abstract | CrossRef Full Text | Google Scholar

Blom, N., and Sygusch, J. (1997). Product binding and role of the C-terminal region in class I D-fructose 1, 6-bisphosphate aldolase. Nat. Struct. Biol. 4, 36–39. doi:10.1038/nsb0197-36

PubMed Abstract | CrossRef Full Text | Google Scholar

Bracher, A., Whitney, S. M., Hartl, F. U., and Hayer-Hartl, M. (2017). Biogenesis and metabolic maintenance of Rubisco. Annu. Rev. Plant Biol. 68, 29–60. doi:10.1146/annurev-arplant-043015-111633

PubMed Abstract | CrossRef Full Text | Google Scholar

Bryant, B. A. (2000). The production and analysis of tobacco plants with reduced levels of the Calvin cycle enzyme sedoheptulose-1, 7-bisphosphatase. Doctoral dissertation, England: University of Essex.

Google Scholar

Cai, B., Li, Q., Liu, F., Bi, H., and Ai, X. (2018). Decreasing fructose-1,6-bisphosphate aldolase activity reduces plant growth and tolerance to chilling stress in tomato seedlings. Physiol. Plant 163, 247–258. doi:10.1111/ppl.12682

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, B., Li, Q., Xu, Y., Yang, L., Bi, H., and Ai, X. (2016). Genome-wide analysis of the fructose 1, 6-bisphosphate aldolase (FBA) gene family and functional characterization of FBA7 in tomato. Plant Physiol. Biochem. 108, 251–265. doi:10.1016/j.plaphy.2016.07.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Capodagli, G. C., Sedhom, W. G., Jackson, M., Ahrendt, K. A., and Pegan, S. D. (2014). A noncompetitive inhibitor for Mycobacterium tuberculosis's class IIa fructose 1, 6-bisphosphate aldolase. Biochemistry 53, 202–213. doi:10.1021/bi401022b

PubMed Abstract | CrossRef Full Text | Google Scholar

Carrera, D. Á., George, G. M., Fischer-Stettler, M., Galbier, F., Eicke, S., Truernit, E., et al. (2021). Distinct plastid fructose bisphosphate aldolases function in photosynthetic and non-photosynthetic metabolism in Arabidopsis. J. Exp. Bot. 72, 3739–3755. doi:10.1093/jxb/erab099

PubMed Abstract | CrossRef Full Text | Google Scholar

Chatterjee, M., Banerjee, A. K., and Hannapel, D. J. (2007). A BELL1-like gene of potato is light activated and wound inducible. Plant Physiol. 145, 1435–1443. doi:10.1104/pp.107.105924

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, C., Chen, H., Zhang, Y., Thomas, H. R., Frank, M. H., He, Y. H., et al. (2020). TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 13, 1194–1202. doi:10.1016/j.molp.2020.06.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Chourasia, K. N., Lal, M. K., Tiwari, R. K., Dev, D., Kardile, H. B., Patil, V. U., et al. (2021). Salinity stress in potato: understanding physiological, biochemical and molecular responses. Life (Basel). 11, 545. doi:10.3390/life11060545

PubMed Abstract | CrossRef Full Text | Google Scholar

Deblonde, P. M. K., and Jean-François, L. (2001). Effects of moderate drought conditions on green leaf number, stem height, leaf length and tuber yield of potato cultivars. Eur. J. Agron. 14, 31–41. doi:10.1016/S1161-0301(00)00081-2

CrossRef Full Text | Google Scholar

Fan, W., Zhang, Z., and Zhang, Y. (2009). Cloning and molecular characterization of fructose-1, 6-bisphosphate aldolase gene regulated by high-salinity and drought in Sesuvium portulacastrum. Plant Cell Rep. 28, 975–984. doi:10.1007/s00299-009-0702-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Fernie, A. R., and Willmitzer, L. (2001). Molecular and biochemical triggers of potato tuber development. Plant Physiol. 127, 1459–1465. doi:10.1104/pp.127.4.1459

PubMed Abstract | CrossRef Full Text | Google Scholar

Flechner, A., Gross, W., Martin, W. F., and Schnarrenberger, C. (1999). Chloroplast class I and class II aldolases are bifunctional for fructose-1, 6-biphosphate and sedoheptulose-1, 7-biphosphate cleavage in the Calvin cycle. FEBS Lett. 447, 200–202. doi:10.1016/s0014-5793(99)00285-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Frenkel, Z. M., and Trifonov, E. N. (2005). Closed loops of TIM barrel protein fold. J. Biomol. Struct. Dyn. 22, 643–656. doi:10.1080/07391102.2005.10507032

PubMed Abstract | CrossRef Full Text | Google Scholar

Galkin, A., Kulakova, L., Melamud, E., Li, L., Wu, C., Mariano, P., et al. (2007). Characterization, kinetics, and crystal structures of fructose-1, 6-bisphosphate aldolase from the human parasite, Giardia lamblia. J. Biol. Chem. 282, 4859–4867. doi:10.1074/jbc.M609534200

PubMed Abstract | CrossRef Full Text | Google Scholar

Gross, W., Lenze, D., Nowitzki, U., Weiske, J., and Schnarrenberger, C. (1999). Characterization, cloning, and evolutionary history of the chloroplast and cytosolic class I aldolases of the red alga Galdieria sulphuraria. Gene 230, 7–14. doi:10.1016/s0378-1119(99)00059-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Haake, V., Geiger, M., Walch-Liu, P., Of Engels, C., Zrenner, R., and Stitt, M. (1999). Changes in aldolase activity in wild-type potato plants are important for acclimation to growth irradiance and carbon dioxide concentration, because plastid aldolase exerts control over the ambient rate of photosynthesis across a range of growth conditions. Plant J. 17, 479–489. doi:10.1046/j.1365-313X.1999.00391.x

CrossRef Full Text | Google Scholar

Haake, V., Zrenner, R., Sonnewald, U., and Stitt, M. (1998). A moderate decrease of plastid aldolase activity inhibits photosynthesis, alters the levels of sugars and starch, and inhibits growth of potato plants. Plant J. 14, 147–157. doi:10.1046/j.1365-313x.1998.00089.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Henkes, S., Sonnewald, U., Badur, R., Flachmann, R., and Stitt, M. (2001). A small decrease of plastid transketolase activity in antisense tobacco transformants has dramatic effects on photosynthesis and phenylpropanoid metabolism. Plant Cell 13, 535–551. doi:10.1105/tpc.13.3.535

PubMed Abstract | CrossRef Full Text | Google Scholar

Hester, G., Brenner-Holzach, O., Rossi, F. A., Struck-Donatz, M., WinterhalterWinterhalter, K. H., Smit, J. D., et al. (1991). The crystal structure of fructose-1, 6-bisphosphate aldolase from Drosophila melanogaster at 2.5 A˚ resolution. FEBS Lett. 292, 237–242. doi:10.1016/0014-5793(91)80875-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, B., Jin, J., Guo, A. Y., Zhang, H., Luo, J., and Gao, G. (2015). GSDS 2.0: an upgraded gene feature visualization server. Bioinformatics 31, 1296–1297. doi:10.1093/bioinformatics/btu817

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, X., Yang, L., Jin, Y., Lin, J., and Liu, F. (2017). Generation, annotation, and analysis of a large-scale expressed sequence tag library from Arabidopsis pumila to explore salt-responsive genes. Front. Plant Sci. 8, 955. doi:10.3389/fpls.2017.00955

PubMed Abstract | CrossRef Full Text | Google Scholar

Izard, T., and Sygusch, J. (2004). Induced fit movements and metal cofactor selectivity of class II aldolases: structure of Thermus aquaticus fructose-1, 6-bisphosphate aldolase. J. Biol. Chem. 279, 11825–11833. doi:10.1074/jbc.M311375200

PubMed Abstract | CrossRef Full Text | Google Scholar

Jacquot, J. P., Lopez-Jaramillo, J., Chueca, A., Cherfils, J., Lemaire, S., Chedozeau, B., et al. (1995). High-level expression of recombinant pea chloroplast fructose-1, 6-bisphosphatase and mutagenesis of its regulatory site. Eur. J. Biochem. 229, 675–681. doi:10.1111/j.1432-1033.1995.tb20513.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, Y., Liu, F., Huang, W., Sun, Q., and Huang, X. (2019). Identification of reliable reference genes for qRT-PCR in the ephemeral plant Arabidopsis pumila based on full-length transcriptome data. Sci. Rep. 9, 8408. doi:10.1038/s41598-019-44849-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Jordan, D. B., and Ogren, W. L. (1984). The CO₂/O₂ specificity of ribulose 1, 5-bisphosphate carboxylase/oxygenase: dependence on ribulose bisphosphate concentration, pH and temperature. Planta 161, 308–313. doi:10.1007/BF00398720

PubMed Abstract | CrossRef Full Text | Google Scholar

Jumper, J., Evans, R., Pritzel, A., Green, T., Figurnov, M., Ronneberger, O., et al. (2021). Highly accurate protein structure prediction with AlphaFold. Nature 596, 583–589. doi:10.1038/s41586-021-03819-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Kong, H., Landherr, L. L., Frohlich, M. W., Leebens-Mack, J., Ma, H., and DePamphilis, C. W. (2007). Patterns of gene duplication in the plant SKP1 gene family in angiosperms: evidence for multiple mechanisms of rapid gene birth. Plant J. 50, 873–885. doi:10.1111/j.1365-313X.2007.03097.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Konishi, H., Yamane, H., Maeshima, M., and Komatsu, S. (2004). Characterization of fructose-bisphosphate aldolase regulated by gibberellin in roots of rice seedling. Plant Mol. Biol. 56, 839–848. doi:10.1007/s11103-004-5920-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Kossmann, J., Sonnewald, U., and Willmitzer, L. (1995). Reduction of the chloroplastic fructose-1,6-bisphosphatase in transgenic potato plants impairs photosynthesis and plant growth. Plant J. 6, 637–650. doi:10.1046/j.1365-313X.1994.6050637.x

CrossRef Full Text | Google Scholar

Kuhlman, B., and Bradley, P. (2019). Advances in protein structure prediction and design. Nat. Rev. Mol. Cell Biol. 20, 681–697. doi:10.1038/s41580-019-0163-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Lebherz, H. G., Leadbetter, M. M., and Bradshaw, R. A. (1984). Isolation and characterization of the cytosolic and chloroplast forms of spinach leaf fructose diphosphate aldolase. J. Biol. Chem. 259, 1011–1017. doi:10.1016/S0021-9258(17)43558-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Lefebvre, S., Lawson, T., Zakhleniuk, O. V., Lloyd, J. C., Raines, C. A., Fryer, M., et al. (2005). Increased sedoheptulose-1, 7-bisphosphatase activity in transgenic tobacco plants stimulates photosynthesis and growth from an early stage in development. Plant Physiol. 138, 451–460. doi:10.1104/pp.104.055046

PubMed Abstract | CrossRef Full Text | Google Scholar

Lescot, M., Déhais, P., Thijs, G., Marchal, K., Moreau, Y., Van de Peer, Y., et al. (2002). PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 30, 325–327. doi:10.1093/nar/30.1.325

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, R. N., You, J., Miao, C., Kong, L., Long, J. H., Yan, Y. Z., et al. (2020). Monochromatic lights regulate the formation, growth, and dormancy of in vitro-grown Solanum tuberosum L. microtubers. Sci. Hortic. 261, 108947. doi:10.1016/j.scienta.2019.108947

CrossRef Full Text | Google Scholar

Li, Z. Q., Zhang, Y., Li, H., Su, T. T., Liu, C. G., Han, Z. C., et al. (2021). Genome-wide characterization and expression analysis provide basis to the biological function of cotton FBA genes. Front. Plant Sci. 12, 696698. doi:10.3389/fpls.2021.696698

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, X., Fu, L., Qin, P., Sun, Y., Liu, J., and Wang, X. (2019). Overexpression of the wheat trehalose 6-phosphate synthase 11 gene enhances cold tolerance in Arabidopsis thaliana. Gene 710, 210–217. doi:10.1016/j.gene.2019.06.006

PubMed Abstract | 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, 402–408. doi:10.1006/meth.2001.1262

PubMed Abstract | CrossRef Full Text | Google Scholar

Long, S. P., Ainsworth, E. A., Rogers, A., and Ort, D. R. (2004). Rising atmospheric carbon dioxide: plants FACE the future. Annu. Rev. Plant Biol. 55, 591–628. doi:10.1146/annurev.arplant.55.031903.141610

PubMed Abstract | CrossRef Full Text | Google Scholar

Lorentzen, E., Pohl, E., Zwart, P., Stark, A., Russell, R. B., Knura, T., et al. (2003). Crystal structure of an archaeal class I aldolase and the evolution of (betaalpha)8 barrel proteins. J. Biol. Chem. 278, 47253–47260. doi:10.1074/jbc.M305922200

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, W., Tang, X., Huo, Y., Xu, R., Qi, S., Huang, J., et al. (2012). Identification and characterization of fructose 1, 6-bisphosphate aldolase genes in Arabidopsis reveal a gene family with diverse responses to abiotic stresses. Gene 503, 65–74. doi:10.1016/j.gene.2012.04.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Lv, G. Y., Guo, X. G., Xie, L. P., Xie, C. G., Zhang, X. H., Yang, Y., et al. (2017). Molecular characterization, gene evolution, and expression analysis of the fructose-1, 6-bisphosphate aldolase (FBA) gene family in wheat (Triticum aestivum L.). Front. Plant Sci. 8, 1030. doi:10.3389/fpls.2017.01030

PubMed Abstract | CrossRef Full Text | Google Scholar

Marsh, J. J., and Lebherz, H. G. (1992). Fructose-bisphosphate aldolases: an evolutionary history. Trends Biochem. Sci. 17, 110–113. doi:10.1016/0968-0004(92)90247-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Mutuku, J. M., and Nose, A. (2012). Changes in the contents of metabolites and enzyme activities in rice plants responding to Rhizoctonia solani kuhn infection: activation of glycolysis and connection to phenylpropanoid pathway. Plant Cell Physio 53, 1017–1032. doi:10.1093/pcp/pcs047

PubMed Abstract | CrossRef Full Text | Google Scholar

Obiadalla-Ali, H., Fernie, A. R., Lytovchenko, A., Kossmann, J., and Lloyd, J. R. (2004). Inhibition of chloroplastic fructose 1, 6-bisphosphatase in tomato fruits leads to decreased fruit size, but only small changes in carbohydrate metabolism. Planta 219, 533–540. doi:10.1007/s00425-004-1257-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Oelze, M. L., Muthuramalingam, M., Vogel, M. O., and Dietz, K. J. (2014). The link between transcript regulation and de novo protein synthesis in the retrograde high light acclimation response of Arabidopsis thaliana. BMC Genomics 15, 320. doi:10.1186/1471-2164-15-320

PubMed Abstract | CrossRef Full Text | Google Scholar

Olas, J. J., Apelt, F., Annunziata, M. G., John, S., Richard, S. I., Gupta, S., et al. (2021). Primary carbohydrate metabolism genes participate in heat-stress memory at the shoot apical meristem of Arabidopsis thaliana. Mol. Plant 14, 1508–1524. doi:10.1016/j.molp.2021.05.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Osakabe, Y., Maruyama, K., Seki, M., Satou, M., Shinozaki, K., and Yamaguchi-Shinozaki, K. (2005). Leucine-rich repeat receptor-like kinase1 is a key membrane-bound regulator of abscisic acid early signaling in Arabidopsis. Plant Cell 17, 1105–1119. doi:10.1105/tpc.104.027474

PubMed Abstract | CrossRef Full Text | Google Scholar

Paterson, A. H., Freeling, M., Tang, H., and Wang, X. (2010). Insights from the comparison of plant genome sequences. Annu. Rev. Plant Biol. 61, 349–372. doi:10.1146/annurev-arplant-042809-112235

PubMed Abstract | CrossRef Full Text | Google Scholar

Pelzer-Reith, B., Penger, A., and Schnarrenberger, C. (1993). Plant aldolase: cDNA and deduced amino-acid sequences of the chloroplast and cytosol enzyme from spinach. Plant Mol. Biol. 21, 331–340. doi:10.1007/BF00019948

PubMed Abstract | CrossRef Full Text | Google Scholar

Perham, R. N. (1990). The fructose-1, 6-bisphosphate aldolases: same reaction, different enzymes. Biochem. Soc. Trans. 18, 185–187. doi:10.1042/bst0180185

PubMed Abstract | CrossRef Full Text | Google Scholar

Pham, G. M., Hamilton, J. P., Wood, J. C., Burke, J. T., Zhao, H. N., Vaillancourt, B., et al. (2020). Construction of a chromosome-scale long-read reference genome assembly for potato. Gigascience 9, giaa100. doi:10.1093/gigascience/giaa100

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiu, Z., Bai, M., Kuang, H., Wang, X., Yu, X., Zhong, X., et al. (2023). Cytosolic fructose-1, 6-bisphosphate aldolases modulate primary metabolism and phytohormone homeostasis in soybean. Agronomy 13, 1383. doi:10.3390/agronomy13051383

CrossRef Full Text | Google Scholar

Raines, C. A. (2003). The Calvin cycle revisited. Photosynth Res. 75, 1–10. doi:10.1023/A:1022421515027

PubMed Abstract | CrossRef Full Text | Google Scholar

Schnarrenberger, C., Jacobshagen, S., Müller, B., and Krüger, I. (1990). Evolution of isozymes of sugar phosphate metabolism in green algae. Prog. Clin. Biol. Res. 344, 743–764.

PubMed Abstract | Google Scholar

Simkin, A. J., Lopez-Calcagno, P. E., Davey, P. A., Headland, L. R., Lawson, T., Timm, S., et al. (2017). Simultaneous stimulation of sedoheptulose 1, 7-bisphosphatase, fructose 1, 6-bisphophate aldolase and the photorespiratory glycine decarboxylase-H protein increases CO₂ assimilation, vegetative biomass and seed yield in Arabidopsis. Plant Biotechnol. J. 15, 805–816. doi:10.1111/pbi.12676

PubMed Abstract | CrossRef Full Text | Google Scholar

Sonnewald, U., Lerchl, J., Zrenner, R., and Frommer, W. (1994). Manipulation of sink-source relations in transgenic plants. Plant Cell Environ. 17, 649–658. doi:10.1111/j.1365-3040.1994.tb00156.x

CrossRef Full Text | Google Scholar

Tang, X., Zhang, N., Si, H., and Calderón-Urrea, A. (2017). Selection and validation of reference genes for RT-qPCR analysis in potato under abiotic stress. Plant Methods 13, 85. doi:10.1186/s13007-017-0238-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Thompson, J. D., Gibson, T. J., and Higgins, D. G. (2002). Multiple sequence alignment using ClustalW and ClustalX. Curr. Protoc. Bioinforma. 2, Unit 2.3. Chapter. doi:10.1002/0471250953.bi0203s00

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsutsumi, K., Kagaya, Y., Hidaka, S., Suzuki, J., Tokairin, Y., Hirai, T., et al. (1994). Structural analysis of the chloroplastic and cytoplasmic aldolase-encoding genes implicated the occurrence of multiple loci in rice. Gene 141, 215–220. doi:10.1016/0378-1119(94)90574-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Uematsu, K., Suzuki, N., Iwamae, T., Inui, M., and Yukawa, H. (2012). Increased fructose 1, 6-bisphosphate aldolase in plastids enhances growth and photosynthesis of tobacco plants. J. Exp. Bot. 63, 3001–3009. doi:10.1093/jxb/ers004

PubMed Abstract | CrossRef Full Text | Google Scholar

Van der Linde, K., Gutsche, N., Leffers, H. M., Lindermayr, C., Müller, B., Holtgrefe, S., et al. (2011). Regulation of plant cytosolic aldolase functions by redox-modifications. Plant Physiol. Biochem. 49, 946–957. doi:10.1016/j.plaphy.2011.06.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Tang, H., Debarry, J. D., Tan, X., Li, J., Wang, X., et al. (2012). MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 40, e49. doi:10.1093/nar/gkr1293

PubMed Abstract | CrossRef Full Text | Google Scholar

Wildman, S. G. (2002). Along the trail from Fraction I protein to Rubisco (ribulose bisphosphate carboxylase-oxygenase). Photosynth Res. 73, 243–250. doi:10.1023/A:1020467601966

PubMed Abstract | CrossRef Full Text | Google Scholar

Yamada, S., Komori, T., Hashimoto, A., Kuwata, S., Imaseki, H., and Kubo, T. (2000). Differential expression of plastidic aldolase genes in Nicotiana plants under salt stress. Plant Sci. 154, 61–69. doi:10.1016/s0168-9452(00)00188-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng, Y., Tan, X., Zhang, L., Jiang, N., and Cao, H. (2014). Identification and expression of fructose-1, 6-bisphosphate aldolase genes and their relations to oil content in developing seeds of tea oil tree (Camellia oleifera). PLoS One. 9, e107422. doi:10.1371/journal.pone.0107422

PubMed Abstract | CrossRef Full Text | Google Scholar

Zgiby, S. M., Thomson, G. J., Qamar, S., and Berry, A. (2000). Exploring substrate binding and discrimination in fructose1, 6-bisphosphate and tagatose 1, 6-bisphosphate aldolases. Eur. J. Biochem. 267, 1858–1868. doi:10.1046/j.1432-1327.2000.01191.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, C., Yang, Z., Tang, D., Zhu, Y., Wang, P., Li, D., et al. (2021). Genome design of hybrid potato. Cell 184, 3873–3883.e12. doi:10.1016/j.cell.2021.06.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, F., Zhang, P., Zhang, Y., Wang, S., Qu, L., Liu, X., et al. (2016). Identification of a peroxisomal-targeted aldolase involved in chlorophyll biosynthesis and sugar metabolism in rice. Plant Sci. 250, 205–215. doi:10.1016/j.plantsci.2016.06.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, H., Xu, F., Wu, Y., Hu, H. H., and Dai, X. F. (2017). Progress of potato staple food research and industry development in China. J. Integr. Agric. 16, 2924–2932. doi:10.1016/S2095-3119(17)61736-2

CrossRef Full Text | Google Scholar

Zhao, W., Liu, H., Zhang, L., Hu, Z., Liu, J., Hua, W., et al. (2019). Genome-wide identification and characterization of FBA gene family in polyploid crop Brassica napus. Int. J. Mol. Sci. 20, 5749. doi:10.3390/ijms20225749

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Y., Jiao, F., Tang, H., Xu, H., Zhang, L., and Wu, H. (2021). Genome-wide characterization, evolution, and expression profiling of FBA gene family in response to light treatments and abiotic stress in Nicotiana tabacum. Plant Signal Behav. 16, 1938442. doi:10.1080/15592324.2021.1938442

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Q., Tang, D., Huang, W., Yang, Z. M., Zhang, Y., Hamilton, J. P., et al. (2020). Haplotype-resolved genome analyses of a heterozygous diploid potato. Nat. Genet. 52, 1018–1023. doi:10.1038/s41588-020–0699-x

PubMed Abstract | CrossRef Full Text | Google Scholar