Multiscale Metabolic Modeling: Dynamic Flux Balance Analysis on a Whole-Plant Scale

Plant metabolism is characterized by a unique complexity on the cellular, tissue, and organ levels. On a whole-plant scale, changing source and sink relations accompanying plant development add another level of complexity to metabolism. With the aim of achieving a spatiotemporal resolution of source-sink interactions in crop plant metabolism, a multiscale metabolic modeling (MMM) approach was applied that integrates static organ-specific models with a whole-plant dynamic model. Allowing for a dynamic flux balance analysis on a whole-plant scale, the MMM approach was used to decipher the metabolic behavior of source and sink organs during the generative phase of the barley (Hordeum vulgare) plant. It reveals a sink-to-source shift of the barley stem caused by the senescence-related decrease in leaf source capacity, which is not sufficient to meet the nutrient requirements of sink organs such as the growing seed. The MMM platform represents a novel approach for the in silico analysis of metabolism on a whole-plant level, allowing for a systemic, spatiotemporally resolved understanding of metabolic processes involved in carbon partitioning, thus providing a novel tool for studying yield stability and crop improvement.Plants are of vital significance as a source of food (Grusak and DellaPenna, 1999; Rogalski and Carrer, 2011), feed (Lu et al., 2011), energy (Tilman et al., 2006; Parmar et al., 2011), and feedstocks for the chemical industry (Metzger and Bornscheuer, 2006; Kinghorn et al., 2011). Given the close connection between plant metabolism and the usability of plant products, there is a growing interest in understanding and predicting the behavior and regulation of plant metabolic processes. In order to increase crop quality and yield, there is a need for methods guiding the rational redesign of the plant metabolic network (Schwender, 2009).Mathematical modeling of plant metabolism offers new approaches to understand, predict, and modify complex plant metabolic processes. In plant research, the issue of metabolic modeling is constantly gaining attention, and different modeling approaches applied to plant metabolism exist, ranging from highly detailed quantitative to less complex qualitative approaches (for review, see Giersch, 2000; Morgan and Rhodes, 2002; Poolman et al., 2004; Rios-Estepa and Lange, 2007).A widely used modeling approach is flux balance analysis (FBA), which allows the prediction of metabolic capabilities and steady-state fluxes under different environmental and genetic backgrounds using (non)linear optimization (Orth et al., 2010). Assuming steady-state conditions, FBA has the advantage of not requiring the knowledge of kinetic parameters and, therefore, can be applied to model detailed, large-scale systems. In recent years, the FBA approach has been applied to several different plant species, such as maize (Zea mays; Dal’Molin et al., 2010; Saha et al., 2011), barley (Hordeum vulgare; Grafahrend-Belau et al., 2009b; Melkus et al., 2011; Rolletschek et al., 2011), rice (Oryza sativa; Lakshmanan et al., 2013), Arabidopsis (Arabidopsis thaliana; Poolman et al., 2009; de Oliveira Dal’Molin et al., 2010; Radrich et al., 2010; Williams et al., 2010; Mintz-Oron et al., 2012; Cheung et al., 2013), and rapeseed (Brassica napus; Hay and Schwender, 2011a, 2011b; Pilalis et al., 2011), as well as algae (Boyle and Morgan, 2009; Cogne et al., 2011; Dal’Molin et al., 2011) and photoautotrophic bacteria (Knoop et al., 2010; Montagud et al., 2010; Boyle and Morgan, 2011). These models have been used to study different aspects of metabolism, including the prediction of optimal metabolic yields and energy efficiencies (Dal’Molin et al., 2010; Boyle and Morgan, 2011), changes in flux under different environmental and genetic backgrounds (Grafahrend-Belau et al., 2009b; Dal’Molin et al., 2010; Melkus et al., 2011), and nonintuitive metabolic pathways that merit subsequent experimental investigations (Poolman et al., 2009; Knoop et al., 2010; Rolletschek et al., 2011). Although FBA of plant metabolic models was shown to be capable of reproducing experimentally determined flux distributions (Williams et al., 2010; Hay and Schwender, 2011b) and generating new insights into metabolic behavior, capacities, and efficiencies (Sweetlove and Ratcliffe, 2011), challenges remain to advance the utility and predictive power of the models.Given that many plant metabolic functions are based on interactions between different subcellular compartments, cell types, tissues, and organs, the reconstruction of organ-specific models and the integration of these models into interacting multiorgan and/or whole-plant models is a prerequisite to get insight into complex plant metabolic processes organized on a whole-plant scale (e.g. source-sink interactions). Almost all FBA models of plant metabolism are restricted to one cell type (Boyle and Morgan, 2009; Knoop et al., 2010; Montagud et al., 2010; Cogne et al., 2011; Dal’Molin et al., 2011), one tissue or one organ (Grafahrend-Belau et al., 2009b; Hay and Schwender, 2011a, 2011b; Pilalis et al., 2011; Mintz-Oron et al., 2012), and only one model exists taking into account the interaction between two cell types by specifying the interaction between mesophyll and bundle sheath cells in C4 photosynthesis (Dal’Molin et al., 2010). So far, no model representing metabolism at the whole-plant scale exists.Considering whole-plant metabolism raises the problem of taking into account temporal and environmental changes in metabolism during plant development and growth. Although classical static FBA is unable to predict the dynamics of metabolic processes, as the network analysis is based on steady-state solutions, time-dependent processes can be taken into account by extending the classical static FBA to a dynamic flux balance analysis (dFBA), as proposed by Mahadevan et al. (2002). The static (SOA) and dynamic optimization approaches introduced in this work provide a framework for analyzing the transience of metabolism by integrating kinetic expressions to dynamically constrain exchange fluxes. Due to the requirement of knowing or estimating a large number of kinetic parameters, so far dFBA has only been applied to a plant metabolic model once, to study the photosynthetic metabolism in the chloroplasts of C3 plants by a simplified model of five biochemical reactions (Luo et al., 2009). Integrating a dynamic model into a static FBA model is an alternative approach to perform dFBA.In this study, a multiscale metabolic modeling (MMM) approach was applied with the aim of achieving a spatiotemporal resolution of cereal crop plant metabolism. To provide a framework for the in silico analysis of the metabolic dynamics of barley on a whole-plant scale, the MMM approach integrates a static multiorgan FBA model and a dynamic whole-plant multiscale functional plant model (FPM) to perform dFBA. The performance of the novel whole-plant MMM approach was tested by studying source-sink interactions during the seed developmental phase of barley plants.     

