A Dibasic Amino Acid Pair Conserved in the Activation Loop Directs Plasma Membrane Localization and Is Necessary for Activity of Plant Type I/II Phosphatidylinositol Phosphate Kinase

Phosphatidylinositol phosphate kinase (PIPK) is an enzyme involved in the regulation of cellular levels of phosphoinositides involved in various physiological processes, such as cytoskeletal organization, ion channel activation, and vesicle trafficking. In animals, research has focused on the modes of activation and function of PIPKs, providing an understanding of the importance of plasma membrane localization. However, it still remains unclear how this issue is regulated in plant PIPKs. Here, we demonstrate that the carboxyl-terminal catalytic domain, which contains the activation loop, is sufficient for plasma membrane localization of PpPIPK1, a type I/II B PIPK from the moss Physcomitrella patens. The importance of the carboxyl-terminal catalytic domain for plasma membrane localization was confirmed with Arabidopsis (Arabidopsis thaliana) AtPIP5K1. Our findings, in which substitution of a conserved dibasic amino acid pair in the activation loop of PpPIPK1 completely prevented plasma membrane targeting and abolished enzymatic activity, demonstrate its critical role in these processes. Placing our results in the context of studies of eukaryotic PIPKs led us to conclude that the function of the dibasic amino acid pair in the activation loop in type I/II PIPKs is plant specific.Phosphoinositides (PIs) are minor lipids found in membrane fractions but implicated in a wide variety of physiological regulations in eukaryotes (Di Paolo and De Camilli, 2006; Zonia and Munnik, 2006). Phosphatidylinositol-4,5-bisphosphate [PtdIns(4,5)P₂] is a major PI in animal plasma membranes, affecting the localization and activity of various kinds of proteins carrying phosphatidylinositol-binding domains, which in turn affect the regulation of cytoskeletal organization, vesicle trafficking, cell proliferation, and cell growth during development and stress responses (Doughman et al., 2003; Downes et al., 2005; Di Paolo and De Camilli, 2006; Zonia and Munnik, 2006; Heck et al., 2007). In addition, PtdIns(4,5)P₂ is also a well-known substrate of phospholipase C, producing second messengers such as diacylglycerol, phosphatidic acid (PA), and inositol-1,4,5-trisphosphate, which are involved in the activation of intracellular signal transduction pathways (Zonia and Munnik, 2006). Transient accumulation of PtdIns(4,5)P₂ has also been observed under various kinds of environmental stress (Pical et al., 1999; DeWald et al., 2001), suggesting an important role of this lipid in the regulation of stress signal transduction pathways also in plants. These findings indicate that PtdIns(4,5)P₂ is multifunctional and involved in a variety of cellular processes. Therefore, elucidation of the mechanisms controlling the cellular levels of PtdIns(4,5)P₂ is important in understanding the significance of PI signaling in eukaryotes.PtdIns(4,5)P₂ is synthesized by phosphatidylinositol phosphate kinases (PIPKs; Anderson et al., 1999; Doughman et al., 2003; Heck et al., 2007). Physiological roles of several plant PIPKs have been reported. In Arabidopsis (Arabidopsis thaliana), AtPIP5K3 is an essential regulator of tip growth of root hairs (Kusano et al., 2008; Stenzel et al., 2008), while AtPIPK4 and AtPIPK5 are essential for pollen germination and pollen tube elongation (Ischebeck et al., 2008; Sousa et al., 2008). In addition, AtPIP5K9 was shown to interact with the cytosolic invertase CINV1 to regulate sugar-mediated root cell elongation negatively (Lou et al., 2007). Rice (Oryza sativa) OsPIPK1 is proposed to be involved in shoot growth and floral initiation through the regulation of floral induction genes (Ma et al., 2004). In animals, membrane-associated type I PIPK mainly phosphorylates the D-5 hydroxyl group of PtdIns4P to produce PtdIns(4,5)P₂ but also produces PtdIns(3,4)P₂ and PtdIns(3,5)P₂ from PtdIns3P with 5- and 4-kinase activity (Anderson et al., 1999; Heck et al., 2007), whereas type II PIPK prefers the D-4 position of PtdIns5P, producing PtdIns(4,5)P₂ in the nucleus and at the endoplasmic reticulum (Clarke et al., 2007). Thus, in animals, type I and II PIPKs are involved in the generation of PtdIns(4,5)P₂ via different pathways. Molecular biological analysis of plant PIPKs was initiated with AtPIP5K1 from Arabidopsis (Mikami et al., 1998), which phosphorylates PtdIns3P, PtdIns4P, and PtdIns(4,5)P₂ to produce PtdIns(3,4)P₂, PtdIns(4,5)P₂, and PtdIns (3,4,5)P₃, respectively, with D-4- and D-5-kinase activity (Elge et al., 2001; Westergren et al., 2001; Im et al., 2007). Similar enzymatic activity was also reported for other PIPKs from Arabidopsis (Ischebeck et al., 2008; Kusano et al., 2008; Stenzel et al., 2008). In addition, a PIPK from the moss Physcomitrella patens (designated as PpPIPK1) preferred PtdIns4P, PtdIns3P, and PtdIns(3,4)P₂ as substrates, but not PtdIns5P, producing PtdIns(4,5)P₂, PtdIns(3,4)P₂, and PtdIns(3,4,5)P₃, respectively (Saavedra et al., 2009). These findings indicate that the substrate specificity of plant PIPKs is essentially the same as that of type I PIPKs. However, AtPIP5K1 has yet to be classified as either type I or type II based on sequence comparisons of the catalytic domain (CD; Mikami et al., 1998). This was confirmed by a genome-wide analysis of PIPK genes in Arabidopsis in which all 11 PIPKs were classified as type I/II based on sequence comparisons of the CDs, which were further subdivided into subtypes A and B (Mueller-Roeber and Pical, 2002). Therefore, it is suggested that typical type I and II PIPKs are absent in plants, although further confirmation is needed.The conserved PIPK CD contains a short highly conserved region near its C-terminal end, designated the activation loop, which acts as the substrate-binding site and is responsible for the differences in substrate specificity and subcellular localization between animal type I and type II PIPKs (Kunz et al., 2000, 2002). Substrate specificities of animal type I and type II PIPKs, for example, are determined by a respective Glu and Ala at the corresponding positions in the activation loop. Moreover, it has been established that substitution of Glu to Ala results in a swap of substrate specificity and subcellular localization between the two types (Kunz et al., 2000, 2002). In contrast to animal PIPKs, a substitution in the activation loop of PpPIPK1 from Glu to Ala resulted in a nearly complete loss of type I/II activity; however, such a mutation did not fully convert the substrate specificity, although an enhancement of type II versus type I activity was observed (Saavedra et al., 2009). Since the corresponding amino acid residue is Glu in all plant PIPKs so far reported, it is suggested that there also is a plant-specific mode of substrate specificity regulation in plant type I/II PIPKs. However, enzymatic activity appears to be modified in similar ways between plant type I/II and animal type I PIPKs; that is, phosphorylation- and PA-dependent activation of PIPKs has been observed in both animals and plants (Moritz et al., 1992; Jenkins et al., 1994; Pical et al., 1999; Westergren et al., 2001; Perera et al., 2005; Saavedra et al., 2009).The regulation of plasma membrane localization of mammalian type I PIPKs remains confusing. In addition to the involvement of a Glu residue as mentioned above, the substitution of two Lys residues in the activation loop to Asn residues changes the subcellular localization from the plasma membrane to the cytosol (Kunz et al., 2000, 2002). However, Arioka et al. (2004) also showed that the plasma membrane localization of type I PIPKs is regulated by another basic amino acid pair localized downstream of the activation loop in the CD, which is not found in type II PIPKs. Interestingly, the mechanism behind plasma membrane localization of plant PIPKs seems to differ significantly from the animal one. The obvious structural feature of plant PIPKs is the presence of a repetition of membrane occupation and recognition nexus (MORN) motifs at the N-terminal half, which is conserved across the B subfamily of plant type I/II PIPKs (Mueller-Roeber and Pical, 2002). The MORN motif was first identified in mammalian junctophilin, an endoplasmic reticulum-membrane-bound component of the junctional complex between the plasma membrane and the endoplasmic reticulum (Takeshima et al., 2000). Since MORN motifs are not found in PIPKs from nonplant organisms, a plant-specific mode of PIPK activation is speculated. Indeed, a regulatory role of the MORN domain was reported in the enzymatic activation of AtPIP5K1 (Im et al., 2007) and in root hair formation, but not in enzymatic activation, of AtPIP5K3 (Stenzel et al., 2008). Moreover, the MORN domain may play a role in the plasma membrane localization of OsPIPK1 from rice and AtPIP5K1 and AtPIP5K3 from Arabidopsis (Ma et al., 2006; Im et al., 2007; Kusano et al., 2008). However, stable transformation of tobacco (Nicotiana tabacum) cells to express an AtPIP5K1 MORN domain-GFP fusion did not allow visualization of the plasma membrane localization of this protein (Im et al., 2007). Thus, it is not clear if the MORN domain functions as a plasma membrane-targeting module.Given the sequence conservation of the CD among eukaryotic PIPKs (Saavedra et al., 2009), we hypothesize that the CD is responsible for the plasma membrane localization of plant PIPKs. Thus, to gain further insight into the mechanisms regulating this issue, we dissected PpPIPK1 to determine the molecular determinants of plasma membrane localization. Here, we show that the MORN domain is not involved in the plasma membrane localization of PpPIPK1 and AtPIP5K1 in P. patens protoplasts and onion (Allium cepa) epidermal cells. We further demonstrate that two basic amino acids, but not Glu, conserved in the activation loop of the CD are required for plasma membrane localization. These findings demonstrate that the activation mode of type I/II PIPKs is plant specific and differs from that of the membrane-localized animal type I PIPKs.     

