Comparative transcriptome analysis coupled to X-ray CT reveals sucrose supply and growth velocity as major determinants of potato tuber starch biosynthesis

Abstract

Background

Even though the process of potato tuber starch biosynthesis is well understood, mechanisms regulating biosynthesis are still unclear. Transcriptome analysis provides valuable information as to how genes are regulated. Therefore, this work aimed at investigating transcriptional regulation of starch biosynthetic genes in leaves and tubers of potato plants under various conditions. More specifically we looked at gene expression diurnally in leaves and tubers, during tuber induction and in tubers growing at different velocities. To determine velocity of potato tuber growth a new method based on X-ray Computed Tomography (X-ray CT) was established.

Results

Comparative transcriptome analysis between leaves and tubers revealed striking similarities with the same genes being differentially expressed in both tissues. In tubers, oscillation of granule bound starch synthase (GBSS) expression) was observed which could be linked to sucrose supply from source leaves. X-ray CT was used to determine time-dependent changes in tuber volume and the growth velocity was calculated. Although there is not a linear correlation between growth velocity and expression of starch biosynthetic genes, there are significant differences between growing and non-growing tubers. Co-expression analysis was used to identify transcription factors positively correlating with starch biosynthetic genes possibly regulating starch biosynthesis.

Conclusion

Most starch biosynthetic enzymes are encoded by gene families. Co-expression analysis revealed that the same members of these gene families are co-regulated in leaves and tubers. This suggests that regulation of transitory and storage starch biosynthesis in leaves and tubers, respectively, is surprisingly similar. X-ray CT can be used to monitor growth and development of belowground organs and allows to link tuber growth to changes in gene expression. Comparative transcriptome analysis provides a useful tool to identify transcription factors possibly involved in the regulation of starch biosynthesis.

Background

Starch is not only the most important carbohydrate source to the human diet, but has major industrial applications. It consists of two major fractions, amylose and amylopectin. Amylose is essentially a linear polymer of glucose units linked with alpha (1,4) bonds whilst amylopectin has a higher percentage of branched alpha (1,6) bonds. Due to its importance, crop plants producing starch in large quantities have been extensively researched. The fourth most important crop in the world in terms of total biomass produced is potato and this is due to its starchy tuber which can store up to 80% of its dry weight as starch [1].

In potato starch is either accumulated transiently in leaves or as storage starch in tubers. In leaves starch is synthesised in the chloroplast from triose-phosphates produced during photosynthesis. After several intermediate steps, glucose 6-phosphate is converted to glucose 1-phosphate by the enzyme phosphoglucomutase (PGM). Glucose 1-phosphate, along with ATP serves as substrate for ADP-glucose production by ADP-glucose pyrophosphorylase (AGPase), which is the first reaction committed to starch biosynthesis. ADP-glucose is then the glycosyl donor for the various starch synthases forming linear glucans which are subsequently branched by branching enzymes to produce starch.

For starch production in the tuber photoassimilates, in the form of sucrose, must first be imported from photoautotrophic tissue via the phloem [2]. There is conflicting evidence to whether carbon accumulation is sink or source limited. It has been shown that under normal conditions flux control of tuber starch biosynthesis is mostly source regulated [3]. In another study, reduced photosynthetic activity of potato leaves by silencing cytosolic fructose 1, 6 bisphosphatase did not have an effect on tuber yield or plant growth and it was concluded that tuber starch biosynthesis is not source limited [4]. Both studies however emphasize the importance of sucrose supply to the tuber. This idea is further supported by constitutive [5] or phloem-specific [6] antisense inhibition of the sucrose transporter. Sucrose transport activity is essential for apoplastic phloem loading, hence, silencing its expression lead to a reduced phloem loading of sucrose which was accompanied by decreased photosynthetic rates and reduced tuber yield of transgenic potato plants. In another study [7], over-expression of a sucrose transporter from spinach in potato lead to a reduced sucrose level in leaves and an increased sucrose content of tubers. However, this had no effect on tuber starch content. Even though this did not lead to an increased starch content of tubers, it did provide evidence as to the important role of sucrose as regulator of carbon metabolism [7]. Of particular note was the decrease in plastidial amino acid synthesis even though sucrose was still in abundant supply.

There is also conflicting evidence as to whether the supply rate of sucrose from the leaves to tubers is constant throughout the diurnal cycle or whether there is a diurnal rhythm. It has been shown that source to sink carbon flux is constant [8]. It can also be argued that for starch biosynthesis in tubers to be unaffected by limitations in photosynthetic capacity [4], sucrose supply to the tuber should be constant. It has been shown that there are significant differences of tuber sucrose content at the start and end of the light period and that this has an effect on tuber metabolism [9].