Disruption of Fumarylacetoacetate Hydrolase Causes Spontaneous Cell Death under Short-Day Conditions in Arabidopsis

Fumarylacetoacetate hydrolase (FAH) hydrolyzes fumarylacetoacetate to fumarate and acetoacetate, the final step in the tyrosine (Tyr) degradation pathway that is essential to animals. Deficiency of FAH in animals results in an inborn lethal disorder. However, the role for the Tyr degradation pathway in plants remains to be elucidated. In this study, we isolated an Arabidopsis (Arabidopsis thaliana) short-day sensitive cell death1 (sscd1) mutant that displays a spontaneous cell death phenotype under short-day conditions. The SSCD1 gene was cloned via a map-based cloning approach and found to encode an Arabidopsis putative FAH. The spontaneous cell death phenotype of the sscd1 mutant was completely eliminated by further knockout of the gene encoding the putative homogentisate dioxygenase, which catalyzes homogentisate into maleylacetoacetate (the antepenultimate step) in the Tyr degradation pathway. Furthermore, treatment of Arabidopsis wild-type seedlings with succinylacetone, an abnormal metabolite caused by loss of FAH in the Tyr degradation pathway, mimicked the sscd1 cell death phenotype. These results demonstrate that disruption of FAH leads to cell death in Arabidopsis and suggest that the Tyr degradation pathway is essential for plant survival under short-day conditions.Programmed cell death (PCD) has been defined as a sequence of genetically regulated events that lead to the elimination of specific cells, tissues, or whole organs (Lockshin and Zakeri, 2004). In plants, PCD is essential for developmental processes and defense responses (Dangl et al., 1996; Greenberg, 1996; Durrant et al., 2007). One well-characterized example of plant PCD is the hypersensitive response occurring during incompatible plant-pathogen interactions (Lam, 2004), which results in cell death to form visible lesions at the site of infection by an avirulent pathogen and consequently limits the pathogen spread (Morel and Dangl, 1997).To date, a large number of mutants that display spontaneous cell death lesions have been identified in barley (Hordeum vulgare), maize (Zea mays), rice (Oryza sativa), and Arabidopsis (Arabidopsis thaliana; Marchetti et al., 1983; Wolter et al., 1993; Dietrich et al., 1994; Gray et al., 1997). Because lesions form in the absence of pathogen infection, these mutants have been collectively termed as lesion-mimic mutants. Many genes with regulatory roles in PCD and defense responses, including LESION SIMULATING DISEASE1, ACCELERATED CELL DEATH11, and VASCULAR ASSOCIATED DEATH1, have been cloned and characterized (Dietrich et al., 1997; Brodersen et al., 2002; Lorrain et al., 2004).The appearance of spontaneous cell death lesions in some lesion-mimic mutants is dependent on photoperiod. For example, the Arabidopsis mutant lesion simulating disease1 and myoinositol-1-phosphate synthase1 show lesions under long days (LD; Dietrich et al., 1994; Meng et al., 2009), whereas the lesion simulating disease2, lesion initiation1, enhancing RPW8-mediated HR-like cell death1, and lag one homolog1 display lesions under short days (SD; Dietrich et al., 1994; Ishikawa et al., 2003; Wang et al., 2008; Ternes et al., 2011).Blockage of some metabolic pathways in plants may cause cell death and result in lesion formation. For example, the lesion-mimic phenotypes in the Arabidopsis mutants lesion initiation2 and accelerated cell death2 and the maize mutant lesion mimic22 result from an impairment of porphyrin metabolism (Hu et al., 1998; Ishikawa et al., 2001; Mach et al., 2001). Deficiency in fatty acid, sphingolipid, and myoinositol metabolism also causes cell death in Arabidopsis (Mou et al., 2000; Liang et al., 2003; Wang et al., 2008; Meng et al., 2009; Donahue et al., 2010; Berkey et al., 2012).Tyr degradation is an essential five-step pathway in animals (Lindblad et al., 1977). First, Tyr aminotransferase catalyzes the conversion of Tyr into 4-hydroxyphenylpyruvate, which is further transformed into homogentisate by 4-hydroxyphenylpyruvate dioxygenase. Through the sequential action of homogentisate dioxygenase (HGO), maleylacetoacetate isomerase (MAAI), and fumarylacetoacetate hydrolase (FAH), homogentisate is catalyzed to generate fumarate and acetoacetate (Lindblad et al., 1977). Blockage of this pathway in animals results in metabolic disorder diseases (Lindblad et al., 1977; Ruppert et al., 1992; Grompe et al., 1993). For example, human FAH deficiency causes hereditary tyrosinemia type I (HT1), an inborn lethal disease (St-Louis and Tanguay, 1997). Although the homologous genes putatively encoding these enzymes exist in plants (Dixon et al., 2000; Lopukhina et al., 2001; Dixon and Edwards, 2006), it is unclear whether this pathway is essential for plant growth and development.In this study, we report the isolation and characterization of a recessive short-day sensitive cell death1 (sscd1) mutant in Arabidopsis. Map-based cloning of the corresponding gene revealed that SSCD1 encodes the Arabidopsis putative FAH. Further knockout of the gene encoding the Arabidopsis putative HGO completely eliminated the spontaneous cell death phenotype in the sscd1 mutant. Furthermore, we found that treatment of Arabidopsis wild-type seedlings with succinylacetone, an abnormal metabolite caused by loss of FAH in the Tyr degradation pathway (Lindblad et al., 1977), is able to mimic the sscd1 cell death phenotype. These results demonstrate that disruption of FAH leads to cell death in Arabidopsis and suggest that the Tyr degradation pathway is essential for plant survival under SD.     

