BACKGROUND
the transition from active growth to dormancy in woody perennials of the temperate and boreal regions involves molecular, cellular, and whole-plant responses.
triggered by gradual decreases in temperature and day length, this remodeling reflects a complex array of adaptations to the physical and metabolic stress imposed by freezing temperatures.
substantial changes in the transcriptome and metabolome accompany this transition and enable the plant to avoid cellular damage resulting from direct thermal effects of freezing temperatures on macromolecules, extracellular ice formation, and the generation of reactive oxygen species .
in deciduous angiosperms and evergreen gymnosperms, the light- and hormone-regulated transcriptional and metabolic changes associated with cambial and vegetative bud dormancy and freezing tolerance are initiated by prolonged nights long prior to the actual freezing  <cit> – <cit> .
because most conifers retain their leaves  during winter, protection from freezing has to extend to these photosynthetic organs.
the photosynthetic apparatus is prone to oxidative damage as temperatures drop due to energetic and metabolic imbalance  <cit> .
this imbalance is a result of the temperature dependency of metabolic processes, including photochemistry.
light is absorbed by photosystems regardless of the temperature, but enzymes are inhibited at low temperatures, resulting in over-reduction of the photosystems and photoinhibition.
high-energy electrons produced in excess reduce molecular oxygen to generate ros causing photosystem damage  <cit> .
in evergreen conifers, such oxygen reduction in a mehler-type reaction is coupled to antioxidant defenses and represents an important mechanism of dissipating excess energy absorbed by photosystem i at low temperatures  <cit> .
photosynthetic acclimation to low temperatures and freezing in evergreen gymnosperms also involves other mechanisms, including the re-localization of chloroplasts, reorganization and aggregation of photosystem antennae, chlororespiration, and non-photochemical quenching leading to conversion of absorbed light energy to heat rather than to reductant for co <dig> fixation and growth  <cit> – <cit> .
lodgepole pine  was shown to reduce the antenna size and the number of reaction centers in photosystem ii during winter hardening to minimize light absorption  <cit> .
additional mechanisms that facilitate the transition from energy harvesting to energy dissipation that occurs with winter hardening include changes in abundance of specific thylakoid proteins  <cit> , non-photochemical quenching via the xanthophyll cycle, and antenna protonation via thylakoid ph changes to dissipate the excess energy to heat  <cit> .
a second, zeaxanthin-independent, quenching mechanism has also been described, involving charge recombination between photosystem ii reaction center components  <cit> .
photosystem i is less sensitive to low temperatures than photosystem ii and supports xanthophyll-mediated non-photochemical quenching, while maintaining active cyclic electron transport for atp synthesis at low temperatures  <cit> .
chlororespiration and cyclic electron transport enable the dissipation of overproduced excitation energy via alternative electron acceptors, including the terminal oxidase in the plastid , a bi-functional protein involved in carotenoid biosynthesis and oxidation of plastoquinol produced by over-reduced electron transport chain   <cit> – <cit> .
potential excess nadph produced as a result of linear photosynthetic electron transport in the almost complete absence of carbon fixation in the stroma can also be dissipated with or without altering atp synthesis through the action of the alternate etc.
in mitochondria via operation of the plastidic malate/oxaloacetate shuttle  <cit> .
the alternate mitochondrial etc.
also provides means of preventing ros formation during mitochondrial aerobic respiration  <cit> .
if ros scavenging does not occur rapidly enough, under such conditions as extreme biotic or abiotic stresses including low temperatures, damage to mitochondrial and other cellular components and even programmed cell death may occur  <cit> – <cit> .
ros production can also be prevented by lowering o <dig> levels and/or redirecting respiratory electron flux to other substrates.
these alternate routes for electrons during respiration involve alternative oxidases, external and internal alternative nadh dehydrogenases, and uncoupling proteins, that are especially active upon imposition of abiotic stress that leads to ros production  <cit> – <cit> .
these enzymes play a key role in plant acclimation and tolerance to low temperatures by preventing ros formation when the traditional route involving complexes iii and iv is inhibited  <cit> – <cit> .
these alternate routes bring about a lowering of ros levels by maintaining active electron flow and preventing over-reduction of electron transport components, or redox imbalance  <cit> .
alternative oxidases compete for electrons with cytochrome c oxidase , using oxygen as a substrate.
in addition to maintaining mitochondrial electron flow, they can also prevent ros production directly by lowering the levels of oxygen, which can react with electrons to generate ros, up to 60% during abiotic stresses  <cit> .
rotenone-insensitive alternative internal and external nadh dehydrogenases interfere with mitochondrial proton transport along with alternative oxidases  <cit> – <cit> .
these dehydrogenases bypass complex i, while alternative oxidases compete for the electrons with complex iii, thus preventing a total of three protons from being transported from the matrix to the inter-membrane space by nadh dehydrogenase and cytochrome c reductase and oxidase  <cit> .
unlike alternative oxidases and dehydrogenases, uncoupling proteins diminish the electrochemical potential directly by transporting protons back to the matrix with the result that less atp is synthesized.
uncoupling of the etc.
from oxidative phosphorylation stimulates fluxes through the etc.
and subsequently substrate oxidation during carbon and energy imbalance that occurs during a variety of stresses.
as such, uncoupling proteins are also involved in controlling ros levels and promoting tolerance of plants to stresses  <cit> – <cit> .
the transition from energy harvesting to dissipation of the photosynthetic apparatus in needles that takes place during winter hardening is a major adaptive response.
however, freezing acclimation also involves a wide range of other cellular adaptations, which requires a substantial source of metabolic precursors and chemical energy that, given the circumstances, clearly cannot originate from linear photosynthetic electron transport and co <dig> fixation.
for example, detailed bioinformatic analyses of transcriptomic and metabolic remodeling during winter hardening in sitka spruce  revealed the presence of highly complex regulatory networks involving ethylene-mediated signaling, active trafficking of cell wall carbohydrates through the endomembrane system, and significant remodeling of cell walls  <cit> .
here we implement systems biology tools to further interrogate these datasets, which involve five time points sampled from late summer to early winter  <cit> .
specifically, the current paper focuses on adaptation of the photosynthetic apparatus to low temperature, with the focus on whole processes relevant to energy metabolism from the global perspective of systems biology rather than individual genes.
RESULTS
photosynthesis and respiration
the vast majority of responsive genes associated with chloroplast function were down-regulated after the time point  <dig>  and remained in that state for the duration of the time course .
notable exceptions were ptox/immutans, a nuclear gene encoding a plastid terminal oxidase, and psb, a chloroplast gene encoding the cp <dig> subunit of the reaction center of photosystem ii , in addition to the plastid-targeted antioxidant gene products discussed below .
in contrast, respiration-related genes were up-regulated .
for example, atmg <dig>  which encodes a subunit of mitochondrial nadh dehydrogenase , showed a consistent about 3-fold up-regulation throughout the time-course, and other components of the etc.
exhibited low, but consistent up-regulation, and clustered with the genes encoding alternative mitochondrial etc.
components.
however, there were also corresponding sets of genes involved in respiration  that were down-regulated .
antioxidant function
genes encoding proteins related to antioxidant activity in various organelles responded in diverse ways during the time course .
massive increases in expression levels of peroxisomal catalase  <dig> were the most striking results in this overall category.
increases in the expression of genes encoding differentially targeted thioredoxin-related proteins were also observed, including three genes encoding protein disulfide isomerases  targeted to the er, and genes homologous to atero, which interacts with pdi in the er.
genes encoding proteins that function in the ubiquitous ros-scavenging ascorbate/glutathione pathway showed varying responses.
glutathione reductase  <dig>  encoding a gene product targeted to the peroxisome and the cytosol, was up-regulated.
two monodehydroascorbate reductases, which generate ascorbate from monodehydroascorbate substrates, were up-regulated.
one of these proteins is targeted to the chloroplast stroma.
two genes encoding ascorbate peroxidases were also up-regulated.
after cat <dig>  thylakoid-associated l-ascorbate peroxidase was among the strongest up-regulated antoxidant genes.
re-reduction of ascorbate using reduced glutathione as a substrate was evident by the up-regulation of a dehydroascorbatereductase , as well as an intrinsic trans-membrane cytochrome b <dig>  and a membrane-associated heme binding protein that may be involved in intra-vesicular reduction of monodehydroascorbate.
within the glutaredoxin category, a spruce homolog of the disulfide oxidoreductase roxy <dig>  which has a regulatory role in angiosperm flower development, was highly up-regulated at the later time points.
notably up-regulated early in the time course were several genes involved in ascorbate biosynthesis, including an l-galactose dehydrogenase, which catalyzes the conversion of l-galactose to l-galactono- <dig> -lactone, and a gdp-mannose  <dig> -epimerase, which catalyzes an important step in the ascorbate biosynthetic pathway – the conversion of gdp-d-mannose to gdp-l-galactose.
the automated layout pipeline for inferred networks  pipeline was used to visualize documented “guilt by association” relationships for the thylakoid-bound ascorbate peroxidase  .
the “guilt by association” relationships are imported from genemania and drawn from a large dataset containing functional associations to tapx.
these functional associations include protein and genetic interactions, pathways, co-expression, co-localization and protein domain similarity.
the alpine visualization integrates the functional associations of genes, gene ontology  annotations, and subcellular location data into a single unified view.
this view is filtered by user supplied gene expression data to return only gene associations obtained from the input.
this ability to exclusively analyze a particular dataset and return only the gene associations of interest and significance is unique to alpine and to our knowledge, there are no bioinformatics tools available with the alpine’s useful and unique combination of functionalities and capabilities.this strategy revealed clues concerning antioxidant-related events in winter hardening that are not apparent using other singular strategies.
of the  <dig> genes that were found to be associated with tapx, all but five were down-regulated, and all but one gene encoded a plastid-targeted protein.
the plastid-targeted down-regulated genes fell into several categories, for example, photosystem ii, chloroplast-based protein synthesis, protein folding, and rna regulation of transcription.
in contrast, a gene encoding a member of the rhodanase family was up-regulated, as was a plastid-targeted, nadp-linked oxidoreductase.
taken together, tapx may be associated with the regulatory mechanism that brings about the virtual cessation of many chloroplast activities that occur during winter hardening.
all associations were found to be co-expression associations, and all  <dig> genes in the network demonstrate a significant  expression value in t <dig> of the dataset.
alternative etc
components of the alternative electron transport are known to interfere with mitochondrial atp production.
two external alternative nadh dehydrogenases showed a consistent up-regulation in response to freezing temperatures across the time series, while two internal ones from tp <dig> to tp <dig> .
alternative oxidase a <dig> gene  was up-regulated at tp <dig> and  <dig>  while an ortholog of uncoupling protein  <dig>  at tp <dig> only .
carbon and reductant sources for respiration
to determine the source of reductant for mitochondrial electron transport, we assessed changes in gene expression of various genes involved in central carbon metabolic pathways that could provide these substrates .
the only potentially active pathways involved sugars.
transcript levels of four invertase genes, which are involved in starch degradation, increased, including a cell wall invertase, whose expression level increased through tp <dig> .
a homolog of at4g <dig>  a cytosolic invertase, showed large increases through tp <dig>  the expression of homologs to members of the sucrose synthase  gene family showed a somewhat analogous pattern with sus <dig> exhibiting a substantial increase through tp <dig>  sus <dig> also showed an increase, albeit a more modest one, which peaked between tp <dig> and tp <dig>  in contrast, sus <dig> showed a decrease in tp <dig> and no significant changes at other tps.
however, sucrose levels remained constant and above baseline levels during freezing acclimation.
more generally, genes encoding starch biosynthesis enzymes were moderately up-regulated, while those associated with degradation were down-regulated .
adp glucose pyrophosphorylase  involved in starch biosynthesis is a heteromer of two large and two small subunits.
apl <dig> and  <dig>  two of the three genes encoding the large agpase subunits responded and showed increases in the expression from tp <dig> to tp <dig>  a gene encoding the small subunit of agpase  showed a decrease in expression at tp <dig>   <dig> and  <dig> and an increase in tp <dig>  one gene encoding starch synthase showed an initial decrease in tp <dig> and  <dig>  while the other responsive starch synthase gene was down-regulated at tp <dig> .
most genes encoding glycolytic enzymes did not have considerably altered levels of relevant transcripts in response to freezing .
however, some glycolytic genes showed a more substantial up- or down-regulation, including the pyrophosphate-dependent phosphofructokinase  <dig>  that showed a maximal  <dig> -fold increase and a constant up-regulation.
the same step in glycolysis can also be catalyzed by another pyrophosphate-dependent enzyme, fructose-6-p 1-phosphotransferase and the corresponding α - subunit  was correspondingly down-regulated .
fructose- <dig> -bisp aldolase  transcript levels increased consistently  in response at all tps.
in arabidopsis, at5g <dig> encodes a protein similar to cytosolic pyruvate kinase, and the corresponding gene in sitka spruce showed a 3-fold increase in transcript levels.
transcript levels of gluconeogenetic pyruvate pyrophosphate dikinase  decreased at tp <dig> and remained low throughout the time course.
phosphoenolpyruvatecarboxykinase  transcript gradually decreased  <dig>  fold through tp <dig>
dihydrolipoamide dehydrogenase  <dig> , the e <dig> subunit  <dig> of pyruvate dehydrogenase, was moderately , but significantly up-regulated from tp <dig> to tp <dig> .
in contrast, lpd <dig> , the e <dig> subunit  <dig> of pyruvate dehydrogenase, was down-regulated about  <dig> -fold between tp <dig> and tp <dig>  the steady-state transcript levels for other subunits of pyruvate dehydrogenase remained unchanged during freezing.
this was true for most of the transcripts encoding the remaining enzymes of the tca cycle, though they tended to be slightly above the baseline levels  throughout the time course.
a putative cytosolic nadp-dependent isocitrate dehydrogenase  showed the same pattern and degree of increases in gene expression as lpd <dig>  a putative regulatory subunit of the mitochondrial nad-dependent isocitrate dehydrogenase iii  increased 2-fold in tp <dig> and showed a stable  <dig> -fold increase in transcript levels in tp <dig> through tp <dig>  malic enzyme gene  showed only a moderate  <dig> -fold up-regulation in response to freezing, but the consistent 9- to 11-fold induction of the at2g13560-encoded nad-dependent malic enzyme transcript levels in tp <dig>   <dig>  and  <dig> was quite striking.
citrate and succinate levels were moderately increased in a stable manner, while malate levels gradually increased up to  <dig>  fold during freezing .
discussion and
CONCLUSIONS
temperate and boreal conifers must survive harsh winters while retaining their photosynthetic tissue, and have evolved intricate mechanisms to maintain a certain level of active metabolism, signaling, and protective processes at low temperatures in their needles.
the comprehensiveness of our data afforded a global view of winter hardening as it is reflected in numerous processes.
the use of bioinformatics tools enabled the visualization of these changes in a cellular context.
the results point to numerous photoprotective processes that are active during winter, serve to confirm previous reports, and provide new insights into the activity of genes associated with chloroplast function, energy metabolism, antioxidant processes, and possible regulatory mechanisms.
genes related to energy quenching mechanisms including a spruce psbs homolog and plastid terminal oxidase ptox increased their expression, which is in an agreement with previous reports  <cit> .
psbs activity is associated with changes in photosystem ii that occur as part of the photoprotective mechanisms that are set in place as the temperature drops, while ptox represents an important component of the dissipation of excess light energy in photosynthesis as an alternative electron acceptor and o2-utilizing enzyme helping to prevent ros formation in hardening conifer needles  <cit> .
in ranunculus glacialis, an alpine plant growing in high altitudes at high light and low temperatures, ptox protein levels correlated with the use of alternative etc.
and enabled high fluxes through photosynthetic etc.
despite the excess of light energy.
ptox was implicated to serve as a safety valve for the over-reduction of photosystems during high light and low temperatures in alpine plants  <cit> .
ptox has also been implicated in increasing the flux through the linear etc.
and subsequent increased tolerance of photosystems to high light in cold-hardened arabidopsis thaliana <cit> .
this result suggests that the energy dissipation mechanisms set in motion within the chloroplast during winter hardening were effective.
one manifestation of this may have been a relatively low level of ros generation as the temperature dropped.
however, redox effects on cellular responses to the environment are not confined to a simple containment of ros.
much data point to regulatory roles for specific antioxidant proteins in plant cells under stress.
several of these antioxidant, regulatory genes were up-regulated over the winter hardening time course, including catalase cat <dig>  glutaredoxin roxy <dig>  and glutathione reductase gr <dig>  catalases regulate the levels of hydrogen peroxide, which is an important signaling molecule in plant stress responses  <cit> .
roxy1encodes a glutaredoxin with a documented regulatory role in development in arabidopsis  <cit> – <cit> , while gr1exerts a regulatory effect on responses to cold stress, an effect that cannot be replaced by gr <dig>  or thioredoxins  <cit> .
the ros scavenging pathway may have been operational in the chloroplast as winter hardening proceeded, since other components, such as mondehydroascorbate reductase, were also up-regulated .
however, it is also possible that tapx was acting in a different/additional capacity to protect the photosynthetic machinery.
duan et al.
have shown that tapx plays a specific role in alleviating the photosystem i and ii-associated photoinhibition that occurs with chilling stress in arabidopsis  <cit> , and spruce tapx may play an analogous role in the chloroplasts of hardening needles.
interestingly, maruta et al.
provide evidence for a specific regulatory role for tapx in chloroplast responses to stress, resonant with the role shown for gr <dig>  <cit> .
in the case of tapx and the winter hardening data, the guilt by association analysis revealed an association of tapx exclusively with plastid-targeted gene products, in which the great majority of those responsive genes were down-regulated.
it is possible that tapx plays an important role in the adaptation mechanism that results in the onset of chloroplast dormancy in conifer needles.
taken together, these results suggest specific regulatory roles for tapx, roxy <dig>  and gr <dig> in the adaptation mechanism underlying winter hardening.
freezing temperatures impose severe stress on the mitochondrial etc.
a set of genes encoding components of the respiratory etc.
were strongly or mildly up-regulated, suggesting that the mitochondrial etc.
is active.
mitochondrial aerobic respiration is a source of ros formed through the reduction of molecular oxygen by electrons coming from complexes i and iii of the etc.
<cit> .
the ros production can be prevented by redirecting respiratory electron flux to other substrates by involving the major players of the alternate etc.
including alternative oxidases, external and internal alternative nadh dehydrogenases, and uncoupling proteins  <cit> – <cit> .
genes encoding all these components of the alternate etc.
were up-regulated during the time course, suggesting that these alternate pathways are active, along with the conventional etc., during freezing acclimation.
up-regulation of genes encoding the components of the alternate etc.
is typical during cold acclimation in plants in general and is accompanied with metabolic changes.
tobacco plants subjected to cold showed an increase in alternative oxidase levels and sugar content compared to the transgenic plants with suppressed expression of this enzyme  <cit> .
these transgenic plants also showed an increased lipid peroxidation at normal, but not at low, temperatures.
alternative oxidases are regulated at the transcriptional level by abscisic acid insensitive  <dig>  transcription factor and are activated by pyruvate  <cit> – <cit> .
accumulation of pyruvate as the final product of glycolysis indicates an imbalance between glycolytic fluxes and the mitochondrial etc.
pyruvate then activates alternative oxidases to stimulate flux through the mitochondrial etc.
to keep pace with glycolysis  <cit> .
suppressing the external nadh dehydrogenase ndb <dig> resulted in increased levels of ndb <dig> and alternative oxidase 1a in tobacco, which lead to increased salt tolerance and improved growth  <cit> .
glycolysis and the tca cycle as sources of carbon, atp, and reductant
avoiding ros production comes at a price, as every alternate route leads to a decreased atp production by mitochondrial atp synthase.
both normal and alternate etcs were active based on our data, which brings a question of the potential source of reductant for these etcs and source of atp that would supplement atp lost due to the potentially high activities of the alternate routes for electrons.
most glycolytic genes did not show major steady-state transcript differences during freezing acclimation, though persistent up-regulation was observed, suggesting that glycolysis was active.
in plants, glycolysis is regulated at the transcriptional and allosteric levels and in a complex way due to metabolic flexibility  <cit> .
in situations when atp cannot be easily produced in heterotrophic systems , pyrophosphate- instead of atp-utilizing enzyme forms are used to drive glycolysis  <cit> .
one of the genes encoding phosphofructokinase was up-regulated up to  <dig> -fold during winter hardening.
up-regulation of a pyrophosphate-dependent phosphofructokinase isoform suggests the possibility of conservation of atp during freezing acclimation, similar to phosphate starvation stress.
to maximize atp production, alternate glycolytic steps involving substrate-level phosphorylation are preferentially utilized.
the conversion of glyceraldehyde-3p to 3-phosphoglycerate can be done in two ways by:  non-phosphorylating glyceraldehyde-3p dehydrogenase in a single step to only produce nadph and  phosphoglycerate kinase and phosphorylating glyceraldehyde-3p dehydrogenase in two steps involving nadh and atp production  <cit> .
the phosphorylating glyceraldehyde-3p dehydrogenase was moderately, but consistently up-regulated, while the kinase showed a consistent, less than  <dig> -fold down-regulation, suggesting that the two-step reaction producing both the reductant and atp could be active.
atp production by substrate-level phosphorylation appears to be active in general – three transcripts coding for pyruvate kinase showed varying, but consistent up-regulation.
pyruvate kinase catalyzes a key step regulating glycolysis and pyruvate is a known allosteric regulator of a number of glycolytic enzymes including pyruvate kinase, activity of which is induced during hypoxia when glycolysis is the major source of atp  <cit> .
based on these observations, it appears that glycolytic routes conserving or producing atp are preferentially used.
this is in an agreement with the global activation of alternative etc.
routes.
as such, glycolysis could be a source of atp supplementing lower atp levels obtained from respiration and could help with high-energy cofactor balancing.
in conventional glycolysis, pyruvate kinase converts phosphoenolpyruvate to pyruvate and atp is made.
however, mitochondria could oxidize substrates other than pyruvate, such as malate  <cit> , during winter hardening.
in this scenario, the pyruvate kinase step is bypassed by three consecutive steps with no net nadh and atp production or consumption .
however, this bypass allows nadh from the cytosol to be shuttled inside mitochondria where it can be used by the conventional or alternative etcs.
first, cytosolic phosphoenolpyruvate is carboxylated to oxaloacetate by phosphoenolpyruvate carboxylase, a highly regulated enzyme whose activity increases during phosphate starvation  <cit> – <cit> .
second, oxaloacetate is reduced to malate by cytosolic nadh-utilizing malate dehydrogenase, and malate is transported to mitochondria.
the dicarboxylate/tricarboxylate transporter dtc encoded by the arabidopsis at5g <dig> gene ortholog that is capable of transporting malate to mitochondria  <cit>  showed a consistent 5-fold up-regulation.
third, malate either enters the tca cycle or is decarboxylated to pyruvate by nad-dependent malic enzyme in mitochondria.
most genes encoding the enzymes of the tca cycle showed a moderate up-regulation or no significant changes in their steady-state transcript levels, suggesting that the tca cycle is operational.
for example, a putative regulatory subunit of the mitochondrial nad-dependent isocitrate dehydrogenase iii   <cit>  showed a stable  <dig> -fold increase in transcript levels in tp <dig> through tp <dig>  malic enzyme has been implicated in generating nadh and pyruvate needed for lipid synthesis in heterotrophic and photoheterotrophic systems such as plant seeds  <cit> – <cit> .
the physiological significance of this bypass in most in vivo situations is not well understood  <cit> .
however, existing metabolomic and transcriptomic data  <cit>  suggest that the bypass is operational  during freezing acclimation.
malate levels increased significantly at tp <dig> and reached the levels that were nearly 3-fold higher than baseline later in the time course.
two genes encoding phosphoenolpyruvate carboxylase and one gene encoding malic enzyme were moderately up-regulated in response to freezing acclimation, while cytosolic malate dehydrogenase transcript levels were not significantly different in all tps.
however, the sitka spruce homolog of the arabidopsis at2g13560-encoded nad-dependent malic enzyme  <dig>  which is allosterically activated by fumarate  <cit> , showed a 9- to 11-fold up-regulation in its transcript levels during freezing acclimation and could drive significant fluxes through this bypass.
low temperature stress was also shown to induce the expression of a gene encoding malic enzyme in maize  <cit> .
to drive efficient conventional and alternative etcs during winter hardening, this bypass could:  shuttle nadh from cytosol to mitochondria,  bring additional carbon from glycolysis to the tca cycle to stimulate alternative and conventional mitochondrial etcs, and  generate phosphate for atp synthesis.
potential sources of hexoses for glycolysis
the reaction catalyzed by invertase requires  <dig> atp molecules to generate fructose-6p, while the net ntp consumption for the sucrose synthase pathway is zero and the latter is up-regulated during hypoxia to save atp  <cit> .
both invertases and sucrose synthases showed significant increases in the corresponding transcript levels in sitka spruce needles during freezing acclimation.
however, steady-state sucrose levels did not change significantly.
so, if there were an increase in sucrose degradation to provide glucose and fructose for glycolysis, there would be a corresponding increase in sucrose supply to maintain unchanging sucrose levels.
because photosynthetic carbon fixation is close to inactive, one can hypothesize that sucrose would have to be transported from organs that have a great potential to store it or made from other, more condensed form of carbohydrate such as starch that can easily be converted to sucrose.
starch can accumulate transiently in photosynthetically active cells or be stored for a prolonged period of time in tubers and seeds  <cit> – <cit> .
this makes starch an excellent candidate as a storage compound that could be mobilized and used as a source of hexoses during winter hardening.
however, the majority of genes associated with starch degradation were down-regulated, while those involved in its biosynthesis were up-regulated, which would make starch present in needles a seemingly unlikely storage form of carbon.
regulation of starch metabolism is complex, as it involves several levels of regulation, including transcriptional, protein phosphorylation, allosteric, and redox regulation, to contribute to the balance between photosynthate production and demand in source and sink tissues  <cit> – <cit>  and as such, transcript levels alone cannot provide enough information regarding the direction of starch metabolism during winter hardening.
changes in starch levels were not measured during freezing acclimation in sitka spruce needles  <cit> , but fisher &holl  <cit>  showed that starch content progressively decreased in scots pine needles during the winter, suggesting that this polysaccharide could possibly serve as a storage compound for the winter hardening processes also in sitka spruce needles despite the up-regulation in starch synthesis genes.
however, an increase in starch synthesis was also observed in arabidopsis leaves exposed to prolonged moderate osmotic stress  <cit>  that lead to similar global gene expression responses involving cell wall remodeling, vesicular transport, and hormonal signaling as in winter hardening  <cit> , but the reason for these increases in starch synthesis remains to be elucidated.
addressing the question of the storage form of carbon in needles during winter hardening is quite complicated.
starch is a potential candidate, but conifers also accumulate di and triacylglycerols as storage compounds in the needles  <cit> .
unfortunately these compounds were not measured in sitka spruce needles during winter hardening  <cit> , but both types of lipids showed a similar trend as their levels gradually increased during fall and winter in scots pine needles, followed by a rapid decrease in spring prior to bud break, suggesting that lipid degradation provides additional carbon and energy for new buds before photosynthesis is fully operational  <cit> .
it is unclear whether these lipids can also be used as a source of energy during winter hardening.
starch levels decreasing and lipid levels increasing in needles during the winter is consistent with the hypothesis that starch rather than lipids are substrates for respiration and active cell-wall remodeling during winter.
in addition, degradation of oligo and polysaccharides is more energy effective for producing hexoses for substantial cell-wall remodeling than potential gluconeogenesis  from lipids.
a possibility also exists that soluble and insoluble sugars are stored in other parts of the trees and then transported in the form of sucrose to needles to help maintaining active metabolism during prolonged freezing.
all these compounds including di and triacylglycerols are present in the sapwood of scots pine  <cit> , but detailed source-sink relationship studies are needed to confirm the possibility of sucrose transport to needles in winter.

 the transition from active growth to dormancy in woody perennials of the temperate and boreal regions involves molecular, cellular, and whole-plant responses. triggered by gradual decreases in temperature and day length, this remodeling reflects a complex array of adaptations to the physical and metabolic stress imposed by freezing temperatures. substantial changes in the transcriptome and metabolome accompany this transition and enable the plant to avoid cellular damage resulting from direct thermal effects of freezing temperatures on macromolecules, extracellular ice formation, and the generation of reactive oxygen species. in deciduous angiosperms and evergreen gymnosperms, the light- and hormone-regulated transcriptional and metabolic changes associated with cambial and vegetative bud dormancy and freezing tolerance are initiated by prolonged nights long prior to the actual freezing –. because most conifers retain their leaves during winter, protection from freezing has to extend to these photosynthetic organs. the photosynthetic apparatus is prone to oxidative damage as temperatures drop due to energetic and metabolic imbalance. this imbalance is a result of the temperature dependency of metabolic processes, including photochemistry. light is absorbed by photosystems regardless of the temperature, but enzymes are inhibited at low temperatures, resulting in over-reduction of the photosystems and photoinhibition. high-energy electrons produced in excess reduce molecular oxygen to generate ros causing photosystem damage. in evergreen conifers, such oxygen reduction in a mehler-type reaction is coupled to antioxidant defenses and represents an important mechanism of dissipating excess energy absorbed by photosystem i at low temperatures. photosynthetic acclimation to low temperatures and freezing in evergreen gymnosperms also involves other mechanisms, including the re-localization of chloroplasts, reorganization and aggregation of photosystem antennae, chlororespiration, and non-photochemical quenching leading to conversion of absorbed light energy to heat rather than to reductant for co fixation and growth –. lodgepole pine was shown to reduce the antenna size and the number of reaction centers in photosystem ii during winter hardening to minimize light absorption. additional mechanisms that facilitate the transition from energy harvesting to energy dissipation that occurs with winter hardening include changes in abundance of specific thylakoid proteins, non-photochemical quenching via the xanthophyll cycle, and antenna protonation via thylakoid ph changes to dissipate the excess energy to heat. a second, zeaxanthin-independent, quenching mechanism has also been described, involving charge recombination between photosystem ii reaction center components. photosystem i is less sensitive to low temperatures than photosystem ii and supports xanthophyll-mediated non-photochemical quenching, while maintaining active cyclic electron transport for atp synthesis at low temperatures. chlororespiration and cyclic electron transport enable the dissipation of overproduced excitation energy via alternative electron acceptors, including the terminal oxidase in the plastid, a bi-functional protein involved in carotenoid biosynthesis and oxidation of plastoquinol produced by over-reduced electron transport chain –. potential excess nadph produced as a result of linear photosynthetic electron transport in the almost complete absence of carbon fixation in the stroma can also be dissipated with or without altering atp synthesis through the action of the alternate etc. in mitochondria via operation of the plastidic malate/oxaloacetate shuttle. the alternate mitochondrial etc. also provides means of preventing ros formation during mitochondrial aerobic respiration. if ros scavenging does not occur rapidly enough, under such conditions as extreme biotic or abiotic stresses including low temperatures, damage to mitochondrial and other cellular components and even programmed cell death may occur –. ros production can also be prevented by lowering o levels and/or redirecting respiratory electron flux to other substrates. these alternate routes for electrons during respiration involve alternative oxidases, external and internal alternative nadh dehydrogenases, and uncoupling proteins, that are especially active upon imposition of abiotic stress that leads to ros production –. these enzymes play a key role in plant acclimation and tolerance to low temperatures by preventing ros formation when the traditional route involving complexes iii and iv is inhibited –. these alternate routes bring about a lowering of ros levels by maintaining active electron flow and preventing over-reduction of electron transport components, or redox imbalance. alternative oxidases compete for electrons with cytochrome c oxidase, using oxygen as a substrate. in addition to maintaining mitochondrial electron flow, they can also prevent ros production directly by lowering the levels of oxygen, which can react with electrons to generate ros, up to 60% during abiotic stresses. rotenone-insensitive alternative internal and external nadh dehydrogenases interfere with mitochondrial proton transport along with alternative oxidases –. these dehydrogenases bypass complex i, while alternative oxidases compete for the electrons with complex iii, thus preventing a total of three protons from being transported from the matrix to the inter-membrane space by nadh dehydrogenase and cytochrome c reductase and oxidase. unlike alternative oxidases and dehydrogenases, uncoupling proteins diminish the electrochemical potential directly by transporting protons back to the matrix with the result that less atp is synthesized. uncoupling of the etc. from oxidative phosphorylation stimulates fluxes through the etc. and subsequently substrate oxidation during carbon and energy imbalance that occurs during a variety of stresses. as such, uncoupling proteins are also involved in controlling ros levels and promoting tolerance of plants to stresses –. the transition from energy harvesting to dissipation of the photosynthetic apparatus in needles that takes place during winter hardening is a major adaptive response. however, freezing acclimation also involves a wide range of other cellular adaptations, which requires a substantial source of metabolic precursors and chemical energy that, given the circumstances, clearly cannot originate from linear photosynthetic electron transport and co fixation. for example, detailed bioinformatic analyses of transcriptomic and metabolic remodeling during winter hardening in sitka spruce revealed the presence of highly complex regulatory networks involving ethylene-mediated signaling, active trafficking of cell wall carbohydrates through the endomembrane system, and significant remodeling of cell walls. here we implement systems biology tools to further interrogate these datasets, which involve five time points sampled from late summer to early winter. specifically, the current paper focuses on adaptation of the photosynthetic apparatus to low temperature, with the focus on whole processes relevant to energy metabolism from the global perspective of systems biology rather than individual genes.  photosynthesis and respiration the vast majority of responsive genes associated with chloroplast function were down-regulated after the time point and remained in that state for the duration of the time course. notable exceptions were ptox/immutans, a nuclear gene encoding a plastid terminal oxidase, and psb, a chloroplast gene encoding the cp subunit of the reaction center of photosystem ii, in addition to the plastid-targeted antioxidant gene products discussed below. in contrast, respiration-related genes were up-regulated. for example, atmg which encodes a subunit of mitochondrial nadh dehydrogenase, showed a consistent about 3-fold up-regulation throughout the time-course, and other components of the etc. exhibited low, but consistent up-regulation, and clustered with the genes encoding alternative mitochondrial etc. components. however, there were also corresponding sets of genes involved in respiration that were down-regulated. antioxidant function genes encoding proteins related to antioxidant activity in various organelles responded in diverse ways during the time course. massive increases in expression levels of peroxisomal catalase were the most striking results in this overall category. increases in the expression of genes encoding differentially targeted thioredoxin-related proteins were also observed, including three genes encoding protein disulfide isomerases targeted to the er, and genes homologous to atero, which interacts with pdi in the er. genes encoding proteins that function in the ubiquitous ros-scavenging ascorbate/glutathione pathway showed varying responses. glutathione reductase encoding a gene product targeted to the peroxisome and the cytosol, was up-regulated. two monodehydroascorbate reductases, which generate ascorbate from monodehydroascorbate substrates, were up-regulated. one of these proteins is targeted to the chloroplast stroma. two genes encoding ascorbate peroxidases were also up-regulated. after cat thylakoid-associated l-ascorbate peroxidase was among the strongest up-regulated antoxidant genes. re-reduction of ascorbate using reduced glutathione as a substrate was evident by the up-regulation of a dehydroascorbatereductase, as well as an intrinsic trans-membrane cytochrome b and a membrane-associated heme binding protein that may be involved in intra-vesicular reduction of monodehydroascorbate. within the glutaredoxin category, a spruce homolog of the disulfide oxidoreductase roxy which has a regulatory role in angiosperm flower development, was highly up-regulated at the later time points. notably up-regulated early in the time course were several genes involved in ascorbate biosynthesis, including an l-galactose dehydrogenase, which catalyzes the conversion of l-galactose to l-galactono- -lactone, and a gdp-mannose -epimerase, which catalyzes an important step in the ascorbate biosynthetic pathway – the conversion of gdp-d-mannose to gdp-l-galactose. the automated layout pipeline for inferred networks pipeline was used to visualize documented “guilt by association” relationships for the thylakoid-bound ascorbate peroxidase. the “guilt by association” relationships are imported from genemania and drawn from a large dataset containing functional associations to tapx. these functional associations include protein and genetic interactions, pathways, co-expression, co-localization and protein domain similarity. the alpine visualization integrates the functional associations of genes, gene ontology annotations, and subcellular location data into a single unified view. this view is filtered by user supplied gene expression data to return only gene associations obtained from the input. this ability to exclusively analyze a particular dataset and return only the gene associations of interest and significance is unique to alpine and to our knowledge, there are no bioinformatics tools available with the alpine’s useful and unique combination of functionalities and capabilities.this strategy revealed clues concerning antioxidant-related events in winter hardening that are not apparent using other singular strategies. of the genes that were found to be associated with tapx, all but five were down-regulated, and all but one gene encoded a plastid-targeted protein. the plastid-targeted down-regulated genes fell into several categories, for example, photosystem ii, chloroplast-based protein synthesis, protein folding, and rna regulation of transcription. in contrast, a gene encoding a member of the rhodanase family was up-regulated, as was a plastid-targeted, nadp-linked oxidoreductase. taken together, tapx may be associated with the regulatory mechanism that brings about the virtual cessation of many chloroplast activities that occur during winter hardening. all associations were found to be co-expression associations, and all genes in the network demonstrate a significant expression value in t of the dataset. alternative etc components of the alternative electron transport are known to interfere with mitochondrial atp production. two external alternative nadh dehydrogenases showed a consistent up-regulation in response to freezing temperatures across the time series, while two internal ones from tp to tp. alternative oxidase a gene was up-regulated at tp and while an ortholog of uncoupling protein at tp only. carbon and reductant sources for respiration to determine the source of reductant for mitochondrial electron transport, we assessed changes in gene expression of various genes involved in central carbon metabolic pathways that could provide these substrates. the only potentially active pathways involved sugars. transcript levels of four invertase genes, which are involved in starch degradation, increased, including a cell wall invertase, whose expression level increased through tp. a homolog of at4g a cytosolic invertase, showed large increases through tp the expression of homologs to members of the sucrose synthase gene family showed a somewhat analogous pattern with sus exhibiting a substantial increase through tp sus also showed an increase, albeit a more modest one, which peaked between tp and tp in contrast, sus showed a decrease in tp and no significant changes at other tps. however, sucrose levels remained constant and above baseline levels during freezing acclimation. more generally, genes encoding starch biosynthesis enzymes were moderately up-regulated, while those associated with degradation were down-regulated. adp glucose pyrophosphorylase involved in starch biosynthesis is a heteromer of two large and two small subunits. apl and two of the three genes encoding the large agpase subunits responded and showed increases in the expression from tp to tp a gene encoding the small subunit of agpase showed a decrease in expression at tp and and an increase in tp one gene encoding starch synthase showed an initial decrease in tp and while the other responsive starch synthase gene was down-regulated at tp. most genes encoding glycolytic enzymes did not have considerably altered levels of relevant transcripts in response to freezing. however, some glycolytic genes showed a more substantial up- or down-regulation, including the pyrophosphate-dependent phosphofructokinase that showed a maximal -fold increase and a constant up-regulation. the same step in glycolysis can also be catalyzed by another pyrophosphate-dependent enzyme, fructose-6-p 1-phosphotransferase and the corresponding α - subunit was correspondingly down-regulated. fructose- -bisp aldolase transcript levels increased consistently in response at all tps. in arabidopsis, at5g encodes a protein similar to cytosolic pyruvate kinase, and the corresponding gene in sitka spruce showed a 3-fold increase in transcript levels. transcript levels of gluconeogenetic pyruvate pyrophosphate dikinase decreased at tp and remained low throughout the time course. phosphoenolpyruvatecarboxykinase transcript gradually decreased fold through tp dihydrolipoamide dehydrogenase, the e subunit of pyruvate dehydrogenase, was moderately, but significantly up-regulated from tp to tp. in contrast, lpd, the e subunit of pyruvate dehydrogenase, was down-regulated about -fold between tp and tp the steady-state transcript levels for other subunits of pyruvate dehydrogenase remained unchanged during freezing. this was true for most of the transcripts encoding the remaining enzymes of the tca cycle, though they tended to be slightly above the baseline levels throughout the time course. a putative cytosolic nadp-dependent isocitrate dehydrogenase showed the same pattern and degree of increases in gene expression as lpd a putative regulatory subunit of the mitochondrial nad-dependent isocitrate dehydrogenase iii increased 2-fold in tp and showed a stable -fold increase in transcript levels in tp through tp malic enzyme gene showed only a moderate -fold up-regulation in response to freezing, but the consistent 9- to 11-fold induction of the at2g13560-encoded nad-dependent malic enzyme transcript levels in tp and was quite striking. citrate and succinate levels were moderately increased in a stable manner, while malate levels gradually increased up to fold during freezing. discussion and  temperate and boreal conifers must survive harsh winters while retaining their photosynthetic tissue, and have evolved intricate mechanisms to maintain a certain level of active metabolism, signaling, and protective processes at low temperatures in their needles. the comprehensiveness of our data afforded a global view of winter hardening as it is reflected in numerous processes. the use of bioinformatics tools enabled the visualization of these changes in a cellular context. the results point to numerous photoprotective processes that are active during winter, serve to confirm previous reports, and provide new insights into the activity of genes associated with chloroplast function, energy metabolism, antioxidant processes, and possible regulatory mechanisms. genes related to energy quenching mechanisms including a spruce psbs homolog and plastid terminal oxidase ptox increased their expression, which is in an agreement with previous reports. psbs activity is associated with changes in photosystem ii that occur as part of the photoprotective mechanisms that are set in place as the temperature drops, while ptox represents an important component of the dissipation of excess light energy in photosynthesis as an alternative electron acceptor and o2-utilizing enzyme helping to prevent ros formation in hardening conifer needles. in ranunculus glacialis, an alpine plant growing in high altitudes at high light and low temperatures, ptox protein levels correlated with the use of alternative etc. and enabled high fluxes through photosynthetic etc. despite the excess of light energy. ptox was implicated to serve as a safety valve for the over-reduction of photosystems during high light and low temperatures in alpine plants. ptox has also been implicated in increasing the flux through the linear etc. and subsequent increased tolerance of photosystems to high light in cold-hardened arabidopsis thaliana. this result suggests that the energy dissipation mechanisms set in motion within the chloroplast during winter hardening were effective. one manifestation of this may have been a relatively low level of ros generation as the temperature dropped. however, redox effects on cellular responses to the environment are not confined to a simple containment of ros. much data point to regulatory roles for specific antioxidant proteins in plant cells under stress. several of these antioxidant, regulatory genes were up-regulated over the winter hardening time course, including catalase cat glutaredoxin roxy and glutathione reductase gr catalases regulate the levels of hydrogen peroxide, which is an important signaling molecule in plant stress responses. roxy1encodes a glutaredoxin with a documented regulatory role in development in arabidopsis –, while gr1exerts a regulatory effect on responses to cold stress, an effect that cannot be replaced by gr or thioredoxins. the ros scavenging pathway may have been operational in the chloroplast as winter hardening proceeded, since other components, such as mondehydroascorbate reductase, were also up-regulated. however, it is also possible that tapx was acting in a different/additional capacity to protect the photosynthetic machinery. duan et al. have shown that tapx plays a specific role in alleviating the photosystem i and ii-associated photoinhibition that occurs with chilling stress in arabidopsis, and spruce tapx may play an analogous role in the chloroplasts of hardening needles. interestingly, maruta et al. provide evidence for a specific regulatory role for tapx in chloroplast responses to stress, resonant with the role shown for gr. in the case of tapx and the winter hardening data, the guilt by association analysis revealed an association of tapx exclusively with plastid-targeted gene products, in which the great majority of those responsive genes were down-regulated. it is possible that tapx plays an important role in the adaptation mechanism that results in the onset of chloroplast dormancy in conifer needles. taken together, these results suggest specific regulatory roles for tapx, roxy and gr in the adaptation mechanism underlying winter hardening. freezing temperatures impose severe stress on the mitochondrial etc. a set of genes encoding components of the respiratory etc. were strongly or mildly up-regulated, suggesting that the mitochondrial etc. is active. mitochondrial aerobic respiration is a source of ros formed through the reduction of molecular oxygen by electrons coming from complexes i and iii of the etc.. the ros production can be prevented by redirecting respiratory electron flux to other substrates by involving the major players of the alternate etc. including alternative oxidases, external and internal alternative nadh dehydrogenases, and uncoupling proteins –. genes encoding all these components of the alternate etc. were up-regulated during the time course, suggesting that these alternate pathways are active, along with the conventional etc., during freezing acclimation. up-regulation of genes encoding the components of the alternate etc. is typical during cold acclimation in plants in general and is accompanied with metabolic changes. tobacco plants subjected to cold showed an increase in alternative oxidase levels and sugar content compared to the transgenic plants with suppressed expression of this enzyme. these transgenic plants also showed an increased lipid peroxidation at normal, but not at low, temperatures. alternative oxidases are regulated at the transcriptional level by abscisic acid insensitive transcription factor and are activated by pyruvate –. accumulation of pyruvate as the final product of glycolysis indicates an imbalance between glycolytic fluxes and the mitochondrial etc. pyruvate then activates alternative oxidases to stimulate flux through the mitochondrial etc. to keep pace with glycolysis. suppressing the external nadh dehydrogenase ndb resulted in increased levels of ndb and alternative oxidase 1a in tobacco, which lead to increased salt tolerance and improved growth. glycolysis and the tca cycle as sources of carbon, atp, and reductant avoiding ros production comes at a price, as every alternate route leads to a decreased atp production by mitochondrial atp synthase. both normal and alternate etcs were active based on our data, which brings a question of the potential source of reductant for these etcs and source of atp that would supplement atp lost due to the potentially high activities of the alternate routes for electrons. most glycolytic genes did not show major steady-state transcript differences during freezing acclimation, though persistent up-regulation was observed, suggesting that glycolysis was active. in plants, glycolysis is regulated at the transcriptional and allosteric levels and in a complex way due to metabolic flexibility. in situations when atp cannot be easily produced in heterotrophic systems, pyrophosphate- instead of atp-utilizing enzyme forms are used to drive glycolysis. one of the genes encoding phosphofructokinase was up-regulated up to -fold during winter hardening. up-regulation of a pyrophosphate-dependent phosphofructokinase isoform suggests the possibility of conservation of atp during freezing acclimation, similar to phosphate starvation stress. to maximize atp production, alternate glycolytic steps involving substrate-level phosphorylation are preferentially utilized. the conversion of glyceraldehyde-3p to 3-phosphoglycerate can be done in two ways by: non-phosphorylating glyceraldehyde-3p dehydrogenase in a single step to only produce nadph and phosphoglycerate kinase and phosphorylating glyceraldehyde-3p dehydrogenase in two steps involving nadh and atp production. the phosphorylating glyceraldehyde-3p dehydrogenase was moderately, but consistently up-regulated, while the kinase showed a consistent, less than -fold down-regulation, suggesting that the two-step reaction producing both the reductant and atp could be active. atp production by substrate-level phosphorylation appears to be active in general – three transcripts coding for pyruvate kinase showed varying, but consistent up-regulation. pyruvate kinase catalyzes a key step regulating glycolysis and pyruvate is a known allosteric regulator of a number of glycolytic enzymes including pyruvate kinase, activity of which is induced during hypoxia when glycolysis is the major source of atp. based on these observations, it appears that glycolytic routes conserving or producing atp are preferentially used. this is in an agreement with the global activation of alternative etc. routes. as such, glycolysis could be a source of atp supplementing lower atp levels obtained from respiration and could help with high-energy cofactor balancing. in conventional glycolysis, pyruvate kinase converts phosphoenolpyruvate to pyruvate and atp is made. however, mitochondria could oxidize substrates other than pyruvate, such as malate, during winter hardening. in this scenario, the pyruvate kinase step is bypassed by three consecutive steps with no net nadh and atp production or consumption. however, this bypass allows nadh from the cytosol to be shuttled inside mitochondria where it can be used by the conventional or alternative etcs. first, cytosolic phosphoenolpyruvate is carboxylated to oxaloacetate by phosphoenolpyruvate carboxylase, a highly regulated enzyme whose activity increases during phosphate starvation –. second, oxaloacetate is reduced to malate by cytosolic nadh-utilizing malate dehydrogenase, and malate is transported to mitochondria. the dicarboxylate/tricarboxylate transporter dtc encoded by the arabidopsis at5g gene ortholog that is capable of transporting malate to mitochondria showed a consistent 5-fold up-regulation. third, malate either enters the tca cycle or is decarboxylated to pyruvate by nad-dependent malic enzyme in mitochondria. most genes encoding the enzymes of the tca cycle showed a moderate up-regulation or no significant changes in their steady-state transcript levels, suggesting that the tca cycle is operational. for example, a putative regulatory subunit of the mitochondrial nad-dependent isocitrate dehydrogenase iii showed a stable -fold increase in transcript levels in tp through tp malic enzyme has been implicated in generating nadh and pyruvate needed for lipid synthesis in heterotrophic and photoheterotrophic systems such as plant seeds –. the physiological significance of this bypass in most in vivo situations is not well understood. however, existing metabolomic and transcriptomic data suggest that the bypass is operational during freezing acclimation. malate levels increased significantly at tp and reached the levels that were nearly 3-fold higher than baseline later in the time course. two genes encoding phosphoenolpyruvate carboxylase and one gene encoding malic enzyme were moderately up-regulated in response to freezing acclimation, while cytosolic malate dehydrogenase transcript levels were not significantly different in all tps. however, the sitka spruce homolog of the arabidopsis at2g13560-encoded nad-dependent malic enzyme which is allosterically activated by fumarate, showed a 9- to 11-fold up-regulation in its transcript levels during freezing acclimation and could drive significant fluxes through this bypass. low temperature stress was also shown to induce the expression of a gene encoding malic enzyme in maize. to drive efficient conventional and alternative etcs during winter hardening, this bypass could: shuttle nadh from cytosol to mitochondria, bring additional carbon from glycolysis to the tca cycle to stimulate alternative and conventional mitochondrial etcs, and generate phosphate for atp synthesis. potential sources of hexoses for glycolysis the reaction catalyzed by invertase requires atp molecules to generate fructose-6p, while the net ntp consumption for the sucrose synthase pathway is zero and the latter is up-regulated during hypoxia to save atp. both invertases and sucrose synthases showed significant increases in the corresponding transcript levels in sitka spruce needles during freezing acclimation. however, steady-state sucrose levels did not change significantly. so, if there were an increase in sucrose degradation to provide glucose and fructose for glycolysis, there would be a corresponding increase in sucrose supply to maintain unchanging sucrose levels. because photosynthetic carbon fixation is close to inactive, one can hypothesize that sucrose would have to be transported from organs that have a great potential to store it or made from other, more condensed form of carbohydrate such as starch that can easily be converted to sucrose. starch can accumulate transiently in photosynthetically active cells or be stored for a prolonged period of time in tubers and seeds –. this makes starch an excellent candidate as a storage compound that could be mobilized and used as a source of hexoses during winter hardening. however, the majority of genes associated with starch degradation were down-regulated, while those involved in its biosynthesis were up-regulated, which would make starch present in needles a seemingly unlikely storage form of carbon. regulation of starch metabolism is complex, as it involves several levels of regulation, including transcriptional, protein phosphorylation, allosteric, and redox regulation, to contribute to the balance between photosynthate production and demand in source and sink tissues – and as such, transcript levels alone cannot provide enough information regarding the direction of starch metabolism during winter hardening. changes in starch levels were not measured during freezing acclimation in sitka spruce needles, but fisher &holl showed that starch content progressively decreased in scots pine needles during the winter, suggesting that this polysaccharide could possibly serve as a storage compound for the winter hardening processes also in sitka spruce needles despite the up-regulation in starch synthesis genes. however, an increase in starch synthesis was also observed in arabidopsis leaves exposed to prolonged moderate osmotic stress that lead to similar global gene expression responses involving cell wall remodeling, vesicular transport, and hormonal signaling as in winter hardening, but the reason for these increases in starch synthesis remains to be elucidated. addressing the question of the storage form of carbon in needles during winter hardening is quite complicated. starch is a potential candidate, but conifers also accumulate di and triacylglycerols as storage compounds in the needles. unfortunately these compounds were not measured in sitka spruce needles during winter hardening, but both types of lipids showed a similar trend as their levels gradually increased during fall and winter in scots pine needles, followed by a rapid decrease in spring prior to bud break, suggesting that lipid degradation provides additional carbon and energy for new buds before photosynthesis is fully operational. it is unclear whether these lipids can also be used as a source of energy during winter hardening. starch levels decreasing and lipid levels increasing in needles during the winter is consistent with the hypothesis that starch rather than lipids are substrates for respiration and active cell-wall remodeling during winter. in addition, degradation of oligo and polysaccharides is more energy effective for producing hexoses for substantial cell-wall remodeling than potential gluconeogenesis from lipids. a possibility also exists that soluble and insoluble sugars are stored in other parts of the trees and then transported in the form of sucrose to needles to help maintaining active metabolism during prolonged freezing. all these compounds including di and triacylglycerols are present in the sapwood of scots pine, but detailed source-sink relationship studies are needed to confirm the possibility of sucrose transport to needles in winter. 
