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.     

CELLULOSE SYNTHASE-LIKE A2, a Glucomannan Synthase,Is Involved in Maintaining Adherent Mucilage Structure in Arabidopsis Seed

Mannans are hemicellulosic polysaccharides that are considered to have both structural and storage functions in the plant cell wall. However, it is not yet known how mannans function in Arabidopsis (Arabidopsis thaliana) seed mucilage. In this study, CELLULOSE SYNTHASE-LIKE A2 (CSLA2; At5g22740) expression was observed in several seed tissues, including the epidermal cells of developing seed coats. Disruption of CSLA2 resulted in thinner adherent mucilage halos, although the total amount of the adherent mucilage did not change compared with the wild type. This suggested that the adherent mucilage in the mutant was more compact compared with that of the wild type. In accordance with the role of CSLA2 in glucomannan synthesis, csla2-1 mucilage contained 30% less mannosyl and glucosyl content than did the wild type. No appreciable changes in the composition, structure, or macromolecular properties were observed for nonmannan polysaccharides in mutant mucilage. Biochemical analysis revealed that cellulose crystallinity was substantially reduced in csla2-1 mucilage; this was supported by the removal of most mucilage cellulose through treatment of csla2-1 seeds with endo-β-glucanase. Mutation in CSLA2 also resulted in altered spatial distribution of cellulose and an absence of birefringent cellulose microfibrils within the adherent mucilage. As with the observed changes in crystalline cellulose, the spatial distribution of pectin was also modified in csla2-1 mucilage. Taken together, our results demonstrate that glucomannans synthesized by CSLA2 are involved in modulating the structure of adherent mucilage, potentially through altering cellulose organization and crystallization.Mannan polysaccharides are a complex set of hemicellulosic cell wall polymers that are considered to have both structural and storage functions. Based on the particular chemical composition of the backbone and the side chains, mannan polysaccharides are classified into four types: pure mannan, glucomannan, galactomannan, and galactoglucomannan (Moreira and Filho, 2008; Wang et al., 2012; Pauly et al., 2013). Each of these polysaccharides is composed of a β-1,4-linked backbone containing Man or a combination of Glc and Man residues. In addition, the mannan backbone can be substituted with side chains of α-1,6-linked Gal residues. Mannan polysaccharides have been proposed to cross link with cellulose and other hemicelluloses via hydrogen bonds (Fry, 1986; Iiyama et al., 1994; Obel et al., 2007; Scheller and Ulvskov, 2010). Furthermore, it has been reported that heteromannans with different levels of substitution can interact with cellulose in diverse ways (Whitney et al., 1998). Together, these observations indicate the complexity of mannan polysaccharides in the context of cell wall architecture.CELLULOSE SYNTHASE-LIKE A (CSLA) enzymes have been shown to have mannan synthase activity in vitro. These enzymes polymerize the β-1,4-linked backbone of mannans or glucomannans, depending on the substrates (GDP-Man and/or GDP-Glc) provided (Richmond and Somerville, 2000; Liepman et al., 2005, 2007; Pauly et al., 2013). In Arabidopsis (Arabidopsis thaliana), nine CSLA genes have been identified; different CSLAs are responsible for the synthesis of different mannan types (Liepman et al., 2005, 2007). CSLA7 has mannan synthase activity in vitro (Liepman et al., 2005) and has been shown to synthesize stem glucomannan in vivo (Goubet et al., 2009). Disrupting the CSLA7 gene results in defective pollen growth and embryo lethality phenotypes in Arabidopsis, indicating structural or signaling functions of mannan polysaccharides during plant embryo development (Goubet et al., 2003). A mutation in CSLA9 results in the inhibition of Agrobacterium tumefaciens-mediated root transformation in the rat4 mutant (Zhu et al., 2003). CSLA2, CSLA3, and CSLA9 are proposed to play nonredundant roles in the biosynthesis of stem glucomannans, although mutations in CSLA2, CSLA3, or CSLA9 have no effect on stem development or strength (Goubet et al., 2009). All of the Arabidopsis CSLA proteins have been shown to be involved in the biosynthesis of mannan polysaccharides in the plant cell wall (Liepman et al., 2005, 2007), although the precise physiological functions of only CSLA7 and CSLA9 have been conclusively demonstrated.In Arabidopsis, when mature dry seeds are hydrated, gel-like mucilage is extruded to envelop the entire seed. Ruthenium red staining of Arabidopsis seeds reveals two different mucilage layers, termed the nonadherent and the adherent mucilage layers (Western et al., 2000; Macquet et al., 2007a). The outer, nonadherent mucilage is loosely attached and can be easily extracted by shaking seeds in water. Compositional and linkage analyses suggest that this layer is almost exclusively composed of unbranched rhamnogalacturonan I (RG-I) (>80% to 90%), with small amounts of branched RG-I, arabinoxylan, and high methylesterified homogalacturonan (HG). By contrast, the inner, adherent mucilage layer is tightly attached to the seed and can only be removed by strong acid or base treatment, or by enzymatic digestion (Macquet et al., 2007a; Huang et al., 2011; Walker et al., 2011). As with the nonadherent layer, adherent mucilage is also mainly composed of unbranched RG-I, but with small numbers of arabinan and galactan ramifications (Penfield et al., 2001; Willats et al., 2001; Dean et al., 2007; Macquet et al., 2007a, 2007b; Arsovski et al., 2009; Haughn and Western, 2012). There are also minor amounts of pectic HG in the adherent mucilage, with high methylesterified HG in the external domain compared with the internal domain of the adherent layer (Willats et al., 2001; Macquet et al., 2007a; Rautengarten et al., 2008; Sullivan et al., 2011; Saez-Aguayo et al., 2013). In addition, the adherent mucilage contains cellulose (Blake et al., 2006; Macquet et al., 2007a), which is entangled with RG-I and is thought to anchor the pectin-rich mucilage onto seeds (Macquet et al., 2007a; Harpaz-Saad et al., 2011, 2012; Mendu et al., 2011; Sullivan et al., 2011). As such, Arabidopsis seed mucilage is considered to be a useful model for investigating the biosynthesis of cell wall polysaccharides and how this process is regulated in vivo (Haughn and Western, 2012).Screening for altered seed coat mucilage has led to the identification of several genes encoding enzymes that are involved in the biosynthesis or modification of mucilage components. RHAMNOSE SYNTHASE2/MUCILAGE-MODIFIED4 (MUM4) is responsible for the synthesis of UDP-l-Rha (Usadel et al., 2004; Western et al., 2004; Oka et al., 2007). The putative GALACTURONSYLTRANSFERASE11 can potentially synthesize mucilage RG-I or HG pectin from UDP-d-GalUA (Caffall et al., 2009). GALACTURONSYLTRANSFERASE-LIKE5 appears to function in the regulation of the final size of the mucilage RG-I (Kong et al., 2011, 2013). Mutant seeds defective in these genes display reduced thickness of the extruded mucilage layer compared with wild-type Arabidopsis seeds.RG-I deposited in the apoplast of seed coat epidermal cells appears to be synthesized in a branched form that is subsequently modified by enzymes in the apoplast. MUM2 encodes a β-galactosidase that removes Gal residues from RG-I side chains (Dean et al., 2007; Macquet et al., 2007b). β-XYLOSIDASE1 encodes an α-l-arabinfuranosidase that removes Ara residues from RG-I side chains (Arsovski et al., 2009). Disruptions of these genes lead to defective hydration properties and affect the extrusion of mucilage. Furthermore, correct methylesterification of mucilage HG is also required for mucilage extrusion. HG is secreted into the wall in a high methylesterified form that can then be enzymatically demethylesterified by pectin methylesterases (PMEs; Bosch and Hepler, 2005). PECTIN METHYLESTERASE INHIBITOR6 (PMEI6) inhibits PME activities (Saez-Aguayo et al., 2013). The subtilisin-like Ser protease (SBT1.7) can activate other PME inhibitors, but not PMEI6 (Rautengarten et al., 2008; Saez-Aguayo et al., 2013). Disruption of either PMEI6 or SBT1.7 results in the delay of mucilage release.Although cellulose is present at low levels in adherent mucilage, it plays an important adhesive role for the attachment of mucilage pectin to the seed coat epidermal cells. The orientation and amount of pectin associated with the cellulose network is largely determined by cellulose conformation properties (Macquet et al., 2007a; Haughn and Western, 2012). Previous studies have demonstrated that CELLULOSE SYNTHASE A5 (CESA5) is required for the production of seed mucilage cellulose and the adherent mucilage in the cesa5 mutant can be easily extracted with water (Harpaz-Saad et al., 2011, 2012; Mendu et al., 2011; Sullivan et al., 2011).Despite all of these discoveries, large gaps remain in the current knowledge of the biosynthesis and functions of mucilage polysaccharides in seed coats. In this study, we show that CSLA2 is involved in the biosynthesis of mucilage glucomannan. Furthermore, we show that CSLA2 functions in the maintenance of the normal structure of the adherent mucilage layer through modifying the mucilage cellulose ultrastructure.     

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

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

