Quantitative analysis of microtubule orientation in interdigitated leaf pavement cells

Machine learning and feature analysis of the cortical microtubule organization of Arabidopsis cotyledon pavement cells

The measurement of cytoskeletal features can provide valuable insights into cell biology. In recent years, digital image analysis of cytoskeletal features has become an important research tool for quantitative evaluation of cytoskeleton organization. In this study, we examined the utility of a supervised machine learning approach with digital image analysis to distinguish different cellular organizational patterns. We focused on the jigsaw puzzle-shaped pavement cells of Arabidopsis thaliana. Measurements of three features of cortical microtubules in these cells (parallelness, density, and the coefficient of variation of the intensity distribution of fluorescently labeled cytoskeletons [as an indicator of microtubule bundling]) were obtained from microscopic images. A random forest machine learning model was then used with these images to differentiate mutant and wild type, and Taxol-treated and control cells. Using these three metrics, we were able to distinguish wild type from bpp125 triple mutant cells, with approximately 80% accuracy; classification accuracy was 88% for control and Taxol-treated cells. Different features contributed most to the classification, namely, coefficient of variation for the wild-type/mutant cells and parallelness for the Taxol-treated/control cells. The random forest method used enabled quantitative evaluation of the contribution of features to the classification, and partial dependence plots showed the relationships between metric values and classification accuracy. While further improvements to the method are needed, our small-scale analysis shows the potential for this approach in large-scale screening analyses.

Arabidopsis KASH Proteins SINE1 and SINE2 Are Involved in Microtubule Reorganization During ABA-Induced Stomatal Closure

Abscisic acid (ABA) induces stomatal closure by utilizing complex signaling mechanisms, allowing for sessile plants to respond rapidly to ever-changing environmental conditions. ABA regulates the activity of plasma membrane ion channels and calcium-dependent protein kinases, Ca2+ oscillations, and reactive oxygen species (ROS) concentrations. Throughout ABA-induced stomatal closure, the cytoskeleton undergoes dramatic changes that appear important for efficient closure. However, the precise role of this cytoskeletal reorganization in stomatal closure and the nature of its regulation are unknown. We have recently shown that the plant KASH proteins SINE1 and SINE2 are connected to actin organization during ABA-induced stomatal closure but their role in microtubule (MT) organization remains to be investigated. We show here that depolymerizing MTs using oryzalin can restore ABA-induced stomatal closure deficits in sine1-1 and sine2-1 mutants. GFP-MAP4-visualized MT organization is compromised in sine1-1 and sine2-1 mutants during ABA-induced stomatal closure. Loss of SINE1 or SINE2 results in loss of radially organized MT patterning in open guard cells, aberrant MT organization during stomatal closure, and an overall decrease in the number of MT filaments or bundles. Thus, SINE1 and SINE2 are necessary for establishing MT patterning and mediating changes in MT rearrangement, which is required for ABA-induced stomatal closure.

Mechanical Stress Initiates and Sustains the Morphogenesis of Wavy Leaf Epidermal Cells

Pavement cells form wavy interlocking patterns in the leaf epidermis of many plants. We use computational mechanics to simulate the morphogenetic process based on microtubule organization and cell wall chemistry. Based on the in silico simulations and experimental evidence, we suggest that a multistep process underlies the morphogenesis of pavement cells. The in silico model predicts alternatingly located, feedback-augmented mechanical heterogeneity of the periclinal and anticlinal walls. It suggests that the emergence of waves is created by a stiffening of the emerging indented sides, an effect that matches cellulose and de-esterified pectin patterns in the cell wall. Further, conceptual evidence for mechanical buckling of the cell walls is provided, a mechanism that has the potential to initiate wavy patterns de novo and may precede chemical and geometrical symmetry breaking.

Mitochondrial dynamics and segregation during the asymmetric division of Arabidopsis zygotes

The zygote is the first cell of a multicellular organism. In most angiosperms, the zygote divides asymmetrically to produce an embryo-precursor apical cell and a supporting basal cell. Zygotic division should properly segregate symbiotic organelles, because they cannot be synthesized de novo. In this study, we revealed the real-time dynamics of the principle source of ATP biogenesis, mitochondria, in Arabidopsis thaliana zygotes using live-cell observations and image quantifications. In the zygote, the mitochondria formed the extended structure associated with the longitudinal array of actin filaments (F-actins) and were polarly distributed along the apical-basal axis. The mitochondria were then temporally fragmented during zygotic division, and the resulting apical cells inherited mitochondria at higher concentration compared to the basal cells. Further observation of postembryonic organs showed that these mitochondrial behaviours are characteristic of the zygote. Overall, our results showed that the zygote has spatiotemporal regulation that unequally distributes the mitochondria.

