Genome-Wide Association of Carbon and Nitrogen Metabolism in the Maize Nested Association Mapping Population

Carbon (C) and nitrogen (N) metabolism are critical to plant growth and development and are at the basis of crop yield and adaptation. We performed high-throughput metabolite analyses on over 12,000 samples from the nested association mapping population to identify genetic variation in C and N metabolism in maize (Zea mays ssp. mays). All samples were grown in the same field and used to identify natural variation controlling the levels of 12 key C and N metabolites, namely chlorophyll a, chlorophyll b, fructose, fumarate, glucose, glutamate, malate, nitrate, starch, sucrose, total amino acids, and total protein, along with the first two principal components derived from them. Our genome-wide association results frequently identified hits with single-gene resolution. In addition to expected genes such as invertases, natural variation was identified in key C₄ metabolism genes, including carbonic anhydrases and a malate transporter. Unlike several prior maize studies, extensive pleiotropy was found for C and N metabolites. This integration of field-derived metabolite data with powerful mapping and genomics resources allows for the dissection of key metabolic pathways, providing avenues for future genetic improvement.Carbon (C) and nitrogen (N) metabolism are the basis for life on Earth. The production, balance, and tradeoffs of C and N metabolism are critical to all plant growth, yield, and local adaptation (Coruzzi and Bush, 2001; Coruzzi et al., 2007). In plants, there is a critical balance between the tissues that are producing energy (sources) and those using it (sinks), as the identities and locations of these vary through time and developmental stage (Smith et al., 2004). While a great deal of research has focused on the key genes and proteins involved in these processes (Wang et al., 1993; Kim et al., 2000; Takahashi et al., 2009), relatively little is known about the natural variation within a species that fine-tunes these processes in individual plants.In addition, a key aspect of core C metabolism involves the nature of plant photosynthesis. While the majority of plants use standard C₃ photosynthetic pathways, some, including maize (Zea mays) and many other grasses, use C₄ photosynthesis to concentrate CO₂ in bundle sheath cells to avoid wasteful photorespiration (Sage, 2004). Under some conditions (such as drought or high temperatures), C₄ photosynthesis is much more efficient than C₃ photosynthesis. Since these conditions are expected to become more prevalent in the near future due to climate change, various research groups are working to convert C₃ crop species to C₄ metabolism in order to boost crop production and food security (Sage and Zhu, 2011). Beyond this, better understanding of both C₃ and C₄ metabolic pathways will aid efforts to breed crops for superior yield, N-use efficiency, and other traits important for global food production.In the last two decades, quantitative trait locus (QTL) mapping, first with linkage analysis and later with association mapping, has been used to dissect C and N metabolism in several species, including Arabidopsis (Arabidopsis thaliana; Mitchell-Olds and Pedersen, 1998; Keurentjes et al., 2008; Lisec et al., 2008; Sulpice et al., 2009), tomato (Solanum lycopersicum; Schauer et al., 2006), and maize (Hirel et al., 2001; Limami et al., 2002; Zhang et al., 2006, 2010a, 2010b). These studies identified key genetic regions underlying variation in core C and N metabolism, many of which include candidate genes known to be involved in these processes.Previous studies of genetic variation for C and N metabolism are limited by the fact that they identified trait loci only through linkage mapping in artificial families or through association mapping across populations of unrelated individuals. Linkage mapping benefits from high statistical power due to many individuals sharing the same genotype at any given location, but it suffers from low resolution due to the limited number of generations (and hence recombination events) since the initial founders. Association mapping, in turn, enjoys high resolution due to the long recombination histories of natural populations but suffers from low power, since most genotypes occur in only a few individuals. In addition, many of these studies focused on C and N in artificial settings (e.g. greenhouses or growth chambers) instead of field conditions, running the risk that important genetic loci could be missed if the conditions do not include important (and potentially unknown) natural environmental variables.To address these issues and improve our understanding of C and N metabolism in maize, we used a massive and diverse germplasm resource, the maize nested association mapping (NAM) population (Buckler et al., 2009; McMullen et al., 2009), to evaluate genetic variation underlying the accumulation of 12 targeted metabolites in maize leaf tissue under field conditions. This population was formed by mating 25 diverse maize lines to the reference line, B73, and creating a 200-member biparental family from each of these crosses. The entire 5,000-member NAM population thus combines the strengths of both linkage and association mapping (McMullen et al., 2009), and it has been used to identify QTLs for important traits such as flowering time (Buckler et al., 2009), disease resistance (Kump et al., 2011; Poland et al., 2011), and plant architecture (Tian et al., 2011; Peiffer et al., 2013). Most importantly, this combination of power and resolution frequently resolves associations down to the single-gene level, even when using field-based data.The metabolites we profiled are key indicators of photosynthesis, respiration, glycolysis, and protein and sugar metabolism in the plant (Sulpice et al., 2009). By taking advantage of a robotized metabolic phenotyping platform (Gibon et al., 2004), we performed more than 100,000 assays across 12,000 samples, with two independent samples per experimental plot. Raw data and the best linear unbiased predictors (BLUPs) of these data were included as part of a study of general functional variation in maize (Wallace et al., 2014), but, to our knowledge, this is the first in-depth analysis of these metabolic data. We find strong correlations among several of the metabolites, and we also find extensive pleiotropy among the different traits. Many of the top QTLs are also near or within candidate genes relating to C and N metabolism, thus identifying targets for future breeding and selection. These results provide a powerful resource for those working with core C and N metabolism in plants and for improving maize performance in particular.     

Rhamnolipids Elicit Defense Responses and Induce Disease Resistance against Biotrophic, Hemibiotrophic, and Necrotrophic Pathogens That Require Different Signaling Pathways in Arabidopsis and Highlight a Central Role for Salicylic Acid