Arabidopsis 3-Ketoacyl-Coenzyme A Synthase9 Is Involved in the Synthesis of Tetracosanoic Acids as Precursors of Cuticular Waxes,Suberins, Sphingolipids,and Phospholipids

Very-long-chain fatty acids (VLCFAs) with chain lengths from 20 to 34 carbons are involved in diverse biological functions such as membrane constituents, a surface barrier, and seed storage compounds. The first step in VLCFA biosynthesis is the condensation of two carbons to an acyl-coenzyme A, which is catalyzed by 3-ketoacyl-coenzyme A synthase (KCS). In this study, amino acid sequence homology and the messenger RNA expression patterns of 21 Arabidopsis (Arabidopsis thaliana) KCSs were compared. The in planta role of the KCS9 gene, showing higher expression in stem epidermal peels than in stems, was further investigated. The KCS9 gene was ubiquitously expressed in various organs and tissues, including roots, leaves, and stems, including epidermis, silique walls, sepals, the upper portion of the styles, and seed coats, but not in developing embryos. The fluorescent signals of the KCS9::enhanced yellow fluorescent protein construct were merged with those of BrFAD2::monomeric red fluorescent protein, which is an endoplasmic reticulum marker in tobacco (Nicotiana benthamiana) epidermal cells. The kcs9 knockout mutants exhibited a significant reduction in C24 VLCFAs but an accumulation of C20 and C22 VLCFAs in the analysis of membrane and surface lipids. The mutant phenotypes were rescued by the expression of KCS9 under the control of the cauliflower mosaic virus 35S promoter. Taken together, these data demonstrate that KCS9 is involved in the elongation of C22 to C24 fatty acids, which are essential precursors for the biosynthesis of cuticular waxes, aliphatic suberins, and membrane lipids, including sphingolipids and phospholipids. Finally, possible roles of unidentified KCSs are discussed by combining genetic study results and gene expression data from multiple Arabidopsis KCSs.Very-long-chain fatty acids (VLCFAs) are fatty acids of 20 or more carbons in length and are essential precursors of functionally diverse lipids, cuticular waxes, aliphatic suberins, phospholipids, sphingolipids, and seed oils in the Brassicaceae. These lipids are involved in various functions, such as acting as protective barriers between plants and the environment, impermeable barriers to water and ions, energy-storage compounds in seeds, structural components of membranes, and lipid signaling, which is involved in the hypersensitive response (Pollard et al., 2008; Kunst and Samuels, 2009; Franke et al., 2012). VLCFAs are synthesized by the microsomal fatty acid elongase complex, which catalyzes the cyclic addition of a C2 moiety obtained from malonyl-CoA to C16 or C18 acyl-CoA. The fatty acid elongation process has been shown to proceed through a series of four reactions: condensation of the C2 carbon moiety to acyl-CoA by 3-ketoacyl coenzyme A synthase (KCS), reduction of KCS by 3-ketoacyl coenzyme A reductase (KCR), dehydration of 3-hydroxyacyl-CoA by 3-hydroxyacyl-CoA dehydratase (PAS2), and reduction of trans-2,3-enoyl-CoA by trans-2-enoyl-CoA reductase (ECR). Except for KCS isoforms with redundancy, disruption of KCR1, ECR/ECERIFERUM10 (CER10), or PAS2 exhibited severe morphological abnormalities and embryo lethality, suggesting that VLCFA homeostasis is essential for plant developmental processes (Zheng et al., 2005; Bach et al., 2008; Beaudoin et al., 2009).Cuticular waxes that cover plant aerial surfaces are known to be involved in limiting nonstomatal water loss and gaseous exchanges (Boyer et al., 1997; Riederer and Schreiber, 2001), repelling lipophilic pathogenic spores and dust (Barthlott and Neinhuis, 1997), and protecting plants from UV light (Reicosky and Hanover, 1978). VLCFAs that are synthesized in the epidermal cells are either directly used or further modified into aldehydes, alkanes, secondary alcohols, ketones, primary alcohols, and wax esters for the synthesis of cuticular waxes. Reverse genetic analysis and Arabidopsis (Arabidopsis thaliana) epidermal peel microarray analysis (Suh et al., 2005) has enabled the research community to identify the functions of many genes involved in cuticular wax biosynthesis (Kunst and Samuels, 2009): CER1 (Bourdenx et al., 2011; Bernard et al., 2012), WAX2/CER3 (Chen et al., 2003; Rowland et al., 2007; Bernard et al., 2012), and MAH1(Greer et al., 2007; Wen and Jetter, 2009) have been shown to be involved in the decarbonylation pathway to form aldehydes, alkanes, secondary alcohols, and ketones, and acyl-coenzyme A reductase (FAR; Aarts et al., 1997; Rowland et al., 2006) and WSD1 (Li et al., 2008) have been shown to be involved in the decarboxylation pathway for the synthesis of primary alcohols and wax esters. The export of wax precursors to the extracellular space is mediated by a heterodimer of the ATP-binding cassette transporters in the plasma membrane (Pighin et al., 2004; Bird et al., 2007; McFarlane et al., 2010). In addition, glycosylphosphatidylinositol-anchored LTP (LTPG1) and LTPG2 contribute either directly or indirectly to the export of cuticular wax (DeBono et al., 2009; Lee et al., 2009; Kim et al., 2012).VLCFAs that are synthesized in the endodermis of primary roots, seed coats, and the chalaza-micropyle region of seeds are used as precursors for the synthesis of aliphatic suberins. The suberin layer is known to function as a barrier against uncontrolled water, gas, and ion loss and provides protection from environmental stresses and pathogens (Pollard et al., 2008; Franke et al., 2012). For aliphatic suberin biosynthesis, the ω-carbon of the VLCFAs is oxidized by the fatty acyl ω-hydroxylase (Xiao et al., 2004; Li et al., 2007; Höfer et al., 2008; Molina et al., 2008, 2009; Compagnon et al., 2009; Li-Beisson et al., 2009), and the ω-hydroxy VLCFAs are further oxidized into α,ω-dicarboxylic acids by the HOTHEAD-like oxidoreductase (Kurdyukov et al., 2006). α,ω-Dicarboxylic acids are acylated to glycerol-3-P via acyl-CoA:glycerol-3-P acyltransferase (Beisson et al., 2007; Li et al., 2007; Li-Beisson et al., 2009; Yang et al., 2010) or to ferulic acid. In addition, C18, C20, and C22 fatty acids are also reduced by FAR enzymes to primary fatty alcohols, which are a common component in root suberin (Vioque and Kolattukudy, 1997). Finally, the aliphatic suberin precursors are likely to be extensively polymerized and cross linked with the polysaccharides or lignins in the cell wall.In addition, VLCFAs are found in sphingolipids, including glycosyl inositolphosphoceramides, glycosylceramides, and ceramides and phospholipids, such as phosphatidylethanolamine (PE) and phosphatidyl-Ser (PS), which are present in the extraplastidial membrane (Pata et al., 2010; Yamaoka et al., 2011). For sphingolipid biosynthesis, VLCFA-CoAs and Ser are condensed to form 3-keto-sphinganine, which is subsequently reduced to produce sphinganine, a long chain base (LCB). LCBs are known to be further modified by 4-hydroxylation, 4-desaturation, and 8-desaturation (Lynch and Dunn, 2004; Chen et al., 2006, 2012; Pata et al., 2010). The additional VLCFAs are linked with 4-hydroxy LCBs via an amino group to form ceramides (Chen et al., 2008). The presence of VLCFA in sphingolipids may contribute to an increase of their hydrophobicity, membrane leaflet interdigitation, and the transition from a fluid to a gel phase, which is required for microdomain formation. In plants, PS is synthesized from CDP-diacylglycerol and Ser by PS synthase or through an exchange reaction between a phospholipid head group and Ser by a calcium-dependent base-exchange-type PS synthase (Vincent et al., 1999; Yamaoka et al., 2011). PE biosynthesis proceeds through decarboxylation via PS decarboxylase (Nerlich et al., 2007), the phosphoethanolamine transfer from CDP-ethanolamine to diacylglycerol (Kennedy pathway), and the exchange of the head group of PE with Ser via a base-exchange enzyme (Marshall and Kates, 1973). In particular, PS containing a relatively large amount of VLCFAs is enriched in endoplasmic reticulum (ER)-derived vesicles that may function in stabilizing small (70- to 80-nm-diameter) vesicles (Vincent et al., 2001).During the fatty acid elongation process, the first committed step is the condensation of C₂ units to acyl-CoA by KCS. Arabidopsis harbors a large family containing 21 KCS members (Joubès et al., 2008). Characterization of Arabidopsis KCS mutants with defects in VLCFA synthesis revealed in planta roles and substrate specificities (based on differences in carbon chain length and degree of unsaturation) of KCSs. For example, FAE1, a seed-specific condensing enzyme, was shown to catalyze C20 and C22 VLCFA biosynthesis for seed storage lipids (James et al., 1995). KCS6/CER6/CUT1 and KCS5/CER60 are involved in the elongation of fatty acyl-CoAs longer than C28 VLCFA for cuticular waxes in epidermis and pollen coat lipids (Millar et al., 1999; Fiebig et al., 2000; Hooker et al., 2002). KCS20 and KCS2/DAISY are functionally redundant in the two-carbon elongation to C22 VLCFA, which is required for cuticular wax and root suberin biosynthesis (Franke et al., 2009; Lee et al., 2009). When KCS1 and KCS9 were expressed in yeast (Saccharomyces cerevisiae), KCS1 showed broad substrate specificity for saturated and monounsaturated C16 to C24 acyl-CoAs and KCS9 utilized the C16 to C22 acyl-CoAs (Trenkamp et al., 2004; Blacklock and Jaworski, 2006; Paul et al., 2006). Recently, CER2 encoding putative BAHD acyltransferase was reported to be a fatty acid elongase that was involved in the elongation of C28 fatty acids for the synthesis of wax precursors (Haslam et al., 2012).In this study, the expression patterns and subcellular localization of KCS9 were examined, and an Arabidopsis kcs9 mutant was isolated to investigate the roles of KCS9 in planta. Diverse classes of lipids, including cuticular waxes, aliphatic suberins, and sphingolipids, as well as fatty acids in various organs were analyzed from the wild type, the kcs9 mutant, and complementation lines. The combined results of this study revealed that KCS9 is involved in the elongation of C22 to C24 fatty acids, which are essential precursors for the biosynthesis of cuticular waxes, aliphatic suberins, and membrane lipids, including sphingolipids. To the best of our knowledge, this is the first study where a KCS9 isoform involved in sphingolipid biosynthesis was identified.     

