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.     

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

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

Leaf Senescence Signaling: The Ca2+-Conducting Arabidopsis Cyclic Nucleotide Gated Channel2 Acts through Nitric Oxide to Repress Senescence Programming

Ca²⁺ and nitric oxide (NO) are essential components involved in plant senescence signaling cascades. In other signaling pathways, NO generation can be dependent on cytosolic Ca²⁺. The Arabidopsis (Arabidopsis thaliana) mutant dnd1 lacks a plasma membrane-localized cation channel (CNGC2). We recently demonstrated that this channel affects plant response to pathogens through a signaling cascade involving Ca²⁺ modulation of NO generation; the pathogen response phenotype of dnd1 can be complemented by application of a NO donor. At present, the interrelationship between Ca²⁺ and NO generation in plant cells during leaf senescence remains unclear. Here, we use dnd1 plants to present genetic evidence consistent with the hypothesis that Ca²⁺ uptake and NO production play pivotal roles in plant leaf senescence. Leaf Ca²⁺ accumulation is reduced in dnd1 leaves compared to the wild type. Early senescence-associated phenotypes (such as loss of chlorophyll, expression level of senescence-associated genes, H₂O₂ generation, lipid peroxidation, tissue necrosis, and increased salicylic acid levels) were more prominent in dnd1 leaves compared to the wild type. Application of a Ca²⁺ channel blocker hastened senescence of detached wild-type leaves maintained in the dark, increasing the rate of chlorophyll loss, expression of a senescence-associated gene, and lipid peroxidation. Pharmacological manipulation of Ca²⁺ signaling provides evidence consistent with genetic studies of the relationship between Ca²⁺ signaling and senescence with the dnd1 mutant. Basal levels of NO in dnd1 leaf tissue were lower than that in leaves of wild-type plants. Application of a NO donor effectively rescues many dnd1 senescence-related phenotypes. Our work demonstrates that the CNGC2 channel is involved in Ca²⁺ uptake during plant development beyond its role in pathogen defense response signaling. Work presented here suggests that this function of CNGC2 may impact downstream basal NO production in addition to its role (also linked to NO signaling) in pathogen defense responses and that this NO generation acts as a negative regulator during plant leaf senescence signaling.Senescence can be considered as the final stage of a plant’s development. During this process, nutrients will be reallocated from older to younger parts of the plant, such as developing leaves and seeds. Leaf senescence has been characterized as a type of programmed cell death (PCD; Gan and Amasino, 1997; Quirino et al., 2000; Lim et al., 2003). During senescence, organelles such as chloroplasts will break down first. Biochemical changes will also occur in the peroxisome during this process. When the chloroplast disassembles, it is easily observed as a loss of chlorophyll. Mitochondria, the source of energy for cells, will be the last cell organelles to undergo changes during the senescence process (Quirino et al., 2000). At the same time, other catabolic events (e.g. protein and lipid breakdown, etc.) are occurring (Quirino et al., 2000). Hormones may also contribute to this process (Gepstein, 2004). From this information we can infer that leaf senescence is regulated by many signals.Darkness treatment can induce senescence in detached leaves (Poovaiah and Leopold, 1973; Chou and Kao, 1992; Weaver and Amasino, 2001; Chrost et al., 2004; Guo and Crawford, 2005; Ülker et al., 2007). Ca²⁺ can delay the senescence of detached leaves (Poovaiah and Leopold, 1973) and leaf senescence induced by methyl jasmonate (Chou and Kao, 1992); the molecular events that mediate this effect of Ca²⁺ are not well characterized at present.Nitric oxide (NO) is a critical signaling molecule involved in many plant physiological processes. Recently, published evidence supports NO acting as a negative regulator during leaf senescence (Guo and Crawford, 2005; Mishina et al., 2007). Abolishing NO generation in either loss-of-function mutants (Guo and Crawford, 2005) or transgenic Arabidopsis (Arabidopsis thaliana) plants expressing NO degrading dioxygenase (NOD; Mishina et al., 2007) leads to an early senescence phenotype in these plants compared to the wild type. Corpas et al. (2004) showed that endogenous NO is mainly accumulated in vascular tissues of pea (Pisum sativum) leaves. This accumulation is significantly reduced in senescing leaves (Corpas et al., 2004). Corpas et al. (2004) also provided evidence that NO synthase (NOS)-like activity (i.e. generation of NO from l-Arg) is greatly reduced in senescing leaves. Plant NOS activity is regulated by Ca²⁺/calmodulin (CaM; Delledonne et al., 1998; Corpas et al., 2004, 2009; del Río et al., 2004; Valderrama et al., 2007; Ma et al., 2008). These studies suggest a link between Ca²⁺ and NO that could be operating during senescence.In animal cells, all three NOS isoforms require Ca²⁺/CaM as a cofactor (Nathan and Xie, 1994; Stuehr, 1999; Alderton et al., 2001). Notably, animal NOS contains a CaM binding domain (Stuehr, 1999). It is unclear whether Ca²⁺/CaM can directly modulate plant NOS or if Ca²⁺/CaM impacts plant leaf development/senescence through (either direct or indirect) effects on NO generation. However, recent studies from our lab suggest that Ca²⁺/CaM acts as an activator of NOS activity in plant innate immune response signaling (Ali et al., 2007; Ma et al., 2008).Although Arabidopsis NO ASSOCIATED PROTEIN1 (AtNOA1; formerly named AtNOS1) was thought to encode a NOS enzyme, no NOS-encoding gene has yet been identified in plants (Guo et al., 2003; Crawford et al., 2006; Zemojtel et al., 2006). However, the AtNOA1 loss-of-function mutant does display reduced levels of NO generation, and several groups have used the NO donor sodium nitroprusside (SNP) to reverse some low-NO related phenotypes in Atnoa1 plants (Guo et al., 2003; Bright et al., 2006; Zhao et al., 2007). Importantly, plant endogenous NO deficiency (Guo and Crawford, 2005; Mishina et al., 2007) or abscisic acid/methyl jasmonate (Hung and Kao, 2003, 2004) induced early senescence can be successfully rescued by application of exogenous NO. Addition of NO donor can delay GA-elicited PCD in barley (Hordeum vulgare) aleurone layers as well (Beligni et al., 2002).It has been suggested that salicylic acid (SA), a critical pathogen defense metabolite, can be increased in natural (Morris et al., 2000; Mishina et al., 2007) and transgenic NOD-induced senescent Arabidopsis leaves (Mishina et al., 2007). Pathogenesis related gene1 (PR1) expression is up-regulated in transgenic Arabidopsis expressing NOD (Mishina et al., 2007) and in leaves of an early senescence mutant (Ülker et al., 2007).Plant cyclic nucleotide gated channels (CNGCs) have been proposed as candidates to conduct extracellular Ca²⁺ into the cytosol (Sunkar et al., 2000; Talke et al., 2003; Lemtiri-Chlieh and Berkowitz, 2004; Ali et al., 2007; Demidchik and Maathuis, 2007; Frietsch et al., 2007; Kaplan et al., 2007; Ma and Berkowitz, 2007; Urquhart et al., 2007; Ma et al., 2009a, 2009b). Arabidopsis “defense, no death” (dnd1) mutant plants have a null mutation in the gene encoding the plasma membrane-localized Ca²⁺-conducting CNGC2 channel. This mutant also displays no hypersensitive response to infection by some pathogens (Clough et al., 2000; Ali et al., 2007). In addition to involvement in pathogen-mediated Ca²⁺ signaling, CNGC2 has been suggested to participate in the process of leaf development/senescence (Köhler et al., 2001). dnd1 mutant plants have high levels of SA and expression of PR1 (Yu et al., 1998), and spontaneous necrotic lesions appear conditionally in dnd1 leaves (Clough et al., 2000; Jirage et al., 2001). Endogenous H₂O₂ levels in dnd1 mutants are increased from wild-type levels (Mateo et al., 2006). Reactive oxygen species molecules, such as H₂O₂, are critical to the PCD/senescence processes of plants (Navabpour et al., 2003; Overmyer et al., 2003; Hung and Kao, 2004; Guo and Crawford, 2005; Zimmermann et al., 2006). Here, we use the dnd1 mutant to evaluate the relationship between leaf Ca²⁺ uptake during plant growth and leaf senescence. Our results identify NO, as affected by leaf Ca²⁺ level, to be an important negative regulator of leaf senescence initiation. Ca²⁺-mediated NO production during leaf development could control senescence-associated gene (SAG) expression and the production of molecules (such as SA and H₂O₂) that act as signals during the initiation of leaf senescence programs.     

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.     

