In Vivo NADH/NAD1 Biosensing Reveals the Dynamics of Cytosolic Redox Metabolism in Plants[]

Janina Steinbeck,a Philippe Fuchs,a,b Yuri L. Negroni,b,c Marlene Elsässer,a,b,d Sophie Lichtenauer,a Yvonne Stockdreher,b Elias Feitosa-Araujo,a,e Johanna B. Kroll,a Jan-Ole Niemeier,a Christoph Humberg,a Edward N. Smith,f Marie Mai,g Adriano Nunes-Nesi,e Andreas J. Meyer,b Michela Zottini,c Bruce Morgan,g Stephan Wagner,a,b,1 and Markus Schwarzländera,2 a Institute of Plant Biology and Biotechnology (IBBP), Westfälische Wilhelms-Universität Münster, D-48143 Münster, Germany b Institute of Crop Science and Resource Conservation (INRES), Rheinische Friedrich-Wilhelms-Universität Bonn, D-53113 Bonn, Germany cDepartment of Biology, University of Padova, 35131 Padova, Italy d Institute for Cellular and Molecular Botany (IZMB), Rheinische Friedrich-Wilhelms-Universität Bonn, D-53115 Bonn, Germany eDepartamento de Biologia Vegetal, Universidade Federal de Viçosa, 36570-900 Viçosa, Minas Gerais, Brazil f Department of Plant Sciences, University of Oxford, Oxford, OX1 3RB, United Kingdom g Institute of Biochemistry, Zentrum für Humanund Molekularbiologie (ZHMB), Saarland University, D-66123 Saarbrücken, Germany