Plant resistance to phytopathogenic microorganisms mainly relies on the activation of an innate immune response usually launched after recognition by the plant cells of microbe-associated molecular patterns. The plant hormones, salicylic acid (SA), jasmonic acid, and ethylene have emerged as key players in the signaling networks involved in plant immunity. Rhamnolipids (RLs) are glycolipids produced by bacteria and are involved in surface motility and biofilm development. Here we report that RLs trigger an immune response in Arabidopsis (Arabidopsis thaliana) characterized by signaling molecules accumulation and defense gene activation. This immune response participates to resistance against the hemibiotrophic bacterium Pseudomonas syringae pv tomato, the biotrophic oomycete Hyaloperonospora arabidopsidis, and the necrotrophic fungus Botrytis cinerea. We show that RL-mediated resistance involves different signaling pathways that depend on the type of pathogen. Ethylene is involved in RL-induced resistance to H. arabidopsidis and to P. syringae pv tomato whereas jasmonic acid is essential for the resistance to B. cinerea. SA participates to the restriction of all pathogens. We also show evidence that SA-dependent plant defenses are potentiated by RLs following challenge by B. cinerea or P. syringae pv tomato. These results highlight a central role for SA in RL-mediated resistance. In addition to the activation of plant defense responses, antimicrobial properties of RLs are thought to participate in the protection against the fungus and the oomycete. Our data highlight the intricate mechanisms involved in plant protection triggered by a new type of molecule that can be perceived by plant cells and that can also act directly onto pathogens.In their environment, plants are challenged by potentially pathogenic microorganisms. In response, they express a set of defense mechanisms including preformed structural and chemical barriers, as well as an innate immune response quickly activated after microorganism perception (Boller and Felix, 2009). Plant innate immunity is triggered after recognition by pattern recognition receptors of conserved pathogen- or microbe-associated molecular patterns (PAMPs or MAMPs, respectively) or by plant endogenous molecules released by pathogen invasion and called danger-associated molecular patterns (Boller and Felix, 2009; Dodds and Rathjen, 2010). This first step of recognition leads to the activation of MAMP-triggered immunity (MTI). Successful pathogens can secrete effectors that interfere or suppress MTI, resulting in effector-triggered susceptibility. A second level of perception involves the direct or indirect recognition by specific receptors of pathogen effectors leading to effector-triggered immunity (ETI; Boller and Felix, 2009; Dodds and Rathjen, 2010). Whereas MTI and ETI are thought to involve common signaling network, ETI is usually quantitatively stronger than MTI and associated with more sustained and robust immune responses (Katagiri and Tsuda, 2010; Tsuda and Katagiri, 2010).The plant hormones, salicylic acid (SA), jasmonic acid (JA), and ethylene (ET) have emerged as key players in the signaling networks involved in MTI and ETI (Robert-Seilaniantz et al., 2007; Tsuda et al., 2009; Katagiri and Tsuda, 2010; Mersmann et al., 2010; Tsuda and Katagiri, 2010; Robert-Seilaniantz et al., 2011). Interactions between these signal molecules allow the plant to activate and/or modulate an appropriate spectrum of responses, depending on the pathogen lifestyle, necrotroph or biotroph (Glazebrook, 2005; Koornneef and Pieterse, 2008). It is assumed that JA and ET signaling pathways are important for resistance to necrotrophic fungi including Botrytis cinerea and Alternaria brassicicola (Thomma et al., 2001; Ferrari et al., 2003; Glazebrook, 2005). Infection of Arabidopsis (Arabidopsis thaliana) with B. cinerea causes the induction of the JA/ET responsive gene PLANT DEFENSIN1.2 (PDF1.2; Penninckx et al., 1996; Zimmerli et al., 2001). Induction of PDF1.2 by B. cinerea is blocked in ethylene-insensitive2 (ein2) and coronatine-insensitive1 (coi1) mutants that are respectively defective in ET and JA signal transduction pathways. Moreover, ein2 and coi1 plants are highly susceptible to B. cinerea infection (Thomma et al., 1998; Thomma et al., 1999). JA/ET-dependent responses do not seem to be usually induced during resistance to biotrophs, but they can be effective if they are stimulated prior to pathogen challenge (Glazebrook, 2005). Plants impaired in SA signaling are highly susceptible to biotrophic and hemibiotrophic pathogens. Following pathogen infection, SA hydroxylase (NahG), enhanced disease susceptibility5 (eds5), or SA induction-deficient2 (sid2) plants are unable to accumulate high SA levels and they display heightened susceptibility to Pseudomonas syringae pv tomato (Pst), Hyaloperonospora arabidopsidis, or Erysiphe orontii (Delaney et al., 1994; Lawton et al., 1995; Wildermuth et al., 2001; Nawrath et al., 2002; Vlot et al., 2009). Mutants that are insensitive to SA, such as nonexpressor of PATHOGENESIS-RELATED (PR) genes1 (npr1), have enhanced susceptibility to these pathogens (Cao et al., 1994; Glazebrook et al., 1996; Shah et al., 1997; Dong, 2004). According to some reports, plant defense against necrotrophs also involves SA. Arabidopsis plants expressing the nahG gene and infected with B. cinerea show larger lesions compared with wild-type plants (Govrin and Levine, 2002). In tobacco (Nicotiana tabacum), acidic isoforms of PR3 and PR5 gene that are specifically induced by SA (Ménard et al., 2004) are up-regulated after challenge by B. cinerea (El Oirdi et al., 2010). Resistance to some necrotrophs like Fusarium graminearum involves both SA and JA signaling pathways (Makandar et al., 2010). It is assumed that SA and JA signaling can be antagonistic (Bostock, 2005; Koornneef and Pieterse, 2008; Pieterse et al., 2009; Thaler et al., 2012). In Arabidopsis, SA inhibits JA-dependent resistance against A. brassicicola or B. cinerea (Spoel et al., 2007; Koornneef et al., 2008). Recent studies demonstrated that ET modulates the NPR1-mediated antagonism between SA and JA (Leon-Reyes et al., 2009; Leon-Reyes et al., 2010a) and suppression by SA of JA-responsive gene expression is targeted at a position downstream of the JA biosynthesis pathway (Leon-Reyes et al., 2010b). Synergistic effects of SA- and JA-dependent signaling are also well documented (Schenk et al., 2000; van Wees et al., 2000; Mur et al., 2006) and induction of some defense responses after pathogen challenge requires intact JA, ET, and SA signaling pathways (Campbell et al., 2003).Isolated MAMPs trigger defense responses that also require the activation of SA, JA, and ET signaling pathways (Tsuda et al., 2009; Katagiri and Tsuda, 2010). For instance, treatment with the flagellin peptide flg22 induces many SA-related genes including SID2, EDS5, NPR1, and PR1 (Ferrari et al., 2007; Denoux et al., 2008), causes SA accumulation (Tsuda et al., 2008; Wang et al., 2009), and activates ET signaling (Bethke et al., 2009; Mersmann et al., 2010). Local application of lipopolysaccharides elevates the level of SA (Mishina and Zeier, 2007). The oomycete Pep13 peptide induces defense responses in potato (Solanum tuberosum) that require both SA and JA (Halim et al., 2009). Although signaling networks induced by isolated MAMPs are well documented, the contribution of SA, JA, and ET in MAMP- or PAMP-induced resistance to biotrophs and necrotrophs is poorly understood.Rhamnolipids (RLs) are glycolipids produced by various bacteria species including some Pseudomonas and Burkholderia species. They are essential for bacterial surface motility and biofilm development (Vatsa et al., 2010; Chrzanowski et al., 2012). RLs are potent stimulators of animal immunity (Vatsa et al., 2010). They have recently been shown to elicit plant defense responses and to induce resistance against B. cinerea in grapevine (Vitis vinifera; Varnier et al., 2009). They also participate to biocontrol activity of the plant beneficial bacteria Pseudomonas aeruginosa PNA1 against oomycetes (Perneel et al., 2008). However, the signaling pathways used by RLs to stimulate plant innate immunity are not known. To gain more insights into RL-induced MTI, we investigated RL-triggered defense responses and resistance to the necrotrophic fungus B. cinerea, the biotroph oomycete H. arabidopsidis, and the hemibiotroph bacterium Pst in Arabidopsis. Our results show that RLs trigger an innate immune response in Arabidopsis that protects the plant against these different lifestyle pathogens. We demonstrate that RL-mediated resistance involves separated signaling sectors that depend on the type of pathogen. In plants challenged by RLs, SA has a central role and participates to the restriction of the three pathogens. ET is fully involved in RL-induced resistance to the biotrophic oomycete and to the hemibiotrophic bacterium whereas JA is essential for the resistance to the necrotrophic fungus.     

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.     

