Overexpression of Flavodoxin in Bacteroids Induces Changes in Antioxidant Metabolism Leading to Delayed Senescence and Starch Accumulation in Alfalfa Root Nodules

Sinorhizobium meliloti cells were engineered to overexpress Anabaena variabilis flavodoxin, a protein that is involved in the response to oxidative stress. Nodule natural senescence was characterized in alfalfa (Medicago sativa) plants nodulated by the flavodoxin-overexpressing rhizobia or the corresponding control bacteria. The decline of nitrogenase activity and the nodule structural and ultrastructural alterations that are associated with nodule senescence were significantly delayed in flavodoxin-expressing nodules. Substantial changes in nodule antioxidant metabolism, involving antioxidant enzymes and ascorbate-glutathione cycle enzymes and metabolites, were detected in flavodoxin-containing nodules. Lipid peroxidation was also significantly lower in flavodoxin-expressing nodules than in control nodules. The observed amelioration of the oxidative balance suggests that the delay in nodule senescence was most likely due to a role of the protein in reactive oxygen species detoxification. Flavodoxin overexpression also led to high starch accumulation in nodules, without reduction of the nitrogen-fixing activity.Symbiotic nodules have a limited functional life that varies among different legume species. Nodule senescence is the sequence of structural, molecular, biochemical, and physiological events taking place in the process that a mature and functional nodule undergoes leading to the loss of the nitrogen-fixing activity and culminating in cell death of symbiotic tissue (Swaraj and Bishnoi, 1996; Puppo et al., 2005; Van de Velde et al., 2006).Various models have been proposed to explain the mechanisms that trigger the process of natural or stress-induced nodule senescence. However, it is generally accepted that a senescence-inducing signal from the plant causes a decrease in antioxidant levels and thus an increase in reactive oxygen species (ROS) up to a point of no return. Numerous studies have shown that ROS and antioxidant systems are involved in natural (Lucas et al., 1998; Evans et al., 1999; Hernández-Jiménez et al., 2002; Puppo et al., 2005) as well as induced (Dalton et al., 1993; Becana et al., 2000; Hernández-Jiménez et al., 2002; Matamoros et al., 2003) nodule senescence. Nitrogen fixation is very sensitive to ROS, and nitrogenase activity drastically decreases during nodule senescence (Dalton et al., 1986).Antioxidant systems that protect cells from oxidative damage have been described in symbiotic nodules (Dalton et al., 1986, 1993; Evans et al., 1999; Becana et al., 2000; Matamoros et al., 2003; Puppo et al., 2005). These include the enzymes superoxide dismutase (SOD), catalase, and peroxidase. Another enzymatic system associated with ROS detoxification is the ascorbate-glutathione pathway, which includes ascorbate peroxidase (APX), dehydroascorbate reductase (DHAR), monodehydroascorbate reductase (MDHAR), and glutathione reductase (GR; Dalton et al., 1986, 1992; Noctor and Foyer 1998; Becana et al., 2000). Ascorbate and reduced glutathione (GSH) in this pathway can also scavenge superoxide and hydrogen peroxide.During nodule senescence, several ultrastructural alterations in the nodule tissues and cells have been observed (Lucas et al., 1998; Hernández-Jiménez et al., 2002; Puppo et al., 2005, and refs. therein; Van de Velde et al., 2006). Cytosol becomes electron dense, altered vesicles proliferate, and eventually the cytosol undergoes lysis. The number of peroxisomes increases, mitochondria form complex elongated structures, and symbiosomes change in size and shape and fuse during natural and induced senescence of nodules (Hernández-Jiménez et al., 2002). Damage of the symbiosome membrane is also detected (Puppo et al., 2005; Van de Velde et al., 2006).A strategy of delayed nodule senescence could lead to increased nitrogen fixation and legume productivity. Delayed nodule senescence together with enhanced sustainability under field conditions are among the key aims of legume improvement programs (Puppo et al., 2005). An interesting approach proposed to achieve delayed senescence is to induce nodulation in legumes using rhizobial strains with modified redox capacity (Zahran, 2001).The protein flavodoxin contains a FMN group acting as a redox center transferring electrons at low potentials (Pueyo et al., 1991; Pueyo and Gómez-Moreno, 1991). The FMN cofactor of flavodoxin can exist in three different redox states: oxidized, one-electron-reduced semiquinone, and two-electron-reduced hydroquinone. This property confers high versatility to flavodoxins in electron transport systems (Simondsen and Tollin, 1980; McIver et al., 1998). To date, flavodoxin has not been described in plants, as flavodoxin-encoding genes were lost during the transition of algae to plants (Zurbriggen et al., 2007) and, consequently, no homologs have been identified in the sequenced genome of Arabidopsis (Arabidopsis thaliana; Arabidopsis Genome Initiative, 2000). Flavodoxin is present as a constitutive or inducible protein in different microorganisms (Klugkist et al., 1986). In the nitrogen-fixing cyanobacterium Anabaena variabilis PCC 7119, flavodoxin is expressed under conditions of limited iron availability, replacing ferredoxin in the photosynthetic electron transport from PSI to NADP⁺ and in nitrogenase reduction (Sandmann et al., 1990). Reversible electron transfer from flavodoxin to NADP⁺ is catalyzed by ferredoxin NADP⁺ reductase in different pathways of oxidative metabolism (Arakaki et al., 1997). In its reduced state, flavodoxin might be able to react with ROS and revert to its original redox state in the presence of an appropriate electron source. This could probably occur without the associated molecular damage that metallic complexes in catalases or SODs suffer (Keyer et al., 1995). The presence of flavodoxin has not been documented to date in the symbiotic bacterium Sinorhizobium meliloti. In Escherichia coli, however, flavodoxin induction is linked to the oxidative stress-responsive regulon soxRS (Zheng et al., 1999). It has been suggested that flavodoxin and ferredoxin (flavodoxin) NADP⁺ reductase might be induced and have a role in reestablishing the cell redox balance under oxidative stress conditions (Liochev et al., 1994). The properties of flavodoxin suggest that its presence in the cell may have a facilitating effect on ROS detoxification. In fact, an increase in the amount of flavodoxin has been observed in some bacterial species subjected to oxidative stress (Zheng et al., 1999; Yousef et al., 2003; Singh et al., 2004), and transgenic tobacco (Nicotiana tabacum) plants expressing flavodoxin in chloroplasts show enhanced tolerance to a broad range of stresses related to oxidative damage (Tognetti et al., 2006, 2007a, 2007b).In this work, Sinorhizobium meliloti was transformed with the A. variabilis flavodoxin gene and used to nodulate alfalfa (Medicago sativa) plants. The effects of flavodoxin expression on nodulation dynamics, on nodule development and senescence processes, and on nitrogen-fixing activity were analyzed. Mechanistic insights suggesting putative roles for flavodoxin in protection from ROS and the induced delay of nodule senescence are likewise discussed.     