Integrating Image-Based Phenomics and Association Analysis to Dissect the Genetic Architecture of Temporal Salinity Responses in Rice

Salinity affects a significant portion of arable land and is particularly detrimental for irrigated agriculture, which provides one-third of the global food supply. Rice (Oryza sativa), the most important food crop, is salt sensitive. The genetic resources for salt tolerance in rice germplasm exist but are underutilized due to the difficulty in capturing the dynamic nature of physiological responses to salt stress. The genetic basis of these physiological responses is predicted to be polygenic. In an effort to address this challenge, we generated temporal imaging data from 378 diverse rice genotypes across 14 d of 90 mm NaCl stress and developed a statistical model to assess the genetic architecture of dynamic salinity-induced growth responses in rice germplasm. A genomic region on chromosome 3 was strongly associated with the early growth response and was captured using visible range imaging. Fluorescence imaging identified four genomic regions linked to salinity-induced fluorescence responses. A region on chromosome 1 regulates both the fluorescence shift indicative of the longer term ionic stress and the early growth rate decline during salinity stress. We present, to our knowledge, a new approach to capture the dynamic plant responses to its environment and elucidate the genetic basis of these responses using a longitudinal genome-wide association model.Nearly one-third of the 54 million ha of the highly saline soils in the world are located in South and Southeast Asia. Rice (Oryza sativa), which is the primary source of calories and protein for these two regions, is very sensitive to salinity stress, with even moderate salinity levels known to decrease yields by 50% (Zeng et al., 2002). Projected sea level rise due to climate change is expected to increase saltwater ingress in coastal rice-growing regions of South and Southeast Asia. Therefore, development of salt-tolerant rice cultivars is essential to maintain rice productivity in the salinity-affected regions globally.Salt tolerance, defined as the ability to maintain growth and productivity in saline conditions, is a complex polygenic trait that may be influenced by distinct physiological mechanisms (Munns et al., 1982; Munns and Termaat, 1986; Cheeseman, 1988; Munns and Tester, 2008; Horie et al., 2012; for a comprehensive review of genes involved in salinity tolerance in rice, see Negrão et al., 2011) At the cellular level, plants respond to saline conditions in two phases, namely an osmotic (shoot ion independent) and an ionic stress phase, which can occur in an overlapping manner with varying intensity during the course of salinity stress (Munns and Termaat, 1986; Munns, 2002; Munns and James, 2003; Munns and Tester, 2008; Horie et al., 2012). During the osmotic stress phase, which occurs soon after the onset of salinity, the reduction in external osmotic potential disrupts water uptake and impedes cell expansion, which, at the whole plant level, leads to reduced growth rate (Matsuda and Riazi, 1981; Munns and Passioura, 1984; Rawson and Munns, 1984; Azaizeh and Steudle, 1991; Fricke and Peters, 2002; Fricke, 2004; Boursiac et al., 2005). As salinity stress persists over several days and weeks, sodium ions (Na⁺) accumulate to toxic levels, resulting in cell death and precocious leaf senescence (Lutts and Bouharmont, 1996; Munns, 2002; Munns and James, 2003; Ghanem et al., 2008). This is typically observed during the ionic phase of the salinity response (Munns, 2002; Munns and James, 2003; Munns and Tester, 2008). Plants possess distinct mechanisms to adapt to these osmotic and ionic stresses that are controlled by a suite of genes (Apse et al., 1999; Carvajal et al., 1999; Halfter et al., 2000; Ishitani et al., 2000; Shi et al., 2000; Zeng and Shannon, 2000; Rus et al., 2001; Berthomieu et al., 2003; Martínez-Ballesta et al., 2003; Boursiac et al., 2005, 2008; Ren et al., 2005; Huang et al., 2006; Davenport et al., 2007; Obata et al., 2007; Székely et al., 2008; Horie et al., 2011; Rivandi et al., 2011; Asano et al., 2012; Munns et al., 2012; Latz et al., 2013; Schmidt et al., 2013; Campo et al., 2014; Choi et al., 2014; Liu et al., 2014). The genetic basis of temporal adaptive responses to salinity stress remains to be explored in rice and other crops. This is primarily due to challenges in capturing the dynamic physiological responses to salinity for a large number of genotypes in a nondestructive manner. Manual phenotyping to detect incremental changes in growth rate during the osmotic stress and slight shifts in leaf color due to ionic stress is difficult to quantify for a large number of genotypes.In rice, at least one major quantitative trait loci (QTL; saltol) for salinity tolerance has been characterized based on end point measurements of biomass, senescence/injury, and Na⁺ and K⁺ concentrations (Bonilla et al., 2002; Lin et al., 2004; Thomson et al., 2010). SHOOT K⁺ CONTENT1 (SKC1) is the causative gene underlying the saltol region. SKC1 encodes a Na⁺-selective high-affinity potassium transporter that regulates Na⁺/K⁺ homeostasis during salinity stress (Ren et al., 2005). High levels of Na⁺ displace cellular K⁺, an essential element for several enzymatic reactions and physiological processes (Gierth and Mäser, 2007). The ability to maintain cellular K⁺ levels during salinity through the action of Na⁺-selective potassium transporters or Na⁺/H⁺ antiporters is a well-characterized tolerance mechanism in cereals including rice (Ren et al., 2005; Sunarpi et al., 2005; Huang et al., 2006; Møller et al., 2009; Mian et al., 2011; Munns et al., 2012).Numerous studies have utilized conventional linkage mapping to identify QTL for morphological and physiological responses to salinity in rice using discrete end point measurements (Bonilla et al., 2002; Lin et al., 2004; Ren et al., 2005; Negrão et al., 2011; Wang et al., 2012). However, the physiological adaptation to saline conditions is a complex continuous process that develops over time. While some accessions will exhibit similar end point phenotypic values, the genetic and physiological mechanisms giving rise to the similar phenotypes may be very different and the growth trajectories throughout the experiment may be distinct. Although single time point studies have yielded important information regarding the genetic basis of salinity tolerance, such approaches are too simple to reveal the genetic architecture of stress adaptation. With the advent of high-throughput image-based phenotyping platforms, it is now feasible to quantify dynamic responses during the stress treatment for a large number of genotypes (Berger et al., 2010; Golzarian et al., 2011; Chen et al., 2014; Honsdorf et al., 2014).Image-based phenotyping has been combined with genome-wide association studies (GWAS) and linkage mapping to examine the genetic basis of complex developmental processes (Busemeyer et al., 2013; Moore et al., 2013; Topp et al., 2013; Slovak et al., 2014; Würschum et al., 2014; Yang et al., 2014; Bac-Molenaar et al., 2015). Moreover, the introduction of the time axis provides a better understanding of the physiological processes underlying complex stress and developmental responses compared with single end point measurements (Zhang et al., 2012; Moore et al., 2013; Brown et al., 2014; Chen et al., 2014; Slovak et al., 2014; Bac-Molenaar et al., 2015). However, to date, no studies have implemented an association mapping approach using image-derived phenotypes to address the genetic basis of dynamic stress responses in plants. Image-based phenotyping offers several advantages over conventional phenotyping: (1) quantitative measurements can be recorded over discrete time points to capture morphological and physiological responses in a nondestructive manner, and (2) the use of various types of spectral imaging address phenotypes that are not detectable to the human eye such as chlorophyll fluorescence and leaf water content. Integrating dynamic phenotypic data and association mapping has the potential to query genetic diversity across hundreds of accessions for complex traits and provides much higher resolution compared with conventional linkage mapping. Here, we explored the dynamic growth and chlorophyll responses to salinity of a diverse set of rice accessions using high-throughput visible and fluorescence imaging. To assess the genetic basis of plant growth in saline conditions, a logistic model was used to describe the temporal growth responses and was incorporated into the statistical framework necessary for association mapping. Coupled with temporal fluorescence imaging, we present, to our knowledge, new insights into the genetic architecture of osmotic and ionic responses during salinity stress in rice.     

