TY - JOUR
AU - Possell, Malcolm
AB - Background and Aims Isoprene is the most important volatile organic compound emitted by land plants in terms of abundance and environmental effects. Controls on isoprene emission rates include light, temperature, water supply and CO2 concentration. A need to quantify these controls has long been recognized. There are already models that give realistic results, but they are complex, highly empirical and require separate responses to different drivers. This study sets out to find a simpler, unifying principle. Methods A simple model is presented based on the idea of balancing demands for reducing power (derived from photosynthetic electron transport) in primary metabolism versus the secondary pathway that leads to the synthesis of isoprene. This model's ability to account for key features in a variety of experimental data sets is assessed. Key results The model simultaneously predicts the fundamental responses observed in short-term experiments, namely: (1) the decoupling between carbon assimilation and isoprene emission; (2) a continued increase in isoprene emission with photosynthetically active radiation (PAR) at high PAR, after carbon assimilation has saturated; (3) a maximum of isoprene emission at low internal CO2 concentration (ci) and an asymptotic decline thereafter with increasing ci; (4) maintenance of high isoprene emissions when carbon assimilation is restricted by drought; and (5) a temperature optimum higher than that of photosynthesis, but lower than that of isoprene synthase activity. Conclusions A simple model was used to test the hypothesis that reducing power available to the synthesis pathway for isoprene varies according to the extent to which the needs of carbon assimilation are satisfied. Despite its simplicity the model explains much in terms of the observed response of isoprene to external drivers as well as the observed decoupling between carbon assimilation and isoprene emission. The concept has the potential to improve global-scale modelling of vegetation isoprene emission. Isoprene, modelling, electron transport, photosynthesis, temperature, carbon dioxide, isoprene emission, volatile organic compounds INTRODUCTION Isoprene (2-methyl-1,3-butadiene; C5H8) is a highly volatile and reactive unsaturated hydrocarbon that is produced continuously in daylight by many terrestrial plants, and in great abundance by broadleaved trees. On a mass basis, it is the most important biogenic volatile organic compound (BVOC) emitted by vegetation, with an annual global emission of approximately 0·5 × 1015 g C. This is similar in magnitude to the total annual emission of the greenhouse gas methane (CH4) from all natural sources combined (Guenther et al., 1995, 2006; Laothawornkitkul et al., 2009). Although not a greenhouse gas itself, isoprene reacts in the atmosphere with oxidants, including hydroxyl radicals (OH) and ozone (O3) (Fan and Zhang, 2004), and consequently influences the atmospheric lifetime and concentration of CH4 (Poisson et al., 2000; Collins et al., 2002, 2010; Pike and Young, 2009). The influence of isoprene on atmospheric oxidation capacity has been proposed as one of the controls of the glacial–interglacial variations of atmospheric CH4, as recorded in ice cores (Valdes et al., 2005; Singarayer et al., 2011). Isoprene also enhances the production of tropospheric ozone (O3), a potent greenhouse gas and toxic pollutant, under high-NOx conditions (Sanderson et al., 2003; Hauglustaine et al., 2005), and can significantly affect the atmosphere's radiative balance through the generation of secondary organic aerosols (Claeys et al., 2004; Heald et al., 2008; Carlton et al., 2009; Nozière et al., 2011). Isoprene emissions by plants at the leaf scale respond to changes in photosynthetically active radiation (PAR), temperature, ambient CO2 concentration and drought (Sharkey and Yeh, 2001; Laothawornkitkul et al., 2009; Pacifico et al., 2009; Niinemets, 2010). Despite general agreement between models under the present climate, simulations of future isoprene emissions, and their potential impact on atmospheric chemistry, change dramatically depending on the temperature and light responses of the model (Keenan et al., 2009) and whether the model includes a physiological response of isoprene emission to CO2 (Heald et al., 2009; Young et al., 2009; Pacifico et al., 2012). Given the continuously increasing atmospheric CO2 concentration and its impact on future temperature, we need to understand the processes behind observed responses, and use that understanding to build better models. The adaptive significance of isoprene emission is thought to be connected with enhancing membrane stability at high temperatures, and protection against oxidative stress – including that induced by high temperatures (Sharkey and Yeh, 2001; Vickers et al., 2009; Velikova et al., 2011, 2012). On time scales of weeks to years, acclimation mechanisms acting at the level of gene transcription may operate, possibly in such a way as to match isoprene synthase activity to adaptive requirements (Grote and Niinemets, 2008; Monson et al., 2012; Harrison et al., 2013). Here, however, we focus on the immediate responses of isoprene emission to environmental variations, as observed in experiments conducted over a time scale of minutes to hours, and the basic metabolic mechanisms that may be responsible for them. The biosynthesis of isoprene occurs via the chloroplastic methylerythritol 4-phosphate (MEP) pathway (Lichtenthaler, 1999; Logan et al., 2000; Sharkey et al., 2008). 13C labelling experiments have shown that the majority of the C in isoprene comes directly from photosynthesis, with the remainder coming from cytosolic C pools depending upon the environmental conditions (Delwiche and Sharkey, 1993; Kreuzwieser et al., 2002; Karl et al., 2002; Affek and Yakir, 2003; Loreto et al., 2004). However, the metabolic controls of the MEP pathway are only beginning to be elucidated (Li and Sharkey, 2012). With incomplete understanding of the metabolic controls of the pathway, models have been developed on the basis of experimental studies of the relationships between isoprene emission and environmental variables. The approach with the longest pedigree combines empirically derived functions for each environmental effect: this is the principle of the MEGAN model (Guenther et al., 2006), developed from the pioneer work of Guenther et al. (1993). Other approaches have made more direct use of the limited available information at the biochemical process level, e.g. SIM-BIM (Zimmer et al., 2000, 2003) and the models of Niinemets et al. (1999) and Martin et al. (2000). Aside from the model from Martin et al. (2000), which has an ATP limitation for isoprene production at high internal CO2 concentration (ci), all these models need an empirical parameterization to reproduce the observed CO2 response. This is potentially quite a severe limitation because there may be unforeseen interactions between the effects of different environmental drivers. Empirical models such as MEGAN include a multiplicity of functions for each environmental response of isoprene emission. More mechanistic approaches such as SIM-BIM, on the other hand, require information on many parameters. This might also be an issue because there is a generally accepted trade-off between the multiplicity of required parameter values and model robustness. We set out to identify a unifying principle that might transcend these limitations. Our starting point was the model of Niinemets et al. (1999), which is based on quantifying the NADPH requirement of isoprene synthesis. Niinemets et al. (1999) assumed that a certain fraction of the total electron flux generated by Photosystem II is allocated to this function. The model we present here, initially proposed in Harrison et al. (2013), builds on Niinemets' work but differs in one fundamental respect: it links isoprene emission to the electron availability for isoprene emission, relative to the needs of CO2 assimilation. Therefore, the model predicts higher isoprene emissions when absorbed radiant energy (leading to the ‘supply’ of NADPH) exceeds the ‘demand’ for CO2 assimilation. An excess of energy arises because of a mismatch between light availability and carboxylation capacity, which typically occurs daily – especially at high PAR, associated high temperature and under water stress. We compare the model's predictions of observed, published environmental responses of isoprene emission to changes in PAR and the leaf-internal concentration of CO2 (ci) with those obtained with the Guenther et al. (1993) algorithm, hereafter called G93, which is the basis of the widely used MEGAN model (Guenther et al., 2006), and with the model of Niinemets et al. (1999), hereafter called the Niinemets model. We also compare the theoretical temperature responses of our model with G93 and the Niinemets model. We focus on these two models as they have been widely used at the global scale (Guenther et al., 2006; Lathière et al., 2010; Arneth et al., 2011; Pacifico et al., 2012). However, other isoprene models have been developed. Reviews can be found in Arneth et al. (2007a), Grote and Niinemets (2008), and Monson et al. (2012). HYPOTHESIS In isoprene-emitting plants with the C3 pathway of photosynthesis, over 90 % of isoprene production takes place in the chloroplast via the MEP pathway (Lichtenthaler et al., 1997; Sharkey et al., 2008). The final stage is the enzymatic synthesis of isoprene from its precursor, dimethylallyl diphosphate (DMADP). On a per-molecule basis, isoprene synthesis is energetically expensive, and has a high requirement for reducing power (14 NADPH for one molecule of isoprene). For comparison, only six NADPH are needed to synthesize glyceradehyde 3-phosphate (G3P), and only five for pyruvate. NADPH consumption for G3P and pyruvate synthesis takes place within the Calvin cycle and therefore is linked to the electron cost for carbon assimilation. Three additional reducing steps are needed within the MEP pathway to reduce G3P and pyruvate to DMADP. These supplementary reducing steps consume one further NADPH, and two additional reducing equivalents in the form of either NADPH or ferredoxin (Fd) (Charon et al., 1999; Hecht et al., 2001; Seemann et al., 2006; Li and Sharkey, 2012). Our hypothesis focuses on these additional reduction steps, which are directly linked to the production of isoprene. Isoprene production is typically measured in nanomoles per second while photosynthesis and respiration rates (to which G3P and pyruvate production are linked) are measured in micromoles per second. Hence, the major consumption of reducing power takes place within the Calvin cycle and associated photorespiration while the diversion of reducing power to the MEP pathway is very small. Yet there is abundant circumstantial evidence for a link between the availability of reducing power (after the requirements of carbon assimilation have been accounted for) and the magnitude of this diversion. The MEP pathway is tightly linked to the photosynthetic apparatus, involves light-dependent reactions and takes place in the chloroplast. Higher isoprene emission capacity is encountered under conditions when photoinhibition occurs, including high light intensities, low ci and high temperatures. Isoprene emissions decrease if plants are fed with nitrate (note that nitrate reduction to ammonia occurs mainly in the cytoplasm and consumes NAPDH) instead of being fed with ammonia directly (Campbell, 1988; Rosenstiel et al., 2004). Li and Sharkey (2012) measured an extremely high level of the intermediate metabolite, methylerythritol cyclodiphosphate (MEcDP), in an N2 atmosphere, where the carbon assimilation and photorespiration sinks for NADPH are blocked. Thus, it might be that the MEP pathway acts as a ‘branch circuit’ with the amount of NAPDH allocated to it increasing in proportion to the amount of reducing power to spare from other functions. Thus, we hypothesize that isoprene emission is regulated in the short term by variations of the DMADP pool size, linked to the excess or deficit of electrons (and so also reducing power) relative to the needs of carbon assimilation. Figure 1 provides a schematic of the processes involved. When the chloroplast is illuminated, light absorbed by the thylakoids generates the electron flux (Jtot) that finally reduces NADP+ to NADPH. Most of the NADPH is used in the Calvin cycle for carbon fixation, but the total NADPH thus generated (≈ 0·5 Jtot) exceeds the amount consumed in the Calvin Cycle (≈ 0·5 JCO2+O2). When assimilation is light-limited (at high ci and/or low PAR) there is still some NADPH available for other functions, which include nitrate reduction (Canvin and Atkins, 1974; Niinemets, 2004; Eichelmann et al., 2011) and DMADP synthesis. When assimilation is Rubisco-limited (at low ci and/or high PAR) this excess of NADPH becomes larger, allowing more NADPH to be used in DMADP synthesis. Fig. 1. Open in new tabDownload slide Schematic of the processes underlying the model of isoprene emissions. The availability of reducing power (NADPH) for CO2 assimilation is represented by a colour scheme, from dark blue (deficit of NADPH) to dark red (excess of NADPH). Symbols: NADPH and NADP+, nicotinamide adenine dinucleotide phosphate; DMADP, dimethylallyl diphosphate; Pyr, pyruvate; G3P, glyceraldehyde 3-phosphate; DOXP, 1-deoxy-d-xylulose 5-phosphate; MEP, 2-C-methyl-d-erythritol 4-phosphate. Fig. 1. Open in new tabDownload slide Schematic of the processes underlying the model of isoprene emissions. The availability of reducing power (NADPH) for CO2 assimilation is represented by a colour scheme, from dark blue (deficit of NADPH) to dark red (excess of NADPH). Symbols: NADPH and NADP+, nicotinamide adenine dinucleotide phosphate; DMADP, dimethylallyl diphosphate; Pyr, pyruvate; G3P, glyceraldehyde 3-phosphate; DOXP, 1-deoxy-d-xylulose 5-phosphate; MEP, 2-C-methyl-d-erythritol 4-phosphate. This reasoning suggests the following simple model: (1) where Iso is the rate of isoprene emission; f(ci) is a function of internal CO2 concentration; J is an estimate of the total electron flux, taken to be a non-rectangular hyperbolic function of absorbed PAR and the maximum electron flux Jmax, following Farquhar et al. (1980); Jv is the electron flux required to support Rubisco-limited carbon assimilation; and a and b are parameters. The electron flux required to support carbon assimilation is derived as follows (from Farquhar et al., 1980): (2) where Aj is the gross (light-limited) assimilation rate and Γ* is the CO2 compensation point in the absence of dark respiration. Hence (3) When Rubisco limits photosynthesis, then (4) where Av is the gross (Rubisco-limited) assimilation rate, Vcmax is the Rubisco capacity and Km = Kc(1 + [O2]/Ko) where Kc and Ko are the Michaelis coefficients of Rubisco for CO2 and O2 respectively (Farquhar et al., 1980). Substituting this into eqn (3) gives: (5) It should be noted that J in eqn (1) is used as an estimate of Jtot and could be an underestimate (Singsaas et al., 2001; Niinemets, 2004). More details of the photosynthetic model, as used in this paper, can be found in the Supplementary Data. The term aJ in eqn (1) represents a ‘baseline’ of isoprene emission under light-limited conditions under the equilibrium conditions for carbon assimilation (J = Jv, energy supply = Rubisco demand), while b(J – Jv) represents variation in isoprene emission due to the disequilibrium between supply and demand. The function f(ci) in eqn (1) is chosen to take the value ci/Γ* when ci ≤ Γ* and 1 otherwise. Because of this function, the model differs slightly from the one we proposed in Harrison et al. (2013). The function f(ci) reflects the idea that a minimum rate of supply of carbon chains is required for isoprene synthesis, and the common observation that isoprene emission ceases abruptly when ci <Γ* (Wolfertz et al., 2003; Rasulov et al., 2009, 2011; Monson et al., 2012; Sun et al., 2012). This fall-off of isoprene at low ci is not always observed: emission of isoprene in CO2-free air has been reported in a few studies (Monson and Fall, 1989; Affek and Yakir, 2003; Li and Sharkey, 2012), but comparable conditions are not found in natural environments. Although based conceptually on the NADPH requirements of isoprene synthesis and the Farquhar photosynthesis model, our approach differs from that of Niinemets et al. (1999, 2004) in which isoprene production was assumed to be closely linked to the light-limited carbon assimilation rate (Aj). This difference has important consequences, as we will show. TESTS OF THE HYPOTHESIS We consider the observed environmental responses of isoprene emission (Iso) and also the ratio of isoprene emission to carbon gross assimilation (Iso/Agross), which is a sensitive indicator of the allocation of reducing power to the MEP pathway versus the Calvin cycle. We will also consider changes in the ratio of isoprene emission to carbon net assimilation (Iso/Anet). Responses to PAR Equation (1) predicts an increase of isoprene emission with PAR, but also an increase of the ratio Iso/Agross (Fig. 2A). The predicted behaviour of Iso/Anet (Anet = A – Rd, where Rd is mitochondrial respiration) is substantially different at low PAR, as shown in Fig. 2A. At saturating PAR, the difference becomes less important. This is due to the introduction of the Rd term, which affects the assimilation independently from the allocation of reducing power between carbon fixation and secondary metabolism. Most laboratory experiments have reported only Anet; this should be kept in mind while interpreting the results. Fig. 2. Open in new tabDownload slide Modelled responses of the normalized ratio of isoprene to CO2 assimilation to changes in PAR for (A) our model, (B) G93 and (C) the Niinemets model. T = 30 °C, ci = 273 µmol mol−1, Vcmax_25 °C = 70 µmol m−2 s−1, Jmax_25 °C = 130 µmol m−2 s−1 based on values from Arneth et al. (2007a). The solid line represents the ratio of isoprene emission to gross assimilation, the dashed line to net assimilation. Normalized ratio of isoprene to CO2 net assimilation was simulated for two extreme values of dark respiration to illustrate the potential effect of the magnitude of Rd on how Iso/Anet varies with PAR: grey short-dashed line, low Rd;Rd_25 °C = 0·5 µmol m−2 s−1; black long-dashed line, high Rd;Rd_25 °C = 2 µmol m−2 s−1. Isoprene model parameters a and b (eqn 1) are based on data of Possell and Hewitt (2011; fig. 6). Fig. 2. Open in new tabDownload slide Modelled responses of the normalized ratio of isoprene to CO2 assimilation to changes in PAR for (A) our model, (B) G93 and (C) the Niinemets model. T = 30 °C, ci = 273 µmol mol−1, Vcmax_25 °C = 70 µmol m−2 s−1, Jmax_25 °C = 130 µmol m−2 s−1 based on values from Arneth et al. (2007a). The solid line represents the ratio of isoprene emission to gross assimilation, the dashed line to net assimilation. Normalized ratio of isoprene to CO2 net assimilation was simulated for two extreme values of dark respiration to illustrate the potential effect of the magnitude of Rd on how Iso/Anet varies with PAR: grey short-dashed line, low Rd;Rd_25 °C = 0·5 µmol m−2 s−1; black long-dashed line, high Rd;Rd_25 °C = 2 µmol m−2 s−1. Isoprene model parameters a and b (eqn 1) are based on data of Possell and Hewitt (2011; fig. 6). The response of normalized Iso/Agross (and Iso/Anet) with PAR is predicted to take place in three stages (Fig. 2A). Stage 1: light-limited carbon assimilation This stage occurs when PAR absorbed is insufficient to generate an electron flux to satisfy Rubisco capacity. It is characterized by an initial steep increase of Iso/Agross with PAR, becoming gradually less steep at higher PAR. For Iso/Anet the form of the response at low PAR depends on the magnitude of Rd. Stage 2: transition between light- and Rubisco-limited carbon assimilation This stage is characterized by a discontinuity (abrupt increase) in the slope of Iso/Agross (and Iso/Anet) versus PAR. Stage 3: Rubisco-limited carbon assimilation When the electron requirement for carbon assimilation is fully satisfied, the additional reducing power generated by increasing PAR allows Iso/Agross to continue increasing while Agross remains constant. In this stage, Iso/Anet follows a similar pattern of the Iso/Agross and increases with PAR. With still further increases in PAR, Iso/Agross and Iso/Anet eventually saturate, as J tends to its maximal value (Jmax). Note that the PAR flux where the transition between light- and Rubisco-limited assimilation occurs (Stage 2) as well as the rate of increase of Iso/Agross with increasing PAR are dependent on both the photosynthetic and the isoprene model parameters. We also examined the normalized responses of Iso/Agross and Iso/Anet to changes in PAR in the G93 and Niinemets models (Fig. 2B, C). Under light-limited conditions (Stage 1), the picture differs dramatically depending on the model. In the Niinemets model, isoprene emissions are tightly linked to Aj (see Supplementary Data) and therefore the response of Iso to PAR necessarily has the same shape as that of Aj, irrespective of the chosen values of Vcmax or Jmax. As a result, the ratio Iso/Agross in this model is always constant under light-limited conditions, where carbon assimilation is equal to Aj. The ratio Iso/Anet, when simulated with the Niinemets model, always decreases with PAR under light-limited conditions. In G93, by contrast, changes in Iso/Agross with PAR are strongly dependent on Vcmax, Jmax and temperature under light-limited conditions. Consequently, the increase followed by a decrease of Iso/Agross under light-limited conditions, shown in Fig. 2B, is one of the possible responses of Iso/Agross for G93, obtained for the temperature and photosynthetic parameters chosen for this simulation. Changing those parameters changes the shape of the response, and Iso/Agross can decrease or increase at low PAR. Introducing a dark respiration term affects the shape of the response of Iso/Anet with PAR, as represented by the dashed lines in Fig. 2B. Hence, G93 can potentially show an Iso/A response to PAR similar to that of our model. All three models predict increasing Iso/Agross with PAR under Rubisco-limited conditions. Indeed, in the Niinemets model, as isoprene emissions are linked to Aj, they must continue to increase even when carbon assimilation is Rubisco-limited. In that sense, the Niinemets model implicitly allows consumption of extra NADPH above the needs for carbon assimilation (for the PAR response only). For G93, isoprene emission approaches an asymptotic value at high PAR, while the Farquhar model fully saturates under Rubisco-limited conditions at high PAR. Most studies reporting the fraction of assimilated carbon that is re-emitted as isoprene have found that it increases with PAR, in line with our predictions (Sharkey and Loreto, 1993; Harley et al., 1996; Lerdau and Keller, 1997; Niinemets et al., 2010). However, one study (Lerdau and Throop, 1999) found no significant increase in Iso/Anet with PAR for most of the tropical taxa they investigated. Figure 3A compares the relationships of Iso/Anet to PAR in our model and in digitized data from Sharkey and Loreto (1993) on kudzu leaves (Pueraria lobata). Assuming Jv is constant (no variation in ci; Fig. 3B), the observed isoprene emissions show a strong positive linear relationship with J (r2 = 0·97). The model parameters a and b (eqn 1) have been estimated from this linear regression. The Farquhar model parameters were estimated with a best data/model fit by minimizing the residual sum of squares (RSS). The comparison between our model and the data for Iso/Anet shows excellent agreement (r2 = 0·92). In comparison, G93 and the Niinemets model both show poor agreement (r2 = 0·19 and r2 = 0·06, respectively). Yet the three models show a good agreement of the modelled isoprene alone (Iso) with data (our model: r2 = 0·97; G93: r2 = 0·92; Niinemets: r2 = 0·97; results not shown). Fig. 3. Open in new tabDownload slide The relationship between isoprene emission and NADPH availability for carbon assimilation with changing PAR. (A) Increasing values of the isoprene to CO2 net assimilation ratio with increasing PAR, based on data digitized from figure 2 in Sharkey and Loreto (1993). Simulations made with our model, G93, and the Niinemets model are as indicated in the key. (B) The linear regression between isoprene data and the light-limited electron flux (J). Plant-specific isoprene parameters (a, b) are estimated from this linear regression and parameters for assimilation (Vcmax, Jmax) were fitted to the assimilation observations by minimizing the residual sum of squares (RSS). In both panels, the availability of reducing power (NADPH) for CO2 assimilation is illustrated by a colour scheme, from dark blue (deficit) to dark red (excess). Fig. 3. Open in new tabDownload slide The relationship between isoprene emission and NADPH availability for carbon assimilation with changing PAR. (A) Increasing values of the isoprene to CO2 net assimilation ratio with increasing PAR, based on data digitized from figure 2 in Sharkey and Loreto (1993). Simulations made with our model, G93, and the Niinemets model are as indicated in the key. (B) The linear regression between isoprene data and the light-limited electron flux (J). Plant-specific isoprene parameters (a, b) are estimated from this linear regression and parameters for assimilation (Vcmax, Jmax) were fitted to the assimilation observations by minimizing the residual sum of squares (RSS). In both panels, the availability of reducing power (NADPH) for CO2 assimilation is illustrated by a colour scheme, from dark blue (deficit) to dark red (excess). We also compiled data on the response of Iso/Anet to PAR from the limited number of published studies to assess the generality of the pattern (Fig. 4). The publications reported Anet rather than Agross, and did not typically provide measurements of Rd. As the predicted response of Iso/Anet for low PAR is dependent on Rd, it is not surprising to observe an initial decline in Iso/Anet with PAR for some of the 18 experiments. More importantly, the great majority of the studies show increasing Iso/Anet up to the highest PAR fluxes, especially when photosynthesis saturates (open circles in Fig. 4). In some studies a drop in assimilation rate at high PAR contributed to this increase in Iso/Anet at high PAR; this was probably due to stomatal closure at high PAR, resulting in reduced ci. Fig. 4. Open in new tabDownload slide Observed responses of the normalized ratio (isoprene emission/net CO2 assimilation) to changes in PAR. Light-limited and light-saturated carbon assimilation as indicated in the key. Data digitized from (A) Sharkey and Loreto (1993), (B, C) Loreto and Sharkey (1990), (D) Monson and Fall (1989), (E–N) Lerdau and Keller (1997), (O–R) Sun et al. (2012). (O) Plant grown at ambient CO2, chamber [CO2] = 380 µmol mol−1. (P) Plant grown at ambient CO2, chamber [CO2] = 780 µmol mol−1. (Q) Plant grown at elevated CO2, chamber [CO2] = 380 µmol mol−1. (R) Plant grown at elevated CO2, chamber [CO2] = 780 µmol mol−1. Fig. 4. Open in new tabDownload slide Observed responses of the normalized ratio (isoprene emission/net CO2 assimilation) to changes in PAR. Light-limited and light-saturated carbon assimilation as indicated in the key. Data digitized from (A) Sharkey and Loreto (1993), (B, C) Loreto and Sharkey (1990), (D) Monson and Fall (1989), (E–N) Lerdau and Keller (1997), (O–R) Sun et al. (2012). (O) Plant grown at ambient CO2, chamber [CO2] = 380 µmol mol−1. (P) Plant grown at ambient CO2, chamber [CO2] = 780 µmol mol−1. (Q) Plant grown at elevated CO2, chamber [CO2] = 380 µmol mol−1. (R) Plant grown at elevated CO2, chamber [CO2] = 780 µmol mol−1. As shown in Fig. 2, the Niinemets model cannot reproduce the positive response of Iso/Anet to PAR that is generally observed under low PAR. Our model, along with G93, fully captures the shape of the observed response of Iso/Anet to PAR over the full range of PAR. But our model also provides a process-based explanation for this response. We also examined the relationship between Iso/Agross and PAR at the canopy scale, at which isoprene emission is more likely to be controlled by the DMADP pool size than by isoprene synthase activity (Vickers et al., 2010). We used simultaneous CO2 and isoprene flux measurements made at Harvard Forest, Massachusetts, USA (42·54°N, 72·17°W) (Urbanski et al., 2007; McKinney et al., 2011). Data used were obtained during the 2007 growing seasons using eddy covariance, with proton transfer reaction mass spectrometry used to measure the isoprene mixing ratio (McKinney et al., 2011). Daytime data were selected for temperatures above 23 °C where variations in isoprene emission were no longer significantly driven by temperature (Fig. A1). Ecosystem respiration, estimated from night-time CO2 flux measurements, was used to convert the daytime measured net ecosystem CO2 exchanges into canopy-scale gross assimilation rates. Canopy-scale carbon assimilation shows a typical rectangular hyperbolic response to PAR, but the response of isoprene emission to PAR is closer to linearity, and emissions do not saturate at high PAR (Fig. 5A, B). Thus, above a PAR threshold of approx. 300 µmol m−2 s−1, Iso/Agross increases with PAR even at high PAR, where assimilation is light-saturated, consistently with our hypothesis. Fig. 5. Open in new tabDownload slide Above-canopy gross assimilation (A), isoprene emissions (B) and isoprene emission/gross assimilation (C), in relation to PAR. From flux measurements at Harvard Forest. Fig. 5. Open in new tabDownload slide Above-canopy gross assimilation (A), isoprene emissions (B) and isoprene emission/gross assimilation (C), in relation to PAR. From flux measurements at Harvard Forest. Scaling from leaf to canopy involves additional processes, such as within-canopy chemistry and canopy structure effects (Grote, 2007; Keenan et al., 2011; Bryan et al., 2012). Therefore, a canopy model is needed to fully account for these results, especially for low PAR where deposition processes can influence the observed above-canopy isoprene emissions and possibly explain the observed drop in Iso/Agross. Nevertheless these results, along with those of laboratory experiments, corroborate the notion that isoprene emission is related to the availability of electrons generated in photosynthesis, relative to the demand for them to be used in carbon assimilation. Responses to ci Responses of isoprene emission to ambient CO2 concentration have been widely reported. Plants grown at high atmospheric CO2 concentrations generally emit less isoprene than those grown at lower CO2 concentrations. On short time scales, isoprene emission has also been shown to respond strongly and rapidly to ci, with lower emission rates at higher ci (Rosenstiel et al., 2003; Wilkinson et al., 2009; Possell and Hewitt, 2011; Sun et al., 2012). The fact that rapid changes in ci evoke instantaneous responses in isoprene emission suggests that the driving mechanism must be tightly linked to processes in the chloroplast. The mechanisms behind the decoupling between isoprene emission and carbon assimilation in the response to ci are not well established. Niinemets et al. (1999) hypothesized that the dependency of isoprene emission on ci might be due to the partitioning of reducing power and ATP into the MEP pathway. However, the model of Niinemets et al. (1999) does not allow for any greater partitioning of reducing power to the MEP pathway at low ci. Isotopic labelling studies have provided evidence for the existence of extra-chloroplastic sources of carbon to support isoprene production. Hence, competition for phosphoenolpyruvate (PEP) between cytosolic and chloroplastic processes has been proposed as an explanation for the drop in isoprene emission at high ci due to the CO2-dependence of PEP carboxylase activity (Karl et al., 2002; Rosenstiel et al., 2003; Possell and Hewitt, 2011; Trowbridge et al., 2012). But these experiments compared plants grown at different CO2 concentrations. Gene expression involving changes in enzyme quantities cannot explain the observed fast (about 10-min) responses to changes in ci. We focus here only on the short-term responses to ci, in which isoprene emission appears to be tightly coupled to changes in the pool size of DMADP (Rasulov et al., 2009). Specifically, we examine whether the fast responses of isoprene emission to ci could be explained in a simple way by our model, based on the same mechanisms we have proposed to explain the response to PAR. At low ci, carbon fixation is Rubisco-limited, resulting in an excess of NADPH (Figs 1 and 6A). The excess of NADPH can be smaller or larger depending on PAR. This provides a simple explanation for why isoprene responses to changes in ci are light-dependent (Loreto and Sharkey, 1993; Fig. 7D). Moreover, our model can indeed reproduce the isoprene emission response to changes in ci. This is shown in Fig. 6A using data on Acacia nigrescens from Possell and Hewitt (2011). Here, isoprene emission shows a strong negative linear relationship with the Rubisco-limited electron flux, Jv (r2 = 0·70), as shown in Fig. 6B. The parameters a and b (eqn 1) of our model were estimated from this linear regression. When plotted against ci, our model shows a good agreement with the data (r2 = 0·70). Figure 6A also shows the response of the G93 and the Niinemets model with and without a CO2 inhibition effect (Arneth et al., 2007a; Pacifico et al., 2011). It is clear from Fig. 6A that these models do not reproduce the observed response of isoprene emission to ci. Without an additional empirical function for CO2 inhibition, isoprene emissions simulated with the Niinemets model are quite out of range. Instead, the model shows a strong negative correlation with the data (r2 = 0·7). The negative relationship can be explained by the fact that although the PAR and therefore the light-limited electron flux (J) are constant, light-limited assimilation (Aj) is strongly ci-dependent. Adding an empirical function to represent the CO2 inhibition effect, as in Arneth et al. (2007a), changes the shape of the response (allowing a decrease at high ci) but still the simulated emissions agree poorly with the data. Fig. 6. Open in new tabDownload slide The relationship between isoprene emission and NADPH availability for carbon assimilation with changing internal CO2 concentration ci. (A) Decreasing isoprene emissions with increasing leaf-internal CO2 concentration, ci (data from Possell and Hewitt, 2011); T = 30 °C, PAR = 1000 µmol m−2 s−1. Simulations made with our model, G93 and the Niinemets model are as indicated in the key. The dashed black line represent the Niinemets model with an additional CO2 effect represented by f(ci)= ci [ca = 390 µmol mol−1]/ci, where ca is the atmospheric CO2 concentration. The plain grey line represent the Niinemets model with an alternative additional CO2 effect represented by f(ca)= [ca = 390 µmol mol−1]/ca. The terms f(ci) and f(ca) are adapted from Arneth et al. (2007a). Standard isoprene emission factor (Is) is taken as the observed emission at ca = 390 µmol mol−1. (B) The linear regression between isoprene data and the electron flux required for carbon assimilation by Rubisco (Jv). Plant-specific isoprene parameters (a, b) are estimated from this linear regression and parameters for assimilation (Vcmax, Jmax) were fitted to the observations by minimizing the residual sum of squares (RSS). In both panels, the availability of reducing power (NADPH) for CO2 assimilation is represented by a colour scheme, from dark blue (deficit) to dark red (excess). Fig. 6. Open in new tabDownload slide The relationship between isoprene emission and NADPH availability for carbon assimilation with changing internal CO2 concentration ci. (A) Decreasing isoprene emissions with increasing leaf-internal CO2 concentration, ci (data from Possell and Hewitt, 2011); T = 30 °C, PAR = 1000 µmol m−2 s−1. Simulations made with our model, G93 and the Niinemets model are as indicated in the key. The dashed black line represent the Niinemets model with an additional CO2 effect represented by f(ci)= ci [ca = 390 µmol mol−1]/ci, where ca is the atmospheric CO2 concentration. The plain grey line represent the Niinemets model with an alternative additional CO2 effect represented by f(ca)= [ca = 390 µmol mol−1]/ca. The terms f(ci) and f(ca) are adapted from Arneth et al. (2007a). Standard isoprene emission factor (Is) is taken as the observed emission at ca = 390 µmol mol−1. (B) The linear regression between isoprene data and the electron flux required for carbon assimilation by Rubisco (Jv). Plant-specific isoprene parameters (a, b) are estimated from this linear regression and parameters for assimilation (Vcmax, Jmax) were fitted to the observations by minimizing the residual sum of squares (RSS). In both panels, the availability of reducing power (NADPH) for CO2 assimilation is represented by a colour scheme, from dark blue (deficit) to dark red (excess). Fig. 7. Open in new tabDownload slide Observed changes in the ratio of isoprene emission to net carbon assimilation with changes in (A) leaf-internal CO2 concentration (ci), (B) electron excess (J – Jv) (data from Possell and Hewitt, 2011). Observed changes with ci of (C) the ratio of isoprene emission to carbon assimilation and (D) isoprene emission, for three PAR fluxes (μmol m−2 s−1). Data digitized from Loreto and Sharkey (1990). Fig. 7. Open in new tabDownload slide Observed changes in the ratio of isoprene emission to net carbon assimilation with changes in (A) leaf-internal CO2 concentration (ci), (B) electron excess (J – Jv) (data from Possell and Hewitt, 2011). Observed changes with ci of (C) the ratio of isoprene emission to carbon assimilation and (D) isoprene emission, for three PAR fluxes (μmol m−2 s−1). Data digitized from Loreto and Sharkey (1990). Again using the data from Possell and Hewitt (2011), we plotted Iso/Anet versus ci (Fig. 7A) and J – Jv (Fig. 7B). These plots confirm that the fraction of assimilated carbon allocated to isoprene production increases under conditions of NADPH excess. This provides a simple explanation for the response of isoprene emission to ci. The extremely steep rise in Iso/Anet when J – Jv becomes positive is due to the combination of steeply increasing isoprene emission with decreasing assimilation rate as ci declines. Loreto and Sharkey (1990) measured changes in isoprene emission with changing ci at different PAR fluxes in Quercus rubra. Both Iso/Anet and isoprene emission are shown (Fig. 7C, D) to increase with PAR, consistent with a dependence on NADPH availability, at all values of ci. We compared the responses of G93 and the Niinemets model to ci at different PAR fluxes, together with our model (Fig. 8). Note that both G93 and the Niinemets model are applied here in their original formulations (see Supplementary Data for details), and therefore do not include additional parameterizations of the CO2 effect. A number of studies have used these same models with additional empirical functions, introduced specifically to account for the observed CO2 inhibition (Arneth et al., 2007b; Heald et al., 2009; Pacifico et al., 2011). G93 in its original formulation simulates no change at all in isoprene emission with changes in ci, although it has isoprene emission depending on PAR (Fig. 8B). The Niinemets model in its original formulation also simulates increasing isoprene emission with PAR, but here the modelled emissions increase with increasing ci, due to the fact that this model tightly links isoprene emission to light-limited assimilation (Fig. 8C). Thus, additional functions are required in both models to account for the observed effects of varying ci. In contrast, our model (Fig. 8A) can reproduce the form of the ci response shown in the data (Fig. 7D), as well as the effects of combined changes in ci and PAR (Fig. 7D), without the need for any additional function. Fig. 8. Open in new tabDownload slide Modelled responses of isoprene emission versus ci for three PAR fluxes (80, 180 and 700 µmol m−2 s−1). (A) our model, (B) G93 and (C) the Niinemets model. Parameters values as in Fig. 2. Emissions were normalized to a standard emission rate at T = 30 °C, ci = 273 µmol mol−1, PAR = 1000 µmol m−2 s−1. Fig. 8. Open in new tabDownload slide Modelled responses of isoprene emission versus ci for three PAR fluxes (80, 180 and 700 µmol m−2 s−1). (A) our model, (B) G93 and (C) the Niinemets model. Parameters values as in Fig. 2. Emissions were normalized to a standard emission rate at T = 30 °C, ci = 273 µmol mol−1, PAR = 1000 µmol m−2 s−1. Responses to leaf temperature The temperature dependence of isoprene emission differs from that of photosynthesis. Temperature optima for carbon assimilation are usually ≤ 30 °C in C3 plants, while isoprene emission peaks at ≈ 40 °C (Guenther et al., 1993; Niinemets et al., 1999; Sharkey and Yeh, 2001; Pacifico et al., 2009). An increase of Iso/A with temperature is usually observed (Sharkey and Loreto, 1993; Harley et al., 1996; Sharkey et al., 1996; Niinemets et al., 1999; Sharkey and Yeh, 2001). The optimum for isoprene emissions rarely exceeds 40 °C, so the temperature dependence of isoprene emission cannot be fully explained by the temperature dependence of isoprene synthase, which is maximally active between 45 and 48 °C (Monson et al., 1992; Lehning et al., 1999; Niinemets et al., 1999; Rasulov et al., 2010). The decrease in isoprene emissions above 40 °C has long been considered to be linked to the behaviour of the photosynthetic electron transport rate (Guenther et al., 1991; Niinemets et al., 1999). Rasulov et al. (2010) found that this decrease is associated with decline in the DMADP pool size and the energetic status of the leaf. In G93 the temperature dependency of isoprene emission is fixed with a temperature optimum around 38 °C. In the Niinemets model it is assumed to be primarily controlled by IspS activity, with the fraction of electrons used for isoprene production exponentially increasing with temperature. The temperature optimum for isoprene emissions in the Niinemets model is thus close to the optimum for IspS. Some global-scale studies have set an upper limit for the increase of ɛ with temperature, thereby reducing the temperature optimum to a value closer to 38 °C (Pacifico et al., 2011) (Supplementary Data A.2). Our model is based on the hypothesis that the production of DMADP depends on photosynthetic electron flux and variations in electron availability for functions other than carbon assimilation. Thus, our modelled optima for isoprene emissions are primarily driven by the behaviour of the light-limited electron flux. Figures 9 and 10 illustrate how a temperature response arises in our model. Carbon assimilation follows the lower of the temperature response curves of the Rubisco and light-limited assimilation rates. Rubisco activity usually has a higher temperature optimum than electron transport (Crafts-Brandner and Salvucci, 2000; Medlyn et al., 2002; Cen and Sage, 2005; Kattge and Knorr, 2007). At high PAR an excess of NADPH can arise for temperatures below the optimum for electron transport (J), so isoprene emissions increase. Above this optimum (Topt_J), Jv still increases even if assimilation is reduced, due to the higher affinity of Rubisco for O2 at high temperatures. Both J and (J − Jv) decrease for temperatures higher than Topt_J (as illustrated in Fig. 9B). Our model thereby predicts a temperature optimum of isoprene emissions that is closer to the temperature optimum of the electron transport rate. Beyond this optimum, our model predicts a drop in the availability of reducing power, leading to a decrease of DMADP and consequently isoprene emissions. At low PAR (light-limited condition), however, (J − Jv) decreases with increasing temperature, compensating for the increase of J. Predicted emissions are thus almost insensitive to temperature or even decrease with temperature (Fig. 10A, B). This behaviour is not realistic, so the model may be overestimating the effect of (J − Jv ) at low PAR. Fig. 9. Open in new tabDownload slide Explanation of the predicted temperature dependency of isoprene emissions. (A) The responses to temperature of the light-limited Aj, the Rubisco-limited Av and the gross assimilation Agross (as indicated in key); (B) the associated electrons fluxes (left-hand axis) and the associated electron availability (J – Jv) (right hand axis). (C) The normalized responses to temperature of our model, normalized IspS activity and the resulting product (see key). Temperature dependency of IspS is as described in Niinemets et al. (1999). Parameters values of the model are taken as in Fig. 2. Fig. 9. Open in new tabDownload slide Explanation of the predicted temperature dependency of isoprene emissions. (A) The responses to temperature of the light-limited Aj, the Rubisco-limited Av and the gross assimilation Agross (as indicated in key); (B) the associated electrons fluxes (left-hand axis) and the associated electron availability (J – Jv) (right hand axis). (C) The normalized responses to temperature of our model, normalized IspS activity and the resulting product (see key). Temperature dependency of IspS is as described in Niinemets et al. (1999). Parameters values of the model are taken as in Fig. 2. Fig. 10. Open in new tabDownload slide Responses to variation in temperature (°C) and PAR (μmol m−2 s−1) of electron availability (J – Jv) (μmol m−2 s−1) (A), our model (B) and our model simulations multiplied by a normalized function of enzymatic activity (C). Farquhar model parameters are for Quercus robur, as described in Medlyn et al. (2002). Isoprene model parameters a and b (eqn 1) are based on data of Possell and Hewitt (2011) (fig. 6). Model outputs in A and C are normalized to be unity at T = 30 °C and PAR = 1000 µmol m−2 s−1. Fig. 10. Open in new tabDownload slide Responses to variation in temperature (°C) and PAR (μmol m−2 s−1) of electron availability (J – Jv) (μmol m−2 s−1) (A), our model (B) and our model simulations multiplied by a normalized function of enzymatic activity (C). Farquhar model parameters are for Quercus robur, as described in Medlyn et al. (2002). Isoprene model parameters a and b (eqn 1) are based on data of Possell and Hewitt (2011) (fig. 6). Model outputs in A and C are normalized to be unity at T = 30 °C and PAR = 1000 µmol m−2 s−1. We infer that energetic control alone is insufficient to fully explain the observed temperature dependency of isoprene emission. In principle the activities of enzymes along the MEP pathway should also influence the production rate of DMADP, but very little is known about their temperature responses (Zimmer et al., 2000). Temperature optima for isoprene production are shifted toward higher temperature than Topt_J, probably because a decrease in DMADP pool size is compensated for by an increase in IspS activity (Rasulov et al., 2010). Taking into account the temperature response of IspS, we can reproduce this shift (Figs 9C and 10C). So we suggest that temperature effects on enzyme activity may need to be considered, as well as temperature effects on electron availability. A further limitation of our model is the paucity of available information on the temperature responses of Jmax and Vcmax (Medlyn et al., 2002; Kattge and Knorr, 2007). The experiments needed to quantify these responses are time-consuming, and in particular, few studies have included temperatures >40 °C. In general we would expect a decline in DMADP production for temperatures >40 °C due to thylakoid damage, while at temperatures above 45–48 °C irreversible damage to enzyme function will cause isoprene emission to cease. Using data from Medlyn et al. (2000) and references therein, we checked variations with temperature of electron availability among isoprene emitting species at 1000 µmol m−2 s−1 PAR (Fig. 11). We also tested the influence of the temperature response parameterization of Vcmax by contrasting an Arrhenius function with a peak function (Supplementary Data), as described in Medlyn et al. (2000). The temperature optima for the selected species are all higher for Vcmax than Jmax (Medlyn et al., 2002; Kattge and Knorr, 2007). Consequently, we predicted a decline in DMADP pool size above Topt_J, due to the decline of J being accompanied by a decline in (J − Jv), but the shape of the decline depended on the parameterization adopted. Fig. 11. Open in new tabDownload slide Left: temperature responses for different species of the light-limited electron flux (J) (dark grey dotted line), the Rubisco-limited electron flux (Jv) using an Arrhenius function for Vcmax (Jv Arrhenius, black dashed line), and the Rubisco-limited electron flux using a peak function for Vcmax (Jv peak, red dashed line). Right: the resulting temperature responses of (J – Jv), using an Arrhenius function for Vcmax (black solid line), and a peak function for Vcmax (red solid line). Farquhar parameters and calculation of Vcmax are as described in Medlyn et al. (2000). For Quercus robur: GH, greenhouse experiment; ME, mini-ecosystem experiment (Medlyn et al., 2000). Simulations are done for ci = 273 µmol mol−1 and PAR = 1000 µmol m−2 s−1. Fig. 11. Open in new tabDownload slide Left: temperature responses for different species of the light-limited electron flux (J) (dark grey dotted line), the Rubisco-limited electron flux (Jv) using an Arrhenius function for Vcmax (Jv Arrhenius, black dashed line), and the Rubisco-limited electron flux using a peak function for Vcmax (Jv peak, red dashed line). Right: the resulting temperature responses of (J – Jv), using an Arrhenius function for Vcmax (black solid line), and a peak function for Vcmax (red solid line). Farquhar parameters and calculation of Vcmax are as described in Medlyn et al. (2000). For Quercus robur: GH, greenhouse experiment; ME, mini-ecosystem experiment (Medlyn et al., 2000). Simulations are done for ci = 273 µmol mol−1 and PAR = 1000 µmol m−2 s−1. DISCUSSION We have used a conceptual model to ask whether variation in the availability of NADPH in the chloroplast can plausibly account for observed changes in isoprene emission with PAR, ci and leaf temperature. The answer is yes. By modelling isoprene emission as proportional to a simple metric of the excess or deficit of electrons relative to the demands of carbon assimilation, we have provided a unifying explanation for the lack of close coupling of isoprene emission with carbon assimilation, the disparities in carbon allocated to isoprene production, high isoprene emissions at low ci and the shift of the temperature optimum for isoprene emission above that of carbon assimilation but below that of isoprene synthase. To our knowledge, this is the first study that has attempted to model the flux of reducing power into the MEP production pathway based on the idea of a balance between electron supply and demand. Our hypothesis invokes mechanisms that are incompletely understood and thus is to some extent speculative. Nevertheless, it appears to have significant predictive power in explaining the already documented responses of isoprene emission to PAR, ci and (with some caveats) temperature. Moreover, this hypothesis provides a parsimonious explanation for the response of isoprene emission to drought. Under moderate to mild drought where the photosynthetic apparatus is not damaged (Cornic and Briantais, 1991), carbon assimilation is first reduced by stomatal closure (and thus reduced ci). Under higher drought severity, this reduction is greatly increased by decreased ATP in water-deficient leaves (Lawlor and Tezara, 2009), which reduces the photosynthetic metabolic potential (Apot), even if ci increases due to light respiration. The resulting oversupply of reducing power ensures that isoprene emissions continue at a high rate, although carbon assimilation is reduced (Niinemets, 2010). However, a decrease in ATP could also reduce isoprene emissions. Under extreme drought, however, damage to the photosynthesis apparatus eventually results in the cessation of both carbon assimilation and isoprene emission. A strong diurnal cycle is observed in isoprene emission at canopy scales. Low emissions during early morning and late afternoon contrast with high emissions during the midday period (Hewitt et al., 2011; Keenan and Niinemets, 2012). Our hypothesis explains this as a consequence of higher PAR and temperature, and lower ci due to partial stomatal closure associated with higher evaporative demand in the midday period. This simple explanation does not require the intervention of a circadian clock, as had been proposed by Hewitt et al. (2011). Note that we are not advocating a function of isoprene emissions as an ‘electron sink’ as was earlier proposed (e.g. Logan et al., 2000). It is clear from the findings of Li and Sharkey (2012) that the quantity of electrons used in isoprene synthesis is far too small for this function to be plausible. The low affinity of IspS for DMADP already argues strongly against this notion (Silver and Fall, 1995; Schnitzler et al., 1997). Our model implies that the allocation of reducing power to this pathway occurs under those conditions when electron availability is in excess, which fortuitously occurs during stress events when isoprene biosynthesis and emission is advantageous to the plant (e.g. Sharkey et al., 2001; Vickers et al., 2009). Our results provide an alternative, robust approach to modelling isoprene emissions for global change applications. But more work is needed before implementing the model in a global context. Particular attention should be given to the influence of enzymatic activity on temperature responses of the modelled rates of isoprene. The values of the parameters a and b (eqn 1), and their potential species and environmental dependencies, also call for further investigation at several scales: For leaves, by setting up experiments that could test interactions among the short-term responses of isoprene emission to different environmental drivers, and associated variations in the excess of electrons (i.e. isoprene/assimilation responses to ci at different PAR fluxes, together with isoprene/assimilation response to PAR at different ci); and the influence of growth conditions on the parameters. Note that, as most of the process-based models are closely linked to photosynthesis models, information on the values of Vcmax and Jmax associated with the isoprene standard emission rate would be valuable. For ecosystems, by upscaling the model from the leaf scale to the canopy, with particular attention to the response of Iso/Agross. This step would require a representation of the canopy structure and vertical mixing as well as the canopy chemistry accounting for isoprene oxidation, deposition and OH regeneration. At the global scale, with the possibility of using remotely sensed formaldehyde column concentrations to better constrain model parameters for different plant function types and environments. Formaldehyde is a product of isoprene oxidation. As it is observed by satellite, with global coverage, numerous studies have used formaldehyde data to investigate isoprene emission at larger scales (Palmer et al., 2003, 2006; Barkley et al., 2008; Stavrakou et al., 2009; Fortems-Cheiney et al., 2012). A comprehensive approach to isoprene modelling would also have to account for longer-term acclimation over a time scale of weeks to months, including responses to antecedent temperatures (Guenther et al., 2006), phenological stages, and differences between the short-term and acclimated responses to CO2 (Sun et al., 2012), which are presumably mediated by transcriptional control of the MEP pathway enzymes. It would be of particular interest to examine whether these acclimatory changes in isoprene emission are correlated with acclimatory changes in reducing power. However, some acclimatory shifts are unlikely to be explained by a model based on reducing power alone. For example, growth at higher temperatures leads to increased emission rates measured at a common temperature (Pétron et al., 2001; Niinemets et al., 2010), whereas reducing power at a given temperature tends to be reduced by high growth temperatures due to a decline in the Jmax/Vcmax ratio (Hikosaka et al., 1999; Onoda et al., 2005). Longer-term responses of isoprene emission to changes in growth temperature are therefore presumably governed by other factors, including transcriptional control of the MEP pathway enzymes. Conclusions A simple model of the biochemistry and physiology of isoprene emissions has been developed and used to test the hypothesis that the reducing power available to the synthesis pathway for isoprene varies according to demands of carbon assimilation. The model explains the observed response of isoprene production to environment and the coupling/decoupling between carbon assimilation and isoprene emission. The model has the potential to improve global-scale modelling of vegetation isoprene emissions, as well as emissions of isoprenoids that do not originate from storages. ACKNOWLEDGEMENTS We thank Karena McKinney for providing the original isoprene data for the Harvard forest site. We thank Russell Monson and Rüdiger Grote for their helpful and constructive comments on the manuscript. C.M. and I.C.P. have received funding from the European Community's Seventh Framework Programme (FP7 2007–2013) under grant agreement no. 238366. LITERATURE CITED Affek HP Yakir D Natural abundance carbon isotope composition of isoprene reflects incomplete coupling between isoprene synthesis and photosynthetic carbon flow Plant Physiology 2003 131 1727 1736 Google Scholar Crossref Search ADS PubMed WorldCat Arneth A Niinemets Ü Pressley S et al. Process-based estimates of terrestrial ecosystem isoprene emissions: incorporating the effects of a direct CO2-isoprene interaction Atmospheric Chemistry and Physics 2007a 7 31 53 Google Scholar Crossref Search ADS WorldCat Arneth A Miller PA Scholze M et al. CO2 inhibition of global terrestrial isoprene emissions: potential implications for atmospheric chemistry Geophysical Research Letters 2007b 34 L18813 Google Scholar Crossref Search ADS WorldCat Arneth A Schurgers G Lathiere J et al. Global terrestrial isoprene emission models: sensitivity to variability in climate and vegetation Atmospheric Chemistry and Physics 2011 11 8037 8052 Google Scholar Crossref Search ADS WorldCat Barkley MP Palmer PI Kuhn U et al. Net ecosystem fluxes of isoprene over tropical South America inferred from Global Ozone Monitoring Experiment (GOME) observations of HCHO columns Journal of Geophysical Research 2008 113 D20304 Google Scholar OpenURL Placeholder Text WorldCat Bryan AM Bertman SB Carroll MA et al. In-canopy gas-phase chemistry during CABINEX 2009: sensitivity of a 1-D canopy model to vertical mixing and isoprene chemistry Atmospheric Chemistry and Physics 2012 12 8829 8849 Google Scholar Crossref Search ADS WorldCat Campbell WH Nitrate reductase and its role in nitrate assimilation in plants Physiologia Plantarum 1988 74 214 219 Google Scholar Crossref Search ADS WorldCat Canvin DT Atkins CA Nitrate, nitrite and ammonia assimilation by leaves?: effect of light, carbon dioxide and oxygen Planta 1974 116 207 224 Google Scholar Crossref Search ADS PubMed WorldCat Carlton AG Wiedinmyer C Kroll JH A review of secondary organic aerosol (SOA) formation from isoprene Atmospheric Chemistry and Physics 2009 9 4986 5005 Google Scholar Crossref Search ADS WorldCat Cen Y Sage RF The regulation of Rubisco activity in response to variation in temperature and atmospheric CO2 partial pressure in sweet potato Plant Physiology 2005 139 979 990 Google Scholar Crossref Search ADS PubMed WorldCat Charon L Pale-grosdemange C Rohmer M On the reduction steps in the mevalonate independent 2-C-methyi-dhritol 4-phosphate (MEP) pathway for isoprenoid biosynthesis in the bacterium Zymomonas mobilis Tetrahedron Letters 1999 40 7231 7234 Google Scholar Crossref Search ADS WorldCat Claeys M Graham B Vas G et al. Formation of secondary organic aerosols through photooxidation of isoprene Science 2004 303 1173 1176 Google Scholar Crossref Search ADS PubMed WorldCat Crafts-Brandner SJ Salvucci ME Rubisco activase constrains the photosynthetic potential of leaves at high temperature and CO2 Proceedings of the National Academy of Sciences of the United States of America 2000 97 13430 13435 Google Scholar Crossref Search ADS PubMed WorldCat Collins WJ Derwent RG Johnson CE Stevenson DS The oxidation of organic compounds in the troposphere and their global warming potentials Climatic Change 2002 52 453 479 Google Scholar Crossref Search ADS WorldCat Collins WJ Sitch S Boucher O How vegetation impacts affect climate metrics for ozone precursors Journal of Geophysical Research 2010 115 D23308 Google Scholar OpenURL Placeholder Text WorldCat Cornic G Briantais J-M Partitioning of photosynthetic electron flow between CO2 and O2 reduction in a C3 leaf (Phaseolus vulgaris L.) at different CO2 concentrations and during drought stress Planta 1991 183 178 184 Google Scholar Crossref Search ADS PubMed WorldCat Delwiche CF Sharkey TD Rapid appearance of 13C in biogenic isoprene when 13CO2 is fed to intact leaves Plant, Cell and Environment 1993 16 587 591 Google Scholar Crossref Search ADS WorldCat Eichelmann H Oja V Peterson RB Laisk A The rate of nitrite reduction in leaves as indicated by O2 and CO2 exchange during photosynthesis Journal of Experimental Botany 2011 62 2205 2215 Google Scholar Crossref Search ADS PubMed WorldCat Fan J Zhang R Atmospheric oxidation mechanism of isoprene Environmental Chemistry 2004 1 140 149 Google Scholar Crossref Search ADS WorldCat Farquhar GD Von Caemmerer S Berry JA A biochemical model of photosynthetic CO2 assimilation in leaves of C3 species Planta 1980 149 78 90 Google Scholar Crossref Search ADS PubMed WorldCat Fortems-Cheiney A Chevallier F Pison I et al. The formaldehyde budget as seen by a global-scale multi-constraint and multi-species inversion system Atmospheric Chemistry and Physics 2012 12 6699 6721 Google Scholar Crossref Search ADS WorldCat Grote R Sensitivity of volatile monoterpene emission to changes in canopy structure: a model-based exercise with a process-based emission model The New Phytologist 2007 173 550 561 Google Scholar Crossref Search ADS PubMed WorldCat Grote R Niinemets Ü Modeling volatile isoprenoid emissions–a story with split ends Plant Biology 2008 10 8 28 Google Scholar Crossref Search ADS PubMed WorldCat Guenther AB Monson RK Fall R Isoprene and monoterpene emission rate variability: observations with eucalyptus and emission rate algorithm development Journal of Geophysical Research 1991 96 10799 10808 Google Scholar Crossref Search ADS WorldCat Guenther AB Zimmerman PR Harley PC Monson RK Fall R Isoprene and monoterpene emission rate variability: model evaluations and sensitivity analyses Journal of Geophysical Research 1993 98 12609 12617 Google Scholar Crossref Search ADS WorldCat Guenther AB Hewitt CN Erickson D et al. A global model of natural volatile organic compound emissions Journal of Geophysical Research 1995 100 8873 8892 Google Scholar Crossref Search ADS WorldCat Guenther AB Karl T Harley P Wiedinmyer C Palmer PI Geron C Estimates of global terrestrial isoprene emissions using MEGAN (Model of Emissions of Gases and Aerosols from Nature) Atmospheric Chemistry and Physics 2006 6 3181 3210 Google Scholar Crossref Search ADS WorldCat Harley P Guenther AB Zimmerman P Effects of light, temperature and canopy position on net photosynthesis and isoprene emission from sweetgum (Liquidambar styraciflua) leaves Tree Physiology 1996 16 25 32 Google Scholar Crossref Search ADS PubMed WorldCat Harrison SP Morfopoulos C Dani KGS et al. Volatile isoprenoid emissions from plastid to planet The New Phytologist 2013 197 49 57 Google Scholar Crossref Search ADS PubMed WorldCat Hauglustaine DA Lathière J Szopa S Folberth GA Future tropospheric ozone simulated with a climate–chemistry–biosphere model Geophysical Research Letters 2005 32 L24807 Google Scholar Crossref Search ADS WorldCat Heald CL Henze DK Horowitz LW et al. Predicted change in global secondary organic aerosol concentrations in response to future climate, emissions, and land use change Journal of Geophysical Research 2008 113 D05211 Google Scholar Crossref Search ADS WorldCat Heald CL Wilkinson MJ Monson RK Alo CA Wang G Guenther A Response of isoprene emission to ambient CO2 changes and implications for global budgets Global Change Biology 2009 15 1127 1140 Google Scholar Crossref Search ADS WorldCat Hecht S Eisenreich W Adam P et al. Studies on the nonmevalonate pathway to terpenes: the role of the GcpE (IspG) protein Proceedings of the National Academy of Sciences of the United States of America 2001 98 14837 14842 Google Scholar Crossref Search ADS PubMed WorldCat Hewitt CN Ashworth K Boynard A et al. Ground-level ozone influenced by circadian control of isoprene emissions Nature Geoscience 2011 4 671 674 Google Scholar Crossref Search ADS WorldCat Hikosaka K Murakami A Hirose T Balancing carboxylation and regeneration of ribulose-1,5- bisphosphate in leaf photosynthesis: temperature acclimation of an evergreen tree, Quercus myrsinaefolia Plant, Cell and Environment 1999 22 841 849 Google Scholar Crossref Search ADS WorldCat Karl T Fall R Rosenstiel TN et al. On-line analysis of the 13CO2 labeling of leaf isoprene suggests multiple subcellular origins of isoprene precursors Planta 2002 215 894 905 Google Scholar Crossref Search ADS PubMed WorldCat Kattge J Knorr W Temperature acclimation in a biochemical model of photosynthesis: a reanalysis of data from 36 species Plant, Cell & Environment 2007 30 1176 1190 Google Scholar Crossref Search ADS PubMed WorldCat Keenan TF Niinemets Ü Circadian control of global isoprene emissions Nature Geoscience 2012 5 435 Google Scholar Crossref Search ADS WorldCat Keenan TF Niinemets Ü Sabate S Gracia C Peñuelas J Process based inventory of isoprenoid emissions from European forests: model comparisons, current knowledge and uncertainties Atmospheric Chemistry and Physics 2009 9 4053 4076 Google Scholar Crossref Search ADS WorldCat Keenan TF Grote R Sabaté S Overlooking the canopy: the importance of canopy structure in scaling isoprenoid emissions from the leaf to the landscape Ecological Modelling 2011 222 737 747 Google Scholar Crossref Search ADS WorldCat Kreuzwieser J Graus M Wisthaler A Hansel A Rennenberg H Xylem-transported glucose as an additional carbon source for leaf isoprene formation in Quercus robur New Phytologist 2002 156 171 178 Google Scholar Crossref Search ADS WorldCat Laothawornkitkul J Taylor JE Paul ND Hewitt CN Biogenic volatile organic compounds in the Earth system The New Phytologist 2009 183 27 51 Google Scholar Crossref Search ADS PubMed WorldCat Lathière J Hewitt CN Beerling DJ Sensitivity of isoprene emissions from the terrestrial biosphere to 20th century changes in atmospheric CO2 concentration, climate, and land use Global Biogeochemical Cycles 2010 24 GB1004 Google Scholar Crossref Search ADS WorldCat Lawlor DW Tezara W Causes of decreased photosynthetic rate and metabolic capacity in water-deficient leaf cells: a critical evaluation of mechanisms and integration of processes Annals of Botany 2009 103 561 579 Google Scholar Crossref Search ADS PubMed WorldCat Lehning A Zimmer W Steinbrecher R Brüggemann N Schnitzler J-P Isoprene synthase activity and its relation to isoprene emission in Quercus robur L. leaves Plant, Cell and Environment 1999 22 495 504 Google Scholar Crossref Search ADS WorldCat Lerdau MT Keller M Controls on isoprene emission from trees in a subtropical dry forest Plant, Cell and Environment 1997 20 569 578 Google Scholar Crossref Search ADS WorldCat Lerdau MT Throop HL Isoprene emission and photosynthesis in a tropical forest canopy: implications for model development Ecological Applications 1999 9 1109 1117 Google Scholar Crossref Search ADS WorldCat Li Z Sharkey TD Metabolic profiling of the methylerythritol phosphate pathway reveals the source of post-illumination isoprene burst from leaves Plant, Cell & Environment 2012 36 429 437 Google Scholar Crossref Search ADS PubMed WorldCat Lichtenthaler HK The 1-deoxy-d-xylulose-5-phosphate pathway of isoprenoid biosynthesis in plants Annual Review of Plant Physiology and Plant Molecular Biology 1999 50 47 65 Google Scholar Crossref Search ADS PubMed WorldCat Lichtenthaler HK Schwender J Disch A Rohmer M Biosynthesis of isoprenoids in higher plant chloroplasts proceeds via a mevalonate-independent pathway FEBS Letters 1997 400 271 274 Google Scholar Crossref Search ADS PubMed WorldCat Logan BA Monson RK Potosnak MJ Biochemistry and physiology of foliar isoprene production Trends in Plant Science 2000 5 477 481 Google Scholar Crossref Search ADS PubMed WorldCat Loreto F Sharkey TD A gas-exchange study of photosynthesis and isoprene emission in Quercus rubra L Planta 1990 182 523 531 Google Scholar Crossref Search ADS PubMed WorldCat Loreto F Sharkey TD On the relationship between isoprene emission and photosynthetic metabolites under different environmental conditions Planta 1993 189 420 424 Google Scholar Crossref Search ADS PubMed WorldCat Loreto F Pinelli P Brancaleoni E 13C labeling reveals chloroplastic and extrachloroplastic pools of dimethylallyl pyrophosphate and their contribution to isoprene formation Plant Physiology 2004 135 1903 1907 Google Scholar Crossref Search ADS PubMed WorldCat Martin M Stirling CM Humphries SW Long SP A process-based model to predict the effects of climatic change on leaf isoprene emission rates Ecological Modelling 2000 131 161 174 Google Scholar Crossref Search ADS WorldCat McKinney KA Lee BH Vasta A Pho TV Munger JW Emissions of isoprenoids and oxygenated biogenic volatile organic compounds from a New England mixed forest Atmospheric Chemistry and Physics 2011 11 4807 4831 Google Scholar Crossref Search ADS WorldCat Medlyn BE Dreyer E Ellsworth D et al. Temperature response of parameters of a biochemically based model of photosynthesis. II. A review of experimental data Plant, Cell and Environment 2002 25 1167 1179 Google Scholar Crossref Search ADS WorldCat Monson RK Fall R Isoprene emission from aspen leaves?: influence of environment and relation to photosynthesis and photorespiration Plant Physiology 1989 90 267 274 Google Scholar Crossref Search ADS PubMed WorldCat Monson RK Jaeger CH Adams WW Driggers EM Silver GM Fall R Relationships among isoprene emission rate, photosynthesis, and isoprene synthase activity as influenced by temperature Plant Physiology 1992 98 1175 1180 Google Scholar Crossref Search ADS PubMed WorldCat Monson RK Grote R Niinemets Ü Schnitzler J-P Modeling the isoprene emission rate from leaves The New Phytologist 2012 195 541 559 Google Scholar Crossref Search ADS PubMed WorldCat Niinemets Ü Costs of production and physiology of emission of volatile leaf isoprenoids Advances in Plant Physiology 2004 7 233 268 Google Scholar OpenURL Placeholder Text WorldCat Niinemets Ü Mild versus severe stress and BVOCs: thresholds, priming and consequences Trends in Plant Science 2010 15 145 153 Google Scholar Crossref Search ADS PubMed WorldCat Niinemets Ü Tenhunen JD Harley PC Steinbrecher R A model of isoprene emission based on energetic requirements for isoprene synthesis and leaf photosynthetic properties for Liquidambar and Quercus Plant, Cell and Environment 1999 22 1319 1335 Google Scholar Crossref Search ADS WorldCat Niinemets Ü Copolovici L Hüve K High within-canopy variation in isoprene emission potentials in temperate trees: Implications for predicting canopy-scale isoprene fluxes Journal of Geophysical Research 2010 115 G04029 Google Scholar Crossref Search ADS WorldCat Nozière B González NJD Borg-Karlson A-K et al. Atmospheric chemistry in stereo: a new look at secondary organic aerosols from isoprene Geophysical Research Letters 2011 38 L11807 Google Scholar Crossref Search ADS WorldCat Onoda Y Hikosaka K Hirose T The balance between RuBP carboxylation and RuBP regeneration: a mechanism underlying the interspecific variation in acclimation of photosynthesis to seasonal change in temperature Functional Plant Biology. 2005 32 903 910 Google Scholar Crossref Search ADS WorldCat Pacifico F Harrison S. P. Jones CD Sitch S Isoprene emissions and climate Atmospheric Environment 2009 43 6121 6135 Google Scholar Crossref Search ADS WorldCat Pacifico F Harrison SP Jones CD et al. Evaluation of a photosynthesis-based biogenic isoprene emission scheme in JULES and simulation of isoprene emissions under present-day climate conditions Atmospheric Chemistry and Physics 2011 11 4371 4389 Google Scholar Crossref Search ADS WorldCat Pacifico F Folberth GA Jones CD Harrison SP Collins WJ Sensitivity of biogenic isoprene emissions to past, present, and future environmental conditions and implications for atmospheric chemistry Journal of Geophysical Research 2012 117 D22302 Google Scholar OpenURL Placeholder Text WorldCat Palmer PI Jacob DJ M. FA Martin RV Mapping isoprene emissions over North America using formaldehyde column observations from space Journal of Geophysical Research 2003 108 4180 Google Scholar Crossref Search ADS WorldCat Palmer PI Abbot DS Fu T-M et al. Quantifying the seasonal and interannual variability of North American isoprene emissions using satellite observations of the formaldehyde column Journal of Geophysical Research 2006 111 D12315 Google Scholar OpenURL Placeholder Text WorldCat Pétron G Harley P Greenberg J Guenther A Seasonal temperature variations influence isoprene emission Geophysical Research Letters 2001 28 1707 1710 Google Scholar Crossref Search ADS WorldCat Pike RC Young PJ How plants can influence tropospheric chemistry: the role of isoprene emissions from the biosphere Weather 2009 64 332 336 Google Scholar Crossref Search ADS WorldCat Poisson N Kanakidou M Crutzen P. J Impact of non-methane hydrocarbons on tropospheric chemistry and the oxidizing power of the global troposphere?: 3-dimensional modelling results Journal of Atmospheric Chemistry 2000 36 157 230 Google Scholar Crossref Search ADS WorldCat Possell M Hewitt CN Isoprene emissions from plants are mediated by atmospheric CO2 concentrations Global Change Biology 2011 17 1595 1610 Google Scholar Crossref Search ADS WorldCat Rasulov B Hüve K Väalbe M Laisk A Niinemets Ü Evidence that light, carbon dioxide, and oxygen dependencies of leaf isoprene emission are driven by energy status in hybrid aspen Plant Physiology 2009 151 448 460 Google Scholar Crossref Search ADS PubMed WorldCat Rasulov B Hüve K Bichele I Laisk A Niinemets Ü Temperature response of isoprene emission in vivo reflects a combined effect of substrate limitations and isoprene synthase activity: a kinetic analysis Plant Physiology 2010 154 1558 1570 Google Scholar Crossref Search ADS PubMed WorldCat Rasulov B Hüve K Laisk A Niinemets Ü Induction of a longer term component of isoprene release in darkened aspen leaves: origin and regulation under different environmental conditions Plant Physiology 2011 156 816 831 Google Scholar Crossref Search ADS PubMed WorldCat Rosenstiel TN Potosnak MJ Griffin KL Fall R Monson RK Increased CO2 uncouples growth from isoprene emission in an agriforest ecosystem Nature 2003 421 256 259 Google Scholar Crossref Search ADS PubMed WorldCat Rosenstiel TN Ebbets AL Khatri WC Fall R Monson RK Induction of poplar leaf nitrate reductase: a test of extrachloroplastic control of isoprene emission rate Plant Biology 2004 6 12 21 Google Scholar Crossref Search ADS PubMed WorldCat Sanderson MG Jones CD Collins WJ Johnson CE Derwent RG Effect of climate change on isoprene emissions and surface ozone levels Geophysical Research Letters 2003 30 1936 Google Scholar OpenURL Placeholder Text WorldCat Schnitzler J-P Lehning A Steinbrecher R Seasonal pattern of isoprene synthase activity in Quercus robur leaves and its significance for modeling isoprene emission rates Botanica Acta 1997 110 240 243 Google Scholar Crossref Search ADS WorldCat Seemann M Tse Sum Bui B Wolff M Miginiac-Maslow M Rohmer M Isoprenoid biosynthesis in plant chloroplasts via the MEP pathway: direct thylakoid/ferredoxin-dependent photoreduction of GcpE/IspG FEBS Letters 2006 580 1547 1552 Google Scholar Crossref Search ADS PubMed WorldCat Sharkey TD Loreto F Water stress, temperature, and light effects on the capacity for isoprene emission and photosynthesis of kudzu leaves Oecologia 1993 95 328 333 Google Scholar Crossref Search ADS PubMed WorldCat Sharkey TD Yeh S Isoprene emission from plants Annual Review of Plant Physiology and Plant Molecular Biology 2001 52 407 436 Google Scholar Crossref Search ADS PubMed WorldCat Sharkey TD Singsaas EL Vanderveer PJ Geron C Field measurements of isoprene emission from trees in response to temperature and light Tree Physiology 1996 16 649 654 Google Scholar Crossref Search ADS PubMed WorldCat Sharkey TD Chen X Yeh S Isoprene increases thermotolerance of fosmidomycin-fed leaves Plant Physiology 2001 125 2001 2006 Google Scholar Crossref Search ADS PubMed WorldCat Sharkey TD Wiberley AE Donohue AR Isoprene emission from plants: why and how Annals of Botany 2008 101 5 18 Google Scholar Crossref Search ADS PubMed WorldCat Silver GM Fall R Characterization of aspen isoprene synthase, an enzyme responsible for leaf isoprene emission to the atmosphere The Journal of Biological Chemistry 1995 270 13010 13016 Google Scholar Crossref Search ADS PubMed WorldCat Singarayer JS Valdes PJ Friedlingstein P Nelson S Beerling DJ Late Holocene methane rise caused by orbitally controlled increase in tropical sources Nature 2011 470 82 85 Google Scholar Crossref Search ADS PubMed WorldCat Singsaas EL Ort DR DeLucia EH Variation in measured values of photosynthetic quantum yield in ecophysiological studies Oecologia 2001 128 15 23 Google Scholar Crossref Search ADS WorldCat Stavrakou T Müller J-F De Smedt I et al. Evaluating the performance of pyrogenic and biogenic emission inventories against one decade of space-based formaldehyde columns Atmospheric Chemistry and Physics 2009 9 1037 1060 Google Scholar Crossref Search ADS WorldCat Sun Z Niinemets Ü Hüve K et al. Enhanced isoprene emission capacity and altered light responsiveness in aspen grown under elevated atmospheric CO2 concentration Global Change Biology 2012 18 3423 3440 Google Scholar Crossref Search ADS WorldCat Trowbridge AM Asensio D Eller ASD et al. Contribution of various carbon sources toward isoprene biosynthesis in poplar leaves mediated by altered atmospheric CO2 concentrations PloS ONE 2012 7 e32387 Google Scholar OpenURL Placeholder Text WorldCat Urbanski S Barford C Wofsy S et al. Factors controlling CO2 exchange on timescales from hourly to decadal at Harvard Forest Journal of Geophysical Research 2007 112 G02020 Google Scholar OpenURL Placeholder Text WorldCat Valdes PJ Beerling DJ Johnson CE The ice age methane budget Geophysical Research Letters 2005 32 L02704 Google Scholar OpenURL Placeholder Text WorldCat Velikova V Váarkonyi Z Szabóo M et al. Increased thermostability of thylakoid membranes in isoprene-emitting leaves probed with three biophysical techniques Plant Physiology 2011 157 905 916 Google Scholar Crossref Search ADS PubMed WorldCat Velikova V Sharkey TD Loreto F Stabilization of thylakoid membranes in isoprene-emitting plants reduces formation of reactive oxygen species Plant Signaling & Behavior 2012 7 139 141 Google Scholar Crossref Search ADS PubMed WorldCat Vickers CE Gershenzon J Lerdau MT Loreto F A unified mechanism of action for volatile isoprenoids in plant abiotic stress Nature Chemical Biology 2009 5 283 291 Google Scholar Crossref Search ADS PubMed WorldCat Vickers CE Possell M Hewitt CN Mullineaux PM Genetic structure and regulation of isoprene synthase in Poplar (Populus spp.) Plant Molecular Biology 2010 73 547 558 Google Scholar Crossref Search ADS PubMed WorldCat Wilkinson MJ Monson RK Trahan N et al. Leaf isoprene emission rate as a function of atmospheric CO2 concentration Global Change Biology 2009 15 1189 1200 Google Scholar Crossref Search ADS WorldCat Wolfertz M Sharkey TD Boland W Kühnemann F Yeh S Weise SE Biochemical regulation of isoprene emission Plant, Cell and Environment 2003 26 1357 1364 Google Scholar Crossref Search ADS WorldCat Young PJ Arneth A Schurgers G Zeng G Pyle JA The CO2 inhibition of terrestrial isoprene emission significantly affects future ozone projections Atmospheric Chemistry and Physics 2009 9 2793 2803 Google Scholar Crossref Search ADS WorldCat Zimmer W Brüggemann N Emeis S Giersch C Lehning A Steinbrecher R Schnitzler J-P Process-based modelling of isoprene emission by oak leaves Plant, Cell and Environment 2000 23 585 595 Google Scholar Crossref Search ADS WorldCat Zimmer W Steinbrecher R Körner C Schnitzler J-P The process-based SIM – BIM model?: towards more realistic prediction of isoprene emissions from adult Quercus petraea forest trees Atmospheric Environment 2003 37 1665 1671 Google Scholar Crossref Search ADS WorldCat © The Author 2013. Published by Oxford University Press on behalf of the Annals of Botany Company. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited © The Author 2013. Published by Oxford University Press on behalf of the Annals of Botany Company.
TI - A unifying conceptual model for the environmental responses of isoprene emissions from plants
JF - Annals of Botany
DO - 10.1093/aob/mct206
DA - 2013-11-01
UR - https://www.deepdyve.com/lp/oxford-university-press/a-unifying-conceptual-model-for-the-environmental-responses-of-3yxCAMghla
SP - 1223
EP - 1238
VL - 112
IS - 7
DP - DeepDyve
ER -