Phosphorylation of the C Terminus of RHD3 Has a Critical Role in Homotypic ER Membrane Fusion in Arabidopsis

The endoplasmic reticulum (ER) consists of dynamically changing tubules and cisternae. In animals and yeast, homotypic ER membrane fusion is mediated by fusogens (atlastin and Sey1p, respectively) that are membrane-associated dynamin-like GTPases. In Arabidopsis (Arabidopsis thaliana), another dynamin-like GTPase, ROOT HAIR DEFECTIVE3 (RHD3), has been proposed as an ER membrane fusogen, but direct evidence is lacking. Here, we show that RHD3 has an ER membrane fusion activity that is enhanced by phosphorylation of its C terminus. The ER network was RHD3-dependently reconstituted from the cytosol and microsome fraction of tobacco (Nicotiana tabacum) cultured cells by exogenously adding GTP, ATP, and F-actin. We next established an in vitro assay system of ER tubule formation with Arabidopsis ER vesicles, in which addition of GTP caused ER sac formation from the ER vesicles. Subsequent application of a shearing force to this system triggered the formation of tubules from the ER sacs in an RHD-dependent manner. Unexpectedly, in the absence of a shearing force, Ser/Thr kinase treatment triggered RHD3-dependent tubule formation. Mass spectrometry showed that RHD3 was phosphorylated at multiple Ser and Thr residues in the C terminus. An antibody against the RHD3 C-terminal peptide abolished kinase-triggered tubule formation. When the Ser cluster was deleted or when the Ser residues were replaced with Ala residues, kinase treatment had no effect on tubule formation. Kinase treatment induced the oligomerization of RHD3. Neither phosphorylation-dependent modulation of membrane fusion nor oligomerization has been reported for atlastin or Sey1p. Taken together, we propose that phosphorylation-stimulated oligomerization of RHD3 enhances ER membrane fusion to form the ER network.In eukaryotic cells, the endoplasmic reticulum (ER) is the organelle with the largest membrane area. The ER consists of an elaborate network of interconnected membrane tubules and cisternae that is continually moving and being remodeled (Friedman and Voeltz, 2011). In plant cells, ER movement and remodeling is primarily driven by the actin-myosin XI cytoskeleton (Sparkes et al., 2009; Ueda et al., 2010; Yokota et al., 2011; Griffing et al., 2014) and secondarily by the microtubule cytoskeleton (Hamada et al., 2014). Several factors involved in creating the ER architecture have been also identified (Anwar et al., 2012; Chen et al., 2012; Goyal and Blackstone, 2013; Sackmann, 2014; Stefano et al., 2014a; Westrate et al., 2015). Among them, ER membrane-bound GTPases, animal atlastins and yeast Sey1p (Synthetic Enhancement of Yop1), function as ER fusogens to form the interconnected tubular network (Hu et al., 2009; Orso et al., 2009; Anwar et al., 2012). Atlastin molecules on the two opposed membranes have been proposed to transiently dimerize to attract the two membranes to each other (Bian et al., 2011; Byrnes and Sondermann, 2011; Morin-Leisk et al., 2011; Moss et al., 2011; Lin et al., 2012; Byrnes et al., 2013). Closely attracted lipid bilayers are supposed to be destabilized by an amphipathic helical domain at the atlastin C terminus to facilitate membrane fusion (Bian et al., 2011; Liu et al., 2012; Faust et al., 2015). Knockdown of atlastins leads to fragmentation of the ER and unbranched ER tubules, while overexpression of atlastins enhances ER membrane fusion, which enlarges the ER profiles (Hu et al., 2009; Orso et al., 2009).An Arabidopsis (Arabidopsis thaliana) protein, ROOT HAIR DEFECTIVE3 (RHD3), has been proposed as a fusogen because (1) when it is disrupted, the ER network is modified into large cable-like strands of poorly branched membranes (Zheng et al., 2004; Chen et al., 2011; Stefano et al., 2012), (2) it shares sequence similarity with the above-mentioned fusogen Sey1p (Hu et al., 2009), and (3) it has structural similarity to atlastin and Sey1p, with a functional GTPase domain at the N-terminal cytosolic domain (Stefano et al., 2012) followed by two transmembrane domains and a cytosolic tail. RHD3 has a longer cytosolic C-terminal tail than do atlastin and Sey1p (Stefano and Brandizzi, 2014). It contains not only an amphipathic region but also a Ser/Thr-rich C terminus.Arabidopsis has two RHD3 isoforms called RHD3-Like 1 and RHD3-Like 2. Fluorescently tagged RHD3 and RHD3-Like 2 localize to the ER (Chen et al., 2011; Stefano et al., 2012; Lee et al., 2013). RHD3 and the two RHD3-Like proteins likely have redundant roles in ER membrane fusion (Zhang et al., 2013). Overexpression of either RHD3 or RHD3-Like 2 with a defective GTPase domain phenocopies the aberrant ER morphology in rhd3-deficient mutants (Chen et al., 2011; Lee et al., 2013).In this study, we show that the Ser/Thr-rich C terminus enhances ER membrane fusion following phosphorylation of its C terminus. We propose a model in which phosphorylation and oligomerization of RHD3 is required for efficient ER membrane fusion. Our findings clarify the mechanisms that regulate RHD3 and consequently the homeostasis of membrane fusion in the ER.     

Characterization of the Possible Roles for B Class MADS Box Genes in Regulation of Perianth Formation in Orchid