Enhanced Photosynthesis and Growth in atquac1 Knockout Mutants Are Due to Altered Organic Acid Accumulation and an Increase in Both Stomatal and Mesophyll Conductance

Stomata control the exchange of CO₂ and water vapor in land plants. Thus, whereas a constant supply of CO₂ is required to maintain adequate rates of photosynthesis, the accompanying water losses must be tightly regulated to prevent dehydration and undesired metabolic changes. Accordingly, the uptake or release of ions and metabolites from guard cells is necessary to achieve normal stomatal function. The AtQUAC1, an R-type anion channel responsible for the release of malate from guard cells, is essential for efficient stomatal closure. Here, we demonstrate that mutant plants lacking AtQUAC1 accumulated higher levels of malate and fumarate. These mutant plants not only display slower stomatal closure in response to increased CO₂ concentration and dark but are also characterized by improved mesophyll conductance. These responses were accompanied by increases in both photosynthesis and respiration rates, without affecting the activity of photosynthetic and respiratory enzymes and the expression of other transporter genes in guard cells, which ultimately led to improved growth. Collectively, our results highlight that the transport of organic acids plays a key role in plant cell metabolism and demonstrate that AtQUAC1 reduce diffusive limitations to photosynthesis, which, at least partially, explain the observed increments in growth under well-watered conditions.Stomata are functionally specialized microscopic pores that control the essential exchange of CO₂ and H₂O with the environment in land plants. Stomata are found on the surfaces of the majority of the aerial parts of plants, rendering them as the main control point regulating the flow of gases between plants and their surrounding atmosphere. Accordingly, the majority of water loss from plants occurs through stomatal pores, allowing plant transpiration and CO₂ absorption for the photosynthetic process (Bergmann and Sack, 2007; Kim et al., 2010). The maintenance of an adequate water balance through stomatal control is crucial to plants because cell expansion and growth require tissues to remain turgid (Sablowski and Carnier Dornelas, 2014), and minor reductions in cell water volume and turgor pressure will therefore compromise both processes (Thompson, 2005). As a result, the high sensitivity of plant tissues to turgor has prompted the use of reverse genetic studies in attempt to engineer plants with improved performance (Cowan and Troughton, 1971; Xiong et al., 2009; Borland et al., 2014; Franks et al., 2015).In most land plants, not only redox signals invoked by shifts in light quality (Busch, 2014) but also the transport of inorganic ions (e.g. K⁺, Cl⁻, and NO⁻₃) as well as metabolites such as the phytohormone abscisic acid (ABA), Suc, and malate, are important players controlling stomatal movements (Hetherington, 2001; Roelfsema and Hedrich, 2005; Pandey et al., 2007; Blatt et al., 2014; Kollist et al., 2014). In this context, although organic acids in plants is known to support numerous and diverse functions both within and beyond cellular metabolism, only recently have we obtained genetic evidence to support that modulation of guard cell malate and fumarate concentration can greatly influence stomatal movements (Nunes-Nesi et al., 2007; Araújo et al., 2011b; Penfield et al., 2012; Medeiros et al., 2015). Notably malate, in particular, has been considered as a key metabolite and one of the most important organic metabolites involved in guard cell movements (Hedrich and Marten, 1993; Fernie and Martinoia, 2009; Meyer et al., 2010). During stomatal aperture, the flux of malate into guard cells coupled with hexoses generated on starch breakdown lead to decreases in the water potential, and consequently, water uptake by the guard cells ultimately opens the stomata pore (Roelfsema and Hedrich, 2005; Vavasseur and Raghavendra, 2005; Lee et al., 2008). On the other hand, during stomatal closure, malate is believed to be converted into starch, which has no osmotic activity (Penfield et al., 2012) or, alternatively, is released from the guard cells to the surrounding apoplastic space (Lee et al., 2008; Negi et al., 2008; Vahisalu et al., 2008; Meyer et al., 2010).The role of organic acids on the stomatal movements has been largely demonstrated by studies related to malate transport (Lee et al., 2008; Meyer et al., 2010; Sasaki et al., 2010). In the last decade, two protein families were identified and functionally characterized to be directly involved with organic acid transport at the guard cell plasma membrane and to be required for stomatal functioning (Lee et al., 2008; Meyer et al., 2010; Sasaki et al., 2010). In summary, AtABCB14, a member of the ABC (ATP binding cassette) family, which is involved in malate transport from apoplast to guard cells, was described as a negative modulator of stomatal closure induced by high CO₂ concentration; notably, exogenous application of malate minimizes this response (Lee et al., 2008). In addition, members of a small gene family, which encode the anion channels SLAC1 (slow anion channel 1) and four SLAC1-homologs (SLAHs) in Arabidopsis (Arabidopsis thaliana), have been described to be involved in stomatal movements. SLAC1 is a well-documented S-type anion channel that preferentially transports chloride and nitrate as opposed to malate (Vahisalu et al., 2008, 2010; Geiger et al., 2010; Du et al., 2011; Brandt et al., 2012; Kusumi et al., 2012). Lack of SLAC1 in Arabidopsis and rice (Oryza sativa) culminated in a failure in stomatal closure in response to high CO₂ levels, low relative humidity, and dark conditions (Negi et al., 2008; Vahisalu et al., 2008; Kusumi et al., 2012). Although mutations in AtSLAC1 impair S-type anion channel functions as a whole, the R-type anion channel remained functional (Vahisalu et al., 2008). Indeed, a member of the aluminum-activated malate transporter (ALMT) family, AtALMT12, an R-type anion channel, has been demonstrated to be involved in malate transport, particularly at the plasma membrane of guard cells (Meyer et al., 2010; Sasaki et al., 2010). Although AtALMT12 is a member of ALMT family, it is not activated by aluminum, and therefore Meyer et al. (2010) proposed to rename it as AtQUAC1 (quick-activating anion channel 1; Imes et al., 2013; Mumm et al., 2013). Hereafter, we will follow this nomenclature. Deficiency of a functional AtQUAC1 has been documented to lead to changes in stomatal closure in response to high levels of CO₂, dark, and ABA (Meyer et al., 2010). Taken together, these studies have clearly demonstrated that both S- and R-type anion channels are key modulators of stomatal movements in response to several environmental factors.Despite a vast number of studies involving the above-mentioned anion channels, little information concerning the metabolic changes caused by their impairment is currently available. Such information is important to understand stomatal movements, mainly considering that organic acids, especially the levels of malate in apoplastic/mesophyll cells, have been highlighted as of key importance in leaf metabolism (Fernie and Martinoia, 2009; Araújo et al., 2011a, 2011b; Lawson et al., 2014; Medeiros et al., 2015). Here, we demonstrate that a disruption in the expression of AtQUAC1, which leads to impaired stomatal closure (Meyer et al., 2010), was accompanied by increases in mesophyll conductance (g_m), which is defined as the conductance for the transfer of CO₂ from the intercellular airspaces (C_i) to the sites of carboxylation in the chloroplastic stroma (C_c). By further characterization of atquac1 knockout plants, we demonstrated that reduced diffusive limitations resulted in higher photosynthetic rates and altered respiration that, in turn, led to enhanced biomass accumulation. Overall, the results obtained are discussed both in terms of the importance of organic acid transport in plant cell metabolism and with regard to the contribution that it plays in the regulation of both stomatal function and growth.     

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.     