After the formation of tubers they represent the predominant sink in the plant. The onset of tuberisation is marked by many physiological and biochemical changes in the stolons. The major changes which are seen as markers of tuberisation are the switch from apoplastic to symplastic sucrose unloading and the change from hydrolytic to sucrolytic cleavage of sucrose [10]. Sucrose synthase (Susy) activity increases with the development of the tuber and is the major determinant of potato sink strength [11]. After the cleavage of sucrose to UDP-glucose and fructose by Susy, the two molecules must be converted to glucose 6-phosphate via different pathways for starch biosynthesis. The major difference between starch biosynthesis in tubers compared to leaves is that hexose-phosphate, in the form of glucose 6-phosphate, and ATP must be imported from the cytosol. Although there is evidence that glucose 1-phosphate can also be transported across the plastid membrane and directly used for starch biosynthesis via starch phosporylase in potato tuber discs [12], it seems that in vivo glucose 6-phosphate is the predominant form of hexose-phosphate transport across the amyloplast membrane. Transport of glucose 6-phosphate and ATP across the amyloplast membrane is facilitated by the glucose 6-phosphate translocator (GPT) [13], and the plastidial ATP/ADP translocator (NTT) [14], respectively. Simultaneous over-expression of these transporters lead to an increase in total starch yield per plant showing that starch biosynthesis is co-limited by hexose-phosphate and ATP in the plastid [15]. Glucose 6-phosphate is then converted to glucose 1-phosphate by plastidial PGM. Silencing of plastidial PGM lead to dramatic reduction in starch content of the tubers, illustrating not only the importance of this enzyme in starch biosynthesis, but provides compelling evidence that starch is synthesised from imported glucose 6-phosphate [16]. As mentioned earlier, glucose 1-phosphate and ATP serve as substrates for ADP-glucose synthesis by AGPase and further starch biosynthesis is similar to leaves. Since large scale comparative transcriptome analysis comparing leaf and tuber starch biosynthesis have not been done, it is not known whether or not the same isoforms of these genes are active in both tissues. The importance of AGPase in starch biosynthesis has been shown on various occasions and reduced activity lead to low starch, high sucrose containing tubers [17]. As far as increasing starch biosynthesis by over-expression of AGPase is concerned, two studies over-expressing an E. coli AGPase give conflicting results. Stark et al. [18] could increase starch content by over-expression but Sweetlove et al. [19] showed that increased biosynthesis is accompanied by increased breakdown. There are mainly four soluble starch synthase enzymes responsible for biosynthesis, starch synthase I, II, III, IV and one that is exclusively granule bound, named granule bound starch synthase (GBSS). Soluble starch synthases and the branching enzymes are believed to be responsible for amylopectin synthesis. The major starch synthases are starch synthase II and III and silencing of these two isoforms lead to greatly altered amylopectin structure [20]. Granule bound starch synthase is exclusively responsible for the synthesis of amylose and a potato mutant lacking any activity produce amylose free starch [21]. Starch synthase IV seems to be important for starch granule initiation and mutants lacking this isoform do not only show changes in diurnal starch accumulation, but also have fewer, and bigger, starch granules [22, 23]. Despite the many isoforms of branching enzyme, they can be grouped in starch branching enzyme A (SBEA) and B (SBEB). The majority of branching enzyme activity can be attributed to SBEB [24], but silencing of only SBEB did not lead to a change in starch structure. Simultaneous silencing of SBEA and SBEB, however, lead to production of high amylose potato starch [25]. Table 1 contains all the genes involved in starch metabolism which will be further discussed in this paper. The pathways in which these enzymes are involved are shown in figure 1.

Table 1 List of genes involved in starch metabolism discussed in this paper.

Figure 1

Proposed pathway of starch metabolism in photosynthetic and non-photosynthetic tissue.

One of the major problems with potato tuber research is the fact that tuber growth and induction rates are not synchronised. Although tuber developmental stages are well defined [26], this does not mean that tubers in the same stage have similar growth- or biochemical characteristics. Adding to this, tubers are underground organs making in vivo analysis without damaging the plant almost impossible. To date no study has determined the in vivo growth velocity of individual tubers. X-ray computed tomography (X-ray CT) provides the opportunity to determine the velocity of tubers without physically damaging the plant. It has been used to study underground plant organs [27], but to date no study has determined tuber growth velocity using this method.

Even though the pathway of starch biosynthesis is well understood, mechanisms regulating biosynthesis are still unclear. Smith et al. [28] did an extensive study on transcriptional regulation of starch metabolism in Arabidopsis leaves over a 24 hour diurnal period. This study was very useful elucidating the regulating machinery of starch metabolism. The design of a custom microarray has made it possible to do similar types of analysis in potato. The POCI microarray was designed from the largest collection of ESTs from potato yet [29] and can be used to study all aspects of the potato transcriptome, which include starch biosynthesis. Due to the abundance of transcriptome data available, co-expression analysis has become a common method to identify new genes involved in specific metabolic pathways. The principle of the technique is that genes involved in similar processes would be expressed or inhibited at the same time point or under similar conditions. This would not only include structural genes, but possibly also regulatory genes like transcription factors. They are able to bind to specific sequences of various targets and thus have the ability to regulate entire metabolic pathways [30]. Moreover, over-expression of two transcription factors in tomato led to the over-expression of the entire anthocyanin synthesis pathway, producing purple anthocyanin rich tomatoes [31]. These results show that by manipulating transcription factors entire biosynthetic pathways can be influenced.