To investigate sepal/petal/lip formation in Oncidium Gower Ramsey, three paleoAPETALA3 genes, O. Gower Ramsey MADS box gene5 (OMADS5; clade 1), OMADS3 (clade 2), and OMADS9 (clade 3), and one PISTILLATA gene, OMADS8, were characterized. The OMADS8 and OMADS3 mRNAs were expressed in all four floral organs as well as in vegetative leaves. The OMADS9 mRNA was only strongly detected in petals and lips. The mRNA for OMADS5 was only strongly detected in sepals and petals and was significantly down-regulated in lip-like petals and lip-like sepals of peloric mutant flowers. This result revealed a possible negative role for OMADS5 in regulating lip formation. Yeast two-hybrid analysis indicated that OMADS5 formed homodimers and heterodimers with OMADS3 and OMADS9. OMADS8 only formed heterodimers with OMADS3, whereas OMADS3 and OMADS9 formed homodimers and heterodimers with each other. We proposed that sepal/petal/lip formation needs the presence of OMADS3/8 and/or OMADS9. The determination of the final organ identity for the sepal/petal/lip likely depended on the presence or absence of OMADS5. The presence of OMADS5 caused short sepal/petal formation. When OMADS5 was absent, cells could proliferate, resulting in the possible formation of large lips and the conversion of the sepal/petal into lips in peloric mutants. Further analysis indicated that only ectopic expression of OMADS8 but not OMADS5/9 caused the conversion of the sepal into an expanded petal-like structure in transgenic Arabidopsis (Arabidopsis thaliana) plants.The ABCDE model predicts the formation of any flower organ by the interaction of five classes of homeotic genes in plants (Yanofsky et al., 1990; Jack et al., 1992; Mandel et al., 1992; Goto and Meyerowitz, 1994; Jofuku et al., 1994; Pelaz et al., 2000, 2001; Theißen and Saedler, 2001; Pinyopich et al., 2003; Ditta et al., 2004; Jack, 2004). The A class genes control sepal formation. The A, B, and E class genes work together to regulate petal formation. The B, C, and E class genes control stamen formation. The C and E class genes work to regulate carpel formation, whereas the D class gene is involved in ovule development. MADS box genes seem to have a central role in flower development, because most ABCDE genes encode MADS box proteins (Coen and Meyerowitz, 1991; Weigel and Meyerowitz, 1994; Purugganan et al., 1995; Rounsley et al., 1995; Theißen and Saedler, 1995; Theißen et al., 2000; Theißen, 2001).The function of B group genes, such as APETALA3 (AP3) and PISTILLATA (PI), has been thought to have a major role in specifying petal and stamen development (Jack et al., 1992; Goto and Meyerowitz, 1994; Krizek and Meyerowitz, 1996; Kramer et al., 1998; Hernandez-Hernandez et al., 2007; Kanno et al., 2007; Whipple et al., 2007; Irish, 2009). In Arabidopsis (Arabidopsis thaliana), mutation in AP3 or PI caused identical phenotypes of second whorl petal conversion into a sepal structure and third flower whorl stamen into a carpel structure (Bowman et al., 1989; Jack et al., 1992; Goto and Meyerowitz, 1994). Similar homeotic conversions for petal and stamen were observed in the mutants of the AP3 and PI orthologs from a number of core eudicots such as Antirrhinum majus, Petunia hybrida, Gerbera hybrida, Solanum lycopersicum, and Nicotiana benthamiana (Sommer et al., 1990; Tröbner et al., 1992; Angenent et al., 1993; van der Krol et al., 1993; Yu et al., 1999; Liu et al., 2004; Vandenbussche et al., 2004; de Martino et al., 2006), from basal eudicot species such as Papaver somniferum and Aquilegia vulgaris (Drea et al., 2007; Kramer et al., 2007), as well as from monocot species such as Zea mays and Oryza sativa (Ambrose et al., 2000; Nagasawa et al., 2003; Prasad and Vijayraghavan, 2003; Yadav et al., 2007; Yao et al., 2008). This indicated that the function of the B class genes AP3 and PI is highly conserved during evolution.It has been thought that B group genes may have arisen from an ancestral gene through multiple gene duplication events (Doyle, 1994; Theißen et al., 1996, 2000; Purugganan, 1997; Kramer et al., 1998; Kramer and Irish, 1999; Lamb and Irish, 2003; Kim et al., 2004; Stellari et al., 2004; Zahn et al., 2005; Hernandez-Hernandez et al., 2007). In the gymnosperms, there was a single putative B class lineage that duplicated to generate the paleoAP3 and PI lineages in angiosperms (Kramer et al., 1998; Theißen et al., 2000; Irish, 2009). The paleoAP3 lineage is composed of AP3 orthologs identified in lower eudicots, magnolid dicots, and monocots (Kramer et al., 1998). Genes in this lineage contain the conserved paleoAP3- and PI-derived motifs in the C-terminal end of the proteins, which have been thought to be characteristics of the B class ancestral gene (Kramer et al., 1998; Tzeng and Yang, 2001; Hsu and Yang, 2002). The PI lineage is composed of PI orthologs that contain a highly conserved PI motif identified in most plant species (Kramer et al., 1998). Subsequently, there was a second duplication at the base of the core eudicots that produced the euAP3 and TM6 lineages, which have been subject to substantial sequence changes in eudicots during evolution (Kramer et al., 1998; Kramer and Irish, 1999). The paleoAP3 motif in the C-terminal end of the proteins was retained in the TM6 lineage and replaced by a conserved euAP3 motif in the euAP3 lineage of most eudicot species (Kramer et al., 1998). In addition, many lineage-specific duplications for paleoAP3 lineage have occurred in plants such as orchids (Hsu and Yang, 2002; Tsai et al., 2004; Kim et al., 2007; Mondragón-Palomino and Theißen, 2008, 2009; Mondragón-Palomino et al., 2009), Ranunculaceae, and Ranunculales (Kramer et al., 2003; Di Stilio et al., 2005; Shan et al., 2006; Kramer, 2009).Unlike the A or C class MADS box proteins, which form homodimers that regulate flower development, the ability of B class proteins to form homodimers has only been reported in gymnosperms and in the paleoAP3 and PI lineages of some monocots. For example, LMADS1 of the lily Lilium longiflorum (Tzeng and Yang, 2001), OMADS3 of the orchid Oncidium Gower Ramsey (Hsu and Yang, 2002), and PeMADS4 of the orchid Phalaenopsis equestris (Tsai et al., 2004) in the paleoAP3 lineage, LRGLOA and LRGLOB of the lily Lilium regale (Winter et al., 2002), TGGLO of the tulip Tulipa gesneriana (Kanno et al., 2003), and PeMADS6 of the orchid P. equestris (Tsai et al., 2005) in the PI lineage, and GGM2 of the gymnosperm Gnetum gnemon (Winter et al., 1999) were able to form homodimers that regulate flower development. Proteins in the euAP3 lineage and in most paleoAP3 lineages were not able to form homodimers and had to interact with PI to form heterodimers in order to regulate petal and stamen development in various plant species (Schwarz-Sommer et al., 1992; Tröbner et al., 1992; Riechmann et al., 1996; Moon et al., 1999; Winter et al., 2002; Kanno et al., 2003; Vandenbussche et al., 2004; Yao et al., 2008). In addition to forming dimers, AP3 and PI were able to interact with other MADS box proteins, such as SEPALLATA1 (SEP1), SEP2, and SEP3, to regulate petal and stamen development (Pelaz et al., 2000; Honma and Goto, 2001; Theißen and Saedler, 2001; Castillejo et al., 2005).Orchids are among the most important plants in the flower market around the world, and research on MADS box genes has been reported for several species of orchids during the past few years (Lu et al., 1993, 2007; Yu and Goh, 2000; Hsu and Yang, 2002; Yu et al., 2002; Hsu et al., 2003; Tsai et al., 2004, 2008; Xu et al., 2006; Guo et al., 2007; Kim et al., 2007; Chang et al., 2009). Unlike the flowers in eudicots, the nearly identical shape of the sepals and petals as well as the production of a unique lip in orchid flowers make them a very special plant species for the study of flower development. Four clades (1–4) of genes in the paleoAP3 lineage have been identified in several orchids (Hsu and Yang, 2002; Tsai et al., 2004; Kim et al., 2007; Mondragón-Palomino and Theißen, 2008, 2009; Mondragón-Palomino et al., 2009). Several works have described the possible interactions among these four clades of paleoAP3 genes and one PI gene that are involved in regulating the differentiation and formation of the sepal/petal/lip of orchids (Tsai et al., 2004; Kim et al., 2007; Mondragón-Palomino and Theißen, 2008, 2009). However, the exact mechanism that involves the orchid B class genes remains unclear and needs to be clarified by more experimental investigations.O. Gower Ramsey is a popular orchid with important economic value in cut flower markets. Only a few studies have been reported on the role of MADS box genes in regulating flower formation in this plant species (Hsu and Yang, 2002; Hsu et al., 2003; Chang et al., 2009). An AP3-like MADS gene that regulates both floral formation and initiation in transgenic Arabidopsis has been reported (Hsu and Yang, 2002). In addition, four AP1/AGAMOUS-LIKE9 (AGL9)-like MADS box genes have been characterized that show novel expression patterns and cause different effects on floral transition and formation in Arabidopsis (Hsu et al., 2003; Chang et al., 2009). Compared with other orchids, the production of a large and well-expanded lip and five small identical sepals/petals makes O. Gower Ramsey a special case for the study of the diverse functions of B class MADS box genes during evolution. Therefore, the isolation of more B class MADS box genes and further study of their roles in the regulation of perianth (sepal/petal/lip) formation during O. Gower Ramsey flower development are necessary. In addition to the clade 2 paleoAP3 gene OMADS3, which was previously characterized in our laboratory (Hsu and Yang, 2002), three more B class MADS box genes, OMADS5, OMADS8, and OMADS9, were characterized from O. Gower Ramsey in this study. Based on the different expression patterns and the protein interactions among these four orchid B class genes, we propose that the presence of OMADS3/8 and/or OMADS9 is required for sepal/petal/lip formation. Further sepal and petal formation at least requires the additional presence of OMADS5, whereas large lip formation was seen when OMADS5 expression was absent. Our results provide a new finding and information pertaining to the roles for orchid B class MADS box genes in the regulation of sepal/petal/lip formation.     