CPK13, a Noncanonical Ca2+-Dependent Protein Kinase,Specifically Inhibits KAT2 and KAT1 Shaker K+ Channels and Reduces Stomatal Opening

Ca²⁺-dependent protein kinases (CPKs) form a large family of 34 genes in Arabidopsis (Arabidopsis thaliana). Based on their dependence on Ca²⁺, CPKs can be sorted into three types: strictly Ca²⁺-dependent CPKs, Ca²⁺-stimulated CPKs (with a significant basal activity in the absence of Ca²⁺), and essentially calcium-insensitive CPKs. Here, we report on the third type of CPK, CPK13, which is expressed in guard cells but whose role is still unknown. We confirm the expression of CPK13 in Arabidopsis guard cells, and we show that its overexpression inhibits light-induced stomatal opening. We combine several approaches to identify a guard cell-expressed target. We provide evidence that CPK13 (1) specifically phosphorylates peptide arrays featuring Arabidopsis K⁺ Channel KAT2 and KAT1 polypeptides, (2) inhibits KAT2 and/or KAT1 when expressed in Xenopus laevis oocytes, and (3) closely interacts in plant cells with KAT2 channels (Förster resonance energy transfer-fluorescence lifetime imaging microscopy). We propose that CPK13 reduces stomatal aperture through its inhibition of the guard cell-expressed KAT2 and KAT1 channels.Stomata are microscopic organs at the leaf surface, each made of two so-called guard cells forming a pore. Opening or closing these pores is the way through which plants control their gas exchanges with the atmosphere (i.e. carbon dioxide uptake to feed the photosynthetic process and transpirational loss of water vapor). Stomatal movements result from osmotically driven fluxes of water, which follow massive exchanges of solutes, including K⁺ ions, between the guard cells and the surrounding tissues (Hetherington, 2001; Nilson and Assmann, 2007).Both Ca²⁺-dependent and Ca²⁺-independent signaling pathways are known to control stomatal movements (MacRobbie, 1993, 1998; Blatt, 2000; Webb et al., 2001; Mustilli et al., 2002; Israelsson et al., 2006; Marten et al., 2007; Laanemets et al., 2013). In particular, Ca²⁺ signals have been reported to promote stomatal closure through the inhibition of inward K⁺ channels and the activation of anion channels (Blatt, 1991, 1992, 2000; Thiel et al., 1992; Grabov and Blatt, 1999; Schroeder et al., 2001; Hetherington and Brownlee, 2004; Mori et al., 2006; Marten et al., 2007; Geiger et al., 2010; Brandt et al., 2012; Scherzer et al., 2012). However, little is known about the molecular identity of the links between Ca²⁺ events and Shaker K⁺ channel activity. Several kinases and phosphatases are believed to be involved in both the Ca²⁺-dependent and Ca²⁺-independent signaling pathways. Plants express two large kinase families whose activity is related to Ca²⁺ signaling. Firstly, CBL-interacting protein kinases (CIPKs; 25 genes in Arabidopsis [Arabidopsis thaliana]) are indirectly controlled by their interaction with a set of calcium sensors, the calcineurin B-like proteins (CBLs; 10 genes in Arabidopsis). This complex forms a fascinating network of potential Ca²⁺ signaling decoders (Luan, 2009; Weinl and Kudla, 2009), which have been addressed in numerous reports (Xu et al., 2006; Hu et al., 2009; Batistic et al., 2010; Held et al., 2011; Chen et al., 2013). In particular, some CBL-CIPK pairs have been shown to regulate Shaker channels such as Arabidopsis K⁺ Transporter1 (AKT1; Xu et al., 2006; Lan et al., 2011) or AKT2 (Held et al., 2011). Second, Ca²⁺-dependent protein kinases (CPKs) form an even larger family (34 genes in Arabidopsis) of proteins combining a kinase domain with the ability to bind Ca²⁺, thanks to the so-called EF hands (Harmon et al., 2000; Harper et al., 2004). CPKs, which, interestingly, are not found in animal cells, exhibit different calcium dependencies (Boudsocq et al., 2012). With respect to this, three types of CPKs can be considered: strictly Ca²⁺-dependent CPKs, Ca²⁺-stimulated CPKs (with a significant basal activity in the absence of Ca²⁺), and essentially Ca²⁺-insensitive CPKs (however, structurally close to kinases of groups 1 and 2).Pioneering work by Luan et al. (1993) demonstrated in Vicia faba guard cells that inward K⁺ channels were regulated by some Ca²⁺-dependent kinases. Then, such a Ca²⁺-dependent kinase was purified from guard cell protoplasts of V. faba and shown to actually phosphorylate the in vitro-translated KAT1 protein, a Shaker channel subunit natively expressed in Arabidopsis guard cells (Li et al., 1998). KAT1 regulation by CPK was shown by the inhibition of KAT1 currents after the coexpression of KAT1 and CDPK from soybean (Glycine max) in oocytes (Berkowitz et al., 2000). Since then, several cpk mutant lines of Arabidopsis have been shown to be impaired in stomatal movements, for example cpk10 (Ca²⁺ insensitive), cpk4/cpk11 (Ca²⁺ dependent), and cpk3/cpk6/cpk23 (Ca²⁺ dependent; Mori et al., 2006; Geiger et al., 2010; Munemasa et al., 2011; Hubbard et al., 2012).Of the nine genes encoding voltage-dependent K⁺ channels (Shaker) in Arabidopsis (Véry and Sentenac, 2002, 2003; Lebaudy et al., 2007; Hedrich, 2012), six are expressed in guard cells and play a role in stomatal movements: the Gated Outwardly-Rectifying K⁺ (GORK) gene, encoding an outward K⁺ channel subunit, and the AKT1, AKT2, Arabidopsis K⁺ Rectifying Channel1 (AtKC1), KAT1, and KAT2 genes, encoding inward K⁺ channel subunits (Pilot et al., 2001; Szyroki et al., 2001; Hosy et al., 2003; Pandey et al., 2007; Lebaudy et al., 2008a). Shaker channels result from the assembly of four subunits, and it has been shown that inward subunits tend to heterotetramerize, thus potentially widening the functional and regulatory scope of inward K⁺ conductance in guard cells (Xicluna et al., 2007; Jeanguenin et al., 2008; Lebaudy et al., 2008a, 2010). Inhibition of inward K⁺ channels has been shown to reduce stomatal opening (Liu et al., 2000; Kwak et al., 2001). This has grounded a strategy for disrupting inward K⁺ channel conductance in guard cells by expressing a nonfunctional KAT2 subunit (dominant negative mutation) in a kat2 knockout Arabidopsis line. The resulting Arabidopsis lines, named kincless, have no functional inward K⁺ channels and exhibit delayed stomatal opening (Lebaudy et al., 2008b) with, in the long term, a biomass reduction compared with the Arabidopsis wild-type line.Among the CPKs presumably expressed in Arabidopsis guard cells (Leonhardt et al., 2004), we looked for CPK13, which belongs to the atypical Ca²⁺-insensitive type of CPKs (Kanchiswamy et al., 2010; Boudsocq et al., 2012; Liese and Romeis, 2013) and whose role remains unknown in stomatal movements. Here, we confirm first that CPK13 kinase activity is independent of Ca²⁺ and show that CPK13 expression is predominant in Arabidopsis guard cells using CPK13-GUS lines. We then report that overexpression of CPK13 in Arabidopsis induces a dramatic default in stomatal aperture. Based on the previously reported kincless phenotype (Lebaudy et al., 2008b), we propose that CPK13 could reduce the activity of inward K⁺ channels in guard cells, particularly that of KAT2. We confirm this hypothesis by voltage-clamp experiments and show an inhibition of KAT2 and KAT1 activity by CPK13 (but not that of AKT2). In addition, we present peptide array phosphorylation assays showing that CPK13 targets, with some specificity, several KAT2 and KAT1 polypeptides. Finally, we demonstrate that KAT2 and CPK13 interact in planta using Förster resonance energy transfer (FRET)-fluorescence lifetime imaging microscopy (FLIM).     

