CYP86A33-Targeted Gene Silencing in Potato Tuber Alters Suberin Composition,Distorts Suberin Lamellae,and Impairs the Periderm's Water Barrier Function

Suberin is a cell wall lipid polyester found in the cork cells of the periderm offering protection against dehydration and pathogens. Its biosynthesis and assembly, as well as its contribution to the sealing properties of the periderm, are still poorly understood. Here, we report on the isolation of the coding sequence CYP86A33 and the molecular and physiological function of this gene in potato (Solanum tuberosum) tuber periderm. CYP86A33 was down-regulated in potato plants by RNA interference-mediated silencing. Periderm from CYP86A33-silenced plants revealed a 60% decrease in its aliphatic suberin load and greatly reduced levels of C18:1 ω-hydroxyacid (approximately 70%) and α,ω-diacid (approximately 90%) monomers in comparison with wild type. Moreover, the glycerol esterified to suberin was reduced by 60% in the silenced plants. The typical regular ultrastructure of suberin, consisting of dark and light lamellae, disappeared and the thickness of the suberin layer was clearly reduced. In addition, the water permeability of the periderm isolated from CYP86A33-silenced lines was 3.5 times higher than that of the wild type. Thus, our data provide convincing evidence for the involvement of ω-functional fatty acids in establishing suberin structure and function.Periderm, the boundary tissue that replaces the epidermis in the secondary organs of plants, provides efficient protection against dehydration, UV radiation, and pathogens (Esau, 1965). Periderm is regularly found in the bark of woody plants, but herbaceous plants may also form a well-developed periderm in roots, tubers, and the oldest portions of stem. The periderm has been widely studied in potato (Solanum tuberosum) tubers because of the latter''s great agronomic significance (Schmidt and Schönherr, 1982; Vogt et al., 1983; Lulai and Freeman, 2001; Sabba and Lulai, 2002). Shrinkage and flaccidity occur in tubers if the protection afforded by the periderm against water loss is compromised (Lulai et al., 2006). Suberin and waxes embedded into the suberin matrix are considered essential for periderm imperviousness (Franke and Schreiber, 2007). Chemically, suberin is a complex lipid polymer consisting of a fatty acid-derived domain (aliphatic suberin) cross-linked by ester bonds to a polyaromatic lignin-like domain (aromatic suberin; Kolattukudy, 2001; Bernards, 2002; Franke and Schreiber, 2007). Aliphatic suberin has been widely analyzed in potato periderm (Kolattukudy and Agrawal, 1974; Graça and Pereira, 2000; Schreiber et al., 2005). The main monomers released from potato aliphatic suberin are a mixture of ω-hydroxyacids and α,ω-diacids with chain lengths ranging from C16 to C28 (mainly C18), together with glycerol. Small amounts of monofunctional fatty acids, alcohols, and ferulic acid are also released. Waxes are complex mixtures of lipids extractable with chloroform that in potato periderm consist mostly of linear very-long-chain aliphatics up to C32 (Schreiber et al., 2005). Suberin is deposited in the cell wall as a continuous deposit or secondary cell wall that overlays the polysaccharide primary cell wall from the inside (Esau, 1965). These suberin deposits appear under the transmission electron microscope (TEM) as regularly alternating opaque and translucent lamellae (Schmidt and Schönherr, 1982). Although several molecular models for suberin have been proposed (Kolattukudy, 1980; Bernards, 2002; Graça and Santos, 2007), how the suberin and wax components are organized in the lamellated suberin secondary cell wall is a matter of debate (Graça and Santos, 2007). Moreover, to what extent suberin and wax deposition and composition determine sealing properties of periderm still remains unclear (Schreiber et al., 2005). Several studies confirm the importance of waxes for the diffusion barrier (Soliday et al., 1979; Vogt et al., 1983; Schreiber et al., 2005), but the significance of aliphatic suberin has hardly been investigated at all. Interestingly, an Arabidopsis (Arabidopsis thaliana) knockout mutant for the GLYCEROL-3-PHOSPHATE ACYLTRANSFERASE5 gene (GPAT5) with altered suberin showed higher tetrazolium salt permeability in the seed coat (Beisson et al., 2007).ω-Hydroxylation of fatty acids, a reaction carried out in plants by cytochrome P450 monooxygenases, is a crucial step in the biosynthesis of plant biopolymers (Kolattukudy, 1980; Nawrath, 2002). The Arabidopsis mutant lacerata, which shows phenotype defects compatible with a cutin deficiency, is defective in CYP86A8 encoding a fatty acid ω-hydroxylase (Wellesen et al., 2001). The aberrant induction of type three genes1 (att1) mutant, showing an altered cuticle ultrastructure and a higher transpiration rate than wild type, is defective in CYP86A2 and contains reduced amounts of ω-functionalized cutin monomers (Xiao et al., 2004). Moreover, a genome-wide study of cork oak (Quercus suber) bark highlighted a member of the cytochrome P450 of the CYP86A subfamily as a strong candidate gene for aliphatic suberin biosynthesis (Soler et al., 2007); and a role for CYP86A1 in the biosynthesis of suberin has recently been confirmed in the primary root of Arabidopsis knockout mutants (Li et al., 2007; Hofer et al., 2008). However, how the lack of fatty acid ω-hydroxylase activity may affect the structural features and sealing properties of suberin in periderm cell walls has not been documented.To provide experimental evidence of the role of CYP86A genes in periderm fine structure and water transpiration properties, especially quantitative permeability studies so far unexplored in Arabidopsis, we followed a strategy based on the potato (cv Desirée). The potato is a model of choice for such studies because transgenic plants can be produced in reasonable time and sufficient amounts of periderm can easily be obtained from tubers for chemical and physiological studies (Vogt et al., 1983; Schreiber et al., 2005). Based on the CYP86A gene that was highlighted in cork oak periderm as a strong suberin candidate for aliphatic suberin biosynthesis, we isolated the putative ortholog in potato and used a reverse genetic approach to analyze the effects of down-regulation on the chemical and ultrastructural features of suberin and on the physiological properties of the tuber periderm. Our findings indicate that down-regulation of CYP86A33, as this gene was designated, results in a strong decrease of the ω-functionalized monomers in aliphatic suberin, which are necessary for the suberin typical lamellar organization and for the periderm resistance to water loss.     

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.     

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.     