Basic Proline-Rich Protein-Mediated Microtubules Are Essential for Lobe Growth and Flattened Cell Geometry

Complex cell shapes are generated first by breaking symmetry, and subsequent polar growth. Localized bending of anticlinal walls initiates lobe formation in the epidermal pavement cells of cotyledons and leaves, but how the microtubule cytoskeleton mediates local cell growth, and how plant pavement cells benefit from adopting jigsaw puzzle-like shapes, are poorly understood. In Arabidopsis (Arabidopsis thaliana), the basic Pro-rich protein (BPP) microtubule-associated protein family comprises seven members. We analyzed lobe morphogenesis in cotyledon pavement cells of a BPP1;BPP2;BPP5 triple knockout mutant. New image analysis methods (MtCurv and BQuant) showed that anticlinal microtubule bundles were significantly reduced and cortical microtubules that fan out radially across the periclinal wall did not enrich at the convex side of developing lobes. Despite these microtubule defects, new lobes were initiated at the same frequency as in wild-type cells, but they did not expand into well-defined protrusions. Eventually, mutant cells formed nearly polygonal shapes and adopted concentric microtubule patterns. The mutant periclinal cell wall bulged outward. The radius of the calculated inscribed circle of the pavement cells, a proposed proxy for maximal stress in the cell wall, was consistently larger in the mutant cells during cotyledon development, and correlated with an increase in cell height. These bpp mutant phenotypes provide genetic and cell biological evidence that initiation and growth of lobes are distinct morphogenetic processes, and that interdigitated cell geometry effectively suppresses large outward bulging of pavement cells.

Asymmetric distribution of cytokinins determines root hydrotropism in Arabidopsis thaliana

The phenomenon of plant root tips sensing moisture gradient in soil and growing towards higher water potential is designated as root hydrotropism, which is critical for plants to survive when water is a limited factor. Molecular mechanisms regulating such a fundamental process, however, are largely unknown. Here we report our identification that cytokinins are key signaling molecules directing root growth orientation in a hydrostimulation (moisture gradient) condition. Lower water potential side of the root tip shows more cytokinin response relative to the higher water potential side. Consequently, two cytokinin downstream type-A response regulators, ARR16 and ARR17, were found to be up-regulated at the lower water potential side, causing increased cell division in the meristem zone, which allows the root to bend towards higher water potential side. Genetic analyses indicated that various cytokinin biosynthesis and signaling mutants, including the arr16 arr17 double mutant, are significantly less responsive to hydrostimulation. Consistently, treatments with chemical inhibitors interfering with either cytokinin biosynthesis or cell division completely abolished root hydrotropic response. Asymmetrically induced expression of ARR16 or ARR17 effectively led to root bending in both wild-type and miz1, a previously known hydrotropism-defective mutant. These data demonstrate that asymmetric cytokinin distribution is a primary determinant governing root hydrotropism.

Microtubule-associated protein IQ67 DOMAIN5 regulates morphogenesis of leaf pavement cells in Arabidopsis thaliana

Plant microtubules form a highly dynamic intracellular network with important roles for regulating cell division, cell proliferation, and cell morphology. Their organization and dynamics are co-ordinated by various microtubule-associated proteins (MAPs) that integrate environmental and developmental stimuli to fine-tune and adjust cytoskeletal arrays. IQ67 DOMAIN (IQD) proteins recently emerged as a class of plant-specific MAPs with largely unknown functions. Here, using a reverse genetics approach, we characterize Arabidopsis IQD5 in terms of its expression domains, subcellular localization, and biological functions. We show that IQD5 is expressed mostly in vegetative tissues, where it localizes to cortical microtubule arrays. Our phenotypic analysis of iqd5 loss-of-function lines reveals functions of IQD5 in pavement cell (PC) shape morphogenesis. Histochemical analysis of cell wall composition further suggests reduced rates of cellulose deposition in anticlinal cell walls, which correlate with reduced anisotropic expansion. Lastly, we demonstrate IQD5-dependent recruitment of calmodulin calcium sensors to cortical microtubule arrays and provide first evidence for important roles for calcium in regulation of PC morphogenesis. Our work identifies IQD5 as a novel player in PC shape regulation and, for the first time, links calcium signaling to developmental processes that regulate anisotropic growth in PCs.

Microtubule-associated protein IQ67 DOMAIN5 regulates morphogenesis of leaf pavement cells in Arabidopsis thaliana