cold acclimation in woody perennials is a metabolically intensive process, but coincides with environmental conditions that are not conducive to the generation of energy through photosynthesis. the mapman tool revealed enrichment of the expression of genes related to mitochondrial function, antioxidant and associated regulatory activity, while changes in metabolite levels over the time course were consistent with the gene expression patterns observed. while the negative effects of low temperatures on the photosynthetic apparatus during winter have been well studied, less is known about how this is reflected at the level of gene and metabolite expression, nor how the plant generates primary metabolites needed for adaptive processes during autumn. genes related to thylakoid function were down-regulated as expected, with the exception of plastid targeted specific antioxidant gene products such as thylakoid-bound ascorbate peroxidase, components of the reactive oxygen species scavenging cycle, and the plastid terminal oxidase. in contrast, the conventional and alternative mitochondrial electron transport chains, the tricarboxylic acid cycle, and redox-associated proteins providing reactive oxygen species scavenging generated by electron transport chains functioning at low temperatures were all active. a regulatory mechanism linking thylakoid-bound ascorbate peroxidase action with “chloroplast dormancy” is proposed.
 cold acclimation in woody perennials is a metabolically intensive process, but coincides with environmental conditions that are not conducive to the generation of energy through photosynthesis. while the negative effects of low temperatures on the photosynthetic apparatus during winter have been well studied, less is known about how this is reflected at the level of gene and metabolite expression, nor how the plant generates primary metabolites needed for adaptive processes during autumn.  the mapman tool revealed enrichment of the expression of genes related to mitochondrial function, antioxidant and associated regulatory activity, while changes in metabolite levels over the time course were consistent with the gene expression patterns observed. genes related to thylakoid function were down-regulated as expected, with the exception of plastid targeted specific antioxidant gene products such as thylakoid-bound ascorbate peroxidase, components of the reactive oxygen species scavenging cycle, and the plastid terminal oxidase. in contrast, the conventional and alternative mitochondrial electron transport chains, the tricarboxylic acid cycle, and redox-associated proteins providing reactive oxygen species scavenging generated by electron transport chains functioning at low temperatures were all active.  a regulatory mechanism linking thylakoid-bound ascorbate peroxidase action with “chloroplast dormancy” is proposed. most importantly, the energy and substrates required for the substantial metabolic remodeling that is a hallmark of freezing acclimation could be provided by heterotrophic metabolism. microarraysitka sprucecarbon metabolismcell wallsadaptation mechanismsvisualization 