Overexpression of BAX INHIBITOR-1 Links Plasma Membrane Microdomain Proteins to Stress

BAX INHIBITOR-1 (BI-1) is a cell death suppressor widely conserved in plants and animals. Overexpression of BI-1 enhances tolerance to stress-induced cell death in plant cells, although the molecular mechanism behind this enhancement is unclear. We recently found that Arabidopsis (Arabidopsis thaliana) BI-1 is involved in the metabolism of sphingolipids, such as the synthesis of 2-hydroxy fatty acids, suggesting the involvement of sphingolipids in the cell death regulatory mechanism downstream of BI-1. Here, we show that BI-1 affects cell death-associated components localized in sphingolipid-enriched microdomains of the plasma membrane in rice (Oryza sativa) cells. The amount of 2-hydroxy fatty acid-containing glucosylceramide increased in the detergent-resistant membrane (DRM; a biochemical counterpart of plasma membrane microdomains) fraction obtained from BI-1-overexpressing rice cells. Comparative proteomics analysis showed quantitative changes of DRM proteins in BI-1-overexpressing cells. In particular, the protein abundance of FLOTILLIN HOMOLOG (FLOT) and HYPERSENSITIVE-INDUCED REACTION PROTEIN3 (HIR3) markedly decreased in DRM of BI-1-overexpressing cells. Loss-of-function analysis demonstrated that FLOT and HIR3 are required for cell death by oxidative stress and salicylic acid, suggesting that the decreased levels of these proteins directly contribute to the stress-tolerant phenotypes in BI-1-overexpressing rice cells. These findings provide a novel biological implication of plant membrane microdomains in stress-induced cell death, which is negatively modulated by BI-1 overexpression via decreasing the abundance of a set of key proteins involved in cell death.BAX INHIBITOR-1 (BI-1) is an endoplasmic reticulum (ER)-based cell death suppressor widely conserved in plants and animals (Xu and Reed, 1998; Kawai et al., 1999). In plants, BI-1 is considered a stress-associated factor, since its expression is stimulated by various stresses (Sanchez et al., 2000; Kawai-Yamada et al., 2001; Matsumura et al., 2003; Watanabe and Lam, 2006; Isbat et al., 2009). Although plants lack the homolog of animal BAX as an inducer of programmed cell death, loss of BI-1 expression results in a severe cell death phenotype under stress conditions, such as fumonisin B1-induced ER stress and disturbance of ion homeostasis (Watanabe and Lam, 2006; Ihara-Ohori et al., 2007). Conversely, plants overexpressing BI-1 exhibit tolerance to cell death induced by various stresses (Kawai-Yamada et al., 2001, 2004; Matsumura et al., 2003; Ihara-Ohori et al., 2007; Watanabe and Lam, 2008; Ishikawa et al., 2010). Moreover, BI-1 overexpression confers not only tolerance to oxidative stress-mediated cell death but also enhanced metabolic acclimation involved in energy and redox balance (Ishikawa et al., 2010). The results of these studies indicate that plant BI-1 is potentially useful for engineering stress-tolerant plants. However, little is known about the mode of action of BI-1 in the cell death regulatory pathway (Ishikawa et al., 2011). While overexpression systems sometimes include artificial or off-site effects, the observation that BI-1 overexpression improves stress tolerance suggests the importance of dissecting plants overexpressing it to further address the molecular basis of BI-1 function and cell death and stress tolerance management.As another approach to understand the molecular function of BI-1, screening of candidates interacting biochemically or functionally with BI-1 has been performed. First, Arabidopsis (Arabidopsis thaliana) BI-1 was confirmed to bind to calmodulin, like barley (Hordeum vulgare) MLO protein, a membrane-bound cell death regulator (Kim et al., 2002; Ihara-Ohori et al., 2007). Since the calmodulin-binding ability of BI-1 and MLO is necessary for their cell death-suppressing activity, Ca²⁺ signaling is critically involved in BI-1- and MLO-mediated cell death regulation (Kim et al., 2002; Kawai-Yamada et al., 2009). More recently, it was also demonstrated that the cell death suppression by BI-1 is mediated, at least in part, through fatty acid hydroxylase (FAH) in a Saccharomyces cerevisiae ectopic expression system (Nagano et al., 2009). In addition, Arabidopsis FAHs (AtFAH1 and AtFAH2) interact with BI-1 via cytochrome b₅ at the ER, resulting in the accumulation of 2-hydroxy fatty acids (2-HFAs) in Arabidopsis plants overexpressing BI-1. 2-HFAs are typical components of the ceramide backbone of sphingolipids (Imai et al., 1995; Pata et al., 2010). Although many functions of plant sphingolipids remain to be elucidated, accumulating evidence clearly indicates that sphingolipids and their metabolism are closely involved in cell death regulation and various stress responses in plants (Ng et al., 2001; Liang et al., 2003; Townley et al., 2005; Chen et al., 2008, 2012; Wang et al., 2008; Saucedo-García et al., 2011; Dutilleul et al., 2012; Kӧnig et al., 2012; Nagano et al., 2012; Mortimer et al., 2013), implying that BI-1 plays a role in cell death regulation through sphingolipid metabolism. Sphingolipids are major components of membrane lipids and are at particularly high concentrations in membrane microdomains, known as lipid rafts in animal cells, which are essential for membrane-mediated signaling and act as a sorting platform for targeted protein traffic (Simons and Toomre, 2000; Staubach and Hanisch, 2011). In mammalian cells, sphingomyelin metabolism in lipid rafts plays a vital role in the initiation of apoptotic cell death (Milhas et al., 2010). Recent studies have demonstrated the presence of raft-like membrane microdomains in plant cells and a role for them in defense responses and targeted protein sorting (Peskan et al., 2000; Fujiwara et al., 2009; Minami et al., 2009; Melser et al., 2010; Markham et al., 2011).This study focused on membrane microdomains in relation to BI-1-mediated sphingolipid metabolism. Our findings indicated that BI-1 alters sphingolipid composition in membrane microdomains, and this is accompanied by dynamic changes in a number of detergent-resistant membrane (DRM) proteins involved in cell death regulation.     

Calcineurin B-Like Protein-Interacting Protein Kinase CIPK21 Regulates Osmotic and Salt Stress Responses in Arabidopsis