The aim of this work was to investigate the transcriptional regulation of starch biosynthesis in potato under various conditions. More specifically we looked at similarities between gene expression in leaves and tubers. Furthermore we established a new technique using X-ray computed tomography to determine tuber growth velocity in vivo and used it to analyse gene expression in tubers growing at different velocities. Finally comparative analysis of transcription profiles were used to identify transcription factors possibly regulating starch biosynthesis.

Results and discussion

Starch biosynthesis in potato leaves follow carbohydrate accumulation and show similarities to tuber starch biosynthesis

Leaf samples for starch and RNA extractions were taken from potato plants grown under a 14 hour light and 10 hour dark cycle. The light conditions were chosen to ensure diurnal turnover of transitory starch. As was expected, starch levels were the lowest at the start of the light period and accumulated during the day to a high at the end of the light period (Figure 2A). The same pattern was observed for sucrose although the accumulation and decline was more rapid (Figure 2B). The levels of starch and sucrose compared well to what was previously measured in potato leaves [4].

Figure 2

Diurnal starch and sucrose content of leaves over a sixteen hour period. A) Starch and B) sucrose content. Error bars indicate standard deviation (n = 3).

After confirming that there was a turnover of transient starch under the specific light conditions used, gene expression at different times of the day was analysed by microarray. To investigate how much similarity there is between leaf and tuber starch biosynthesis, the expression of genes known to be important for tuber starch biosynthesis was investigated in leaves. Starch biosynthetic genes had a strong diurnal rhythm of expression. Interestingly, genes involved in the import of glucose 6-phosphate and ATP into the plastid are co-regulated with ADP-glucose pyrophosphorylase, suggesting that both genes are also involved in transitory starch biosynthesis in leaves (Figure 3A). To rule out that the observed expression of GPT and NTT in leaf extracts was due to contaminating epidermis rather than mesophyll cells, cell-specific RNA analysis was performed. To this end epidermis stripes were harvested and expression of GPT and NTT were probed by quantitative real-time PCR (qRT-PCR). As shown in figure 4, GPT and NTT expression in epidermis cells was about ten fold lower compared to whole leaf extracts. This supports the hypothesis that GPT and NTT are indeed required for photosynthetic starch biosynthesis in leaves. However, other cell types such as companion and bundle sheath cells could still contribute to the observed expression pattern. It can also not be excluded that the transporters are exporting substrates from the plastid.

Figure 3

Diurnal expression of genes known to be involved in starch biosynthesis. A) Import of glucose 6-phosphate and ATP into the plastid and the conversion thereof to ADP-Glucose. GPT1 (light blue), GPT2 (dark blue), PGM (orange), NTT1 (brown), AGPase LS (light green), AGPase SS (dark green). B) Starch synthases and branching enzymes. SSII (dark grey), SSIII (light grey), GBSS (blue), SSIV (black), SBEA (light purple), SBEB (dark purple), ISA1 (dark red) and ISA2 (pink). C) Photosynthetic and Calvin cycle related genes. Plastocyanin (pink solid), CAB ( pink dotted), Rubisco (blue dotted), TPT (blue solid). D) Sucrose cleavage and phosphorylation. CW-Inv (olive green), FK (blue), Susy (gold), HK1 (dark yellow) and HK2 (green) Values are the mean of two replicates.

Figure 4

Relative expression of GPT and NTT in epidermal and whole leaf tissue. A) GPT and B) NTT. Error bars represent standard deviation (n = 3).

Most starch synthases and branching enzymes had a similar pattern to the above mentioned genes in figure 3A with the exception of GBSS which was highest expressed two hours into the light period (Figure 3B). Smith et al. [28] argue and provide evidence that since the enzyme is present within the granule, the protein is degraded together with starch at night and must very quickly be re-synthesised in the morning. Two genes involved in the light reaction of photosynthesis, plastocyanin and chlorophyll a/b binding protein, were already up-regulated at the first time point which was taken moments after the lights came on and were highest expressed two hours into the light. Ribulose 1,5 bisphosphate carboxylase oxygenase (Rubisco) had an expression pattern similar to that of starch biosynthetic genes (Figure 3C). The sucrose cleavage enzymes cell wall-bound invertase and sucrose synthase had very different expression patterns with Susy 4 being much stronger regulated. Susy 4 increased during the light and went down in the dark. Hexokinase 1 was not diurnally regulated, whilst hexokinase 2 and fructokinase had a similar pattern of increasing early in the morning and declining at the end of the light period (Figure 3D).