Mitochondrial Alternative Oxidase Maintains Respiration and Preserves Photosynthetic Capacity during Moderate Drought in Nicotiana tabacum

The mitochondrial electron transport chain includes an alternative oxidase (AOX) that is hypothesized to aid photosynthetic metabolism, perhaps by acting as an additional electron sink for photogenerated reductant or by dampening the generation of reactive oxygen species. Gas exchange, chlorophyll fluorescence, photosystem I (PSI) absorbance, and biochemical and protein analyses were used to compare respiration and photosynthesis of Nicotiana tabacum ‘Petit Havana SR1’ wild-type plants with that of transgenic AOX knockdown (RNA interference) and overexpression lines, under both well-watered and moderate drought-stressed conditions. During drought, AOX knockdown lines displayed a lower rate of respiration in the light than the wild type, as confirmed by two independent methods. Furthermore, CO₂ and light response curves indicated a nonstomatal limitation of photosynthesis in the knockdowns during drought, relative to the wild type. Also relative to the wild type, the knockdowns under drought maintained PSI and PSII in a more reduced redox state, showed greater regulated nonphotochemical energy quenching by PSII, and displayed a higher relative rate of cyclic electron transport around PSI. The origin of these differences may lie in the chloroplast ATP synthase amount, which declined dramatically in the knockdowns in response to drought. None of these effects were seen in plants overexpressing AOX. The results show that AOX is necessary to maintain mitochondrial respiration during moderate drought. In its absence, respiration rate slows and the lack of this electron sink feeds back on the photosynthetic apparatus, resulting in a loss of chloroplast ATP synthase that then limits photosynthetic capacity.The plant mitochondrial electron transport chain (ETC) is bifurcated such that electrons in the ubiquinone pool partition between the cytochrome (cyt) pathway (consisting of Complex III, cyt c, and Complex IV) and alternative oxidase (AOX; Finnegan et al., 2004; Millar et al., 2011; Vanlerberghe, 2013). AOX directly couples ubiquinol oxidation with O₂ reduction to water. This reduces the energy yield of respiration because, unlike Complexes III and IV, AOX is not proton pumping. Hence, AOX is an electron sink, the capacity of which is little encumbered by rates of ATP turnover. In this way, AOX might be well suited to prevent cellular over-reduction. Supporting this, transgenic Nicotiana tabacum leaves with suppressed amounts of AOX have increased concentrations of mitochondrial-localized superoxide radical (O₂−) and nitric oxide, the products that can arise when an over-reduced ETC results in electron leakage to O₂ or nitrite (Cvetkovska and Vanlerberghe, 2012, 2013).In angiosperms, AOX is encoded by a small gene family (Considine et al., 2002). In Arabidopsis (Arabidopsis thaliana), mutation or knockdown of the stress-responsive AOX1a gene family member dramatically reduces AOX protein and the capacity of the AOX respiration pathway to consume O₂. Several studies have shown that this loss of AOX capacity in Arabidopsis aox1a plants affected processes such as growth, carbon and energy metabolism, and/or the cellular network of reactive oxygen species (ROS) scavengers (Fiorani et al., 2005; Umbach et al., 2005; Watanabe et al., 2008; Giraud et al., 2008; Skirycz et al., 2010). However, in studies in which respiration was measured, it was consistently reported that the lack of AOX capacity had no significant impact on the respiration rate in the dark (R_D; Umbach et al., 2005; Giraud et al., 2008; Strodtkötter et al., 2009; Florez-Sarasa et al., 2011; Yoshida et al., 2011b; Gandin et al., 2012). The exceptions are two reports that R_D was actually higher in aox1a than in the wild type under some conditions (Watanabe et al., 2008; Vishwakarma et al., 2014). To our knowledge, how the lack of AOX affects respiration rate in the light (R_L) is not reported in Arabidopsis or other species.Numerous studies have established the importance of mitochondrial metabolism in the light to optimize photosynthesis (Hoefnagel et al., 1998; Raghavendra and Padmasree, 2003). In recent years, the potential importance of specifically AOX respiration during photosynthesis has been examined using the Arabidopsis aox1a plants (Giraud et al., 2008; Strodtkötter et al., 2009; Zhang et al., 2010; Florez-Sarasa et al., 2011; Yoshida et al., 2011a, 2011b). In general, these studies reported small perturbations of photosynthesis in standard-grown aox1a plants, including slightly lower rates of CO₂ uptake or O₂ release (Gandin et al., 2012; Vishwakarma et al., 2014), slightly higher rates of cyclic electron transport (CET; Yoshida et al., 2011b), and slightly increased susceptibility to photoinhibition after a high light treatment (Florez-Sarasa et al., 2011). Generally, these studies concluded that aox1a plants exhibit a biochemical limitation of photosynthesis, in line with the hypothesis that AOX serves as a sink for excess photogenerated reducing power, with the reductant likely reaching the mitochondrion via the malate valve (Noguchi and Yoshida, 2008; Taniguchi and Miyake, 2012). Similar to these Arabidopsis studies, we recently reported that well-watered N. tabacum AOX knockdowns grown at moderate irradiance display a slight reduced rate of photosynthesis (approximately 10%–15%) when measured at high irradiance. However, we established that the lower photosynthetic rate was the result of a stomatal rather than biochemical limitation of photosynthesis, and provided evidence that this stomatal limitation resulted from disrupted nitric oxide homeostasis within the guard cells of AOX knockdown plants (Cvetkovska et al., 2014).Drought is a common abiotic stress that can substantially curtail photosynthesis because stomatal closure, meant to conserve water, also restricts CO₂ availability to the Calvin cycle. Besides this well established stomatal limitation of photosynthesis, there may also be water deficit-sensitive biochemical components that contribute to the reduction of photosynthesis during drought. However, the nature of this biochemical limitation and the degree to which it contributes to the curtailment of photosynthesis during drought remain areas of active debate (Flexas et al., 2004; Lawlor and Tezara, 2009; Pinheiro and Chaves, 2011). Additional factors, such as patchy stomatal closure (Sharkey and Seemann, 1989; Gunasekera and Berkowitz, 1992) or changes in the conductance to CO₂ of mesophyll cells (Perez-Martin et al., 2009), can further complicate analyses of photosynthesis during drought.Metabolism can experience energy imbalances, when there is a mismatch between rates of synthesis and rates of utilization of ATP and/or NADPH, and the importance of mechanisms to minimize such imbalances has been emphasized (Cruz et al., 2005; Kramer and Evans, 2011; Vanlerberghe, 2013). For example, such imbalances may occur in the chloroplast when the use of ATP and NADPH by the Calvin cycle does not keep pace with the harvesting of light energy (Hüner et al., 2012). This can result in excess excitation energy that can damage photosynthetic components, perhaps through the generation of ROS (Asada, 2006; Noctor et al., 2014). Such a scenario has been hypothesized to underlie the development of the biochemical limitations of photosynthesis reported during drought (Lawlor and Tezara, 2009).In this study, we find that N. tabacum AOX knockdowns show a compromised rate of mitochondrial respiration in the light during moderate drought. This corresponds with a strong nonstomatal limitation of photosynthesis in these plants relative to the wild type, and we describe a biochemical basis for this photosynthetic limitation. The results indicate that AOX is a necessary electron sink to support photosynthesis during drought, a condition when the major photosynthetic electron sink, the Calvin cycle, is becoming limited by CO₂ availability.     