The role of calcium-mediated signaling has been extensively studied in plant responses to abiotic stress signals. Calcineurin B-like proteins (CBLs) and CBL-interacting protein kinases (CIPKs) constitute a complex signaling network acting in diverse plant stress responses. Osmotic stress imposed by soil salinity and drought is a major abiotic stress that impedes plant growth and development and involves calcium-signaling processes. In this study, we report the functional analysis of CIPK21, an Arabidopsis (Arabidopsis thaliana) CBL-interacting protein kinase, ubiquitously expressed in plant tissues and up-regulated under multiple abiotic stress conditions. The growth of a loss-of-function mutant of CIPK21, cipk21, was hypersensitive to high salt and osmotic stress conditions. The calcium sensors CBL2 and CBL3 were found to physically interact with CIPK21 and target this kinase to the tonoplast. Moreover, preferential localization of CIPK21 to the tonoplast was detected under salt stress condition when coexpressed with CBL2 or CBL3. These findings suggest that CIPK21 mediates responses to salt stress condition in Arabidopsis, at least in part, by regulating ion and water homeostasis across the vacuolar membranes.Drought and salinity cause osmotic stress in plants and severely affect crop productivity throughout the world. Plants respond to osmotic stress by changing a number of cellular processes (Xiong et al., 1999; Xiong and Zhu, 2002; Bartels and Sunkar, 2005; Boudsocq and Lauriére, 2005). Some of these changes include activation of stress-responsive genes, regulation of membrane transport at both plasma membrane (PM) and vacuolar membrane (tonoplast) to maintain water and ionic homeostasis, and metabolic changes to produce compatible osmolytes such as Pro (Stewart and Lee, 1974; Krasensky and Jonak, 2012). It has been well established that a specific calcium (Ca²⁺) signature is generated in response to a particular environmental stimulus (Trewavas and Malhó, 1998; Scrase-Field and Knight, 2003; Luan, 2009; Kudla et al., 2010). The Ca²⁺ changes are primarily perceived by several Ca²⁺ sensors such as calmodulin (Reddy, 2001; Luan et al., 2002), Ca²⁺-dependent protein kinases (Harper and Harmon, 2005), calcineurin B-like proteins (CBLs; Luan et al., 2002; Batistič and Kudla, 2004; Pandey, 2008; Luan, 2009; Sanyal et al., 2015), and other Ca²⁺-binding proteins (Reddy, 2001; Shao et al., 2008) to initiate various cellular responses.Plant CBL-type Ca²⁺ sensors interact with and activate CBL-interacting protein kinases (CIPKs) that phosphorylate downstream components to transduce Ca²⁺ signals (Liu et al., 2000; Luan et al., 2002; Batistič and Kudla, 2004; Luan, 2009). In several plant species, multiple members have been identified in the CBL and CIPK family (Luan et al., 2002; Kolukisaoglu et al., 2004; Pandey, 2008; Batistič and Kudla, 2009; Weinl and Kudla, 2009; Pandey et al., 2014). Involvement of specific CBL-CIPK pair to decode a particular type of signal entails the alternative and selective complex formation leading to stimulus-response coupling (D’Angelo et al., 2006; Batistič et al., 2010).Several CBL and CIPK family members have been implicated in plant responses to drought, salinity, and osmotic stress based on genetic analysis of Arabidopsis (Arabidopsis thaliana) mutants (Zhu, 2002; Cheong et al., 2003, 2007; Kim et al., 2003; Pandey et al., 2004, 2008; D’Angelo et al., 2006; Qin et al., 2008; Tripathi et al., 2009; Held et al., 2011; Tang et al., 2012; Drerup et al., 2013; Eckert et al., 2014). A few CIPKs have also been functionally characterized by gain-of-function approach in crop plants such as rice (Oryza sativa), pea (Pisum sativum), and maize (Zea mays) and were found to be involved in osmotic stress responses (Mahajan et al., 2006; Xiang et al., 2007; Yang et al., 2008; Tripathi et al., 2009; Zhao et al., 2009; Cuéllar et al., 2010).In this report, we examined the role of the Arabidopsis CIPK21 gene in osmotic stress response by reverse genetic analysis. The loss-of-function mutant plants became hypersensitive to salt and mannitol stress conditions, suggesting that CIPK21 is involved in the regulation of osmotic stress response in Arabidopsis. These findings are further supported by an enhanced tonoplast targeting of the cytoplasmic CIPK21 through interaction with the vacuolar Ca²⁺ sensors CBL2 and CBL3 under salt stress condition.     

Patterns of Metabolite Changes Identified from Large-Scale Gene Perturbations in Arabidopsis Using a Genome-Scale Metabolic Network

Metabolomics enables quantitative evaluation of metabolic changes caused by genetic or environmental perturbations. However, little is known about how perturbing a single gene changes the metabolic system as a whole and which network and functional properties are involved in this response. To answer this question, we investigated the metabolite profiles from 136 mutants with single gene perturbations of functionally diverse Arabidopsis (Arabidopsis thaliana) genes. Fewer than 10 metabolites were changed significantly relative to the wild type in most of the mutants, indicating that the metabolic network was robust to perturbations of single metabolic genes. These changed metabolites were closer to each other in a genome-scale metabolic network than expected by chance, supporting the notion that the genetic perturbations changed the network more locally than globally. Surprisingly, the changed metabolites were close to the perturbed reactions in only 30% of the mutants of the well-characterized genes. To determine the factors that contributed to the distance between the observed metabolic changes and the perturbation site in the network, we examined nine network and functional properties of the perturbed genes. Only the isozyme number affected the distance between the perturbed reactions and changed metabolites. This study revealed patterns of metabolic changes from large-scale gene perturbations and relationships between characteristics of the perturbed genes and metabolic changes.Rational and quantitative assessment of metabolic changes in response to genetic modification (GM) is an open question and in need of innovative solutions. Nontargeted metabolite profiling can detect thousands of compounds, but it is not easy to understand the significance of the changed metabolites in the biochemical and biological context of the organism. To better assess the changes in metabolites from nontargeted metabolomics studies, it is important to examine the changed metabolites in the context of the genome-scale metabolic network of the organism.Metabolomics is a technique that aims to quantify all the metabolites in a biological system (Nikolau and Wurtele, 2007; Nicholson and Lindon, 2008; Roessner and Bowne, 2009). It has been used widely in studies ranging from disease diagnosis (Holmes et al., 2008; DeBerardinis and Thompson, 2012) and drug discovery (Cascante et al., 2002; Kell, 2006) to metabolic reconstruction (Feist et al., 2009; Kim et al., 2012) and metabolic engineering (Keasling, 2010; Lee et al., 2011). Metabolomic studies have demonstrated the possibility of identifying gene functions from changes in the relative concentrations of metabolites (metabotypes or metabolic signatures; Ebbels et al., 2004) in various species including yeast (Saccharomyces cerevisiae; Raamsdonk et al., 2001; Allen et al., 2003), Arabidopsis (Arabidopsis thaliana; Brotman et al., 2011), tomato (Solanum lycopersicum; Schauer et al., 2006), and maize (Zea mays; Riedelsheimer et al., 2012). Metabolomics has also been used to better understand how plants interact with their environments (Field and Lake, 2011), including their responses to biotic and abiotic stresses (Dixon et al., 2006; Arbona et al., 2013), and to predict important agronomic traits (Riedelsheimer et al., 2012). Metabolite profiling has been performed on many plant species, including angiosperms such as Arabidopsis, poplar (Populus trichocarpa), and Catharanthus roseus (Sumner et al., 2003; Rischer et al., 2006), basal land plants such as Selaginella moellendorffii and Physcomitrella patens (Erxleben et al., 2012; Yobi et al., 2012), and Chlamydomonas reinhardtii (Fernie et al., 2012; Davis et al., 2013). With the availability of whole genome sequences of various species, metabolomics has the potential to become a useful tool for elucidating the functions of genes using large-scale systematic analyses (Fiehn et al., 2000; Saito and Matsuda, 2010; Hur et al., 2013).Although metabolomics data have the potential for identifying the roles of genes that are associated with metabolic phenotypes, the biochemical mechanisms that link functions of genes with metabolic phenotypes are still poorly characterized. For example, we do not yet know the principles behind how perturbing the expression of a single gene changes the metabolic system as a whole. Large-scale metabolomics data have provided useful resources for linking phenotypes to genotypes (Fiehn et al., 2000; Roessner et al., 2001; Tikunov et al., 2005; Schauer et al., 2006; Lu et al., 2011; Fukushima et al., 2014). For example, Lu et al. (2011) compared morphological and metabolic phenotypes from more than 5,000 Arabidopsis chloroplast mutants using gas chromatography (GC)- and liquid chromatography (LC)-mass spectrometry (MS). Fukushima et al. (2014) generated metabolite profiles from various characterized and uncharacterized mutant plants and clustered the mutants with similar metabolic phenotypes by conducting multidimensional scaling with quantified metabolic phenotypes. Nonetheless, representation and analysis of such a large amount of data remains a challenge for scientific discovery (Lu et al., 2011). In addition, these studies do not examine the topological and functional characteristics of metabolic changes in the context of a genome-scale metabolic network. To understand the relationship between genotype and metabolic phenotype, we need to investigate the metabolic changes caused by perturbing the expression of a gene in a genome-scale metabolic network perspective, because metabolic pathways are not independent biochemical factories but are components of a complex network (Berg et al., 2002; Merico et al., 2009).Much progress has been made in the last 2 decades to represent metabolism at a genome scale (Terzer et al., 2009). The advances in genome sequencing and emerging fields such as biocuration and bioinformatics enabled the representation of genome-scale metabolic network reconstructions for model organisms (Bassel et al., 2012). Genome-scale metabolic models have been built and applied broadly from microbes to plants. The first step toward modeling a genome-scale metabolism in a plant species started with developing a genome-scale metabolic pathway database for Arabidopsis (AraCyc; Mueller et al., 2003) from reference pathway databases (Kanehisa and Goto, 2000; Karp et al., 2002; Zhang et al., 2010). Genome-scale metabolic pathway databases have been built for several plant species (Mueller et al., 2005; Zhang et al., 2005, 2010; Urbanczyk-Wochniak and Sumner, 2007; May et al., 2009; Dharmawardhana et al., 2013; Monaco et al., 2013, 2014; Van Moerkercke et al., 2013; Chae et al., 2014; Jung et al., 2014). Efforts have been made to develop predictive genome-scale metabolic models using enzyme kinetics and stoichiometric flux-balance approaches (Sweetlove et al., 2008). de Oliveira Dal’Molin et al. (2010) developed a genome-scale metabolic model for Arabidopsis and successfully validated the model by predicting the classical photorespiratory cycle as well as known key differences between redox metabolism in photosynthetic and nonphotosynthetic plant cells. Other genome-scale models have been developed for Arabidopsis (Poolman et al., 2009; Radrich et al., 2010; Mintz-Oron et al., 2012), C. reinhardtii (Chang et al., 2011; Dal’Molin et al., 2011), maize (Dal’Molin et al., 2010; Saha et al., 2011), sorghum (Sorghum bicolor; Dal’Molin et al., 2010), and sugarcane (Saccharum officinarum; Dal’Molin et al., 2010). These predictive models have the potential to be applied broadly in fields such as metabolic engineering, drug target discovery, identification of gene function, study of evolutionary processes, risk assessment of genetically modified crops, and interpretations of mutant phenotypes (Feist and Palsson, 2008; Ricroch et al., 2011).Here, we interrogate the metabotypes caused by 136 single gene perturbations of Arabidopsis by analyzing the relative concentration changes of 1,348 chemically identified metabolites using a reconstructed genome-scale metabolic network. We examine the characteristics of the changed metabolites (the metabolites whose relative concentrations were significantly different in mutants relative to the wild type) in the metabolic network to uncover biological and topological consequences of the perturbed genes.     