Evolution of Conifer Diterpene Synthases: Diterpene Resin Acid Biosynthesis in Lodgepole Pine and Jack Pine Involves Monofunctional and Bifunctional Diterpene Synthases

Abscisic acid (ABA) induces stomatal closure and inhibits light-induced stomatal opening. The mechanisms in these two processes are not necessarily the same. It has been postulated that the ABA receptors involved in opening inhibition are different from those involved in closure induction. Here, we provide evidence that four recently identified ABA receptors (PYRABACTIN RESISTANCE1 [PYR1], PYRABACTIN RESISTANCE-LIKE1 [PYL1], PYL2, and PYL4) are not sufficient for opening inhibition in Arabidopsis (Arabidopsis thaliana). ABA-induced stomatal closure was impaired in the pyr1/pyl1/pyl2/pyl4 quadruple ABA receptor mutant. ABA inhibition of the opening of the mutant’s stomata remained intact. ABA did not induce either the production of reactive oxygen species and nitric oxide or the alkalization of the cytosol in the quadruple mutant, in accordance with the closure phenotype. Whole cell patch-clamp analysis of inward-rectifying K⁺ current in guard cells showed a partial inhibition by ABA, indicating that the ABA sensitivity of the mutant was not fully impaired. ABA substantially inhibited blue light-induced phosphorylation of H⁺-ATPase in guard cells in both the mutant and the wild type. On the other hand, in a knockout mutant of the SNF1-related protein kinase, srk2e, stomatal opening and closure, reactive oxygen species and nitric oxide production, cytosolic alkalization, inward-rectifying K⁺ current inactivation, and H⁺-ATPase phosphorylation were not sensitive to ABA.The phytohormone abscisic acid (ABA), which is synthesized in response to abiotic stresses, plays a key role in the drought hardiness of plants. Reducing transpirational water loss through stomatal pores is a major ABA response (Schroeder et al., 2001). ABA promotes the closure of open stomata and inhibits the opening of closed stomata. These effects are not simply the reverse of one another (Allen et al., 1999; Wang et al., 2001; Mishra et al., 2006).A class of receptors of ABA was identified (Ma et al., 2009; Park et al., 2009; Santiago et al., 2009; Nishimura et al., 2010). The sensitivity of stomata to ABA was strongly decreased in quadruple and sextuple mutants of the ABA receptor genes PYRABACTIN RESISTANCE/PYRABACTIN RESISTANCE-LIKE/REGULATORY COMPONENT OF ABSCISIC ACID RECEPTOR (PYR/PYL/RCAR; Nishimura et al., 2010; Gonzalez-Guzman et al., 2012). The PYR/PYL/RCAR receptors are involved in the early ABA signaling events, in which a sequence of interactions of the receptors with PROTEIN PHOSPHATASE 2Cs (PP2Cs) and subfamily 2 SNF1-RELATED PROTEIN KINASES (SnRK2s) leads to the activation of downstream ABA signaling targets in guard cells (Cutler et al., 2010; Kim et al., 2010; Weiner et al., 2010). Studies of Commelina communis and Vicia faba suggested that the ABA receptors involved in stomatal opening are not the same as the ABA receptors involved in stomatal closure (Allan et al., 1994; Anderson et al., 1994; Assmann, 1994; Schwartz et al., 1994). The roles of PYR/PYL/RCAR in either stomatal opening or closure remained to be elucidated.Blue light induces stomatal opening through the activation of plasma membrane H⁺-ATPase in guard cells that generates an inside-negative electrochemical gradient across the plasma membrane and drives K⁺ uptake through voltage-dependent inward-rectifying K⁺ channels (Assmann et al., 1985; Shimazaki et al., 1986; Blatt, 1987; Schroeder et al., 1987; Thiel et al., 1992). Phosphorylation of the penultimate Thr of the plasma membrane H⁺-ATPase is a prerequisite for blue light-induced activation of the H⁺-ATPase (Kinoshita and Shimazaki, 1999, 2002). ABA inhibits H⁺-ATPase activity through dephosphorylation of the penultimate Thr in the C terminus of the H⁺-ATPase in guard cells, resulting in prevention of the opening (Goh et al., 1996; Zhang et al., 2004; Hayashi et al., 2011). Inward-rectifying K⁺ currents (I_Kin) of guard cells are negatively regulated by ABA in addition to through the decline of the H⁺ pump-driven membrane potential difference (Schroeder and Hagiwara, 1989; Blatt, 1990; McAinsh et al., 1990; Schwartz et al., 1994; Grabov and Blatt, 1999; Saito et al., 2008). This down-regulation of ion transporters by ABA is essential for the inhibition of stomatal opening.A series of second messengers has been shown to mediate ABA-induced stomatal closure. Reactive oxygen species (ROS) produced by NADPH oxidases play a crucial role in ABA signaling in guard cells (Pei et al., 2000; Zhang et al., 2001; Kwak et al., 2003; Sirichandra et al., 2009; Jannat et al., 2011). Nitric oxide (NO) is an essential signaling component in ABA-induced stomatal closure (Desikan et al., 2002; Guo et al., 2003; Garcia-Mata and Lamattina, 2007; Neill et al., 2008). Alkalization of cytosolic pH in guard cells is postulated to mediate ABA-induced stomatal closure in Arabidopsis (Arabidopsis thaliana) and Pisum sativum and Paphiopedilum species (Irving et al., 1992; Gehring et al., 1997; Grabov and Blatt, 1997; Suhita et al., 2004; Gonugunta et al., 2008). These second messengers transduce environmental signals to ion channels and ion transporters that create the driving force for stomatal movements (Ward et al., 1995; MacRobbie, 1998; Garcia-Mata et al., 2003).In this study, we examined the mobilization of second messengers, the inactivation of I_Kin, and the suppression of H⁺-ATPase phosphorylation evoked by ABA in Arabidopsis mutants to clarify the downstream signaling events of ABA signaling in guard cells. The mutants included a quadruple mutant of PYR/PYL/RCARs, pyr1/pyl1/pyl2/pyl4, and a mutant of a SnRK2 kinase, srk2e.     

