Spore Density Determines Infection Strategy by the Plant Pathogenic Fungus Plectosphaerella cucumerina

Necrotrophic and biotrophic pathogens are resisted by different plant defenses. While necrotrophic pathogens are sensitive to jasmonic acid (JA)-dependent resistance, biotrophic pathogens are resisted by salicylic acid (SA)- and reactive oxygen species (ROS)-dependent resistance. Although many pathogens switch from biotrophy to necrotrophy during infection, little is known about the signals triggering this transition. This study is based on the observation that the early colonization pattern and symptom development by the ascomycete pathogen Plectosphaerella cucumerina (P. cucumerina) vary between inoculation methods. Using the Arabidopsis (Arabidopsis thaliana) defense response as a proxy for infection strategy, we examined whether P. cucumerina alternates between hemibiotrophic and necrotrophic lifestyles, depending on initial spore density and distribution on the leaf surface. Untargeted metabolome analysis revealed profound differences in metabolic defense signatures upon different inoculation methods. Quantification of JA and SA, marker gene expression, and cell death confirmed that infection from high spore densities activates JA-dependent defenses with excessive cell death, while infection from low spore densities induces SA-dependent defenses with lower levels of cell death. Phenotyping of Arabidopsis mutants in JA, SA, and ROS signaling confirmed that P. cucumerina is differentially resisted by JA- and SA/ROS-dependent defenses, depending on initial spore density and distribution on the leaf. Furthermore, in situ staining for early callose deposition at the infection sites revealed that necrotrophy by P. cucumerina is associated with elevated host defense. We conclude that P. cucumerina adapts to early-acting plant defenses by switching from a hemibiotrophic to a necrotrophic infection program, thereby gaining an advantage of immunity-related cell death in the host.Plant pathogens are often classified as necrotrophic or biotrophic, depending on their infection strategy (Glazebrook, 2005; Nishimura and Dangl, 2010). Necrotrophic pathogens kill living host cells and use the decayed plant tissue as a substrate to colonize the plant, whereas biotrophic pathogens parasitize living plant cells by employing effector molecules that suppress the host immune system (Pel and Pieterse, 2013). Despite this binary classification, the majority of pathogenic microbes employ a hemibiotrophic infection strategy, which is characterized by an initial biotrophic phase followed by a necrotrophic infection strategy at later stages of infection (Perfect and Green, 2001). The pathogenic fungi Magnaporthe grisea, Sclerotinia sclerotiorum, and Mycosphaerella graminicola, the oomycete Phytophthora infestans, and the bacterial pathogen Pseudomonas syringae are examples of hemibiotrophic plant pathogens (Perfect and Green, 2001; Koeck et al., 2011; van Kan et al., 2014; Kabbage et al., 2015).Despite considerable progress in our understanding of plant resistance to necrotrophic and biotrophic pathogens (Glazebrook, 2005; Mengiste, 2012; Lai and Mengiste, 2013), recent debate highlights the dynamic and complex interplay between plant-pathogenic microbes and their hosts, which is raising concerns about the use of infection strategies as a static tool to classify plant pathogens. For instance, the fungal genus Botrytis is often labeled as an archetypal necrotroph, even though there is evidence that it can behave as an endophytic fungus with a biotrophic lifestyle (van Kan et al., 2014). The rice blast fungus Magnaporthe oryzae, which is often classified as a hemibiotrophic leaf pathogen (Perfect and Green, 2001; Koeck et al., 2011), can adopt a purely biotrophic lifestyle when infecting root tissues (Marcel et al., 2010). It remains unclear which signals are responsible for the switch from biotrophy to necrotrophy and whether these signals rely solely on the physiological state of the pathogen, or whether host-derived signals play a role as well (Kabbage et al., 2015).The plant hormones salicylic acid (SA) and jasmonic acid (JA) play a central role in the activation of plant defenses (Glazebrook, 2005; Pieterse et al., 2009, 2012). The first evidence that biotrophic and necrotrophic pathogens are resisted by different immune responses came from Thomma et al. (1998), who demonstrated that Arabidopsis (Arabidopsis thaliana) genotypes impaired in SA signaling show enhanced susceptibility to the biotrophic pathogen Hyaloperonospora arabidopsidis (formerly known as Peronospora parastitica), while JA-insensitive genotypes were more susceptible to the necrotrophic fungus Alternaria brassicicola. In subsequent years, the differential effectiveness of SA- and JA-dependent defense mechanisms has been confirmed in different plant-pathogen interactions, while additional plant hormones, such as ethylene, abscisic acid (ABA), auxins, and cytokinins, have emerged as regulators of SA- and JA-dependent defenses (Bari and Jones, 2009; Cao et al., 2011; Pieterse et al., 2012). Moreover, SA- and JA-dependent defense pathways have been shown to act antagonistically on each other, which allows plants to prioritize an appropriate defense response to attack by biotrophic pathogens, necrotrophic pathogens, or herbivores (Koornneef and Pieterse, 2008; Pieterse et al., 2009; Verhage et al., 2010).In addition to plant hormones, reactive oxygen species (ROS) play an important regulatory role in plant defenses (Torres et al., 2006; Lehmann et al., 2015). Within minutes after the perception of pathogen-associated molecular patterns, NADPH oxidases and apoplastic peroxidases generate early ROS bursts (Torres et al., 2002; Daudi et al., 2012; O’Brien et al., 2012), which activate downstream defense signaling cascades (Apel and Hirt, 2004; Torres et al., 2006; Miller et al., 2009; Mittler et al., 2011; Lehmann et al., 2015). ROS play an important regulatory role in the deposition of callose (Luna et al., 2011; Pastor et al., 2013) and can also stimulate SA-dependent defenses (Chaouch et al., 2010; Yun and Chen, 2011; Wang et al., 2014; Mammarella et al., 2015). However, the spread of SA-induced apoptosis during hyperstimulation of the plant immune system is contained by the ROS-generating NADPH oxidase RBOHD (Torres et al., 2005), presumably to allow for the sufficient generation of SA-dependent defense signals from living cells that are adjacent to apoptotic cells. Nitric oxide (NO) plays an additional role in the regulation of SA/ROS-dependent defense (Trapet et al., 2015). This gaseous molecule can stimulate ROS production and cell death in the absence of SA while preventing excessive ROS production at high cellular SA levels via S-nitrosylation of RBOHD (Yun et al., 2011). Recently, it was shown that pathogen-induced accumulation of NO and ROS promotes the production of azelaic acid, a lipid derivative that primes distal plants for SA-dependent defenses (Wang et al., 2014). Hence, NO, ROS, and SA are intertwined in a complex regulatory network to mount local and systemic resistance against biotrophic pathogens. Interestingly, pathogens with a necrotrophic lifestyle can benefit from ROS/SA-dependent defenses and associated cell death (Govrin and Levine, 2000). For instance, Kabbage et al. (2013) demonstrated that S. sclerotiorum utilizes oxalic acid to repress oxidative defense signaling during initial biotrophic colonization, but it stimulates apoptosis at later stages to advance necrotrophic colonization. Moreover, SA-induced repression of JA-dependent resistance not only benefits necrotrophic pathogens but also hemibiotrophic pathogens after having switched from biotrophy to necrotrophy (Glazebrook, 2005; Pieterse et al., 2009, 2012).Plectosphaerella cucumerina ((P. cucumerina, anamorph Plectosporum tabacinum) anamorph Plectosporum tabacinum) is a filamentous ascomycete fungus that can survive saprophytically in soil by decomposing plant material (Palm et al., 1995). The fungus can cause sudden death and blight disease in a variety of crops (Chen et al., 1999; Harrington et al., 2000). Because P. cucumerina can infect Arabidopsis leaves, the P. cucumerina-Arabidopsis interaction has emerged as a popular model system in which to study plant defense reactions to necrotrophic fungi (Berrocal-Lobo et al., 2002; Ton and Mauch-Mani, 2004; Carlucci et al., 2012; Ramos et al., 2013). Various studies have shown that Arabidopsis deploys a wide range of inducible defense strategies against P. cucumerina, including JA-, SA-, ABA-, and auxin-dependent defenses, glucosinolates (Tierens et al., 2001; Sánchez-Vallet et al., 2010; Gamir et al., 2014; Pastor et al., 2014), callose deposition (García-Andrade et al., 2011; Gamir et al., 2012, 2014; Sánchez-Vallet et al., 2012), and ROS (Tierens et al., 2002; Sánchez-Vallet et al., 2010; Barna et al., 2012; Gamir et al., 2012, 2014; Pastor et al., 2014). Recent metabolomics studies have revealed large-scale metabolic changes in P. cucumerina-infected Arabidopsis, presumably to mobilize chemical defenses (Sánchez-Vallet et al., 2010; Gamir et al., 2014; Pastor et al., 2014). Furthermore, various chemical agents have been reported to induce resistance against P. cucumerina. These chemicals include β-amino-butyric acid, which primes callose deposition and SA-dependent defenses, benzothiadiazole (BTH or Bion; Görlach et al., 1996; Ton and Mauch-Mani, 2004), which activates SA-related defenses (Lawton et al., 1996; Ton and Mauch-Mani, 2004; Gamir et al., 2014; Luna et al., 2014), JA (Ton and Mauch-Mani, 2004), and ABA, which primes ROS and callose deposition (Ton and Mauch-Mani, 2004; Pastor et al., 2013). However, among all these studies, there is increasing controversy about the exact signaling pathways and defense responses contributing to plant resistance against P. cucumerina. While it is clear that JA and ethylene contribute to basal resistance against the fungus, the exact roles of SA, ABA, and ROS in P. cucumerina resistance vary between studies (Thomma et al., 1998; Ton and Mauch-Mani, 2004; Sánchez-Vallet et al., 2012; Gamir et al., 2014).This study is based on the observation that the disease phenotype during P. cucumerina infection differs according to the inoculation method used. We provide evidence that the fungus follows a hemibiotrophic infection strategy when infecting from relatively low spore densities on the leaf surface. By contrast, when challenged by localized host defense to relatively high spore densities, the fungus switches to a necrotrophic infection program. Our study has uncovered a novel strategy by which plant-pathogenic fungi can take advantage of the early immune response in the host plant.     