For analysis of gene expression during tuber induction, tuber developmental stages were designated according to Kloosterman et al. [26] and the expression of genes were studied by microarray analysis. Stage 1, representing an unswollen stolon, showed very low expression for all starch biosynthetic genes except starch synthase IV which was highest expressed at stage 1, fitting with its proposed role in starch granule initiation [22, 23]. From stages 3-5 there was an increase in the expression of all starch biosynthetic genes with the exception of starch synthase IV (Figure 5A and 5B). As evidence that stage 1 was before the onset of tuberisation, cell wall-bound invertase expression was very high and sucrose synthase 4, a tuber expressed isoform, low. From stages 3-5 the expression cell wall-bound invertase went down and that of Susy 4 increased (Figure 5C). The expression patterns compared well to earlier experiments [29]

Figure 5

Expression of genes known to be involved in starch biosynthesis during tuber induction. A) Import of glucose 6-phosphate and ATP into the plastid and the conversion thereof to ADP-Glucose. GPT1 (light blue), GPT2 (dark blue), PGM (orange), NTT1 (brown), AGPase LS (light green), AGPase SS (dark green). B) Starch synthases and branching enzymes. SSII (dark grey), SSIII (light grey), GBSS (blue), SSIV (black), SBEA (light purple), SBEB (dark purple), ISA1 (dark red) and ISA2 (pink). C) Sucrose cleavage and phosphorylation. Cw-Inv (Black), FK (blue), Susy (gold), HK1 (Dark yellow) and HK2 (green). Values are the mean of two replicates.

qRT-PCR is a well accepted method for verifying microarray data and was used to validate the expression patterns of a selected number of genes in independent samples. The expression patterns and levels of Susy 4, GPT and GBSS compared well to the microarray data in both leaves and tubers and confirmed the accuracy of the microarray results (Figure 6).

Figure 6

Quantitative real time PCR confirmation of microarray results. A-C) Relative expression of GPT2, Susy4 and GBSS diurnally in leaves. D-F) Relative expression of GPT2, Susy4 and GBSS in tubers. Error bars represent standard deviation (n = 3).

In contrast to starch biosynthesis, the pathway of starch degradation in potato tubers is not well understood. Currently it is believed that the major pathway of starch degradation in Arabidopsis and other organisms proceeds via ί-amylase [32]. Starch must first be phosphorylated by glucan-water dikinase (GWD), and to a lesser extent phospho-glucan water dikinase [33, 34]. It is then either directly degraded to maltose by ί-amylase or first to linear glucans by Isoamylase 3. Before complete hydrolysis of phosphorylated glucans can occur, phosphate groups must be removed by specific a phosphatase, known as SEX4 phosphoglucan phosphatase [35]. The linear glucans are then hydrolysed to maltose by ί-amylase. Maltose is subsequently exported to the cytosol via the maltose transporter (MEX1) [36] and it is believed that the disproportioning enzyme 2 (DPE2) plays an important role in its further degradation [37]. Maltose can also be degraded by disproportioning enzyme 1 (DPE1) inside the plastid. For potato evidence has been provided that DPE2 might be present inside the plastid [38], which questions the role of MEX1. If the expression data is compared to that of Arabidopsis, indications are there that a similar pathway for starch degradation probably occurs in potato leaves, except that isoamylase 3 was not differentially expressed. The other genes mentioned were all differentially co-expressed in leaves, and expression patterns followed that of starch accumulation in the leaf (Figure 7). The entire microarray dataset for both the diurnal leaf experiment and the tuber induction experiment have been deposited on ArrayExpress (accession numbers E-MEXP-2481 Ferreira et al. Potato diurnal leaf time-course and E-MEXP-2482 Ferreira et al. Potato tuber induction). POCI sequence and annotation data are available through the POCI online tool http://pgrc.ipk-gatersleben.de/poci.

Figure 7

Diurnal expression of genes known to be involved in starch degradation. GWD (black), DPE1 (green), DPE2 (light blue), MEX1 (red), Alpha amylase (grey), Isoamylase 3 (dark blue), PCT-BMY (brown) and SEX4 (orange). Values are the mean of two replicates.

Diurnal oscillation of GBSS in potato tubers can be linked to differences in sucrose supply