Cavitation Resistance in Seedless Vascular Plants: The Structure and Function of Interconduit Pit Membranes

Plant water transport occurs through interconnected xylem conduits that are separated by partially digested regions in the cell wall known as pit membranes. These structures have a dual function. Their porous construction facilitates water movement between conduits while limiting the spread of air that may enter the conduits and render them dysfunctional during a drought. Pit membranes have been well studied in woody plants, but very little is known about their function in more ancient lineages such as seedless vascular plants. Here, we examine the relationships between conduit air seeding, pit hydraulic resistance, and pit anatomy in 10 species of ferns (pteridophytes) and two lycophytes. Air seeding pressures ranged from 0.8 ± 0.15 MPa (mean ± sd) in the hydric fern Athyrium filix-femina to 4.9 ± 0.94 MPa in Psilotum nudum, an epiphytic species. Notably, a positive correlation was found between conduit pit area and vulnerability to air seeding, suggesting that the rare-pit hypothesis explains air seeding in early-diverging lineages much as it does in many angiosperms. Pit area resistance was variable but averaged 54.6 MPa s m⁻¹ across all surveyed pteridophytes. End walls contributed 52% to the overall transport resistance, similar to the 56% in angiosperm vessels and 64% in conifer tracheids. Taken together, our data imply that, irrespective of phylogenetic placement, selection acted on transport efficiency in seedless vascular plants and woody plants in equal measure by compensating for shorter conduits in tracheid-bearing plants with more permeable pit membranes.Water transport in plants occurs under tension, which renders the xylem susceptible to air entry. This air seeding may lead to the rupture of water columns (cavitation) such that the air expands within conduits to create air-vapor embolisms that block further transport. (Zimmermann and Tyree, 2002). Excessive embolism such as that which occurs during a drought may jeopardize leaf hydration and lead to stomatal closure, overheating, wilting, and possibly death of the plant (Hubbard et al., 2001; Choat et al., 2012; Schymanski et al., 2013). Consequently, strong selection pressure resulted in compartmentalized and redundant plant vascular networks that are adapted to a species habitat water availability by way of life history strategy (i.e. phenology) or resistance to air seeding (Tyree et al., 1994; Mencuccini et al., 2010; Brodersen et al., 2012). The spread of drought-induced embolism is limited primarily by pit membranes, which are permeable, mesh-like regions in the primary cell wall that connect two adjacent conduits. The construction of the pit membrane is such that water easily moves across the membrane between conduits, but because of the small membrane pore size and the presence of a surface coating on the membrane (Pesacreta et al., 2005; Lee et al., 2012), the spread of air and gas bubbles is restricted up to a certain pressure threshold known as the air-seeding pressure (ASP). When xylem sap tension exceeds the air-seeding threshold, air can be aspirated from an air-filled conduit into a functional water-filled conduit through perhaps a large, preexisting pore or one that is created by tension-induced membrane stress (Rockwell et al., 2014). Air seeding leads to cavitation and embolism formation, with emboli potentially propagating throughout the xylem network (Tyree and Sperry, 1988; Brodersen et al., 2013). So, on the one hand, pit membranes are critical to controlling the spread of air throughout the vascular network, while on the other hand, they must facilitate the efficient flow of water between conduits (Choat et al., 2008; Domec et al., 2008; Pittermann et al., 2010; Schulte, 2012). Much is known about such hydraulic tradeoffs in the pit membranes of woody plants, but comparatively little data exist on seedless vascular plants such as ferns and lycophytes. Given that seedless vascular plants may bridge the evolutionary transition from bryophytes to woody plants, the lack of functional data on pit membrane structure in early-derived tracheophytes is a major gap in our understanding of the evolution of plant water transport.In woody plants, pit membranes fall into one of two categories: the torus-margo type found in most gymnosperms and the homogenous pit membrane characteristic of angiosperms (Choat et al., 2008; Choat and Pittermann, 2009). In conifers, water moves from one tracheid to another through the margo region of the membrane, with the torus sealing the pit aperture should one conduit become embolized. Air seeding occurs when water potential in the functional conduit drops low enough to dislodge the torus from its sealing position, letting air pass through the pit aperture into the water-filled tracheid (Domec et al., 2006; Delzon et al., 2010; Pittermann et al., 2010; Schulte, 2012; but see Jansen et al., 2012). Across north-temperate conifer species, larger pit apertures correlate with lower pit resistance to water flow (r_pit; MPa s m⁻¹), but it is the ratio of torus-aperture overlap that sets a species cavitation resistance (Pittermann et al., 2006, 2010; Domec et al., 2008; Hacke and Jansen, 2009). A similar though mechanistically different tradeoff exists in angiosperm pit membranes. Here, air seeding reflects a probabilistic relationship between membrane porosity and the total area of pit membranes present in the vessel walls. Specifically, the likelihood of air aspirating into a functional conduit is determined by the combination of xylem water potential and the diameter of the largest pore and/or the weakest zone in the cellulose matrix in the vessel’s array of pit membranes (Wheeler et al., 2005; Hacke et al., 2006; Christman et al., 2009; Rockwell et al., 2014). As it has come to be known, the rare-pit hypothesis suggests that the infrequent, large-diameter leaky pore giving rise to that rare pit reflects some combination of pit membrane traits such as variation in conduit membrane area (large or small), membrane properties (tight or porous), and hydrogel membrane chemistry (Hargrave et al., 1994; Choat et al., 2003; Wheeler et al., 2005; Hacke et al., 2006; Christman et al., 2009; Lee et al., 2012; Plavcová et al., 2013; Rockwell et al., 2014). The maximum pore size is critical because, per the Young-Laplace law, the larger the radius of curvature, the lower the air-water pressure difference under which the contained meniscus will fail (Jarbeau et al., 1995; Choat et al., 2003; Jansen et al., 2009). Consequently, angiosperms adapted to drier habitats may exhibit thicker, denser, smaller, and less abundant pit membranes than plants occupying regions with higher water availability (Wheeler et al., 2005; Hacke et al., 2007; Jansen et al., 2009; Lens et al., 2011; Scholz et al., 2013). However, despite these qualitative observations, there is no evidence that increased cavitation resistance arrives at the cost of higher r_pit. Indeed, the bulk of the data suggest that prevailing pit membrane porosity is decoupled from the presence of the single largest pore that allows air seeding to occur (Choat et al., 2003; Wheeler et al., 2005 Hacke et al., 2006, 2007).As water moves from one conduit to another, pit membranes offer considerable hydraulic resistance throughout the xylem network. On average, r_pit contributes 64% and 56% to transport resistance in conifers and angiosperms, respectively (Wheeler et al., 2005; Pittermann et al., 2006; Sperry et al., 2006). In conifers, the average r_pit is estimated at 6 ± 1 MPa s m⁻¹, almost 60 times lower than the 336 ± 81 MPa s m⁻¹ computed for angiosperms (Wheeler et al., 2005; Hacke et al., 2006; Sperry et al., 2006). Presumably, the high porosity of conifer pits compensates for the higher transport resistance offered by a vascular system composed of narrow, short, single-celled conduits (Pittermann et al., 2005; Sperry et al., 2006).Transport in seedless vascular plants presents an interesting conundrum because, with the exception of a handful of species, their primary xylem is composed of tracheids, the walls of which are occupied by homogenous pit membranes (Gibson et al., 1985; Carlquist and Schneider, 2001, 2007; but see Morrow and Dute, 1998, for torus-margo membranes in Botrychium spp.). At first pass, this combination of traits appears hydraulically maladaptive, but several studies have shown that ferns can exhibit transport capacities that are on par with more recently evolved plants (Wheeler et al., 2005; Watkins et al., 2010; Pittermann et al., 2011, 2013; Brodersen et al., 2012). Certainly, several taxa possess large-diameter, highly overlapping conduits, some even have vessels such as Pteridium aquilinum and many species have high conduit density, all of which could contribute to increased hydraulic efficiency (Wheeler et al., 2005; Pittermann et al., 2011, 2013). But how do the pit membranes of seedless vascular plants compare? Scanning electron micrographs of fern and lycopod xylem conduits suggest that they are thin, diaphanous, and susceptible to damage during specimen preparation (Carlquist and Schneider 2001, 2007). Consistent with such observations, two estimates of r_pit imply that r_pit in ferns may be significantly lower than in angiosperms; Wheeler et al. (2005) calculated r_pit in the fern Pteridium aquilinum at 31 MPa s m⁻¹, while Schulte et al. (1987) estimated r_pit at 1.99 MPa s m⁻¹ in the basal fern Psilotum nudum. The closest structural analogy to seedless vascular plant tracheids can be found in the secondary xylem of the early-derived vesselless angiosperms, in which tracheids possess homogenous pit membranes with r_pit values that at 16 MPa s m⁻¹ are marginally higher than those of conifers (Hacke et al., 2007). Given that xylem in seedless vascular plants is functionally similar to that in vesselless angiosperms, we expected convergent r_pit values in these two groups despite their phylogenetic distance. We tested this hypothesis, as well as the intrinsic cavitation resistance of conduits in seedless vascular plants, by scrutinizing the pit membranes of ferns and fern allies using the anatomical and experimental approaches applied previously to woody taxa. In particular, we focused on the relationship between pit membrane traits and cavitation resistance at the level of the individual conduit.