Differential Impact of Lipoxygenase 2 and Jasmonates on Natural and Stress-Induced Senescence in Arabidopsis

Jasmonic acid and related oxylipins are controversially discussed to be involved in regulating the initiation and progression of leaf senescence. To this end, we analyzed profiles of free and esterified oxylipins during natural senescence and upon induction of senescence-like phenotypes by dark treatment and flotation on sorbitol in Arabidopsis (Arabidopsis thaliana). Jasmonic acid and free 12-oxo-phytodienoic acid increased during all three processes, with the strongest increase of jasmonic acid after dark treatment. Arabidopside content only increased considerably in response to sorbitol treatment. Monogalactosyldiacylglycerols and digalactosyldiacylglycerols decreased during these treatments and aging. Lipoxygenase 2-RNA interference (RNAi) plants were generated, which constitutively produce jasmonic acid and 12-oxo-phytodienoic acid but do not exhibit accumulation during natural senescence or upon stress treatment. Chlorophyll loss during aging and upon dark incubation was not altered, suggesting that these oxylipins are not involved in these processes. In contrast, lipoxygenase 2-RNAi lines and the allene oxid synthase-deficient mutant dde2 were less sensitive to sorbitol than the wild type, indicating that oxylipins contribute to the response to sorbitol stress.Senescence is an important, highly regulated process at the end of development. Senescence is characterized by breakdown of organelles and molecules, export and transport of these nutrients to other organs/parts of the organism, and finally programmed cell death of the senescing organ.The process of senescence has been intensively studied in leaves, and morphological as well as molecular changes in senescing leaves have been described. Yellowing as a consequence of chlorophyll and chloroplast degradation is the most obvious process during natural leaf senescence. In addition, gene expression changes dramatically during senescence. Some senescence-associated genes (SAG, SEN) have been reported that are induced during this process, and several of the encoded proteins function in macromolecule degradation, detoxification and defense metabolism, or signal transduction (Gepstein et al., 2003). Based on the degradation of chloroplasts and macromolecules, leaf metabolism changes from carbon assimilation to catabolism (Lim et al., 2007).The initiation and progression of senescence is regulated by endogenous as well as exogenous factors. Among the endogenous factors, the developmental status of the organ and of the whole plant (e.g. age and progress in flowering and seed production) has a great impact on the process of senescence. Different stress factors such as pathogen attack, drought, osmotic stress, heat, cold, ozone, UV light, and shading can induce or accelerate senescence (Quirino et al., 2000). Phytohormones are very important regulators that integrate information about the developmental status and the environmental factors. Cytokinins are antagonistic signals and delay senescence. Endogenous levels of cytokinins decrease during senescence, and exogenous application and transgenic approaches, enhancing endogenous levels of these compounds, lead to delayed senescence (Gan and Amasino, 1995). In contrast, the gaseous phytohormone ethylene is known to induce and accelerate senescence (John et al., 1995). There are also several indications that abscisic acid modulates senescence (van der Graaff et al., 2006). The roles of other phytohormones/signaling compounds such as auxin, salicylic acid, and jasmonates are less clear (Lim et al., 2007).Jasmonates are oxylipin signaling molecules derived from linolenic acid. The term jasmonates comprises 12-oxo-phytodienoic acid (OPDA), jasmonic acid (JA), and derivatives such as the methyl ester and amino acid conjugates of JA. One of the first biological activities described for these compounds was the promotion of senescence in oat (Avena sativa) leaves by methyl jasmonate (MeJa) isolated from Artemisia absinthium (Ueda and Kato, 1980). Later on, the induction of senescence-like phenotypes by exogenous application of MeJa was also found in other plant species (Ueda and Kato, 1980; Weidhase et al., 1987a; He et al., 2002). On the molecular level, this senescence-promoting effect of MeJa is accompanied by chlorophyll loss and decreases in Rubisco and photosynthesis (Weidhase et al., 1987a, 1987b). In addition, expression of some senescence-up-regulated genes is also responsive to JA; examples are SEN1, SEN4, SEN5, SAG12, SAG14, and SAG15 (Park et al., 1998; Schenk et al., 2000; He et al., 2002). Due to the results described above, jasmonates have been described for decades as compounds with senescence-promoting activities, while the function of these compounds in natural senescence in planta was critically discussed (Parthier, 1990; Sembdner and Parthier, 1993; Creelman and Mullet, 1997; Wasternack, 2007; Balbi and Devoto, 2008; Reinbothe et al., 2009). Additional indications for a role of jasmonates in regulating senescence are the transient up-regulation of expression of some enzymes involved in JA biosynthesis, such as allene oxide synthase (AOS) and OPDA reductase 3 (OPR3), and the increase in JA levels during natural senescence (He et al., 2002; van der Graaff et al., 2006). Furthermore, alterations in natural and induced senescence have been reported for some mutants with defects in the JA pathway. The mutant coi1, which is impaired in JA signaling, exhibited delayed chlorophyll loss upon dark incubation of detached leaves (Castillo and Leon, 2008). Plants with reduced expression of the 3-ketoacyl-CoA-thiolase KAT2, which is involved in β-oxidation and JA production, showed delayed yellowing during natural senescence and upon dark incubation of detached leaves (Castillo and Leon, 2008).However, there are also several reports that cast doubt on an important function of JA in senescence. For most mutants in JA biosynthesis or signaling, no differences in natural senescence are apparent (He et al., 2002; Schommer et al., 2008). In addition, mutants defective in the expression of AOS or OPR3 do not show altered senescence-like phenotypes upon dark treatment (Schommer et al., 2008; Kunz et al., 2009). It has to be taken into consideration that the knockout in these mutants has pleiotrophic effects during whole plant development. For example, the leaves of plants with reduced expression of the lipase DGL or of OPR3 are larger (Hyun et al., 2008). In addition, several knockout mutants defective in JA biosynthesis or signaling do not produce fertile flowers (Feys et al., 1994; McConn and Browse, 1996; Sanders et al., 2000; Stintzi and Browse, 2000; Ishiguro et al., 2001; von Malek et al., 2002). These changes in development might affect other developmental processes such as senescence.To investigate the function of jasmonates in senescence in more detail, we compared the oxylipin profile of wild-type leaves during natural senescence and upon stress induction of senescence-like phenotypes. The analysis of lipoxygenase 2 (LOX2)-RNA interference (RNAi) plants, which produce low basal levels of oxylipins but are impaired in the accumulation of OPDA and JA during senescence or in response to stress, indicates that 13-LOX products are not necessary for natural senescence or dark-induced chlorophyll loss but are involved in the response to sorbitol.     