Polyamine Oxidase5 Regulates Arabidopsis Growth through Thermospermine Oxidase Activity

The major plant polyamines (PAs) are the tetraamines spermine (Spm) and thermospermine (T-Spm), the triamine spermidine, and the diamine putrescine. PA homeostasis is governed by the balance between biosynthesis and catabolism; the latter is catalyzed by polyamine oxidase (PAO). Arabidopsis (Arabidopsis thaliana) has five PAO genes, AtPAO1 to AtPAO5, and all encoded proteins have been biochemically characterized. All AtPAO enzymes function in the back-conversion of tetraamine to triamine and/or triamine to diamine, albeit with different PA specificities. Here, we demonstrate that AtPAO5 loss-of-function mutants (pao5) contain 2-fold higher T-Spm levels and exhibit delayed transition from vegetative to reproductive growth compared with that of wild-type plants. Although the wild type and pao5 are indistinguishable at the early seedling stage, externally supplied low-dose T-Spm, but not other PAs, inhibits aerial growth of pao5 mutants in a dose-dependent manner. Introduction of wild-type AtPAO5 into pao5 mutants rescues growth and reduces the T-Spm content, demonstrating that AtPAO5 is a T-Spm oxidase. Recombinant AtPAO5 catalyzes the conversion of T-Spm and Spm to triamine spermidine in vitro. AtPAO5 specificity for T-Spm in planta may be explained by coexpression with T-Spm synthase but not with Spm synthase. The pao5 mutant lacking T-Spm oxidation and the acl5 mutant lacking T-Spm synthesis both exhibit growth defects. This study indicates a crucial role for T-Spm in plant growth and development.Polyamines (PAs) are low-molecular mass aliphatic amines that are present in almost all living organisms. Cellular PA concentrations are governed primarily by the balance between biosynthesis and catabolism. In plants, the major PAs are the diamine putrescine (Put), the triamine spermidine (Spd), and the tetraamines spermine (Spm) and thermospermine (T-Spm; Kusano et al., 2008; Alcázar et al., 2010; Mattoo et al., 2010; Takahashi and Kakehi, 2010; Tiburcio et al., 2014). Put is synthesized from Orn by Orn decarboxylase and/or from Arg by three sequential reactions catalyzed by Arg decarboxylase (ADC), agmatine iminohydrolase, and N-carbamoylputrescine amidohydrolase. Arabidopsis (Arabidopsis thaliana) does not contain an ORNITHINE DECARBOXYLASE gene (Hanfrey et al., 2001) and synthesizes Put from Arg via the ADC pathway. Put is further converted to Spd via an aminopropyltransferase reaction catalyzed by spermidine synthase (SPDS). In this reaction, an aminopropyl residue is transferred to Put from decarboxylated S-adenosyl-Met, which is synthesized by S-adenosyl-Met decarboxylase (SAMDC; Kusano et al., 2008). Spd is then converted to Spm or T-Spm, reactions catalyzed in Arabidopsis by spermine synthase (SPMS; encoded by SPMS) or thermospermine synthase (encoded by Acaulis5 [ACL5]), respectively (Hanzawa et al., 2000; Knott et al., 2007; Kakehi et al., 2008; Naka et al., 2010). A recent review reports that T-Spm is ubiquitously present in the plant kingdom (Takano et al., 2012).The PA catabolic pathway has been extensively studied in mammals. Spm and Spd acetylation by Spd/Spm-N¹-acetyltransferase (Enzyme Commission no. 2.3.1.57) precedes the catabolism of PAs and is a rate-limiting step in the catabolic pathway (Wallace et al., 2003). A mammalian polyamine oxidase (PAO), which requires FAD as a cofactor, oxidizes N¹-acetyl Spm and N¹-acetyl Spd at the carbon on the exo-side of the N⁴-nitrogen to produce Spd and Put, respectively (Wang et al., 2001; Vujcic et al., 2003; Wu et al., 2003; Cona et al., 2006). Mammalian spermine oxidases (SMOs) perform oxidation of the carbon on the exo-side of the N⁴-nitrogen to produce Spd, 3-aminopropanal, and hydrogen peroxide (Vujcic et al., 2002; Cervelli et al., 2003; Wang et al., 2003). Thus, mammalian PAOs and SMOs are classified as back-conversion (BC)-type PAOs.In plants, Spm, T-Spm, and Spd are catabolized by PAO. Plant PAOs derived from maize (Zea mays) and barley (Hordeum vulgare) catalyze terminal catabolism (TC)-type reactions (Tavladoraki et al., 1998). TC-type PAOs oxidize the carbon at the endo-side of the N⁴-nitrogen of Spm and Spd to produce N-(3-aminopropyl)-4-aminobutanal and 4-aminobutanal, respectively, plus 1,3-diaminopropane and hydrogen peroxide (Cona et al., 2006; Angelini et al., 2008, 2010). The Arabidopsis genome contains five PAO genes, designated as AtPAO1 to AtPAO5. Four recombinant AtPAOs, AtPAO1 to AtPAO4, have been homogenously purified and characterized (Tavladoraki et al., 2006; Kamada-Nobusada et al., 2008; Moschou et al., 2008; Takahashi et al., 2010; Fincato et al., 2011, 2012). AtPAO1 to AtPAO4 possess activities that convert Spm (or T-Spm) to Spd, called partial BC, or they convert Spm (or T-Spm) first to Spd and subsequently to Put, called full BC. Ahou et al. (2014) report that recombinant AtPAO5 also catalyzes a BC-type reaction. Therefore, all Arabidopsis PAOs are BC-type enzymes (Kamada-Nobusada et al., 2008; Moschou et al., 2008; Takahashi et al., 2010; Fincato et al., 2011, 2012; Ahou et al., 2014). Four of the seven PAOs in rice (Oryza sativa; OsPAO1, OsPAO3, OsPAO4, and OsPAO5) catalyze BC-type reactions (Ono et al., 2012; Liu et al., 2014a), whereas OsPAO7 catalyzes a TC-type reaction (Liu et al., 2014b). OsPAO2 and OsPAO6 remain to be characterized, but may catalyze TC-type reactions based on their structural similarity with OsPAO7. Therefore, plants possess both TC-type and BC-type PAOs.PAs are involved in plant growth and development. Recent molecular genetic analyses in Arabidopsis indicate that metabolic blocks at the ADC, SPDS, or SAMDC steps lead to embryo lethality (Imai et al., 2004; Urano et al., 2005; Ge et al., 2006). Potato (Solanum tuberosum) plants with suppressed SAMDC expression display abnormal phenotypes (Kumar et al., 1996). It was also reported that hydrogen peroxide derived from PA catabolism affects root development and xylem differentiation (Tisi et al., 2011). These studies indicate that flux through metabolic and catabolic PA pathways is required for growth and development. The Arabidopsis acl5 mutant, which lacks T-Spm synthase activity, displays excessive differentiation of xylem tissues and a dwarf phenotype, especially in stems (Hanzawa et al., 2000; Kakehi et al., 2008, 2010). An allelic ACL5 mutant (thickvein [tkv]) exhibits a similar phenotype as that of acl5 (Clay and Nelson, 2005). These results indicate that T-Spm plays an important role in Arabidopsis xylem differentiation (Vera-Sirera et al., 2010; Takano et al., 2012).Here, we demonstrate that Arabidopsis pao5 mutants contain 2-fold higher T-Spm levels and exhibit aerial tissue growth retardation approximately 50 d after sowing compared with that of wild-type plants. Growth inhibition of pao5 stems and leaves at an early stage of development is induced by growth on media containing low T-Spm concentrations. Complementation of pao5 with AtPAO5 rescues T-Spm-induced growth inhibition. We confirm that recombinant AtPAO5 catalyzes BC of T-Spm (or Spm) to Spd. Our data strongly suggest that endogenous T-Spm levels in Arabidopsis are fine tuned, and that AtPAO5 regulates T-Spm homeostasis through a T-Spm oxidation pathway.     

