Dynamics and Organization of Cortical Microtubules as Revealed by Superresolution Structured Illumination Microscopy

Plants employ acentrosomal mechanisms to organize cortical microtubule arrays essential for cell growth and differentiation. Using structured illumination microscopy (SIM) adopted for the optimal documentation of Arabidopsis (Arabidopsis thaliana) hypocotyl epidermal cells, dynamic cortical microtubules labeled with green fluorescent protein fused to the microtubule-binding domain of the mammalian microtubule-associated protein MAP4 and with green fluorescent protein-fused to the alpha tubulin6 were comparatively recorded in wild-type Arabidopsis plants and in the mitogen-activated protein kinase mutant mpk4 possessing the former microtubule marker. The mpk4 mutant exhibits extensive microtubule bundling, due to increased abundance and reduced phosphorylation of the microtubule-associated protein MAP65-1, thus providing a very useful genetic tool to record intrabundle microtubule dynamics at the subdiffraction level. SIM imaging revealed nano-sized defects in microtubule bundling, spatially resolved microtubule branching and release, and finally allowed the quantification of individual microtubules within cortical bundles. Time-lapse SIM imaging allowed the visualization of subdiffraction, short-lived excursions of the microtubule plus end, and dynamic instability behavior of both ends during free, intrabundle, or microtubule-templated microtubule growth and shrinkage. Finally, short, rigid, and nondynamic microtubule bundles in the mpk4 mutant were observed to glide along the parent microtubule in a tip-wise manner. In conclusion, this study demonstrates the potential of SIM for superresolution time-lapse imaging of plant cells, showing unprecedented details accompanying microtubule dynamic organization.Plant cell growth and differentiation depend on dynamic cortical microtubule organization mechanisms (Ehrhardt, 2008). Such mechanisms include branched microtubule formation and release (Murata et al., 2005; Nakamura et al., 2010; Fishel and Dixit, 2013), microtubule-templated microtubule growth (Chan et al., 2009), angle-of-contact microtubule bundling or catastrophe induction (Dixit and Cyr, 2004; Tulin et al., 2012), severing at microtubule crossovers (Wightman and Turner, 2007), and unique dynamic behavior between steady-state treadmilling and dynamic instability (Shaw et al., 2003).Cortical microtubule dynamics have been studied in vivo and in vitro with total internal reflection microscopy (TIRFM; Vizcay-Barrena et al., 2011), variable-angle emission microscopy (VAEM; Wan et al., 2011), spinning-disc microscopy (SD; Shaw and Lucas, 2011), and confocal laser scanning microscopy (CLSM; Shaw et al., 2003). TIRFM and VAEM provide sufficient resolution and speed but at limited depth of imaging (approximately 200 nm; Martin-Fernandez et al., 2013) and inevitably a very narrow field of view when used for in vivo studies (Mattheyses et al., 2010). Dynamic CLSM imaging suffers from field-of-view limitations while also introducing phototoxicity to the imaged sample. Furthermore, CLSM is based on a speed-to-resolution tradeoff that will necessitate computational extrapolation to bring resolution to affordable levels (Rosero et al., 2014). Finally, SD can provide sufficient depth and speed but otherwise poor resolution, owing to aberrations arising from the sample and the properties of the optics commonly used (Shaw and Ehrhardt, 2013).Microtubule research evolved concomitant with optical microscopy and the development of fluorescent proteins markers, allowing the resolution of microtubule dynamics and organization at video rates (Marc et al., 1998; Shaw and Ehrhardt, 2013). However, the bulk of plant cells organized in tissues and the optical properties of cell walls hamper microscopic observations, so that the delineation of fine details of microtubule organization still relies on laborious transmission electron microscopy (Kang, 2010).Alternatively, in vitro assays using total internal reflection (TIRFM) or Allen’s video-enhanced contrast-differential interference contrast microscopy (Allen et al., 1981) with purified components have advanced the understanding of microtubule-microtubule-associated protein (MAP) interactions while providing mechanistic insight on the function of MAP65 proteins (Tulin et al., 2012; Portran et al., 2013; Stoppin-Mellet et al., 2013), kinesin motors (Song et al., 1997), katanin-mediated microtubule severing (Stoppin-Mellet et al., 2007), and microtubule dynamics (Moore et al., 1997). However, it is explicitly acknowledged that such in vitro assays should be addressed in biologically coherent systems with physiological relevance to microtubule dynamics (Gardner et al., 2013; Zanic et al., 2013). Thus, an ideal approach would be to address microtubule dynamics in the complex cellular environment at spatiotemporal resolutions achieved by in vitro assays.Subdiffraction optical microscopy techniques allow subcellular observations below Abbe’s resolution threshold (Verdaasdonk et al., 2014), complementing the use of transmission electron microscopy. Such approaches permit dynamic subcellular tracking of appropriately tagged structures within living cells (Tiwari and Nagai, 2013). Practically, two superresolution strategies exist. The first involves patterned light illumination, allowing superresolution acquisitions by two fundamentally different methods, stimulated emission depletion (STED; Hell, 2007) and structured illumination microscopy (SIM; Gustafsson, 2000). The second interrogates the precision of fluorophore localization and includes stochastic optical reconstruction microscopy (STORM; Kamiyama and Huang, 2012) and photoactivation localization microscopy (PALM; Sengupta et al., 2012). The above regimes differ in translational and axial resolution, and their temporal efficiency depends on the size of the imaged area. SIM is probably the best compromise for superresolution live imaging, as it offers reasonable lateral resolution (approximately 100 nm; Gustafsson, 2000), which may be reduced to 50 nm (Rego et al., 2012), and sufficient depth of imaging combined with a reasonable axial resolution (approximately 200 nm). SIM allows dynamic imaging in a broader field of view than STED, at biologically meaningful rates compared with PALM and STORM (Kner et al., 2009), and with deeper imaging capacity compared with other superresolution regimes and with TIRFM/VAEM (Leung and Chou, 2011). Superresolution approaches have received limited attention in the plant cell biology field (Fitzgibbon et al., 2010; Kleine-Vehn et al., 2011), and their resolution potential during live imaging was not quantified previously.Here, high-numerical aperture (NA) objectives were combined with SIM for the acquisition and systematic quantification of subdiffraction details of cortical microtubules labeled either with GFP fused to the microtubule-binding domain of mammalian MAP4 (GFP-MBD; Marc et al., 1998) or with GFP fused to alpha tubulin6 (GFP-TUA6; Shaw et al., 2003). For such studies, wild-type plants and a mitogen-activated protein kinase4 (mpk4) mutant, exhibiting extensive microtubule bundling due to the overexpression and underphosphorylation of MAP65-1 (Beck et al., 2010), were used.     

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

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

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

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