The Arabidopsis Endosomal Sorting Complex Required for Transport III Regulates Internal Vesicle Formation of the Prevacuolar Compartment and Is Required for Plant Development

We have established an efficient transient expression system with several vacuolar reporters to study the roles of endosomal sorting complex required for transport (ESCRT)-III subunits in regulating the formation of intraluminal vesicles of prevacuolar compartments (PVCs)/multivesicular bodies (MVBs) in plant cells. By measuring the distributions of reporters on/within the membrane of PVC/MVB or tonoplast, we have identified dominant negative mutants of ESCRT-III subunits that affect membrane protein degradation from both secretory and endocytic pathways. In addition, induced expression of these mutants resulted in reduction in luminal vesicles of PVC/MVB, along with increased detection of membrane-attaching vesicles inside the PVC/MVB. Transgenic Arabidopsis (Arabidopsis thaliana) plants with induced expression of ESCRT-III dominant negative mutants also displayed severe cotyledon developmental defects with reduced cell size, loss of the central vacuole, and abnormal chloroplast development in mesophyll cells, pointing out an essential role of the ESCRT-III complex in postembryonic development in plants. Finally, membrane dissociation of ESCRT-III components is important for their biological functions and is regulated by direct interaction among Vacuolar Protein Sorting-Associated Protein20-1 (VPS20.1), Sucrose Nonfermenting7-1, VPS2.1, and the adenosine triphosphatase VPS4/SUPPRESSOR OF K⁺ TRANSPORT GROWTH DEFECT1.Endomembrane trafficking in plant cells is complicated such that secretory, endocytic, and recycling pathways are usually integrated with each other at the post-Golgi compartments, among which, the trans-Golgi network (TGN) and prevacuolar compartment (PVC)/multivesicular body (MVB) are best studied (Tse et al., 2004; Lam et al., 2007a, 2007b; Müller et al., 2007; Foresti and Denecke, 2008; Hwang, 2008; Otegui and Spitzer, 2008; Robinson et al., 2008; Richter et al., 2009; Ding et al., 2012; Gao et al., 2014). Following the endocytic trafficking of a lipophilic dye, FM4-64, the TGN and PVC/MVB are sequentially labeled and thus are defined as the early and late endosome, respectively, in plant cells (Lam et al., 2007a; Chow et al., 2008). While the TGN is a tubular vesicular-like structure that may include several different microdomains and fit its biological function as a sorting station (Chow et al., 2008; Kang et al., 2011), the PVC/MVB is 200 to 500 nm in size with multiple luminal vesicles of approximately 40 nm (Tse et al., 2004). Membrane cargoes destined for degradation are sequestered into these tiny luminal vesicles and delivered to the lumen of the lytic vacuole (LV) via direct fusion between the PVC/MVB and the LV (Spitzer et al., 2009; Viotti et al., 2010; Cai et al., 2012). Therefore, the PVC/MVB functions between the TGN and LV as an intermediate organelle and decides the fate of membrane cargoes in the LV.In yeast (Saccharomyces cerevisiae), carboxypeptidase S (CPS) is synthesized as a type II integral membrane protein and sorted from the Golgi to the lumen of the vacuole (Spormann et al., 1992). Genetic analyses on the trafficking of CPS have led to the identification of approximately 17 class E genes (Piper et al., 1995; Babst et al., 1997, 2002a, 2002b; Odorizzi et al., 1998; Katzmann et al., 2001) that constitute the core endosomal sorting complex required for transport (ESCRT) machinery. The evolutionarily conserved ESCRT complex consists of several functionally different subcomplexes, ESCRT-0, ESCRT-I, ESCRT-II, and ESCRT-III and the ESCRT-III-associated/Vacuolar Protein Sorting4 (VPS4) complex. Together, they form a complex protein-protein interaction network that coordinates sorting of cargoes and inward budding of the membrane on the MVB (Hurley and Hanson, 2010; Henne et al., 2011). Cargo proteins carrying ubiquitin signals are thought to be passed from one ESCRT subcomplex to the next, starting with their recognition by ESCRT-0 (Bilodeau et al., 2002, 2003; Hislop and von Zastrow, 2011; Le Bras et al., 2011; Shields and Piper, 2011; Urbé, 2011). ESCRT-0 recruits the ESCRT-I complex, a heterotetramer of VPS23, VPS28, VPS37, and MVB12, from the cytosol to the endosomal membrane (Katzmann et al., 2001, 2003). The C terminus of VPS28 interacts with the N terminus of VPS36, a member of the ESCRT-II complex (Kostelansky et al., 2006; Teo et al., 2006). Then, cargoes passed from ESCRT-I and ESCRT-II are concentrated in certain membrane domains of the endosome by ESCRT-III, which includes four coiled-coil proteins and is sufficient to induce the membrane invagination (Babst et al., 2002b; Saksena et al., 2009; Wollert et al., 2009). Finally, the ESCRT components are disassociated from the membrane by the adenosine triphosphatase (ATPase) associated with diverse cellular activities (AAA) VPS4/SUPPRESSOR OF K+ TRANSPORT GROWTH DEFECT1 (SKD1) before releasing the internal vesicles (Babst et al., 1997, 1998).Putative homologs of ESCRT-I–ESCRT-III and ESCRT-III-associated components have been identified in plants, except for ESCRT-0, which is only present in Opisthokonta (Winter and Hauser, 2006; Leung et al., 2008; Schellmann and Pimpl, 2009). To date, only a few plant ESCRT components have been studied in detail. The Arabidopsis (Arabidopsis thaliana) AAA ATPase SKD1 localized to the PVC/MVB and showed ATPase activity that was regulated by Lysosomal Trafficking Regulator-Interacting Protein5, a plant homolog of Vps Twenty Associated1 Protein (Haas et al., 2007). Expression of the dominant negative form of SKD1 caused an increase in the size of the MVB and a reduction in the number of internal vesicles (Haas et al., 2007). This protein also contributes to the maintenance of the central vacuole and might be associated with cell cycle regulation, as leaf trichomes expressing its dominant negative mutant form lost the central vacuole and frequently contained multiple nuclei (Shahriari et al., 2010). Double null mutants of CHARGED MULTIVESICULAR BODY PROTEIN, chmp1achmp1b, displayed severe growth defects and were seedling lethal. This may be due to the mislocalization of plasma membrane (PM) proteins, including those involved in auxin transport such as PINFORMED1, PINFORMED2, and AUXIN-RESISTANT1, from the vacuolar degradation pathway to the tonoplast of the LV (Spitzer et al., 2009).Plant ESCRT components usually contain several homologs, with the possibility of functional redundancy. Single mutants of individual ESCRT components may not result in an obvious phenotype, whereas knockout of all homologs of an ESCRT component by generating double or triple mutants may be lethal to the plant. As a first step to carry out systematic analysis on each ESCRT complex in plant cells, here, we established an efficient analysis system to monitor the localization changes of four vacuolar reporters that accumulate either in the lumen (LRR84A-GFP, EMP12-GFP, and aleurain-GFP) or on the tonoplast (GFP-VIT1) of the LV and identified several ESCRT-III dominant negative mutants. We reported that ESCRT-III subunits were involved in the release of PVC/MVB’s internal vesicles from the limiting membrane and were required for membrane protein degradation from secretory and endocytic pathways. In addition, transgenic Arabidopsis plants with induced expression of ESCRT-III dominant negative mutants showed severe cotyledon developmental defects. We also showed that membrane dissociation of ESCRT-III subunits was regulated by direct interaction with SKD1.     