Metabolic Flux Analysis of Plastidic Isoprenoid Biosynthesis in Poplar Leaves Emitting and Nonemitting Isoprene

The plastidic 2-C-methyl-d-erythritol-4-phosphate (MEP) pathway is one of the most important pathways in plants and produces a large variety of essential isoprenoids. Its regulation, however, is still not well understood. Using the stable isotope ¹³C-labeling technique, we analyzed the carbon fluxes through the MEP pathway and into the major plastidic isoprenoid products in isoprene-emitting and transgenic isoprene-nonemitting (NE) gray poplar (Populus × canescens). We assessed the dependence on temperature, light intensity, and atmospheric [CO₂]. Isoprene biosynthesis was by far (99%) the main carbon sink of MEP pathway intermediates in mature gray poplar leaves, and its production required severalfold higher carbon fluxes compared with NE leaves with almost zero isoprene emission. To compensate for the much lower demand for carbon, NE leaves drastically reduced the overall carbon flux within the MEP pathway. Feedback inhibition of 1-deoxy-d-xylulose-5-phosphate synthase activity by accumulated plastidic dimethylallyl diphosphate almost completely explained this reduction in carbon flux. Our data demonstrate that short-term biochemical feedback regulation of 1-deoxy-d-xylulose-5-phosphate synthase activity by plastidic dimethylallyl diphosphate is an important regulatory mechanism of the MEP pathway. Despite being relieved from the large carbon demand of isoprene biosynthesis, NE plants redirected only approximately 0.5% of this saved carbon toward essential nonvolatile isoprenoids, i.e. β-carotene and lutein, most probably to compensate for the absence of isoprene and its antioxidant properties.Isoprenoids represent the largest and most diverse group (over 50,000) of natural compounds and are essential in all living organisms (Gershenzon and Dudareva, 2007; Thulasiram et al., 2007). They are economically important for humans as flavor and fragrance, cosmetics, drugs, polymers for rubber, and precursors for the chemical industry (Chang and Keasling, 2006). The broad variety of isoprenoid products is formed from two building blocks, dimethylallyl diphosphate (DMADP) and isopentenyl diphosphate (IDP). In plants, the plastidic 2-C-methyl-d-erythritol-4-phosphate (MEP) pathway (Zeidler et al., 1997) produces physiologically and ecologically important volatile organic compounds (VOCs), the carotenoids (tetraterpenes; Giuliano et al., 2008; Cazzonelli and Pogson, 2010), diterpenes, the prenyl side-chains of chlorophylls (Chls) and plastoquinones, isoprenylated proteins, the phytohormones gibberellins, and side-chain of cytokinins (for review, see Dudareva et al., 2013; Moses et al., 2013). Industrially important prokaryotes (e.g. Escherichia coli) also use the MEP pathway for the biosynthesis of isoprenoids (Vranová et al., 2012), and there is an increasing interest in manipulating the MEP pathway of engineered microbes to increase production of economically relevant isoprenoids (Chang and Keasling, 2006). To achieve this, a mechanistic understanding of the regulation of the MEP pathway is needed (Vranová et al., 2012).Some plants, including poplars (Populus spp.), produce large amounts of the hemiterpene VOC isoprene. Worldwide isoprene emissions from plants are estimated to be 600 teragrams per year and to account for one-third of all hydrocarbons emitted to the atmosphere (Arneth et al., 2008; Guenther, 2013). Isoprene has strong effects on air chemistry and climate by participating in ozone formation reactions (Fuentes et al., 2000), by prolonging the lifespan of methane, a greenhouse gas (Poisson et al., 2000; Archibald et al., 2011), and by taking part in the formation of secondary organic aerosols (Kiendler-Scharr et al., 2012).Poplar leaves invest a significant amount of recently fixed carbon in isoprene biosynthesis (Delwiche and Sharkey, 1993; Schnitzler et al., 2010; Ghirardo et al., 2011) to cope with abiotic stresses (Sharkey, 1995; Velikova and Loreto, 2005; Behnke et al., 2007, 2010b, 2013; Vickers et al., 2009; Loreto and Schnitzler, 2010; Sun et al., 2013b), although there are indications that other protective mechanisms can partially compensate the lack of isoprene emission in genetically transformed poplars (Behnke et al., 2012; Way et al., 2013). It has been suggested that in isoprene-emitting (IE) species, most of the carbon that passes through the MEP pathway is used for isoprene biosynthesis (Sharkey and Yeh, 2001). However, a recent study using pulse-chase labeling with ¹⁴C has shown continuous synthesis and degradation of carotenes and Chl a in mature leaves of Arabidopsis (Arabidopsis thaliana; Beisel et al., 2010), and the amount of flux diverted to carotenoid and Chl synthesis compared with isoprene biosynthesis in poplar leaves is not known.Isoprene emission is temperature, light, and CO₂ dependent (Schnitzler et al., 2005; Rasulov et al., 2010; Way et al., 2011; Monson et al., 2012; Li and Sharkey, 2013a). It has been demonstrated that isoprene biosynthesis depends on the activities of IDP isomerase (EC 5.3.3.2), isoprene synthase (ISPS; EC 4.2.3.27), and the amount of ISPS substrate, DMADP (Brüggemann and Schnitzler, 2002a, 2002b; Schnitzler et al., 2005; Rasulov et al., 2009b). In turn, DMADP concentration has been hypothesized to act as a feedback regulator of the MEP pathway by inhibiting 1-deoxy-d-xylulose-5-phosphate synthase (DXS; EC 2.2.1.7), the first enzyme of the MEP pathway (Banerjee et al., 2013). Understanding the controlling mechanism of isoprene biosynthesis is not only of fundamental relevance, but also necessary for engineering the MEP pathway in various organisms and for accurate simulation of isoprene emissions by plants in predicting atmospheric reactivity (Niinemets and Monson, 2013).There is ample evidence that silencing the ISPS in poplar has a broad effect on the leaf metabolome (Behnke et al., 2009, 2010a, 2013; Way et al., 2011; Kaling et al., 2014). While some of those changes (e.g. ascorbate and α-tocopherol) are compensatory mechanisms to cope with abiotic stresses, others (e.g. shikimate pathway and phenolic compounds) might be related to the alteration of the MEP pathway (Way et al., 2013; Kaling et al., 2014). The perturbation of these metabolic pathways can be attributed to the removal of a major carbon sink of the MEP pathway and the resulting change in the energy balance within the plant cell (Niinemets et al., 1999; Ghirardo et al., 2011). In this work, we analyzed the carbon fluxes through the MEP pathway into the main plastidic isoprenoid products.We used the ¹³C-labeling technique as a tool to measure the carbon fluxes through the MEP pathway at different temperatures, light intensities, and CO₂ concentrations in mature leaves of IE and transgenic, isoprene-nonemitting (NE) gray poplar (Populus × canescens). Isoprene emission was drastically reduced in the transgenic trees through knockdown of PcISPS gene expression by RNA interference, resulting in plants with only 1% to 5% of isoprene emission potential compared with wild-type plants (Behnke et al., 2007).We measured the appearance of ¹³C in the isoprenoid precursors 2-C-methyl-d-erythritol-2,4-cyclodiphosphate (MEcDP) and DMADP as well as isoprene and the major downstream products of the MEP pathway, i.e. carotenoids and Chls. To reliably detect de novo synthesis of the pigments, which occur at very low rates (Beisel et al., 2010), we used isotope ratio mass spectrometry (IRMS).Here, (1) we quantify the effect of isoprene biosynthesis on the MEP pathway in poplar, and (2) we show that suppression of isoprene biosynthesis negatively affects the carbon flux through the MEP pathway by accumulating plastidic DMADP, which feeds back to inhibit PcDXS, leading to (3) a slight increase of carbon flux toward production of greater chain-length isoprenoids and (4) a strong decrease in the overall isoprenoid carbon fluxes to compensate for the much lower MEP pathway demand for carbon. This study strongly supports the hypothesis that an important regulatory mechanism of the MEP pathway is the feedback regulation of plastidic DMADP on DXS. The large carbon flux through the MEP pathway of IE poplar plastids demonstrates the potential of transgenically altered IE plant species to produce economically valuable isoprenoids at high rates in, for instance, industrial applications.     

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 Endocytosis of Cellulose Synthase in Arabidopsis Is Dependent on μ2, a Clathrin-Mediated Endocytosis Adaptin