Nontransgenic Genome Modification in Plant Cells

Zinc finger nucleases (ZFNs) are a powerful tool for genome editing in eukaryotic cells. ZFNs have been used for targeted mutagenesis in model and crop species. In animal and human cells, transient ZFN expression is often achieved by direct gene transfer into the target cells. Stable transformation, however, is the preferred method for gene expression in plant species, and ZFN-expressing transgenic plants have been used for recovery of mutants that are likely to be classified as transgenic due to the use of direct gene-transfer methods into the target cells. Here we present an alternative, nontransgenic approach for ZFN delivery and production of mutant plants using a novel Tobacco rattle virus (TRV)-based expression system for indirect transient delivery of ZFNs into a variety of tissues and cells of intact plants. TRV systemically infected its hosts and virus ZFN-mediated targeted mutagenesis could be clearly observed in newly developed infected tissues as measured by activation of a mutated reporter transgene in tobacco (Nicotiana tabacum) and petunia (Petunia hybrida) plants. The ability of TRV to move to developing buds and regenerating tissues enabled recovery of mutated tobacco and petunia plants. Sequence analysis and transmission of the mutations to the next generation confirmed the stability of the ZFN-induced genetic changes. Because TRV is an RNA virus that can infect a wide range of plant species, it provides a viable alternative to the production of ZFN-mediated mutants while avoiding the use of direct plant-transformation methods.Methods for genome editing in plant cells have fallen behind the remarkable progress made in whole-genome sequencing projects. The availability of reliable and efficient methods for genome editing would foster gene discovery and functional gene analyses in model plants and the introduction of novel traits in agriculturally important species (Puchta, 2002; Hanin and Paszkowski, 2003; Reiss, 2003; Porteus, 2009). Genome editing in various species is typically achieved by integrating foreign DNA molecules into the target genome by homologous recombination (HR). Genome editing by HR is routine in yeast (Saccharomyces cerevisiae) cells (Scherer and Davis, 1979) and has been adapted for other species, including Drosophila, human cell lines, various fungal species, and mouse embryonic stem cells (Baribault and Kemler, 1989; Venken and Bellen, 2005; Porteus, 2007; Hall et al., 2009; Laible and Alonso-González, 2009; Tenzen et al., 2009). In plants, however, foreign DNA molecules, which are typically delivered by direct gene-transfer methods (e.g. Agrobacterium and microbombardment of plasmid DNA), often integrate into the target cell genome via nonhomologous end joining (NHEJ) and not HR (Ray and Langer, 2002; Britt and May, 2003).Various methods have been developed to indentify and select for rare site-specific foreign DNA integration events or to enhance the rate of HR-mediated DNA integration in plant cells. Novel T-DNA molecules designed to support strong positive- and negative-selection schemes (e.g. Thykjaer et al., 1997; Terada et al., 2002), altering the plant DNA-repair machinery by expressing yeast chromatin remodeling protein (Shaked et al., 2005), and PCR screening of large numbers of transgenic plants (Kempin et al., 1997; Hanin et al., 2001) are just a few of the experimental approaches used to achieve HR-mediated gene targeting in plant species. While successful, these approaches, and others, have resulted in only a limited number of reports describing the successful implementation of HR-mediated gene targeting of native and transgenic sequences in plant cells (for review, see Puchta, 2002; Hanin and Paszkowski, 2003; Reiss, 2003; Porteus, 2009; Weinthal et al., 2010).HR-mediated gene targeting can potentially be enhanced by the induction of genomic double-strand breaks (DSBs). In their pioneering studies, Puchta et al. (1993, 1996) showed that DSB induction by the naturally occurring rare-cutting restriction enzyme I-SceI leads to enhanced HR-mediated DNA repair in plants. Expression of I-SceI and another rare-cutting restriction enzyme (I-CeuI) also led to efficient NHEJ-mediated site-specific mutagenesis and integration of foreign DNA molecules in plants (Salomon and Puchta, 1998; Chilton and Que, 2003; Tzfira et al., 2003). Naturally occurring rare-cutting restriction enzymes thus hold great promise as a tool for genome editing in plant cells (Carroll, 2004; Pâques and Duchateau, 2007). However, their wide application is hindered by the tedious and next to impossible reengineering of such enzymes for novel DNA-target specificities (Pâques and Duchateau, 2007).A viable alternative to the use of rare-cutting restriction enzymes is the zinc finger nucleases (ZFNs), which have been used for genome editing in a wide range of eukaryotic species, including plants (e.g. Bibikova et al., 2001; Porteus and Baltimore, 2003; Lloyd et al., 2005; Urnov et al., 2005; Wright et al., 2005; Beumer et al., 2006; Moehle et al., 2007; Santiago et al., 2008; Shukla et al., 2009; Tovkach et al., 2009; Townsend et al., 2009; Osakabe et al., 2010; Petolino et al., 2010; Zhang et al., 2010). Here too, ZFNs have been used to enhance DNA integration via HR (e.g. Shukla et al., 2009; Townsend et al., 2009) and as an efficient tool for the induction of site-specific mutagenesis (e.g. Lloyd et al., 2005; Zhang et al., 2010) in plant species. The latter is more efficient and simpler to implement in plants as it does not require codelivery of both ZFN-expressing and donor DNA molecules and it relies on NHEJ—the dominant DNA-repair machinery in most plant species (Ray and Langer, 2002; Britt and May, 2003).ZFNs are artificial restriction enzymes composed of a fusion between an artificial Cys₂His₂ zinc-finger protein DNA-binding domain and the cleavage domain of the FokI endonuclease. The DNA-binding domain of ZFNs can be engineered to recognize a variety of DNA sequences (for review, see Durai et al., 2005; Porteus and Carroll, 2005; Carroll et al., 2006). The FokI endonuclease domain functions as a dimer, and digestion of the target DNA requires proper alignment of two ZFN monomers at the target site (Durai et al., 2005; Porteus and Carroll, 2005; Carroll et al., 2006). Efficient and coordinated expression of both monomers is thus required for the production of DSBs in living cells. Transient ZFN expression, by direct gene delivery, is the method of choice for targeted mutagenesis in human and animal cells (e.g. Urnov et al., 2005; Beumer et al., 2006; Meng et al., 2008). Among the different methods used for high and efficient transient ZFN delivery in animal and human cell lines are plasmid injection (Morton et al., 2006; Foley et al., 2009), direct plasmid transfer (Urnov et al., 2005), the use of integrase-defective lentiviral vectors (Lombardo et al., 2007), and mRNA injection (Takasu et al., 2010).In plant species, however, efficient and strong gene expression is often achieved by stable gene transformation. Both transient and stable ZFN expression have been used in gene-targeting experiments in plants (Lloyd et al., 2005; Wright et al., 2005; Maeder et al., 2008; Cai et al., 2009; de Pater et al., 2009; Shukla et al., 2009; Tovkach et al., 2009; Townsend et al., 2009; Osakabe et al., 2010; Petolino et al., 2010; Zhang et al., 2010). In all cases, direct gene-transformation methods, using polyethylene glycol, silicon carbide whiskers, or Agrobacterium, were deployed. Thus, while mutant plants and tissues could be recovered, potentially without any detectable traces of foreign DNA, such plants were generated using a transgenic approach and are therefore still likely to be classified as transgenic. Furthermore, the recovery of mutants in many cases is also dependent on the ability to regenerate plants from protoplasts, a procedure that has only been successfully applied in a limited number of plant species. Therefore, while ZFN technology is a powerful tool for site-specific mutagenesis, its wider implementation for plant improvement may be somewhat limited, both by its restriction to certain plant species and by legislative restrictions imposed on transgenic plants.Here we describe an alternative to direct gene transfer for ZFN delivery and for the production of mutated plants. Our approach is based on the use of a novel Tobacco rattle virus (TRV)-based expression system, which is capable of systemically infecting its host and spreading into a variety of tissues and cells of intact plants, including developing buds and regenerating tissues. We traced the indirect ZFN delivery in infected plants by activation of a mutated reporter gene and we demonstrate that this approach can be used to recover mutated plants.     

