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

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

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

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

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

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.     

Genetic Interactions between PEROXIN12 and Other Peroxisome-Associated Ubiquitination Components

Most eukaryotic cells require peroxisomes, organelles housing fatty acid β-oxidation and other critical metabolic reactions. Peroxisomal matrix proteins carry peroxisome-targeting signals that are recognized by one of two receptors, PEX5 or PEX7, in the cytosol. After delivering the matrix proteins to the organelle, these receptors are removed from the peroxisomal membrane or matrix. Receptor retrotranslocation not only facilitates further rounds of matrix protein import but also prevents deleterious PEX5 retention in the membrane. Three peroxisome-associated ubiquitin-protein ligases in the Really Interesting New Gene (RING) family, PEX2, PEX10, and PEX12, facilitate PEX5 retrotranslocation. However, the detailed mechanism of receptor retrotranslocation remains unclear in plants. We identified an Arabidopsis (Arabidopsis thaliana) pex12 Glu-to-Lys missense allele that conferred severe peroxisomal defects, including impaired β-oxidation, inefficient matrix protein import, and decreased growth. We compared this pex12-1 mutant to other peroxisome-associated ubiquitination-related mutants and found that RING peroxin mutants displayed elevated PEX5 and PEX7 levels, supporting the involvement of RING peroxins in receptor ubiquitination in Arabidopsis. Also, we observed that disruption of any Arabidopsis RING peroxin led to decreased PEX10 levels, as seen in yeast and mammals. Peroxisomal defects were exacerbated in RING peroxin double mutants, suggesting distinct roles of individual RING peroxins. Finally, reducing function of the peroxisome-associated ubiquitin-conjugating enzyme PEX4 restored PEX10 levels and partially ameliorated the other molecular and physiological defects of the pex12-1 mutant. Future biochemical analyses will be needed to determine whether destabilization of the RING peroxin complex observed in pex12-1 stems from PEX4-dependent ubiquitination on the pex12-1 ectopic Lys residue.Oilseed plants obtain energy for germination and early development by utilizing stored fatty acids (Graham, 2008). This β-oxidation of fatty acids to acetyl-CoA occurs in peroxisomes, organelles that also house other important metabolic reactions, including the glyoxylate cycle, several steps in photorespiration, and phytohormone production (Hu et al., 2012). For example, indole-3-butyric acid (IBA) is β-oxidized into the active auxin indole-3-acetic acid (IAA) in peroxisomes (Zolman et al., 2000, 2007, 2008; Strader et al., 2010; Strader and Bartel, 2011). Many peroxisomal metabolic pathways generate reactive oxygen species (Inestrosa et al., 1979; Hu et al., 2012), and peroxisomes also house antioxidative enzymes, like catalase and ascorbate peroxidase, to detoxify hydrogen peroxide (Wang et al., 1999; Mhamdi et al., 2012).Peroxisomes can divide by fission or be synthesized de novo from the endoplasmic reticulum (ER). Preperoxisomes with peroxisomal membrane proteins bud from the ER and fuse, allowing matrix proteins to be imported to form mature peroxisomes (van der Zand et al., 2012; Mayerhofer, 2016). Peroxin (PEX) proteins facilitate peroxisome biogenesis and matrix protein import. Most peroxins are involved in importing proteins destined for the peroxisome matrix, which are imported after recognition of a type 1 or type 2 peroxisome-targeting signal (PTS). The PTS1 is a tripeptide located at the C terminus of most peroxisome-bound proteins (Gould et al., 1989; Chowdhary et al., 2012). The less common PTS2 is a nonapeptide usually located near the N terminus (Swinkels et al., 1991; Reumann, 2004). PTS1 proteins are recognized by PEX5 (van der Leij et al., 1993; Zolman et al., 2000), PTS2 proteins are recognized by PEX7 (Marzioch et al., 1994; Braverman et al., 1997; Woodward and Bartel, 2005), and PEX7 binds to PEX5 to allow matrix protein delivery in plants and mammals (Otera et al., 1998; Hayashi et al., 2005; Woodward and Bartel, 2005). The cargo-receptor complex docks with the membrane peroxins PEX13 and PEX14 (Urquhart et al., 2000; Otera et al., 2002; Woodward et al., 2014), and PEX5 assists cargo translocation into the peroxisomal matrix (Meinecke et al., 2010) before dissociating from its cargo (Freitas et al., 2011).After cargo delivery, PEX5 is recycled to enable further rounds of cargo recruitment (Thoms and Erdmann, 2006). This process requires a set of peroxins that is implicated in ubiquitinating PEX5 so that it can be retrotranslocated back to the cytosol. PEX5 ubiquitination is best understood in yeast. In Saccharomyces cerevisiae, Pex5 is monoubiquitinated through the action of the peroxisome-tethered ubiquitin-conjugating enzyme Pex4 and the peroxisomal ubiquitin-protein ligase Pex12 (Platta et al., 2009) and returned to the cytosol with the assistance of a peroxisome-tethered ATPase complex containing Pex1 and Pex6 (Grimm et al., 2012). S. cerevisiae Pex5 also can be polyubiquitinated and targeted for proteasomal degradation (Kiel et al., 2005). The cytosolic ubiquitin-conjugating enzyme Ubc4 cooperates with the peroxisomal ubiquitin-protein ligase Pex2 to polyubiquitinate Pex5 (Platta et al., 2009). Pex10 has ubiquitin-protein ligase activity (Williams et al., 2008; Platta et al., 2009; El Magraoui et al., 2012), but whether Pex10 directly ubiquitinates Pex5 is controversial. Pex10 promotes Ubc4-dependent Pex5 polyubiquitination when Pex4 is absent (Williams et al., 2008); however, Pex10 is not essential for Pex5 mono- or polyubiquitination (Platta et al., 2009), but rather enhances both Pex4/Pex12- and Ubc4/Pex2-mediated ubiquitination (El Magraoui et al., 2012). Recycling of the PTS2 receptor PEX7 is less understood, although the Pex5 recycling pathways are implicated in shuttling and degrading Pex7 in Pichia pastoris (Hagstrom et al., 2014).Although PEX5 ubiquitination has not been directly demonstrated in plants, the implicated peroxins are conserved in Arabidopsis, and several have been connected to PEX5 retrotranslocation. The PEX4 ubiquitin-conjugating enzyme binds to PEX22, which is predicted to be a peroxisomal membrane protein based on ability to restore peroxisome function to yeast mutants (Zolman et al., 2005). The pex4-1 mutant displays increased membrane-associated PEX5 (Ratzel et al., 2011; Kao and Bartel, 2015), suggesting that ubiquitin supplied by PEX4 promotes PEX5 retrotranslocation. PEX1 and PEX6 are members of the ATPases associated with diverse cellular activities (AAA) family and are tethered to peroxisomes by the peroxisomal membrane protein PEX26 (Goto et al., 2011; Li et al., 2014). The pex6-1 mutant displays PTS1 import defects and decreased PEX5 levels (Zolman and Bartel, 2004), suggesting that impaired PEX5 recycling can lead to increased PEX5 degradation. Indeed, pex4-1 restores PEX5 levels in the pex6-1 mutant (Ratzel et al., 2011), suggesting that Arabidopsis PEX4 also is involved in PEX5 ubiquitination and degradation when retrotranslocation is impeded.In addition to allowing for further rounds of PTS1 cargo import, several lines of evidence suggest that in the absence of efficient retrotranslocation, PEX5 retention in the peroxisomal membrane impairs peroxisome function. Slightly reducing levels of the PEX13 docking peroxin ameliorates the physiological defects of pex4-1 without restoring matrix protein import (Ratzel et al., 2011), presumably because decreasing PEX5 docking reduces its accumulation in the peroxisomal membrane. In addition, overexpressing PEX5 exacerbates rather than ameliorates the peroxisomal defects of pex4-1 (Kao and Bartel, 2015), suggesting that pex4-1 defects are linked to excessive PEX5 lingering in the peroxisome membrane rather than a lack of PEX5 available for import.The three Really Interesting New Gene (RING) peroxins (PEX2, PEX10, and PEX12) from Arabidopsis each possesses in vitro ubiquitin-protein ligase activity (Kaur et al., 2013). Null mutations in the RING peroxin genes confer embryo lethality in Arabidopsis (Hu et al., 2002; Schumann et al., 2003; Sparkes et al., 2003; Fan et al., 2005; Prestele et al., 2010), necessitating other approaches to study the in vivo functions of these peroxins. Expressing RING peroxins with mutations in the C-terminal zinc-binding RING domains (ΔZn) confers matrix protein import defects for PEX2-ΔZn and photorespiration defects for PEX10-ΔZn but no apparent defects for PEX12-ΔZn (Prestele et al., 2010). Targeting individual RING peroxins using RNAi confers β-oxidation deficiencies and impairs PTS1 cargo import (Fan et al., 2005; Nito et al., 2007). A screen for delayed matrix protein degradation (Burkhart et al., 2013) uncovered a missense pex2-1 mutant and a splicing pex10-2 mutant that both display PTS1 import defects (Burkhart et al., 2014), suggesting roles in regulating the PTS1 receptor, PEX5. A missense pex12 mutant (aberrant peroxisome morphology 4, apm4) has defects in β-oxidation and PTS1 import and increased membrane-associated PEX5 (Mano et al., 2006). These findings highlight the essential roles of the RING peroxins in Arabidopsis development and peroxisomal functions, but the RING peroxin interactions and the individual roles of the RING peroxins in PEX5 retrotranslocation remain incompletely understood.In this study, we describe a missense pex12-1 mutant recovered from a forward genetic screen for β-oxidation deficient mutants. The pex12-1 mutant displayed severe peroxisomal defects, including reduced growth, β-oxidation deficiencies, matrix protein import defects, and inefficient processing of PTS2 proteins. Comparing single and double mutants with impaired RING peroxins revealed that each RING peroxin contributes to complex stability and influences PEX5 accumulation. Furthermore, decreasing PEX4 function ameliorated pex12-1 defects, suggesting that the Glu-to-Lys substitution in pex12-1 lures ubiquitination, perhaps by pex12-1 itself, leading to PEX4-dependent degradation of the mutant protein.     

Nontransgenic Genome Modification in Plant Cells

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