Clathrin-mediated endocytosis (CME) is the best-characterized type of endocytosis in eukaryotic cells. Plants appear to possess all of the molecular components necessary to carry out CME; however, functional characterization of the components is still in its infancy. A yeast two-hybrid screen identified μ2 as a putative interaction partner of CELLULOSE SYNTHASE6 (CESA6). Arabidopsis (Arabidopsis thaliana) μ2 is homologous to the medium subunit 2 of the mammalian ADAPTOR PROTEIN COMPLEX2 (AP2). In mammals, the AP2 complex acts as the central hub of CME by docking to the plasma membrane while concomitantly recruiting cargo proteins, clathrin triskelia, and accessory proteins to the sites of endocytosis. We confirmed that μ2 interacts with multiple CESA proteins through the μ-homology domain of μ2, which is involved in specific interactions with endocytic cargo proteins in mammals. Consistent with its role in mediating the endocytosis of cargos at the plasma membrane, μ2-YELLOW FLUORESCENT PROTEIN localized to transient foci at the plasma membrane, and loss of μ2 resulted in defects in bulk endocytosis. Furthermore, loss of μ2 led to increased accumulation of YELLOW FLUORESCENT PROTEIN-CESA6 particles at the plasma membrane. Our results suggest that CESA represents a new class of CME cargo proteins and that plant cells might regulate cellulose synthesis by controlling the abundance of active CESA complexes at the plasma membrane through CME.Cellulose microfibrils, as the major load-bearing polymers in plant cell walls, are the predominant component that enforces asymmetric cell expansion (Green, 1962). In higher plants, cellulose is synthesized by multimeric rosettes, which are also referred to as cellulose synthase complexes (CSCs; Kimura et al., 1999). Genetic and coimmunoprecipitation studies have indicated that CELLULOSE SYNTHASE1 (CESA1), CESA3, and CESA6-like (CESA6, CESA2, CESA5, and CESA9) isoforms are constituents of CSCs during primary cell wall synthesis (Persson et al., 2005; Desprez et al., 2007; Persson et al., 2007; Wang et al., 2008), whereas CESA4, CESA7, and CESA8 are implicated in the cellulose synthesis of secondary cell walls (Taylor et al., 1999, 2003; Brown et al., 2005). Knowledge about cellulose synthesis has recently been enhanced by the development of a system whereby the dynamics of CESA can be imaged in living cells (Paredez et al., 2006; Desprez et al., 2007). In agreement with earlier transmission electron microscopy studies in which rosettes were visualized in Golgi cisternae, vesicles, and at the plasma membrane (Haigler and Brown, 1986), fluorescent protein tagging of CESA has identified CESA localization at the plasma membrane, in Golgi bodies, and in small intracellular compartments (Paredez et al., 2006; Desprez et al., 2007; Crowell et al., 2009; Gutierrez et al., 2009; Gu et al., 2010; Lei et al., 2012; Li et al., 2012b).Assuming that cellulose synthesis occurs solely at the plasma membrane, the trafficking of CSCs to and from the plasma membrane may act as a significant regulatory mechanism. Although the mechanistic details of CESA trafficking are lacking, live cell imaging has shown that CESA localizes to various subcellular compartments. A subset of CESAs colocalize with markers of the trans-Golgi network (TGN)/early endosome (EE), an organelle that is part of both the secretory and endocytic pathways in Arabidopsis (Arabidopsis thaliana; Dettmer et al., 2006; Lam et al., 2007; Crowell et al., 2009, 2010; Viotti et al., 2010). CESAs also localize to microtubule-associated cellulose synthase compartments (MASCs) and small CESA-containing compartments (SmaCCs). The exact function of SmaCCs/MASCs is unknown, but it has been proposed that SmaCCs/MASCs might result from the internalization of CSCs or might act in the delivery of CSCs to the plasma membrane (Crowell et al., 2009, 2010; Gutierrez et al., 2009).Clathrin-mediated endocytosis (CME) has been shown to be a major endocytic pathway in Arabidopsis (Holstein, 2002; Samaj et al., 2005; Dhonukshe et al., 2007; Kleine-Vehn and Friml, 2008; Chen et al., 2011; Beck et al., 2012; Wang et al., 2013), although there is also evidence of clathrin-independent endocytosis mechanisms (Bandmann and Homann, 2012). The function of many CME proteins has been extensively characterized in mammals (McMahon and Boucrot, 2011), and homologs of many CME components are encoded by the Arabidopsis genome, including multiple copies of clathrin H chain and clathrin light chain (CLC), all four subunits of the heterotetrameric ADAPTOR PROTEIN COMPLEX2 (AP2) complex, dynamin-related proteins, and accessory proteins such as AP180 (Holstein, 2002; Chen et al., 2011); however, many CME components have yet to be characterized in plants.It has been suggested that CME might also function in controlling cell wall metabolism. For example, dividing and growing cells internalize cross-linked cell wall pectins, which might allow for cell wall remodeling (Baluska et al., 2002, 2005; Samaj et al., 2004). Moreover, the importance of endocytosis for cell wall morphogenesis is apparent from the functional characterization of proteins involved in CME. A dynamin-related protein, DRP1A, plays a significant role in endocytosis and colocalizes with CLC (Collings et al., 2008; Konopka and Bednarek, 2008). Defective endocytosis in RADIAL SWELLING9 (rsw9) plants, which contain a mutation in DRP1A, results in cellulose deficiency and defects in cell elongation (Collings et al., 2008). A mutation in rice, brittle culm3 (bc3), was mapped to the dynamin-related gene OsDRP2A, which has been proposed to function in CME. The brittle-culm phenotype in this mutant was attributed to cellulose deficiency (Xiong et al., 2010). Although the abundance of OsCESA4 was also altered in bc3, it remains unclear whether the cellulose deficiency of either bc3 or rsw9 results directly from perturbations in CESA trafficking.To identify proteins involved in the regulation of cellulose biosynthesis, a yeast two-hybrid (Y2H) screen was performed in which the central domain of CESA6 (CESA6CD) was used as bait to screen an Arabidopsis complementary DNA library for potential interaction partners of CESA6 (Gu et al., 2010; Gu and Somerville, 2010). The Y2H screen identified μ2 as a putative interaction partner of CESA6CD. The mammalian homolog of μ2 is the medium subunit of the tetrameric AP2, which acts as the core of the CME machinery by docking to the plasma membrane while concomitantly recruiting cargo proteins, clathrin triskelia, and accessory proteins to the sites of endocytosis (Jackson et al., 2010; McMahon and Boucrot, 2011; Cocucci et al., 2012). In this study, we provide evidence that μ2 plays a role in CME in Arabidopsis, that CESAs are a new set of CME cargo proteins, and that plant cells might regulate cellulose synthesis by controlling the abundance of CSCs at the plasma membrane through CME. To our knowledge, this study is the first to show the affect of an AP2 complex component on endocytosis in plants and the first to visualize an AP2 complex component in living plant cells. Furthermore, our data suggest that the role of AP2 in plants may differ from what has been shown in animals.     