Plasma Membrane Localization of Solanum tuberosum Remorin from Group 1, Homolog 3 Is Mediated by Conformational Changes in a Novel C-Terminal Anchor and Required for the Restriction of Potato Virus X Movement]

The formation of plasma membrane (PM) microdomains plays a crucial role in the regulation of membrane signaling and trafficking. Remorins are a plant-specific family of proteins organized in six phylogenetic groups, and Remorins of group 1 are among the few plant proteins known to specifically associate with membrane rafts. As such, they are valuable to understand the molecular bases for PM lateral organization in plants. However, little is known about the structural determinants underlying the specific association of group 1 Remorins with membrane rafts. We used a structure-function approach to identify a short C-terminal anchor (RemCA) indispensable and sufficient for tight direct binding of potato (Solanum tuberosum) REMORIN 1.3 (StREM1.3) to the PM. RemCA switches from unordered to α-helical structure in a nonpolar environment. Protein structure modeling indicates that RemCA folds into a tight hairpin of amphipathic helices. Consistently, mutations reducing RemCA amphipathy abolished StREM1.3 PM localization. Furthermore, RemCA directly binds to biological membranes in vitro, shows higher affinity for Detergent-Insoluble Membranes lipids, and targets yellow fluorescent protein to Detergent-Insoluble Membranes in vivo. Mutations in RemCA resulting in cytoplasmic StREM1.3 localization abolish StREM1.3 function in restricting potato virus X movement. The mechanisms described here provide new insights on the control and function of lateral segregation of plant PM.Protein-lipid interactions are increasingly recognized as key regulatory processes for signal perception and cellular signaling cascades (Cho and Stahelin, 2005). During signal transduction and trafficking, a number of soluble proteins dynamically associate with plasma membranes (PMs) to deliver their cargo and to recruit pathway components to the sites of action (Seong et al., 2011). For such proteins, membrane association can be critical for function (Porter and Koelle, 2010).PM targeting of peripheral proteins is achieved through (1) binding to integral membrane proteins, (2) posttranslational modifications, or (3) directly by intrinsic membrane anchor domains. Posttranslational modifications function as auxiliary modifications for transient or weak association of soluble proteins to the intracellular face of the PM. In plants, these include N-myristoylation, S-palmitoylation, prenylation by farnesyl or geranylgeranyl moieties, or attachment of glycosylphosphatidylinositol (GPI) anchors (Thompson and Okuyama, 2000). GPI anchors, for example, tightly associate proteins to the extracellular face of PMs by interaction of the inositol head group of the membrane lipid phosphatidylinositol with a glucosamine residue linked to the C-terminal amino acid of the protein (Paulick and Bertozzi, 2008). As an alternative mechanism, globular structures either recognize phospholipids in a stereospecific manner or associate with membranes by their biophysical properties (for review, see Lemmon, 2008). Other proteins expose unstructured clusters of basic and hydrophobic residues to mediate PM binding (McLaughlin et al., 2002; McLaughlin and Murray, 2005).Selective recognition of membrane compartments or domains by protein anchors can be critical in triggering the appropriate downstream trafficking and signaling events (for review, see Gruenberg, 2003; De Matteis and Godi, 2004). Membrane domain selectivity can be specified by the anchoring posttranslational modification or by a protein anchor domain. For instance, proteins carrying GPI anchors are overrepresented in membrane rafts, indicating that addition of this lipid anchor directs proteins to these microdomains (Cordy et al., 2003; Kierszniowska et al., 2009). Membrane rafts are enriched in highly saturated long-chain sphingolipids, sterols, and saturated phospholipids, creating tightly packed domains, designated as “liquid ordered.” These lipids display a stronger affinity to saturated acyl chains as found in GPI-anchored and acylated proteins (Brown, 2006). The composition of membrane rafts also prevents solubilization by detergent at low temperature with nonionic detergent and allows the partial purification of rafts in so-called Detergent-Insoluble Membrane (DIM) fractions, which are supposedly biochemical counterparts of membrane rafts. Many signaling proteins are found in membrane rafts, supporting the hypothesis that they serve as key platforms for cellular signal transduction and cell-to-cell communication (Lingwood and Simons, 2010; Simon-Plas et al., 2011). For example, in human (Homo sapiens) cells, key soluble signaling components such as the Ser/Thr kinase Akt (protein kinase B) are recruited to membrane rafts where they activate signal transduction cascades (Lasserre et al., 2008). Nevertheless, few protein motifs were described to contribute to raft targeting (Rossin et al., 2010), although the 6-amino-acid-long raft target signal from human Tyr phosphatase Src homology 2-containing phosphatase1 is the only known motif sufficient to anchor soluble proteins specifically to domains of the intracellular face of the PM (Sankarshanan et al., 2007). However, the anchoring mechanism itself remains to be unraveled in plants even though DIMs also exist (Mongrand et al., 2010) and functional PM domains have been reported (Bhat et al., 2005). The molecular basis for specific targeting and binding of proteins to membrane rafts has never been described.Remorins form a diverse family of plant-specific proteins organized in six distinct phylogenetic groups (Raffaele et al., 2007). Remorins from group 1 have been reported to localize to the PM despite their overall hydrophilic nature (Reymond et al., 1996; Raffaele et al., 2007). Moreover, group 1 Remorins almost exclusively associate to DIMs and localize to membrane microdomains in a sterol-dependent manner (Lefebvre et al., 2007; Kierszniowska et al., 2009; Raffaele et al., 2009). The function of Remorins is mostly unknown, but we showed in a previous study that StREM1.3 (for potato (Solanum tuberosum) Remorin from group 1, homolog 3; initially described in Reymond et al., 1996) regulates cell-to-cell propagation of the potato virus X (PVX), likely by directly interacting with the viral movement protein Triple Gene Block protein1 (TGBp1; Raffaele et al., 2009). StREM1.3 localizes to the inner leaflet of PMs and along plasmodesmata, bridges connecting neighbor cells essential for cell-to-cell communication in plants (Maule, 2008). Other members of the Remorin family group 1 are likely involved in innate immune responses (Liu et al., 2009; Widjaja et al., 2009; Keinath et al., 2010). Remorins from group 2 are involved in the control of infection by symbiotic bacteria at nodular infection threads and the peribacteroid membrane (Lefebvre et al., 2010) These data suggest general roles for Remorins in regulating signaling in plant-microbe interactions (Jarsch and Ott, 2011).Elucidating the mechanisms driving StREM1.3 association with PM microdomains therefore provides a unique opportunity for understanding the regulation and function of membrane lateral segregation in plants. StREM1.3 does not contain predictable transmembrane or membrane-associated domains. The bases for its association to PMs and selective targeting to DIMs are unknown. Here we identified a novel membrane anchor domain required for StREM1.3 tight and direct association with the detergent-insoluble fraction of the PM. We combined biophysics, in silico analysis, and directed mutagenesis to unravel the molecular bases of StREM1.3 membrane binding and its biological significance in the control of PVX propagation.     