Cortical Microtubule Arrays Are Initiated from a Nonrandom Prepattern Driven by Atypical Microtubule Initiation

The ordered arrangement of cortical microtubules in growing plant cells is essential for anisotropic cell expansion and, hence, for plant morphogenesis. These arrays are dismantled when the microtubule cytoskeleton is rearranged during mitosis and reassembled following completion of cytokinesis. The reassembly of the cortical array has often been considered as initiating from a state of randomness, from which order arises at least partly through self-organizing mechanisms. However, some studies have shown evidence for ordering at early stages of array assembly. To investigate how cortical arrays are initiated in higher plant cells, we performed live-cell imaging studies of cortical array assembly in tobacco (Nicotiana tabacum) Bright Yellow-2 cells after cytokinesis and drug-induced disassembly. We found that cortical arrays in both cases did not initiate randomly but with a significant overrepresentation of microtubules at diagonal angles with respect to the cell axis, which coincides with the predominant orientation of the microtubules before their disappearance from the cell cortex in preprophase. In Arabidopsis (Arabidopsis thaliana) root cells, recovery from drug-induced disassembly was also nonrandom and correlated with the organization of the previous array, although no diagonal bias was observed in these cells. Surprisingly, during initiation, only about one-half of the new microtubules were nucleated from locations marked by green fluorescent protein-γ-tubulin complex protein2-tagged γ-nucleation complexes (γ-tubulin ring complex), therefore indicating that a large proportion of early polymers was initiated by a noncanonical mechanism not involving γ-tubulin ring complex. Simulation studies indicate that the high rate of noncanonical initiation of new microtubules has the potential to accelerate the rate of array repopulation.Higher plant cells feature ordered arrays of microtubules at the cell cortex (Ledbetter and Porter, 1963) that are essential for cell and tissue morphogenesis, as revealed by disruption of cortical arrays by drugs that cause microtubule depolymerization (Green, 1962) or stabilization (Weerdenburg and Seagull, 1988) and by loss-of-function mutations in a wide variety of microtubule-associated proteins (Baskin, 2001; Whittington et al., 2001; Buschmann and Lloyd, 2008; Lucas et al., 2011). The structure of these arrays is thought to control the pattern of cell growth primarily by its role in the deposition of cellulose microfibrils, the load-bearing component of the cell wall (Somerville, 2006). Functional relations between cortical microtubules and cellulose microfibrils have been proposed since the early sixties, even before cortical microtubules had been visualized (Green, 1962). Recent live-cell imaging studies have confirmed that cortical microtubules indeed guide the movement of cellulose synthase complexes that produce cellulose microfibrils (Paredez et al., 2006) and have shown further that microtubules position the insertion of most cellulose synthase complexes into the plasma membrane (Gutierrez et al., 2009). These activities of ordered cortical microtubules are proposed to facilitate the organization of cell wall structure, creating material properties that underlie cell growth anisotropy.While organization of the interphase cortical array appears to be essential for cell morphogenesis, this organization is disrupted during the cell cycle as microtubules are rearranged to create the preprophase band, spindle, and phragmoplast during mitosis and cytokinesis (for review, see Wasteneys, 2002). Upon completion of cytokinesis, an organized interphase cortical array is regenerated, but the pathway for this reassembly is not well understood.The plant interphase microtubule array is organized and maintained without centrosomes as organizing centers (for review, see Wasteneys, 2002; Bartolini and Gundersen, 2006; Ehrhardt and Shaw, 2006), and microtubule self-organization is proposed to play an important role in cortical microtubule array ordering (Dixit and Cyr, 2004). In electron micrographs, microtubules have been observed to be closely associated with the plasma membrane (Hardham and Gunning, 1978), and live-cell imaging provides evidence for attachment of microtubules to the cell cortex (Shaw et al., 2003; Vos et al., 2004). The close association to the plasma membrane restricts the cortical microtubules to a quasi two-dimensional plane where they interact through polymerization-driven collisions (Shaw et al., 2003; Dixit and Cyr, 2004). Microtubule encounters at shallow angles (<40°) have a high probability of leading to bundling, while microtubule encounters at steeper angles most likely result in induced catastrophes or microtubule crossovers (Dixit and Cyr, 2004). Several computational modeling studies have since shown that these types of interactions between surface-bound dynamical microtubules can indeed explain spontaneous coalignment of microtubules (Allard et al., 2010; Eren et al., 2010; Hawkins et al., 2010; Tindemans et al., 2010).The question of how the orientation of the cortical array is established with respect to the cell axis is less well understood. One possibility is that microtubules are selectively destabilized with respect to cellular coordinates (Ehrhardt and Shaw, 2006). Indeed, recent results from biological observations and modeling suggest that catastrophic collisions induced at the edges between cell faces or heighted catastrophe rates in cell caps could be sufficient to selectively favor microtubules in certain orientation and hence determine the final orientation of the array (Allard et al., 2010; Eren et al., 2010; Ambrose et al., 2011; Dhonukshe et al., 2012).To date, all models of cortical array assembly assume random initial conditions. However, experimental work by Wasteneys and Williamson (1989a, 1989b) in Nitella tasmanica showed that, during array reassembly after drug-induced disruption, microtubules were initially transverse. This was followed by a less ordered phase and later by the acquisition of the final transverse order. A nonrandom initial ordering was also observed in tobacco (Nicotiana tabacum) Bright Yellow-2 (BY-2) cells by Kumagai et al. (2001), who concluded that the process of transverse array establishment starts with longitudinal order but did not provide quantitative data for the process of array assembly. The initial conditions for the cortical microtubule array formation are important to consider, as they may strongly influence the speed at which order is established and could even prevent it from being established over a biologically relevant time scale.In this study, we used live-cell imaging to follow and record the whole transition from the cortical microtubule-free state to the final transverse array and used digital tracking algorithms to quantify the microtubule order. Nucleation stands out as a central process to characterize during array initiation. Lacking a central body to organize microtubule nucleations, the higher plant cell has dispersed nucleation complexes (Wasteneys and Williamson, 1989a, 1989b; Chan et al., 2003; Shaw et al., 2003; Murata et al., 2005; Pastuglia et al., 2006; Nakamura et al., 2010). Therefore, we performed high time resolution observations to quantify nucleation complex recruitment, nucleation rates, and microtubule nucleation angles. We found evidence for a highly nonrandom initial ordering state that features diagonal microtubule orientation and an atypical microtubule initiation mechanism. Simulation analysis indicates that these atypical nucleations have the potential to accelerate the recovery of cortical array density.     

