- Research article
- Open Access
Identification and biochemical characterization of the fructokinase gene family in Arabidopsis thaliana
BMC Plant Biologyvolume 17, Article number: 83 (2017)
Fructose is an abundant sugar in plants as it is a breakdown product of both major sucrose-cleaving enzymes. To enter metabolism, fructose is phosphorylated by a fructokinase (FRK). Known FRKs are members of a diverse family of carbohydrate/purine kinases known as the phosphofructokinase B (pfkB) family. The complete complement of active fructokinases has not been reported for any plant species.
Protein sequence analysis of the 22 Arabidopsis thaliana pfkB members identified eight highly related predicted proteins, including one with previously demonstrated FRK activity. For one, At1g50390, the predicted open reading frame is half the size of active FRKs, and only incompletely spliced RNAs were identified, which led to a premature stop codon, both indicating that this gene does not produce active FRK. The remaining seven proteins were expressed in E. coli and phosphorylated fructose specifically in vitro leading us to propose a unifying nomenclature (FRK1–7). Substrate inhibition was observed for fructose in all FRKs except FRK1. Fructose binding was on the same order of magnitude for FRK1–6, between 260 and 480 μM. FRK7 was an outlier with a fructose Km of 12 μM. ATP binding was similar for all FRKs and ranged between 52 and 280 μM. YFP-tagged AtFRKs were cytosolic, except plastidic FRK3. T-DNA alleles with non-detectable wild-type RNAs in five of the seven active FRK genes produced no overt phenotype. We extended our sequence comparisons to include putative FRKs encoded in other plant sequenced genomes. We observed that different subgroups expanded subsequent to speciation.
Arabidopsis thaliana as well as all other plant species analyzed contain multiple copies of genes encoding FRK activity. Sequence comparisons among multiple species identified a minimal set of three distinct FRKs present on all species investigated including a plastid-localized form. The selective expansion of specific isozymes results in differences in FRK gene number among species. AtFRKs exhibit substrate inhibition, typical of their mammalian counterparts with the single AtFRK1 lacking this property, suggesting it may have a distinct in vivo role. Results presented here provide a starting point for the engineering of specific FRKs to affect biomass production.
Fructose is an abundant sugar found in plants, generated via the breakdown of sucrose either by invertases, which hydrolytically produce fructose and glucose, or by sucrose synthases, reversibly producing fructose and NDP-glucose using a nucleoside diphosphate. To enter metabolism, fructose must first be phosphorylated most typically to fructose-6-phosphate (F6P) by a fructokinase (FRK, EC 22.214.171.124). Curiously, FRKs are members of a family of diverse kinases known as the phosphofructokinase B (pfkB) family based on their sequence similarity to the founding member, E. coli phosphofructokinase 2 (Pfk-2) which is the minor fructokinase isoform in this species . The pfkB family of enzymes across kingdoms includes ribokinases, phosphofructokinases, adenosine kinases (ADKs) and several others with diverse substrate specificities .
Typical PfkB protein structure is a large domain consisting of a β-sheet sandwiched between several α-helices, and a smaller domain, known as the lid domain, comprised of another β-sheet attached to the larger domain by short loops that act as a hinge [3, 4]. The active site lies in a cleft between the two domains. Substrate binding induces a conformational change in which the lid closes over the substrates . Catalysis occurs while the protein is in the closed state. After catalysis, the protein returns to an open state and the products are released . Binding of substrates to pfkB proteins follows ordered bi-bi kinetics where the carbohydrate enters first, followed by ATP .
PfkB proteins possess two signature motifs: a di-gly (GG) motif in the N-terminal region and a G/AXGD motif in the C-terminal region . The role of each motif has been identified through mutational and structural analyses. The GG motif provides flexibility in the hinge region that connects the lid and the large domain. Substitution of the second glycine to aspartate in the GG motif in the Leishmania donovani pfkB member ADK, decreased enzyme activity to less than 1% of that of the wild-type enzyme and reduced substrate binding affinity, presumably due to the inability of the mutant protein to adopt the closed conformation that supports catalysis . The aspartate in the G/AXGD motif acts as a base during catalysis and activates the C6 fructose hydroxyl group for nucleophilic attack on the γ-phosphate in ATP [2, 4]. Mutation of the codon for aspartate to asparagine in E. coli Pfk-2 significantly abrogated enzymatic activity, but substrate binding was only mildly affected , demonstrating a key catalytic role for aspartate in the G/AXGD motif.
Many pfkB proteins are active as dimers generated through interactions between the lid domains of two monomers [3, 4]. Interactions between β-sheets are mainly hydrophobic and, interestingly, both β-sheets contribute a strand to the other β-sheet in what has been called a β-clasp . ATP has been shown to lead to substrate inhibition of pfkB family members. Excess ATP binds an allosteric site and leads to the formation of tetramers, which are inactive, though it has been shown that substrate inhibition by ATP can still occur in proteins containing substitutions that prohibit formation of tetramers so the regulatory role of ATP is not completely understood [8,9,10]. Aside from their regulation by ATP, pfkB family enzymes are activated by both monovalent and divalent cationic cofactors. Potassium is thought to be the physiological monovalent cofactor. The binding of a potassium ion in the active site activates pfkB enzymes via a conformational shift that results in an anion hole . The ATP-magnesium chelate is considered the actual substrate of pfkB enzymes and magnesium is postulated to aid in catalysis [12, 13].
Despite renewed interest in understanding regulation of carbon flux within the plant and overall plant biomass accumulation for energy purposes, biochemical studies on FRK activities have been characterized only in a few species, mostly tomato, potato and Arabidopsis thaliana. Among these species, FRK activities are best characterized in tomato and four tomato FRK isozymes have been studied functionally and biochemically . As GFP fusions, tomato FRK3 localizes to plastid stroma and the other three FRKs localize to the cytosol . Tomato plants with antisense RNAs targeting FRK1 or FRK2 had markedly less FRK activity than control plants . Surprisingly, fructose levels were lower in all FRK knock-down (KD) plant lines compared to control plants. While the growth of all FRK KD plant lines analyzed in the aforementioned study were stunted, KD of FRK2 affected plant size more than FRK1. Subsequently FRK2 was shown to be involved with both xylem and phloem development . Vascular cells were smaller in FRK2 KD plants than control plants and vessels in the stems were thinner, leading to impaired water conductance. RNAi-mediated KD of tomato FRK3 correlated with diminished stem xylem and reduced water conductance. Simultaneous KD of both FRK2–3 led to severe defects in plant growth that were more drastic than seen in FRK2 KD plants alone . In experiments using the tomato FRK4 promoter to drive expression of diphtheria toxin A in Arabidopsis to assess the tissue-specific expression profile of FRK4, David-Schwartz et al. found that the lethal effects of the toxin were limited to pollen and developing anthers suggesting FRK4 expression is restricted to these organs . The differences in tissue-specific expression and the differential phenotypes seen upon KD of different FRKs suggests non-overlapping roles for different FRKs in tomato.
The growth and vasculature defects in the tomato KD plants suggest that FRKs could be involved in cell wall development, which is a logical possibility for carbohydrate metabolizing enzymes. Sucrose synthase generates fructose and UDP-glucose, the latter of which is a precursor of cellulose, a major constituent of plant cell walls, and fructose is a feedback inhibitor of sucrose synthase. Since FRKs phosphorylate fructose, and thereby reduce intracellular fructose pools, they are hypothesized to indirectly affect cellulose production through modulation of sucrose synthase activity and UDP-glucose production. This hypothesis implicates FRKs as relevant to crop engineering for lignocellulosic biofuel production [20, 21]. In aspen, RNAi targeting FRK2 led to a reduction in the UDP-glucose pool, which was accompanied by reduced cell wall fiber thickness and a lower proportion of cellulose in cell walls . These data together with results from experiments in tomato suggest that further studies of FRKs may lead to a better understanding of, and advances in, the generation of cellulosic biomass.
In Arabidopsis, two electrophoretically distinct FRK activities have been identified and characterized biochemically for substrate specificities and sensitivities to ions and metabolites, though the identities of protein or proteins in each band are unknown . In control experiments while studying proteins related to FRKs, Arsova et al. demonstrated FRK activity for a plastidic Arabidopsis pfkB protein, At1g66430, and denoted it as FRK3 due to sequence similarity to tomato FRK3 . To explore the diversity of FRKs within a single species, we report here the identification of the FRK family in the model plant Arabidopsis thaliana, biochemical characterization of seven active members, and the identification of one pseudogene. All seven members phosphorylated fructose specifically, and henceforth will be referred to as FRKs. Fluorescently tagged Arabidopsis FRK enzymes localized to the cytosol, with the exception of At1g66430/FRK3, which was plastid-localized. Sequence analysis using predicted FRKs from other sequenced genomes support a minimal set of three types of FRKs in plants. Analysis of other plant species also suggests that some of the isozymes may have expanded after speciation. Arabidopsis lines with T-DNA-mediated disruption of single FRKs did not noticeably affect plant growth. The data presented in this manuscript represents the first steps in defining the biological and biochemical roles for the individual FRK enzymes in Arabidopsis.
The Arabidopsis genome encodes seven functional FRKs
Arsova et al.  determined that one member of the pfkB family of proteins in Arabidopsis thaliana ecotype Col-0, encoded in At1g66430, therein named FRK3, was an active FRK after expression of the recombinant protein in E. coli. To determine whether other pfkB enzymes in Arabidopsis are active FRKs, we first generated a phylogenetic tree based on an alignment of the amino acid sequences of the 22 putative pfkB family members in Arabidopsis to determine whether and which other pfkB proteins clade with FRK3 (Fig. 1a). Several clades with 1–2 members and one large sub-family were apparent. One clade contained the two characterized adenosine kinases, ADK1 (At3g09820) and ADK2 (At5g03300) . Approximately half of the Arabidopsis pfkB proteins are un-characterized, so the functional significance of these phylogenetic groupings remains uncertain. The one large sub-family contained sequences from nine proteins divided into two branches. One branch within this clade comprised the two fructokinase-like proteins FLN1 (At3g59480) and FLN2 (At1g69200), for which no kinase activity has been detected . The other branch comprised seven proteins including FRK3. An eighth putative pfkB family member was excluded from these analyses because it was found to be a pseudogene (discussed in the following paragraph). Two peptide sequences that were present in the FRK3 inclusive branch, At1g66430 and At4g10260, were both designated FRK3 in different publications [23, 25]. We have chosen to denote At1g66430 as FRK3 because of its similarity to the tomato FRK3 (Table 1) as previously noted . At4g10260 is designated FRK4 here . FRK1 (At5g51830) and FRK2 (At2g31390) were previously annotated due to their similarity to the tomato isozymes . To visualize specific regions of highest identity we aligned the peptide sequences of the seven putative Arabidopsis FRK proteins (Fig. 1b). FRK3 was predicted to contain a chloroplast transit peptide (cTP) by the chloroP software  and the predicted cTP was removed for the alignment. 47% of the amino acids are identical between all seven members of the putative FRK branch. Both the GG and G/AXGD motifs, underlined in red and green respectively (Fig. 1b), were conserved in all members, and, notably, the G/AXGD motif was in the form GAGD in all members.
Another predicted pfkB protein, At1g50390, claded with the FRK branch, but the predicted protein is only 16 kDa, while other FRKs range from ~34–41 kDa. The 16 kDa protein is missing at least half of the lid domain and alpha helices in the large domain. In addition, there is no evidence that this gene is expressed based on data from TAIR and Genvestigator . According to the TAIR11 gene model, At1g50390 includes three introns (Fig. 2a), but two independent cDNAs synthesized from total Arabidopsis ecotype Col-0 seedling RNA contained sequences identical to the predicted first intron (Fig. 2a-b). As a result, the cDNA contains an in-frame stop codon encoded in the unspliced first intron. If this mRNA were translated, the protein would be 7.85 kDa and lack a functional lid and the catalytic residues. Based on these data we conclude that At1g50390 does not express a functional protein in Col-0 and, thus, have not included it in our phylogenetic or biochemical analyses.
cDNAs for the remaining seven predicted proteins were either isolated from Col-0 cDNA or obtained from the Arabidopsis Biological Resource Center and sequenced to verify the predicted ORFs. All seven ORFs matched the TAIR11 representative gene models, and contained both signature motifs characteristic of pfkB proteins.
All seven putative FRKs are active and phosphorylate fructose specifically
To verify whether the seven members of the putative FRK clade exhibit FRK activity, we expressed each protein in E. coli with 6xHis and FLAG epitope tags and purified them using Ni-agarose resin (Fig. 3). Since the substrate specificity of known pfkB family proteins varies greatly , we measured NADH oxidation as a proxy for FRK activity in the presence of an ATP regeneration system (see Experimental Procedures) in pools of carbohydrate substrates. No change in NADH concentration was seen for any of the proteins in a reaction mixture containing fructose-6-phosphate, sucrose, and N-acetyl-mannosamine suggesting that none of the putative FRKs could phosphorylate any of these carbohydrates (Fig. 4a). Similarly, the proteins were not active in the presence of ribose, xylose, and arabinose (Fig. 4b). In contrast, NADH oxidation was observed in the presence of glucose, fructose, fucose, and tagatose (Fig. 4c), suggesting that one or more of these sugars is a substrate. When these carbohydrates were parsed, fucose, glucose or tagatose alone produced no change in NADH concentration (Fig. 4d-f). However, all seven enzymes were active in the presence of fructose (Fig. 4g), suggesting that each of the seven enzymes possesses FRK activity and that their activity is specific for fructose. The two proteins that comprise the most closely related phylogenetic branch, FLN1 and FLN2, lack detectable fructokinase activity , and the only biochemically characterized enzyme in the unresolved polytomy of proteins related to the FRK clade, At1g17160, specifically phosphorylates ribose , suggesting that the seven FRKs identified herein likely comprise the Arabidopsis FRK family. As such, we propose the nomenclature in Table 1 based on previously defined nomenclature [23, 25] and degree of similarities among the Arabidopsis proteins, and we will use it in the remainder of this manuscript.
The same NADH oxidation assay was used to determine the kinetic parameters of each of the FRK enzymes for their substrates fructose and ATP (summarized in Table 2). The Km values were similar for FRK1–6 as they were all between 470 ± 77 μM and 230 ± 56 μM (Fig. 5a-g, Table 2), whereas FRK7, the least active enzyme exhibited the strongest binding affinity (Km = 12 ± 8 μM, Fig. 5g). Substrate inhibition was observed for all FRKs at levels of fructose above 1 mM with the exception of FRK1 and FRK7 (Fig. 5). The activity of FRK7 was much lower than the other FRKs and, as such, the activity measurements produced more experimental variation, which may have masked its inhibition by fructose.
As demonstrated for other plant FRKs [22, 29], most of the Arabidopsis FRK enzymes were inhibited by higher concentrations of fructose (Fig. 5a-h, Table 2). KmATP values were also very similar and inhibition by ATP was not apparent for any of the FRK enzymes (Fig. 6a-h, Table 2). As such we calculated Vmax and kcat values under varied concentrations of ATP under a constant concentration of fructose (Table 2). Calculated Km and kcat values were also used to calculate the catalytic efficiency of the FRKs (Table 2). Two groups with different amounts of activity were apparent (Fig. 6h). The group with higher activity included FRK1–4 and FRK6; while the group with lower activity included FRK5 and FRK7 (Fig. 6h). The group with higher activity had similar kcat values, which varied between 10.3 and 14.3 s−1, whereas the group with lower activity had kcat values of 4.8 s−1 and 1.5 s−1 for FRK5 and FRK7, respectively (Table 2).
We then assessed the cofactor requirements of a subset of Arabidopsis FRKs, FRK1–3. Unlike other pfkB proteins, most notably the ribokinases and adenosine kinases [30, 31], Arabidopsis FRK1–3 were not significantly sensitive to the absence of inorganic phosphate from the reaction mixture (Fig. 7). Because potassium ions activate pfkB family members via creation of an anion hole in their active site , we tested whether Arabidopsis FRK1–3 were activated to the same extent upon the substitution of sodium for potassium as the monovalent cation in the reaction mixture. Replacement with sodium resulted in the reduced activity of FRK1–3 by roughly 60% of their respective activities in the presence of potassium (Fig. 7). Since magnesium is known to be required for ATP-dependent phosphorylation reactions we next, tested whether or not magnesium is required for FRK activity. Magnesium was omitted from the reaction mixture and EDTA was added to chelate any magnesium that may have co-purified with the enzymes. FRK1–3 activity was completely lost in the reaction mixture containing EDTA (Fig. 7). Other FRK enzymes can use other nucleotide triphosphates besides ATP as kinase substrates, albeit with decreased effectiveness . Accordingly, FRK1 activity decreased by 35% activity in the presence of 1 mM GTP (Fig. 7), while FRK2–3 activities were approximately halved (Fig. 7). Substitution of ATP with 1 mM UTP resulted in the loss of 85% and approximately 75% of FRK1 and FRK2–3 activities, respectively (Fig. 7).
Single T-DNA insertional mutants exhibit no overt growth phenotype
To further our understanding of the in vivo role of FRK enzymes in Arabidopsis, we acquired Arabidopsis plants harboring T-DNA-mediated gene disruption mutants for five of the seven FRKs (T-DNA insertional mutants for FRK4 and FRK6 were unavailable). Each mutant line was confirmed homozygous for the T-DNA insertion in the respective gene (data not shown) and none of the T-DNA insertional mutants made authentic transcripts for their respective disrupted FRK gene (Fig. 8a-e). After 21 days of growth on soil there was no difference in rosette diameter among any of the single FRK mutants or WT control plants (Fig. 8f), suggesting a lack of vegetative phenotype and possible redundancy between some of the FRKs.
Subcellular localization of Arabidopsis FRKs
To gain further insight into the biological roles for Arabidopsis FRK enzymes, we visualized the intracellular localization of FRKs expressed as C-terminal yellow fluorescent protein (YFP) fusions transiently in tobacco leaves (N. benthamiana) and examined by confocal microscopy. The majority of a plant cell’s cytoplasm comprises the vacuole, such that the cytosol is compressed between the vacuolar and plasma membranes. Cytosolic proteins therefore appear around the edges of cells by microscopy, as seen for a known GFP-tagged cytosolic protein, β-glucuronidase (Fig. 9a). FRK1-YFP was observed in the cytosol as the YFP signal was most abundant around the edges of cells with apparent cytosolic bridges connecting parts of the cell (Fig. 9b). The same was true for FRK2-YFP, and FRK4–7-YFP (Fig. 9c, e-h). As mentioned previously, FRK3 was predicted to contain a cTP and, therefore, predicted to be chloroplast-localized. We found FRK3-YFP fluorescence to be coincident with chlorophyll auto-fluorescence, which confirms its predicted chloroplast localization (Fig. 9d).
Other plant species also encode multiple FRKs
In order to determine whether orthologous FRKs are present in other plant species for all seven Arabidopsis FRKs, we compared the amino acid sequences of the seven active Arabidopsis FRKs with those of other plant species with sequenced genomes including Medicago truncatula, Oryza sativa (rice), Solanum lycospersicum (tomato), Populus trichocarpa (poplar), Zea mays (maize), Brachypodium distaychon, and Physcomitrella patens (moss). We then created a phylogenetic tree using all protein sequences along with all Arabidopsis pfkB sequences. We included all proteins that were within one branch from any of the seven active Arabidopsis FRKs. This yielded 51 protein sequences comprised of seven from Arabidopsis, 11 from Medicago, three from rice, four from tomato, four from Populus, four from Brachypodium, eight from maize, and eight from Physcomitrella (moss).
Several interesting relationships were apparent from the resulting phylogenetic tree comprised of the 51 confirmed and putative FRKs (Fig. 10). All Physcomitrella sequences were confined to two clades and, interestingly, those clades included Physcomitrella sequences only, indicating that Physcomitrella FRK sequences are more similar to each other than they are to FRKs from other species. There were two clades that included a single sequence from all other species excluding Physcomitrella (Fig. 10, green shading). One clade contains AtFRK3 and the second includes AtFRK4. Another large clade includes four Arabidopsis proteins, FRK2, FRK5–7, as well as sequences from poplar, tomato, and Medicago, but does not include any sequences from a monocot species. Similarly, another clade contains FRK1, and a tomato sequence, as well as two poplar sequences, and six Medicago sequences that are more distantly related, and no monocot proteins. There is one weakly separated monocot-specific clade with four maize sequences and one each from rice and Brachypodium. A few monocot proteins seem equidistant from all other pfkBs.
Fructose is one of the major sugars found in plant cells as it is a common product of the two major sucrose cleaving enzymes, and it must be phosphorylated to enter metabolism as F6P. Here we provide evidence that there are at least seven active fructokinase isozymes in the model plant Arabidopsis. This trend of multiple FRKs also appears to extend to other plant species as there were at least three putative FRKs encoded by the genome of each species included in our analysis. The Physcomitrella genome encodes eight putative FRKs, which clade together but group away from all other species in our analysis. The fewest FRKs were identified in rice and we speculate that the three putative rice FRK enzymes may represent the minimum FRK cohort that supports plant survival. The groups of FRK paralogs in Arabidopsis, Physcomitrella, and Medicago that are more closely related to each other than those of any of the other species in our analysis suggest that each has a common ancestral gene that then may have duplicated multiple times after speciation.
The similarity between the Arabidopsis FRK protein sequences suggested that their biochemical activity would likely be very similar. Indeed, each of the Arabidopsis FRKs in our analysis phosphorylated fructose specifically. Five Arabidopsis FRKs were roughly equivalent in activity while the other two were significantly lower. The lower activities of FRK5 and FRK7 is surprising given the sequence similarities between these two and the other FRKs. Upon closer inspection of the peptide sequences there are no amino acid or motif differences that are shared between only FRK5 and FRK7, leaving the molecular explanation for the difference in activity as compared to the other FRKs unresolved. The regulation of most Arabidopsis FRKs, with the exception of FRK1 that lacks apparent substrate inhibition by fructose, also appears to be very similar, though the mechanism of substrate inhibition in FRK enzymes by fructose remains unknown. The E. coli PFK-2 shows substrate inhibition by ATP, which involves an allosteric ATP binding site made up partly by the first ATP molecule bound to the protein . The structures of several pfkB family members have been determined and, consistent with the known ordered bi-bi reaction mechanism, the binding site for the carbohydrate ligand is deep in the cleft between the lid and αβα domains and would be occluded if ATP, which binds closer to the surface, bound first . Allosteric substrate inhibition by ATP is a property of other pfkB enzymes such as ribokinases and PFK-2 [8,9,10], though residues that interact with the allosteric ATP in PFK-2  are not conserved in any of the Arabidopsis FRKs studied herein. This finding explains why the Arabidopsis FRKs do not appear to be inhibited by ATP. However, it is possible that an allosteric fructose binds in an analogous way to the allosteric ATP in PFK-2, though currently there is no structural or biochemical data to support this hypothesis.
The ionic cofactor requirements for the tested Arabidopsis FRKs are comparable to those of characterized pfkB family enzymes from other species, which further supports that potassium and magnesium are general cofactors for pfkB family enzymes [11,12,13]. Inorganic phosphate activates ribokinase (RBSK) enzymes [13, 30, 31] and was included in our reactions to determine the kinetic parameters of the Arabidopsis FRKs. However, inorganic phosphate was not required for FRK activity, which suggests that activation by inorganic phosphate may be specific to RBSK enzymes rather than a general characteristic of the pfkB family. FRK activity was diminished but detectable upon the replacement of ATP with either GTP or UTP. This demonstrates that, like previously characterized FRK enzymes, Arabidopsis FRKs are able to utilize other nucleoside triphosphates as phosphate donors , though ATP is likely their preferred substrate.
An important point to consider as it pertains to the possible redundancy between the Arabidopsis FRKs is their tissue-specific expression profiles. In their manuscript describing the tomato FRK4, David-Schwartz et al.  discussed publicly available Arabidopsis microarray data for four of the FRKs characterized in this manuscript. Briefly, FRK1 and FRK3 were expressed at moderate to high levels in all plant tissues with the exception of stamens where their expression was low. FRK7 was expressed at high levels in pollen and in roots, and was also present in siliques and stems. For FRK4, which they report to be the likely functional Arabidopsis homolog of tomato FRK4, expression was confined to mature pollen. Examining the tissue specific expression profiles from the database used in the aforementioned study, Arabidopsis FRK5 and FRK6 are both expressed only at low levels. FRK5 is expressed mainly in roots, flowers, and seeds, while FRK6 is expressed at similar low levels in roots, flowers, and pollen. FRK2, on the other hand, was expressed at moderate levels throughout the plant and expressed highly in hypocotyls, roots, and stems. The ubiquitous and generally moderate to high levels of expression for FRK1, FRK2, and FRK3 further suggests that these enzymes may exhibit functional redundancy. The low level of expression of FRK5 and FRK6 in any tissue suggests that they play minor roles in plant metabolism. We did not examine the effects of the loss of FRK4 as no T-DNA lines were available when we began these studies. However, we speculate that the loss of FRK4 may result in pollen deficiencies due to the tissue specific expression of its putative tomato ortholog .
The biochemical data presented here combined with the similar subcellular localization for six of the seven Arabidopsis FRKs further suggest possible redundancies between two or more Arabidopsis FRKs. Indeed, we observed no discernable differences in the vegetative growth upon T-DNA-mediated inactivation of Arabidopsis FRK1, FRK2, FRK5 or FRK7. FRK3 is the only isozyme that is differentially localized, so we expected to see a phenotypic difference upon its inactivation. However, growth of FRK3 T-DNA line was also not different than wild type.
The role of plastidic FRK is unknown. In prokaryotes it is well established that fructose in the unphosphorylated form will more easily pass through cell membranes. So, unphosphorylated plastidic fructose may also freely cross the chloroplast membrane making it available for phosphorylation by a cytosolic FRK. This scenario explains why the FRK3 mutant may not exhibit a growth phenotype but also does not explain the possible importance of plastidic FRK.
Currently, there is no evidence that FRKs, themselves, are signal transducers (reviewed in [32, 33]). However, there is evidence of fructose-mediated signaling, which could implicate FRKs in plant carbohydrate signaling. The addition of psicose, a C3 epimer of fructose that can be phosphorylated by FRKs, to growth media inhibited root growth in both lettuce and Arabidopsis [34, 35]. The inhibition of root growth by psicose was similar to inhibition by mannose through the known hexokinase (HXK) pathway . Interestingly, psicose-mediated inhibition of root growth was shown to act independently of the HXK pathway . Furthermore, the addition of fructose resulted in the loss of psicose-mediated inhibition . Fructose, itself, acts as a signaling molecule through FRUCTOSE INSENSITIVE 1 (FINS1) and the transcription factor FSQ6/ANAC089, both of which, in turn, are linked to abscisic acid signaling in plants [36, 37]. Whether acting directly through an as yet unknown mechanism or indirectly via their effect on intracellular fructose pools, it is clear that in addition to the assimilation of fructose produced from sucrose, FRKs have some role in plant signaling.
Knockdown experiments in tomato and aspen suggest that FRKs could participate in regulating cell wall development, likely by modulating the flux of cellulose precursors. This implicates FRKs as possible targets for biofuel crop engineering. Mukherjee et al. (19) sought to increase cotton fiber yield via the overexpression of tomato FRK1. They noted an increase in seed cotton production, though much of the increase was likely due to an increase in seed number. The transgenic cotton had reduced sucrose levels, and, curiously, increased stem diameters, though they speculate that other factors could have affected stem diameter . Although an increase in seed number did not fulfill the initial goal, this result suggests that characterizing FRKs in dual-use food and biomass crops such as sorghum could be of particular interest. Our analysis could provide a starting point for engineering of specific FRKs for utilization in such experiments.
This manuscript describes the identification and initial biochemical characterization of seven active fructokinases (FRKs) in the model plant Arabidopsis thaliana. FRKs represent an important group of enzymes that maintain the metabolic flux of fructose, one of the breakdown products of sucrose, the mobile sugar in plants. A large subfamily of pfkB family genes with high peptide sequence similarity was identified that included eight putative enzymes. One was found to likely encode a pseudogene, as its cDNA product contained an unspliced intron that contained a premature stop codon. The remaining seven were bacterially expressed and found to phosphorylate fructose specifically with roughly similar kinetic parameters. All were cytosolic except one, which was plastidic. Inactivation of any of five of the seven for which T-DNA insertion mutants were available resulted in no overt phenotype. Sequence analysis, when extended to other species, suggests that the pattern of encoding several FRKs is a general property of plants. Our study represents the first side-by-side characterization of the FRK family as a whole in a plant species and provides the basis for further studies linking FRK activity through sucrose metabolism to the agronomically important area of cellulosic biomass production.
All chemicals were obtained from Sigma-Aldrich unless otherwise noted.
Plant material and growth conditions
Arabidopsis thaliana Col-0 ecotype seeds (originally from the Arabidopsis Biological Resource Center (ABRC) and propagated in the laboratory) were surface sterilized in a solution of 30% bleach and 0.1% Triton X-100 for 10 min, rinsed with water, stratified for 48 h at 4 °C and grown for 10 days on solid (0.8% bacto-agar, Difco) germination media (GM) (4.3 g/l Murashige and Skoog (MS) basal salts, 2.5 mm MES, 1× B vitamins (0.5 μg/ml nicotinic acid, 1.0 μg/ml thiamine·HCl, 0.5 μg/ml pyroxidine·Cl, 0.1 μg/ml myo-inositol), 1% w/v sucrose (Fisher), pH 5.7). All plate-grown seedlings were grown at 20 °C under constant white light at 40–50 μ mol sec−1 m−2 and used for RNA isolation.
For phenotypic analysis of T-DNA insertional mutants, FRK1, SALK_046463 (frk1–1); two alleles for FRK2, SALKseq_17726 (frk2–1) and SALK_114786 (frk2–2); two alleles for FRK3, SALK_044085 (frk3–1) and SALK_035386 (frk3–2); two alleles for FRK5, SALK_027635 (frk5–1) and SALK_057002 (frk5–2); and one allele for FRK7, GABI_253H07 (frk7–1), seeds were obtained from ABRC and sown directly on soil, thinned to 1 plant per 8 cm by 8 cm square pot, with multiple pots for each genotype distributed among 4–7 trays and grown at 22 °C in cycles of 16 h light and 8 h dark (average 116 μmol⋅sec−1⋅m−2). Untransformed Col-0 was grown as control. The growth experiment was performed twice with identical results. Rosette size at 3 weeks was measured with Image J software from photographs, graphed and analyzed by one-way ANOVA with post-hoc tests in Prism (GraphPad). Genomic DNA isolation and PCR genotypic were performed as described previously . Primers used are listed in Table 4.
Nicotiana benthamiana seeds were obtained from Bo Liu (Department of Plant Biology, from Valerie Williamson, Department of Plant Pathology and Nematology UC-Davis ) and grown in soil at 20 °C 18h day/night for 4 weeks. Fully expanded leaves were used.
Protein sequence analysis
Full-length protein sequences of Arabidopsis FRKs were aligned using MUSCLE  and visualized using Boxshade. Gene identifiers for all proteins analyzed are listed in Table 3. Phylogenetic relationships between the 22 Arabidopsis pfkB-type proteins or sequences of putative FRKs from other species were evaluated by maximum parsimony in PAUP*4.0b10  using an alignment of pfkB protein sequences constructed with MUSCLE  and manually adjusted to remove phylogenetically uninformative sequences. Statistical support for relationships was evaluated by 100 bootstrap replicates for the phylogram in Fig. 1, with ten random addition heuristic searches using the branch and bound algorithm for each bootstrap replicate using PAUP*4.0b10. For the generation of the phylogram in Fig. 10, Phytozome was used to search the genomes of Medicago truncatula, Oryza sativa, Solanum lycopersicum, Populus trichocarpa, Physcomitrella patens, Zea mays, and Brachypodium distachyon for the keyword “pfkB”. All resulting sequences were downloaded and aligned along with the peptide sequences from seven Arabidopsis FRKs in MUSCLE . Then a phylogenetic tree was calculated based on that alignment in PhyML 3.0  and visualized in iTOL . Sequences that were within one branch of Arabidopsis FRKs were considered to be putative FRKs. All putative FRKs were realigned in MUSCLE  and the alignment was manually adjusted using jalview . The phylogenetic tree in Fig. 10 was calculated from the trimmed alignment of putative and known FRKs from all species in PhyML 3.0 using default settings. Statistical support for relationships was evaluated by 100 bootstrap replicates in PhyML 3.0 and visualized in iTOL. Branches with >50% bootstrap support are depicted on the phylogram as colored edges where red represents those close to 50% and green represents those approaching 100%.
For FRK2, FRK3, FRK4, FRK5, and FRK7, RNA was extracted from 7-day-old Arabidopsis seedlings with the RNeasy Plant Mini Kit (Qiagen) and used to generate cDNA using the SuperScript III First-Strand Supermix (Invitrogen). pENTR223 and pDONR221 based clones for FRK1 (G09822) and FRK5 (DQ056446), respectively, were obtained from the ABRC and were recombined directly into destination vectors for bacterial expression. Clones for transient expression in tobacco leaves and for Arabidopsis transformation had their stop codon mutated and Gateway (Invitrogen) recombination sites added via PCR using the primers listed in Table 4. All PCR for cloning was carried out using Phusion high fidelity DNA polymerase (Fisher) according to the manufacturer’s instructions. PCR products for all clones were recombined into pDONR201 or pDONR207 (Invitrogen) via the BP clonase (Invitrogen) reaction and transformed into E. coli strain DH5α (New England Biolabs) for propagation of plasmid DNA. Sequences were verified at the UC Davis DNA sequencing facility. Clones were then recombined into pEAK2  for bacterial expression; or pEARLEYGATE101  or pGWB441  via the LR clonase (Invitrogen) reaction for transient expression in tobacco leaves.
Preparation of recombinant proteins
Proteins were expressed in E. coli strain BL21-pLys-S (New England Biolabs) and purified via Ni Sepharose (GE Healthcare) affinity chromatography. Bacterial cells were lysed by sonication in lysis buffer (50 mM Tris, 100 mM NaCl, 50 mM imidazole, pH 7.5, complete mini protease inhibitor cocktail (Roche)), mixed with resin and the beads washed with at least 10 column volumes of wash buffer (50 mM Tris, 1.5 M NaCl, 50 mM imidazole, pH 7.5) and eluted in elution buffer (50 mM Tris, 100 mM NaCl, 300 mM imidazole, pH 7.5). Eluate was brought to 20% glycerol and was flash frozen and stored at −80 °C. Concentration of recombinant proteins was determined using Protein Assay Reagent (BioRad) and BSA to generate standard curves. Purity was determined by coomassie staining of an 10% SDS-acrylamide gel.
Enzymatic assays were carried out in assays coupling the production of ADP to NADH oxidation and monitored by A340. The reaction mixture consisted of 50 mM Tris pH 7.5, 100 mM KCl, 10 mM MgCl2, 0.1 mM ATP, 1 mM phosphoenolpyruvate, 10 mM inorganic phosphate, 1 mM NADH, 1 mM fructose, 2 U/mL lactate dehydrogenase, 2 U/mL pyruvate kinase. Reactions were carried out in 100 μl volumes in microplate format (NUNC) and A340 was monitored using a Spectra Max 340 microplate reader (Molecular Devices). For the determination of Km values of ribose and ATP, initial reaction rates were measured in concentrations of ribose or ATP, respectively, varied between 2 and 0.0625 mM with all other parameters held constant. Data were plotted and biochemical parameters determined via non-linear curve fitting to either the Michaelis Menten equation or substrate inhibition equation as appropriate in Prism 7 (GraphPad). To determine the effect of monovalent cations, initial reaction rates were determined in the presence of 100 mM KCl, or NaCl and normalized to that of KCl. To determine the effect of divalent cations, initial reaction rates were determined in the presence of 10 mM MgCl2 or EDTA, and normalized to that of MgCl2. Inorganic phosphate was omitted from the reaction to test its requirement. Nucleotide preference was tested by substituting the 0.1 mM ATP with 0.1 mM GTP or UTP, and normalized to that of 0.1 mM ATP. All assays to determine Vmax and Km were performed in triplicate. Assays to determine cofactor requirements were performed in quadruplicate and expressed as scatterplots showing all data points. Statistical differences were determined using two-way analysis of variance and post hoc Dunnett’s multiple comparison test in Prism (GraphPad).
Agrobacteria strain AGL1 (gift from Charles Gasser, Department of Molecular and Cellular Biology, UC-Davis ) carrying vectors to express C-terminally YFP-tagged FRKs were co-infiltrated with AGL Agrobacteria carrying P19 (gift from Richard Michelmore, Plant Sciences Department, UC-Davis ) into Nicotiana benthamiana leaves. After 48 h, a small leaf disc close to the site of infiltration was excised and then prepared and mounted similarly to what was shown in Littlejohn and Love . Briefly, leaf discs were soaked in perfluorodecalin for 5 min and then mounted on slides in a small well created with Carolina Observation Gel (Carolina Biological Supply Company), covered with a coverslip, and imaged on an Olympus FV1000 confocal microscope using a 40× objective and the appropriate filters.
Columbia ecotype Arabidopsis
Chloroplast transit peptide
E. coli phosphofructokinase-2
Yellow fluorescent protein
Morrissey AT, Fraenkel DG. Suppressor of phosphofructokinase mutations of Escherichia coli. J Bacteriol. 1972;112(1):183–7.
Cabrera R, Babul J, Guixe V. Ribokinase family evolution and the role of conserved residues at the active site of the PfkB subfamily representative, Pfk-2 from Escherichia coli. Arch Biochem Biophys. 2010;502(1):23–30.
Chua TK, Seetharaman J, Kasprzak JM, Ng C, Patel BK, Love C, Bujnicki JM, Sivaraman J. Crystal structure of a fructokinase homolog from Halothermothrix orenii. J Struct Biol. 2010;171(3):397–401.
Sigrell JA, Cameron AD, Jones TA, Mowbray SL. Structure of Escherichia coli ribokinase in complex with ribose and dinucleotide determined to 1.8 A resolution: insights into a new family of kinase structures. Structure. 1998;6(2):183–93.
Campos G, Guixe V, Babul J. Kinetic mechanism of phosphofructokinase-2 from Escherichia coli. A mutant enzyme with a different mechanism. J Biol Chem. 1984;259(10):6147–52.
Park J, Gupta RS. Adenosine kinase and ribokinase--the RK family of proteins. Cell Mol Life Sci. 2008;65(18):2875–96.
Datta R, Das I, Sen B, Chakraborty A, Adak S, Mandal C, Datta AK. Mutational analysis of the active-site residues crucial for catalytic activity of adenosine kinase from Leishmania donovani. Biochem J. 2005;387(Pt 3):591–600.
Baez M, Merino F, Astorga G, Babul J. Uncoupling the MgATP-induced inhibition and aggregation of Escherichia coli phosphofructokinase-2 by C-terminal mutations. FEBS Lett. 2008;582(13):1907–12.
Cabrera R, Ambrosio AL, Garratt RC, Guixe V, Babul J. Crystallographic structure of phosphofructokinase-2 from Escherichia coli in complex with two ATP molecules. Implications for substrate inhibition. J Mol Biol. 2008;383(3):588–602.
Cabrera R, Baez M, Pereira HM, Caniuguir A, Garratt RC, Babul J. The crystal complex of phosphofructokinase-2 of Escherichia coli with fructose-6-phosphate: kinetic and structural analysis of the allosteric ATP inhibition. J Biol Chem. 2011;286(7):5774–83.
Andersson CE, Mowbray SL. Activation of ribokinase by monovalent cations. J Mol Biol. 2002;315(3):409–19.
Quiroga-Roger D, Babul J, Guixe V. Role of monovalent and divalent metal cations in human ribokinase catalysis and regulation. Biometals. 2015;28(2):401–13.
Parducci RE, Cabrera R, Baez M, Guixe V. Evidence for a catalytic Mg2+ ion and effect of phosphate on the activity of Escherichia coli phosphofructokinase-2: regulatory properties of a ribokinase family member. Biochemistry. 2006;45(30):9291–9.
Granot D. Role of tomato hexose kinases. Funct Plant Biol. 2007;34(6):564–70.
Damari-Weissler H, Kandel-Kfir M, Gidoni D, Mett A, Belausov E, Granot D. Evidence for intracellular spatial separation of hexokinases and fructokinases in tomato plants. Planta. 2006;224(6):1495–502.
Odanaka S, Bennett AB, Kanayama Y. Distinct physiological roles of fructokinase isozymes revealed by gene-specific suppression of Frk1 and Frk2 expression in tomato. Plant Physiol. 2002;129(3):1119–26.
Damari-Weissler H, Rachamilevitch S, Aloni R, German MA, Cohen S, Zwieniecki MA, Michele Holbrook N, Granot D. LeFRK2 is required for phloem and xylem differentiation and the transport of both sugar and water. Planta. 2009;230(4):795–805.
Stein O, Damari-Weissler H, Secchi F, Rachamilevitch S, German MA, Yeselson Y, Amir R, Schaffer A, Holbrook NM, Aloni R, et al. The tomato plastidic fructokinase SlFRK3 plays a role in xylem development. New Phytol. 2016;209(4):1484–95.
David-Schwartz R, Weintraub L, Vidavski R, Zemach H, Murakhovsky L, Swartzberg D, Granot D. The SlFRK4 promoter is active only during late stages of pollen and anther development. Plant Sci. 2013;199-200:61–70.
Mukherjee T, Ivanova M, Dagda M, Kanayama Y, Granot D, Holaday AS. Constitutively overexpressing a tomato fructokinase gene (LeFRK1) in cotton (Gossypium hirsutum L. cv. Coker 312) positively affects plant vegetative growth, boll number and seed cotton yield. Funct Plant Biol. 2015;42(9):899–908.
Roach M, Gerber L, Sandquist D, Gorzsas A, Hedenstrom M, Kumar M, Steinhauser MC, Feil R, Daniel G, Stitt M, et al. Fructokinase is required for carbon partitioning to cellulose in aspen wood. Plant J. 2012;70(6):967–77.
Gonzali S, Pistelli L, De Bellis L, Alpi A. Characterization of two Arabidopsis thaliana fructokinases. Plant Sci. 2001;160(6):1107–14.
Arsova B, Hoja U, Wimmelbacher M, Greiner E, Ustun S, Melzer M, Petersen K, Lein W, Bornke F. Plastidial thioredoxin z interacts with two fructokinase-like proteins in a thiol-dependent manner: evidence for an essential role in chloroplast development in Arabidopsis and Nicotiana benthamiana. Plant Cell. 2010;22(5):1498–515.
Moffatt BA, Wang L, Allen MS, Stevens YY, Qin W, Snider J, von Schwartzenberg K. Adenosine kinase of Arabidopsis. Kinetic properties and gene expression. Plant Physiol. 2000;124(4):1775–85.
Pego JV, Smeekens SC. Plant fructokinases: a sweet family get-together. Trends Plant Sci. 2000;5(12):531–6.
Emanuelsson O, Nielsen H, von Heijne G. ChloroP, a neural network-based method for predicting chloroplast transit peptides and their cleavage sites. Protein Sci. 1999;8(5):978–84.
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 Bioinforma. 2008;2008:420747.
Riggs JW, Rockwell NC, Cavales PC, Callis J. Identification of the Plant Ribokinase and Discovery of a Role for Arabidopsis Ribokinase in Nucleoside Metabolism. J Biol Chem. 2016;291(43):22572–82.
Martinez-Barajas E, Krohn BM, Stark DM, Randall DD. Purification and characterization of recombinant tomato fruit (Lycopersicon esculentum Mill.) fructokinase expressed in Escherichia coli. Protein Expr Purif. 1997;11(1):41–6.
Maj MC, Singh B, Gupta RS. Pentavalent ions dependency is a conserved property of adenosine kinase from diverse sources: identification of a novel motif implicated in phosphate and magnesium ion binding and substrate inhibition. Biochemistry. 2002;41(12):4059–69.
Park J, van Koeverden P, Singh B, Gupta RS. Identification and characterization of human ribokinase and comparison of its properties with E. coli ribokinase and human adenosine kinase. FEBS Lett. 2007;581(17):3211–6.
Granot D, Kelly G, Stein O, David-Schwartz R. Substantial roles of hexokinase and fructokinase in the effects of sugars on plant physiology and development. J Exp Bot. 2014;65(3):809–19.
Li L, Sheen J. Dynamic and diverse sugar signaling. Curr Opin Plant Biol. 2016;33:116–25.
Kato-Noguchi H, Okada K. Phosphorylation and growth inhibitory activity of all ketohexose analogs. Weed Biol Manag. 2013;13(3):114–8.
Kato-Noguchi H, Takaoka T, Izumori K. Psicose inhibits lettuce root growth via a hexokinase-independent pathway. Physiol Plant. 2005;125(3):293–8.
Cho YH, Yoo SD. Signaling Role of Fructose Mediated by FINS1/FBP in Arabidopsis thaliana. PLoS Genet. 2011;7(1):e1001263.
Li P, Wind JJ, Shi X, Zhang H, Hanson J, Smeekens SC, Teng S. Fructose sensitivity is suppressed in Arabidopsis by the transcription factor ANAC089 lacking the membrane-bound domain. Proc Natl Acad Sci U S A. 2011;108(8):3436–41.
Gilkerson J, Perez-Ruiz JM, Chory J, Callis J. The plastid-localized pfkB-type carbohydrate kinases FRUCTOKINASE-LIKE 1 and 2 are essential for growth and development of Arabidopsis thaliana. BMC Plant Biol. 2012;12:102.
Goodin MM, Zaitlin D, Naidu RA, Lommel SA. Nicotiana benthamiana: Its History and Future as a Model for Plant-Pathogen Interactions. Mol Plant-Microbe Interact. 2015;2015(1):28–39.
Edgar RC. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004;5:113.
Rogers JS, Swofford DL. Multiple local maxima for likelihoods of phylogenetic trees: a simulation study. Mol Biol Evol. 1999;16(8):1079–85.
Dereeper A, Guignon V, Blanc G, Audic S, Buffet S, Chevenet F, Dufayard JF, Guindon S, Lefort V, Lescot M, et al. Phylogeny.fr: robust phylogenetic analysis for the non-specialist. Nucleic Acids Res. 2008;36(Web Server issue):W465–9.
Letunic I, Bork P. Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016;44(W1):W242–5.
Waterhouse AM, Procter JB, Martin DM, Clamp M, Barton GJ. Jalview Version 2--a multiple sequence alignment editor and analysis workbench. Bioinformatics. 2009;25(9):1189–91.
Kraft E. An investigation of the ubiquitin conjugating enzymes and RING E3 ligases in Arabidopsis thaliana. PhD thesis. Davis: University of California; 2007.
Earley KW, Haag JR, Pontes O, Opper K, Juehne T, Song K, Pikaard CS. Gateway-compatible vectors for plant functional genomics and proteomics. Plant J. 2006;45(4):616–29.
Nakagawa T, Suzuki T, Murata S, Nakamura S, Hino T, Maeo K, Tabata R, Kawai T, Tanaka K, Niwa Y, et al. Improved Gateway binary vectors: high-performance vectors for creation of fusion constructs in transgenic analysis of plants. Biosci Biotechnol Biochem. 2007;71(8):2095–100.
Lazo GR, Stein PA, Ludwig RA. A DNA transformation-competent Arabidopsis genomic library in Agrobacterium. Biotechnology (N Y). 1991;9(10):963–7.
Guerra D, Chapiro SM, Pratelli R, Yu S, Jia W, Leary J, Pilot G, Callis J. Control of Amino Acid Homeostasis by a Ubiquitin Ligase-Coactivator Protein Complex. J Biol Chem. 2017;292(9):3827–40.
Littlejohn GR, Love J. A simple method for imaging Arabidopsis leaves using perfluorodecalin as an infiltrative imaging medium. J Vis Exp. 2012;59:e3394.
We thank Drs. Kentaro Inoue and Li Tian (UC-Davis) for helpful discussions and Li Tian for her helpful comments on the manuscript. We thank Rachel Carr for help with the enzymology and Frank Lin for assistance with analysis of insertional mutants. We thank Dr. Damian Guerra and Katrina Linden for their helpful comments on the manuscript and Charles Gasser for assistance with phylogenetic analyses. Microscopy was performed at the UCD Light Microscopy Imaging Facility and plants were grown in the Controlled Environment Facility, both at UCD.
This work was funded by the Division of Chemical Sciences, Geosciences, and Biosciences, Office of Basic Energy Sciences of the U.S. Department of Energy through Grants DE-FG02-09ER16077 and DE- SC0002175 to JC and the Eric and Louise Conn Graduate Student Award in Plant Biochemistry to JR.
Availability of data and materials
Phylogenetic data generated and analyzed during the current study are available in the TreeBASE repository, [submission #20814 and #20818 for Figs. 1 and 10, respectively]. URLs to access Figs. 1 and 10 in TreeBASE are http://purl.org/phylo/treebase/phylows/study/TB2:S20814 and http://purl.org/phylo/treebase/phylows/study/TB2:S20818, respectively.
Any other datasets used and/or analyzed during the current study available from the corresponding author on reasonable request.
Both JR and JC conceived the project, analyzed data, wrote and edited the manuscript. JR, PC, JC and SC performed experiments. JC secured funding. This work was performed in partial fulfillment of the requirements for the degree of PhD by JR All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.