Plant microtubules form a highly dynamic intracellular network with important roles for regulating cell division, cell proliferation, and cell morphology. Their organization and dynamics are co-ordinated by various microtubule-associated proteins (MAPs) that integrate environmental and developmental stimuli to fine-tune and adjust cytoskeletal arrays. IQ67 DOMAIN (IQD) proteins recently emerged as a class of plant-specific MAPs with largely unknown functions. Here, using a reverse genetics approach, we characterize Arabidopsis IQD5 in terms of its expression domains, subcellular localization, and biological functions. We show that IQD5 is expressed mostly in vegetative tissues, where it localizes to cortical microtubule arrays. Our phenotypic analysis of iqd5 loss-of-function lines reveals functions of IQD5 in pavement cell (PC) shape morphogenesis. Histochemical analysis of cell wall composition further suggests reduced rates of cellulose deposition in anticlinal cell walls, which correlate with reduced anisotropic expansion. Lastly, we demonstrate IQD5-dependent recruitment of calmodulin calcium sensors to cortical microtubule arrays and provide first evidence for important roles for calcium in regulation of PC morphogenesis. Our work identifies IQD5 as a novel player in PC shape regulation and, for the first time, links calcium signaling to developmental processes that regulate anisotropic growth in PCs.

Quantitative and Comparative Analysis of Global Patterns of (Microtubule) Cytoskeleton Organization with CytoskeletonAnalyzer2D

The microtubule cytoskeleton plays important roles in cell morphogenesis. To investigate the mechanisms of cytoskeletal organization, for example, during growth or development, in genetic studies, or in response to environmental stimuli, image analysis tools for quantitative assessment are needed. Here, we present a method for texture measure-based quantification and comparative analysis of global microtubule cytoskeleton patterns and subsequent visualization of output data. In contrast to other approaches that focus on the extraction of individual cytoskeletal fibers and analysis of their orientation relative to the growth axis, CytoskeletonAnalyzer2D quantifies cytoskeletal organization based on the analysis of local binary patterns. CytoskeletonAnalyzer2D thus is particularly well suited to study cytoskeletal organization in cells where individual fibers are difficult to extract or which lack a clearly defined growth axis, such as leaf epidermal pavement cells. The tool is available as ImageJ plugin and can be combined with publicly available software and tools, such as R and Cytoscape, to visualize similarity networks of cytoskeletal patterns.

The Microtubule-Associated Protein IQ67 DOMAIN5 Modulates Microtubule Dynamics and Pavement Cell Shape

The dynamic arrangement of cortical microtubules (MTs) plays a pivotal role in controlling cell growth and shape formation in plants, but the mechanisms by which cortical MTs are organized to regulate these processes are not well characterized. In particular, the dynamic behavior of cortical MTs is critical for their spatial organization, yet the molecular mechanisms controlling MT dynamics remain poorly understood. In this study, we used the puzzle piece-shaped pavement cells of Arabidopsis (Arabidopsis thaliana) leaves as a model system in which to study cortical MT organization. We isolated an ethyl methanesulfonate mutant with reduced interdigitation of pavement cells in cotyledons. This line carried a mutation in IQ67 DOMAIN5 (IQD5), which encodes a member of the plant-specific IQ motif protein family. Live-cell imaging and biochemical analyses demonstrated that IQD5 binds to MTs and promotes MT assembly. MT-depolymerizing drug treatment and in vivo MT dynamics assays suggested that IQD5 functions to stabilize MTs. Hence, our findings provide genetic, cell biological, and biochemical evidence that IQD5 regulates MT dynamics that affect MT organization and subsequent cell shape formation.

A new kymogram-based method reveals unexpected effects of marker protein expression and spatial anisotropy of cytoskeletal dynamics in plant cell cortex

BACKGROUND

Cytoskeleton can be observed in live plant cells in situ with high spatial and temporal resolution using a combination of specific fluorescent protein tag expression and advanced microscopy methods such as spinning disc confocal microscopy (SDCM) or variable angle epifluorescence microscopy (VAEM). Existing methods for quantifying cytoskeletal dynamics are often either based on laborious manual structure tracking, or depend on costly commercial software. Current automated methods also do not readily allow separate measurements of structure lifetime, lateral mobility, and spatial anisotropy of these parameters.

RESULTS

We developed a new freeware-based, operational system-independent semi-manual technique for analyzing VAEM or SDCM data, QuACK (Quantitative Analysis of Cytoskeletal Kymograms), and validated it on data from Arabidopsis thaliana fh1 formin mutants, previously shown by conventional methods to exhibit altered actin and microtubule dynamics compared to the wild type. Besides of confirming the published mutant phenotype, QuACK was used to characterize surprising differential effects of various fluorescent protein tags fused to the Lifeact actin probe on actin dynamics in A. thaliana cotyledon epidermis. In particular, Lifeact-YFP slowed down actin dynamics compared to Lifeact-GFP at marker expression levels causing no macroscopically noticeable phenotypic alterations, although the two fluorophores are nearly identical. We could also demonstrate the expected, but previously undocumented, anisotropy of cytoskeletal dynamics in elongated epidermal cells of A. thaliana petioles and hypocotyls.