A +TIP for a smooth trip

Is there a cellular mechanism for preventing a depolymerizing microtubule track from “slipping out from under” its cargo? A recent study in budding yeast indicates that when a chromosome is transported to the minus end of a spindle microtubule, its kinetochore-bound microtubule plus end–tracking protein (+TIP) Stu2 may move to the plus end to promote rescue; i.e., to switch the depolymerizing end to a polymerizing end. The possibility that other +TIPs may play a similar role in sustaining a microtubule track during vesicular transport deserves investigation.Microtubule motor proteins such as dynein and kinesins are responsible for transporting cellular cargos along microtubule tracks (Vale, 2003). The net direction and speed of cargo movement, however, are likely to be regulated in a very complicated fashion, especially when a cargo is bound to multiple motors with opposite directionalities (Vale 2003; Mallik and Gross 2004; Levi et al., 2006). The fact that the microtubule track is not very stable further complicates matters. The plus ends of microtubules, which face the cell periphery in most cell types, are highly dynamic, exhibiting alternating periods of polymerization (growth) and depolymerization (shrinkage; Desai and Mitchison, 1997). Such plus end dynamics may be useful for searching and capturing relatively stationary cargos near the cell periphery that need to be transported inward (Vaughan et al., 2002). However, the dynamic nature of the track can also create an obvious problem for the transport process. If a microtubule''s rate of shrinkage is greater than the rate of cargo transport, then the microtubule may shrink past an attached cargo, causing its dissociation from the track. Does this happen in cells, or do cells have a mechanism to prevent it?Although this question has never been directly addressed, a recent study on budding yeast chromosome segregation has shed new light on the issue (Tanaka et al., 2005). In this study, the authors took advantage of a strategy that allowed them to specifically shut off the function of a single kinetochore, thereby preventing it from attaching to a spindle microtubule while, at the same time, permitting other kinetochores to attach to the spindle. After the function of this single kinetochore was switched back on, the behavior of its associated chromosome on a spindle microtubule was subjected to a detailed image analysis. Several important insights from this study on chromosome–microtubule interactions during mitosis have been recently reviewed (Bloom 2005), and, thus, only those observations that pertain to cargo transport will be highlighted here. The chromosome was first seen to undergo a lateral interaction with the microtubule followed by minus end–directed transport toward the pole. The mechanism of the minus end–directed transport is not entirely clear, although a member of the kinesin-14 family, Kar3, may be one of the players in this process (Tanaka et al., 2005). During transport, the attached microtubule can undergo shrinkage with a rate higher than that of the minus end–directed chromosome movement (Tanaka et al., 2005); however, it never shrank beyond the position of the cargo. Such exquisite control over the extent of shrinkage appears to rely on a conversation between the cargo and the plus end of the microtubule that is mediated by the microtubule plus end–tracking protein Stu2 (Bloom 2005; Tanaka et al., 2005).Microtubule plus end–tracking proteins (+TIPs) are a class of proteins that use different structural motifs or specific targeting mechanisms to localize to the dynamic plus ends of microtubules (Carvalho et al., 2003; Akhmanova and Hoogenraad, 2005). Although most +TIPs associate with only the growing ends of microtubules, several +TIPs also localize to the shrinking ends (Carvalho et al., 2003; 2004; Akhmanova and Hoogenraad, 2005; Mennella et al., 2005; Sproul et al., 2005; Molk et al., 2006; Wu et al., 2006). Many +TIPs have been found to impact microtubules by either promoting their growth or promoting dynamic behavior. Stu2 is a member of the XMAP215/TOG/Dis1/DdCP224 family of proteins that have been shown to affect microtubule dynamics in multiple ways depending on different experimental conditions (Ohkura et al., 2001; Popov and Karsenti 2003; Holmfeldt et al., 2004; Akhmanova and Hoogenraad, 2005). In vitro, Stu2 binds to the plus ends of preformed microtubules and promotes catastrophe, which is a switch from growth to shrinkage (van Breugel et al., 2003). In vivo studies using mutants of Stu2, however, indicate that Stu2 promotes microtubule growth (Severin et al., 2001) and the dynamics of both kinetochore and cytoplasmic microtubules (Kosco et al., 2001; Pearson et al., 2003). During anaphase B spindle elongation, Stu2 may antagonize the function of Kip3 (a kinesin-13 family member) to promote the plus end polymerization of overlapping microtubules (Severin et al., 2001). Although it is not fully understood how or why Stu2 is so versatile, it is well recognized that the in vivo interactions among +TIPs are very complicated, and the loss of function of a +TIP in vivo may decrease or increase the accumulation of other +TIPs that also regulate microtubule dynamics (Carvalho et al., 2003; 2004; Lansbergen et al., 2004; Akhmanova and Hoogenraad, 2005; Galjart 2005; Komarova et al., 2005). Tanaka et al. (2005) identified Stu2 as a rescue (a switch from shrinkage to growth) factor based not on phenotypic studies of Stu2 mutants but, instead, on a direct observation of the relationship between microtubule plus end behavior and Stu2 localization. They found that Stu2 was localized at the plus ends of microtubules emanating from the spindle pole body, and, during periods of microtubule shrinkage, Stu2 levels at the plus ends were decreased. Interestingly, Stu2 was also localized at the unbound kinetochore. When the kinetochore subsequently attached laterally to a spindle microtubule and underwent minus end–directed transport, the Stu2 proteins were transported from the kinetochore to the microtubule plus end. The arrival of Stu2 at the plus end closely correlated to the rescue of the shrinking microtubule (Tanaka et al., 2005). These observations strongly suggest that the Stu2 carried by the kinetochore may serve as a rescue factor for the microtubule track, preventing it from vanishing before the migrating chromosome.Could such a scenario exist during microtubule-dependent transport of nonchromosomal cargoes during interphase? We do not yet know the answer. However, based on published studies, it seems reasonable to hypothesize that other +TIPs, especially the cytoplasmic linker protein CLIP-170, may function in a manner similar to yeast Stu2 to ensure a safe trip for a minus end–directed cargo. CLIP-170 contains CAP-Gly microtubule-binding motifs at its NH₂ terminus and was initially identified as a protein required for linking endocytic vesicles to microtubules in vitro (Pierre et al., 1992; Rickard and Kreis 1996). Later, CLIP-170 was identified as a founding member of the microtubule plus end–tracking proteins (Perez et al., 1999). The connection between CLIP-170''s in vitro endosome–microtubule linking property and its in vivo plus end tracking behavior has not been clearly made. Could an endocytic vesicle use its bound CLIP-170 as a rescue factor to prevent the disappearance of the track on which it is traveling?CLIP-170 is indeed considered to be a rescue factor in mammalian cells (Komarova et al., 2002a). Komarova et al. (2002b) have found that in cultured CHO and NRK cells, microtubule dynamics seem to be controlled spatially; catastrophe and rescue occur frequently only near the cell periphery. Although the mechanisms behind catastrophe and rescue are not fully understood, protein factors are required for regulating both events in vivo (Desai and Mitchison, 1997). In CHO cells, a dominant-negative form of CLIP-170 that displaces the endogenous CLIP-170 from microtubule plus ends severely reduces the rescue frequency so that microtubules are more likely to shrink all the way back to the microtubule-organizing center (Komarova et al., 2002a). Moreover, both in vivo and in vitro studies suggest that the rescue activity of CLIP-170 is localized to the NH₂ terminus containing the CAP-Gly motifs (Komarova et al., 2002a; Arnal et al., 2004). How CLIP-170 rescues a shrinking end is not clear. CLIP-170 can promote tubulin oligomerization (Diamantopoulos et al., 1999; Arnal et al., 2004), and it is likely that this property serves to increase the local concentration of tubulin substrate, thereby lowering the entropic barrier for the polymerization reaction. CLIP-170 in mammalian cells has only been found at growing plus ends, most likely as a result of copolymerization with tubulin subunits followed by its release from older segments (Diamantopoulos et al., 1999; Perez et al., 1999; Folker et al., 2005). When a microtubule end shrinks, CLIP-170 falls off. Is there a mechanism to get CLIP-170 close to the depolymerizing end and facilitate its function as a rescue factor? Given the proposed function of Stu2 as a rescue factor for spindle microtubules, one may easily imagine a similar scenario in which vesicle-bound CLIP-170 may be transported to the approaching microtubule end to rescue it from further shrinkage.If vesicle-bound CLIP-170 is transported to the plus end in a manner similar to Stu2, could such transport be mediated by plus end–directed kinesins? Although the kinesin involved in transporting Stu2 toward the microtubule plus end still needs to be identified, detailed image analyses have revealed a role for the Kip2/Tea2 kinesins (members of the kinesin-7 family) in transporting CLIP-170 homologues in fungi (Busch et al., 2004; Carvalho et al., 2004). Bik1 and Tip1 are the CLIP-170 homologues in budding and fission yeasts, respectively, and these proteins are found at microtubule plus ends, where they act as growth-promoting factors or anticatastrophe factors (Berlin et al., 1990; Brunner and Nurse 2000; Carvalho et al., 2004). In both yeasts, the Kip2/Tea2 kinesins bind to and comigrate with the CLIP-170 homologues along the microtubule toward the plus end (Busch et al., 2004; Carvalho et al., 2004). Kinesins have also been implicated in targeting other +TIPs to microtubule plus ends (Jimbo et al., 2002; Maekawa et al., 2003; Zhang et al., 2003; Wu et al., 2006). For example, the mammalian tumor suppressor protein APC (adenomatous polyposis coli) may be targeted to the plus end by KIF3A/KIF3B (a heterotrimeric kinesin II in the kinesin-2 family) as well as by other mechanisms (Jimbo et al., 2002; Nathke 2004; Slep et al., 2005). It will be interesting to see whether a similar transport process for CLIP-170 exists in higher eukaryotic cells. It is possible that such a mechanism would deliver just enough CLIP-170 to the shrinking plus end to initiate rescue. When microtubule growth is resumed, CLIP-170''s intrinsic higher affinity for tubulin subunits and lower affinity for the microtubule wall may allow these proteins to “treadmill” on the growing end (Perez et al., 1999; Folker et al., 2005).The regulation of CLIP-170 activity appears to be rather complex. CLIP-170 is most likely phosphorylated by multiple kinases, including FKBP12–rapamycin-associated protein (mTOR; Choi et al., 2002). Although phosphorylation by mTOR/FKBP 12–rapamycin-associated protein may stimulate CLIP-170''s microtubule binding, phosphorylation by other kinases may cause CLIP-170 to dissociate from microtubules (Rickard and Kreis 1996; Choi et al., 2002). In vivo, CLIP-170 has a closed conformation that is presumably inactive and an open conformation that may interact with microtubules and dynein regulators such as dynactin (Schroer 2004) and LIS1 (Morris et al., 1998; Lansbergen et al., 2004). It is possible that phosphorylation may regulate the conversion between these two forms, but the specific mechanism and the spatial regulation for this conversion have yet to be resolved. If CLIP-170 is indeed released from a membranous cargo to move to the plus end in order to serve as a rescue factor, it would be interesting to know when and/or where such a conformational switch occurs. Finally, other proteins may play redundant roles with CLIP-170 in vesicular trafficking, which may explain why a dramatic defect in vesicle/organelle distribution is not detected when the CLIP-170 level is lowered or when the gene is knocked out (Lansbergen et al., 2004; Akhmanova et al., 2005).+TIPs other than CLIP-170 may play a similar role in rescuing shrinking microtubule tracks. For example, the dynactin complex that links dynein to membranous cargoes and promotes the processive motion of dynein (Schroer 2004) may act as a rescue factor. The p150^Glued subunit of dynactin and CLIP-170 both contain CAP-Gly microtubule-binding motifs at their NH2 termini, although p150^Glued contains one, whereas CLIP-170 contains two such motifs. Dynactin has been shown to behave as a +TIP facilitating the capture of vesicular cargo for minus end–directed transport (Vaughan et al., 1999; 2002). The head domain of the p150^Glued subunit containing the CAP-Gly motif has been shown to promote rescue in vivo in the absence of endogenous CLIP-170, although the effect was much weaker than that caused by the exogenous CLIP-170 head domain (Kamarova et al., 2002a). In vitro studies showed that dynactin may promote nucleation during microtubule assembly (Ligon et al., 2003), which is consistent with it being a potential rescue factor. As shown with CLIP-170, this capacity to bring multiple tubulins together may help to overcome the entropic barrier of the polymerization reaction. Finally, cargo-bound dynactin may also use kinesin to get to the plus end. The p150^Glued subunit of dynactin has been shown to interact directly with the COOH terminus of KAP3, a subunit of the heterotrimeric kinesin II (a member of the kinesin-2 family) that also binds to APC (Jimbo et al., 2002; Deacon et al., 2003; Dell 2003). Although this binding is implicated in dynactin''s role as a cargo adaptor for kinesin II, it is possible, in theory, that a small amount of dynactin may use this connection to move to the plus end.The proposed hypothesis that +TIPs may be released from a membranous cargo to rescue a shrinking microtubule track may apply to both minus and plus end–directed transport. In addition, it is important to point out that this hypothesis does not exclude other mechanisms for rescuing long microtubule tracks. Rescue may occur stochastically, and, sometimes, +TIPs may participate in other ways such as mediating microtubule capture by the actin-rich cortex to stabilize the track (Wen et al., 2004; Galjart 2005). In some situations, microtubule dynamics are modulated by the direct binding of membranous cargo to the growing or shrinking plus ends of microtubules (Waterman-Storer and Salmon, 1998).Currently, the ability of CLIP-170 or other +TIPs to be released from a membranous cargo and to act as a rescue factor for a shrinking microtubule is just a hypothesis. Nevertheless, searching for proteins involved in the communication between a cargo and the approaching shrinking end of its microtubule track is clearly an endeavor worth pursuing.     

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.     