Control of Tobacco mosaic virus Movement Protein Fate by CELL-DIVISION-CYCLE Protein48

Like many other viruses, Tobacco mosaic virus replicates in association with the endoplasmic reticulum (ER) and exploits this membrane network for intercellular spread through plasmodesmata (PD), a process depending on virus-encoded movement protein (MP). The movement process involves interactions of MP with the ER and the cytoskeleton as well as its targeting to PD. Later in the infection cycle, the MP further accumulates and localizes to ER-associated inclusions, the viral factories, and along microtubules before it is finally degraded. Although these patterns of MP accumulation have been described in great detail, the underlying mechanisms that control MP fate and function during infection are not known. Here, we identify CELL-DIVISION-CYCLE protein48 (CDC48), a conserved chaperone controlling protein fate in yeast (Saccharomyces cerevisiae) and animal cells by extracting protein substrates from membranes or complexes, as a cellular factor regulating MP accumulation patterns in plant cells. We demonstrate that Arabidopsis (Arabidopsis thaliana) CDC48 is induced upon infection, interacts with MP in ER inclusions dependent on the MP N terminus, and promotes degradation of the protein. We further provide evidence that CDC48 extracts MP from ER inclusions to the cytosol, where it subsequently accumulates on and stabilizes microtubules. We show that virus movement is impaired upon overexpression of CDC48, suggesting that CDC48 further functions in controlling virus movement by removal of MP from the ER transport pathway and by promoting interference of MP with microtubule dynamics. CDC48 acts also in response to other proteins expressed in the ER, thus suggesting a general role of CDC48 in ER membrane maintenance upon ER stress.Plant viruses are obligate intracellular pathogens that replicate in association with host membranes (Laliberté and Sanfaçon, 2010) and subvert host intra- and intercellular trafficking pathways to achieve cell-to-cell and systemic spread (Harries and Ding, 2011; Niehl and Heinlein, 2011). In the case of the well-studied Tobacco mosaic virus (TMV), viral replication factories form on membranes of the endoplasmic reticulum (ER; Heinlein et al., 1995, 1998). As the plant ER is continuous between cells through plasmodesmata (PD; Ding et al., 1992), this membrane network provides a direct pathway for the spread of replicated virus from the replication sites in infected cells into the ER network of noninfected cells. The spread of plant viruses depends on virus-encoded movement proteins (MPs; Deom et al., 1987; Lucas, 2006). The MP of TMV facilitates the cell-to-cell passage of the infectious particle by forming a ribonucleoprotein complex with the viral RNA (Citovsky et al., 1990) and by increasing the size exclusion limit of PD (Wolf et al., 1989).During the course of infection, as well as when ectopically expressed, the MP associates with PD, the ER/actin network, and microtubules (Heinlein et al., 1995, 1998; Reichel and Beachy, 1998; Wright et al., 2007; Sambade et al., 2008; Hofmann et al., 2009; Boutant et al., 2010; Peña and Heinlein, 2012; Supplemental Fig. S1). Shortly after infection of a new cell, the MP localizes to small, mobile, ER-associated particles proposed to play a role in PD targeting of the viral RNA (Boyko et al., 2007; Sambade et al., 2008). Similar small, mobile MP particles are observed early upon ectopic expression of the protein. These particles colocalize with RNA and undergo stop-and-go movements in association with the ER (Sambade et al., 2008). The particle movements pause at microtubule proximal sites and their detachment requires microtubule polymerization (Sambade et al., 2008). These observations suggest that the interaction with the microtubule system plays a critical role in the maturation and ER-mediated delivery of infectious viral RNA particles to PD during early infection stages. Consistently, tobacco (Nicotiana tabacum) mutants with reduced microtubule dynamics exhibit reduced TMV movement (Ouko et al., 2010). Following virus movement, the previously infected cell further accumulates MP at the ER, a process that coincides with the formation of large ER inclusions that contain viral replicase and viral RNA in addition to MP and likely function as virus factories (Heinlein et al., 1998; Más and Beachy, 1999). In mature form, these inclusions may represent the so-called viroplasms or X-bodies described in the classical literature (Bawden and Sheffield, 1939; Esau and Cronshaw, 1967; Hills et al., 1987). Their formation is associated with rearrangements of the ER membrane and likely mediated by the accumulated MP since the inclusions diminish and reconstitute a native ER structure when MP becomes degraded by the 26S proteasome (Reichel and Beachy, 1998, 2000). Transfected cells accumulate MP in similar inclusions as those formed during infection, indicating that accumulated MP is indeed necessary and sufficient to form inclusions in association with the ER (Reichel and Beachy, 1998; Supplemental Fig. S1). Following accumulation of MP in virus factories, the infected cells accumulate the MP also along microtubules (Heinlein et al., 1998). The accumulation of MP in virus factories and on microtubules in cells behind the leading front of infection is dispensable for virus movement (Heinlein et al., 1998; Boyko et al., 2000a). At these late infection stages, the virus factories may enable the virus to produce high virion titers (Laliberté and Sanfaçon, 2010; Tilsner et al., 2012), and the subsequent accumulation along microtubules may play a role in withdrawing MP from the cell-to-cell communication pathway (Curin et al., 2007) and in stockpiling MP prior to degradation (Padgett et al., 1996; Gillespie et al., 2002).The molecular mechanisms that guide the MP to the ER and subsequently to microtubules during infection are not known. The MP is a hydrophobic protein that behaves like a membrane-integral or tightly membrane-associated protein in differential fractionation experiments and contains two predicted transmembrane domains (Reichel and Beachy, 1998; Brill et al., 2000, 2004) involved in ER association (Fujiki et al., 2006). The association with microtubules depends on MP amino acids 1 to 213 required for MP function (Kahn et al., 1998; Boyko et al., 2000b,Boyko et al., 2000c, 2002; Kotlizky et al., 2001). Moreover, certain amino acid exchange mutations known to affect the function of MP in virus movement in a temperature-sensitive manner also affect the ability of MP to interact with microtubules (Boyko et al., 2007,Boyko et al., 2000b). Interestingly, these mutations cluster together in a short domain of 25 amino acids showing a structural similarity with the M-loop of tubulin involved in tubulin-tubulin interactions (Boyko et al., 2000b; Waigmann et al., 2007). Importantly, this M-loop similarity domain overlaps with the predicted transmembrane domain (Brill et al., 2000, 2004) thus suggesting that the association of MP with membranes or microtubules is an alternative event that may depend on specific posttranslational modifications or specific folds of MP. However, although the different subcellular localizations of MPs during the course of infection indicate directional transport of MP from the ER to microtubules and may indicate different folds and functions of the protein when associated with these different subcellular components, the mechanism that controls the subcellular localization and, thus, the fate and function of MP is not known.Here, we identify CELL-DIVISION-CYCLE protein48 (CDC48), named p97/VCP (Valosin-containing protein) in mammals and Cdc48p in yeast (Saccharomyces cerevisiae), as a cellular factor regulating MP subcellular accumulation patterns. CDC48 functions are well characterized in mammalian and yeast systems but remain poorly investigated in plants. Yeast and mammalian CDC48s are essential, conserved chaperones involved in diverse cellular processes by controlling protein fate through extraction of substrates from membranes or complexes (Tsai et al., 2002; Meusser et al., 2005; Römisch, 2005; Rumpf and Jentsch, 2006; Schrader et al., 2009; Eisele et al., 2010; Meyer et al., 2012; Yamanaka et al., 2012). We show that virus infection leads to the induction of Arabidopsis (Arabidopsis thaliana) CDC48 isoforms and demonstrate a function of CDC48 in ER maintenance upon ER stress conditions. We further demonstrate that CDC48 interacts with MP and that CDC48 activity is required for MP degradation. Interaction of CDC48 with MP depends on the MP N terminus, which is required for degradation of the protein, for PD localization and microtubule accumulation of MP, and for function of MP in cell-to-cell transport of the viral RNA. Overexpressed CDC48 shifts MP subcellular localization from ER inclusions to microtubules, suggesting that CDC48 extracts the MP from ER-associated inclusions, where it accumulates in midstages of infection, to the cytosol, where it accumulates along microtubules during late infection stages. Moreover, overexpression of active, but not inactive, CDC48 inhibits virus movement. Our data demonstrate that a CDC48-dependent pathway leading to the clearance of ER-associated protein inclusions exists in plants, that plant viral MPs are substrates for this pathway, and that this pathway determines viral protein fate during infection. We suggest that CDC48-mediated extraction of MP from the ER is part of a plant defense response to remove MP from the ER, the compartment the virus uses for replication and movement.     