Suberin-Associated Fatty Alcohols in Arabidopsis: Distributions in Roots and Contributions to Seed Coat Barrier Properties

Suberin is found in a variety of tissues, such as root endoderms and periderms, storage tuber periderms, tree cork layer, and seed coats. It acts as a hydrophobic barrier to control the movement of water, gases, and solutes as well as an antimicrobial barrier. Suberin consists of polymerized phenolics, glycerol, and a variety of fatty acid derivatives, including primary fatty alcohols. We have conducted an in-depth analysis of the distribution of the C18:0 to C22:0 fatty alcohols in Arabidopsis (Arabidopsis thaliana) roots and found that only 20% are part of the root suberin polymer, together representing about 5% of its aliphatic monomer composition, while the remaining 80% are found in the nonpolymeric (soluble) fraction. Down-regulation of Arabidopsis FATTY ACYL REDUCTASE1 (FAR1), FAR4, and FAR5, which collectively produce the fatty alcohols found in suberin, reduced their levels by 70% to 80% in (1) the polymeric and nonpolymeric fractions from roots of tissue culture-grown plants, (2) the suberin-associated root waxes from 7-week-old soil-grown plants, and (3) the seed coat suberin polymer. By contrast, the other main monomers of suberin were not altered, indicating that reduced levels of fatty alcohols did not influence the suberin polymerization process. Nevertheless, the 75% reduction in total fatty alcohol and diol loads in the seed coat resulted in increased permeability to tetrazolium salts and a higher sensitivity to abscisic acid. These results suggest that fatty alcohols and diols play an important role in determining the functional properties of the seed coat suberin barrier.Suberin is a cell wall-linked polymeric barrier that plays a critical role in the survival of plants by protecting them against various biotic and abiotic stresses. It primarily acts as a hydrophobic barrier to control the movement of water, gases, and solutes, but also contributes to the strength of the cell wall (Ranathunge et al., 2011). Suberin is deposited at the inner face of primary cell walls next to the plasma membrane (Kolattukudy, 1980; Franke and Schreiber, 2007). It is typically found as lamellae (alternating dark and light bands when viewed by transmission electron microscopy) in the endodermis, exodermis, and peridermis of roots, as well as in the peridermis of underground storage tubers (Bernards, 2002). Suberin is also found in shoot periderms of trees (i.e. cork layer) and in seed coats (Molina et al., 2006, 2008) and is deposited in response to wounding (Kolattukudy, 2001).Suberin is a polymer consisting of aliphatics (fatty acid derivatives), phenolics, and glycerol. The predominant aliphatic components of suberin are ω-hydroxy fatty acids, α,ω-dicarboxylic acids, very-long-chain fatty acids, and primary fatty alcohols, while the major phenolic components are p-hydroxycinnamic acids, especially ferulic acid (Kolatukudy, 1980; Bernards et al., 1995; Pollard et al., 2008; Ranathunge et al., 2011). In the periderm of underground storage organs, suberin is found in association with waxes, which are isolated either by extensive extraction in solvent (Soliday et al., 1979; Serra et al., 2009) or by brief immersion of tubers in chloroform (Espelie et al., 1980). These suberin-associated waxes consist of linear aliphatics (e.g. alkanes, fatty acids, and fatty alcohols), which are similar to cuticular wax components of aerial tissues but generally of shorter chain lengths (Espelie et al., 1980). In waxes extracted from 3-week-old wounded potato (Solanum tuberosum) periderms, alkyl ferulates (i.e. ferulic acid linked by an ester bond to a C16:0–C32:0 fatty alcohol) represent up to 60% of the total wax load (Schreiber et al., 2005). Root waxes are also found in 6- to 7-week-old mature taproots of Arabidopsis (Arabidopsis thaliana) with a fully developed periderm (Li et al., 2007; Kosma et al., 2012). They are enriched in alkyl hydroxycinnamates (AHCs) made of C18:0 to C22:0 fatty alcohols esterified with coumaric, caffeic, or ferulic acids (Kosma et al., 2012). The monomer composition (in terms of major chemical species and chain length) of both suberin and suberin-associated waxes varies considerably between plant species, tissues, and developmental stages. Aliphatic suberin and suberin-associated waxes are considered the major contributors to the diffusion resistance of suberized cell walls to radial transport of water and solutes (Soliday et al., 1979; Espelie et al., 1980; Zimmermann et al., 2000; Ranathunge and Schreiber, 2011). The organization of suberin components in the lamellated structure as well as how waxes may be associated with the polymer is a matter of debate (Graça and Santos, 2007).Primary fatty alcohols are long-chain hydrocarbons containing a single hydroxyl group at the terminal position. They are ubiquitously detected as components of the suberin polymer, representing 1% to 10% of the total monomer mass recovered after transesterification (Holloway, 1983; Bernards, 2002; Pollard et al., 2008). Primary fatty alcohols are also typical components of suberin-associated waxes, where they can be found either in free form or linked by an ester bond with a hydroxycinnamic acid (i.e. as AHCs; Soliday et al., 1979; Espelie et al., 1980; Bernards and Lewis 1992; Li et al., 2007; Kosma et al., 2012). In mechanically isolated endodermis of soybean (Glycine max) roots, fatty alcohols represent about 1.5% and 0.2% of the total aliphatics found in suberin-associated waxes and suberin polymer, respectively (Thomas et al., 2007). In onion (Allium cepa) root exodermis, fatty alcohols (C14:0–C28:0) account for 7% to 12% of the soluble fraction, while the suberin fraction contains only C22:0 fatty alcohol, which makes up 3% of the suberin fraction across all exodermal maturation zones (Meyer et al., 2011). In suberizing potato periderms 7 d post wounding, C16:0 to C28:0 fatty alcohols represent about 10% and 18% of the total aliphatics in the insoluble poly(aliphatic) domain (suberin polymer) and in the soluble (nonpolymeric) fraction, respectively (Yang and Bernards, 2006). A similar study on native periderms from 21-d-stored potato (Serra et al., 2009) reported that fatty alcohols represent about 20% of the total aliphatic components found in the suberin polyester, while unlinked fatty alcohols and alkyl ferulates accounted for about 23% and 44% of the total aliphatics in the soluble waxes.In Arabidopsis, C18:0, C20:0, and C22:0 fatty alcohols account for slightly less than 3% of the polymerized aliphatics in roots of soil-grown plants (Domergue et al., 2010), but as much as 36% [w/w] of the soluble wax load (Li et al., 2007). Arabidopsis fatty acyl reductases FAR1 (At5g22500), FAR4 (At3g44540), and FAR5 (At3g44550) generate, respectively, the C22:0, C20:0, and C18:0 fatty alcohol present in the suberin of root, seed coat, and wounded leaf tissues (Domergue et al., 2010). These three enzymes also generate the C18:0 to C22:0 fatty alcohol components that make up AHCs of root waxes (Kosma et al., 2012). Although one particular chain length of primary alcohol was reduced in each far single mutant line (C18:0-OH, C20:0-OH, and C22:0-OH in far5, far4, and far1, respectively), the total fatty alcohol load of the suberin polymer and its composition were only slightly affected and mutant plants had no obvious developmental or physiological defects (Domergue et al., 2010). In this study, we report on the distribution of primary fatty alcohols in the soluble (nonpolymeric) and insoluble (suberin polymer) fractions from mature roots of Arabidopsis. We report that far double and triple mutants have highly reduced fatty alcohol levels, in a chain length-specific manner, in both fractions as well as in the seed coat suberin polymer. The significant reductions in total fatty alcohol and diol levels in the seed coat of these mutants lead to increased permeability and higher sensitivity to abscisic acid (ABA), bringing to light insights on the roles of fatty alcohols and diols in determining functional properties of suberin.     

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.     

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.     

