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.     

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

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

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

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

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

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

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.     

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

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

Surveying Rubisco Diversity and Temperature Response to Improve Crop Photosynthetic Efficiency

The threat to global food security of stagnating yields and population growth makes increasing crop productivity a critical goal over the coming decades. One key target for improving crop productivity and yields is increasing the efficiency of photosynthesis. Central to photosynthesis is Rubisco, which is a critical but often rate-limiting component. Here, we present full Rubisco catalytic properties measured at three temperatures for 75 plants species representing both crops and undomesticated plants from diverse climates. Some newly characterized Rubiscos were naturally “better” compared to crop enzymes and have the potential to improve crop photosynthetic efficiency. The temperature response of the various catalytic parameters was largely consistent across the diverse range of species, though absolute values showed significant variation in Rubisco catalysis, even between closely related species. An analysis of residue differences among the species characterized identified a number of candidate amino acid substitutions that will aid in advancing engineering of improved Rubisco in crop systems. This study provides new insights on the range of Rubisco catalysis and temperature response present in nature, and provides new information to include in models from leaf to canopy and ecosystem scale.In a changing climate and under pressure from a population set to hit nine billion by 2050, global food security will require massive changes to the way food is produced, distributed, and consumed (Ort et al., 2015). To match rising demand, agricultural production must increase by 50 to 70% in the next 35 years, and yet the gains in crop yields initiated by the green revolution are slowing, and in some cases, stagnating (Long and Ort, 2010; Ray et al., 2012). Among a number of areas being pursued to increase crop productivity and food production, improving photosynthetic efficiency is a clear target, offering great promise (Parry et al., 2007; von Caemmerer et al., 2012; Price et al., 2013; Ort et al., 2015). As the gatekeeper of carbon entry into the biosphere and often acting as the rate-limiting step of photosynthesis, Rubisco, the most abundant enzyme on the planet (Ellis, 1979), is an obvious and important target for improving crop photosynthetic efficiency.Rubisco is considered to exhibit comparatively poor catalysis, in terms of catalytic rate, specificity, and CO₂ affinity (Tcherkez et al., 2006; Andersson, 2008), leading to the suggestion that even small increases in catalytic efficiency may result in substantial improvements to carbon assimilation across a growing season (Zhu et al., 2004; Parry et al., 2013; Galmés et al., 2014a; Carmo-Silva et al., 2015). If combined with complimentary changes such as optimizing other components of the Calvin Benson or photorespiratory cycles (Raines, 2011; Peterhansel et al., 2013; Simkin et al., 2015), optimized canopy architecture (Drewry et al., 2014), or introducing elements of a carbon concentrating mechanism (Furbank et al., 2009; Lin et al., 2014a; Hanson et al., 2016; Long et al., 2016), Rubisco improvement presents an opportunity to dramatically increase the photosynthetic efficiency of crop plants (McGrath and Long, 2014; Long et al., 2015; Betti et al., 2016). A combination of the available strategies is essential for devising tailored solutions to meet the varied requirements of different crops and the diverse conditions under which they are typically grown around the world.Efforts to engineer an improved Rubisco have not yet produced a “super Rubisco” (Parry et al., 2007; Ort et al., 2015). However, advances in engineering precise changes in model systems continue to provide important developments that are increasing our understanding of Rubisco catalysis (Spreitzer et al., 2005; Whitney et al., 2011a, 2011b; Morita et al., 2014; Wilson et al., 2016), regulation (Andralojc et al., 2012; Carmo-Silva and Salvucci, 2013; Bracher et al., 2015), and biogenesis (Saschenbrecker et al., 2007; Whitney and Sharwood, 2008; Lin et al., 2014b; Hauser et al., 2015; Whitney et al., 2015).A complementary approach is to understand and exploit Rubisco natural diversity. Previous characterization of Rubisco from a limited number of species has not only demonstrated significant differences in the underlying catalytic parameters, but also suggests that further undiscovered diversity exists in nature and that the properties of some of these enzymes could be beneficial if present in crop plants (Carmo-Silva et al., 2015). Recent studies clearly illustrate the variation possible among even closely related species (Galmés et al., 2005, 2014b, 2014c; Kubien et al., 2008; Andralojc et al., 2014; Prins et al., 2016).Until recently, there have been relatively few attempts to characterize the consistency, or lack thereof, of temperature effects on in vitro Rubisco catalysis (Sharwood and Whitney, 2014), and often studies only consider a subset of Rubisco catalytic properties. This type of characterization is particularly important for future engineering efforts, enabling specific temperature effects to be factored into any attempts to modify crops for a future climate. In addition, the ability to coanalyze catalytic properties and DNA or amino acid sequence provides the opportunity to correlate sequence and biochemistry to inform engineering studies (Christin et al., 2008; Kapralov et al., 2011; Rosnow et al., 2015). While the amount of gene sequence information available grows rapidly with improving technology, knowledge of the corresponding biochemical variation resulting has yet to be determined (Cousins et al., 2010; Carmo-Silva et al., 2015; Sharwood and Whitney, 2014; Nunes-Nesi et al., 2016).This study aimed to characterize the catalytic properties of Rubisco from diverse species, comprising a broad range of monocots and dicots from diverse environments. The temperature dependence of Rubisco catalysis was evaluated to tailor Rubisco engineering for crop improvement in specific environments. Catalytic diversity was analyzed alongside the sequence of the Rubisco large subunit gene, rbcL, to identify potential catalytic switches for improving photosynthesis and productivity. In vitro results were compared to the average temperature of the warmest quarter in the regions where each species grows to investigate the role of temperature in modulating Rubisco catalysis.     