It is known that starch biosynthetic genes are diurnally regulated by several factors, with sucrose and the circadian clock seemingly being the most important. Bläsing et al. [39] showed that between 30-50% of genes in Arabidopsis rosettes show diurnal changes in their transcripts and that this was especially true for genes involved in redox regulation, nutrient acquisition and assimilation, and starch and sucrose metabolism. Comparative analysis of nutrient feeding and diurnal transcription profiles indicate that sugars make a major contribution to diurnal regulation. Furthermore, Osuna et al. [40] analysed gene expression in carbon-deprived Arabidopsis seedlings after the addition of sucrose. Genes involved in central carbon metabolism, and more specifically starch biosynthesis, showed a response to sucrose and this did lead to an increase in starch content. A second major regulator of diurnal gene expression seems to be the circadian clock. Bläsing et al. [39] identified a subset of 373 genes known to be circadian regulated [41]. The gene set was used in a principle component analysis which showed that sucrose and the circadian clock are the predominant factors in regulating diurnal gene expression and that light, nitrogen and water deficiency makes a smaller contribution. Starch degradation related genes also have a strong diurnal rhythm in Arabidopsis which was maintained under continuous light, but not continuous darkness [42].

Since GBSS is one of the strongest diurnally regulated genes involved in starch biosynthesis [28], it was investigated whether this gene is also diurnally regulated in tubers. GBSS expression was measured at different times of the day and under different light regimes. Plants were grown in normal long day conditions and samples were taken from stolon-ends and tubers at different time-points of the diurnal cycle. GBSS expression in tubers oscillated during the day and was highest at the end of the light period and lowest 2 hours into the next light period (Figure 8A). Expression was also significantly lower in tubers from plants kept in twenty four hours of darkness compared to tubers from plants kept in a normal light/dark cycle (Figure 8B).

Figure 8

GBSS relative expression and stolon sucrose content at different time-points of the day. A) Diurnal expression of GBSS in tubers. B) GBSS expression at the 24 hour time-point from plant grown in light/dark cycle and from plants kept in twenty four hours of darkness. C) Stolon sucrose content at different time-points of the day (dark grey bars) and at the same time-points from plant kept in darkness from 0 hours onward (dark grey bars). Two values were not determined (n.d.). Error bars represent standard deviation (n = 3-7).

As mentioned, sucrose and the circadian clock seem to be important regulators of diurnal gene expression and this is especially true for GBSS. Tenioro et al. [43] showed that GBSS is strongly regulated by the circadian clock and that expression is markedly lower in mutants lacking clock genes LHY and CCA-OX respectively. In a detailed study done on GBSS in snapdragon (Antirrhinum majus) it was shown that GBSS is diurnally regulated in leaves even under continuous light and it was concluded that the regulation is due to the circadian clock. This was not the case in snapdragon roots though, where expression was the same in the middle of the day and in the middle of the night [44]. Also in rice leaves GBSS continues its diurnal cycling under continuous light suggesting circadian regulation, but expression can be induced by nitrogen starvation or sucrose feeding and repressed by darkness, indicating the importance of sucrose in its regulation [45]. Moreover, sucrose floating experiments with potato leaves show that GBSS expression can be induced by sucrose [46].

To investigate whether oscillation of GBSS might be due to changes in sucrose supply from the source, sucrose content of stolon-ends attached to a tuber was measured. Phloem sucrose content is very high compared to those tissues surrounding it [47] and phloem signifies a large proportion of the total stolon tissue [48]. This makes it possible to determine sucrose import to the tuber by measuring sucrose content of the stolon-end [3]. Sucrose content differed significantly during the day and was highest at the end of the light and lowest at the end of the dark period. When plants were kept in constant darkness, sucrose content of stolon-ends declined linearly over time indicating that changes in sucrose supply from the leaves contributes to the oscillating expression of GBSS in tubers (Figure 8C).

Although GBSS expression in tubers can be linked to diurnal changes in sucrose supply from the source, caution should be exercised in the interpretation of the significance of this in terms of enzyme activity. As mentioned earlier, the level of GBSS protein does change substantially during the day in Arabidopsis leaves and the reason for this is probably the location of GBSS within the granule [28]. This is also true for algae Chlamydomonas reinhardtii suspension cultures where GBSS expression and enzyme activity correlates with starch levels. The authors state that the correlation is probably due to the fact that the analysis was conducted in suspension cultures, where new cells are produced constantly leading to continuous production of new GBSS protein [49]. This however is not necessarily true for GBSS in other tissues or for other enzymes of starch biosynthesis and several studies have shown that diurnal changes in expression do not lead to changes in protein levels [9, 42]. However, it still remains an interesting finding that GBSS expression in tubers follows a diurnal rhythm which declines when sucrose supply from the leaves are reduced.

Starch biosynthetic gene expression is influenced by tuber growth velocity