The K-Segment of Maize DHN1 Mediates Binding to Anionic Phospholipid Vesicles and Concomitant Structural Changes

Dehydrins (DHNs; late embryogenesis abundant D11 family) are a family of intrinsically unstructured plant proteins that accumulate in the late stages of seed development and in vegetative tissues subjected to water deficit, salinity, low temperature, or abscisic acid treatment. We demonstrated previously that maize (Zea mays) DHNs bind preferentially to anionic phospholipid vesicles; this binding is accompanied by an increase in α-helicity of the protein, and adoption of α-helicity can be induced by sodium dodecyl sulfate. All DHNs contain at least one “K-segment,” a lysine-rich 15-amino acid consensus sequence. The K-segment is predicted to form a class A2 amphipathic α-helix, a structural element known to interact with membranes and proteins. Here, three K-segment deletion proteins of maize DHN1 were produced. Lipid vesicle-binding assays revealed that the K-segment is required for binding to anionic phospholipid vesicles, and adoption of α-helicity of the K-segment accounts for most of the conformational change of DHNs upon binding to anionic phospholipid vesicles or sodium dodecyl sulfate. The adoption of structure may help stabilize cellular components, including membranes, under stress conditions.When plants encounter environmental stresses such as drought or low temperature, various responses take place to adapt to these conditions. Typical responses include increased expression of chaperones, signal transduction pathway and late embryogenesis abundant (LEA) proteins, osmotic adjustment, and induction of degradation and repair systems (Ingram and Bartels, 1996).Dehydrins (DHNs; LEA D11 family) are a subfamily of group 2 LEA proteins that accumulate to high levels during late stages of seed development and in vegetative tissues subjected to water deficit, salinity, low temperature, or abscisic acid (ABA) treatment (Svensson et al., 2002). Some DHNs are expressed constitutively during normal growth (Nylander et al., 2001; Rorat et al., 2004, 2006; Rodriguez et al., 2005). DHNs exist in a wide range of photosynthetic organisms, including angiosperms, gymnosperms, algae, and mosses (Svensson et al., 2002). DHNs are encoded by a dispersed multigene family and are differentially regulated, at least in higher plants. For example, 13 Dhn genes have been identified in barley (Hordeum vulgare), dispersed over seven genetic map locations (Choi et al., 1999; Svensson et al., 2002) and regulated variably by drought, low temperature, and embryo development (Tommasini et al., 2008). DHNs are localized in various subcellular compartments, including cytosol (Roberts et al., 1993), nucleus (Houde et al., 1995), chloroplast (Artus et al., 1996), vacuole (Heyen et al., 2002), and proximal to the plasma membrane and protein bodies (Asghar et al., 1994; Egerton-Warburton et al., 1997; Puhakainen et al., 2004). Elevated expression of Dhn genes generally has been correlated with the acquisition of tolerance to abiotic stresses such as drought (Whitsitt et al., 1997), salt (Godoy et al., 1994; Jayaprakash et al., 1998), chilling (Ismail et al., 1999a), or freezing (Houde et al., 1995; Danyluk et al., 1998; Fowler et al., 2001). The differences in expression and tissue location suggest that individual members of the Dhn multigene family have somewhat distinct biological functions (Close, 1997; Zhu et al., 2000; Nylander et al., 2001). Many studies have observed a positive correlation between the accumulation of DHNs and tolerance to abiotic stresses (Svensson et al., 2002). However, overexpression of a single DHN protein has not, in general, been sufficient to confer stress tolerance (Puhakainen et al., 2004).DHNs are subclassified by sequence motifs referred to as the K-segment (Lys-rich consensus sequence), the Y-segment (N-terminal conserved sequence), the S-segment (a tract of Ser residues), and the φ-segment (Close, 1996). Because of high hydrophilicity, high content of Gly (>20%), and the lack of a defined three-dimensional structure in the pure form (Lisse et al., 1996), DHNs have been categorized as “intrinsically disordered/unstructured proteins” or “hydrophilins” (Wright and Dyson, 1999; Garay-Arroyo et al., 2000; Tompa, 2005; Kovacs et al., 2008). On the basis of compositional and biophysical properties and their link to abiotic stresses, several functions of DHNs have been proposed, including ion sequestration (Roberts et al., 1993), water retention (McCubbin et al., 1985), and stabilization of membranes or proteins (Close, 1996, 1997). Observations from in vitro experiments include DHN binding to lipid vesicles (Koag et al., 2003; Kovacs et al., 2008) or metals (Svensson et al., 2000; Heyen et al., 2002; Kruger et al., 2002; Alsheikh et al., 2003; Hara et al., 2005), protection of membrane lipid against peroxidation (Hara et al., 2003), retention of hydration or ion sequestration (Bokor et al., 2005; Tompa et al., 2006), and chaperone activity against the heat-induced inactivation and aggregation of various proteins (Kovacs et al., 2008).Intrinsically disordered/unstructured proteins that lack a well-defined three-dimensional structure have recently been recognized to be prevalent in prokaryotes and eukaryotes (Oldfield et al., 2005). They fulfill important functions in signal transduction, gene expression, and binding to targets such as protein, RNA, ions, and membranes (Wright and Dyson, 1999; Tompa, 2002; Dyson and Wright, 2005). The disorder confers structural flexibility and malleability to adapt to changes in the protein environment, including water potential, pH, ionic strength, and temperature, and to undergo structural transition when complexed with ligands such as other proteins, DNA, RNA, or membranes (Prestrelski et al., 1993; Uversky, 2002). Structural changes from disorder to ordered functional structure also can be induced by the folding of a partner protein (Wright and Dyson, 1999; Tompa, 2002; Mouillon et al., 2008).The idea that DHNs interact with membranes is consistent with many immunolocalization studies, which have shown that DHNs accumulate near the plasma membrane or membrane-rich areas surrounding lipid and protein bodies (Asghar et al., 1994; Egerton-Warburton et al., 1997; Danyluk et al., 1998; Puhakainen et al., 2004). The K-segment is predicted to form a class A2 amphipathic α-helix, in which hydrophilic and hydrophobic residues are arranged on opposite faces (Close, 1996). The amphipathic α-helix is a structural element known to interact with membranes and proteins (Epand et al., 1995). Also, in the presence of helical inducers such as SDS and trifluoroethanol (Dalal and Pio, 2006), DHNs take on α-helicity (Lisse et al., 1996; Ismail et al., 1999b). We previously examined the binding of DHN1 to liposomes and found that DHNs bind preferentially to anionic phospholipids and that this binding is accompanied by an increase in α-helicity of the protein (Koag et al., 2003). Similarly, a mitochondrial LEA protein, one of the group III LEA proteins, recently has been shown to interact with and protect membranes subjected to desiccation, coupled with the adoption of amphipathic α-helices (Tolleter et al., 2007).Here, we explore the basis of DHN-vesicle interaction using K-segment deletion proteins. This study reveals that the K-segment is necessary and sufficient for binding to anionic phospholipid vesicles and that the adoption of α-helicity of DHN proteins can be attributed mainly to the K-segment.     

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.