SALT-OVERLY SENSITIVE5 Mediates Arabidopsis Seed Coat Mucilage Adherence and Organization through Pectins

Interactions between cell wall polymers are critical for establishing cell wall integrity and cell-cell adhesion. Here, we exploit the Arabidopsis (Arabidopsis thaliana) seed coat mucilage system to examine cell wall polymer interactions. On hydration, seeds release an adherent mucilage layer strongly attached to the seed in addition to a nonadherent layer that can be removed by gentle agitation. Rhamnogalacturonan I (RG I) is the primary component of adherent mucilage, with homogalacturonan, cellulose, and xyloglucan constituting minor components. Adherent mucilage contains rays composed of cellulose and pectin that extend above the center of each epidermal cell. CELLULOSE SYNTHASE5 (CESA5) and the arabinogalactan protein SALT-OVERLY SENSITIVE5 (SOS5) are required for mucilage adherence through unknown mechanisms. SOS5 has been suggested to mediate adherence by influencing cellulose biosynthesis. We, therefore, investigated the relationship between SOS5 and CESA5. cesa5-1 seeds show reduced cellulose, RG I, and ray size in adherent mucilage. In contrast, sos5-2 seeds have wild-type levels of cellulose but completely lack adherent RG I and rays. Thus, relative to each other, cesa5-1 has a greater effect on cellulose, whereas sos5-2 mainly affects pectin. The double mutant cesa5-1 sos5-2 has a much more severe loss of mucilage adherence, suggesting that SOS5 and CESA5 function independently. Double-mutant analyses with mutations in MUCILAGE MODIFIED2 and FLYING SAUCER1 that reduce mucilage release through pectin modification suggest that only SOS5 influences pectin-mediated adherence. Together, these findings suggest that SOS5 mediates adherence through pectins and does so independently of but in concert with cellulose synthesized by CESA5.Cellulosic cell walls are a defining feature of land plants. Primary cell walls are composed of three major classes of polysaccharides: cellulose, hemicelluloses, and pectins. In addition, approximately 10% of the primary cell wall is composed of protein (Burton et al., 2010). Cell walls provide mechanical support for the cell, and cell wall carbohydrates in the middle lamellae mediate cell-cell adhesion (Caffall and Mohnen, 2009). Current models of cell wall structure depict a cellulose-hemicellulose network embedded in an independent pectin gel (for review, see Albersheim et al., 2011). These components are believed to interact through both covalent and noncovalent bonds to provide structure and strength to the cell wall, although the relative importance of pectin and its interactions with the hemicellulose-cellulose network remain unclear (for review, see Cosgrove, 2005).Another gap in our understanding of cell wall structure and assembly is the role of arabinogalactan proteins (AGPs). AGPs are a family of evolutionarily conserved secreted proteins highly glycosylated with type II arabinogalactans, and they can be localized to the plasma membrane by a C-terminal glycophosphatidylinositol (GPI) lipid anchor (for review, see Schultz et al., 2000; Showalter, 2001; Johnson et al., 2003; Seifert and Roberts, 2007; Ellis et al., 2010). AGPs can be extensively modified in the cell wall; many glycosyl hydrolases can affect AGP function by cleaving their glycosyl side chains (Sekimata et al., 1989; Cheung et al., 1995; Wu et al., 1995; Kotake et al., 2005). The GPI anchor can also be cleaved, releasing the AGPs from the membrane into the cell wall (Schultz et al., 2000). Although their exact roles are still unclear, AGPs have been proposed to interact with cell wall polysaccharides, initiate intracellular signaling cascades, and influence a wide variety of biological processes (for review, see Seifert and Roberts, 2007; Ellis et al., 2010; Tan et al., 2013).Many fasciclin-like AGPs (FLAs), which contain at least one fasciclin domain (FAS) associated with protein-protein interactions, have been suggested to influence cellulose biosynthesis or organization (Seifert and Roberts, 2007; Li et al., 2010; MacMillan et al., 2010). FLA3 RNA interference lines have reduced intine cell wall biosynthesis and loss of Calcofluor white (a fluorescent dye specific for glycan molecules) staining in aborted pollen grains (Li et al., 2010). A fla11 fla12 double mutant was shown to have reduced cellulose deposition, altered cellulose microfibril angle, and reduced cell wall integrity (MacMillan et al., 2010). The fla11 fla12 double mutant also had reductions in arabinans, galactans, and rhamnose (MacMillan et al., 2010). FLA4/SALT-OVERLY SENSITIVE5 (SOS5) was identified in a screen for salt sensitivity in roots. The SOS5 gene encodes an FLA protein with a GPI anchor, two AGP-like domains, and two FAS domains (Shi et al., 2003). Plants homozygous for the loss-of-function conditional allele sos5-1 have thinner root cell walls that appear less organized (Shi et al., 2003). The presence of the two FAS domains has led to the suggestion that SOS5 may interact with other proteins, forming a network that strengthens the cell wall (Shi et al., 2003). SOS5 is involved in regulation of cell wall rheology through a pathway involving two Leu-rich repeat receptor-like kinases, FEI1 and FEI2 (Xu et al., 2008). SOS5 and FEI2 are also required for normal seed coat mucilage adherence and hypothesized to do so by influencing cellulose biosynthesis (Harpaz-Saad et al., 2011, 2012).Arabidopsis (Arabidopsis thaliana) seed coat mucilage is a powerful model for studying cell wall biosynthesis and polysaccharide interactions (Arsovski et al., 2010; Haughn and Western, 2012). Seed coat epidermal cells sequentially produce two distinct types of secondary cell walls with unique morphologies and properties (Western et al., 2000; Windsor et al., 2000). Between approximately 5 and 9 d approximate time of fertilization (DPA), seed coat epidermal cells synthesize mucilage and deposit it in the apoplast, creating a donut-shaped mucilage pocket that surrounds a central cytoplasmic column (Western et al., 2000, 2004; Haughn and Chaudhury, 2005). From 9 to 13 DPA, the cytoplasmic column is gradually replaced by a cellulose-rich, volcano-shaped secondary cell wall called the columella (Beeckman et al., 2000; Western et al., 2000; Windsor et al., 2000; Stork et al., 2010; Mendu et al., 2011).Seed mucilage is composed primarily of relatively unbranched rhamnogalacturonan I (RG I) with minor amounts of homogalacturonan (HG), cellulose, and hemicelluloses (for review, see Haughn and Western, 2012). When mucilage is hydrated, it expands rapidly from the apoplastic pocket, forming a halo that surrounds the seed. Mucilage separates into two fractions: a loose nonadherent fraction and an inner adherent fraction that can only be released by vigorous shaking, strong bases, or glycosidases (for review, see North et al., 2014). Galactans and arabinans are also present in mucilage, and their regulation by glycosidases is required for correct mucilage hydration (Dean et al., 2007; Macquet et al., 2007b; Arsovski et al., 2009). For example, β-XYLOSIDASE1 encodes a bifunctional β-d-xylosidase/α-l-arabinofuranosidase required for arabinan modification in mucilage, and β-xylosidase1 mutant seeds have a delayed mucilage release phenotype (Arsovski et al., 2009). MUCILAGE MODIFIED2 (MUM2) encodes a β-d-galactosidase, and mum2 seeds fail to release mucilage when hydrated in water (Dean et al., 2007; Macquet et al., 2007b). MUM2 is believed to modify RG I galactan side chains but may also affect the galactan component of other mucilage components (Dean et al., 2007; Macquet et al., 2007b). Galactans are capable of binding to cellulose in vitro and could affect mucilage hydration through pectin-cellulose interactions (Zykwinska et al., 2005, 2007a, 2007b; Dick-Pérez et al., 2011; Wang et al., 2012), although carbohydrate linkage analysis suggests that the galactan side chains are very short.Several studies indicate that seed mucilage extrusion and expansion are also influenced by methylesterification of HG. For example, both SUBTILISIN-LIKE SER PROTEASE1.7 and PECTIN METHYLESTERASE INHIBITOR6 are required for proper methyl esterification of mucilage (Rautengarten et al., 2008; Saez-Aguayo et al., 2013). Mutations in another gene, FLYING SAUCER1 (FLY1; a transmembrane E3 ubiquitin ligase), reduce the degree of pectin methylesterification in mucilage and cause increased mucilage adherence and defective mucilage extrusion (Voiniciuc et al., 2013). fly1 seeds have disc-like structures at the edge of the mucilage halo, which are outer primary cell wall fragments that detach from the columella during extrusion and are difficult to separate from the adherent mucilage (Voiniciuc et al., 2013).Recently, CELLULOSE SYNTHASE5 (CESA5) and SOS5 were proposed to facilitate cellulose-mediated mucilage adherence (Harpaz-Saad et al., 2011; Mendu et al., 2011; Sullivan et al., 2011). A simple hypothesis for the role of CESA5 in mucilage adherence is that it synthesizes cellulose, which interacts with the mucilage pectin to mediate adherence. Loss of CESA5 function results in a reduction of mucilage cellulose biosynthesis and a less adherent mucilage cell wall matrix (Mendu et al., 2011; Sullivan et al., 2011). The role of SOS5 in mucilage adherence is more difficult to explain. SOS5 null mutations cause a loss-of-adherence phenotype similar to cesa5-1 seeds, suggesting that SOS5 may regulate mucilage adherence by influencing CESA5 function (Harpaz-Saad et al., 2011). However, the mechanism through which SOS5 could influence CESA5 and/or cellulose biosynthesis is not clear.To better understand the role of SOS5 in mucilage adherence and its relationship to CESA5, we thoroughly investigated the seed coat epidermal cell phenotypes of the cesa5-1 and sos5-2 single mutants as well as those of the cesa5-1 sos5-2 double mutant. We also investigated how cellulose, SOS5, and pectin interact to mediate mucilage adherence by constructing double mutants with either cesa5-1 or sos5-2 together with either mum2-1 or fly1. Our results suggest that SOS5 mediates mucilage adherence independently of CESA5. Furthermore, compared with CESA5, SOS5 has a greater influence on mucilage pectin structure, suggesting that SOS5 mediates mucilage adherence through pectins, not cellulose.