The Role of a Potassium Transporter OsHAK5 in Potassium Acquisition and Transport from Roots to Shoots in Rice at Low Potassium Supply Levels

In plants, K transporter (KT)/high affinity K transporter (HAK)/K uptake permease (KUP) is the largest potassium (K) transporter family; however, few of the members have had their physiological functions characterized in planta. Here, we studied OsHAK5 of the KT/HAK/KUP family in rice (Oryza sativa). We determined its cellular and tissue localization and analyzed its functions in rice using both OsHAK5 knockout mutants and overexpression lines in three genetic backgrounds. A β-glucuronidase reporter driven by the OsHAK5 native promoter indicated OsHAK5 expression in various tissue organs from root to seed, abundantly in root epidermis and stele, the vascular tissues, and mesophyll cells. Net K influx rate in roots and K transport from roots to aerial parts were severely impaired by OsHAK5 knockout but increased by OsHAK5 overexpression in 0.1 and 0.3 mm K external solution. The contribution of OsHAK5 to K mobilization within the rice plant was confirmed further by the change of K concentration in the xylem sap and K distribution in the transgenic lines when K was removed completely from the external solution. Overexpression of OsHAK5 increased the K-sodium concentration ratio in the shoots and salt stress tolerance (shoot growth), while knockout of OsHAK5 decreased the K-sodium concentration ratio in the shoots, resulting in sensitivity to salt stress. Taken together, these results demonstrate that OsHAK5 plays a major role in K acquisition by roots faced with low external K and in K upward transport from roots to shoots in K-deficient rice plants.Potassium (K) is one of the three most important macronutrients and the most abundant cation in plants. As a major osmoticum in the vacuole, K drives the generation of turgor pressure, enabling cell expansion. In the vascular tissue, K is an important participant in the generation of root pressure (for review, see Wegner, 2014 [including his new hypothesis]). In the phloem, K is critical for the transport of photoassimilates from source to sink (Marschner, 1996; Deeken et al., 2002; Gajdanowicz et al., 2011). In addition, enhancing K absorption and decreasing sodium (Na) accumulation is a major strategy of glycophytes in salt stress tolerance (Maathuis and Amtmann, 1999; Munns and Tester, 2008; Shabala and Cuin, 2008).Plants acquire K through K-permeable proteins at the root surface. Since available K concentration in the soil may vary by 100-fold, plants have developed multiple K uptake systems for adapting to this variability (Epstein et al., 1963; Grabov, 2007; Maathuis, 2009). In a classic K uptake experiment in barley (Hordeum vulgare), root K absorption has been described as a high-affinity and low-affinity biphasic transport process (Epstein et al., 1963). It is generally assumed that the low-affinity transport system (LATS) in the roots mediates K uptake in the millimolar range and that the activity of this system is insensitive to external K concentration (Maathuis and Sanders, 1997; Chérel et al., 2014). In contrast, the high-affinity transport system (HATS) was rapidly up-regulated when the supply of exogenous K was halted (Glass, 1976; Glass and Dunlop, 1978).The membrane transporters for K flux identified in plants are generally classified into three channels and three transporter families based on phylogenetic analysis (Mäser et al., 2001; Véry and Sentenac, 2003; Lebaudy et al., 2007; Alemán et al., 2011). For K uptake, it was predicted that, under most circumstances, K transporters function as HATS, while K-permeable channels mediate LATS (Maathuis and Sanders, 1997). However, a root-expressed K channel in Arabidopsis (Arabidopsis thaliana), Arabidopsis K Transporter1 (AKT1), mediates K absorption over a wide range of external K concentrations (Sentenac et al., 1992; Lagarde et al., 1996; Hirsch et al., 1998; Spalding et al., 1999), while evidence is accumulating that many K transporters, including members of the K transporter (KT)/high affinity K transporter (HAK)/K uptake permease (KUP) family, are low-affinity K transporters (Quintero and Blatt, 1997; Senn et al., 2001), implying that functions of plant K channels and transporters overlap at different K concentration ranges.Out of the three families of K transporters, cation proton antiporter (CPA), high affinity K/Na transporter (HKT), and KT/HAK/KUP, CPA was characterized as a K⁺(Na⁺)/H⁺ antiporter, HKT may cotransport Na and K or transport Na only (Rubio et al., 1995; Uozumi et al., 2000), while KT/HAK/KUP were predicted to be H⁺-coupled K⁺ symporters (Mäser et al., 2001; Lebaudy et al., 2007). KT/HAK/KUP were named by different researchers who first identified and cloned them (Quintero and Blatt, 1997; Santa-María et al., 1997). In plants, the KT/HAK/KUP family is the largest K transporter family, including 13 members in Arabidopsis and 27 members in the rice (Oryza sativa) genome (Rubio et al., 2000; Mäser et al., 2001; Bañuelos et al., 2002; Gupta et al., 2008). Sequence alignments show that genes of this family share relatively low homology to each other. The KT/HAK/KUP family was divided into four major clusters (Rubio et al., 2000; Gupta et al., 2008), and in cluster I and II, they were further separated into A and B groups. Genes of cluster I or II likely exist in all plants, cluster III is composed of genes from both Arabidopsis and rice, while cluster IV includes only four rice genes (Grabov, 2007; Gupta et al., 2008).The functions of KT/HAK/KUP were studied mostly in heterologous expression systems. Transporters of cluster I, such as AtHAK5, HvHAK1, OsHAK1, and OsHAK5, are localized in the plasma membrane (Kim et al., 1998; Bañuelos et al., 2002; Gierth et al., 2005) and exhibit high-affinity K uptake in the yeast Saccharomyces cerevisiae (Santa-María et al., 1997; Fu and Luan, 1998; Rubio et al., 2000) and in Escherichia coli (Horie et al., 2011). Transporters of cluster II, like AtKUP4 (TINY ROOT HAIRS1, TRH1), HvHAK2, OsHAK2, OsHAK7, and OsHAK10, could not complement the K uptake-deficient yeast (Saccharomyces cerevisiae) but were able to mediate K fluxes in a bacterial mutant; they might be tonoplast transporters (Senn et al., 2001; Bañuelos et al., 2002; Rodríguez-Navarro and Rubio, 2006). The function of transporters in clusters III and IV is even less known (Grabov, 2007).Existing data suggest that some KT/HAK/KUP transporters also may respond to salinity stress (Maathuis, 2009). The cluster I transporters of HvHAK1 mediate Na influx (Santa-María et al., 1997), while AtHAK5 expression is inhibited by Na (Rubio et al., 2000; Nieves-Cordones et al., 2010). Expression of OsHAK5 in tobacco (Nicotiana tabacum) BY2 cells enhanced the salt tolerance of these cells by accumulating more K without affecting their Na content (Horie et al., 2011).There are only scarce reports on the physiological function of KT/HAK/KUP in planta. In Arabidopsis, mutation of AtKUP2 (SHORT HYPOCOTYL3) resulted in a short hypocotyl, small leaves, and a short flowering stem (Elumalai et al., 2002), while a loss-of-function mutation of AtKUP4 (TRH1) resulted in short root hairs and a loss of gravity response in the root (Rigas et al., 2001; Desbrosses et al., 2003; Ahn et al., 2004). AtHAK5 is the only system currently known to mediate K uptake at concentrations below 0.01 mm (Rubio et al., 2010) and provides a cesium uptake pathway (Qi et al., 2008). AtHAK5 and AtAKT1 are the two major physiologically relevant molecular entities mediating K uptake into roots in the range between 0.01 and 0.05 mm (Pyo et al., 2010; Rubio et al., 2010). AtAKT1 may contribute to K uptake within the K concentrations that belong to the high-affinity system described by Epstein et al. (1963).Among all 27 members of the KT/HAK/KUP family in rice, OsHAK1, OsHAK5, OsHAK19, and OsHAK20 were grouped in cluster IB (Gupta et al., 2008). These four rice HAK members share 50.9% to 53.4% amino acid identity with AtHAK5. OsHAK1 was expressed in the whole plant, with maximum expression in roots, and was up-regulated by K deficiency; it mediated high-affinity K uptake in yeast (Bañuelos et al., 2002). In this study, we examined the tissue-specific localization and the physiological functions of OsHAK5 in response to variation in K supply and to salt stress in rice. By comparing K uptake and translocation in OsHAK5 knockout (KO) mutants and in OsHAK5-overexpressing lines with those in their respective wild-type lines supplied with different K concentrations, we found that OsHAK5 not only mediates high-affinity K acquisition but also participates in root-to-shoot K transport as well as in K-regulated salt tolerance.