As mentioned earlier, tuber initiation and growth rates are not synchronised. Because analysis of tubers growing at different velocities would provide important information on factors determining growth, it was thought necessary to develop a method to determine growth velocity before harvesting the tubers. For this, X-ray CT was used to determine tuber volume over a time course and then calculate the growth velocity. An overview of the method for calculating tuber volume using X-ray CT is given in figure 9. The accuracy of the growth velocity determination depends largely on the accuracy of volume calculation and it was crucial to establish whether tuber volume could be accurately calculated. Figure 10A illustrate that potato tubers were clearly distinguishable from the surrounding soil, and a histogram demonstrated that segmentation of tuber material was possible. The software used for segmentation also allowed for manual correction of possible mistakes made during the automatic segmentation process. Tuber growth velocities were calculated for seventeen tubers. To validate the method, calculated tuber volumes were compared to measured tuber volumes after harvest. X-ray CT calculated volumes showed a correlation coefficient with experimental volume measurements of 0.986 showing that X-ray CT calculated volumes were accurate (Figure 10B).

Figure 9

Schematic scheme illustrating how tuber volume is measured using X-ray CT. Potato plants are scanned with X-ray in a chamber containing an X-ray beam and a two dimensional detector. After projecting X-ray images of the potato plant on the detector at different angles, the projections are reconstructed in silico to create a three dimensional image. From this image the tuber volume can be calculated.

Figure 10

Potato tuber segmentation and volume calculation. A) Two dimensional X-ray images illustrating that potato tubers can be distinguished from the surrounding soil. The histogram shows that segmentation is possible. The red bar indicates the grey level threshold selected. B) Linear regression of X-ray calculated and real volume measurements confirming the accuracy of X-ray CT calculated volumes.

Six tubers with different growth velocities (Figure 11) were selected for microarray hybridisation and the gene expression patterns were analysed. Interestingly no linear correlation between the expression of starch biosynthetic genes and growth velocity could be observed. The expression of genes involved in starch biosynthesis is shown in figure 12. The expression of these genes, and also other genes involved in starch biosynthesis, did not show large differences in expression levels between tubers that were still growing, albeit at very different velocities. Starch biosynthetic gene expression was however much lower in tubers that have virtually stopped growing and it seems as though the relationship between growth velocity and starch biosynthetic gene expression is qualitative rather than quantitative. This was especially true for Susy, a major determinant of tuber sink strength [11].

Figure 11

Estimated growth velocity of tubers in cubic centimetre volume increase per day. Tubers chosen for microarray hybridisation are marked by arrow heads

Figure 12

Relative expression of starch biosynthetic genes in tubers growing at different velocities. GPT1 (pink), GPT2 (yellow), ISA2 (orange), AGPase LS (purple), AGPase SS (brown) and Susy4 (green). Values are the mean of two replicates.

These data indicate that tubers that look visually similar have large differences in gene expression depending on their growth stage. Moreover, transcription profiles reveal many genes that are differentially expressed between growing and non-growing tubers which provide important information towards the identification of factors determining tuber growth. The entire microarray dataset for the tuber growth velocity experiment has been deposited on ArrayExpress (accession number E-MEXP-2484 Ferreira et al. tuber growth velocity).

Comparative analysis of transcription profiles reveals genes possibly regulating starch biosynthesis

The similarities observed in terms of starch biosynthetic gene expression between leaves and tubers suggests that the same regulators might be involved in both processes. In an attempt to identify these regulators, comparative analyses of various transcription profiles were conducted. Firstly, features differentially expressed between a stolon/swollen stolon, diurnally in leaves and between growing/non-growing tubers were selected. To further decrease the number of features, only those features differentially expressed under all conditions were selected which reduced the number to 1662 (Figure 13A). K-means clustering analysis was performed on these features and most starch biosynthetic genes were present in three clusters. These clusters also showed the expected expression patterns for starch biosynthetic genes (Figure 13B). A functional assignment on the 913 features present in these three clusters was performed which showed that 37% of the features were involved in metabolism and a further 9% storage protein related (Figure 13C). Transcription factors made up 6% of the features and these were further investigated.

Figure 13

Comparative analysis of transcription profiles. A) Venn diagramm showing 1662 features differentially expressed under all conditions selected B) K means clustering showing 5 clusters (A-E) with three clusters having the desired pattern of expression for starch biosynthetic genes (B-D). C) Functional assignment of features present in clusters B-D).

Arabidopsis orthologs for these transcription factors were identified and used to gather more information on their possible function. Six of these transcription factors showed increased expression in Arabidopsis suspension culture fed with sucrose [50] and most of these genes were also down-regulated in leaves when the dark period was extended [51]. Figure 14 contain the expression profiles of these transcription factors including the identifiers of their orthologs in Arabidopsis. All of these transcription factors, with the exception of AT4G00870 and AT4G32730, were also down-regulated in the 35S::amiR-white-1(MIR172a) knockout [52]. The construct targets the GUN4 allele which is important for plastid to nucleus signalling and when this is interrupted it affects plastidic expression of several nuclear encoded genes. Furthermore, the expression of GPT2, plastidial PGM, GBSS and SBEB was also significantly down in these plants. This, along with the fact that these transcription factors were co-expressed in the same clusters as starch biosynthetic genes, indicate that they might be important in regulating starch accumulation.