Large-Scale Phosphoprotein Analysis in Medicago truncatula Roots Provides Insight into in Vivo Kinase Activity in Legumes

Nitrogen fixation in legumes requires the development of root organs called nodules and their infection by symbiotic rhizobia. Over the last decade, Medicago truncatula has emerged as a major model plant for the analysis of plant-microbe symbioses and for addressing questions pertaining to legume biology. While the initiation of symbiosis and the development of nitrogen-fixing root nodules depend on the activation of a protein phosphorylation-mediated signal transduction cascade in response to symbiotic signals produced by the rhizobia, few sites of in vivo phosphorylation have previously been identified in M. truncatula. We have characterized sites of phosphorylation on proteins from M. truncatula roots, from both whole cell lysates and membrane-enriched fractions, using immobilized metal affinity chromatography and tandem mass spectrometry. Here, we report 3,457 unique phosphopeptides spanning 3,404 nonredundant sites of in vivo phosphorylation on 829 proteins in M. truncatula Jemalong A17 roots, identified using the complementary tandem mass spectrometry fragmentation methods electron transfer dissociation and collision-activated dissociation. With this being, to our knowledge, the first large-scale plant phosphoproteomic study to utilize electron transfer dissociation, analysis of the identified phosphorylation sites revealed phosphorylation motifs not previously observed in plants. Furthermore, several of the phosphorylation motifs, including LxKxxs and RxxSxxxs, have yet to be reported as kinase specificities for in vivo substrates in any species, to our knowledge. Multiple sites of phosphorylation were identified on several key proteins involved in initiating rhizobial symbiosis, including SICKLE, NUCLEOPORIN133, and INTERACTING PROTEIN OF DMI3. Finally, we used these data to create an open-access online database for M. truncatula phosphoproteomic data.Medicago truncatula has become a model for studying the biology of leguminous plants such as soybean (Glycine max), alfalfa (Medicago sativa), and clover (Trifolium spp.; Singh et al., 2007). Most members of this vast family have the ability to fix atmospheric nitrogen by virtue of an endosymbiotic association with rhizobial bacteria, through which legumes undergo nodulation, the process of forming root nodules (Jones et al., 2007). Legumes are central to modern agriculture and civilization because of their ability to grow in nitrogen-depleted soils and replenish nitrogen through crop rotation. Consequently, there is great interest in understanding the molecular events that allow legumes to recognize their symbionts, develop root nodules, and fix nitrogen. Nod factors are lipochitooligosaccharidic signals secreted by the rhizobia and are required, in most legumes, for intracellular infection and nodule development. In recent decades, an elegant combination of genetics, biochemistry, and cell biology has shown that Nod factors activate intricate signaling events within cells of legume roots, including protein phosphorylation cascades and intracellular ion fluxes (Oldroyd and Downie, 2008).Protein phosphorylation is a central mechanism of signal transfer in cells (Laugesen et al., 2006; Peck, 2006; Huber, 2007). Several characterized protein kinases are required for symbiosis signal transduction in M. truncatula roots (Lévy et al., 2004; Yoshida and Parniske, 2005; Smit et al., 2007). A recent antibody-based study of cultured M. truncatula cells observed protein phosphorylation changes at the proteomic level in response to fungal infection (Trapphoff et al., 2009); however, the target residues of the phosphorylation events were not determined. A variety of studies have determined in vitro phosphorylation sites on legume proteins and demonstrated the biological importance of the target residues by mutagenesis (Yoshida and Parniske, 2005; Arrighi et al., 2006; Lima et al., 2006; Miyahara et al., 2008; Yano et al., 2008). To our knowledge, only six sites of in vivo protein phosphorylation have been detected for M. truncatula (Laugesen et al., 2006; Lima et al., 2006; Wienkoop et al., 2008), demonstrating the need for the identification of endogenous protein phosphorylation sites in legume model organisms on a proteome-wide scale.While considerable advancements have been made in the global analysis of protein phosphorylation (Nita-Lazar et al., 2008; Macek et al., 2009; Piggee, 2009; Thingholm et al., 2009), phosphoproteomics in plants has lagged years behind that of the mammalian systems (Kersten et al., 2006, 2009; Peck, 2006), which have more fully sequenced genomes and better annotated protein predictions. Arabidopsis (Arabidopsis thaliana), the first plant genome sequenced (Arabidopsis Genome Initiative, 2000), is now predicted to have over 1,000 protein kinases (Finn et al., 2008), approximately twice as many as in human (Manning et al., 2002). Because many of the kinases in the commonly studied mammalian systems are not conserved in the plant kingdom, there is significant need for large-scale phosphoproteomic technologies to discern the intricacies of phosphorylation-mediated cell signaling in plants. With the high mass accuracy afforded by the linear ion trap-orbitrap hybrid mass spectrometer (Makarov et al., 2006; Yates et al., 2006), recent studies in Arabidopsis have reported 2,597 phosphopeptides from suspension cell culture (Sugiyama et al., 2008) and 3,029 phosphopeptides from seedlings (Reiland et al., 2009).All previous large-scale plant phosphoproteomic studies have relied solely on collision-activated dissociation (CAD) during tandem mass spectrometry (MS/MS) and have not taken advantage of the more recently developed methods (Kersten et al., 2009) electron capture dissociation (Kelleher et al., 1999) or electron transfer dissociation (ETD; Coon et al., 2004; Syka et al., 2004). Mapping sites of posttranslational modifications, such as phosphorylation, is often more straightforward using electron-based fragmentation methods, as they frequently produce a full spectrum of sequence-informative ions without causing neutral loss of the modifying functional groups (Meng et al., 2005; Chi et al., 2007; Khidekel et al., 2007; Molina et al., 2007; Wiesner et al., 2008; Chalkley et al., 2009; Swaney et al., 2009). With an ETD-enabled hybrid orbitrap mass spectrometer (McAlister et al., 2007, 2008), we previously compared the performance of CAD and ETD tandem MS for large-scale identification of phosphopeptides (Swaney et al., 2009). ETD identified a greater percentage of unique phosphopeptides and more frequently localized phosphorylation sites. Still, the low overlap of identified phosphopeptides indicates that the two methods are highly complementary. With this in mind, we recently developed a decision tree-driven tandem MS algorithm to select the optimal fragmentation method for each precursor (Swaney et al., 2008).Here, we utilize this technology to map sites of in vivo protein phosphorylation in roots of M. truncatula Jemalong A17 plants. Phosphoproteins, from both whole-cell lysate and membrane-enriched fractions, were analyzed after digestion with a variety of different enzymes individually. Utilizing the complementary fragmentation methods of ETD and CAD, we report 3,404 nonredundant phosphorylation sites at an estimated false discovery rate (FDR) of 1%. Analysis of these data revealed several phosphorylation motifs not previously observed in plants. The phosphorylation sites identified provide insight into the potential regulation of key proteins involved in rhizobial symbiosis, potential consensus sequences by which kinases recognize their substrates, and critical phosphorylation events that are conserved between plant species.