PLASTID MOVEMENT IMPAIRED1 and PLASTID MOVEMENT IMPAIRED1-RELATED1 Mediate Photorelocation Movements of Both Chloroplasts and Nuclei

Organelle movement and positioning play important roles in fundamental cellular activities and adaptive responses to environmental stress in plants. To optimize photosynthetic light utilization, chloroplasts move toward weak blue light (the accumulation response) and escape from strong blue light (the avoidance response). Nuclei also move in response to strong blue light by utilizing the light-induced movement of attached plastids in leaf cells. Blue light receptor phototropins and several factors for chloroplast photorelocation movement have been identified through molecular genetic analysis of Arabidopsis (Arabidopsis thaliana). PLASTID MOVEMENT IMPAIRED1 (PMI1) is a plant-specific C2-domain protein that is required for efficient chloroplast photorelocation movement. There are two PLASTID MOVEMENT IMPAIRED1-RELATED (PMIR) genes, PMIR1 and PMIR2, in the Arabidopsis genome. However, the mechanism in which PMI1 regulates chloroplast and nuclear photorelocation movements and the involvement of PMIR1 and PMIR2 in these organelle movements remained unknown. Here, we analyzed chloroplast and nuclear photorelocation movements in mutant lines of PMI1, PMIR1, and PMIR2. In mesophyll cells, the pmi1 single mutant showed severe defects in both chloroplast and nuclear photorelocation movements resulting from the impaired regulation of chloroplast-actin filaments. In pavement cells, pmi1 mutant plants were partially defective in both plastid and nuclear photorelocation movements, but pmi1pmir1 and pmi1pmir1pmir2 mutant lines lacked the blue light-induced movement responses of plastids and nuclei completely. These results indicated that PMI1 is essential for chloroplast and nuclear photorelocation movements in mesophyll cells and that both PMI1 and PMIR1 are indispensable for photorelocation movements of plastids and thus, nuclei in pavement cells.In plants, organelles move within the cell and become appropriately positioned to accomplish their functions and adapt to the environment (for review, see Wada and Suetsugu, 2004). Light-induced chloroplast movement (chloroplast photorelocation movement) is one of the best characterized organelle movements in plants (Suetsugu and Wada, 2012). Under weak light conditions, chloroplasts move toward light to capture light efficiently (the accumulation response; Zurzycki, 1955). Under strong light conditions, chloroplasts escape from light to avoid photodamage (the avoidance response; Kasahara et al., 2002; Sztatelman et al., 2010; Davis and Hangarter, 2012; Cazzaniga et al., 2013). In most green plant species, these responses are induced primarily by the blue light receptor phototropin (phot) in response to a range of wavelengths from UVA to blue light (approximately 320–500 nm; for review, see Suetsugu and Wada, 2012; Wada and Suetsugu, 2013; Kong and Wada, 2014). Phot-mediated chloroplast movement has been shown in land plants, such as Arabidopsis (Arabidopsis thaliana; Jarillo et al., 2001; Kagawa et al., 2001; Sakai et al., 2001), the fern Adiantum capillus-veneris (Kagawa et al., 2004), the moss Physcomitrella patens (Kasahara et al., 2004), and the liverwort Marchantia polymorpha (Komatsu et al., 2014). Two phots in Arabidopsis, phot1 and phot2, redundantly mediate the accumulation response (Sakai et al., 2001), whereas phot2 primarily regulates the avoidance response (Jarillo et al., 2001; Kagawa et al., 2001; Luesse et al., 2010). M. polymorpha has only one phot that mediates both the accumulation and avoidance responses (Komatsu et al., 2014), although two or more phots mediate chloroplast photorelocation movement in A. capillus-veneris (Kagawa et al., 2004) and P. patens (Kasahara et al., 2004). Thus, duplication and functional diversification of PHOT genes have occurred during land plant evolution, and plants have gained a sophisticated light sensing system for chloroplast photorelocation movement.In general, movements of plant organelles, including chloroplasts, are dependent on actin filaments (for review, see Wada and Suetsugu, 2004). Most organelles common in eukaryotes, such as mitochondria, peroxisomes, and Golgi bodies, use the myosin motor for their movements, but there is no clear evidence that chloroplast movement is myosin dependent (for review, see Suetsugu et al., 2010a). Land plants have innovated a novel actin-based motility system that is specialized for chloroplast movement as well as a photoreceptor system (for review, see Suetsugu et al., 2010a; Wada and Suetsugu, 2013; Kong and Wada, 2014). Chloroplast-actin (cp-actin) filaments, which were first found in Arabidopsis, are short actin filaments specifically localized around the chloroplast periphery at the interface between the chloroplast and the plasma membrane (Kadota et al., 2009). Strong blue light induces the rapid disappearance of cp-actin filaments and then, their subsequent reappearance preferentially at the front region of the moving chloroplasts. This asymmetric distribution of cp-actin filaments is essential for directional chloroplast movement (Kadota et al., 2009; Kong et al., 2013a). The greater the difference in the amount of cp-actin filaments between the front and rear regions of chloroplasts becomes, the faster the chloroplasts move, in which the magnitude of the difference is determined by fluence rate (Kagawa and Wada, 2004; Kadota et al., 2009; Kong et al., 2013a). Strong blue light-induced disappearance of cp-actin filaments is regulated in a phot2-dependent manner before the intensive polymerization of cp-actin filaments at the front region occurs (Kadota et al., 2009; Ichikawa et al., 2011; Kong et al., 2013a). This phot2-dependent response contributes to the greater difference in the amount of cp-actin filaments between the front and rear regions of chloroplasts. Similar behavior of cp-actin filaments has also been observed in A. capillus-veneris (Tsuboi and Wada, 2012) and P. patens (Yamashita et al., 2011).Like chloroplasts, nuclei also show light-mediated movement and positioning (nuclear photorelocation movement) in land plants (for review, see Higa et al., 2014b). In gametophytic cells of A. capillus-veneris, weak light induced the accumulation responses of both chloroplasts and nuclei, whereas strong light induced avoidance responses (Kagawa and Wada, 1993, 1995; Tsuboi et al., 2007). However, in mesophyll cells of Arabidopsis, strong blue light induced both chloroplast and nuclear avoidance responses, but weak blue light induced only the chloroplast accumulation response (Iwabuchi et al., 2007, 2010; Higa et al., 2014a). In Arabidopsis pavement cells, small numbers of tiny plastids were found and showed autofluorescence under the confocal laser-scanning microscopy (Iwabuchi et al., 2010; Higa et al., 2014a). Hereafter, the plastid in the pavement cells is called the pavement cell plastid. Strong blue light-induced avoidance responses of pavement cell plastids and nuclei were induced in a phot2-dependent manner, but the accumulation response was not detected for either organelle (Iwabuchi et al., 2007, 2010; Higa et al., 2014a). In both Arabidopsis and A. capillus-veneris, phots mediate nuclear photorelocation movement, and phot2 mediates the nuclear avoidance response (Iwabuchi et al., 2007, 2010; Tsuboi et al., 2007). The nuclear avoidance response is dependent on actin filaments in both mesophyll and pavement cells of Arabidopsis (Iwabuchi et al., 2010). Recently, it was shown that the nuclear avoidance response relies on cp-actin-dependent movement of pavement cell plastids, where nuclei are associated with pavement cell plastids of Arabidopsis (Higa et al., 2014a). In mesophyll cells, nuclear avoidance response is likely dependent on cp-actin filament-mediated chloroplast movement, because the mutants deficient in chloroplast movement were also defective in nuclear avoidance response (Higa et al., 2014a). Thus, phots mediate both chloroplast (and pavement cell plastid) and nuclear photorelocation movement by regulating cp-actin filaments.Molecular genetic analyses of Arabidopsis mutants deficient in chloroplast photorelocation movement have identified many molecular factors involved in signal transduction and/or motility systems as well as those involved in the photoreceptor system for chloroplast photorelocation movement (and thus, nuclear photorelocation movement; for review, see Suetsugu and Wada, 2012; Wada and Suetsugu, 2013; Kong and Wada, 2014). CHLOROPLAST UNUSUAL POSITIONING1 (CHUP1; Oikawa et al., 2003) and KINESIN-LIKE PROTEIN FOR ACTIN-BASED CHLOROPLAST MOVEMENT (KAC; Suetsugu et al., 2010b) are key factors for generating and/or maintaining cp-actin filaments. Both proteins are highly conserved in land plants and essential for the movement and attachment of chloroplasts to the plasma membrane in Arabidopsis (Oikawa et al., 2003, 2008; Suetsugu et al., 2010b), A. capillus-veneris (Suetsugu et al., 2012), and P. patens (Suetsugu et al., 2012; Usami et al., 2012). CHUP1 is localized on the chloroplast outer membrane and binds to globular and filamentous actins and profilin in vitro (Oikawa et al., 2003, 2008; Schmidt von Braun and Schleiff, 2008). Although KAC is a kinesin-like protein, it lacks microtubule-dependent motor activity but has filamentous actin binding activity (Suetsugu et al., 2010b). An actin-bundling protein THRUMIN1 (THRUM1) is required for efficient chloroplast photorelocation movement (Whippo et al., 2011) and interacts with cp-actin filaments (Kong et al., 2013a). chup1 and kac mutant plants were shown to lack detectable cp-actin filaments (Kadota et al., 2009; Suetsugu et al., 2010b; Ichikawa et al., 2011; Kong et al., 2013a). Similarly, cp-actin filaments were rarely detected in thrum1 mutant plants (Kong et al., 2013a), indicating that THRUM1 also plays an important role in maintaining cp-actin filaments.Other proteins J-DOMAIN PROTEIN REQUIRED FOR CHLOROPLAST ACCUMULATION RESPONSE1 (JAC1; Suetsugu et al., 2005), WEAK CHLOROPLAST MOVEMENT UNDER BLUE LIGHT1 (WEB1; Kodama et al., 2010), and PLASTID MOVEMENT IMPAIRED2 (PMI2; Luesse et al., 2006; Kodama et al., 2010) are involved in the light regulation of cp-actin filaments and chloroplast photorelocation movement. JAC1 is an auxilin-like J-domain protein that mediates the chloroplast accumulation response through its J-domain function (Suetsugu et al., 2005; Takano et al., 2010). WEB1 and PMI2 are coiled-coil proteins that interact with each other (Kodama et al., 2010). Although web1 and pmi2 were partially defective in the avoidance response, the jac1 mutation completely suppressed the phenotype of web1 and pmi2, suggesting that the WEB1/PMI2 complex suppresses JAC1 function (i.e. the accumulation response) under strong light conditions (Kodama et al., 2010). Both web1 and pmi2 showed impaired disappearance of cp-actin filaments in response to strong blue light (Kodama et al., 2010). However, the exact molecular functions of these proteins are unknown.In this study, we characterized mutant plants deficient in the PMI1 gene and two homologous genes PLASTID MOVEMENT IMPAIRED1-RELATED1 (PMIR1) and PMIR2. PMI1 was identified through molecular genetic analyses of pmi1 mutants that showed severe defects in chloroplast accumulation and avoidance responses (DeBlasio et al., 2005). PMI1 is a plant-specific C2-domain protein (DeBlasio et al., 2005; Zhang and Aravind, 2010), but its roles and those of PMIRs in cp-actin-mediated chloroplast and nuclear photorelocation movements remained unclear. Thus, we analyzed chloroplast and nuclear photorelocation movements in the single, double, and triple mutants of pmi1, pmir1, and pmir2.     

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.     

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.