Figure 14

Expression profiles of transcription factors possibly regulating starch biosynthesis. A) Diurnally in leaves, B) during tuber induction and C) growing and non-growing tubers. Micro.8007 (AT2g40820), orange, Micro.7865.c1 (AT3G16280), blue, Micro.5635.c1 (AT4G37750), black, Micro.5579.c2 (AT4G34590), light blue, Micro.15471.c1 (AT4G00870), brown and Micro.4326.c1 (AT4G32730), red.

Conclusion

Comparative transcriptome analysis between leaves and tubers indicate that transient and storage starch biosynthesis might not be all that different with the same isoforms being differentially expressed in both tissues. There was also a diurnal rhythm of GBSS expression in tubers which could be correlated to sucrose supply from the leaves. This provided evidence not only of the diurnal regulation of starch biosynthetic gene expression in tubers, but also showed the importance of sucrose supply in regulating gene expression in tubers.

Since tuber initiation and growth is not synchronised, it was important to determine the growth velocity of individual tubers. To this end X-ray CT was used to determine the volume of individual tubers at different time points and calculate the growth velocity. This was the first time that the growth velocities of tubers were determined in a natural environment. Tuber growth velocity could not be correlated to starch biosynthetic gene expression, although it was clear that gene expression is different between growing and non-growing tubers. The relationship between gene expression and growth velocity seems to be qualitative rather than quantitative and the data provides important information towards the identification of factors determining potato tuber growth.

Comparative analysis made it possible to select for genes differentially regulated under various conditions. Since the microarray experiments performed were set up to comprise conditions of active starch biosynthesis, it was believed that comparative analysis of transcription profiles would select for genes involved in this process. Cluster analysis of differentially expressed genes revealed clusters containing genes known to be involved in starch biosynthesis, and further analysis revealed transcription factors which could be used to influence starch biosynthesis. Closer analysis reveal that orthologs of these transcription factors in Arabidopsis are positively regulated by sucrose indicating that they could be interesting targets for influencing starch biosynthesis.

Methods

Plants and growth conditions

Solanum tuberosum (cv Solara) were propagated in tissue culture on MS medium [53] containing 2% sucrose. To obtain tubers, plants were transferred to soil and cultivated until harvest in individual pots in the greenhouse or growth chambers. For transcriptional analysis of potato leaves over a diurnal period, plants were grown in a growth chamber under a 14 hour light and 10 hour dark cycle. Plants for the tuber induction and X-ray CT studies were grown in the greenhouse under normal long day conditions. Plants used for analysing for the diurnal rhythm in tuber gene expression were grown in growth chambers under normal long day conditions.

RNA isolation, cRNA synthesis and Cy3-labeling

Isolation of total RNA was performed as described previously [54].

POCI array and database

The construction of the POCI array and database has been described in detail [29].

Sample preparation and microarray hybridization

For each hybridisation at least two biological replicates were included. For diurnal leaf time-points five different leaves from five plants were sampled for each biological replicate. For tuber induction samples stolon or tuber material from 40 plants were pooled into two groups according to the developmental stages. For the hybridization comparing fast and slow growing tubers, two independent samples were taken from each tuber and treated as biological replicates. RNA purity was measured by the ND-1000 Spectrophotometer (NanoDrop Technologies). To check for RNA degradation two μg of total RNA were separated on 1.5% formaldehyde containing agarose gel. Total RNA was purified using RNeasy Mini Spin Columns (QIAGEN, Valencia, CA). Afterwards RNA quality and quantity was tested using the Agilent 2100 BioAnalyzer (vB.02.03 BSI307) as recommended by manufacturer's protocol (Agilent RNA 6000 Nano Assay Protocol2). Synthesis of cDNA and cRNA was performed as described in the one-color microarray-based gene expression analysis protocol provided by Agilent including the one-color RNA spike-in kit (v5.0.1, 2006; Agilent Technologies, Santa Clara). After fragmentation Cy3-labelled samples were loaded on the array and hybridised over night (17h/65°C). Slides were washed as recommended in the manufacturer's protocol and scanned on the Agilent Microarray Scanner with extended dynamic range (XDR) at high resolution (5 μm). Data sets were extracted by using the feature extraction software (v9.5.3.1/Agilent Technologies) using a standard protocol.

Array data analysis

