|
|
||||||||
|
First published online September 10, 2008; 10.1104/pp.108.125880 Plant Physiology 148:1324-1341 (2008) © 2008 American Society of Plant Biologists OPEN ACCESS ARTICLE
Complex I Dysfunction Redirects Cellular and Mitochondrial Metabolism in Arabidopsis1,[W],[OA]Australian Research Council Centre of Excellence in Plant Energy Biology M316 (M.G., A.J.C., E.D., C.V., I.D.S., A.H.M.) and School of Biomedical, Biomolecular, and Chemical Sciences M310 (A.J.C.), University of Western Australia, Crawley, Western Australia 6009, Australia; and Australian Research Council Centre of Excellence in Plant Energy Biology, School of Biological Sciences, University of Sydney, New South Wales 2006, Australia (D.A.D.)
Mitochondrial complex I is a major avenue for reduced NAD oxidation linked to oxidative phosphorylation in plants. However, the plant enzyme has structural and functional features that set it apart from its counterparts in other organisms, raising questions about the physiological significance of this complex in plants. We have developed an experimental model in which rotenone, a classic complex I inhibitor, has been applied to Arabidopsis (Arabidopsis thaliana) cell suspension cultures in order to dissect early metabolic adjustments involved in cell acclimation to mitochondrial dysfunction. Rotenone induced a transitory decrease in cellular respiration (0–4 h after treatment). Cell respiration then progressively recovered and reached a steady state at 10 to 12 h after treatment. Complex I inhibition by rotenone did not induce obvious oxidative stress or cell death but affected longer term cell growth. Integrated analyses of gene expression, the mitochondrial proteome, and changes in primary metabolism indicated that rotenone treatment caused changes in mitochondrial function via alterations in specific components. A physical disengagement of glycolytic activities associated with the mitochondrial outer membrane was observed, and the tricarboxylic acid cycle was altered. Amino acid and organic acid pools were also modified by rotenone treatment, with a marked early decrease of 2-oxoglutarate, aspartate, and glutamine pools. These data demonstrate that, in Arabidopsis cells, complex I inhibition by rotenone induces significant remodeling of metabolic pathways involving the mitochondria and other compartments and point to early metabolic changes in response to mitochondrial dysfunction.
Complex I (NADH:ubiquinone oxidoreductase; EC 1.6.5.3) is a major entry point into the mitochondrial electron transport chain (ETC) of reductant generated within the mitochondrial matrix. Concomitant proton translocation by complex I is coupled to mitochondrial oxidative phosphorylation, generating ATP. Complex I is composed of over 40 different subunits in eukaryotes with a native mass of 800 to 1,000 kD (Brandt et al., 2003
As complex I represents the sole entry point of electrons from matrix NADH to the electron transport chain in mammals, mutation, damage, or inhibition of the complex has a profound effect on cellular energetics (DiMauro and Schon, 2003
Key differences exist between animal and plant complex I, as the role of the complex in plants is complicated by a series of specific peripheral enzymatic associations and the presence of complex I respiratory bypasses. For instance, the final enzyme of ascorbate synthesis in plants, converting galactolactone to ascorbate, is physically associated with complex I (Heazlewood et al., 2003
Plants containing complex I mutations exist and are viable, presumably due to the presence and activation of the type II bypass dehydrogenases that allow alternative means of oxidizing the matrix NADH pool, albeit with a lowered efficiency of coupled oxidative phosphorylation. The best studied case has been the CMSII (for cytoplasmic male sterile) mutant of tobacco (Nicotiana sylvestris), which harbors a deletion in the mitochondrial nad7 gene encoding the NAD7 subunit of complex I (Gutierres et al., 1997 However, it is extremely difficult to use the comparison of wild-type and mutant phenotypes to understand the short-term events resulting from changes in complex I activity in plants. Mutants and wild types represent distinct steady states separated by an unknown series of events and exhibit pleiotropic effects that do not necessarily portray direct roles of complex I or specific responses to complex I loss. Study of the early events associated with loss of complex I activity not only presents an opportunity to understand the transition from one metabolic state to another but also to uncover elements in the pathway of mitochondrion to nucleus signaling that must occur to initiate these events. In this study, we have used rotenone to inhibit complex I function in an Arabidopsis (Arabidopsis thaliana) cell suspension and have followed cell responses. By combining proteomic, transcriptomic, and metabolomic analyses, we provide a detailed view of how cells manage mitochondrial dysfunction. Complex I inhibition by rotenone led to the induction of alternative respiratory pathways; thus, overall cellular respiration was maintained. Transcript, protein, and metabolite analyses revealed complex metabolic adjustments, with a disengagement of mitochondria from glycolysis.
Rotenone Transiently Inhibits Cell Respiration and Affects Long-Term Cell Growth A heterotrophic cell suspension culture of Arabidopsis was grown in the dark until its midlog growth phase and treated with either 40 µM rotenone dissolved in methanol (0.25%, v/v) or methanol (to 0.25%, v/v) as a control. At specific time intervals, cell culture samples were frozen and kept for protein, RNA, and metabolite analyses, whereas fresh samples were taken for respiration, cell viability, and growth analyses. Rotenone significantly affected cell respiration over the first 4 h of treatment, decreasing respiration rates by 25% to 45% (Fig. 1A , phase I). Respiration recovered to the control value over the following 4 h, and no significant change was observed over the next 20 h (Fig. 1A, phase II). To assess the longevity of the impact of rotenone on respiration, the respiratory rate of cells treated for 48 h with rotenone was measured with or without a further rotenone addition (Supplemental Fig. S1, A and B), revealing that the pretreated cells were resistant to further rotenone treatments. To determine if the cell supernatant still contained rotenone for a substantial period, supernatant at 6, 24, and 48 h after rotenone treatment was harvested and used to inhibit fresh cells that had not been previously rotenone treated. In each case, the supernatants inhibited respiration, and this rate could not be further inhibited by further rotenone additions (Supplemental Fig. S1C). The longer term impact of rotenone treatment on cell growth and viability was also assessed over 7 d. Rotenone inhibited cell growth by nearly 10% after 16 h and by a further 20% after 48 h and up to 7 d (Fig. 1B). Measurement of cell viability by propidium iodide staining showed no significant decrease in cell viability at 7 d after rotenone treatment when compared with methanol treatment (Fig. 1C). This suggests that the cells were growing, albeit more slowly, after rotenone treatment but that they were not senescing or undergoing programmed cell death at a significantly higher rate than methanol-treated or untreated cells. These conclusions were reinforced by the analysis of the ATH1 GeneChip data from the samples, which showed significant decreases in transcripts encoding histones (MapMan BIN 28.1.3), cell cycle and cell division proteins (BINs 31.3 and 31.2), and ribosomal proteins (BIN 29.2.2), while no significant modifications (P > 0.05) in transcript levels of BINs for stress-related genes were detected (BIN 20; Supplemental Table S1). Thus, although whole cell respiration rapidly recovered from phase I inhibition, rotenone treatment affected long-term cell culture growth. Respiration as well as transcriptional and metabolic modifications induced by rotenone were then further dissected in phases I and II of the response.
Respiration Acclimation to Rotenone (Phase II) Is Associated with the Induction of Alternative Respiratory Pathways
We have shown previously that transcripts for respiratory bypasses are induced in the first hours following rotenone treatment of these cell cultures (Clifton et al., 2005
After 16 h, when the respiratory rate had recovered and reached an equilibrium (Fig. 1A, phase II), whole cell respiration was further dissected using respiratory inhibitors (Supplemental Fig. S2A). Total oxygen consumption rate, after addition of carbonyl cyanide m-chlorophenylhydrazone to remove adenylate control, was unchanged by rotenone treatment when expressed on a cell dry weight basis. Addition of 1 mM n-propyl gallate (nPG) to inhibit the AOX pathway showed that maximal respiration through the COX pathway was unchanged by rotenone treatment (Supplemental Fig. S2A). Blue-native PAGE (BN-PAGE) separation of respiratory complexes and detection of NADH dehydrogenase activity revealed that complex I was present and active in mitochondria isolated from rotenone-treated cells (Supplemental Fig. S2B), consistent with the evidence that rotenone is a reversible inhibitor and is lost during mitochondrial isolation (Singer, 1979
Addition of 1 mM KCN to whole cells to inhibit the COX pathway revealed that the maximal respiratory rate via the AOX pathway had doubled during rotenone treatment. Succinate-dependent respiratory rates of mitochondria isolated from the rotenone-treated and methanol control cells mirrored this selective induction of AOX by rotenone (Supplemental Fig. S2C). Using malate and Glu to provide substrates for complex I and internal NADH dehydrogenases through the tricarboxylic acid (TCA) cycle, we observed no change in the respiratory capacity and no change in the degree of rotenone sensitivity (Supplemental Fig. S2D), consistent with the relatively reversible nature of rotenone inhibition in vivo and its removal during mitochondrial isolation. Rates of external NADH oxidation by mitochondria isolated from rotenone-treated cells were not significantly different from those in control cell mitochondria, but again they showed that AOX activity was greatly enhanced by rotenone treatment (Fig. 2B). Western blots of isolated mitochondria using antibodies to AOX confirmed that AOX protein accumulated in mitochondria after rotenone treatment (Fig. 2C), consistent with the increased capacity of the alternative pathway. The in vivo redox state of AOX was investigated by performing rapid total membrane isolations from control and rotenone-treated cells according to Noguchi et al. (2005) Altogether, these data confirm that rotenone treatment induced an up-regulation of complex I respiratory bypasses [internal and external NAD(P)H dehydrogenases] and the AOX pathway, but without affecting the cytochrome pathway. Transcriptional changes for these components were significant at 12 h after treatment, suggesting that alternative pathway induction is involved in cell respiratory acclimation to rotenone during phase II (Fig. 1A).
Complexes I and III are primary sites of mitochondrial reactive oxygen species (mtROS) formation. In Arabidopsis cells, inhibition of complex III by antimycin A affected the redox status of mitochondria and led to a global oxidative stress (Sweetlove et al., 2002
Rotenone Induces Changes in Mitochondrial Biogenesis In order to study the effect of complex I inhibition on other mitochondrial functions, we performed a quantitative comparison of mitochondrial protein samples isolated from control and rotenone-treated cells using DIGE CyDye fluorophore technology (Fig. 4 ). Spots that increased or decreased more than 1.5-fold reproducibly over three independent experiments were selected and analyzed using liquid chromatography-tandem mass spectrometry, and the identities of the proteins are summarized in Table II . We also looked at the responses of transcripts for a broad set of 556 nucleus-encoded components of the mitochondrial proteome to 3 and 12 h of rotenone treatment (Table I; Supplemental Table S2).
Two mitochondrial proteins increased in abundance following rotenone treatment: formate dehydrogenase and arginase (Table II). Three well-known mitochondrial proteins decreased in abundance following rotenone treatment: two belong to the respiratory chain (matrix processing peptidase [MPP] β-subunit and SDH6) and one is part of the translational apparatus (elongation factor Tu). SDH6 is a plant-specific isoform of a subunit of respiratory complex II, but a catalytic role has not been identified for this protein (Millar et al., 2004 The microarray analysis showed that none of 556 known mitochondrial components changed transcript abundance significantly at 3 h (P > 0.05 after false discovery rate correction); however, by 12 h, 98 of the 556 had changed, 89 by increasing in abundance (Supplemental Table S2). Analysis of these changes in functional BINs showed a significant induction of components of the alternative pathway, protein import and fate, signaling and structure, general metabolism, heat shock proteins, mitochondrial DNA replication and transcription, and MAM33 glycoproteins after 12 h. Transcripts of genes encoding proteins involved in translation and C1 metabolism, on the other hand, were slightly down at 3 h (Table I).
Seven other protein spots that decreased in abundance by 1.5- to 5.5-fold in isolated mitochondria following rotenone treatment (Fig. 4; Table II) were identified as major cytosolic enzymes of primary carbon metabolism. Notably among these were isoforms of aldolase, glyceraldehyde-3-P dehydrogenase (GAPDH), enolase, and triosephosphate isomerase, all of which have been shown to selectively bind to plant mitochondria as a functional glycolytic unit on the outer membrane (Giegé et al., 2003
Three other cytosolic carbon metabolism enzymes were also decreased in mitochondrial extracts: NAD-malate dehydrogenase, NADP-malic enzyme (ME), and NADP-isocitrate dehydrogenase (ICDH; Table II); all three of these are directly linked to the import and export of carbon skeletons in mitochondria. Cytosolic malate dehydrogenase and NAD-ME are involved in the connections between malate, OAA, and pyruvate provision to mitochondria and thus can be considered as extensions of the glycolytic pathway, while cytosolic ICDH is involved in converting isocitrate exported from mitochondria to 2-oxoglutarate as a carbon skeleton for nitrogen assimilation. By analogy with the claims made for the functional association of classical glycolytic enzymes with mitochondria, these proteins may also be associated in a functional manner, connecting mitochondrial metabolism with the broader cellular metabolism. The apparent dissolution of these linkages may indicate a degree of disengagement of mitochondria from cytosolic carbon metabolism.
Quantitative profiles of the abundance of the major metabolites from rotenone-treated cells over the first 24 h of treatment were generated using gas chromatography-mass spectrometry (GC-MS) analysis of derivatized compounds from methanol-soluble cell extracts. A series of 48 significant changes were recorded in response to the treatment in a time-dependent fashion. These included 18 metabolites that significantly increased and 30 metabolites that significantly decreased in abundance during at least one point of the time series. In Figure 6
, a heat map of the abundance of major compounds unambiguously identified by comparison with library standards and classified into amino acids, organic acids, sugars, and other compounds is presented. Within each grouping, the metabolites are ordered according to their degree of change at the 16-h time point, where most of the measurements reported in Figures 2 to 5
The first significant changes recorded were the approximate halving of 2-oxoglutarate concentration and the approximate 2.5-fold increases in Trp and Lys within 1 h of rotenone treatment. Interestingly, the decrease in 2-oxoglutarate was sustained for at least 16 h, while Trp and Lys returned to control levels within 3 and 6 h, respectively. These changes were rapidly followed by decreases in Gln, Glu, Asp, and Asn, all significantly reduced within 3 to 6 h. These rapid changes are consistent with a slowing in TCA cycle function. At 12 to 16 h, there was evidence of decreases in a variety of TCA cycle intermediates from the decarboxylating portion of the cycle, namely citrate, aconitate, isocitrate, and of course 2-oxoglutarate, while there were trends to increases in abundance for fumarate and malate in the other side of the cycle. The activity of TCA cycle enzymes, on the other hand, showed either no decrease or a slight increase (Table III ), suggesting that the slowing of the TCA cycle could be caused by an increase in the redox poise of the matrix NADH pool upon inhibition of complex I by rotenone. Pyruvate, 2-oxoglutarate, isocitrate, and malate dehydrogenases are all sensitive to the redox poise of the NADH pool and dramatically decrease their activity as the NADH to NAD ratio increases (for review, see Noctor et al., 2007
Accumulation of carbon intermediates at the end of the glycolytic pathway led to the accumulation of Ala (first recorded at 6 h) and the pyruvate-derived branched chain amino acids Leu and Val (first recorded at 12 h). By 6 h, there were already net increases in the pyruvate-derived fermentation product lactate. This is consistent with transcriptional increases in the components of the glycolytic pathway (BINs 4.9 and 5) and branched chain amino acid degradation pathways (BINs 13.2.4.1 and 13.2.4.4; Supplemental Table S1) and increased whole cell glycolytic enzymatic activities at 16 h (Fig. 5). Phosphoenolpyruvate-derived amino acids and compounds in the shikimate pathway decreased from 3 h (Fig. 6) with the exception of Trp, which first increased at 1 h before dropping over the time course. Despite the decrease in Asp and Asn, increases were noted in the homoserine branch of the Asp family of amino acids, notably in Ile, Met, and 1-aminocyclopropane carboxylic acid, while aminoadipic acid in Lys catabolism decreased by half within 16 h. Losses in both the Asp and Glu pools might be expected to alter the biosynthetic functions of the cell, and this can be seen in the decrease in allantoin from the purine synthesis pathway. Decreases in Glu were consistent with decreases in Orn, and the concomitant loss of Asp appears to have lowered the urea cycle and decreased urea abundance by 24 h. As Glu is also the precursor for Pro and thus 4-Hyp, the decrease of the latter compound was also likely to have stemmed from the lowering of TCA cycle activity and 2-oxoglutarate availability (Fig. 7). Late decreases (12–24 h) were also recorded in hexose phosphate pools that either fed, or were fed by, changes in several Glc-derived sugars. However, hexoses have largely reverted to control levels by 24 h.
The aim of this work was to identify the early responses of Arabidopsis cells involved in acclimation to complex I dysfunction. The long-term consequences of complex I loss have been studied in several mutant plants, but it is difficult to unravel the initial events involved in the modification of primary metabolism in these plants.
We showed that rotenone treatment induced a two-phase response of Arabidopsis cells: phase I, from 0 to 4 h after treatment, characterized by a strong inhibition of cellular respiration by rotenone; and phase II, from 4 to 32 h, characterized by a progressive recovery of cell respiration to initial rates (Fig. 1A). Phase II was associated with transcriptional and posttranscriptional induction of alternative respiratory pathways, complex I bypasses [NAD(P)H dehydrogenases], and AOX, suggesting that these processes were part of cell acclimation to complex I inhibition. Up-regulation of alternative pathways is a classic adaptive response of plants lacking functional complex I, such as the tobacco CMSII mutant or the Arabidopsis otp43 mutant (Sabar et al., 2000
The use of chemical inhibitors raises the possibility of pleiotrophic effects, as both mitochondria and chloroplasts can contain rotenone targets. We used a concentration of 40 µM because previous experiments have optimized this concentration for cell culture growth and survival (Lister et al., 2004 From our data, changes in enzyme and protein were correlated with changes in metabolite and transcriptional profiles over the first 24 h of rotenone treatment. These changes included alterations of specific pathways of mitochondrial electron transport, the source of reductant in the matrix and cytosol, and the disengagement of mitochondria from glycolysis and related carbon metabolism. Although most of the changes in transcript, protein, and metabolite levels were observed during phase II "acclimation" (6–24 h), some changes occurred in phase I (1–3 h; Fig. 6), preceding measurable changes in nuclear gene transcription.
Loss of complex I appeared to lead to a degree of disconnection between glycolysis, the TCA cycle, and the electron transport chain. We saw this through a series of cellular changes that reflect this disconnection and that could also be interpreted as attempts to reconcile this problem through the enhancement of alternative metabolic pathways. Even with the availability of time series data, it has been relatively difficult to neatly place these events in series, but it is clear that progressive changes occur in only a few hours that lead to a new steady state appearing within 10 to 12 h. The early decreases in 2-oxoglutarate and Gln/Glu are consistent with slowing of the TCA cycle due to a loss of a major entry point for NADH, complex I. This appears to set in motion a series of changes in metabolism that alter the carbon sources used to drive respiration and also increase capacity to bypass complex I as an entry point, in a complex and interrelated manner.
The loss of glycolytic enzymes from mitochondria (Table II; Fig. 5) suggests a degree of physical disengagement of mitochondria from glycolysis. The observed accumulation of glycolytic end products and associated amino acid pools, together with decreases in early TCA cycle-derived products, suggests that this physical disengagement leads to a degree of functional disengagement. This is consistent with the recent evidence provided by Graham et al. (2007) The substrates for these compensatory pathways may, in fact, be provided through the disconnection of glycolysis from mitochondria. As pyruvate cannot enter mitochondria or is not rapidly used in mitochondria, it can be used for the synthesis of Leu and Val, which are key substrates of the branched chain degradation pathway. There is clear transcriptional evidence for the induction of branched chain catabolism components located both within the mitochondria (Supplemental Table S2, BIN 12) and more generally in the cell (Supplemental Table S1, BIN 13) and clear metabolomic data for the increased availability of these amino acids (Figs. 6 and 7). Increased glycolytic flux disconnected from the mitochondria also leads to fermentation and lactate formation, shown by both the transcript and metabolite analysis (Fig. 6; Supplemental Table S1). This raises the potential for enhanced cytosolic NADH levels that could act as the substrate for the external rotenone-insensitive NADH dehydrogenases that are induced by rotenone treatment.
Formate oxidation generates NADH in the matrix, and formate could be provided as a consequence of glycolytic interruption, either from a fermentation aldehyde product or via the predicted pyruvate-lyase reaction. Alternatively, it could be generated via the Met salvage pathway/Yang cycle, as Met and 1-aminocyclopropane carboxylic acid are two of the more rapidly increased metabolites and Met was the only metabolite to significantly increase and then significantly decrease during the 24-h time course (Fig. 6). Interestingly, formate dehydrogenase is also elevated after disruption of a plant-specific subunit of complex I, which leads to a substantial loss of assembly of the respiratory complex (Perales et al., 2005
Details of the signal transduction pathways that alter nuclear gene expression upon mitochondrial dysfunction remain elusive in plants. Mitochondrial dysfunction leads to changes in a large number of metabolites, and a single one or a combination of several might be involved in the signaling process.
Redox-based signaling is often raised as a key component in mitochondria-nucleus communication (Noctor et al., 2007
Interestingly, many of the components transcriptionally induced by rotenone belong to the machinery involved in mitochondrial transcription, translation, ETC organization/assembly, and protein fate (Table I; Supplemental Table S2), showing that cells respond by increasing their capacity for mitochondrial biogenesis and ETC protection. Rotenone also induced HSP60s, HSP70s, and HSP90-related proteins, which act as chaperones and ensure the correct function of proteins by preventing the aggregation of denatured proteins, by refolding of stress-denatured proteins, and by assisting in the rapid assembly of the oligomeric protein structures. HSP60 and HSP70 have also been reported to be involved in Arabidopsis cell tolerance to heat stress and to respiration deficiency (Kuzmin et al., 2004
Many stress-responsive genes are repressed during normal plant development, and stress-induced expression could be due to the unbinding of a repressor rather than the binding of an activator molecule. Given this, the initial signal might also be loss of a normally present component rather than either the extraordinary accumulation of a primary metabolite or the generation of a specific "stress signaling" component. In this context, 2-oxoglutarate, Asp, and Gln pools rapidly responded to rotenone inhibition, decreasing significantly by 1 h (Fig. 6) and preceding transcriptional responses of genes for nucleus-encoded mitochondrial components (Fig. 2A; Table I). Consistent and prolonged elevation of other metabolic components (Fig. 6) did not occur significantly until some hours later, and in many instances this could have been caused by the loss of these primary metabolites. In yeast, evidence has accumulated that 2-oxoglutarate and Glu are probable signals in retrograde regulation of mitochondrial function. Glu is a potent repressor of the retrograde-dependent expression of genes that leads to 2-oxoglutarate formation (notably peroxisomal and mitochondrial citrate synthase and mitochondrial aconitase and ICDH), so loss of Glu leads to enhanced expression of these enzymes in a pathway allowing truncated TCA cycle function during respiratory deficiency (Liu and Butow, 1999
Complex I dysfunction has been reported to have various effects on plant growth and development. Rotenone treatment did not alter cell viability but reduced growth (Fig. 1). A similar growth reduction was seen in suspension cell cultures generated from Arabidopsis plants with a knockout of a
As observed after rotenone treatment, AOX and/or alternative dehydrogenase pathways are increased in leaves of different complex I mutant plants, at the level of gene expression, protein accumulation, or respiration rates (Sabar et al., 2000
The metabolite profile of CMSII leaves is enriched in amino acids with low carbon to nitrogen ratios and depleted in both starch and 2-oxoglutarate (Dutilleul et al., 2005
Using rotenone inhibition in Arabidopsis cell cultures, we have shown that plant cells respond to the loss of complex I function by readjusting their electron transport properties and general cellular carbon and nitrogen metabolism to minimize oxidative stress and allow cell survival. This provides insights into the molecular and metabolic phenotypes of complex I mutants and also highlights very early changes in metabolite concentrations that may function in retrograde signaling to the nucleus upon changes in mitochondrial function.
Arabidopsis Cell Culture, Plant Growth, and Treatments
A heterotrophic cell suspension culture of Arabidopsis (Arabidopsis thaliana ecotype Landsberg erecta) was maintained on Murashige and Skoog basal salt medium as described by Sweetlove et al. (2002)
Total RNA from treated and control samples was extracted using the Qiagen Plant RNeasy kit, and genomic DNA was removed using DNA-free DNase (Ambion). Complete removal of both mitochondrial and nuclear DNA was checked by PCR on diluted RNA prior to RT with random primers (SuperScript III; Invitrogen). Quantitative PCR was performed with the LightCycler 480 real-time PCR system (Roche) using the LightCycler 480 SYBR Green 1 Master Mix and the primer sets detailed in Supplemental Table S3, originally developed by Clifton et al. (2005)
Mitochondria were isolated and purified by differential centrifugations followed by two Percoll step gradients according to the method of Millar et al. (2001)
Oxygen consumption was measured using a Clark-type oxygen electrode (Hansatech Instrument) in 1 mL of reaction medium. Respiration on whole Arabidopsis cell culture was performed at 22°C by suspending 300 µL of cell culture in 700 µL of fresh cell medium. Total respiration was measured after the addition of the uncoupler carbonyl cyanide m-chlorophenylhydrazone (4 µM). Subsequent addition of 1 mM nPG to inhibit the AOX pathway or 1 mM KCN to inhibit the COX pathway allowed measurements of maximal activity of the COX and AOX pathways, respectively. Respiration rates on isolated mitochondria were determined at 25°C on 200 µg of mitochondrial proteins suspended in the reaction medium (0.3 M Suc, 5 mM KH2PO4, 10 mM TES, 10 mM NaCl, 2 mM MgSO4, and 0.1% bovine serum albumin [w/v], pH 7.2). Different substrates, cofactors, and inhibitors were successively added to the reaction medium to modulate oxygen consumption by mitochondria. Total respiration was measured in the presence of succinate (10 mM) and ATP (0.5 mM). Maximal capacity of the AOX pathway was measured after successive additions of the complex III inhibitor myxothiazol (2.5 µM), activators of AOX, pyruvate (10 mM), and the reductant DTT (5 mM), and finally the AOX inhibitor nPG (0.5 mM). Respiration rate through the COX pathway was measured in the presence of nPG (0.5 mM) and the complex IV inhibitor KCN (1 mM). Total respiration was also measured in the presence of Glu (10 mM) and malate (10 mM), which provide substrates for complex I and internal NADH dehydrogenases through the TCA cycle. In this case, the pH of the reaction medium was adjusted to 6.8. To Glu and malate were added coenzyme A (12 µM), thiamine pyrophosphate (0.2 mM), NAD+ (2 mM), and ADP (0.3 mM) in order to initiate electron transport. Rotenone (5 mM), nPG (0.5 mM), and KCN (1 mM) were then added. The capacity of mitochondria to use external NADH was measured by additions of NADH (1 mM), CaCl2 (0.1 mM), and ADP (0.3 mM) and then myxothiazol (2.5 µM), DTT (5 mM)/pyruvate (10 mM), and nPG (0.5 mM). Mitochondrial integrity was measured by following oxygen uptake after the addition of ascorbate (10 mM), cytochrome c (25 µM), and Triton X-100 (0.05%, w/v). In all experiments, it was consistently better than 90%.
In order to study AOX content and redox state in treated cells, one-dimensional SDS-PAGE and western blot analysis were performed on mitochondrial proteins or on membrane proteins quickly extracted from whole cells. For mitochondrial proteins, 50 µg of proteins was solubilized in sample buffer (2% SDS, 125 mM Tris-HCl, 10% glycerol, 10% β-mercaptoethanol, and 0.002% bromphenol blue, pH 6.8) and heated at 95°C for 3 min before loading onto a polyacrylamide gel. For isolated cell membranes, 0.5 g of cells was ground at 4°C in 1 mL of an extraction buffer (0.45 M mannitol, 50 mM TES, 0.5% bovine serum albumin, 0.5% polyvinylpyrrolidone-40, 2 mM EGTA, 20 mM ascorbate, and a proteinase inhibitor tablet [Roche], pH 8.0) as described by Noguchi et al. (2005)
BN-PAGE was performed according to the method described by Jansch et al. (1996)
Complex I (NADH dehydrogenase) activity was measured after protein separation by BN-PAGE according to Zerbetto et al. (1997)
The two-dimensional fluorescence difference gel electrophoresis (2D-DIGE) technology was used to allow comparison on one gel of mitochondrial protein contents from two different samples: mitochondrial proteins from control cells (treated with methanol) and mitochondrial proteins from rotenone-treated cells. Proteins of each sample were first individually labeled with CyDye fluorophores according to the manufacturer's instructions (GE Healthcare). Proteins were precipitated by incubating purified mitochondria overnight at –20°C in acetone. After centrifugation at 20,000g, 4°C for 20 min, pellets were resuspended in 10 µL of DIGE lysis buffer (8 M urea, 4% CHAPS, and 40 mM Tris, pH 8.5) and centrifuged again at 12,000g, 4°C for 10 min in order to remove unsolubilized material. Fifty micrograms of proteins of each sample (proteins from methanol- and rotenone-treated cells) were stained with a different CyDye (Cy-2 or Cy-5) by addition of 400 pmol of dye (freshly diluted in dimethylformamide) and incubation on ice in the dark for 30 min. The reaction was stopped by incubating samples for 10 min on ice in the dark with 10 nmol of Lys. An equal volume (12 µL) of DIGE lysis buffer with 22 mM DTT was added to each sample. Each of the labeled protein samples was mixed, and rehydration buffer (8 M urea, 4% CHAPS, 0.5% 3–10 nonlinear immobilized pH gradient buffer, 18 mM DTT, and 0.001% [w/v] bromphenol blue) was added to give a final volume of 450 µL. The mix was loaded onto a strip of 24 cm in length with an immobilized nonlinear pH gradient of 3 to 10 (Immobiline DryStrip; GE Healthcare). Rehydration of the strips and the first IEF dimension electrophoresis were performed on an IPGphor Unit (GE Healthcare) using the following settings: 12 h at 30 V (rehydration step), 1 h at 500 V, 1-h gradient from 500 to 1,000 V, 1-h gradient from 1,000 to 3,000 V, 2-h gradient from 3,000 to 8,000 V, and 5 h at 8,000 V. After IEF, strips were incubated for 15 min in an equilibration buffer (6 M urea, 2% SDS, 26% glycerol, 65 mM DTT, 0.001% [w/v] bromphenol blue, and 50 mM Tris-HCl, pH 8.8) and then for 15 min in an equilibration buffer containing iodoacetamide (6 M urea, 2% SDS, 26% glycerol, 135 mM DTT, 0.001% [w/v] bromphenol blue, and 50 mM Tris-HCl, pH 8.8). The equilibrated strips were then loaded on top of a 12% acrylamide gel. Following separation, gels were scanned using the Typhoon Trio Variable Mode Imager at a resolution of 100 (pixel size) with the photomultiplier tube set to 500 V. Proteins were visualized using the Typhoon Scanner Control software (version 5.0) and processed (quantification) using DeCyder 2-D Differential Analysis software version 6.5 (GE Healthcare). In order to get statistical significance from these experiments, three sets of proteins (control and "rotenone" proteins) from three independent experiments were labeled and submitted to electrophoresis. Standard gels were also used, and precipitation, IEF, and SDS-PAGE were run in parallel with labeled samples. These standard gels were loaded with a mix of 150 µg of each protein sample (control and rotenone). After electrophoresis, proteins were visualized by colloidal Coomassie Brilliant Blue G250 staining (Sweetlove et al., 2002
Gel spots to be analyzed were cut from colloidal Coomassie Brilliant Blue-stained gels. Proteins were digested with trypsin (trypsin sequencing grade; Roche Diagnostic) according to Sweetlove et al. (2001)
Activities were measured on purified mitochondria according to methods optimized for plant mitochondrial activities: aconitase (EC 4.2.1.3), NAD-ME (EC 1.1.1.39), pyruvate dehydrogenase complex (PDC), and 2-oxoglutarate dehydrogenase complex (OGDC) as described by Millar et al. (1999)
Metabolites were extracted from frozen Arabidopsis cells according to a method adapted from Roessner-Tunali et al. (2003)
Derivatized metabolite samples were analyzed on an Agilent GC/MSD system composed of an Agilent GC 6890N gas chromatograph fitted with a 7683B Automatic Liquid Sampler and 5975B Inert MSD quadrupole MS detector (Agilent Technologies). The gas chromatograph was fitted with a 0.25-mm i.d., 0.25-µm film thickness, 30-m Varian FactorFour VF-5ms capillary column with 10 m integrated guard column (Varian). GC-MS run conditions were essentially as described for GC-quadrupole-MS metabolite profiling on the Golm Metabolome Database Web site (http://csbdb.mpimp-golm.mpg.de/csbdb/gmd/analytic/gmd_meth.html; Kopka et al., 2005
After collection, cells were snap frozen in liquid nitrogen and stored at –80°C. Frozen cells were ground to a fine powder with liquid nitrogen and homogenized in 1 M perchloric acid (with a cell biomass:perchloric acid ratio of 1:5 [w/v]). After centrifugation for 15 min at 15,000g and 4°C, a volume of 0.1 mL of 0.12 M NaH2PO4 (pH 5.6) was added to 0.5 mL of supernatant and the pH was adjusted to a value between 5 and 6 with 2.5 M K2CO3. Insoluble KClO4 was removed from samples by centrifugation at 15,000g for 5 min at 4°C. Ascorbate and glutathione were determined in the same supernatant by spectrophotometry as described by Dutilleul et al. (2003b)
ANOVA was performed with StatBox 6.6 software (Grimmersoft). Data are expressed as means ± SE or SD as specified.
The following materials are available in the online version of this article.
From the Australian Research Council Centre of Excellence in Plant Energy Biology, we thank Holger Eubel, Alison M. Winger, and Etienne H. Meyer for assistance with BN-PAGE and 2D-DIGE technologies and Joshua L. Heazlewood for advice and assistance with peptide mass spectrometry. Received July 5, 2008; accepted September 5, 2008; published September 10, 2008.
1 This work was supported by the Australian Research Council (grant no. CE0561495 to A.H.M., I.D.S., and D.A.D.) through the Centre of Excellence program. M.G. was the recipient of an Australian Research Council Linkage International fellowship (no. LX0560236), A.J.C. was the recipient of a Grains Research and Development Corporation Postgraduate Award, I.D.S. was supported as a Western Australia Premier's Fellow, and A.H.M. is an Australian Research Council Australian Professorial Fellow (award no. DP0771156).
2 Present address: Institut de Biotechnologie des Plantes, Université Paris-Sud 11, CNRS, UMR 8618, Bâtiment 630, 91405 Orsay cedex, France. The author responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (www.plantphysiol.org) is: A. Harvey Millar (hmillar{at}cyllene.uwa.edu.au).
[W] The online version of this article contains Web-only data.
[OA] Open Access articles can be viewed online without a subscription. www.plantphysiol.org/cgi/doi/10.1104/pp.108.125880 * Corresponding author; e-mail hmillar{at}cyllene.uwa.edu.au.
Ayala A, Venero JL, Cano J, Machado A (2007) Mitochondrial toxins and neurodegenerative diseases. Front Biosci 12: 986–1007[CrossRef][Web of Science][Medline] Bardey V, Vallet C, Robas N, Charpentier B, Thouvenot B, Mougin A, Hajnsdorf E, Regnier P, Springer M, Branlant C (2005) Characterization of the molecular mechanisms involved in the differential production of erythrose-4-phosphate dehydrogenase, 3-phosphoglycerate kinase and class II fructose-1,6-bisphosphate aldolase in Escherichia coli. Mol Microbiol 57: 1265–1287[CrossRef][Web of Science][Medline] Betarbet R, Sherer TB, MacKenzie G, Garcia-Osuna M, Panov AV, Greenamyre JT (2000) Chronic systemic pesticide exposure reproduces features of Parkinson's disease. Nat Neurosci 3: 1301–1306[CrossRef][Web of Science][Medline] Brandt U, Kerscher S, Dröse S, Zwicker K, Zickermann V (2003) Proton pumping by NADH:ubiquinone oxidoreductase: a redox driven conformational change mechanism? FEBS Lett 54: 9–17 Clifton R, Lister R, Parker KL, Sappl PG, Elhafez D, Millar AH, Day DA, Whelan J (2005) Stress-induced co-expression of alternative respiratory chain components in Arabidopsis thaliana. Plant Mol Biol 58: 193–212[CrossRef][Web of Science][Medline] Corneille S, Cournac L, Guedeney G, Havaux M, Peltier G (1998) Reduction of the plastoquinone pool by exogenous NADH and NADPH in higher plant chloroplasts: characterization of a NAD(P)H-plastoquinone oxidoreductase activity. Biochim Biophys Acta 1363: 59–69[Medline] DiMauro S, Schon EA (2003) Mitochondrial respiratory-chain diseases. N Engl J Med 348: 2656–2668 Djajanegara I, Finnegan PM, Mathieu C, McCabe T, Whelan J, Day DA (2002) Regulation of alternative oxidase gene expression in soybean. Plant Mol Biol 50: 735–742[CrossRef][Web of Science][Medline] Dudkina NV, Eubel H, Keegstra W, Boekema EJ, Braun HP (2005) Structure of a mitochondrial supercomplex formed by respiratory-chain complexes I and III. Proc Natl Acad Sci USA 102: 3225–3229 Dutilleul C, Driscoll S, Cornic G, De Paepe R, Foyer CH, Noctor G (2003a) Functional mitochondrial complex I is required by tobacco leaves for optimal photosynthetic performance in photorespiratory conditions and during transients. Plant Physiol 131: 264–275 Dutilleul C, Garmier M, Noctor G, Mathieu C, Chetrit P, Foyer CH, De Paepe R (2003b) Leaf mitochondria modulate whole cell redox homeostasis, set antioxidant capacity, and determine stress resistance through altered signaling and diurnal regulation. Plant Cell 15: 1212–1226 Dutilleul C, Lelarge C, Prioul JL, De Paepe R, Foyer CH, 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 Elthon TE, Nickels RL, McIntosh L (1989) Monoclonal antibodies to the alternative oxidase of higher plant mitochondria. Plant Physiol 89: 1311–1317 Eubel H, Heinemeyer J, Braun HP (2004) Identification and characterization of respirasomes in potato mitochondria. Plant Physiol 134: 1450–1459 Falcon de Longevialle AF, Meyer EH, Andrés C, Taylor NL, Lurin C, Millar AH, Small ID (2007) The pentatricopeptide repeat gene OTP43 is required for trans-splicing of the mitochondrial nad1 intron 1 in Arabidopsis thaliana. Plant Cell 19: 3256–3265 Giegé P, Heazlewood JL, Roessner-Tunali U, Millar AH, Fernie AR, Leaver CJ, Sweetlove LJ (2003) Enzymes of glycolysis are functionally associated with the mitochondrion in Arabidopsis cells. Plant Cell 15: 2140–2151 Graham JW, Williams TC, Morgan M, Fernie AR, Ratcliffe RG, Sweetlove LJ (2007) Glycolytic enzymes associate dynamically with mitochondria in response to respiratory demand and support substrate channeling. Plant Cell 19: 3723–3728 Gutierres S, Sabar M, Lelandais C, Chetrit P, Diolez P, Degand H, Boutry M, Vedel F, de Kouchkovsky Y, De Paepe R (1997) Lack of mitochondrial and nuclear-encoded subunits of complex I and alteration of the respiratory chain in Nicotiana sylvestris mitochondrial deletion mutants. Proc Natl Acad Sci USA 94: 3436–3441 Heazlewood JL, Howell KA, Millar AH (2003) Mitochondrial complex I from Arabidopsis and rice: orthologs of mammalian and fungal components coupled with plant-specific subunits. Biochim Biophys Acta 1604: 159–169[Medline] Howell KA, Cheng K, Murcha MW, Jenkin LE, Millar AH, Whelan J (2007) Oxygen initiation of respiration and mitochondrial biogenesis in rice. J Biol Chem 282: 15619–15631 Igambardiev AU, Zhou G, Malmberg G, Gardestrom P (1997) Respiration of barley protoplasts before and after illumination. Physiol Plant 99: 15–22[CrossRef] Ikezawa N, Ifuku K, Endo T, Sato F (2002) Inhibition of photosystem II of spinach by the respiration inhibitors piericidin A and thenoyltrifluoroacetone. Biosci Biotechnol Biochem 66: 1925–1929[CrossRef][Medline] Ishizaki K, Schauer N, Larson TR, Graham IA, Fernie AR, Leaver CJ (2006) The mitochondrial electron transfer flavoprotein complex is essential for survival of Arabidopsis in extended darkness. Plant J 47: 751–760[CrossRef][Web of Science][Medline] Jansch L, Kruft V, Schmitz UK, Braun HP (1996) New insights into the composition, molecular mass and stoichiometry of the protein complexes of plant mitochondria. Plant J 9: 357–368[CrossRef][Web of Science][Medline] Karpova OV, Kuzmin EV, Elthon TE, Newton KJ (2002) Differential expression of alternative oxidase genes in maize mitochondrial mutants. Plant Cell 14: 3271–3284 Kopka J, Schauer N, Krueger S, Birkemeyer C, Usadel B, Bergmüller E, Dörmann P, Weckwerth W, Gibon Y, Stitt M, et al (2005) GMD{at}CSB.DB: the Golm Metabolome Database. Bioinformatics 21: 1635–1638 Kussmaul L, Hirst J (2006) The mechanism of superoxide production by NADH:ubiquinone oxidoreductase (complex I) from bovine heart mitochondria. Proc Natl Acad Sci USA 103: 7607–7612 Kuzmin EV, Karpova OV, Elthon TE, Newton KJ (2004) Mitochondrial respiratory deficiencies signal up-regulation of genes for heat shock proteins. J Biol Chem 279: 20672–20677 Lister R, Chew O, Lee MN, Heazlewood JL, Clifton R, Parker KL, Millar AH, Whelan J (2004) A transcriptomic and proteomic characterization of the Arabidopsis mitochondrial protein import apparatus and its response to mitochondrial dysfunction. Plant Physiol 134: 777–789 Liu Z, Butow RA (1999) A transcriptional switch in the expression of yeast tricarboxylic acid cycle genes in response to a reduction or loss of respiratory function. Mol Cell Biol 19: 6720–6728 Liu Z, Butow RA (2006) Mitochondrial retrograde signaling. Annu Rev Genet 40: 159–185[CrossRef][Web of Science][Medline] Maxwell DP, Wang Y, McIntosh L (1999) The alternative oxidase lowers mitochondrial reactive oxygen production in plant cells. Proc Natl Acad Sci USA 96: 8271–8276 Melo AMP, Roberts TH, Møller AM (1996) Evidence for the presence of two rotenone-insensitive NAD(P)H dehydrogenases on the inner surface of the inner membrane of potato tuber mitochondria. Biochim Biophys Acta 1276: 133–139 Meyer EH, Heazlewood JL, Millar AH (2007) Mitochondrial acyl carrier proteins in Arabidopsis thaliana are predominantly soluble matrix proteins and none can be confirmed as subunits of respiratory complex I. Plant Mol Biol 64: 319–327[CrossRef][Web of Science][Medline] Millar AH, Eubel H, Jänsch L, Kruft V, Heazlewood JL, Braun HP (2004) Mitochondrial cytochrome c oxidase and succinate dehydrogenase complexes contain plant specific subunits. Plant Mol Biol 56: 77–90[CrossRef][Web of Science][Medline] Millar AH, Hill SA, Leaver CJ (1999) Plant mitochondrial 2-oxoglutarate dehydrogenase complex: purification and characterization in potato. Biochem J 343: 327–334[CrossRef][Web of Science][Medline] Millar AH, Mittova V, Kiddle G, Heazlewood JL, Bartoli CG, Theodoulou FL, Foyer CH (2003) Control of ascorbate synthesis by respiration and its implications for stress responses. Plant Physiol 133: 443–447 Millar AH, Sweetlove LJ, Giegé P, Leaver CJ (2001) Analysis of the Arabidopsis mitochondrial proteome. Plant Physiol 127: 1711–1727 Nettleton D (2006) A discussion of statistical methods for design and analysis of microarray experiments for plant scientists. Plant Cell 18: 2112–2121 Noctor G, De Paepe R, Foyer CH (2007) Mitochondrial redox biology and homeostasis in plants. Trends Plant Sci 12: 125–134[CrossRef][Web of Science][Medline] Noguchi K, Taylor NL, Millar AH, Lambers H, Day DA (2005) Response of mitochondria to light intensity in the leaves of sun and shade species. Plant Cell Environ 28: 760–771[CrossRef] Perales M, Eubel H, Heinemeyer J, Colaneri A, Zabaleta E, Braun HP (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[CrossRef][Web of Science][Medline] Pineau B, Mathieu C, Gerard-Hirne C, De Paepe R, Chetrit P (2005) Targeting the NAD7 subunit to mitochondria restores a functional complex I and a wild type phenotype in the Nicotiana sylvestris CMS II mutant lacking nad7. J Biol Chem 15: 25994–256001 Rasmusson AG, Møller IM (1990) NADP-utilizing enzymes in the matrix of plant mitochondria. Plant Physiol 94: 1012–1018 Rasmusson AG, Møller IM (1991) Effects of calcium ions and inhibitors on internal NAD(P)H dehydrogenases in plant mitochondria. Eur J Biochem 202: 617–623[Web of Science][Medline] Rasmusson AG, Soole KL, Elthon TE (2004) Alternative NAD(P)H dehydrogenases of plant mitochondria. Annu Rev Plant Biol 55: 23–39[CrossRef][Medline] Rhoads DM, Subbaiah CC (2007) Mitochondrial retrograde regulation in plants. Mitochondrion 7: 177–194[CrossRef][Web of Science][Medline] Rikhvanov EG, Gamburg KZ, Varakina NN, Rusaleva TM, Fedoseeva IV, Tauson EL, Stupnikova IV, Stepanov AV, Borovskii GB, Voinikov VK (2007) Nuclear-mitochondrial cross-talk during heat shock in Arabidopsis cell culture. Plant J 52: 763–778[CrossRef][Web of Science][Medline] Roberts TH, Fredlund KM, Møller IM (1995) Direct evidence for the presence of two external NAD(P)H dehydrogenases coupled to the electron transport chain in plant mitochondria. FEBS Lett 373: 307–309[CrossRef][Web of Science][Medline] Roessner-Tunali U, Urbanczyk-Wochniak E, Czechowski T, Kolbe A, Willmitzer L, Fernie AR (2003) De novo amino acid biosynthesis in potato tubers is regulated by sucrose levels. Plant Physiol 133: 683–692 Sabar M, De Paepe R, de Kouchkovsky Y (2000) Complex I impairment, respiratory compensations, and photosynthetic decrease in nuclear and mitochondrial male sterile mutants of Nicotiana sylvestris. Plant Physiol 124: 1239–1250 Saisho D, Nambara E, Naito S, Tsutsumi N, Hirai A, Nakazono M (1997) Characterization of the gene family for alternative oxidase from Arabidopsis thaliana. Plant Mol Biol 35: 585–596[CrossRef][Web of Science][Medline] Singer TP (1979) Mitochondrial electron-transport inhibitors. Methods Enzymol 55: 454–462[Medline] Sweetlove L, Mowday B, Hebestreit HF, Leaver CJ, Millar AH (2001) Nucleoside disphosphate kinase III is localised to the inter-membrane space in plant mitochondria. FEBS Lett 508: 272–276[CrossRef][Web of Science][Medline] Sweetlove LJ, Heazlewood JL, Herald V, Holtzapffel R, Day DA, Leaver CJ, Millar AH (2002) The impact of oxidative stress on Arabidopsis mitochondria. Plant J 32: 891–904[CrossRef][Web of Science][Medline] Taylor NL, Heazlewood JL, Day DA, Millar AH (2005) Differential impact of environmental stresses on the pea mitochondrial proteome. Mol Cell Proteomics 4: 1122–1133 Taylor NL, Rudhe C, Hulett JM, Lithgow T, Glaser E, Day DA, Millar AH, Whelan J (2003) Environmental stresses inhibit and stimulate different protein import pathways in plant mitochondria. FEBS Lett 547: 125–130[CrossRef][Web of Science][Medline] Thimm O, Blaesing O, Gibon Y, Nagel A, Meyer S, Krüger P, Selbig J, Müller LA, Rhee SY, Stitt M (2004) MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J 37: 914–939[CrossRef][Web of Science][Medline] Vidal G, Ribas-Carbo M, Garmier M, Dubertret G, Rasmusson AG, Mathieu C, Foyer CH, De Paepe R (2007) Lack of respiratory chain complex I impairs alternative oxidase engagement and modulates redox signaling during elicitor-induced cell death in tobacco. Plant Cell 19: 640–655 Zerbetto E, Vergani L, Dabbeni-Sala F (1997) Quantification of muscle mitochondria oxidative phosphorylation enzymes via histochemical staining of blue native polyacrylamide gels. Electrophoresis 18: 2059–2064[CrossRef][Web of Science][Medline] This article has been cited by other articles:
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| ASPB Publications | PLANT PHYSIOLOGY® | THE PLANT CELL | |
|---|---|---|---|