INTRODUCTION
Nicotinamide adenine dinucleotide (NAD) acts as a central cofactor of redox metabolism across organisms, including plants (Noctor, 2006;Gakière et al., 2018). The oxidized and reduced forms of NAD, NAD 1 , and NADH, are involved in the electron transfer reactions of a large number of dehydrogenase enzymes, including those of glycolysis, the tricarboxylic acid (TCA) cycle, and the mitochondrial electron transport chain (mETC; Fernie et al., 2004). NAD also regulates enzyme activity, and NAD 1 is required for the de-acetylation of histones and metabolic proteins (Møller et al., 2020). Modification of cellular NAD homeostasis can have pronounced, but remarkably specific, effects on plant hormone regulation, reproductive efficiency, and guard cell density, providing a link between redox metabolism and plant developmental programs (Hashida et al., 2007(Hashida et al., , 2010de Souza Chaves et al., 2019;Feitosa-Araujo et al., 2020;Wei et al., 2020). Phosphorylation of NAD 1 by NAD 1 kinase produces NADP 1 and in analogy, phosphorylation of NADH by NADH kinase produces NADPH (Li et al., 2018). NADP has vital roles in anabolic metabolism, such as photosynthetic carbon fixation, and it is an important electron donor for thiol-based redox regulation and antioxidant defense (Møller et al., 2020). Consequently, NAD constitutes a backbone of cellular redox metabolism.
Since the NAD redox state is directly coupled to the flux through many metabolic reactions and has fundamental impacts on fluxes throughout the entire plant metabolic network, extremes of NAD oxidation and reduction can disrupt metabolism. When oxidative phosphorylation capacity is limited, for example, under hypoxia, the resulting Pasteur effect includes a boost in glycolytic flux to maintain ATP production. However, this leads to a high rate of NAD 1 reduction and requires the activation of fermentation of pyruvate to ethanol or lactate (Tadege et al., 1999;António et al., 2016). Fermentation allows NADH to be re-oxidized to avoid largescale disruption of metabolic function, albeit at the major cost of inefficient carbon and energy use, highlighting the vital significance of NAD redox maintenance. A similar function is fulfilled by the alternative respiratory pathways in plant mitochondria, including the alternative NADH dehydrogenases and alternative oxidase (van Dongen et al., 2011). Although oxidation of NADH via the alternative respiratory pathway still requires oxygen as a terminal electron acceptor, it provides a mechanism for maintaining NAD redox balance independent of the capacity of the classical respiratory pathway, ATP-synthase activity, and adenylate charge (Wagner et al., 2018).
The redox state of the cellular NAD pools can fluctuate; such changes have been associated with a variety of physiological conditions and stress situations (Hashida et al., 2007;Hunt et al., 2007;Schippers et al., 2008;Pétriacq et al., 2016). NAD redox management by the plant cell is likely to include moving major quantities of reducing equivalents between different cell compartments. While NAD 1 can be transported across the membranes of different cellular organelles, including those of the chloroplast, mitochondrion, and peroxisome (Palmieri et al., 2009;Bernhardt et al., 2012;de Souza Chaves et al., 2019;Luo et al., 2019;Feitosa-Araujo et al., 2020), there is no evidence for any efficient direct intracellular transmembrane transport of NADH. This highlights the need for metabolite shuttles to move redox equivalents between intracellular compartments. Among different mechanisms, the malate valves that operate between the cytosol and the plastid stroma, mitochondrial matrix, and peroxisomal lumen have attracted major attention due to their potential to reach major flux rates in plant cells (Selinski and Scheibe, 2019). The enzymes involved and their regulation have been studied intensely, and a set of dicarboxylate carriers present in the inner mitochondrial membrane, the inner plastid envelope, and the peroxisome membrane have been implicated in the required transmembrane transport (Weber and Flügge, 2002;Renné et al., 2003;Lee and Millar, 2016). The directionality of shuttling and the flux rates involved in vivo have remained less clear, however. Both parameters strongly depend on the exact metabolic situation, making it likely that they are highly variable depending on developmental and environmental conditions. Even assuming a metabolic steady state, any reliable estimation of flux rates through the redox shuttle systems in the living cell remains extremely difficult, even at a qualitative level. Major advances have recently been made by metabolic modeling (Shameer et al., 2019), while a reliable means for direct experimental assessment remains to be established.
Despite the central role of plant NAD in metabolism and its emerging roles in regulation and signaling, the current understanding of NAD redox dynamics in vivo, NAD compartmentation, and redox exchange between the cellular compartments remains surprisingly sketchy. It has been technically challenging to monitor NAD redox dynamics in vivo. Typically, NADH/NAD 1 ratios are determined by biochemical analyses of whole-cell extracts, in which metabolite pools from different cell compartments have been mixed. Such analyses therefore offer little insight into subcellular compartment-specific NADH/NAD 1 dynamics inside living cells. Hans Krebs established a method for indirectly deducing cytosolic free NADH/NAD 1 ratios from metabolite ratios more than 50 years ago. Such analysis is typically based on measuring pyruvate and lactate concentrations and the assumption of their equilibrium in situ (Williamson et al., 1967;Sun et al., 2012). This method has been widely exploited, but dynamic in vivo measurements have remained nonaccessible for many systems since they are technically highly demanding (Christensen et al., 2014). Although dynamic in vivo measurements of NAD(P)H concentrations are possible based on the autofluorescence of the reduced species, as proposed nearly 60 years ago (Chance et al., 1962), for plants, this remains limited to isolated organelles due to interfering autofluorescence from other cellular components and the high intensity of UV illumination required (Kasimova et al., 2006). Even in systems that allow measurement of NAD(P)H autofluorescence, NAD(P)H in the mitochondrial matrix dominates the cellular signals, thereby impeding attempts to assess NAD(P)H dynamics in other cell compartments. Furthermore, NADH shows the same spectral properties as NADPH, resulting in a mixed readout of both independent pools. An elegant approach that was developed to separate both species using their different fluorescence lifetimes (Blacker et al., 2014) offers several advantages, but it does not overcome the hurdle of distinguishing NAD redox status between cell compartments and has not been optimized for plants.
Several fluorescent protein-based NAD biosensors were recently developed, providing novel insight into NAD dynamics in mammalian, fungal, and bacterial systems (Hung et al., 2011;Zhao et al., 2011Zhao et al., , 2015Bhat et al., 2016;Cambronne et al., 2016;Ishikawa et al., 2017;Tejwani et al., 2017;Hartmann et al., 2018). One of these sensors, Peredox-mCherry, was designed using a T-Rex NAD 1 -and NADH binding protein linked to the blue-green fluorescent protein tSapphire (tS; Hung et al., 2011). The binding of NADH to this sensor causes a steric rearrangement that leads to an increase in fluorescence intensity. NAD 1 competes for the binding site, but its binding has a neutral impact on the fluorescence emission intensity. As a result, tS fluorescence intensity depends on the ambient NADH/NAD 1 ratio. For internal normalization C-terminal fusion with mCherry (mC), which is inert to the NADH/NAD 1 ratio, has been used to generate a ratiometric sensor construct. The pH sensitivity of the sensor was minimized by targeted mutagenesis (Hung et al., 2011).
Alternative NAD sensors include Frex, FrexH, RexYFP, and SoNar (Zhao et al., 2011(Zhao et al., , 2015Bilan et al., 2014; for a review, see Bilan and Belousov, 2017). Frex and FrexH provide a readout for absolute NADH levels but show pronounced pH sensitivity; SoNar allows NADH/NAD 1 measurements with a strongly increased spectroscopic response range but is also highly pH sensitive. A sensor for NAD 1 concentrations based on cpVenus fused with a bipartite NAD 1 binding site of a bacterial DNA ligase was recently developed and used in the mitochondrial matrix (Cambronne et al., 2016), and a semisynthetic fluorescent NAD 1 biosensor, NAD-Snifit, was subsequently introduced (Sallin et al., 2018). Very recently, a sensor named FiNAD was developed that responds to NAD 1 /AXP ratios, where AXP is the pool of ATP and ADP (Zou et al., 2020). Despite the important advances that NAD biosensing has delivered in mammalian and fungal systems, and the potential to gain insight into the flexibility of subcellular plant metabolism, these NAD biosensors have not been established in plants. One recent study made use of NAD sensing as part of a multiparameter approach for analyzing hypoxia physiology, but without specific development of the methodology (Wagner et al., 2019). During the revision of this work, a study was published exploiting SoNar in combination with a pH correction regime for NAD redox monitoring in the cytosol and the chloroplast stroma of Arabidopsis (Arabidopsis thaliana; Lim et al., 2020).
Here, based on the properties of the currently available biosensors, we focused on Peredox-mCherry to establish a method for sensing the NAD redox status in living plants. Peredox-mCherry was the only sensor that (1) provides an integrated measurement of NAD redox status, as opposed to NAD 1 or NADH concentration; (2) can be measured in a self-normalizing ratiometric manner; and (3) was reported to be largely pH inert. First, we expressed Peredox-mCherry in the cytosol to generate stable Arabidopsis biosensor lines. Second, we validated the biophysical and biochemical characteristics of Peredox-mCherry in vitro under conditions that are expected to be relevant in planta. Third, we imaged the NAD redox landscape in living Arabidopsis seedlings to reveal tissuespecific differences and gradients in cytosolic NAD redox status. Analysis of etiolated seedlings demonstrated that the establishment of these gradients is flexibly modified depending on developmental status and the environment. Fourth, we investigated the impact of short-term transitions in photosynthetic and respiratory activity, as well as external carbon availability and elicitor exposure, on cytosolic NAD redox status and resolved the kinetic responses of NAD redox metabolism. Finally, we engineered a sensor variant for use in more reducing NAD redox environments to provide a usable toolkit for NAD redox monitoring in plants.

Arabidopsis Lines Expressing an NAD Redox Sensor in the Cytosol
To perform NAD redox measurements in living plants, we generated Arabidopsis lines that expressed Peredox-mCherry constitutively under the control of the UBIQUITIN10 promoter in the cytosol. We selected three independent lines that showed intense fluorescence signal, which localized to the peripheral cytosol, trans-vacuolar strands, and nucleoplasm (Figures 1A to 1E;Supplemental Figures 1A and 1B). When we attempted to generate analogous sensor lines for the mitochondrial matrix, selection for transformed individuals yielded no results. This is in line with the previous observation that the import of large biosensor constructs into mitochondria compromises plant development (De Col et al., 2017). Even though severe growth phenotypes can also occur for specific biosensors expressed in the cytosol, we observed no phenotypic differences at the whole-plant level between the homozygous Peredox-mCherry lines and the wild type without sensor expression. Quantitative phenotyping of rosette size, inflorescence height, and root length (Figures 1F to 1I;Supplemental Figures 1C and 1D) confirmed that sensor expression did not induce apparent phenotypes. Homozygous line #1 was an exception, as it displayed slightly increased inflorescence height and root growth compared to the wild-type control, perhaps due to the specific location of the T-DNA insertion rather than sensor expression. We did not observe any obvious reproductive phenotypes, but the potential impact on a specific parameter of interest should be assessed when the lines are used in future studies.
pH Stability and Specificity for NAD Make Peredox-mCherry Useful for in Planta Measurements To allow for reliable interpretation of in planta sensor measurements, we sought to obtain additional insight into the specificity of the sensor response, complementing previous in vitro analyses of the Peredox-mCherry protein (Hung et al., 2011;Tejwani et al., 2017;Hartmann et al., 2018). In particular, we were interested in sensor specificity to NAD versus NADP, as well as changes in ambient pH and sensor response kinetics, to appropriately consider their influence on in planta measurements. The model of Peredox-mCherry function suggests that the homodimeric Rex domain undergoes a conformational change in response to NADH binding, thereby increasing the fluorescence intensity of the circularly permuted GFP-derived protein tS, which is integrated between the two Rex monomers (Figure 2A; Hung et al., 2011). The fluorescence spectra of the C-terminal mC remain unaffected, allowing the tS signal to be internally normalized. Since NAD 1 binds to the sensor, but does not affect its fluorescence, the tS/mC fluorescence ratio reports on the NADH/NAD 1 ratio rather than NADH concentration (Hung et al., 2011).
To obtain purified recombinant Peredox-mCherry sensor protein, we optimized the expression conditions, Escherichia coli strain, and protein storage conditions (Supplemental Figure 2). Measuring the emission spectrum of the purified sensor protein confirmed a ratiometric shift in response to NADH, with a dynamic range (d) of 2.2 ( Figure 2B; Supplemental Figure 3). tS fluorescence increased in response to NADH, but not to NAD 1 (excitation at 400 nm, emission peak at 513 nm), while mC fluorescence remained unchanged (excitation at 580 nm, emission peak at 605 nm). Medium NAD 1 concentrations of 500 mM were chosen to approximate the lower end of the range of cytosolic NAD 1 concentrations previously measured from extracts of mesophyll protoplasts (Igamberdiev and Gardeström, 2003;Gakière et al., 2018). Adjusting NADH concentrations from 0 to 100 mM in the presence of 500 mM NAD 1 led to a sigmoidal increase in the tS/mC fluorescence ratio, with a K d -value for NADH of 1.2 6 0.15 mM ( Figure 2C; Supplemental Figures 3A and 3B). Cytosolic NADH concentrations of pea (Pisum sativum) leaf protoplasts in the dark were previously measured at 18 6 14 mM NADH (Igamberdiev and Gardeström, 2003). A response to NADPH in the presence of 150 mM NADP 1 (previously estimated as the concentration of the cytosolic NADP 1 pool; Igamberdiev and Gardeström, 2003) was only observable at concentrations two orders of magnitude higher (K d for NADPH, 531 6 12 mM; Figure 2C; Supplemental Figures 3C and 3D). As a result, the in planta response of Peredox-mCherry is likely to be specific to NAD redox status and will most likely not be affected by the NADP pool under physiological conditions. tS/mC was largely stable within the physiologically relevant pH range from 6 to 9 ( Figure 2D), as assessed for three different conditions of NAD 1 and NADH binding ( Figure 2C): occupied by NAD 1 (0.1 mM NADH, 500 mM NAD 1 ), occupied by NADH (5 mM NADH, 500 mM NAD 1 ), and at similar occupancy by NAD 1 and NADH (1.2 mM NADH, 500 mM NAD 1 ; Figure 2D). The residual pH dependency observed was minor compared to the NAD redox response. The absence of direct effects of pH on the fluorophores is particularly important, since pH can fluctuate substantially in the plant cytosol (e.g., Felle, 2001;Behera et al., 2018), and pH stability represents a decisive advantage over several other NAD sensors (Zhao et al., 2011(Zhao et al., , 2015. Finally, we investigated the binding and dissociation kinetics of NADH to and from the sensor protein in the presence of NAD 1 as a competitor for binding ( Figure 2E). The tS/mC ratio values were also log 10 transformed, allowing us to accurately represent the ratio variance and to restore symmetry of the ratio data; all subsequent data are presented in this manner. Equilibration of the sensor protein with NAD 1 prior to NADH addition ensured that no NADH was bound when the assay was started. At NADH addition (1 mM), tS/mC increased to half of the full response amplitude in 73 6 1 s, whereas the dissociation of NADH, as triggered by the addition of NAD 1 (500 mM), was slower and required 211 6 7 s to reach half of the response amplitude (a more detailed analysis of the binding and dissociation kinetics for different NADH concentrations is (I) Vertical primary root length on half-strength MS medium supplemented with 1% (w/v) Suc and 1% (w/v) agar (n 5 25, means 6 SE). Two-way ANOVA followed by Dunnett's multiple comparisons test for data in (G) to (I). Note that line #1 showed significantly increased inflorescence height and root length (P < 0.001) compared to Col-0; rosette size of line #1 and all comparisons of lines #2 and #3 were not significantly different from the Col-0 control (P > 0.05; Supplemental Data Sets 1A to 1C).  (A) Hypothetical structural model of Peredox-mCherry. tS (Protein Data Bank [PDB]: 3evp, green) is interposed between two Rex subunits (SU; PDB: 2vt2 and 1xcb, gray). The binding of NADH (red, right) increases tS fluorescence, while NAD 1 binding does not (red, left). tS signal is normalized by the mC (PDB: 2 h5q, orange) signal, which is unresponsive to NAD status. mC is fused to the C terminus of the second Rex domain. (B) Emission spectra of tS (excitation at 400 6 5 nm, left y axis) and mC (excitation at 580 6 5 nm, right y axis) in arbitrary units (AU) 3 10 3 after the addition of 500 mM NAD 1 (blue) or 500 mM NAD 1 plus 100 mM NADH (gray). The spectroscopic dynamic range (d) was calculated from the tS/mC ratios at the respective emission maxima of the two fluorophores (tS, 513 6 5 nm; mC, 610 6 5 nm) for the NAD 1 -bound and NADH-bound sensor. Sensor protein concentration: 0.025 mg/mL in Tris-HCl, pH 7.5 (also for subsequent measurements unless otherwise stated). (C) Ratiometric response of Peredox-mCherry to different NADH and NADPH concentrations in the presence of 500 mM NAD 1 and 150 mM NADP 1 , respectively (n 5 3, mean 6 SD), pH 7.5. The ratio values were normalized to the highest ratio (set to 1) to allow comparison of NADH and NADPH binding (corresponding tS and mC spectra are presented in Supplemental Figures 3A to 3D). Data were fitted to a sigmoidal interpolation (Supplemental Data Set 1D). Dotted lines indicate K d -values for NAD(P)H. (D) Response of Peredox-mCherry to different pH values. Sensor protein was equilibrated with three different NADH concentrations (0.1, 1.2, and 5 mM) in the presence of 500 mM NAD 1 to achieve low, mean, and high NADH/NAD 1 binding at different pH values (6.0 to 9.0, n 5 3, mean 6 SD). The tS/mC fluorescence ratios were determined as in (C) at the indicated pH of the buffer (Bis-Tris for pH 6.0 to 7.0 and Tris-HCl for pH 7.5 to 9.0). presented in Supplemental Figures 3E and 3F). Consequently, the binding and dissociation kinetics may limit the ability of the sensor to monitor very rapid NAD redox dynamics. While the delayed sensor response represents a disadvantage that needs to be considered when interpreting the temporal responses of in vivo experiments, this behavior will also extend the time window in which very rapid metabolic transitions between a more oxidized and a more reduced NAD pool (and vice versa) can be resolved. Endogenous proteins can also show low rates in NAD binding and NAD dissociation, which has given rise to the concept of bound and free NAD pools. Here, Peredox-mCherry makes no exception and provides insight into NAD redox dynamics, as occurs for endogenous proteins with similar binding and dissociation rates.

Peredox-mCherry Reports on Changes in Net NAD Reduction and Oxidation by Mitochondrial Malate Metabolism ex Situ
Since NADH and NAD 1 represent a central redox couple of respiratory metabolism, we aimed to explore how mitochondrial metabolic activities can affect NAD redox status depending on substrate availability. We recently established an ex situ strategy to assess mitochondria-mediated changes in ATP (De Col et al., 2017), and we have adapted that strategy here to monitor NAD redox dynamics as mediated by malate metabolism. We supplemented suspensions of Arabidopsis seedling mitochondria, purified Peredox-mCherry protein, and NAD 1 (500 mM) with different concentrations of malate ( Figure 3A) to test the ability of the NAD redox sensor to report NAD redox changes from a complex organelle system. Malate is readily imported into the matrix of plant mitochondria, where it is metabolized by malate dehydrogenase (MDH) and malic enzyme at high capacity coupled to the reduction of NAD 1 (Zoglowek et al., 1988;Fuchs et al., 2020). Any efficient direct export of NADH from intact mitochondria into the medium is unlikely (Douce et al., 1973;Palmieri et al., 2009). We reasoned that the rupture of a fraction of the mitochondria would bring matrix enzymes in direct contact with NAD 1 and the sensor in the medium; alternatively, MDH activity is typically also present at the outer mitochondrial membrane of isolated mitochondria. Once generated, NADH can be re-oxidized by the NADH dehydrogenases of the mETC (Rasmusson et al., 2004).
The addition of malate at 0.5, 1.5, and 10 mM indeed resulted in a dose-dependent reduction of the NAD pool, as indicated by the increase in the log 10 (tS/mC) signal ratio by 0.10, 0.26, and 0.31, respectively ( Figure 3B), with a corresponding spectroscopic dynamic range of 2.04, that is, similar to the dynamic range determined by the external addition of NADH ( Figure 2C). MDH activity likely accounted for most of the reduction of NAD 1 to NADH, linked to malate oxidation, since the assay pH is far above the pH optimum of malic enzyme. The re-oxidation of the NAD pool, as indicated by a decreasing tS/mC, can be explained by the notion that the reaction catalyzed by MDH reached equilibrium (note that the reaction has a very low equilibrium constant favoring malate and NAD 1 over oxaloacetate and NADH; Coleman and Palmer 1972). Since further malate oxidation and NAD 1 reduction are readily inhibited through the accumulation of oxaloacetate and NADH, net oxidation of NADH can occur via a reaction mediated by the activity of mETC NADH-dehydrogenases, including complex I and the alternative NADH dehydrogenases (Rasmusson et al., 2004). Even though malate will not be fully depleted, the added malate concentration determines how long it will take until equilibrium will be reached and net NADH oxidation can commence. The NAD redox status is therefore set by the rates of NAD 1 reduction and NADH oxidation in the ex situ model system. Since net reduction and net oxidation dominate at different time windows of the assay, which can be modified by changing the amount of added substrate, the dose-dependent dynamics verify the ability of Peredox-mCherry to report NAD redox dynamics, like those that can be induced by changes in the fluxes through respiratory carbon metabolism. (E) NADH binding to and dissociation from Peredox-mCherry in vitro. Time resolved tS/mC (gray, left y axis) and corresponding log 10 -transformed ratio changes (black, right y axis) in response to NADH (1.2 mM, first arrow) in the presence of 500 mM NAD 1 and to further NAD 1 addition (500 mM, second arrow). Dotted lines indicate time to reach half of the response amplitude at NADH or NAD 1 addition (linear regression analysis of log 10 -transformed ratios (onephase decay fit; Supplemental Data Set 1E). Fluorescence emission was measured by multiwell plate reader-based fluorimetry in each well every 13 s. Excitation, 400 6 10 nm (tS) and 570 6 10 nm (mC); emission, 520 6 5 nm (tS) and 610 6 5 nm (mC). The assay was independently repeated three times with consistent results. (A) Mitochondria isolated from the 16-d-old Arabidopsis wild-type seedlings were mixed with purified recombinant Peredox-mCherry protein (0.5 mM) and NAD 1 (500 mM) in basic incubation medium, pH 7.5, and supplemented with malate. The sensor protein provides a readout for redox changes in the NAD pool in the medium. ETC, electron transport chain; metab., metabolism. (B) NAD redox dynamics in response to the addition of malate at three different concentrations compared to the control (arrow). tS/mC ratio indicates redox changes in the NAD pool in the medium. Fluorescence emission was measured in each well of a multiwell plate every 149 s and log 10 transformed. Averaged sensor ratios are plotted after correction of background signal. Excitation, 400 6 10 nm (tS) and 570 6 10 nm (mC); emission, 515 6 7.5 nm (tS) and 610 6 5 nm (mC). n 5 4, mean 1 SD.

A Seedling-Scale Atlas of the Cytosolic NAD Redox Landscape
We next sought to explore how heterogeneous the cytosolic NAD redox status is across tissues. This issue has been unclear because NAD measurements in tissue extracts can only provide an average value across cell compartments, cell types, and tissues, resulting in the loss of information about specific tissues and cell types. To address this issue, we performed confocal laser scanning microscopy analysis of cytosolic NAD redox state across entire 5-d-old Arabidopsis seedlings ( Figure 4A). The log 10 (tS/mC) ratios ranged from 0.38 6 0.06 (indicative of a highly reduced NAD state) in the root epidermis down to 0.005 6 0.049 (indicative of a highly oxidized NAD state) in the area of the leaf primordia in the vicinity of the shoot apical meristem (SAM) region ( Figures 4A and 4B). The resulting log 10 -response range of 0.37 (d 5 2.35) equals the theoretical range of the sensor determined in vitro ( Figure 2C), suggesting that the sensor affinities for NADH and NAD 1 are well matched with the cytosolic NAD redox environment in planta, but also that the NAD pool in the tissues with the lowest and highest values may adopt redox states beyond the sensitivity range of the sensor. To control for potential artifacts of steady state measurements through the contributions of tissue autofluorescence, we also imaged wild-type seedlings without sensor as a control (Supplemental Figure 4A). The average contribution was minor, suggesting that meaningful measurements can be obtained without further correction (autofluorescence contributions empirically determined for representative seedlings of line #3 were 2.9% for tS and 0.4% for mC signal; Supplemental Figure 4A). We then explored potential differences in sensor status between the cytosol and the nucleoplasm. The assessment of guard cells, pavement cells, and hypocotyl epidermal cells did not reveal any difference, which is in line with the finding that the sensor protein showed efficient diffusion across the nuclear pores (Supplemental Figure 4B).
The tS/mC ratios of cotyledons and hypocotyl were low (more oxidized NAD pool) relative to the higher tS/mC ratios in the root (more reduced NAD pool; Figures 4A and 4B). The observed gradient of increasing NAD reduction from shoot to root negatively correlates with the gradient in the concentrations of cytosolic MgATP 2-, which are high in shoots and lower in roots, as detected with the biosensor ATeam 1.03 nD/nA (De Col et al. 2017). The epidermis displayed a more reduced NAD pool compared to the vasculature in both cotyledons (epidermis, 0.13 6 0.03; vasculature, 0.06 6 0.03) and roots (epidermis, 0.21 6 0.09; vasculature, 0.17 6 0.08). No such differences were detected in the hypocotyl (epidermis, 0.14 6 0.05; vasculature, 0.10 6 0.07; Figures 4A and 4B).
By contrast, the NAD pools in 5-d-old etiolated seedlings grown in the dark under otherwise identical conditions showed a general shift toward oxidation. This shift was most pronounced in root tissues, which markedly decreased the gradient in the whole plant from shoot and root (Figures 4B and 4C). This reveals that even during early seedling establishment, when the photoautotrophic lifestyle is not yet fully functional, the status of photomorphogenesis and/or establishment of the photosynthetic machinery of the aerial organs has a strong impact on the metabolism of the heterotrophic root tissues.

Establishment of Parallelized Fluorimetry to Monitor Cytosolic NAD Redox Dynamics
The tissue-specific differences in NADH/NAD 1 ratios suggest that different ratios can be adopted by cells without crossing the limits of normal NAD redox physiology. To explore how volatile cytosolic NAD redox status really is, we aimed to study responses to external changes in a dynamic fashion using leaf tissue as a model. Plate reader-based multiwell fluorimetry allows reversible biosensor responses to be recorded over time, with parallelized controls and replicates, and is particularly well suited for monitoring the impact of external treatments at the tissue and organ level (De Col et al., 2017;Nietzel et al., 2019;Wagner et al., 2019). The ratiometric nature of the sensor is critical, since it provides internal normalization for variability in sensor expression levels as well as tissue amounts and positioning, but it requires the sensor signal to be sufficiently high in order to be reliably discriminated from the background. To test whether the Peredox-mCherry sensor lines pass this criterion, we assessed the degree of interference from endogenous fluorescence with the sensor signals in situ in leaf discs of 5-week-old Arabidopsis plants. Twodimensional spectra of excitation and emission from the three independent Peredox-mCherry lines compared to the wild type revealed the characteristic spectra of tS and mC (Supplemental Figure 5), which clearly resembled the in vitro spectra of the sensor ( Figure 2B). The sensor maxima were absent in the wild-type control, further indicating that autofluorescence was low in the relevant spectral regions. Importantly, chlorophyll fluorescence was spectrally well separated, and the emission peaks of tS and mC generally showed a favorable signal-to-noise ratio (e.g., 7

In Vivo Monitoring of NAD Redox Dynamics in Response to Changes in Illumination, Respiration, and Carbon Supply
We next aimed to explore the responses of cytosolic NAD redox status to external transitions with known effects on intracellular energy and redox physiology. To obtain insight into intercompartment redox exchange, we made use of the exclusive cytonuclear sensor localization and applied stimuli with primary effects on the chloroplasts and mitochondria, respectively ( Figure 5A). We first used illumination as a stimulus and performed dark-light transition experiments to trigger photosynthesis ( Figure 5B). Since fluorescence readings are not possible during illumination, we recorded tS/mC in the periods before and after light exposure. Interestingly, illumination led to a pronounced reduction of the cytosolic NAD pool, as indicated by the recovery of tS/mC (re-oxidation of the NAD pool) directly after the light treatment. The relatively slow dissociation kinetics of NADH from the Peredox-mCherry protein ( Figure 2E) likely underestimate the speed of the redox changes but is helpful for monitoring the transitions, where the true starting values in the light are not accessible. The presence of the photosystem II inhibitor 3-(3,4-dichlorophenyl)-1,1-dimethylurea (DCMU; Wilcox and Moreland, 1969) abolished the reduction of the cytosolic NAD pool by light, indicating that active photosynthetic electron transport in the (A) tS and mC fluorescence of 5-d-old seedlings was recorded by confocal laser scanning microscopy, and the ratio is plotted in false colors (sequential excitation, 405 and 570 nm; emission recorded at 517.5 6 7.5 and 612.5 6 7.5 nm, respectively), where high tS/mC (red) corresponds to reduced NADH/ NAD 1 steady state levels. The seedlings were kept in the dark for 120 min before image acquisition to avoid potential effects of active photosynthesis. A representative seedling of line #3 is shown. In the close-ups, tS/mC ratios were analyzed in the indicated ROIs, with 1 out of 10 randomly distributed circular ROIs across a cotyledon shown to facilitate visibility of the actual differences. Enlarged are I, cotyledon; II, hypocotyl and shoot apical meristem (SAM) region; III, root; and IV, root tip. Note that the SAM region was not exposed in the etiolated seedlings and was hence not analyzed to avoid mechanically invasive sample preparation. Bars 5 1 mm. cotyl., cotyledon; epid., epidermis; hypoc., hypocotyl; vasc., vasculature.
(B) Comparison of the tS/mC ratios of light-grown (A) versus etiolated seedlings in the dark (C), and comparison of tissue-specific differences in the seedlings under each condition. Seedlings of line #3 were used for the analysis; n 5 6, whiskers represent minimum and maximum values. Statistical analysis with a twoway ANOVA followed by the Tukey test; different letters indicate significant differences between ROIs of seedlings under the indicated growth condition (P # 0.05), asterisks indicate significant differences between light-grown versus etiolated seedlings for the indicated ROI (***P # 0.001; Supplemental Data Set 1F). (C) Confocal microscopy analysis of a representative 5-d-old seedling etiolated in the dark. For statistical comparison of the NAD redox state to light-grown seedlings (B), ROIs were selected analogously to (A). Bars 5 1 mm. thylakoids is causative for NAD reduction in the cytosol and implying that reducing equivalents are exported from the chloroplast to the cytosol.
We then investigated the effect of carbon status (in the form of free Suc) on cytosolic NAD redox status, as we hypothesized that the tissue-specific differences in the cytosolic NAD redox state ( Figure 4) between the vasculature, where sensor fluorescence predominantly comes from the phloem, and the epidermis were at least partly due to differences in sugar status. External supplementation with Suc induced a strong increase in the tS/mC ratio indicative of a reduction of the NAD pool ( Figure 5C), implying that the Suc was taken up, imported into the cell, and metabolized. The (A) to (F) Overview of the stimuli used to explore cytosolic NAD redox dynamics: illumination to activate photosynthesis including the chloroplastic electron transport chain (cETC) and cETC inhibition with photosystem II inhibitor DCMU (B), external Suc supplementation (C), inhibition of the mETC at complex III by AA (E), or combined stimuli (see [D] and [F]). Fluorescence emission was measured in each well, log 10 transformed, and averaged sensor ratios were plotted after autofluorescence correction using Col-0 leaves.
(B) to (F) Time-course measurements of Peredox-mCherry ratio (log 10 (tS/mC)) of leaf discs from 5-week-old plants. Arrows indicate treatment and control treatment (mock) application after at least 90-min dark adaption. Excitation at 400 6 10 nm (tS) and 570 6 10 nm (mC), emission at 515 6 7.5 nm (tS) and 610 6 5 nm (mC). n $ 5, means 1 SD. reductive impact of Suc feeding on the cytosolic NAD pool, as demonstrated here for leaf discs, may be accounted for by a boost in glycolysis and the TCA cycle. While photosynthetic activity and Suc feeding both boost the rate of reductant influx into the cytosolic NAD pool, we next tested the impact of decreasing the rate of reductant efflux by inhibiting mitochondrial electron transport as a major reductant sink. Inhibiting complex III using antimycin A (AA) led to the gradual reduction of NAD in the cytosol ( Figure 5E). The observation that inhibiting mitochondrial function leads to pronounced responses in the cytosol supports the concept that cytosolic NAD redox state is dependent on mitochondrial electron transport activity, both directly via external NADH dehydrogenases and indirectly by metabolite shuttles (Heber, 1974;Fernie et al., 2004;Selinski and Scheibe, 2019).
Similar NAD redox responses to those observed in leaf discs were also observed in 5-d-old seedlings (Supplemental Figures 6A  and 6B) and were reproducible between the different independent sensor lines (Supplemental Figure 6C), highlighting the robustness of the metabolic responses. The dynamics of the individual fluorescence signals of tS and mC further validate this technique. As expected from the sensor design concept, the mC signal remained largely stable over time, while tS was responsible for the transitions observed in the tS/mC ratio (Supplemental Figure 6D).
Combining light treatment with either Suc supplementation ( Figure 5D) or inhibition of mitochondrial electron transport ( Figure 5F) had a cumulative impact. Suc and AA induced gradual NAD reduction, which was only mildly increased further by illumination. Similarly, re-oxidation after illumination was abolished. Since this was the case at different stages of the gradual NAD reduction ( Figure 5D), it is unlikely that the abolished re-oxidation was caused by sensor saturation.
The reversible responses to external stimuli are suggestive of a remarkable degree of physiological flexibility of the cytosolic NAD pool, which is now accessible for specific measurement. Exclusive cytosolic NAD redox monitoring after the application of organelle-targeted stimuli demonstrates efficient reductant exchange between the subcellular compartments, which is likely due to the tight association of NAD redox status to the metabolic network and mediated by active metabolite shuttle systems.

Engineering a Sensor Variant for Monitoring More Reduced NAD Pools
The pronounced responses of Peredox-mCherry indicate that it is well suited to cover the physiologically relevant range of NAD redox states in the cytosol. However, the responses in leaves covered a large proportion of the full spectroscopic dynamic range of the sensor, indicating that it is likely to reach saturation at its maximal amplitude. Furthermore, reductive responses in tissues with already more reduced NAD pools are expected to be limited by the sensitivity of the sensor. To address this limitation and to expand the range in which NAD redox measurements can be performed, we aimed to engineer a sensor variant with a lower binding affinity for NADH for measuring higher NADH/NAD 1 ratios ( Figure 6). We chose a structure-informed rationale to destabilize NADH binding by replacing an Asp (D) with a Ser (S) in both Rex domains ( Figure 6A). This D is conserved in NAD-dependent oxidoreductases (Laurino et al., 2016) and stabilizes NAD binding through hydrogen bonding with a hydroxyl-group of the NAD Rib moiety. To generate the "Peredox-mCherry DS" variant, we used targeted mutagenesis of the sensor sequence and expressed the sensor variant recombinantly to assess its biochemical properties.
The DS variant showed decreased binding affinity to NADH, with a K d -value for NADH of 31.4 6 4.4 mM compared to 1.2 6 0.15 mM for Peredox-mCherry in the presence of 500 mM NAD 1 ( Figure 6B; Supplemental Figure 7A). The K d -value for NADPH (measured in the presence of 150 mM NADP 1 ) remained one order of magnitude higher at 324 6 18 mM NADPH. The spectroscopic dynamic range (d) was 2.88, that is, slightly higher than that of Peredox-mCherry (d 5 2.2). While the pH stability of the Peredox-mCherry sensor was maintained ( Figure 6C), the NADH binding and NADH dissociation kinetics were accelerated in the DS variant. Following the rationale of assessing the kinetics of Peredox-mCherry ( Figure 2E), we determined the binding kinetics of the DS variant at 500 mM NAD 1 and 30 mM NADH to ensure that both species bind to the sensor at a ratio of ;1:1. At NADH binding, tS/ mC of the DS variant required 26.5 6 4.5 s to increase to 50% of the full response amplitude (2.75 times faster than Peredox-mCherry at equivalent sensor occupation; Figure 6D). Replacement of NADH through the addition of NAD 1 (500 mM) took 1038 6 217 s to reach 50% of the response amplitude (4.9 times slower than for Peredox-mCherry, one-phase decay fit; Supplemental Data Set 1H).
The decreased (albeit not absent) NAD 1 response of the NADHbound DS variant in terms of both speed and impact on sensor ratio prompted us to assess the NAD 1 dependence of both Peredox-mCherry and Peredox-mCherry DS in further detail (Supplemental Figures 7B to 7E). Interestingly, these titrations revealed that NAD 1 concentration had a much more minor impact on the sensor response to NADH of the DS variant compared to Peredox-mCherry. The reduced NAD 1 sensitivity may be explained by the replacement of the negatively charged D by a neutral S at the binding pocket, by which charge-based stabilization of NAD 1 binding appears to be lost. As a result, the DS variant is mainly responsive to NADH concentration instead of NADH/NAD 1 ratio. Under physiologically meaningful conditions, where NAD 1 concentrations are stable compared to NADH, the increased K d remains the dominant change in the DS variant compared to Peredox-mCherry. The decreased affinity of Peredox-mCherry DS is expected to allow more reduced NAD pools to be investigated compared to Peredox-mCherry, which can be critical depending on external conditions, tissue, and subcellular compartment. Specifically, with respect to the data presented here, we asked whether the DS variant may allow for discrimination between biological and technical causes when tS/ mC ratio kinetics reach a plateau at high ratios (reduced NAD pool), as illustrated in Figures 5C and 5E.

Peredox-mCherry DS Expands the Response Range of in Vivo NAD Redox Sensing
To explore the response characteristics of the Peredox-mCherry DS variant, we generated Arabidopsis lines constitutively expressing Peredox-mCherry DS in the cytosol and nuclei (Supplemental Figures 8A to 8I). We selected three independent lines based on their bright fluorescence signal. Two-dimensional excitation-emission scans confirmed that the signals of tS and mC were sufficiently intense for measurements (Supplemental Figures  8J to 8L). All three sensor lines (#7 and #14 homozygous, #1 segregating) showed a wild-type-like phenotype (Supplemental Figure 8M). We then compared the responses of Arabidopsis leaf discs from 5-week-old plants expressing either Peredox-mCherry or Peredox-mCherry DS to illumination, the inhibition of respiration, and external Suc supplementation side by side by multiwell plate reader-based fluorimetry (Figure 7). Prior to any treatment, a significant difference in tS/mC steady state was consistently observed between tissues expressing the two sensor variants (Figure 7; Supplemental Figure 9A). A shift to lower ratios by the Peredox-mCherry DS variant could be expected due to the lower NADH affinity of the sensor ( Figure 6B). The tS/mC ratio of Peredox-mCherry DS showed a diminished response to illumination ( Figure 7A). With increasing light exposure times, a minor post-illumination reduction transient was observed. By contrast, inhibiting respiration with AA resulted in an increase in the tS/mC ratio of Peredox-mCherry DS, albeit at a lower rate than Peredox-mCherry ( Figure 7B). As a result, the tS/mC values of DS variant only reached a plateau after 2.3 h, that is, 2.7-fold slower than Peredox-mCherry. This difference is in line with the higher K d of the DS variant and suggests that the plateau reached by Peredox-mCherry represented saturation of the sensor rather than physiological stabilization of the NADH/NAD 1 ratio.
As an alternative approach to interfere with NADH oxidation by the mitochondria, we used low oxygen (2.5%) to induce hypoxia ( Figure 7C; Wagner et al., 2019). At decreasing oxygen concentrations, the cytosolic NAD redox state, as monitored by Peredox-mCherry, rapidly shifted toward reduction. Once the oxygen concentration was stabilized at 2.5%, reduction progressed more slowly. Re-establishment of normoxia triggered the gradual reoxidation of NAD. The response dynamics of Peredox-mCherry and Peredox-mCherry DS were very similar at the qualitative level, (B) Ratiometric response of the Peredox-mCherry DS variant compared to Peredox-mCherry (log 10 (tS/mC)) at different NADH and NADPH concentrations in the presence of 500 mM NAD 1 and 150 mM NADP 1 , respectively. Sensor protein concentration: 0.025 mg/mL, n 5 3, mean 6 SD, pH 7.5. The ratio values were normalized to the highest ratio (set to 1). Dotted lines indicate K d -values (sigmoidal curve fitting with details provided in Supplemental Data Set 1G). (C) Response of the DS sensor variant to different pH values. Sensor protein was equilibrated with three different NADH concentrations (1.5, 30, and 200 mM) in the presence of 500 mM NAD 1 to achieve low, mean, and high NADH occupancy at different pH values (6.0 to 9.0, n 5 3, mean 6 SD; Bis-Tris for pH 6.0 to 7.0 and Tris-HCl for pH 7.5 to 9.0). (D) NADH binding to and dissociation from Peredox-mCherry DS in vitro. Time resolved tS/mC (magenta, left y axis) and corresponding log 10 ratio (dark purple, right y axis) changes in recombinant sensor protein (0.025 mg/mL, Tris-HCl, pH 7.5) in response to NADH (30 mM, first arrow) in the presence of 500 mM NAD 1 and to further NAD 1 addition (500 mM, second arrow). Dotted lines indicate time to reach half of the response amplitude at NADH or NAD 1 addition (one-phase decay fit; Supplemental Data Set 1H). Fluorescence emission was measured every 13 s. Excitation, 400 6 10 nm (tS) and 570 6 10 nm (mC); emission, 520 6 5 nm (tS) and 610 6 5 nm (mC). The assay was independently repeated three times with consistent results. but the relative tS/mC changes were less strong for DS, as expected since it has a lower K d . Both sensors resolved the NAD redox response dynamics when NAD was shifted toward reduction, without sensor saturation. Comparing the impact of Suc supplementation on the in vivo responses of Peredox-mCherry ( Figure 7D) and Peredox-mCherry DS ( Figure 7E) further validated the characteristics of the DS variant. Suc induced a concentrationdependent increase in the tS/mC of both Peredox-mCherry and Peredox-mCherry DS, with a log 10 -response range of ;0.4 (d of ;2.5) for both two sensor variants, covering their entire spectroscopic dynamic range. The rate of tS/mC increase was concentration dependent up to 10% (w/v) Suc. The plateau of saturated tS/mC values was reached by Peredox-mCherry after 11 h (1%), 2.7 h (5%), and 1.3 h (10%), and reaching the plateau took ;2.7 times longer for the DS variant (28.2 h for 1%, 7.5 h for 5%, 3.6 h for 10%, one-phase decay fit of ratio values from time point of treatment; Supplemental Data Set 1I).
Unexpectedly, we also observed a pronounced impact of elicitor exposure on cytosolic NAD redox dynamics ( Figure 7F). Monitoring cytosolic NAD redox dynamics of leaf discs in response to the bacterial elicitor flg22 revealed a rapid and transient reduction, the timing of which appeared to correlate with stimulation of plant NADPH oxidase (respiratory burst oxidase homolog [RBOH]) activity (Nietzel et al., 2019). The characteristic and rapid dynamics allowed for a side-by-side comparison of Peredox-mCherry and the DS variant in vivo. Peredox-mCherry responded with an initial linear increase in log 10 (tS/mC) ratio of 0.51 6 0.11 per hour, whereas the response rate by Peredox-mCherry DS was nearly twice as high (0.95 6 0.31 per hour; Figure 7Fi). That is in line with the faster in vitro NADH binding rate of the DS variant. The amplitude of the DS variant was much lower than that of Peredox-mCherry (0.105 6 0.011 versus 0.177 6 0.034; Figure 7Fii), which is consistent with a lower affinity for NADH. Strikingly, the DS variant reached the maximum of the transient much earlier than Peredox-mCherry (18.6 6 2.02 versus 41.38 6 1.48 min; Figure 7Fiii), which probably results from a combination of the different binding and dissociation kinetics as well as the different affinities for NADH and the differential NAD 1 sensitivities. Redox responses were reproducible between the different independent DS sensor lines (Supplemental Figures 9B to 9D), highlighting the robustness of the metabolic responses. The similarity of the dynamics of Peredox-mCherry DS and Peredox-mCherry at the qualitative level ( Figures 7A to 7E) validate faithful responses with the same in vivo specificity by both sensors.

Understanding in Vivo Metabolism Requires a Dynamic Understanding of NAD Redox Status
Monitoring intracellular NAD redox dynamics in planta has been a long-standing challenge. Established approaches to measure NADH and NAD 1 concentrations or NADH/NAD 1 ratios by performing enzymatic assays, analytical methods from whole-cell extracts, inference from marker metabolites, or in vivo using NAD(P)H autofluorescence are limited in their spatial or temporal resolution (Williamson et al., 1967;Queval and Noctor, 2007;Blacker et al., 2014). Analytical methods such as HPLC analysis can offer high accuracy and absolute quantitation, at the cost of mixing subcellular pools (Lu et al., 2018). An exception has been subcellular fractionation of pea leaf protoplasts (Igamberdiev and Gardeström, 2003), which indeed allowed subcellular NAD pools to be resolved and their changes measured in response to modified light and CO 2 supply. While this method remains technically challenging, which has prevented its broader adoption, the recent refinement of nonaqueous fractionation techniques may hold promise for performing analytical NAD measurements in nonprotoplast systems (Fürtauer et al., 2016;Beshir et al., 2019;Medeiros et al., 2019). Although rapid changes are an integral feature of plant metabolism, including redox metabolism (Eberhard et al., 2008), it has remained extremely challenging to capture NAD redox changes with a temporal resolution that matches these metabolic fluctuations. As a consequence, direct evidence remains scarce about how different cell compartments are differentially and dynamically affected by metabolic fluctuations as induced, for instance, by light-to-dark transitions. Recent NAD redox measurements in rapidly isolated mitochondria allowed absolute as well as relative NAD redox states to be determined, but rapid organelle isolation approaches cannot be analogously implemented for the cytosol (Luo et al., 2020). Furthermore, it is difficult to quantify the impact resulting from the organelle isolation procedure by removing organelles from the intact plant system. In vivo autofluorescence measurements of NADH and NADPH allow for temporal resolution but lack the specificity to distinguish NADH from NADPH. In addition, these measurements are technically challenging in plants due to interfering autofluorescence from other cellular compounds in the short wavelength range (Shaw and Ehrhardt, 2013). Yet, it is evident that an understanding of NAD redox dynamics, as they occur in the living plant at a tissue, cellular, and subcellular scales, will be critical if we are to integrate our mechanistic understanding about NAD biochemistry into a meaningful physiological framework in planta.

Peredox-mCherry Sensor Lines Allow NAD Redox Dynamics to Be Analyzed in the Plant Cytosol
The Arabidopsis sensor lines that we introduce here show constitutive cytosolic expression ( Figure 1) and are not affected by silencing effects, as judged from the regular use of these lines for experimentation for up to 7 years. Unambiguous targeting and stable expression is in contrast to observations for several other sensors associated with severe silencing and subcellular mistargeting (Albrecht et al., 2011;Behera et al., 2015;Schwarzländer et al., 2016). Despite the successful expression in the plant cytosol and nucleoplasm, efforts to target Peredox-mCherry into mitochondria with three different constructs did not lead to the generation of any transformants, suggesting that the size (103.5 kD) and structure of the biosensor interfere with mitochondrial protein import, which leads to embryo lethality. This interpretation remains speculative but is consistent with the phenotypes of plants expressing other mitochondrial sensors (De Col et al., 2017;Nietzel et al., 2019). Peredox-mCherry is significantly larger than any of the sensors that we have found to cause phenotypes, and it contains two fluorescent protein domains, which may be able to arrest import if folding and/or chromophore maturation takes   Side-by-side multiwell plate reader-based time-course measurements of Peredox-mCherry and Peredox-mCherry DS ratio (log 10 (tS/mC)) of leaf discs from 5-week-old Arabidopsis plants dark adapted for at least 90 min prior to treatments. Excitation at 400 6 10 nm (tS) and 570 6 10 nm (mC); emission at 515 6 7.5 nm (tS) and 610 6 5 nm (mC). n $ 5, means 1 SD. place before import is completed. Importantly, expression of Peredox-mCherry in the plant cytosol did not induce growth phenotypes (Figure 1). This is consistent with the idea that NAD binding to the sensor in addition to the many endogenous NAD binding proteins does not disrupt cytosolic NAD physiology.

Biochemical Characteristics of Peredox-mCherry Fit Plant Cytosolic Physiology
We selected Peredox-mCherry over other NAD biosensors, as it is the only sensor that reports changes in the NAD redox status rather than changes in NADH or NAD 1 concentration, while providing a ratiometric readout and being largely pH inert (Figure 2; Hung et al., 2011). In vivo, the NAD redox state is predominantly set by the activity of the metabolic oxidoreductases, some of which can carry large fluxes and use NAD as a cofactor. De novo NAD 1 synthesis, subcellular transport, recycling through the salvage pathway, degradation and utilization as a substrate, and the conversion of NAD(H) to NADP(H) by NAD(H) kinases (and vice versa by phosphatases) are generally less rapid and are expected to have a much smaller impact on NAD redox status. Instead, they can adjust subcellular NAD pool size Kawai and Murata, 2008;Waller et al., 2010). As such, it is the NADH/ NAD 1 ratio, as opposed to the concentrations of NADH or NAD 1 , that is the decisive parameter to measure when exploring dynamics in redox metabolism. Since both NADH and NAD 1 bind to Peredox-mCherry and Peredox-mCherry DS, but only NADH binding affects tS fluorescence, both sensors report NADH/NAD 1 ratios, while changes in total NAD 1 and NADH concentrations cannot be monitored (Hung et al., 2011). Previous estimations based on rapid subcellular fractionation of pea leaf protoplasts indicated an NADH/NAD 1 ratio between 0.03 to 0.2, with NAD 1 at ;500 mM and NADH at ;20 mM (Igamberdiev and Gardeström, 2003). Considering the difference in concentration between both NAD species, any ratio change by NAD 1 -NADH interconversion will be dominated by the change in NADH concentration, while NAD 1 concentration will have only a minor influence. Moreover, the sensor provides a readout only for the free NAD pool, which remains largely unchanged regardless of the content of NAD bound to proteins (Kasimova et al., 2006). A substantial fraction of the total NAD pool was estimated to be bound to proteins in plant mitochondria (e.g., 30 to 50% of NAD 1 and 75% of NADH; Kasimova et al., 2006), similar to mammalian erythrocytes (50% free NAD; Canepa et al., 1991). Although the NAD redox state is independent of total NAD pool size (in contrast to, e.g., glutathione), preferential sequestration or release of NAD 1 or NADH between the bound and free NAD pools may still lead to considerable redox changes of the free pool. Combinatorial analyses with sensors that monitor NAD 1 concentration (Cambronne et al., 2016) have the potential to further uncover the interplay between in vivo NAD pool sizes and redox status in the future.
An important characteristic of Peredox-mCherry is the overall stability of its readout under physiological pH changes (Figure 2). The impact of pH changes on the response of fluorescent biosensors remains a widespread problem that has previously led to ambiguous interpretations (Schwarzländer et al., 2011(Schwarzländer et al., , 2012(Schwarzländer et al., , 2014. The issue is increasingly recognized and has stimulated the development of pH-stable sensor variants and dedicated correction protocols (Tao et al., 2017;Nietzel et al., 2019;Pak et al., 2020). Also, the development of Peredox-mCherry included a dedicated mutagenesis and screening step to identify a pH inert sensor (Hung et al., 2011), while alternative sensors for NAD redox dynamics, such as SoNar, show pronounced pH dependence (Zhao et al., 2015). pH stability is particularly relevant for the use of sensors in plants under conditions known to trigger intracellular pH changes, such as hypoxia, under which cytosolic pH can drop by more than 1 pH unit (Geigenberger, 2003;Wagner et al., 2019). However, pH changes are also relevant under nonstress conditions and require careful consideration. Physiological events, such as Ca 21 transients or change in illumination status, can trigger pH changes in the range of 0.1 to 0.5 units (Yin et al., 1990;Mühling et al., 1995;Behera et al., 2018), which bears a considerable risk of generating pH-related sensing artifacts when pHresponsive biosensor variants are applied without caution and robust controls. Although cytosolic pH changes are highly likely to occur in response to the illumination, low oxygen, and respiratory inhibition treatments applied here, these changes are unlikely to have a pronounced impact of the response of the Peredox-mCherry sensors.

Toward an Integrated Understanding of Tissue-Specific NAD Redox Metabolism
Direct in vivo mapping of the cytosolic NAD redox state across an entire Arabidopsis seedling by confocal microscopy resolved a remarkable degree of heterogeneity between different tissues and cells (Figure 4). Green tissues showed more oxidized NAD redox states compared to heterotrophic tissues with more reduced NAD redox states. Since the seedlings were dark adapted before imaging, the more oxidized NAD pool in photosynthetically active tissues may reflect its capacity to accept electrons exported by the chloroplast in the light, and indeed illumination of leaf discs led to NAD pool reduction in the cytosol ( Figure 5B; Figure 7A). Interestingly, the shoot-to-root redox gradient was markedly decreased in etiolated seedlings grown in the dark, mainly due to a more oxidized NAD pool in the root tissues ( Figure 4). This suggests that photo-morphogenesis affects the redox metabolism of heterotrophic root tissues in particular. Superimposing the NAD redox seedling atlas with a recently reported cytosolic MgATP 2map (De Col et al., 2017) revealed a correlation between ATP concentration and NAD oxidation, which were generally both high in shoot tissues and both low in root tissues. Also, the shoot-root gradients of both ATP and NAD redox state were flattened in etiolated seedlings, raising the possibility that the correlation is underpinned mechanistically by metabolism, most likely respiration. This overall correlation was not observed in the SAM region, which displayed both low ATP levels and a particularly oxidized NAD pool. Low ATP levels are consistent with respiratory inhibition resulting from a hypoxic environment, which was recently found to regulate the development of leaf primordia (Weits et al., 2019). However, hypoxia is typically associated with NAD reduction, rather than oxidation (Figure 7; Chance et al., 1962;Wagner et al., 2019), raising intriguing questions about the specific metabolic programs used by cells in the meristem and leaf primordia.
Furthermore, cytosolic NAD generally appeared to be more oxidized in the vasculature than in the surrounding tissues. The sensor fluorescence signal in the vasculature most likely originates mainly from the phloem, which shows elevated Suc levels, high rates of oxygen consumption, and few intercellular air spaces and is prone to hypoxia (van Dongen et al., 2003). Like the shoot meristem region, a particularly reduced rather than oxidized NAD pool may be expected in the vasculature, considering that both increased Suc supply and hypoxia induced NAD reduction in leaf discs (Figures 5 and 7;Zhang et al., 2017). The possibility of a highly specialized phloem metabolism deserves further attention in the future, aided by sensor expression under cell type-specific promoters, such as the SUCROSE-PROTON SYMPORTER2 promoter for phloem companion cells, in analogy to recent advances in cell type-specific transcriptome analyses and even organelle isolations (Boussardon et al., 2020;Rich-Griffin et al., 2020).

Cytosolic NAD Redox Dynamics Are Remarkably Flexible and Integrate Organelle Redox Metabolism
In addition to the dependence of NAD redox status on developmental status, online sensor monitoring revealed remarkable flexibility of the cytosolic NAD redox state, as apparent from the study of dynamic NAD redox responses in leaf tissue (Figures 5  and 7). While large fluctuations, rather than strict buffering, may seem counterintuitive for a central hub of the metabolic network, changes in cytosolic NADH/NAD 1 ratio are in general agreement with the results of previous analytical single time point measurements, for example, under different light and CO 2 conditions (Igamberdiev and Gardeström, 2003).
The reduction of cytosolic NAD in response to external Suc ( Figures 5C, 7, and 7E) is in line with the efficient uptake of Suc into the cytosol and its breakdown by glycolysis and the mitochondrial TCA cycle, which is connected to the cytosol by metabolite shuttles. Even though the rate of NADH production likely increases in response to external Suc, the reduction of the NAD pool results from a shift in the balance between the rates of NADH production and NADH oxidation, suggesting that mitochondrial electron transport is unable to fully compensate for the increased flux through primary metabolism. This interpretation is validated by the increased NADH/NAD 1 ratio when mitochondrial electron transport capacity is constrained, as is the case under AA inhibition or hypoxia ( Figure 5E, 7B, and 7C). The cytosolic response further emphasizes the notion that the mETC acts as major electron sink for cytosolic NADH, and the activity of metabolite shuttles and the external NADH dehydrogenases represent obvious links (Geigenberger, 2003;Rasmusson et al., 2004).
The responses of the sensors to illumination and respiratory inhibition provide direct in vivo evidence that the cytosol is central to integrating redox metabolism from the organelles. Since the in vivo measurements integrate all operational crosstalk mechanisms-known and unknown-and consider their relative quantitative contributions, these observations not only validate the wealth of previous biochemical measurements but also allow comparative measurements in response to different and combined external stimuli (Figures 5 and 7).
The activation of photosynthesis by illumination leads to a reduction of the cytosolic NAD pool, as indicated by the return of the sensor ratio to steady state in the dark (Figures 5B and 7A). The export of photosynthetic reductant from the chloroplast via the chloroplastic malate valve and triose phosphate-3-phosphoglycerate shuttle has been studied intensely and can prevent the photoinhibition of photosynthetic electron transport during illumination (Heber, 1974;Taniguchi and Miyake, 2012;Selinski and Scheibe, 2019;Shameer et al., 2019). Efficient electron shuttling through the malate valve depends on malate-oxaloacetate exchange across the chloroplast inner envelope, the light/thioredoxin-activated plastidic NADP-MDH, and cytosolic NAD-MDH (Wolosiuk et al., 1977;Scheibe and Jacquot, 1983;Hatch et al., 1984;Kinoshita et al., 2011). As such, the exported reductant from the NADP pool in the chloroplast stroma can enter the NAD pool in the cytosol, emphasizing the need to monitor cytosolic NAD rather than NADP in the context of photosynthetic redox metabolism. It should be emphasized, however, that the sensor response indicates net changes in NAD redox status, while the underlying gross electron flux rates cannot be resolved. It would be desirable to investigate NAD and NADP redox dynamics simultaneously in the future. A family of NADPH sensors was recently introduced, but caution should be taken before applying these sensors to studying plant NADPH dynamics and photosynthesis due to their pronounced inherent pH sensitivity (Tao et al., 2017).
While this work was under revision, Lim et al. (2020) reported the use of the NAD redox sensor SoNar and the NADPH sensor family iNap to assess the subcellular responses of redox metabolism to photosynthetic activity in the cytosol, chloroplast stroma, and peroxisomes of Arabidopsis cotyledons. The authors performed pH correction using control sensor constructs without pyrimidine nucleotide binding capacity expressed in separate control plants measured in separate experiments. Consistent with our observations based on Peredox-mCherry, the SoNar sensor lines did not show any obvious developmental phenotypes, and the measurements using the cytosolic SoNar sensor indicated a more oxidized NAD pool in roots than in shoots, as well as pronounced NAD pool reduction in response to mitochondrial electron transport inhibitors, such as AA. While the pH-corrected SoNar data indicated strong and reversible reduction of the stromal NAD pool in response to illumination, cytosolic NAD reduction was detectable, but the sensor response was subtle at light intensities that were only slightly lower than those applied here (296 versus 400 mmol m 22 s 21 ). This apparent difference may be due to the different NAD binding and NAD response characteristics between Peredox-mCherry and SoNar, the required pH correction regime for SoNar, the difference in light intensity, and/or the use of cotyledons by Lim et al. (2020) compared to the true leaf tissue of 5week-old Arabidopsis rosettes investigated here.
The mechanisms that underpin photosynthetic metabolism and redox dynamics between cell compartments likely differ between cotyledons and true leaves. The cumulative effects of illumination and respiratory inhibition clearly demonstrate that redox crosstalk between chloroplasts and mitochondria is integrated by the cytosolic NAD pool in true leaf tissue from 5-week-old rosettes ( Figure 5F). This system provides a platform to experimentally dissect the relative significance of major processes in organelle energy metabolism and their impact on inter-organelle redox crosstalk, including that of cyclic electron flow, the photorespiratory pathway (as initiated by Lim et al. [2020] using a pharmacological approach), and the alternative respiratory pathway. We have initiated reverse genetic approaches to manipulate chloroplast-cytosol-mitochondria interactions at various steps of metabolism and transport and expect online NAD redox monitoring to provide novel insight into their operation under changing external conditions in the future.
While the observed NAD redox dynamics in response to Suc, respiratory inhibition, hypoxia, and illumination ( Figures 5 and 7) validate established mechanisms of intracellular redox metabolism, the observation of an elicitor-induced NAD redox transient came as a surprise and currently lacks a defined mechanistic foundation ( Figure 7F). Since the transient is reminiscent of Ca 21 and reactive oxygen species transients that have been extensively investigated in the context of elicitor-induced NADPH oxidase/ RBOH activity, the mechanism of NAD reduction may be directly or indirectly linked to RBOH activity as well. Recent work has implicated a calmodulin/Ca 21 -dependent NAD kinase in providing NADP 1 to support RBOH activity (Dell'Aglio et al., 2019). The phosphorylation of large amounts of NAD 1 at a high rate may indeed result in a reductive shift of the cytosolic NAD pool, even though any quantitative estimation of the actual impact remains highly speculative. Consistently, a knockdown mutant impaired in the first step of de novo NAD biosynthesis shows deficient flg22triggered ROS production by RBOHD and reduced stomatal immunity (Macho et al., 2012). Reduction of the NAD pool by NAD 1 depletion may further be due to the elicitor-induced NADase activity of the Toll/interleukin-1 receptor homology (TIR) domains of TNL proteins, which degrade NAD 1 as an essential step in signaling between TNLs and ENHANCED DISEASE SUSCEPTI-BILITY1 (EDS1; Wiermer et al. 2005;Wagner et al. 2013) to induce cell death (Horsefield et al. 2019;Wan et al. 2019). Moreover, a recent study established a link between impaired NAD biosynthesis and abscisic acid sensitivity, as mediated by SnRK2 protein kinases, RBOHF, and increased rates of reactive oxygen species production (Hong et al., 2020); high NADH/NAD 1 ratios in an NAD biosynthesis mutant caused increased RBOHF activity. Similarly, we observed a reduction of the NAD pool in response to the flg22-induced RBOHD/F-mediated respiratory burst (Kadota et al., 2014). In line with these results, defense-induced RBOHF activity has been implicated in pronounced metabolomic changes (Chaouch et al., 2012). Exactly how changes in NAD redox state and RBOH activity are mechanistically linked deserves further attention.

Peredox-mCherry DS Variant Expands the Measurement Range for NAD Redox Dynamics
The pronounced flexibility of NAD redox dynamics poses the challenge that a single sensor variant cannot cover the full spectrum of physiological responses. By performing structureinformed targeted mutagenesis, we generated a Peredox-mCherry DS variant that is able to resolve a higher NADH concentration range while showing decreased NAD 1 binding ( Figure 6; Supplemental Figure 7). Even though the DS variant indicates NADH concentration rather than NADH/NAD 1 , the physiological responses of both sensors can be compared, and the DS variant reports NAD redox in vivo under the assumption that changes in NAD 1 are minor compared with changes in NADH (which is justified in the cytosol under physiological conditions). When starting from a comparatively oxidized NAD pool, as the case in Arabidopsis leaves, the response to NAD reduction was slower for the DS variant compared to Peredox-mCherry, while the tS/mC values of the DS variant continued to increase in the presence of a more reduced NAD pool, when those of Peredox-mCherry have already reached saturation ( Figures 7B to 7E). This behavior is in line with the lower NADH binding affinity of the DS variant, as observed in vitro ( Figure 6B). On the other hand, the DS variant responded to elicitor exposure at a higher rate than Peredox-mCherry, indicating that in the initial phase of the rapid transient, its faster NADH binding kinetics appear to dominate the response rate ( Figure 7F). As Peredox-mCherry can reach saturation under conditions that are physiologically meaningful in the plant cytosol (Figures 4 and 5), the DS variant is able to respond to more reduced NAD redox states. Conversely, Peredox-mCherry DS is limited in resolving the dynamics of more oxidized NAD redox states, as indicated by the light response ( Figure 7A) and the lower amplitude of the flg22-induced transient ( Figure 7F), suggesting that the DS variant is fully devoid of bound NADH in the dark. As a result, the combination of both sensors expands the range of NAD redox dynamics than can be resolved, and their respective use can be matched to the expected NAD redox dynamics of the specific physiological scenario to be investigated. If the problem of mitochondrial targeting can be resolved in the future, the DS variant might be a good candidate to monitor NAD redox dynamics in mitochondria, which are estimated to display 100-to 1000-fold higher NADH/NAD 1 ratios than the cytosol (Williamson et al., 1967;Igamberdiev and Gardeström, 2003;Hung et al., 2011;Luo et al., 2020).
In principle, in vivo monitoring using two sensors with different affinities and their side-by-side calibration by driving the in situ NAD redox states to their extremes would allow spectroscopic sensor ratios to be converted to absolute NADH/NAD 1 ratios (and NADH concentrations) and NAD redox potential, as implemented, for example, for Ca 21 and glutathione redox sensing (Palmer and Tsien, 2006;Schwarzländer et al., 2008;Aller et al., 2013). However, any absolute quantitation of fluorescent protein sensor data should be interpreted with the utmost caution, since the conversion requires additional assumptions that are prone to error, resulting in the misleading representation of sensor data. Nevertheless, with this in mind, plausible estimates can be made. Previous estimations based on rapid subcellular fractionation of pea leaf protoplasts have indicated an NADH/NAD 1 ratio between 0.03 to 0.1, with NAD 1 at ;500 mM and NADH at ;20 mM (Igamberdiev and Gardeström, 2003); similar NADH/NAD 1 ratios were obtained from cellular extracts of cucumber (Cucumis sativus) leaves (Szal et al., 2008). Ratios of ;0.001, that is, at least one order of magnitude lower, are regularly determined in the cytosol of mammalian tissues (Williamson et al., 1967), suggesting that plants and mammals maintain a fundamentally different cytosolic NAD redox status or that the measurements in subcellular extracts still overestimate even lower in vivo ratios. The finding that Peredox-mCherry (K d for NADH of 1.2 mM in the presence of 500 mM NAD 1 , with a resulting K ratio of 0.0024) responds to in vivo NAD reduction in the Arabidopsis cytosol, rather than being saturated, strongly suggests that NADH/NAD 1 ratios are in the order of 0.001, that is, similar to those determined in the mammalian cytosol. Assuming that cytosolic NAD 1 concentrations from subcellular fractionation of pea leaf protoplasts (570 6 30 mM in the dark; Igamberdiev and Gardeström, 2003) are also a good approximation for the Arabidopsis leaf (mesophyll) cytosol, free NADH concentrations of ;0.5 to 1 mM are estimated at steady state in the dark based on the biosensors (Figures 6 and 7), as opposed to 18 6 14 mM NADH in cytosolic fractions from pea mesophyll protoplasts in the dark. These lower NADH values are consistent with other determinations from plant cell extracts, although these analyses lacked subcellular resolution (Heineke et al., 1991;Wigge et al., 1993).

Future Perspectives for NAD Redox Sensing in Plants
This work explored in vivo NAD redox dynamics in Arabidopsis. We carefully established an approach using a NAD sensor and showcased the power of this approach for plant research, while clearly just scratching the surface of its full potential. The ability to monitor NAD redox dynamics live and in cells of the intact living organism opens a methodological door to answering several longstanding and pressing questions. For instance, insights into the relationship between cellular redox metabolism and developmental programs will be in reach (de Souza Chaves et al., 2019;Feitosa-Araujo et al., 2020). Furthermore, the spatiotemporal signatures of NAD redox dynamics will shed light on the role of NAD under stress conditions or during immune responses (Zhang and Mou, 2009;Hashida et al., 2010;Pétriacq et al., 2012Pétriacq et al., , 2016Horsefield et al., 2019;Wan et al., 2019). Well-established concepts can also be put to the test in vivo, including the significance of alternative respiratory activity for subcellular redox metabolism (Michalecka et al., 2003;Rasmusson et al., 2008) or the differential roles of different metabolic redox shuttle systems (Selinski and Scheibe, 2019;Shameer et al., 2019). Biosensor expression driven by cell type-specific promoters will facilitate the analysis of cell type-specific metabolic modes and refine the NAD redox atlas generated here by global sensor expression. A particularly attractive route will be to obtain a systems perspective that links metabolic heterogeneity and plant development by integrating the cell-specific NAD redox data with data from other biosensors (e.g., for ATP or plant hormone distribution; Brunoud et al., 2012;Jones et al., 2014;Waadt et al., 2014;De Col et al., 2017;Rizza et al., 2017), but also with cell-specific microgenomic data, as well as single cell proteome and metabolite profiles once they become available (Bailey-Serres, 2013;Li et al., 2016). Finally, the optimization of NAD sensing for the apoplast, that is, the generation of high-affinity sensor variants that remain usable at very low pH, will open new avenues to understand the physiological context of extracellular NAD signaling (Wang et al., 2017(Wang et al., , 2019.

Generation of Expression Constructs and Plant Lines
The coding sequence of Peredox-mCherry was PCR amplified from GW1-Peredox-mCherry (Hung et al., 2011), inserted into pDONR207 (Invitrogen), and transferred to pSS02 (derivative of pMDC32; Curtis and Grossniklaus, 2003) as well as pETG10A by Gateway cloning (Invitrogen). pSS02 allows for the constitutive expression of genes in plants under the control of the UBIQUITIN10 promoter; pETG10A allows for inducible expression of recombinant protein in Escherichia coli cells. The plasmid GW1-Peredox-mCherry (Hung et al., 2011) was obtained from Addgene (plasmid no. 32380; http://www.addgene.org/32380/). The DS variant was generated by site-directed mutagenesis. Primers used to change the two codons are listed in the Supplemental Table. Agrobacterium (Agrobacterium tumefaciens)-mediated transformation of Arabidopsis (Arabidopsis thaliana; accession Col-0) was performed by the floral dip method (Clough and Bent, 1998). Transformants and homozygous lines were selected based on the resistance marker for hygromycin and sensor fluorescence.

Plant Culture and Phenotyping
Arabidopsis seedlings were grown from surface-sterilized seeds on vertical plates containing half-strength Murashige and Skoog (MS) medium (Murashige and Skoog, 1962) with 10 mM 2-(N-morpholino)ethanesulfonic acid (MES buffer), pH 5.8 (KOH), supplemented with 1% (w/v) Suc and 1% (w/v) agar and cultivated under long-day conditions (16 h, 75 to 100 mmol photons m 22 s 21 using Osram L 18W 840 Lumilux cool-white or Philips TL5 HO 49W 840 [MASTER] cool-white tubes at 22°C, 8 h dark at 18°C) after stratification for 2 d at 4°C in the dark. Primary root length was documented from seedlings growing side by side with their corresponding controls and measured using ImageJ (n 5 25). To analyze leaf rosette development and measure primary inflorescence height, plants were individually transferred to soil after 5 d and regularly documented photographically. Leaf rosette size was analyzed using the custom Leaf Lab tool (version 1.41) as previously described by Wagner et al. (2015), and height of the primary inflorescence was quantified using ImageJ (n 5 19 to 20). To prepare leaf discs, plants were germinated and grown in soil under long-day conditions.

Purification of Recombinant Sensor Proteins
Protein expression and purification were performed according to De Col et al. (2017) with the following modifications: E. coli strain BL21 (DE3) ArcticExpress cells carrying pETG10A:Peredox-mCherry or pETG10A:-Peredox-mCherry DS were selected for overexpression based on an optimization assay (Supplemental Figure 2). Overexpression was induced by isopropyl b-D-1-thiogalactopyranoside at a final concentration of 0.1 mM. Cells were collected by centrifugation at 5000g for 15 min at 4°C, and the pellet was resuspended in Tris buffer (50 mM Tris-HCl, pH 7.5, 100 mM NaCl, and 0.5 mM MgCl 2 ) supplemented with 1 mg/mL lysozyme, 0.1 mg/ mL DNaseI (Roche), and complete protease inhibitor cocktail (Roche). After a 30-min incubation on ice, the cells were sonicated (3 3 2 min, 50% power output, 50% duty cycle). The lysate was centrifuged at 40,000g for 30 min at 4°C, and the supernatant was loaded onto a Ni 21 nitrilotriacetic acid Hi-sTrap column (GE Healthcare). Proteins were eluted with an imidazole gradient (10 to 250 mM in Tris buffer) using an ÄKTAprime plus system (GE Healthcare). Fractions containing recombinant protein were pooled, concentrated by ultrafiltration, and applied to a HiLoad 16/600 Superdex 200 column (GE Healthcare) pre-equilibrated with Tris buffer. Concentrated Peredox-mCherry and Peredox-mCherry DS supplemented with 20% (v/v) glycerol were stored at -80°C.

Characterization of Sensor Protein in Vitro
The concentration of purified recombinant sensor protein was quantified according to Bradford (1976). Protein at a final concentration of 0.25 mM was mixed with Tris buffer for all in vitro assays, except for the pH assays, in which Tris was replaced by Bis-Tris for pH 6.0 to 7.0. An FP-8300 spectrofluorometer (Jasco) was used at 25°C to excite tS at 400 6 5 nm with emission spectra recorded between 450 and 600 6 5 nm and to excite mC at 580 6 5 nm with emission spectra recorded between 600 and 650 6 5 nm (Supplemental Figure 3). tS/mC ratios were calculated from the fluorescence maxima empirically determined at 513 and 604 nm, respectively, after NAD(P)H autofluorescence correction (Supplemental Figure 3C). Sensor proteins were allowed to equilibrate with 500 mM NAD 1 (or 150 mM NADP 1 ) for 30 min prior to measurements to ensure that no residual NADH was bound to the sensor. To compare the ratiometric responses of sensor proteins to NADH versus NADPH in the presence of 500 mM NAD 1 and 150 mM NADP 1 , respectively, ratio values were normalized to the highest ratio (set to 1). K d -values were determined by sigmoidal curve fitting using Prism software (version 8.0.1; GraphPad).

Isolation of Mitochondria
Mitochondria were isolated from the 16-d-old wild-type Arabidopsis seedlings grown in hydroponic pots under long-day conditions as described previously by Sweetlove et al. (2007) and Schwarzländer et al. (2011).

Multiwell Plate Reader-Based Fluorimetry
tS and mC were excited at 400 6 5 and 540 6 10 nm, respectively, and emission was collected at 520 6 5 and 615 6 9 nm in a CLARIOstar microplate reader (BMG Labtech). For in vivo experiments with Arabidopsis leaf discs or seedlings, samples were monitored in a total volume of 200 mL of assay medium (10 mM MES, 5 mM KCl, 10 mM MgCl 2 , and 10 mM CaCl 2 , pH 5.8) using top optics with adjusted focal height to 6 to 8 mm, 25°C incubation temperature, well-multichromatic monitoring, and 40 flashes per well and measuring cycle with orbital averaging to minimize local fluorescence artifacts. Sample sizes represent total number of leaf discs, which were cut from three to eight different plants per line cultured side by side under the same conditions. A maximum of 12 leaf discs per plant were cut and equally distributed to the treatment and mock control groups. Experiments were independently repeated with three different sensor lines and different plant batches of the same sensor line to validate the reproducibility of the results and to justify the presentation of representative data sets. Plates were kept in the dark for at least 90 min before recording to minimize potential effects of active photosynthesis. For the dark-light-dark transition experiments, measurements were paused, the plates were moved to be illuminated from the top (400 mmol s 21 m 22 ) for 5, 10, or 30 min, respectively, after which measurements were continued. For in vitro monitoring of mitochondrial NADH oxidation, assays were adapted from De Col et al. (2017). Proteins (10 mg) from freshly isolated mitochondria were supplemented with 500 mM NAD 1 and 0.5 mM purified recombinant Peredox-mCherry in a total volume of 200 mL per well. Samples were imaged using top optics with a focal height of 8.0 mm, 50 flashes per well and measuring cycle (orbital averaging), and double orbital shaking at 400 rpm for 10 s before each measurement cycle (excitation at 400 6 10 and 570 6 10 nm, emission collected at 515 6 7.5 and 610 6 5 nm, respectively, for tS and mC).

Confocal Laser Scanning Microscopy
Confocal imaging of 5-d-old Arabidopsis seedlings was performed using a DMi8 inverted microscope equipped with a TCS SP8 laser scanning device (Leica Microsystems) and a 310 (0.40 numerical aperture) or 340 (1.10 numerical aperture, water immersion) lens as previously described by Wagner et al. (2015). tS fluorescence excited at 405 nm was recorded at 517.5 6 7.5 nm and mC fluorescence excited at 575 nm (White Light Laser at 85%) was recorded at 612.5 6 7.5 nm, with sequential excitation and an open pinhole (5 airy units for seedling maps, 310 lens; 3 airy units for subcellular localization, 340 lens). Chlorophyll fluorescence was collected at 650 to 700 nm after excitation at 575 nm. To analyze cytosolic and nuclear sensor ratios in regions of interest (ROIs), images were recorded using an LSM980 inverted microscope (Carl Zeiss Microscopy) equipped with a 340 lens (C-Apochromat, 1.20 numerical aperture, water immersion). tS fluorescence excited at 405 nm was recorded at 517 6 9 nm and mC fluorescence excited at 561 nm was recorded at 615 6 9 nm, with sequential excitation and the pinhole set to 2 airy units. Plants were dark adapted for at least 90 min before image acquisition. Single plane images were processed with custom MATLAB-based software (MathWorks; Fricker, 2016) using x,y noise filtering and fluorescence background subtraction.

Statistical Analysis
Statistical analysis and curve fitting were performed as described in the figure legends using the software Prism (version 8.0.1; GraphPad; Supplemental Data Set 1). To investigate statistically significant phenotypic differences between genotypes, and tissue-specific differences in the NAD redox state between different growth conditions, two-way ANOVA followed by the Tukey test was performed. Differences in sensor-specific K d -values in response to NADH and NADPH were determined by sigmoidal curve fitting. Kinetics were analyzed by nonlinear or linear regression analysis as specified in Supplemental Data Set 1.

Supplemental Data
Supplemental Figure 1. Phenotypic characterization of plants of three independent Arabidopsis biosensor lines expressing Peredox-mCherry in the cytosol compared with the wild-type control.
Supplemental Figure 2. Optimization of protein overexpression and purification of Peredox-mCherry.
Supplemental Figure 5. In vivo fluorescence spectra of three independent Peredox-mCherry lines.
Supplemental Figure 6. In vivo live interventions of cytosolic NAD redox dynamics in seedlings and leaf discs of independent sensor lines. Figure 7. In vitro responses of Peredox-mCherry DS to NADH and NADH/NAD 1 .
Supplemental Figure 9. In vivo NAD redox dynamics of three independent Peredox-mCherry DS lines compared with Peredox-mCherry in response to changes in respiratory activity and external sugar supply.
Supplemental Table. Oligonucleotides used for cloning and mutagenesis.
Supplemental Data Set 1. Statistical analyses.

ACKNOWLEDGMENTS
We thank the Imaging Network of the University of Münster (RI_00497) and Jörg Kudla for access to the Leica SP8 confocal microscope system; Sonja Raetz, Martin Bonin, and Eva Regel for help with protein purification; and Thomas Nietzel and Richard Jaeger for preliminary optimizations of the flg22 experiments. We further thank the reviewers for adding important thoughts to this article and Alberto P. Macho for stimulating exchange on the role of NAD in immunity signaling. This work was supported by the