Focus Issue on Calcium Signaling: Identification of Cyclic GMP-Activated Nonselective Ca2+-Permeable Cation Channels and Associated CNGC5 and CNGC6 Genes in Arabidopsis Guard Cells

Cytosolic Ca²⁺ in guard cells plays an important role in stomatal movement responses to environmental stimuli. These cytosolic Ca²⁺ increases result from Ca²⁺ influx through Ca²⁺-permeable channels in the plasma membrane and Ca²⁺ release from intracellular organelles in guard cells. However, the genes encoding defined plasma membrane Ca²⁺-permeable channel activity remain unknown in guard cells and, with some exceptions, largely unknown in higher plant cells. Here, we report the identification of two Arabidopsis (Arabidopsis thaliana) cation channel genes, CNGC5 and CNGC6, that are highly expressed in guard cells. Cytosolic application of cyclic GMP (cGMP) and extracellularly applied membrane-permeable 8-Bromoguanosine 3′,5′-cyclic monophosphate-cGMP both activated hyperpolarization-induced inward-conducting currents in wild-type guard cells using Mg²⁺ as the main charge carrier. The cGMP-activated currents were strongly blocked by lanthanum and gadolinium and also conducted Ba²⁺, Ca²⁺, and Na⁺ ions. cngc5 cngc6 double mutant guard cells exhibited dramatically impaired cGMP-activated currents. In contrast, mutations in CNGC1, CNGC2, and CNGC20 did not disrupt these cGMP-activated currents. The yellow fluorescent protein-CNGC5 and yellow fluorescent protein-CNGC6 proteins localize in the cell periphery. Cyclic AMP activated modest inward currents in both wild-type and cngc5cngc6 mutant guard cells. Moreover, cngc5 cngc6 double mutant guard cells exhibited functional abscisic acid (ABA)-activated hyperpolarization-dependent Ca²⁺-permeable cation channel currents, intact ABA-induced stomatal closing responses, and whole-plant stomatal conductance responses to darkness and changes in CO₂ concentration. Furthermore, cGMP-activated currents remained intact in the growth controlled by abscisic acid2 and abscisic acid insensitive1 mutants. This research demonstrates that the CNGC5 and CNGC6 genes encode unique cGMP-activated nonselective Ca²⁺-permeable cation channels in the plasma membrane of Arabidopsis guard cells.Plants lose water via transpiration and take in CO₂ for photosynthesis through stomatal pores. Each stomatal pore is surrounded by two guard cells, and stomatal movements are driven by the change of turgor pressure in guard cells. The intracellular second messenger Ca²⁺ functions in guard cell signal transduction (Schroeder and Hagiwara, 1989; McAinsh et al., 1990; Webb et al., 1996; Grabov and Blatt, 1998; Allen et al., 1999; MacRobbie, 2000; Mori et al., 2006; Young et al., 2006; Siegel et al., 2009; Chen et al., 2010; Hubbard et al., 2012). Plasma membrane ion channel activity and gene expression in guard cells are finely regulated by the intracellular free calcium concentration ([Ca2+]_cyt; Schroeder and Hagiwara, 1989; Webb et al., 2001; Allen et al., 2002; Siegel et al., 2009; Kim et al., 2010; Stange et al., 2010). Ca²⁺-dependent protein kinases (CPKs) function as targets of the cytosolic Ca²⁺ signal, and several members of the CPK family have been shown to function in stimulus-induced stomatal closing, including the Arabidopsis (Arabidopsis thaliana) CPK3, CPK4, CPK6, CPK10, and CPK11 proteins (Mori et al., 2006; Zhu et al., 2007; Zou et al., 2010; Brandt et al., 2012; Hubbard et al., 2012). Further research found that several CPKs could activate the S-type anion channel SLAC1 in Xenopus laevis oocytes, including CPK21, CPK23, and CPK6 (Geiger et al., 2010; Brandt et al., 2012). At the same time, the Ca²⁺-independent protein kinase Open Stomata1 mediates stomatal closing and activates the S-type anion channel SLAC1 (Mustilli et al., 2002; Yoshida et al., 2002; Geiger et al., 2009; Lee et al., 2009; Xue et al., 2011), indicating that both Ca²⁺-dependent and Ca²⁺-independent pathways function in guard cells.Multiple essential factors of guard cell abscisic acid (ABA) signal transduction function in the regulation of Ca²⁺-permeable channels and [Ca2+]_cyt elevations, including Abscisic Acid Insensitive1 (ABI1), ABI2, Enhanced Response to Abscisic Acid1 (ERA1), the NADPH oxidases AtrbohD and AtrbohF, the Guard Cell Hydrogen Peroxide-Resistant1 (GHR1) receptor kinase, as well as the Ca²⁺-activated CPK6 protein kinase (Pei et al., 1998; Allen et al., 1999, 2002; Kwak et al., 2003; Miao et al., 2006; Mori et al., 2006; Hua et al., 2012). [Ca2+]_cyt increases result from both Ca²⁺ release from intracellular Ca²⁺ stores (McAinsh et al., 1992) and Ca²⁺ influx across the plasma membrane (Hamilton et al., 2000; Pei et al., 2000; Murata et al., 2001; Kwak et al., 2003; Hua et al., 2012). Electrophysiological analyses have characterized nonselective Ca²⁺-permeable channel activity in the plasma membrane of guard cells (Schroeder and Hagiwara, 1990; Hamilton et al., 2000; Pei et al., 2000; Murata et al., 2001; Köhler and Blatt, 2002; Miao et al., 2006; Mori et al., 2006; Suh et al., 2007; Vahisalu et al., 2008; Hua et al., 2012). However, the genetic identities of Ca²⁺-permeable channels in the plasma membrane of guard cells have remained unknown despite over two decades of research on these channel activities.The Arabidopsis genome includes 20 genes encoding cyclic nucleotide-gated channel (CNGC) homologs and 20 genes encoding homologs to animal Glu receptor channels (Lacombe et al., 2001; Kaplan et al., 2007; Ward et al., 2009), which have been proposed to function in plant cells as cation channels (Schuurink et al., 1998; Arazi et al., 1999; Köhler et al., 1999). Recent research has demonstrated functions of specific Glu receptor channels in mediating Ca²⁺ channel activity (Michard et al., 2011; Vincill et al., 2012). Previous studies have shown cAMP activation of nonselective cation currents in guard cells (Lemtiri-Chlieh and Berkowitz, 2004; Ali et al., 2007). However, only a few studies have shown the disappearance of a defined plasma membrane Ca²⁺ channel activity in plants upon mutation of candidate Ca²⁺ channel genes (Ali et al., 2007; Michard et al., 2011; Laohavisit et al., 2012; Vincill et al., 2012). Some CNGCs have been found to be involved in cation nutrient intake, including monovalent cation intake (Guo et al., 2010; Caballero et al., 2012), salt tolerance (Guo et al., 2008; Kugler et al., 2009), programmed cell death and pathogen responses (Clough et al., 2000; Balagué et al., 2003; Urquhart et al., 2007; Abdel-Hamid et al., 2013), thermal sensing (Finka et al., 2012; Gao et al., 2012), and pollen tube growth (Chang et al., 2007; Frietsch et al., 2007; Tunc-Ozdemir et al., 2013a, 2013b). Direct in vivo disappearance of Ca²⁺ channel activity in cngc disruption mutants has been demonstrated in only a few cases thus far (Ali et al., 2007; Gao et al., 2012). In this research, we show that CNGC5 and CNGC6 are required for a cyclic GMP (cGMP)-activated nonselective Ca²⁺-permeable cation channel activity in the plasma membrane of Arabidopsis guard cells.     

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.     