Data were imported into GeneSpring GX 7.3.1 (Silicon Genetics, Palo Alto, CA, USA) and additionally stored on a local server. A three step normalization was applied: (1) values less than five were set to five, (2) per chip normalization to 50^th percentile and (3) subsequently the signal for each feature was normalised to the median of its value across the entire dataset. For comparative analysis of various transcription profiles, a volcano plot was applied to select for features more than two-fold differentially expressed between two conditions including the Benjamini-Hochberg multiple test correction for the four replicate experiment comparing slow and fast growing tubers. No multiple test correction was employed for the experiments comparing stolon and swollen stolon, and the diurnal leaf experiment, where only two replicates per data-point were used. Features commonly differentially expressed in all experiments were identified using a Venn diagram. K-means clustering was performed using Pearson correlation to split the selected features into five clusters from where features were chosen to conduct a functional assignment.

Starch and sucrose measurements

Sucrose and starch were measured according to a modified method of Müller-Röber et al [17]. To determine sucrose supply rate to a tuber, sucrose content of stolon-ends was determined. A stolon-end is defined as the 15 millimetres of a stolon directly above the connection to the tuber.

qPCR analysis

For all qPCR analysis at least three biological repeats were used unless stated otherwise. Expression levels of genes were determined by real-time quantitative RT-PCR and the corresponding primers for the amplification of targets between 75 and 150 bp were designed using Primer3plus software [55]. Total RNA (five μg) from each of the developmental time points was treated with DNaseI (Fermentas GmbH) before undergoing reverse transcription using oligo d(T) primers and RevertAid™ H minus first strand cDNA synthesis kit (Fermentas GmbH) to generate a first strand cDNA template. Potato ubiquitin primers were used as a control as described previously [26]. One μl of 1:10 diluted cDNA for each time point were amplified with gene-specific primers in three technical replicates on a Mx3000P Q-PCR system (Stratagene) in combination with the Brilliant II SYBR Green Q-PCR Master Mix Kit (Stratagene). The thermal profile was as follows: 1 cycle 10 min at 95°C for DNA polymerase activation followed by 35 cycles of 10 s at 95°C, 15 s 60°C and 20 s 72°C. The primer sequences were as follows: ubi3 (L22576) forward primer, 5'-TTCCGACACCAT CGACAATGT-3'; reverse primer, 5'-CGACCATCCTCAAGCTGCTT-3'. For GBSS the primer sequence was based on POCI feature micro.920.c2. The forward primer was designated GBSS_920.c2 F (5' - CAGACTTGAGGAGCAGAAAGG - 3') and the reverse primer GBSS_920.c2 R (5' - GTGAGCCAAAGGGACATTGA - 3'). For GPT the primer sequence was based on POCI feature micro.1076.c1. The forward primer was designated GPT_1076.c1 F (5' -CCTTGTTTCCTGTTGCTGTG- 3') and the reverse primer GBSS_1076.c1 R (5' -AAAGCAGGCTCTCCACTCTT- 3'). For Susy the primer sequence was based on POCI feature micro.196.c8. The forward primer was designated Susy_196.c8 F (5' -CTGCTGTTTATGGGTTCTGG- 3') and the reverse primer Susy_196.c1 R (5' -GGCACACCTTCATTCACTCA- 3'). For NTT the primer sequence was based on POCI feature micro.1831.c2. The forward primer was designated NTT_1831.c2 F (5' -GAGCAGCAGCCAAGATAACAC- 3') and the reverse primer GBSS_1831.c2 R (5' -GTTCTGCATTGCACCCACA- 3'). Relative gene expression was calculated using the Pfaffl method [56].

Description of X-ray CT

With X-ray CT the 3D volume information of objects can be reconstructed using X-ray projections of the object from different aspects. The geometry used for the investigation was the axial 3D-CT, where a conical X-ray beam projects the object onto a flat 2D image detector. Using axial 3D-CT, projections of the object are taken under different viewing angles, rotating the object perpendicular to the central X-ray beam. The reconstructed volume data set consists of volumetric elements, called voxels, containing grey levels which represent information about the X-ray attenuation characteristics depending on the mass attenuation coefficient and the density distribution of the material [57, 58]. The mass attenuation coefficient itself is dependent on the applied X-ray spectrum and the effective atomic number of the X-rayed material. The calculation of tuber volumes of potted potato plants embedded in soil requires the segmentation of tubers from other materials in the X-ray CT volume data. Therefore a careful selection of exposure conditions is necessary to achieve sufficient data quality. This comprises X-ray parameters and filters as well as the condition of the soil. The parameters were defined as such: The X-ray source was FXE-225.45, accelerating voltage 200 kV, total emission 200 μA, the filter 1 mm Cu, detector Perkin Elmer, 1024 × 1024 pixel 200 μm, scan period 50 minutes and a resolution of 141 μm. After segmentation voxel elements, which have a known volume and specific grey level for tubers, were used to determine tuber volume by calculating how many voxels are present in a reconstructed tuber image. The lower and upper threshold grey values were set at 2938 and 3963 after which mistakes were corrected for manually. Tuber volume calculations were performed using Image J software http://rsbweb.nih.gov/ij.