CONCLUSIONS

Our new method for evaluating plant cytoskeletal dynamics has several advantages over existing techniques. It is intuitive, rapid compared to fully manual approaches, based on the free ImageJ software (including macros we provide here for download), and allows measurement of multiple parameters. Our approach was already used to document unexpected differences in actin mobility in transgenic A. thaliana expressing Lifeact fusion proteins with different fluorophores, highlighting the need for cautious interpretation of experimental results, as well as to reveal hitherto uncharacterized anisotropy of cytoskeletal mobility in elongated plant cells.

Cell wall accumulation of fluorescent proteins derived from a trans-Golgi cisternal membrane marker and paramural bodies in interdigitated Arabidopsis leaf epidermal cells

In most dicotyledonous plants, leaf epidermal pavement cells develop jigsaw puzzle-like shapes during cell expansion. The rapid growth and complicated cell shape of pavement cells is suggested to be achieved by targeted exocytosis that is coordinated with cytoskeletal rearrangement to provide plasma membrane and/or cell wall materials for lobe development during their morphogenesis. Therefore, visualization of membrane trafficking in leaf pavement cells should contribute an understanding of the mechanism of plant cell morphogenesis. To reveal membrane trafficking in pavement cells, we observed monomeric red fluorescent protein-tagged rat sialyl transferases, which are markers of trans-Golgi cisternal membranes, in the leaf epidermis of Arabidopsis thaliana. Quantitative fluorescence imaging techniques and immunoelectron microscopic observations revealed that accumulation of the red fluorescent protein occurred mostly in the curved regions of pavement cell borders and guard cell ends during leaf expansion. Transmission electron microscopy observations revealed that apoplastic vesicular membrane structures called paramural bodies were more frequent beneath the curved cell wall regions of interdigitated pavement cells and guard cell ends in young leaf epidermis. In addition, pharmacological studies showed that perturbations in membrane trafficking resulted in simple cell shapes. These results suggested possible heterogeneity of the curved regions of plasma membranes, implying a relationship with pavement cell morphogenesis.

Exogenous Cellulase Switches Cell Interdigitation to Cell Elongation in an RIC1-dependent Manner in Arabidopsis thaliana Cotyledon Pavement Cells

Pavement cells in cotyledons and true leaves exhibit a jigsaw puzzle-like morphology in most dicotyledonous plants. Among the molecular mechanisms mediating cell morphogenesis, two antagonistic Rho-like GTPases regulate local cell outgrowth via cytoskeletal rearrangements. Analyses of several cell wall-related mutants suggest the importance of cell wall mechanics in the formation of interdigitated patterns. However, how these factors are integrated is unknown. In this study, we observed that the application of exogenous cellulase to hydroponically grown Arabidopsis thaliana cotyledons switched the interdigitation of pavement cells to the production of smoothly elongated cells. The cellulase-induced inhibition of cell interdigitation was not observed in a RIC1 knockout mutant. This gene encodes a Rho-like GTPase-interacting protein important for localized cell growth suppression via microtubule bundling on concave cell interfaces. Additionally, to characterize pavement cell morphologies, we developed a mathematical model that considers the balance between cell and cell wall growth, restricted global cell growth orientation, and regulation of local cell outgrowth mediated by a Rho-like GTPase-cytoskeleton system. Our computational simulations fully support our experimental observations, and suggest that interdigitated patterns form because of mechanical buckling in the absence of Rho-like GTPase-dependent regulation of local cell outgrowth. Our model clarifies the cell wall mechanics influencing pavement cell morphogenesis.

A Theoretical Model of Jigsaw-Puzzle Pattern Formation by Plant Leaf Epidermal Cells

Plant leaf epidermal cells exhibit a jigsaw puzzle-like pattern that is generated by interdigitation of the cell wall during leaf development. The contribution of two ROP GTPases, ROP2 and ROP6, to the cytoskeletal dynamics that regulate epidermal cell wall interdigitation has already been examined; however, how interactions between these molecules result in pattern formation remains to be elucidated. Here, we propose a simple interface equation model that incorporates both the cell wall remodeling activity of ROP GTPases and the diffusible signaling molecules by which they are regulated. This model successfully reproduces pattern formation observed in vivo, and explains the counterintuitive experimental results of decreased cellulose production and increased thickness. Our model also reproduces the dynamics of three-way cell wall junctions. Therefore, this model provides a possible mechanism for cell wall interdigitation formation in vivo.