Focus on Ethylene: Ethylene Biosynthesis Is Promoted by Very-Long-Chain Fatty Acids during Lysigenous Aerenchyma Formation in Rice Roots

In rice (Oryza sativa) roots, lysigenous aerenchyma, which is created by programmed cell death and lysis of cortical cells, is constitutively formed under aerobic conditions, and its formation is further induced under oxygen-deficient conditions. Ethylene is involved in the induction of aerenchyma formation. reduced culm number1 (rcn1) is a rice mutant in which the gene encoding the ATP-binding cassette transporter RCN1/OsABCG5 is defective. Here, we report that the induction of aerenchyma formation was reduced in roots of rcn1 grown in stagnant deoxygenated nutrient solution (i.e. under stagnant conditions, which mimic oxygen-deficient conditions in waterlogged soils). 1-Aminocyclopropane-1-carboxylic acid synthase (ACS) is a key enzyme in ethylene biosynthesis. Stagnant conditions hardly induced the expression of ACS1 in rcn1 roots, resulting in low ethylene production in the roots. Accumulation of saturated very-long-chain fatty acids (VLCFAs) of 24, 26, and 28 carbons was reduced in rcn1 roots. Exogenously supplied VLCFA (26 carbons) increased the expression level of ACS1 and induced aerenchyma formation in rcn1 roots. Moreover, in rice lines in which the gene encoding a fatty acid elongase, CUT1-LIKE (CUT1L; a homolog of the gene encoding Arabidopsis CUT1, which is required for cuticular wax production), was silenced, both ACS1 expression and aerenchyma formation were reduced. Interestingly, the expression of ACS1, CUT1L, and RCN1/OsABCG5 was induced predominantly in the outer part of roots under stagnant conditions. These results suggest that, in rice under oxygen-deficient conditions, VLCFAs increase ethylene production by promoting 1-aminocyclopropane-1-carboxylic acid biosynthesis in the outer part of roots, which, in turn, induces aerenchyma formation in the root cortex.Aerenchyma formation is a morphological adaptation of plants to complete submergence and waterlogging of the soil, and facilitates internal gas diffusion (Armstrong, 1979; Jackson and Armstrong, 1999; Colmer, 2003; Voesenek et al., 2006; Bailey-Serres and Voesenek, 2008; Licausi and Perata, 2009; Sauter, 2013; Voesenek and Bailey-Serres, 2015). To adapt to waterlogging in soil, rice (Oryza sativa) develops lysigenous aerenchyma in shoots (Matsukura et al., 2000; Colmer and Pedersen, 2008; Steffens et al., 2011) and roots (Jackson et al., 1985b; Justin and Armstrong, 1991; Kawai et al., 1998), which is formed by programmed cell death and subsequent lysis of some cortical cells (Jackson and Armstrong, 1999; Evans, 2004; Yamauchi et al., 2013). In rice roots, lysigenous aerenchyma is constitutively formed under aerobic conditions (Jackson et al., 1985b), and its formation is further induced under oxygen-deficient conditions (Colmer et al., 2006; Shiono et al., 2011). The former and latter are designated constitutive and inducible lysigenous aerenchyma formation, respectively (Colmer and Voesenek, 2009). The gaseous plant hormone ethylene regulates adaptive growth responses of plants to submergence (Voesenek and Blom, 1989; Voesenek et al., 1993; Visser et al., 1996a,b; Lorbiecke and Sauter, 1999; Hattori et al., 2009; Steffens and Sauter, 2009; van Veen et al., 2013). Ethylene also induces lysigenous aerenchyma formation in roots of some gramineous plants (Drew et al., 2000; Shiono et al., 2008). The treatment of roots with ethylene or its precursor (1-aminocyclopropane-1-carboxylic acid [ACC]) stimulates aerenchyma formation in rice (Justin and Armstrong, 1991; Colmer et al., 2006; Yukiyoshi and Karahara, 2014), maize (Zea mays; Drew et al., 1981; Jackson et al., 1985a; Takahashi et al., 2015), and wheat (Triticum aestivum; Yamauchi et al., 2014a,b). Moreover, treatment of roots with inhibitors of ethylene action or ethylene biosynthesis effectively blocks aerenchyma formation under hypoxic conditions in maize (Drew et al., 1981; Konings, 1982; Jackson et al., 1985a; Rajhi et al., 2011).Ethylene biosynthesis is accomplished by two main successive enzymatic reactions: conversion of S-adenosyl-Met to ACC by 1-aminocyclopropane-1-carboxylic acid synthase (ACS), and conversion of ACC to ethylene by 1-aminocyclopropane-1-carboxylic acid oxidase (ACO; Yang and Hoffman, 1984). The activities of both enzymes are enhanced during aerenchyma formation under hypoxic conditions in maize root (He et al., 1996). Since the ACC content in roots of maize is increased by oxygen deficiency and is strongly correlated with ethylene production (Atwell et al., 1988), ACC biosynthesis is essential for ethylene production during aerenchyma formation in roots. In fact, exogenously supplied ACC induced ethylene production in roots of maize (Drew et al., 1979; Konings, 1982; Atwell et al., 1988) and wheat (Yamauchi et al., 2014b), even under aerobic conditions. Ethylene production in plants is inversely related to oxygen concentration (Yang and Hoffman, 1984). Under anoxic conditions, the oxidation of ACC to ethylene by ACO, which requires oxygen, is almost completely repressed (Yip et al., 1988; Tonutti and Ramina, 1991). Indeed, anoxic conditions stimulate neither ethylene production nor aerenchyma formation in maize adventitious roots (Drew et al., 1979). Therefore, it is unlikely that the root tissues forming inducible aerenchyma are anoxic, and that the ACO-mediated step is repressed. Moreover, aerenchyma is constitutively formed in rice roots even under aerobic conditions (Jackson et al., 1985b), and thus, after the onset of waterlogging, oxygen can be immediately supplied to the apical regions of roots through the constitutively formed aerenchyma.Very-long-chain fatty acids (VLCFAs; ≥20 carbons) are major constituents of sphingolipids, cuticular waxes, and suberin in plants (Franke and Schreiber, 2007; Kunst and Samuels, 2009). In addition to their structural functions, VLCFAs directly or indirectly participate in several physiological processes (Zheng et al., 2005; Reina-Pinto et al., 2009; Roudier et al., 2010; Ito et al., 2011; Nobusawa et al., 2013; Tsuda et al., 2013), including the regulation of ethylene biosynthesis (Qin et al., 2007). During fiber cell elongation in cotton ovules, ethylene biosynthesis is enhanced by treatment with saturated VLCFAs, especially 24-carbon fatty acids, and is suppressed by an inhibitor of VLCFA biosynthesis (Qin et al., 2007). The first rate-limiting step in VLCFA biosynthesis is condensation of acyl-CoA with malonyl-CoA by β-ketoacyl-CoA synthase (KCS; Joubès et al., 2008). KCS enzymes are thought to determine the substrate and tissue specificities of fatty acid elongation (Joubès et al., 2008). The Arabidopsis (Arabidopsis thaliana) genome has 21 KCS genes (Joubès et al., 2008). In the Arabidopsis cut1 mutant, which has a defect in the gene encoding CUT1 that is required for cuticular wax production (i.e. one of the KCS genes), the expression of AtACO genes and growth of root cells were reduced when compared with the wild type (Qin et al., 2007). Furthermore, expression of the AtACO genes was rescued by exogenously supplied saturated VLCFAs (Qin et al., 2007). These observations imply that VLCFAs or their derivatives work as regulatory factors for gene expression during some physiological processes in plants.reduced culm number1 (rcn1) was first identified as a rice mutant with a low tillering rate in a paddy field (Takamure and Kinoshita, 1985; Yasuno et al., 2007). The rcn1 (rcn1-2) mutant has a single nucleotide substitution in the gene encoding a member of the ATP-binding cassette (ABC) transporter subfamily G, RCN1/OsABCG5, causing an Ala-684Pro substitution (Yasuno et al., 2009). The mutation results in several mutant phenotypes, although the substrates of RCN1/OsABCG5 have not been determined (Ureshi et al., 2012; Funabiki et al., 2013; Matsuda et al., 2014). We previously found that the rcn1 mutant has abnormal root morphology, such as shorter root length and brownish appearance of roots, under stagnant (deoxygenated) conditions (which mimics oxygen-deficient conditions in waterlogged soils). We also found that the rcn1 mutant accumulates less of the major suberin monomers originating from VLCFAs in the outer part of adventitious roots, and this results in a reduction of a functional apoplastic barrier in the root hypodermis (Shiono et al., 2014a).The objective of this study was to elucidate the molecular basis of inducible aerenchyma formation. To this end, we examined lysigenous aerenchyma formation and ACC, ethylene, and VLCFA accumulation and their biosyntheses in rcn1 roots. Based on the results of these studies, we propose that VLCFAs are involved in inducible aerenchyma formation through the enhancement of ethylene biosynthesis in rice roots.