Focus on Ethylene: Ethylene and the Regulation of Physiological and Morphological Responses to Nutrient Deficiencies

To cope with nutrient deficiencies, plants develop both morphological and physiological responses. The regulation of these responses is not totally understood, but some hormones and signaling substances have been implicated. It was suggested several years ago that ethylene participates in the regulation of responses to iron and phosphorous deficiency. More recently, its role has been extended to other deficiencies, such as potassium, sulfur, and others. The role of ethylene in so many deficiencies suggests that, to confer specificity to the different responses, it should act through different transduction pathways and/or in conjunction with other signals. In this update, the data supporting a role for ethylene in the regulation of responses to different nutrient deficiencies will be reviewed. In addition, the results suggesting the action of ethylene through different transduction pathways and its interaction with other hormones and signaling substances will be discussed.When plants suffer from a mineral nutrient deficiency, they develop morphological and physiological responses (mainly in their roots) aimed to facilitate the uptake and mobilization of the limiting nutrient. After the nutrient has been acquired in enough quantity, these responses need to be switched off to avoid toxicity and conserve energy. In recent years, different plant hormones (e.g. ethylene, auxin, cytokinins, jasmonic acid, abscisic acid, brassinosteroids, GAs, and strigolactones) have been implicated in the regulation of these responses (Romera et al., 2007, 2011, 2015; Liu et al., 2009; Rubio et al., 2009; Kapulnik et al., 2011; Kiba et al., 2011; Iqbal et al., 2013; Zhang et al., 2014).Before the 1990s, there were several publications relating ethylene and nutrient deficiencies (cited in Lynch and Brown [1997] and Romera et al. [1999]) without establishing a direct implication of ethylene in the regulation of nutrient deficiency responses. In 1994, Romera and Alcántara (1994) published an article in Plant Physiology suggesting a role for ethylene in the regulation of Fe deficiency responses. In 1999, Borch et al. (1999) showed the participation of ethylene in the regulation of P deficiency responses. Since then, evidence has been accumulating in support of a role for ethylene in the regulation of both Fe (Romera et al., 1999, 2015; Waters and Blevins, 2000; Lucena et al., 2006; Waters et al., 2007; García et al., 2010, 2011, 2013, 2014; Yang et al., 2014) and P deficiency responses (Kim et al., 2008; Lei et al., 2011; Li et al., 2011; Nagarajan and Smith, 2012; Wang et al., 2012, 2014c). Both Fe and P may be poorly available in most soils, and plants develop similar responses under their deficiencies (Romera and Alcántara, 2004; Zhang et al., 2014). More recently, a role for ethylene has been extended to other deficiencies, such as K (Shin and Schachtman, 2004; Jung et al., 2009; Kim et al., 2012), S (Maruyama-Nakashita et al., 2006; Wawrzyńska et al., 2010; Moniuszko et al., 2013), and B (Martín-Rejano et al., 2011). Ethylene has also been implicated in both N deficiency and excess (Tian et al., 2009; Mohd-Radzman et al., 2013; Zheng et al., 2013), and its participation in Mg deficiency has been suggested (Hermans et al., 2010).In this update, we will review the information supporting a role for ethylene in the regulation of different nutrient deficiency responses. For information relating ethylene to other aspects of plant mineral nutrition, such as N₂ fixation and responses to excess of nitrate or essential heavy metals, the reader is referred to other reviews (for review, see Maksymiec, 2007; Mohd-Radzman et al., 2013; Steffens, 2014).     

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.