Autophagy Plays a Role in Chloroplast Degradation during Senescence in Individually Darkened Leaves

Chloroplasts contain approximately 80% of total leaf nitrogen and represent a major source of recycled nitrogen during leaf senescence. While bulk degradation of the cytosol and organelles in plants is mediated by autophagy, its role in chloroplast catabolism is largely unknown. We investigated the effects of autophagy disruption on the number and size of chloroplasts during senescence. When leaves were individually darkened, senescence was promoted similarly in both wild-type Arabidopsis (Arabidopsis thaliana) and in an autophagy-defective mutant, atg4a4b-1. The number and size of chloroplasts decreased in darkened leaves of wild type, while the number remained constant and the size decrease was suppressed in atg4a4b-1. When leaves of transgenic plants expressing stroma-targeted DsRed were individually darkened, a large accumulation of fluorescence in the vacuolar lumen was observed. Chloroplasts exhibiting chlorophyll fluorescence, as well as Rubisco-containing bodies, were also observed in the vacuole. No accumulation of stroma-targeted DsRed, chloroplasts, or Rubisco-containing bodies was observed in the vacuoles of the autophagy-defective mutant. We have succeeded in demonstrating chloroplast autophagy in living cells and provide direct evidence of chloroplast transportation into the vacuole.Chloroplasts contain 75% to 80% of total leaf nitrogen mainly as proteins (Makino and Osmond, 1991). During leaf senescence, chloroplast proteins are gradually degraded as a major source of nitrogen for new growth (Wittenbach, 1978; Friedrich and Huffaker, 1980; Mae et al., 1984), correlating with a decline in photosynthetic activity, while chloroplasts gradually shrink and transform into gerontoplasts, characterized by the disintegration of the thylakoid membranes and accumulation of plastoglobuli (for a recent review, see Krupinska, 2006). Concomitantly, a decline in the cellular population of chloroplasts is also evident in many cases, for example, during natural (Kura-Hotta et al., 1990; Inada et al., 1998), dark-induced (Wittenbach et al., 1982), and nutrient-limited senescence (Mae et al., 1984; Ono et al., 1995), suggesting the existence of a whole chloroplast degradation system. Some electron microscopic studies have shown whole chloroplasts in the central vacuole, which is rich in lytic hydrolases (Wittenbach et al., 1982; Minamikawa et al., 2001). However, there is no direct evidence of chloroplasts moving into the vacuole in living cells and the mechanism of transport is not yet understood (Hörtensteiner and Feller, 2002; Krupinska, 2006).The most abundant chloroplast protein is Rubisco (EC 4.1.1.39), comprising approximately 50% of the soluble protein (Wittenbach, 1978). The amount of Rubisco decreases rapidly in the early phase of leaf senescence, although more slowly in the later phase (Friedrich and Huffaker, 1980; Mae et al., 1984). In contrast, the chloroplast number remains relatively constant, making it impossible to explain Rubisco loss solely by whole chloroplast degradation. However, the mechanism of intrachloroplastic Rubisco degradation is still unknown (for review, see Feller et al., 2008). Using immunoelectron microscopy, we previously demonstrated in naturally senescing wheat (Triticum aestivum) leaves that Rubisco is released from chloroplasts into the cytoplasm and transported to the vacuole for subsequent degradation in small spherical bodies, named Rubisco-containing bodies (RCBs; Chiba et al., 2003). Similar chloroplast-derived structures were also subsequently confirmed in senescent leaves of soybean (Glycine max) and/or Arabidopsis (Arabidopsis thaliana) by electron microscopy (Otegui et al., 2005), and recently in tobacco (Nicotiana tabacum) leaves by immunoelectron microscopy, although the authors gave them a different name, Rubisco vesicular bodies (Prins et al., 2008). RCBs have double membranes, which seem to be derived from the chloroplast envelope; thus, the RCB-mediated degradation of stromal proteins represents a potential mechanism for chloroplast shrinkage during senescence. We recently demonstrated that Rubisco and stroma-targeted fluorescent proteins can be mobilized to the vacuole by ATG-dependent autophagy via RCBs, using leaves treated with concanamycin A, a vacuolar H⁺-ATPase inhibitor (Ishida et al., 2008). To investigate further, we wished to observe chloroplast autophagy and degradation directly in living cells to determine whether autophagy is responsible for chloroplast shrinkage and whether it is involved in the vacuolar degradation of whole chloroplasts during leaf senescence.Autophagy is known to be a major system for the bulk degradation of intracellular proteins and organelles in the vacuole in yeast and plants, or the lysosome in animals (for detailed mechanisms, see reviews by Ohsumi, 2001; Levine and Klionsky, 2004; Thompson and Vierstra, 2005; Bassham et al., 2006). In those systems, a portion of the cytoplasm, including entire organelles, is engulfed in membrane-bound vesicles and delivered to the vacuole/lysosome. A recent genome-wide search confirmed that Arabidopsis has many genes homologous to the yeast autophagy genes (ATGs; Doelling et al., 2002; Hanaoka et al., 2002; for detailed functions of ATGs, see the reviews noted above). Using knockout mutants of ATGs and a monitoring system with an autophagy marker, GFP-ATG8, numerous studies have demonstrated the presence of the autophagy system in plants and its importance in several biological processes (Yoshimoto et al., 2004; Liu et al., 2005; Suzuki et al., 2005; Thompson et al., 2005; Xiong et al., 2005, 2007; Fujiki et al., 2007; Phillips et al., 2008). These articles suggest that autophagy plays an important role in nutrient recycling during senescence, especially in nutrient-starved plants. The atg mutants exhibited an accelerated loss of some chloroplast proteins, but not all, under nutrient-starved conditions and during senescence, suggesting that autophagy is not the sole mechanism for the degradation of chloroplast proteins; other, as yet unidentified systems must be responsible for the degradation of chloroplast contents when the ATG system is compromised (Levine and Klionsky, 2004; Bassham et al., 2006). However, it still remains likely that autophagy is responsible for the vacuolar degradation of chloroplasts in wild-type plants.Prolonged observation is generally required to follow leaf senescence events in naturally aging leaves and senescence-associated processes tend to become chaotic over time. To observe chloroplast degradation over a short period, and to draw clear conclusions, a suitable experimental model of leaf senescence is required. Weaver and Amasino (2001) reported that senescence is rapidly induced in individually darkened leaves (IDLs) of Arabidopsis, but retarded in plants subjected to full darkness. In addition, Keech et al. (2007) observed a significant decrease of both the number and size of chloroplasts in IDLs within 6 d.In this study, using IDLs as a senescence model, we aimed to investigate the involvement of autophagy in chloroplast degradation. We show direct evidence for the transport of whole chloroplasts and RCBs to the vacuole by autophagy.     

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.     

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).     

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.     

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.