- 1 Departamento de Biologia Vegetal, Universidade Federal de Viçosa, Viçosa, Brazil
- 2 Max-Planck Partner Group, Departamento de Biologia Vegetal, Universidade Federal de Viçosa, Viçosa, Brazil
The application of post-genomic techniques in plant respiration studies has greatly improved our ability to assign functions to gene products. In addition it has also revealed previously unappreciated interactions between distal elements of metabolism. Such results have reinforced the need to consider plant respiratory metabolism as part of a complex network and making sense of such interactions will ultimately require the construction of predictive and mechanistic models. Transcriptomics, proteomics, metabolomics, and the quantification of metabolic flux will be of great value in creating such models both by facilitating the annotation of complex gene function, determining their structure and by furnishing the quantitative data required to test them. In this review, we highlight how these experimental approaches have contributed to our current understanding of plant respiratory metabolism and its interplay with associated process (e.g., photosynthesis, photorespiration, and nitrogen metabolism). We also discuss how data from these techniques may be integrated, with the ultimate aim of identifying mechanisms that control and regulate plant respiration and discovering novel gene functions with potential biotechnological implications.
Introduction
Metabolism represents one of the most important and most likely the best characterized network within biological systems. From pioneering studies defining key metabolic pathways, subsequent decades of enzymology characterizing the catalytic and regulatory properties of enzymes, through to more recent genetic studies of metabolism, there is an unprecedented density of both mechanistic and descriptive data relating to metabolic behavior. Advances in the understanding of metabolic regulation and control, however, still suffer from insufficient research concerning the mode of operation of metabolic networks, and plant respiratory metabolism provides an illustrative example of this problem. Due to the crucial roles of mitochondria in meeting the energy demands of the cell whilst simultaneously being involved in amino acid metabolism and photorespiration, precise control of mitochondrial metabolism and function is crucial for cellular homeostasis, and unsurprisingly mitochondrial dysfunction leads to diverse metabolic and phenotypic consequences (Nunes-Nesi et al., 2011). However, despite the fact that the major respiratory pathways in plants were elucidated decades ago (Beevers, 1961), our current knowledge about their regulation and control, and how this relates to whole plant physiology, is still limited (Fernie et al., 2004). Unsurprisingly then, considerable research effort is being devoted to elucidating the metabolic basis of the regulation of the tricarboxylic acid (TCA) cycle and the mitochondrial electron transport chain (Sweetlove et al., 2007, 2010; Fahnenstich et al., 2008; Fukushima et al., 2009; Araújo et al., 2010; Tomaz et al., 2010; Zell et al., 2010) as well as their interactions with photosynthesis (Nunes-Nesi et al., 2007, 2011), photorespiration (Bauwe et al., 2010, 2012), and nitrate assimilation (Foyer et al., 2011).
Given the complex patterns of regulation and the interactions between respiratory metabolism and other processes uncovered by such studies, it is likely that only through the use of model-based approaches we will be able to understand how the individual components of metabolism (e.g., enzymes and metabolites) work together to produce a functioning system, and in this way truly understand the function of individual components of that system. In this context systems biology offers the promise of a better understanding of respiratory metabolism and its role within the metabolic network.
Systems biology aims to understand how populations of molecules, cells, and organisms interact to give rise to complex biological processes, including cell division, growth, development, metabolism, and behavioral and ecological patterns (Karsenti, 2012). Although initially applied exclusively to mathematical modeling strategies (Edwards and Palsson, 1999) systems biology may now incorporate data from the functional genomics tools of transcriptomics, proteomics, and metabolomics as well as genomics (Somerville et al., 2004; Baginsky and Fernie, 2007), and has become a truly interdisciplinary area of research, some basic and applied aspects of which have been expertly reviewed elsewhere (Ideker et al., 2001; Oltvai and Barabási, 2002; Kitano, 2004; Somerville et al., 2004; Breitling, 2010). Whilst an accurate understanding of metabolism as a system requires far more than enumeration of its components (Kleessen et al., 2012), it is nonetheless the case that the development of mathematical models of plant metabolism in general will require accurate gene annotations, functional characterization of enzymes, and the acquisition and statistical analysis of quantitative data.
In the case of plant respiratory metabolism, whilst the core components, including enzymes of the TCA cycle and components of the mitochondrial electron transport chain, have been identified and characterized in numerous species (Millar et al., 2011; Araújo et al., 2012a), how they are regulated in response to changing environmental conditions and developmental programs remains less clear, and extent to which respiratory metabolism interacts with other parts of the metabolic network is only now beginning to emerge (Balmer et al., 2004; Sweetlove et al., 2006; Timm et al., 2008, 2011; Fukushima et al., 2009; Tomaz et al., 2010; Foyer et al., 2011). Moreover, the functions of more peripheral components of plant respiratory metabolism remain to be investigated. Genes for many of the mitochondrial transporters that are hypothesized to exist have yet to be identified, or have only been assigned on the basis of homology with non-plant species (Linka and Weber, 2010; Palmieri et al., 2011), whilst it is only recently that several alternative electron donors to the mitochondrial electron transport chain have been characterized (Ishizaki et al., 2005, 2006; Araújo et al., 2010). Although here we focus mainly on function and interactions of respiratory enzymes, it is important to note that the construction of predictive models of mitochondrial metabolism will also require a better understanding of electron transport processes. Thus, both structural biology and measurements of electron transport rates remain exciting topics for further study.
The availability of comprehensive plant genome information, coupled with the integration of large-scale unbiased molecular profiling technologies such as whole-genome microarrays, quantitative proteomics, and metabolite profiling (Table 1; Figure 1) has increased our ability to annotate gene functions within plant respiratory metabolism, obtain quantitative data and uncover unanticipated relationships between respiration and other cellular processes such as photosynthesis, photorespiration, redox regulation, and signaling (Sweetlove et al., 2002; Heazlewood et al., 2004; Dutilleul et al., 2005; Urbanczyk-Wochniak et al., 2006). This information will ultimately facilitate the application of systems biology to plant respiratory metabolism (Fernie, 2012) and progress in these areas as they relate to research on plant respiration and the annotation of gene function will be discussed in the following sections.
TABLE 1. Omics fields associated with illustrative and available Arabidopsis resources on web.
FIGURE 1. Schematic representation of a synergetic integration of multiple omics approaches. Within individual experiments, a range of data is generated. This data can be associated with different levels of omics and the possible interactions are illustrated. The usage of modeling approaches linking the interactions between the results obtained using different post-genomic techniques is expected. Additionally, the data obtained can be deposited in specific database contributing to the construction of general prediction databases (e.g. coexpression database). This also allows the generation of further testable hypothesis and to adjust such hypothesis for individual experiments.
Applying Transcriptomics to the Study of Mitochondrial Function
Determining the complete transcriptional capability of the cell, including large and small RNAs, novel transcripts from unannotated genes, splicing isoforms, and gene-fusion transcripts, serves as the foundation for a comprehensive study of the transcriptome. Recent advances in DNA sequencing technology have greatly increased both sequencing scale and throughput and led to the development of RNA-Seq (Liu et al., 2011; Martin and Wang, 2011), a technique capable of complementing microarray-based transcriptomics. In recent years transcriptomic techniques have been extensively used to reveal interactions between plant respiratory metabolism and other important processes such as seed germination (Howell et al., 2009), stress tolerance (Giraud et al., 2008; Lehmann et al., 2009; Meyer et al., 2009), hypoxia (Narsai et al., 2009), and the operation of the circadian clock (Fukushima et al., 2009).
The combination of a comprehensive list of putatively mitochondrially expressed genes with a large number of publicly available stress-related microarray datasets was used to define the mitochondrial stress response in Arabidopsis (van Aken et al., 2009). This complex analysis indicated that the mitochondrial stress response extends far beyond the alternative electron transport chain components, with, most prominently, mitochondrial substrate carrier proteins and heat shock proteins also showing extensive stress responsiveness. It is important to note, however, that many of the stress-responsive genes discovered by this approach are of unknown function and this indicates an attractive avenue for research aiming to identify new mitochondrial functions and targets for engineering of stress tolerance.
Transcriptomic studies can also provide details of expression patterns during plant development. Recent work indicated that BCS1, a likely AAA-ATPase, displayed similar levels of expression throughout development (Duncan et al., 2011), despite the fact that previous studies have shown that this gene responds to various stresses (van Aken et al., 2009). Based on these patterns of expression the authors propose that the protein encoded by this gene may be involved in protein repair in response to stress (Xu et al., 2011). This study also indicated that the transcript abundance for genes encoding outer and inner membrane mitochondrial proteins was maximal during seed germination, and identified a subset of genes that were highly expressed in root tissues (Duncan et al., 2011).
Responses to stress conditions, including the effect of absence of alternative oxidase1a (aox1a) expression in Arabidopsis, has also been investigated using transcriptomics. Such plants, when grown in moderate light under drought conditions, displayed 10-fold increases in leaf anthocyanin levels, together with alterations in photosynthetic efficiency, increased superoxide radical production, and reduced root growth (Giraud et al., 2008). Furthermore, microarray and quantitative reverse transcription polymerase chain reaction (qRT-PCR) analysis indicated that even under normal growth conditions genes normally induced under stress conditions were expressed (Giraud et al., 2008), including those involved in defense against reactive oxygen species (ROS) and in stress signaling. This study was incorporated into a recent meta-analysis of transcriptomic studies that aimed to uncover targets of mitochondrial retrograde signaling (Schwarzländer et al., 2012), and construct a model for how mitochondrial dysfunction may affect nuclear gene expression. Transcripts encoding proteins involved in photosynthesis, protein synthesis, and plant–pathogen interactions were revealed as the major targets of retrograde regulation.
A combination of transcript and metabolite profiling methods has been used to investigate the molecular and physiological responses following root hypoxia caused by flooding in gray poplar (Populus × canescens; Kreuzwieser et al., 2009). Interestingly, changes in metabolite levels occurred in both roots and leaves, whilst changes in transcript abundance were restricted to the roots, the actual site of hypoxia. The general pattern of metabolite and transcript abundance suggests that the response to hypoxia comprises both the repression of energy demanding processes such as cell wall biosynthesis, and the reconfiguration of carbohydrate and nitrogen metabolism to ensure sufficient substrate supply during long periods of root hypoxia (Kreuzwieser et al., 2009).
Proteomic Approaches for the Discovery of Gene Function in Plants
Proteomic approaches have been used to generate great insights into plant metabolism in general and into mitochondrial respiratory metabolism in particular. Recently, the utilization of multiple fluorescent dyes in single gels (DIGE) has renewed the use of 2-D gels for quantitative comparisons of the mitochondrial proteome. Studies using DIGE have investigated knockout mutants for carbonic anhydrase-like proteins (Perales et al., 2005), complex I (Meyer et al., 2009), and malate dehydrogenase (Tomaz et al., 2010), as well as differences between tissues (Lee et al., 2008), the impact of rotenone (Garmier et al., 2008), and changes associated with the diurnal cycle (Lee et al., 2010). Shotgun proteomic techniques utilizing Liquid chromatography coupled tandem mass spectrometry (LCMS/MS) of trypsin-digested samples, without the use of gels, have also been used for in-depth studies of the mitochondrial proteome in Arabidopsis (Brugiere et al., 2004; Heazlewood et al., 2004) and rice (Heazlewood et al., 2003b; Huang et al., 2009). These studies greatly increased the number of proteins identified within plant mitochondrial isolates that are below the level of detection by techniques based on gel separation and staining.
The mitochondrial proteome of Arabidopsis in particular has been extensively analyzed by both gel-based and gel-free strategies (Kruft et al., 2001; Millar et al., 2001; Heazlewood et al., 2004; Lee et al., 2008, 2011), and may contain as many as 2000–3000 different proteins, each of which may be subject to post-translational modification (Millar et al., 2005; Taylor et al., 2011). Recent proteomic studies have been able to identify more than 500 proteins, including subunits of mitochondrial respiratory complexes, supercomplexes, phosphorylated proteins, and oxidized proteins (Millar et al., 2005; Klodmann and Braun, 2011). However, even these identified proteins have not all been functionally characterized. Moreover, targeting prediction software tools assign more than 1,500 proteins encoded by the Arabidopsis genome to this subcellular compartment and it therefore seems likely that the function of most mitochondrial proteins, especially those of low abundance and/or high hydrophobicity (Klodmann and Braun, 2011), remains to be discovered.
In addition to analysis of the mitochondrial proteome as a whole, proteomic approaches have been used extensively to dissect mitochondrial protein complexes, in particular complex I of the respiratory electron transport chain (Leterme and Boutry, 1993; Herz et al., 1994; Klodmann and Braun, 2011), and the importance of this characterization using different proteomic approaches has recently been expertly reviewed (Klodmann and Braun, 2011). These methods have revealed that the complex I in plants has an extremely sophisticated configuration and although some of the postulated functions of the characterized extra subunits remain to be fully investigated, proteomic approaches have already contributed significantly to our current understanding of this complex (Heazlewood et al., 2003a; Cardol et al., 2004; Klodmann et al., 2010, 2011). For example, recent work has revealed that plant complex I is especially large and includes 15 extra unique subunits. Five of these subunits represent proteins resembling carbonic anhydrases and one represents a L-galactono-1,4-lactone dehydrogenase (GLDH), introducing side activities to plant complex I, while all the remaining nine subunits are rather small (7–12 kDa) and of so far unknown function (Klodmann et al., 2010; Klodmann and Braun, 2011).
Many mitochondrial proteins in fact form part of protein complexes and it has recently been shown that the formation of supercomplexes, “complexes of complexes,” could have a regulatory function in guiding electrons through alternative respiratory pathways, particularly under variations in oxygen levels (Ramírez-Aguilar et al., 2011). The extent to which this phenomenon occurs is currently unclear though as only the most abundant complexes have so far been functionally characterized.
Proteomics may also help to provide details of the mechanisms that regulate mitochondrial metabolism, and recent studies have identified potential thioredoxin-linked proteins in mitochondria isolated from autotrophic and heterotrophic plant tissues (Balmer et al., 2004; Martí et al., 2009). Fifty mitochondrial thioredoxin-linked proteins were identified including six TCA cycle enzymes [aconitase, succinyl-CoA ligase, isocitrate, malate, pyruvate, and succinate dehydrogenases (SDH); Balmer et al., 2004]. These findings suggest mitochondrial enzymes as potential targets for redox regulation through thioredoxins. However, despite the fact that these studies, amongst others, indicate that by sensing redox state thioredoxins enable mitochondria to communicate with other organelles in photosynthetic tissues, considerable experimental effort is still required to provide functional details.
Insights Into plant Respiratory Metabolism Obtained from Metabolite Profiling
The term metabolome can be defined as the total small-molecule complement of a cell, and metabolomics is therefore the study of all small molecules or metabolites presents in a cell or organism (Oliver et al., 1998; Tweeddale et al., 1998). Although no single analytical system is ever likely to cover the whole metabolome technological developments have considerably extended our ability to analyze complex biological systems, facilitating the simultaneous detection of different compound classes with diverse chemical properties (Osorio et al., 2012). Gas chromatography coupled to mass spectrometry (GC-MS) is a versatile and widely applied technique in modern metabolomic studies and allows the identification and quantification of a relatively broad range of compounds, including organic and amino acids, sugars, sugar alcohols, phosphorylated intermediates, and lipophilic compounds (Fiehn et al., 2000; Roessner et al., 2001; Weckwerth, 2003; Lisec et al., 2006; Schauer et al., 2006), resulting in fairly comprehensive coverage of the central pathways of primary metabolism. Metabolite profiling using LC-MS and capillary electrophoresis-MS can complement the analysis by GC-MS for certain classes of metabolite. Of particular note, tandem MS methods for the quantification of phosphorylated sugars (Lunn et al., 2006; Luo et al., 2007; Hasunuma et al., 2010), Calvin–Benson cycle intermediates (Arrivault et al., 2009) and cell wall precursors (Alonso et al., 2010), which are difficult to analyze by GC-MS, have recently been developed.
Due in part to these technological advances application of metabolite profiling to plant biology is becoming increasingly common, including diagnostic and descriptive analyses of metabolic responses to genetic and/or environmental perturbations, the use in annotation of gene function and in systems biology (Schauer and Fernie, 2006; Guy et al., 2008; Fernie and Schauer, 2009). Recently, metabolite profiling has, in particular, been used to generate significant insight into plant respiratory metabolism and its interactions with photosynthesis. Despite recognition of the importance of mitochondrial metabolism our understanding of the complex pathways through which organic acids are metabolized as well as how these pathways are regulated in vivo remains far from complete (Fernie et al., 2004; Sweetlove et al., 2007; Fait et al., 2008; Møller and Sweetlove, 2010). As part of attempts to remedy this issue, the operation and role of the TCA cycle in illuminated leaves has received particular interest (Krömer, 1995; Raghavendra and Padmasree, 2003; Noguchi and Yoshida, 2008; Nunes-Nesi et al., 2011). The precise mode of operation of the TCA cycle in the light remains somewhat obscure, in part because whilst TCA cycle flux in the light is reduced through inhibition of pyruvate dehydrogenase carbon skeletons deriving from TCA cycle intermediates are still required for nitrogen assimilation. In attempt to address this apparent contradiction, a process of systematic suppression of TCA cycle enzymes through reverse genetics was undertaken. Metabolite profiling of the transgenic lines generated through this process indicated the importance of the TCA cycle in metabolism in illuminated leaves (Nunes-Nesi et al., 2008, 2011). However, these experiments have also revealed a surprising complexity in the response to these manipulations, with suppression of some enzymes leading to increased photosynthesis (Carrari et al., 2003; Nunes-Nesi et al., 2005; Araújo et al., 2011b), others leading to decreased photosynthesis (Nunes-Nesi et al., 2007) and yet others having no effect (Studart-Guimarães et al., 2007; Sienkiewicz-Porzucek et al., 2008, 2010; Sulpice et al., 2010a; Araújo et al., 2012b). Several of the more recent studies from this investigation, highlighting the importance of information obtained from metabolite profiling experiments, are discussed here.
Antisense inhibition of the 2-oxoglutarate dehydrogenase (OGDH) complex in tomato led to alterations in whole plant development that were linked to reductions in total amino acid and nitrate pools despite unaltered photosynthesis (Araújo et al., 2012b). These results clearly imply that OGDH plays an important role in both the TCA cycle and nitrogen assimilation and suggest a novel role for this enzyme in whole plant development and fruit maturation. Furthermore, as was reported for succinyl-CoA ligase antisense plants (Studart-Guimarães et al., 2007), the OGDH antisense plants displayed increased GABA shunt flux, presumably in compensation for decreased succinate production via the TCA cycle (Araújo et al., 2012b). However, it is important to stress that such a compensatory response also caused significant shifts in cellular pools of amino acids and nitrate that were detected by metabolite profiling.
Antisense inhibition of the iron–sulfur subunit of the SDH in tomato plants, on the other hand, resulted in increased photosynthesis and whole plant biomass via an organic acid-mediated effect on stomatal aperture. These results contrast with those obtained for antisense inhibition of fumarase, where decreased photosynthesis and biomass were observed (Nunes-Nesi et al., 2007; Araújo et al., 2011b). Furthermore, measurement of apoplastic organic acid levels in SDH and fumarase antisense plants, revealed a negative correlation between the levels of fumarate and stomatal conductance, though the influence of fumarate appears to be weaker than that of malate (Araújo et al., 2011b). These results provided strong evidence to support that modulation of malate and fumarate concentration can greatly influence stomatal function. Thus, metabolite profiling has aided in the identification of novel interactions between the TCA cycle and photosynthesis, stomatal function and nitrogen metabolism. Further work is clearly required in order to establish the regulatory mechanisms involved in such responses.
Protein degradation during plant development and substrate deficiency can be an important source of substrate for respiratory metabolism, however, our current understanding of the regulation of the classical and alternative pathways of respiration under these conditions is still limited. The combination of genetic approaches with GC-MS-based metabolite profiling has greatly improved our understanding of the complex metabolic interactions observed during dark-induced senescence (Ishizaki et al., 2005, 2006; Araújo et al., 2010, 2011a). These studies have demonstrated that during dark-induced senescence there is a significant accumulation of amino acids and TCA cycle intermediates. Moreover, it has also been demonstrated that both isovaleryl-CoA dehydrogenase and D-2-hydroxyglutarate dehydrogenase provide electrons to the plant ubiquinol pool via the electron transfer flavoprotein (ETF)-ETF:ubiquinone oxidoreductase (ETF/ETFQO) complex (Engqvist et al., 2009; Araújo et al., 2010). Given that the chlorophyll breakdown intermediate phytanoyl-CoA accumulates dramatically both in knockout mutants of the ETF/ETFQO complex and of isovaleryl-CoA dehydrogenase following growth in extended dark periods it was suggested that chlorophyll breakdown could be important for the supply of carbon and electrons during this process. However, metabolic analyses of phytanoyl-CoA 2-hydroxylase knockout mutants under the same extended darkness regime as previously used suggest that phytol and phytanoyl-CoA enzyme does not primarily function as a substrate of the ETF/ETFQO pathway (Araújo et al., 2011a). These studies also showed that these mutants were not compromised in their ability to withstand significant extension of the dark period but do accumulate phytanoyl-CoA and to a lesser extent 2-hydroxyglutarate as well as sharing some of the other metabolic features of mutants of the ETF/ETFQO complex following dark-induced senescence treatment. In summary, the results obtained through this work indicated that both isovaleryl-CoA dehydrogenase and 2-hydroxyglutarate dehydrogenase essentially account for the entire electron input via the ETF complex.
Plant metabolism is also reorganized under a range of different stress conditions including salt, cold, drought, and oxidative stress (Kaplan et al., 2004, 2007; Gong et al., 2005; Sanchez et al., 2008, 2012; Alet et al., 2012; Siahpoosh et al., 2012), allowing plants to continue to produce indispensable metabolites whilst preventing the accumulation of ROS. Metabolite profiling has proven a powerful tool to gain an overview of such reorganizations in response to stressful conditions (Shulaev et al., 2008). Menadione, a quinone which causes ROS generation from both mitochondrial and plastidal electron transport chains, has been used to study metabolic and transcriptomic responses to ROS in Arabidopsis roots (Lehmann et al., 2009). Detailed metabolic studies revealed a down-regulation of glycolysis and TCA cycle and the redirection of carbon from glycolysis to the oxidative pentose phosphate pathway in menadione-treated roots suggesting a reorganization of central carbon metabolism under oxidative stress conditions that was undone rapidly after the removal of menadione (Lehmann et al., 2012).
Contributions of Metabolic Flux Analysis to our Understanding of Plant Respiration
Metabolic flux ultimately underpins both plant growth and development and for this reason experiments that provide insight into metabolic flux, as well as its regulation and control, can be valuable tools to both improve our understanding of biological systems and aid in the discovery of gene function. Feeding experiments using both 13C- and 14C-labeled precursors have been used extensively to delineate metabolic pathways, but may also be used to investigate and quantify metabolic flux. Labeling experiments with 14C offer high sensitivity and fractionation of labeled metabolites and biomass components can readily indicate the fate of metabolized radiolabel. Accordingly, feeding experiments with 14CO2 demonstrated the effect of reduced malate dehydrogenase activity on photosynthetic carbon assimilation in tomato leaves, and revealed that ascorbate feeding led to increased photosynthesis and altered assimilate partitioning in these transgenic lines (Nunes-Nesi et al., 2005). Quantifying the release of 14CO2 from the metabolism of different isotopomers of 14C-glucose can also provide information about the relative activity of different primary metabolic pathways (Nunes-Nesi et al., 2005). This method has been used extensively to identify changes in metabolic flux in plants with altered TCA cycle enzyme activities; for example, decreased carbon entry into the TCA cycle was identified in tomato plants with decreased SDH and OGDH activity (Araújo et al., 2011b, 2012b), which are in good agreement with measurements of respiration and the results of metabolite profiling experiments.
Labeling experiments with stable isotopes can be used for similar ends, however since detection of label incorporation is in these cases carried out by mass spectrometry or Nuclear Magnetic Resonance (NMR) spectroscopy it is possible to determine both enrichment and location of incorporation of label in individual metabolites without laborious chemical cleavage. This method was used to investigate the fate of 13C-labeled glutamate in leaves of tomato plants with decreased succinyl-CoA ligase, and revealed that TCA cycle flux was maintained in these plants by diversion of carbon through the GABA shunt (Studart-Guimarães et al., 2007). Additionally, in vivo NMR spectroscopy can be used to monitor such labeling experiments in real-time. For example, tracing the metabolism of 13C-pyruvate by mitochondria isolated from Arabidopsis plants with reduced manganese superoxide dismutase activity revealed decreased TCA cycle flux (Morgan et al., 2008), most likely as a result of decreased TCA cycle enzyme activity caused by oxidative damage.
Whilst the experiments described above may provide valuable qualitative information about metabolic flux, the number of fluxes that can be quantified using these approaches is typically quite limited. Steady state metabolic flux analysis (Wiechert, 2001; Ratcliffe and Shachar-Hill, 2006) on the other hand allows the absolute quantification of metabolic fluxes in medium size networks, permitting the degree to which a particular environmental condition or gene product controls or regulates multiple fluxes to be determined. In this way, steady state flux analysis has been able to provide important insights into the relationship between photosynthesis and respiration in developing seeds. Labeling experiments and metabolic modeling indicated that developing green embryos are able to decrease losses of fixed carbon during oil synthesis by both refixing CO2 released by respiration (Schwender et al., 2006), and reducing the need for TCA cycle flux through the use of photosynthesis to meet energy demands (Allen et al., 2009). Although these experiments effectively revealed a new function for RuBisCO, steady state analysis can also be used to reveal the effect of alterations in gene expression. This method has been widely applied for the discovery of gene function in microorganisms (Blank et al., 2005; Fischer and Sauer, 2005), and its feasibility in plants was demonstrated by experiments using Arabidopsis embryos (Lonien and Schwender, 2009) that revealed how flux is rerouted in embryos deficient in two plastidic pyruvate kinase isoforms.
The lack of wider adoption of steady state flux analysis as a tool in plant systems biology is in part due to the strict requirements for an experimental system that can obtain metabolic and isotopic steady state, which limits such experiments to tissue culture and isolated organs. Therefore, the so called instationary flux analysis (Wiechert and Nöh, 2005; Young et al., 2008) has the potential to allow quantification of metabolic flux in networks that do not reach isotopic steady state, by incorporating measurements of labeling dynamics and metabolite pool sizes into the modeling process. Recent work has demonstrated the feasibility of this approach for flux quantification in photoautotrophic organisms (Young et al., 2011), whilst the necessary technological approaches for carrying out this work in plants have already been developed (Arrivault et al., 2009; Hasunuma et al., 2010), suggesting that this method may be used in the discovery of gene function in higher plants in the relatively near future.
Whilst empirical determination of fluxes using labeling experiments can generate important biological insight, constraints based modeling (Price et al., 2003) can also be used to investigate the behavior of metabolic networks. The construction of genome scale metabolic models (Durot et al., 2009; Fell et al., 2010) and their analysis using flux balance analysis, provides both a way to exploit the information pertaining to metabolism contained within a sequenced genome and models now exist for several plant species including Arabidopsis (Poolman et al., 2009; Mintz-Oron et al., 2012), maize (Saha et al., 2011), sugarcane (de Oliveira Dal’Molin et al., 2010), and Brassica napus (Hay and Schwender, 2011). Recent work has highlighted the degree to which flux quantification using constraints based modeling agrees with empirical results (Williams et al., 2010; Beurton-Aimar et al., 2011; Chen et al., 2011), but perhaps a more important role of genome scale modeling is as a means to generate hypotheses regarding the capacities of large metabolic networks.
In this regard, flux balance analysis has already provided interesting insights into plant respiratory metabolism. For example, work using a model of Arabidopsis carbon metabolism revealed the large demand for energy made by cell maintenance processes (Poolman et al., 2009), whilst a study of metabolism in barley seeds (Grafahrend-Belau et al., 2009) indicated how fluxes in central carbon metabolism and the TCA cycle may be affected by the hypoxic conditions that can occur within the developing seed. Such models also provide an excellent means for the discovery of new functions for previously annotated genes though the prediction of the consequences of gene knockouts in silico. Such methods are expected to form an important part of methods for metabolic engineering in species with biotechnological importance. Indeed, models have already been created for the model algal species Chlamydomonas reinhardtii (Boyle and Morgan, 2009; Chang et al., 2011; de Oliveira Dal’Molin et al., 2011) that can be further used to explore the potential of microalgae for biofuel production.
Integration of Profiling Data: from Diagnostics to Systems Biology
Whilst investigation of the transcriptome, proteome, and metabolome alone may be highly informative, cellular physiology and the regulation of plant respiratory metabolism is the result of complex interactions between transcripts, proteins, and metabolites. Bearing that in mind straightforward relationships between these components of the metabolic network are therefore not to be expected. Thus, despite the apparent easy understanding of simple relationships between changes in transcripts, proteins, metabolite, and downstream biological functions, quite often there is a major discordance. This apparent discrepancy can often be explained by quantitative analysis of the data sets and information about translational regulation and protein turnover. These incoherencies are not only due to technical issues but also related to the complexity and structure of metabolic networks. Accordingly, some priorities in metabolic research to reduce these problems have recently been proposed (Fernie and Stitt, 2012). For this reason recent research into plant metabolism has often included the usage of modeling and correlation based approaches, and interactions between the results obtained using different post-genomic techniques (Figure 1 and Table 1) are expected.
Such studies have the ability to reveal previously unappreciated and complex interactions between transcripts, proteins, metabolites, and fluxes. In this way, a recent study was able to uncover a potential signaling role for leucine; a combination of extensive metabolite profiling and transcript profiling demonstrated a correlation between expression of a subset of transcripts and the levels of this amino acid (Hannah et al., 2010). An ability to simultaneously quantify components of the metabolic network can also provide information about the relationship between different components of central carbon metabolism and resource allocation. For example, by exploiting the natural variation present in different accessions of Arabidopsis (Sulpice et al., 2009, 2010b), it was possible to detect correlations between enzyme activities, metabolite abundances, and biomass. These studies indicate the presence of a strong negative correlation between starch and biomass, whilst the relative proportion of protein invested in enzymes is positively correlated to biomass accumulation. The relationship between transcripts, metabolites, and fluxes in plant respiratory metabolism was investigated in Arabidopsis cell suspension cultures subjected to menadione treatment (Baxter et al., 2007). The combinatory analysis of these parameters identified a transient reconfiguration of metabolism, involving down-regulation of the TCA cycle and amino acid biosynthesis, apparently working to avoid wasting energy under conditions of oxidative stress. Such methods may also have biotechnological applications. For example, metabolite:transcript correlations identified in potato tubers have revealed target genes for metabolic engineering of sucrose metabolism (Urbanczyk-Wochniak et al., 2003).
Whilst the unbiased nature of such approaches makes them excellent discovery tools, their capacity to improve our understanding of the regulation and control of respiratory metabolism could be further improved by the use of a mechanistic framework that can relate the components of respiratory metabolism and their properties to the functions of this system. In this regard, a strategy recently employed in yeast linking transcriptional regulation to metabolic fluxes (Moxley et al., 2009) may facilitate the integration of data from multiple post-genomic platforms to improve our understanding of plant respiratory metabolism. Whilst direct mRNA measurements were only poorly correlated with metabolic flux, the incorporation of a parameter representing regulation of enzyme activity by metabolite abundance greatly improved this correlation and the resulting model could be used successfully to predict flux changes occurring in amino acid metabolism in yeast. Genome scale metabolic models also provide a promising means via which such integration can take place (Durot et al., 2009; Poolman et al., 2009; Fell et al., 2010; Mintz-Oron et al., 2012). Since such models are based on annotated genome information, transcript profiling, and proteomic data which are already linked to specific genes can potentially be incorporated and used to constrain the behavior of the network. Such an approach has been shown to be useful in microorganisms (Covert and Palsson, 2002) and might also be successfully applied in plants.
Integration may also take place through the use of kinetic models of metabolism (Schallau and Junker, 2010; Rohwer, 2012). Such models are inevitably on a smaller scale than genome scale models, but compensate by providing a mechanistic link between the kinetic and regulatory properties of enzymes and the behavior of the system under study. This approach was used to study metabolism of the aspartate family of amino acids in Arabidopsis (Curien et al., 2009), and revealed how the allosteric properties of enzymes in this pathway permit fluxes in different branches to vary independently of one another. An additional advantage of such models is that once constructed they immediately allow the application of metabolic control analysis, and this has been exploited to investigate the distribution of control over sucrose cycling amongst the enzymes of sucrose synthesis in sugarcane (Rohwer and Botha, 2001; Uys et al., 2007). Whilst the mitochondrial TCA cycle in plants has been the subject of structural modeling (Steuer et al., 2007) and included in models of wider scope (Leduc et al., 2006), a detailed kinetic model based on experimental data has yet to be produced. Profiling techniques could aid in the provision of data for the construction of such a model, but more importantly, a kinetic model of the plant TCA cycle would provide a medium through which the results of profiling experiments could be better understood.
Overall, whilst post-genomic technologies have already aided our understanding of respiratory metabolism at a systems level through additions to the “parts list” of the metabolic network (Sweetlove and Fernie, 2005), they are poised to become more useful still through the discovery of new interactions and the provision of the quantitative data required for the construction of predictive models of metabolism. The further development and combination of many analytical techniques (Figure 1) will additionally allow a fuller description of the metabolic status of a plant. When this is achieved, global analyses of RNA, protein, and metabolites will allow us to obtain a full picture of the complexity of the system under study. Although the application of these techniques requires substantial financial investment, it is likely to bring returns in the form of an improved ability to carry out rational engineering of the plant metabolic network.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
Financial support was provided by the National Council for Scientific and Technological Development CNPq-Brazil (grant number 478261/2010-1 to Adriano Nunes-Nesi and 472787/2011-0 to Wagner L. Araújo). Thomas C. R. Williams was supported by a Petrobras (Petroleo Brasileiro S.A.) post-doctoral grant.
References
Alet, A. I., Sanchez, D. H., Cuevas, J. C., Marina, M., Carrasco, P., Altabella, T., Tiburcio, A. F., and Ruiz, O. A. (2012). New insights into the role of spermine in Arabidopsis thaliana under long-term salt stress. Plant Sci. 182, 94–100.
Allen, D. K., Ohlrogge, J. B., and Shachar-Hill, Y. (2009). The role of light in soybean seed filling metabolism. Plant J. 58, 220–234.
Alonso, A. P., Piasecki, R. J., Wang, Y., Laclair, R. W., and Shachar-Hill, Y. (2010). Quantifying the labeling and the levels of plant cell wall precursors using ion chromatography tandem mass spectrometry. Plant Physiol. 153, 915–924.
Araújo, W. L., Ishizaki, K., Nunes-Nesi, A., Larson, T. R., Tohge, T., Krahnert, I., Witt, S., Obata, T., Schauer, N., Graham, I. A., Leaver, C. J., and Fernie, A. R. (2010). Identification of the 2-hydroxyglutarate and isovaleryl-CoA dehydrogenases as alternative electron donors linking lysine catabolism to the electron transport chain of Arabidopsis mitochondria. Plant Cell 22, 1549–1563.
Araújo, W. L., Ishizaki, K., Nunes-Nesi, A., Tohge, T., Larson, T. R., Krahnert, I., Balbo, I., Witt, S., Doermann, P., Graham, I. A., Leaver, C. J., and Fernie, A. R. (2011a). Analysis of a range of catabolic mutants provides evidence that phytanoyl-coenzyme A does not act as a substrate of the electron-transfer flavoprotein/electron-transfer flavoprotein:ubiquinone oxidoreductase complex in Arabidopsis during dark-induced senescence. Plant Physiol. 157, 55–69.
Araújo, W. L., Nunes-Nesi, A., Osorio, S., Usadel, B., Fuentes, D., Nagy, R., Balbo, I., Lehmann, M., Studart-Witkowski, C., Tohge, T., Martinoia, E., Jordana, X., Damatta, F. M., and Fernie, A. R. (2011b). Antisense inhibition of the iron-sulphur subunit of succinate dehydrogenase enhances photosynthesis and growth in tomato via an organic acid-mediated effect on stomatal aperture. Plant Cell 23, 600–627.
Araújo, W. L., Nunes-Nesi, A., Nikoloski, Z., Sweetlove, L. J., and Fernie, A. R. (2012a). Metabolic control and regulation of the tricarboxylic acid cycle in photosynthetic and heterotrophic plant tissues. Plant Cell Environ. 35, 1–21.
Araújo, W. L., Tohge, T., Osorio, S., Lohse, M., Balbo, I., Krahnert, I., Sienkiewicz-Porzucek, A., Usadel, B., Nunes-Nesi, A., and Fernie, A. R. (2012b). Antisense inhibition of the 2-oxoglutarate dehydrogenase complex in tomato demonstrates its importance for plant respiration and during leaf senescence and fruit maturation. Plant Cell 24, 2328–2351.
Arrivault, S., Guenther, M., Ivakov, A., Feil, R., Vosloh, D., Van Dongen, J. T., Sulpice, R., and Stitt, M. (2009). Use of reverse-phase liquid chromatography, linked to tandem mass spectrometry, to profile the Calvin cycle and other metabolic intermediates in Arabidopsis rosettes at different carbon dioxide concentrations. Plant J. 59, 824–839.
Balmer, Y., Vensel, W. H., Tanaka, C. K., Hurkman, W. J., Gelhaye, E., Rouhier, N., Jacquot, J.-P., Manieri, W., Schürmann, P., Droux, M., and Buchanan, B. B. (2004). Thioredoxin links redox to the regulation of fundamental processes of plant mitochondria. Proc. Nat. Acad. Sci. U.S.A. 101, 2642–2647.
Bauwe, H., Hagemann, M., and Fernie, A. R. (2010). Photorespiration: players, partners and origin. Trends Plant Sci. 15, 330–336.
Bauwe, H., Hagemann, M., Kern, R., and Timm, S. (2012). Photorespiration has a dual origin and manifold links to central metabolism. Curr. Opin. Plant Biol. 15, 269–275.
Baxter, C. J., Redestig, H., Schauer, N., Repsilber, D., Patil, K. R., Nielsen, J., Selbig, J., Liu, J., Fernie, A. R., and Sweetlove, L. J. (2007). The metabolic response of heterotrophic Arabidopsis cells to oxidative stress. Plant Physiol. 143, 312–325.
Beurton-Aimar, M., Beauvoit, B., Monier, A., Vallee, F., Dieuaide-Noubhani, M., and Colombie, S. (2011). Comparison between elementary flux modes analysis and 13C-metabolic fluxes measured in bacterial and plant cells. BMC Syst. Biol. 5, 95. doi: 10.1186/1752-0509-5-95
Blank, L., Kuepfer, L., and Sauer, U. (2005). Large-scale 13C-flux analysis reveals mechanistic principles of metabolic network robustness to null mutations in yeast. Genome Biol. 6, R49.
Boyle, N., and Morgan, J. (2009). Flux balance analysis of primary metabolism in Chlamydomonas reinhardtii. BMC Syst. Biol. 3, 4. doi: 10.1186/1752-0509-3-4
Brugiere, S., Kowalski, S., Ferro, M., Seigneurin-Berny, D., Miras, S., Salvi, D., Ravanel, S., D’Herin, P., Garin, J., Bourguignon, J., Joyard, J., and Rolland, N. (2004). The hydrophobic proteome of mitochondrial membranes from Arabidopsis cell suspensions. Phytochemistry 65, 1693–1707.
Cardol, P., Vanrobaeys, F., Devreese, B., Van Beeumen, J., Matagne, R. F., and Remacle, C. (2004). Higher plant-like subunit composition of mitochondrial complex I from Chlamydomonas reinhardtii: 31 conserved components among eukaryotes. Biochim. Biophys. Acta 1658, 212–224.
Carrari, F., Nunes-Nesi, A., Gibon, Y., Lytovchenko, A., Loureiro, M. E., and Fernie, A. R. (2003). Reduced expression of aconitase results in an enhanced rate of photosynthesis and marked shifts in carbon partitioning in illuminated leaves of wild species tomato. Plant Physiol. 133, 1322–1335.
Chang, R. L., Ghamsari, L., Manichaikul, A., Hom, E. F. Y., Balaji, S., Fu, W., Shen, Y., Hao, T., Palsson, B. O., Salehi-Ashtiani, K., and Papin, J. A. (2011). Metabolic network reconstruction of Chlamydomonas offers insight into light-driven algal metabolism. Mol. Syst. Biol. 7, 518.
Chen, X., Alonso, A. P., Allen, D. K., Reed, J. L., and Shachar-Hill, Y. (2011). Synergy between 13C-metabolic flux analysis and flux balance analysis for understanding metabolic adaption to anaerobiosis in E. coli. Metab. Eng. 13, 38–48.
Covert, M. W., and Palsson, B. Ø. (2002). Transcriptional regulation in constraints-based metabolic models of Escherichia coli. J. Biol. Chem. 277, 28058–28064.
Curien, G., Bastien, O., Robert-Genthon, M., Cornish-Bowden, A., Cardenas, M. L., and Dumas, R. (2009). Understanding the regulation of aspartate metabolism using a model based on measured kinetic parameters. Mol. Syst. Biol. 5, 271.
de Oliveira Dal’Molin, C. G., Quek, L.-E., Palfreyman, R., and Nielsen, L. (2011). AlgaGEM – a genome-scale metabolic reconstruction of algae based on the Chlamydomonas reinhardtii genome. BMC Genomics 12, S5. doi: 10.1186/1471-2164-12-S4-S5
de Oliveira Dal’Molin, C. G., Quek, L.-E., Palfreyman, R. W., Brumbley, S. M., and Nielsen, L. K. (2010). C4GEM, a genome-scale metabolic model to study C4 plant metabolism. Plant Physiol. 154, 1871–1885.
Duncan, O., Taylor, N. L., Carrie, C., Eubel, H., Kubiszewski-Jakubiak, S., Zhang, B., Narsai, R., Millar, A. H., and Whelan, J. (2011). Multiple lines of evidence localize signaling, morphology, and lipid biosynthesis machinery to the mitochondrial outer membrane of Arabidopsis. Plant Physiol. 157, 1093–1113.
Durot, M., Bourguignon, P.-Y., and Schachter, V. (2009). Genome-scale models of bacterial metabolism: reconstruction and applications. FEMS Microbiol. Rev. 33, 164–190.
Dutilleul, C., Lelarge, C., Prioul, J.-L., De Paepe, R., Foyer, C. H., and Noctor, G. (2005). Mitochondria-driven changes in leaf NAD status exert a crucial influence on the control of nitrate assimilation and the integration of carbon and nitrogen metabolism. Plant Physiol. 139, 64–78.
Edwards, J. S., and Palsson, B. O. (1999). Systems properties of the Haemophilus influenzae Rd metabolic genotype. J. Biol. Chem. 274, 17410–17416.
Engqvist, M., Drincovich, M. F., Flügge, U.-I., and Maurino, V. G. (2009). Two D-2-hydroxy-acid dehydrogenases in Arabidopsis thaliana with catalytic capacities to participate in the last reactions of the methylglyoxal and β-oxidation pathways. J. Biol. Chem. 284, 25026–25037.
Fahnenstich, H., Scarpeci, T. E., Valle, E. M., Fluegge, U.-I., and Maurino, V. G. (2008). Generation of hydrogen peroxide in chloroplasts of Arabidopsis overexpressing glycolate oxidase as an inducible system to study oxidative stress. Plant Physiol. 148, 719–729.
Fait, A., Fromm, H., Walter, D., Galili, G., and Fernie, A. R. (2008). Highway or byway: the metabolic role of the GABA shunt in plants. Trends Plant Sci. 13, 14–19.
Fell, D. A., Poolman, M. G., and Gevorgyan, A. (2010). Building and analysing genome-scale metabolic models. Biochem. Soc. Trans. 38, 1197–1201.
Fernie, A. (2012). Grand challenges in plant systems biology: closing the circle(s). Front. Plant Sci. 3:35. doi: 10.3389/fpls.2012.00035
Fernie, A. R., Carrari, F., and Sweetlove, L. J. (2004). Respiratory metabolism: glycolysis, the TCA cycle and mitochondrial electron transport. Curr. Opin. Plant Biol. 7, 254–261.
Fernie, A. R., and Schauer, N. (2009). Metabolomics-assisted breeding: a viable option for crop improvement? Trends Genet. 25, 39–48.
Fernie, A. R., and Stitt, M. (2012). On the discordance of metabolomics with proteomics and transcriptomics: coping with increasing complexity in logic, chemistry, and network interactions. Plant Physiol. 158, 1139–1145.
Fiehn, O., Kopka, J., Dörmann, P., Altmann, T., Trethewey, R. N., and Willmitzer, L. (2000). Metabolite profiling for plant functional genomics. Nat. Biotechnol. 18, 1157.
Fischer, E., and Sauer, U. (2005). Large-scale in vivo flux analysis shows rigidity and suboptimal performance of Bacillus subtilis metabolism. Nat. Genet. 37, 636–640.
Foyer, C. H., Noctor, G., and Hodges, M. (2011). Respiration and nitrogen assimilation: targeting mitochondria-associated metabolism as a means to enhance nitrogen use efficiency. J. Exp. Bot. 62, 1467–1482.
Fukushima, A., Kusano, M., Nakamichi, N., Kobayashi, M., Hayashi, N., Sakakibara, H., Mizuno, T., and Saito, K. (2009). Impact of clock-associated Arabidopsis pseudo-response regulators in metabolic coordination. Proc. Nat. Acad. Sci. U.S.A. 106, 7251–7256.
Garmier, M., Carroll, A. J., Delannoy, E., Vallet, C., Day, D. A., Small, I. D., and Millar, A. H. (2008). Complex I dysfunction redirects cellular and mitochondrial metabolism in Arabidopsis. Plant Physiol. 148, 1324–1341.
Giraud, E., Ho, L. H. M., Clifton, R., Carroll, A., Estavillo, G., Tan, Y.-F., Howell, K. A., Ivanova, A., Pogson, B. J., Millar, A. H., and Whelan, J. (2008). The Absence of ALTERNATIVE OXIDASE1a in Arabidopsis results in acute sensitivity to combined light and drought stress. Plant Physiol. 147, 595–610.
Gong, Q., Li, P., Ma, S., Indu Rupassara, S., and Bohnert, H. J. (2005). Salinity stress adaptation competence in the extremophile Thellungiella halophila in comparison with its relative Arabidopsis thaliana. Plant J. 44, 826–839.
Grafahrend-Belau, E., Schreiber, F., Koschützki, D., and Junker, B. H. (2009). Flux balance analysis of barley seeds: a computational approach to study systemic properties of central metabolism. Plant Physiol. 149, 585–598.
Guy, C., Kopka, J., and Moritz, T. (2008). Plant metabolomics coming of age. Physiol. Plant. 132, 113–116.
Hannah, M. A., Caldana, C., Steinhauser, D., Balbo, I., Fernie, A. R., and Willmitzer, L. (2010). Combined transcript and metabolite profiling of Arabidopsis grown under widely variant growth conditions facilitates the identification of novel metabolite-mediated regulation of gene expression. Plant Physiol. 152, 2120–2129.
Hasunuma, T., Harada, K., Miyazawa, S.-I., Kondo, A., Fukusaki, E., and Miyake, C. (2010). Metabolic turnover analysis by a combination of in vivo 13C-labelling from 13CO2 and metabolic profiling with CE-MS/MS reveals rate-limiting steps of the C3 photosynthetic pathway in Nicotiana tabacum leaves. J. Exp. Bot. 61, 1041–1051.
Hay, J., and Schwender, J. (2011). Metabolic network reconstruction and flux variability analysis of storage synthesis in developing oilseed rape (Brassica napus L.) embryos. Plant J. 67, 526–541.
Heazlewood, J. L., Howell, K. A., and Millar, A. H. (2003a). Mitochondrial complex I from Arabidopsis and rice: orthologs of mammalian and fungal components coupled with plant-specific subunits. Biochim. Biophys. Acta 1604, 159–169.
Heazlewood, J. L., Howell, K. A., Whelan, J., and Millar, A. H. (2003b). Towards an analysis of the rice mitochondrial proteome. Plant Physiol. 132, 230–242.
Heazlewood, J. L., Tonti-Filippini, J. S., Gout, A. M., Day, D. A., Whelan, J., and Millar, A. H. (2004). Experimental analysis of the Arabidopsis mitochondrial proteome highlights signaling and regulatory components, provides assessment of targeting prediction programs, and indicates plant-specific mitochondrial proteins. Plant Cell 16, 241–256.
Herz, U., Schröder, W., Liddell, A., Leaver, C. J., Brennicke, A., and Grohmann, L. (1994). Purification of the NADH:ubiquinone oxidoreductase (complex I) of the respiratory chain from the inner mitochondrial membrane of Solanum tuberosum. J. Biol. Chem. 269, 2263–2269.
Howell, K. A., Narsai, R., Carroll, A., Ivanova, A., Lohse, M., Usadel, B., Millar, A. H., and Whelan, J. (2009). Mapping metabolic and transcript temporal switches during germination in rice highlights specific transcription factors and the role of RNA instability in the germination process. Plant Physiol. 149, 961–980.
Huang, S., Taylor, N. L., Narsai, R., Eubel, H., Whelan, J., and Millar, A. H. (2009). Experimental analysis of the rice mitochondrial proteome, its biogenesis, and heterogeneity. Plant Physiol. 149, 719–734.
Ideker, T., Galitski, T., and Hood, L. (2001). A new approach to decoding life: systems biology. Ann. Rev. Genomics Hum. Genet. 2, 343–372.
Ishizaki, K., Larson, T. R., Schauer, N., Fernie, A. R., Graham, I. A., and Leaver, C. J. (2005). The critical role of Arabidopsis electron-transfer flavoprotein:ubiquinone oxidoreductase during dark-induced starvation. Plant Cell 17, 2587–2600.
Ishizaki, K., Schauer, N., Larson, T. R., Graham, I. A., Fernie, A. R., and Leaver, C. J. (2006). The mitochondrial electron transfer flavoprotein complex is essential for survival of Arabidopsis in extended darkness. Plant J. 47, 751–760.
Kaplan, F., Kopka, J., Haskell, D. W., Zhao, W., Schiller, K. C., Gatzke, N., Sung, D. Y., and Guy, C. L. (2004). Exploring the temperature-stress metabolome of Arabidopsis. Plant Physiol. 136, 4159–4168.
Kaplan, F., Kopka, J., Sung, D. Y., Zhao, W., Popp, M., Porat, R., and Guy, C. L. (2007). Transcript and metabolite profiling during cold acclimation of Arabidopsis reveals an intricate relationship of cold-regulated gene expression with modifications in metabolite content. Plant J. 50, 967–981.
Kleessen, S., Araújo, W. L., Fernie, A. R., and Nikoloski, Z. (2012). Model-based confirmation of alternative substrates of mitochondrial electron transport chain. J. Biol. Chem. 287, 11122–11131.
Klodmann, J., and Braun, H.-P. (2011). Proteomic approach to characterize mitochondrial complex I from plants. Phytochemistry 72, 1071–1080.
Klodmann, J., Senkler, M., Rode, C., and Braun, H.-P. (2011). Defining the protein complex proteome of plant mitochondria. Plant Physiol. 157, 587–598.
Klodmann, J., Sunderhaus, S., Nimtz, M., Jänsch, L., and Braun, H.-P. (2010). Internal architecture of mitochondrial complex I from Arabidopsis thaliana. Plant Cell 22, 797–810.
Kreuzwieser, J., Hauberg, J., Howell, K. A., Carroll, A., Rennenberg, H., Millar, A. H., and Whelan, J. (2009). Differential response of gray poplar leaves and roots underpins stress adaptation during hypoxia. Plant Physiol. 149, 461–473.
Krömer, S. (1995). Respiration during photosynthesis. Annu. Rev. Plant Physiol. Plant Mol. Biol. 46, 45–70.
Kruft, V., Eubel, H., Jänsch, L., Werhahn, W., and Braun, H.-P. (2001). Proteomic approach to identify novel mitochondrial proteins in Arabidopsis. Plant Physiol. 127, 1694–1710.
Leduc, M., Tikhomiroff, C., Cloutier, M., Perrier, M., and Jolicoeur, M. (2006). Development of a kinetic metabolic model: application to Catharanthus roseus hairy root. Bioprocess Biosyst. Eng. 28, 295–313.
Lee, C. P., Eubel, H., and Millar, A. H. (2010). Diurnal changes in mitochondrial function reveal daily optimization of light and dark respiratory metabolism in Arabidopsis. Mol. Cell. Proteomics 9, 2125–2139.
Lee, C. P., Eubel, H., O’Toole, N., and Millar, A. H. (2008). Heterogeneity of the mitochondrial proteome for photosynthetic and non-photosynthetic Arabidopsis metabolism. Mol. Cell. Proteomics 7, 1297–1316.
Lee, C. P., Eubel, H., O’Toole, N., and Millar, A. H. (2011). Combining proteomics of root and shoot mitochondria and transcript analysis to define constitutive and variable components in plant mitochondria. Phytochemistry 72, 1092–1108.
Lehmann, M., Laxa, M., Sweetlove, L., Fernie, A., and Obata, T. (2012). Metabolic recovery of Arabidopsis thaliana roots following cessation of oxidative stress. Metabolomics 8, 143–153.
Lehmann, M., Schwarzländer, M., Obata, T., Sirikantaramas, S., Burow, M., Olsen, C. E., Tohge, T., Fricker, M. D., Møller, B. L., Fernie, A. R., Sweetlove, L. J., and Laxa, M. (2009). The metabolic response of Arabidopsis roots to oxidative stress is distinct from that of heterotrophic cells in culture and highlights a complex relationship between the levels of transcripts, metabolites, and flux. Mol. Plant 2, 390–406.
Leterme, S., and Boutry, M. (1993). Purification and preliminary characterization of mitochondrial complex I (NADH:Ubiquinone reductase) from broad bean (Vicia faba L.). Plant Physiol. 102, 435–443.
Linka, N., and Weber, A. P. M. (2010). Intracellular metabolite transporters in plants. Mol. Plant 3, 21–53.
Lisec, J., Schauer, N., Kopka, J., Willmitzer, L., and Fernie, A. R. (2006). Gas chromatography mass spectrometry-based metabolite profiling in plants. Nat. Protoc. 1, 387–396.
Liu, S., Lin, L., Jiang, P., Wang, D., and Xing, Y. (2011). A comparison of RNA-Seq and high-density exon array for detecting differential gene expression between closely related species. Nucleic Acids Res. 39, 578–588.
Lonien, J., and Schwender, J. (2009). Analysis of metabolic flux phenotypes for two Arabidopsis mutants with severe impairment in seed storage lipid synthesis. Plant Physiol. 151, 1617–1634.
Lunn, J. E., Feil, R., Hendriks, J. H. M., Gibon, Y., Morcuende, R., Osuna, D., Scheible, W.-R., Carillo, P., Hajirezaei, M.-R., and Stitt, M. (2006). Sugar-induced increases in trehalose 6-phosphate are correlated with redox activation of ADPglucose pyrophosphorylase and higher rates of starch synthesis in Arabidopsis thaliana. Biochem. J. 397, 139–148.
Luo, B., Groenke, K., Takors, R., Wandrey, C., and Oldiges, M. (2007). Simultaneous determination of multiple intracellular metabolites in glycolysis, pentose phosphate pathway and tricarboxylic acid cycle by liquid chromatography–mass spectrometry. J. Chromatogr. A 1147, 153–164.
Martí, M. C., Olmos, E., Calvete, J. J., Díaz, I., Barranco-Medina, S., Whelan, J., Lázaro, J. J., Sevilla, F., and Jiménez, A. (2009). Mitochondrial and nuclear localization of a novel pea thioredoxin: identification of its mitochondrial target proteins. Plant Physiol. 150, 646–657.
Martin, J. A., and Wang, Z. (2011). Next-generation transcriptome as-sembly. Nat. Rev. Genet. 12, 671–682.
Meyer, E. H., Tomaz, T., Carroll, A. J., Estavillo, G., Delannoy, E., Tanz, S. K., Small, I. D., Pogson, B. J., and Millar, A. H. (2009). Remodeled respiration in ndufs4 with low phosphorylation efficiency suppresses Arabidopsis germination and growth and alters control of metabolism at night. Plant Physiol. 151, 603–619.
Millar, A. H., Heazlewood, J. L., Kristensen, B. K., Braun, H.-P., and Møller, I. M. (2005). The plant mitochondrial proteome. Trends Plant Sci. 10, 36–43.
Millar, A. H., Sweetlove, L. J., Giegé, P., and Leaver, C. J. (2001). Analysis of the Arabidopsis mitochondrial proteome. Plant Physiol. 127, 1711–1727.
Millar, A. H., Whelan, J., Soole, K. L., and Day, D. A. (2011). Organization and regulation of mitochondrial respiration in plants. Annu. Rev. Plant Biol. 62, 79–104.
Mintz-Oron, S., Meir, S., Malitsky, S., Ruppin, E., Aharoni, A., and Shlomi, T. (2012). Reconstruction of Arabidopsis metabolic network models accounting for subcellular compartmentalization and tissue-specificity. Proc. Natl. Acad. Sci. U.S.A. 109, 339–344.
Møller, I. M., and Sweetlove, L. J. (2010). ROS signalling – specificity is required. Trends Plant Sci. 15, 370–374.
Morgan, M. J., Lehmann, M., Schwarzlander, M., Baxter, C. J., Sienkiewicz-Porzucek, A., Williams, T. C. R., Schauer, N., Fernie, A. R., Fricker, M. D., Ratcliffe, R. G., Sweetlove, L. J., and Finkemeier, I. (2008). Decrease in manganese superoxide dismutase leads to reduced root growth and affects tricarboxylic acid cycle flux and mitochondrial redox homeostasis. Plant Physiol. 147, 101–114.
Moxley, J. F., Jewett, M. C., Antoniewicz, M. R., Villas-Boas, S. G., Alper, H., Wheeler, R. T., Tong, L., Hinnebusch, A. G., Ideker, T., Nielsen, J., and Stephanopoulos, G. (2009). Linking high-resolution metabolic flux phenotypes and transcriptional regulation in yeast modulated by the global regulator Gcn4p. Proc. Natl. Acad. Sci. U.S.A. 106, 6477–6482.
Narsai, R., Howell, K. A., Carroll, A., Ivanova, A., Millar, A. H., and Whelan, J. (2009). Defining core metabolic and transcriptomic responses to oxygen availability in rice embryos and young seedlings. Plant Physiol. 151, 306–322.
Noguchi, K., and Yoshida, K. (2008). Interaction between photosynthesis and respiration in illuminated leaves. Mitochondrion 8, 87–99.
Nunes-Nesi, A., Araújo, W. L., and Fernie, A. R. (2011). Targeting mitochondrial metabolism and machinery as a means to enhance photosynthesis. Plant Physiol. 155, 101–107.
Nunes-Nesi, A., Carrari, F., Gibon, Y., Sulpice, R., Lytovchenko, A., Fisahn, J., Graham, J., Ratcliffe, R. G., Sweetlove, L. J., and Fernie, A. R. (2007). Deficiency of mitochondrial fumarase activity in tomato plants impairs photosynthesis via an effect on stomatal function. Plant J. 50, 1093–1106.
Nunes-Nesi, A., Carrari, F., Lytovchenko, A., Smith, A. M. O., Loureiro, M. E., Ratcliffe, R. G., Sweetlove, L. J., and Fernie, A. R. (2005). Enhanced photosynthetic performance and growth as a consequence of decreasing mitochondrial malate dehydrogenase activity in transgenic tomato plants. Plant Physiol. 137, 611–622.
Nunes-Nesi, A., Sulpice, R., Gibon, Y., and Fernie, A. R. (2008). The enigmatic contribution of mitochondrial function in photosynthesis. J. Exp. Bot. 59, 1675–1684.
Oliver, S. G., Winson, M. K., Kell, D. B., and Baganz, F. (1998). Systematic functional analysis of the yeast genome. Trends Biotechnol. 16, 373–378.
Osorio, S., Do, P. T., and Fernie, A. R. (2012). Profiling primary metabolites of tomato fruit with gas chromatography/mass spectrometry. Methods Mol. Biol. 860, 101–109.
Palmieri, F., Pierri, C. L., De Grassi, A., Nunes-Nesi, A., and Fernie, A. R. (2011). Evolution, structure and function of mitochondrial carriers: a review with new insights. Plant J. 66, 161–181.
Perales, M., Eubel, H., Heinemeyer, J., Colaneri, A., Zabaleta, E., and Braun, H. P. (2005). Disruption of a nuclear gene encoding a mitochondrial gamma carbonic anhydrase reduces complex I and supercomplex I + III2 levels and alters mitochondrial physiology in Arabidopsis. J. Mol. Biol. 350, 263–277.
Poolman, M. G., Miguet, L., Sweetlove, L. J., and Fell, D. A. (2009). A genome-scale metabolic model of Arabidopsis and some of its properties. Plant Physiol. 151, 1570–1581.
Price, N. D., Papin, J. A., Schilling, C. H., and Palsson, B. O. (2003). Genome-scale microbial in silico models: the constraints-based approach. Trends Biotechnol. 21, 162–169.
Raghavendra, A. S., and Padmasree, K. (2003). Beneficial interactions of mitochondrial metabolism with photosynthetic carbon assimilation. Trends Plant Sci. 8, 546–553.
Ramírez-Aguilar, S. J., Keuthe, M., Rocha, M., Fedyaev, V. V., Kramp, K., Gupta, K. J., Rasmusson, A. G., Schulze, W. X., and Van Dongen, J. T. (2011). The composition of plant mitochondrial supercomplexes changes with oxygen availability. J. Biol. Chem. 286, 43045–43053.
Ratcliffe, R. G., and Shachar-Hill, Y. (2006). Measuring multiple fluxes through plant metabolic networks. Plant J. 45, 490–511.
Roessner, U., Luedemann, A., Brust, D., Fiehn, O., Linke, T., Willmitzer, L., and Fernie, A. R. (2001). Metabolic profiling allows comprehensive phenotyping of genetically or environmentally modified plant systems. Plant Cell 13, 11–29.
Rohwer, J. M., and Botha, F. C. (2001). Analysis of sucrose accumulation in the sugar cane culm on the basis of in vitro kinetic data. Biochem. J. 358, 437–445.
Saha, R., Suthers, P. F., and Maranas, C. D. (2011). Zea mays RS1563: a comprehensive genome-scale metabolic reconstruction of maize metabolism. PLoS ONE 6, e21784. doi: 10.1371/journal.pone.0021784
Sanchez, D. H., Lippold, F., Redestig, H., Hannah, M. A., Erban, A., Kramer, U., Kopka, J., and Udvardi, M. K. (2008). Integrative functional genomics of salt acclimatization in the model legume Lotus japonicus. Plant J. 53, 973–987.
Sanchez, D. H., Schwabe, F., Erban, A., Udvardi, M. K., and Kopka, J. (2012). Comparative metabolomics of drought acclimation in model and forage legumes. Plant Cell Environ. 35, 136–149.
Schallau, K., and Junker, B. H. (2010). Simulating plant metabolic pathways with enzyme-kinetic models. Plant Physiol. 152, 1763–1771.
Schauer, N., and Fernie, A. R. (2006). Plant metabolomics: towards biological function and mechanism. Trends Plant Sci. 11, 508–516.
Schauer, N., Semel, Y., Roessner, U., Gur, A., Balbo, I., Carrari, F., Pleban, T., Perez-Melis, A., Bruedigam, C., Kopka, J., Willmitzer, L., Zamir, D., and Fernie, A. R. (2006). Comprehensive metabolic profiling and phenotyping of interspecific introgression lines for tomato improvement. Nat. Biotechnol. 24, 447–454.
Schwarzländer, M., König, A.-C., Sweetlove, L. J., and Finkemeier, I. (2012). The impact of impaired mitochondrial function on retrograde signalling: a meta-analysis of transcriptomic responses. J. Exp. Bot. 63, 1735–1750.
Schwender, J., Shachar-Hill, Y., and Ohlrogge, J. B. (2006). Mitochondrial metabolism in developing embryos of Brassica napus. J. Biol. Chem. 281, 34040–34047.
Shulaev, V., Cortes, D., Miller, G., and Mittler, R. (2008). Metabolomics for plant stress response. Physiol. Plant. 132, 199–208.
Siahpoosh, M. R., Sanchez, D. H., Schlereth, A., Scofield, G. N., Furbank, R. T., Van Dongen, J. T., and Kopka, J. (2012). Modification of OsSUT1 gene expression modulates the salt response of rice Oryza sativa cv. Taipei 309. Plant Sci. 182, 101–111.
Sienkiewicz-Porzucek, A., Nunes-Nesi, A., Sulpice, R., Lisec, J., Centeno, D. C., Carillo, P., Leisse, A., Urbanczyk-Wochniak, E., and Fernie, A. R. (2008). Mild reductions in mitochondrial citrate synthase activity result in a compromised nitrate assimilation and reduced leaf pigmentation but have no effect on photosynthetic performance or growth. Plant Physiol. 147, 115–127.
Sienkiewicz-Porzucek, A., Sulpice, R., Osorio, S., Krahnert, I., Leisse, A., Urbanczyk-Wochniak, E., Hodges, M., Fernie, A. R., and Nunes-Nesi, A. (2010). Mild reductions in mitochondrial NAD-dependent isocitrate dehydrogenase activity result in altered nitrate assimilation and pigmentation but do not impact growth. Mol. Plant 3, 156–173.
Somerville, C., Bauer, S., Brininstool, G., Facette, M., Hamann, T., Milne, J., Osborne, E., Paredez, A., Persson, S., Raab, T., Vorwerk, S., and Youngs, H. (2004). Toward a systems approach to understanding plant cell walls. Science 306, 2206–2211.
Steuer, R., Nesi, A. N., Fernie, A. R., Gross, T., Blasius, B., and Selbig, J. (2007). From structure to dynamics of metabolic pathways: application to the plant mitochondrial TCA cycle. Bioinformatics 23, 1378–1385.
Studart-Guimarães, C., Fait, A., Nunes-Nesi, A., Carrari, F., Usadel, B., and Fernie, A. R. (2007). Reduced expression of succinyl-coenzyme A ligase can be compensated for by up-regulation of the γ-aminobutyrate shunt in illuminated tomato leaves. Plant Physiol. 145, 626–639.
Sulpice, R., Pyl, E.-T., Ishihara, H., Trenkamp, S., Steinfath, M., Witucka-Wall, H., Gibon, Y., Usadel, B., Poree, F., Piques, M. C., Von Korff, M., Steinhauser, M. C., Keurentjes, J. J. B., Guenther, M., Hoehne, M., Selbig, J., Fernie, A. R., Altmann, T., and Stitt, M. (2009). Starch as a major integrator in the regulation of plant growth. Proc. Natl. Acad. Sci. U.S.A. 106, 10348–10353.
Sulpice, R., Sienkiewicz-Porzucek, A., Osorio, S., Krahnert, I., Stitt, M., Fernie, A. R., and Nunes-Nesi, A. (2010a). Mild reductions in cytosolic NADP-dependent isocitrate dehydrogenase activity result in lower amino acid contents and pigmentation without impacting growth. Amino Acids 39, 1055–1066.
Sulpice, R., Trenkamp, S., Steinfath, M., Usadel, B., Gibon, Y., Witucka-Wall, H., Pyl, E.-T., Tschoep, H., Steinhauser, M. C., Guenther, M., Hoehne, M., Rohwer, J. M., Altmann, T., Fernie, A. R., and Stitt, M. (2010b). Network analysis of enzyme activities and metabolite levels and their relationship to biomass in a large panel of Arabidopsis accessions. Plant Cell 22, 2872–2893.
Sweetlove, L. J., Beard, K. F. M., Nunes-Nesi, A., Fernie, A. R., and Ratcliffe, R. G. (2010). Not just a circle: flux modes in the plant TCA cycle. Trends Plant Sci. 15, 462–470.
Sweetlove, L. J., Fait, A., Nunes-Nesi, A., Williams, T., and Fernie, A. R. (2007). The mitochondrion: an integration point of cellular metabolism and signalling. Crit. Rev. Plant Sci. 26, 17–43.
Sweetlove, L. J., and Fernie, A. R. (2005). Regulation of metabolic networks: understanding metabolic complexity in the systems biology era. New Phytol. 168, 9–24.
Sweetlove, L. J., Heazlewood, J. L., Herald, V., Holtzapffel, R., Day, D. A., Leaver, C. J., and Millar, A. H. (2002). The impact of oxidative stress on Arabidopsis mitochondria. Plant J. 32, 891–904.
Sweetlove, L. J., Lytovchenko, A., Morgan, M., Nunes-Nesi, A., Taylor, N. L., Baxter, C. J., Eickmeier, I., and Fernie, A. R. (2006). Mitochondrial uncoupling protein is required for efficient photosynthesis. Proc. Natl. Acad. Sci. U.S.A. 103, 19587–19592.
Taylor, N. L., Heazlewood, J. L., and Millar, A. H. (2011). The Arabidopsis thaliana 2-D gel mitochondrial proteome: refining the value of reference maps for assessing protein abundance, contaminants and post-translational modifications. Proteomics 11, 1720–1733.
Timm, S., Florian, A., Jahnke, K., Nunes-Nesi, A., Fernie, A. R., and Bauwe, H. (2011). The hydroxypyruvate-reducing system in Arabidopsis: multiple enzymes for the same end. Plant Physiol. 155, 694–705.
Timm, S., Nunes-Nesi, A., Pärnik, T., Morgenthal, K., Wienkoop, S., Keerberg, O., Weckwerth, W., Kleczkowski, L. A., Fernie, A. R., and Bauwe, H. (2008). A cytosolic pathway for the conversion of hydroxypyruvate to glycerate during photorespiration in Arabidopsis. Plant Cell 20, 2848–2859.
Tomaz, T., Bagard, M., Pracharoenwattana, I., Linden, P., Lee, C. P., Carroll, A. J., Stroeher, E., Smith, S. M., Gardestrom, P., and Millar, A. H. (2010). Mitochondrial malate dehydrogenase lowers leaf respiration and alters photorespiration and plant growth in Arabidopsis. Plant Physiol. 154, 1143–1157.
Tweeddale, H., Notley-Mcrobb, L., and Ferenci, T. (1998). Effect of slow growth on metabolism of Escherichia coli, as revealed by global metabolite pool (“Metabolome”) analysis. J. Bacteriol. 180, 5109–5116.
Urbanczyk-Wochniak, E., Luedemann, A., Kopka, J., Selbig, J., Roessner-Tunali, U., Willmitzer, L., and Fernie, A. R. (2003). Parallel analysis of transcript and metabolic profiles: a new approach in systems biology. EMBO Rep. 4, 989–993.
Urbanczyk-Wochniak, E., Usadel, B., Thimm, O., Nunes-Nesi, A., Carrari, F., Davy, M., Blasing, O., Kowalczyk, M., Weicht, D., Polinceusz, A., Meyer, S., Stitt, M., and Fernie, A. R. (2006). Conversion of MapMan to allow the analysis of transcript data from Solanaceous species: effects of genetic and environmental alterations in energy metabolism in the leaf. Plant Mol. Biol. 60, 773–792.
Uys, L., Botha, F. C., Hofmeyr, J. H. S., and Rohwer, J. M. (2007). Kinetic model of sucrose accumulation in maturing sugarcane culm tissue. Phytochemistry 68, 2375–2392.
van Aken, O., Zhang, B., Carrie, C., Uggalla, V., Paynter, E., Giraud, E., and Whelan, J. (2009). Defining the mitochondrial stress response in Arabidopsis thaliana. Mol. Plant 2, 1310–1324.
Wiechert, W., and Nöh, K. (2005). From stationary to instationary metabolic flux analysis. Adv. Biochem. Eng. Biotechnol. 92, 145–172.
Williams, T. C. R., Poolman, M. G., Howden, A. J. M., Schwarzlander, M., Fell, D. A., Ratcliffe, R. G., and Sweetlove, L. J. (2010). A genome-scale metabolic model accurately predicts fluxes in central carbon metabolism under stress conditions. Plant Physiol. 154, 311–323.
Xu, S., Peng, G., Wang, Y., Fang, S., and Karbowski, M. (2011). The AAA-ATPase p97 is essential for outer mitochondrial membrane protein turnover. Mol. Biol. Cell 22, 291–300.
Young, J. D., Shastri, A. A., Stephanopoulos, G., and Morgan, J. A. (2011). Mapping photoautotrophic metabolism with isotopically nonstationary 13C flux analysis. Metab. Eng. 13, 656–665.
Young, J. D., Walther, J. L., Antoniewicz, M. R., Yoo, H., and Stephanopoulos, G. (2008). An elementary metabolite unit (EMU) based method of isotopically nonstationary flux analysis. Biotechnol. Bioeng. 99, 686–699.
Zell, M. B., Fahnenstich, H., Maier, A., Saigo, M., Voznesenskaya, E. V., Edwards, G. E., Andreo, C., Schleifenbaum, F., Zell, C., Drincovich, M. F., and Maurino, V. G. (2010). Analysis of Arabidopsis with highly reduced levels of malate and fumarate sheds light on the role of these organic acids as storage carbon molecules. Plant Physiol. 152, 1251–1262.
Keywords: functional genomics, systems biology, plant respiration, omics, data integration, modeling approaches
Citation: Araújo WL, Nunes-Nesi A and Williams TCR (2012) Functional geno-mics tools applied to plant metabolism: a survey on plant respiration, its connections and the annotation of complex gene functions. Front. Plant Sci. 3:210. doi: 10.3389/fpls.2012.00210
Received: 25 June 2012; Accepted: 20 August 2012;
Published online: 06 September 2012.
Edited by:
Daniel J. Kliebenstein, University of California, Davis, USAReviewed by:
John A. Morgan, Purdue University, USARonan Sulpice, Max Planck Institute of Molecular Plant Physiology, Germany
Copyright: © 2012 Araújo, Nunes-Nesi and Williams. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in other forums, provided the original authors and source are credited and subject to any copyright notices concerning any third-party graphics etc.
*Correspondence: Wagner L. Araújo, Departamento de Biologia Vegetal, Universidade Federal de Viçosa, 36570-000 Viçosa, Minas Gerais, Brazil. e-mail: wlaraujo@ufv.br