Cooperative Protein Folding by Two Protein Thiol Disulfide Oxidoreductases and ERO1 in Soybean

Most proteins produced in the endoplasmic reticulum (ER) of eukaryotic cells fold via disulfide formation (oxidative folding). Oxidative folding is catalyzed by protein disulfide isomerase (PDI) and PDI-related ER protein thiol disulfide oxidoreductases (ER oxidoreductases). In yeast and mammals, ER oxidoreductin-1s (Ero1s) supply oxidizing equivalent to the active centers of PDI. In this study, we expressed recombinant soybean Ero1 (GmERO1a) and found that GmERO1a oxidized multiple soybean ER oxidoreductases, in contrast to mammalian Ero1s having a high specificity for PDI. One of these ER oxidoreductases, GmPDIM, associated in vivo and in vitro with GmPDIL-2, was unable to be oxidized by GmERO1a. We therefore pursued the possible cooperative oxidative folding by GmPDIM, GmERO1a, and GmPDIL-2 in vitro and found that GmPDIL-2 synergistically accelerated oxidative refolding. In this process, GmERO1a preferentially oxidized the active center in the a′ domain among the a, a′, and b domains of GmPDIM. A disulfide bond introduced into the active center of the a′ domain of GmPDIM was shown to be transferred to the active center of the a domain of GmPDIM and the a domain of GmPDIM directly oxidized the active centers of both the a or a′ domain of GmPDIL-2. Therefore, we propose that the relay of an oxidizing equivalent from one ER oxidoreductase to another may play an essential role in cooperative oxidative folding by multiple ER oxidoreductases in plants.In eukaryotes, many secretory and membrane proteins fold via disulfide bond formation in the endoplasmic reticulum (ER). Seed storage proteins of major crops, such as wheat, corn, rice, and beans, which are important protein sources for humans and domestic animals, are synthesized in the ER of the endosperm or cotyledon. A number of seed storage proteins fold by the formation of intramolecular disulfide bonds (oxidative folding) and are transported to and accumulate in protein bodies (Kermode and Bewley, 1999; Jolliffe et al., 2005). In contrast to normally folded proteins, misfolded and unfolded proteins are retained in the ER and degraded by an ER-associated degradation or vacuolar system (Smith et al., 2011; Pu and Bassham, 2013). Therefore, quick and efficient oxidative folding of nascent seed storage proteins is needed for their accumulation in protein bodies.During this process, protein disulfide isomerase (PDI; EC 5.3.4.1) and other ER protein thiol disulfide oxidoreductases (ER oxidoreductases) are thought to catalyze the formation and isomerization of disulfide bonds in nascent proteins (Hatahet and Ruddock, 2009; Feige and Hendershot, 2011; Lu and Holmgren, 2014). After phylogenetic analysis of the Arabidopsis genome, 10 classes of ER oxidoreductases (classes I–X) were identified (Houston et al., 2005). Among them, class I ER oxidoreductase, a plant PDI ortholog, has been studied in a wide variety of plants. Class I ER oxidoreductases have two catalytically active domains a and a′, containing active centers composed of Cys-Gly-His-Cys and two catalytically inactive domains b and b′. An Arabidopsis ortholog of class I ER oxidoreductases is required for proper seed development and regulates the timing of programmed cell death by chaperoning and inhibiting Cys proteases (Andème Ondzighi et al., 2008). OaPDI, a PDI from Oldenlandia affinis, a coffee family (Rubiaceae) plant, is involved in the folding of knotted circular proteins (Gruber et al., 2007). The rice ortholog (PDIL1-1) was suggested to be involved in the maturation of the major seed storage protein glutelin (Takemoto et al., 2002). Furthermore, rice PDIL1-1 plays a role in regulatory activities for various proteins that are essential for the synthesis of grain components as determined by analysis of a T-DNA insertion mutant (Satoh-Cruz et al., 2010).The oxidative refolding ability of class I ER oxidoreductases was confirmed in recombinant soybean (GmPDIL-1) and wheat proteins produced by an Escherichia coli expression system established from cDNAs (Kamauchi et al., 2008; Kimura et al., 2015).Class II and III ER oxidoreductases have an a–b–b′–a′ domain structure. Class II ER oxidoreductases have an acidic amino acid-rich sequence in the N-terminal region ahead of the a domain. Recombinant soybean (GmPDIL-2) and wheat class II ER oxidoreductases have oxidative refolding activities similar to that of class I (Kamauchi et al., 2008; Kimura et al., 2015). Class III ER oxidoreductases contain the nonclassical redox-center Cys-X-X-Ser/Cys motifs, as opposed to the more traditional CGHC sequence, in the a and a′ domains. Recombinant soybean (GmPDIL-3) and wheat proteins lack oxidative refolding activity in vitro (Iwasaki et al., 2009; Kimura et al., 2015). Class IV ER oxidoreductases are unique to plants and have an a–a′–ERp29 domain structure, which is homologous to the C-terminal domain of mammalian ERp29 (Demmer et al., 1997).Recombinant soybean class IV ER oxidoreductases (GmPDIS-1 and GmPDIS-2) and wheat class IV ER oxidoreductase possess an oxidative refolding activity that is weaker than that of classes I and II (Wadahama et al., 2007; Kimura et al., 2015). Class V ER oxidoreductases are plant orthologs of mammalian P5 and have an a–a′–b domain structure. A rice class V ER oxidoreductase, consisting of PDIL2 and PDIL3, plays an important role in the accumulation of the seed storage protein Cys-rich 10-kD prolamin (crP10; Onda et al., 2011). Recombinant soybean class V ER oxidoreductase, GmPDIM and wheat class V ER oxidoreductase possess an oxidative refolding activity similar to that of class IV (Wadahama et al., 2008; Kimura et al., 2015). In the soybean, GmPDIL-1, GmPDIL-2, GmPDIM, GmPDIS-1, and GmPDIS-2 were found to associate transiently with a seed storage precursor protein, proglycinin, in the ER of the cotyledon by coimmunoprecipitation experiments, suggesting that multiple ER oxidoreductases are involved in the folding of the nascent proglycinin.The disulfide bond in the active center of ER oxidoreductases is reduced as a result of catalyzing disulfide bond formation in an unfolded protein. The reduced active center of PDI was discovered to be oxidized again by ER oxidoreductin-1 (Ero1p) in yeast (Frand and Kaiser, 1998; Pollard et al., 1998). Ero1p orthologs are present universally in eukaryotes. Yeast and flies have a single copy of the ERO1 gene, which is essential for survival (Frand and Kaiser, 1998; Pollard et al., 1998; Tien et al., 2008). Mammals have two genes encoding Ero1-α (Cabibbo et al., 2000) and Ero1-β (Pagani et al., 2000) that function as major disulfide donors to nascent proteins in the ER, but are not critical for survival (Zito et al., 2010). Domain a of yeast PDI is the most favored substrate of yeast Ero1p (Vitu et al., 2010), whereas a′ of human PDI is specifically oxidized by human Ero1-α (Chambers et al., 2010) and Ero1-β (Wang et al., 2011). Electrons from Cys residues of the active centers of PDI are transferred to oxygen by Ero1 (Tu and Weissman, 2004; Sevier and Kaiser, 2008). The reaction mechanisms of yeast Ero1p and human Ero1s have been intensively investigated; their regulation by PDI has been extensively studied as well (Tavender and Bulleid, 2010; Araki and Inaba, 2012; Benham et al., 2013; Ramming et al., 2015). Only rice Ero1 (OsERO1) has been identified as a plant ortholog of Ero1p (Onda et al., 2009). OsERO1 is necessary for disulfide bond formation in rice endosperm. The formation of native disulfide bonds in the major seed storage protein proglutelin was demonstrated to depend upon OsERO1 by RNAi knockdown experiments. However, no plant protein thiol disulfide oxidoreductases that are oxidized by a plant Ero1 ortholog have been identified to date.In this study, we show that multiple soybean ER oxidoreductases can be activated by a soybean Ero1 ortholog (GmERO1a). In addition, we propose a synergistic mechanism by which GmPDIM and GmPDIL-2 cooperatively fold unfolded proteins using oxidizing equivalents provided by GmERO1 in vitro.     

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.     

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.     

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.