A Two-Hybrid-Receptor Assay Demonstrates Heteromer Formation as Switch-On for Plant Immune Receptors

Receptor kinases sense extracellular signals and trigger intracellular signaling and physiological responses. However, how does signal binding to the extracellular domain activate the cytoplasmic kinase domain? Activation of the plant immunoreceptor Flagellin sensing2 (FLS2) by its bacterial ligand flagellin or the peptide-epitope flg22 coincides with rapid complex formation with a second receptor kinase termed brassinosteroid receptor1 associated kinase1 (BAK1). Here, we show that the receptor pair of FLS2 and BAK1 is also functional when the roles of the complex partners are reversed by swapping their cytosolic domains. This reciprocal constellation prevents interference by redundant partners that can partially substitute for BAK1 and demonstrates that formation of the heteromeric complex is the molecular switch for transmembrane signaling. A similar approach with swaps between the Elongation factor-Tu receptor and BAK1 also resulted in a functional receptor/coreceptor pair, suggesting that a “two-hybrid-receptor assay” is of more general use for studying heteromeric receptor complexes.Cell surface receptors are chemical sensors, often with an exquisite specificity and sensitivity, which detect extracellular signals and initiate corresponding intracellular response programs. Many of these receptors are transmembrane proteins with an extracellular ligand-binding domain and an intracellular protein kinase domain. Higher plants, such as Arabidopsis (Arabidopsis thaliana), have several hundred genes encoding receptor like kinases (Shiu and Bleecker, 2001; Shiu and Li, 2004). How are these receptor like kinases activated by their ligands, and how do they initiate a subsequent intracellular signaling cascade? In our work, we used the leucine-rich repeat receptor kinase (LRR-RK) Flagellin sensing2 (FLS2), which specifically detects bacterial flagellin or its peptide epitope flg22 at subnanomolar concentrations (Gomez-Gomez and Boller, 2000; Gómez-Gómez et al., 2001; Chinchilla et al., 2006). FLS2 undergoes heteromeric complex formation with brassinosteroid receptor1 associated kinase1 (BAK1) within seconds after application of the flagellin-derived peptide ligand flg22 (Chinchilla et al., 2007; Heese et al., 2007; Schulze et al., 2010). Thus, BAK1 might act as a coreceptor of FLS2. However, as previously observed (Chinchilla et al., 2007; Roux et al., 2011), FLS2 is still functional in the absence of BAK1, although with a reduced efficiency. This raises the question whether the ligand-induced heteromeric complex has merely an enhancing effect or whether association with BAK1 or a functional substitute acts as the essential switch-on for transmembrane signaling of FLS2. BAK1 is one of the five members that form the somatic embryogenesis receptor kinase (SERK) family (Albrecht et al., 2008), and other members of this family might partially substitute for BAK1 (Roux et al., 2011). However, a rigorous genetic approach to delineate the role of these potential substitutes is not feasible because triple mutants (serk1 serk3 serk4) and quadruple mutants (serk1 serk2 serk3 serk4) exhibit severe general phenotypes of dwarfing or even lethality at the early embryo stage (He et al., 2007; Gou et al., 2012) that might be due to the important role of SERKs in plant developmental processes (Li et al., 2002; Nam and Li, 2002).To address the role of the heteromeric complex with BAK1 in the absence of other interfering SERKs, we took a two-hybrid-receptor approach based on the premise that the apoplastic and cytoplasmic domains of FLS2 and BAK1 function in a modular manner. A heteromeric complex might thus also form and function when the roles of FLS2 and BAK1 are reversed by reciprocal swapping of their cytoplasmic protein kinase domains (Fig. 1A, schematic view).Open in a separate windowFigure 1.Heteromeric complex formation of FLS with BAK1 switches on flagellin-dependent transmembrane signaling. A, Model for the flg22-dependent heteromeric receptor complex. Schematic representation of FLS2 (blue), BAK1 (red), and the chimeric receptor constructs Ftm-B and Btm-F. Ftm-B comprises the extracellular and the transmembrane domains of FLS2 (blue) and the cytoplasmic domain of BAK1 (red). The receptor chimera Btm-F represents the reciprocal construct with the cytoplasmic part of BAK1 replaced by that of FLS2. The cytoplasmic domain of FLS2 was C-terminally tagged with a GFP. B and C, Functional comparison of native FLS2 and BAK1 (B) with the two hybrid receptors Ftm-B and Btm-F (C). The experiments show luciferase activity in Arabidopsis fls2 bak1-4 double mutant protoplasts cotransformed with pFRK1::Luciferase as a reporter and the receptor constructs indicated. At 0 h (dashed line) protoplasts were treated with 100 nm of flg22 (black diamonds) or 100 nm of the inactive analog flg22^Atum (white circles). Light emission of the protoplasts was measured with a luminometer. Values represent averages and sds of three replicates. Data shown are representative for at least three independent repetitions of the experiments with all constructs. D and E, Dose-response relationship for flg22-dependent induction of pFRK1::Luciferase in protoplasts expressing the receptor constructs indicated. Values represent increase in luciferase activity after 5 h of treatment as percentage of the increase observed with FLS2 plus BAK1 treated with saturating doses of greater than or equal to 10 nm flg22. Comparison of (half-maximal stimulation values in D and E shows that cells coexpressing Ftm-B plus Btm-F responded at least as sensitive to flg22 as cells coexpressing FLS2 plus BAK1. A combination of Ftm-B with the kinase dead version Btm-F_KD was not functional, even when treated with 1,000 nm flg22. LU, Light units.     

Ethylene Signaling Regulates Accumulation of the FLS2 Receptor and Is Required for the Oxidative Burst Contributing to Plant Immunity

Reactive oxygen species (ROS) are potent signal molecules rapidly generated in response to stress. Detection of pathogen-associated molecular patterns induces a transient apoplastic ROS through the function of the NADPH respiratory burst oxidase homologs D (RbohD). However, little is known about the regulation of pathogen-associated molecular pattern-elicited ROS or its role in plant immunity. We investigated ROS production triggered by bacterial flagellin (flg22) in Arabidopsis (Arabidopsis thaliana). The oxidative burst was diminished in ethylene-insensitive mutants. Flagellin Sensitive2 (FLS2) accumulation was reduced in etr1 and ein2, indicating a requirement of ethylene signaling for FLS2 expression. Multiplication of virulent bacteria was enhanced in Arabidopsis lines displaying altered ROS production at early but not late stages of infection, suggesting an impairment of preinvasive immunity. Stomatal closure, a mechanism used to reduce bacterial entry into plant tissues, was abolished in etr1, ein2, and rbohD mutants. These results point to the importance of flg22-triggered ROS at an early stage of the plant immune response.A rapid and transient increase in reactive oxygen species (ROS), termed an “oxidative burst,” is often associated with responses to abiotic and biotic stresses and could trigger changes in stomatal aperture or programmed cell death in defense against pathogens (Kwak et al., 2003; Torres and Dangl, 2005). ROS production can occur extracellularly through activities of plasma membrane-resident NADPH oxidases (Kangasjärvi et al., 2005; Torres and Dangl, 2005). In plants, Rboh proteins, which are homologs of mammalian NADPH oxidase 2, were shown to be the predominant mediators of apoplastic ROS production (Torres et al., 1998; Galletti et al., 2008). Respiratory burst oxidase homologs D and F (RbohD and RbohF) were identified by mutation to be the responsible oxidases in Arabidopsis (Arabidopsis thaliana) defense responses (Torres et al., 2002). While most ROS generated in response to avirulent Pseudomonas syringae bacteria and Hyaloperonospora oomycete pathogens depend on RbohD function, the induced cell death response by these pathogens appears to be mostly regulated by RbohF. Cell death provoked upon infection with the necrotizing fungus Alternaria, however, is under the control of RbohD (Pogány et al., 2009). The contribution of NADPH oxidases to plant immunity was also described in barley (Hordeum vulgare) and tobacco (Nicotiana benthamiana), where resistance to powdery mildew fungi and the oomycete Phytophthora infestans, respectively, was dependent on Rboh functions (Yoshioka et al., 2003; Trujillo et al., 2006).An early layer of active plant defense is mediated by pattern recognition receptors, which sense microbes according to conserved constituents, so-called pathogen-associated molecular patterns (PAMPs). These initiate a plethora of defense responses referred to as PAMP-triggered immunity (Boller and Felix, 2009). The Arabidopsis receptor kinase Flagellin Sensitive2 (FLS2) recognizes and physically interacts with flg22, the elicitor-active epitope of bacterial flagellin (Felix et al., 1999; Gomez-Gomez and Boller, 2000; Chinchilla et al., 2006). FLS2 is plasma membrane localized and expressed throughout the plant (Robatzek et al., 2006). FLS2 requires the receptor kinase BRI1-Associated Kinase1 (BAK1), which forms a heteromeric complex upon flg22 binding (Chinchilla et al., 2007). Subsequently, a rapid and transient flg22-stimulated oxidative burst occurs that is dependent on RbohD (Zhang et al., 2007). In addition, flg22 triggers early responses, such as ethylene biosynthesis, activation of mitogen-activated protein (MAP) kinase cascades, and changes in gene expression (Felix et al., 1999; Asai et al., 2002; Zipfel et al., 2004). Late flg22 responses include the accumulation of salicylic acid (SA), callose deposition, and an arrest of seedling growth (Gomez-Gomez et al., 1999; Mischina and Zeier, 2007). This collectively contributes to plant immunity (Zipfel et al., 2004; Melotto et al., 2006).Little is known about the regulatory components of FLS2-activated early flg22 responses and their relevance in plant resistance to pathogens. Here, we investigated flg22-triggered ROS production in Arabidopsis seedlings and have identified ethylene signaling as a critical component of the oxidative burst in response to flg22, partly through promoting the accumulation of FLS2. We further provide evidence that the flg22-triggered oxidative burst is required for resistance to bacterial infection at the point of pathogen entry through stomata.     

Genetics of Sexual Development: An Evolutionary Playground for Fish

Teleost fishes are the most species-rich clade of vertebrates and feature an overwhelming diversity of sex-determining mechanisms, classically grouped into environmental and genetic systems. Here, we review the recent findings in the field of sex determination in fish. In the past few years, several new master regulators of sex determination and other factors involved in sexual development have been discovered in teleosts. These data point toward a greater genetic plasticity in generating the male and female sex than previously appreciated and implicate novel gene pathways in the initial regulation of the sexual fate. Overall, it seems that sex determination in fish does not resort to a single genetic cascade but is rather regulated along a continuum of environmental and heritable factors.IN contrast to mammals and birds, cold-blooded vertebrates, and among them teleost fishes in particular, show a variety of strategies for sexual reproduction (Figure 1), ranging from unisexuality (all-female species) to hermaphroditism (sequential, serial, and simultaneous, including outcrossing and selfing species) to gonochorism (two separate sexes at all life stages). The underlying phenotypes are regulated by a variety of sex determination (SD) mechanisms that have classically been divided into two main categories: genetic sex determination (GSD) and environmental sex determination (ESD) (Figure 2).Open in a separate windowFigure 1Reproductive strategies in fish. Fish can be grouped according to their reproductive strategy into unisexuals, hermaphrodites, and gonochorists. Further subdivisions of these three categories are shown with pictures of species exemplifying the strategies. Fish images: Amphiprion clarkii courtesy of Sara Mae Stieb; Hypoplectrus nigricans courtesy of Oscar Puebla; Scarus ferrugineus courtesy of Moritz Muschick; Astatotilapia burtoni courtesy of Anya Theis; Poecilia formosa and Kryptolebias marmoratus courtesy of Manfred Schartl; Trimma sp. courtesy of Rick Winterbottom [serial hermaphroditism has been described in several species of the genus Trimma (Kuwamura and Nakashima 1998; Sakurai et al. 2009; and references therein)].Open in a separate windowFigure 2Sex-determining mechanisms in fish. Sex-determining systems in fish have been broadly classified into environmental and genetic sex determination. For both classes, the currently described subsystems are shown.Environmental factors impacting sex determination in fish are water pH, oxygen concentration, growth rate, density, social state, and, most commonly, temperature (for a detailed review on ESD see, e.g., Baroiller et al. 2009b and Stelkens and Wedekind 2010). As indicated in Figure 2, GSD systems in fish compose a variety of different mechanisms and have been reviewed in detail elsewhere (e.g., Devlin and Nagahama 2002; Volff et al. 2007).The GSD systems that have received the most scientific attention so far are those involving sex chromosomes, which either may be distinguishable cytologically (heteromorphic) or appear identical (homomorphic). In both cases, one sex is heterogametic (possessing two different sex chromosomes and hence producing two types of gametes) and the other one homogametic (a genotype with two copies of the same sex chromosome, producing only one type of gamete). A male-heterogametic system is called an XX-XY system, and female-heterogametic systems are denoted as ZZ-ZW. Both types of heterogamety exist in teleosts and are even found side by side in closely related species [e.g., tilapias (Cnaani et al. 2008), ricefishes (Takehana et al. 2008), or sticklebacks (Ross et al. 2009)]; for more details on the phylogenetic distribution of GSD mechanisms in teleost fish, see Mank et al. (2006). Note that sex chromosomes in fish are mostly homomorphic and not differentiated (Ohno 1974), which is in contrast to the degenerated Y and W chromosomes in mammals (Graves 2006) and birds (Takagi and Sasaki 1974), respectively. This is one possible explanation for the viable combination of different sex chromosomal systems within a single species or population of fish (Parnell and Streelman 2013) and could be a mechanistic reason why sex chromosome turnovers occur easily and frequently in this group (Mank and Avise 2009). Additionally, fish can have more complex sex chromosomal systems involving more than one chromosome pair (see Figure 2). Even within a single fish species, more than two sex chromosomes may occur at the same time, or more than two types of sex chromosomes may co-exist in the same species (Schultheis et al. 2006; Cioffi et al. 2013), which can sometimes be due to chromosome fusions (Kitano and Peichel 2012).Detailed insights on the gene level for GSD/sex chromosomal systems are currently available for only a limited number of fish species, and all but one of these cases involve a rather simple genetic system with male heterogamety and one major sex determiner (see below). The only exception is the widely used model species zebrafish (Danio rerio), which has a polyfactorial SD system implicating four different chromosomes (chromosomes 3, 4, 5, and 16) (Bradley et al. 2011; Anderson et al. 2012) and also environmental cues (Shang et al. 2006).In this review, we focus on newly described genetic sex-determining systems and possible mechanisms allowing their emergence in fishes, which are the most successful group of vertebrates with ∼30,000 species.     

Conformational States of Melittin at a Bilayer Interface

The distribution of peptide conformations in the membrane interface is central to partitioning energetics. Molecular-dynamics simulations enable characterization of in-membrane structural dynamics. Here, we describe melittin partitioning into dioleoylphosphatidylcholine lipids using CHARMM and OPLS force fields. Although the OPLS simulation failed to reproduce experimental results, the CHARMM simulation reported was consistent with experiments. The CHARMM simulation showed melittin to be represented by a narrow distribution of folding states in the membrane interface.Unstructured peptides fold into the membrane interface because partitioned hydrogen-bonded peptide bonds are energetically favorable compared to free peptide bonds (1–3). This folding process is central to the mechanisms of antimicrobial and cell-penetrating peptides, as well as to lipid interactions and stabilities of larger membrane proteins (4). The energetics of peptide partitioning into membrane interfaces can be described by a thermodynamic cycle (Fig. 1). State A is a theoretical state representing the fully unfolded peptide in water, B is the unfolded peptide in the membrane interface, C is the peptide in water, and D is the folded peptide in the membrane. The population of peptides in solution (State C) is best described as an ensemble of folded and unfolded conformations, whereas the population of peptides in State D generally is assumed to have a single, well-defined helicity, as shown in Fig. 1 A (5). Given that, in principle, folding in solution and in the membrane interface should follow the same basic rules, peptides in state D could reasonably be assumed to also be an ensemble. A fundamental question (5) is therefore whether peptides in state D can be correctly described as having a single helicity. Because differentiating an ensemble of conformations and a single conformation may be an impossible experimental task (5), molecular-dynamics (MD) simulations provide a unique high-resolution view of the phenomenon.Open in a separate windowFigure 1Thermodynamic cycles for peptide partitioning into a membrane interface. States A and B correspond to the fully unfolded peptide in solution and membrane interface, respectively. The folded peptide in solution is best described as an ensemble of unfolded and folded conformations (State C). State D is generally assumed to be one of peptides with a narrow range of conformations, but the state could actually be an ensemble of states as in the case of State C.Melittin is a 26-residue, amphipathic peptide that partitions strongly into membrane interfaces and therefore has become a model system for describing folding energetics (3,6–8). Here, we describe the structural dynamics of melittin in a dioleoylphosphatidylcholine (DOPC) bilayer by means of two extensive MD simulations using two different force fields.We extended a 12-ns equilibrated melittin-DOPC system (9) by 17 μs using the Anton specialized hardware (10) with the CHARMM22/36 protein/lipid force field and CMAP correction (11,12) (see Fig. S1 and Fig. S2 in the Supporting Material). To explore force-field effects, a similar system was simulated for 2 μs using the OPLS force field (13) (see Methods in the Supporting Material). In agreement with x-ray diffraction measurements on melittin in DOPC multilayers (14), melittin partitioned spontaneously into the lipid headgroups at a position below the phosphate groups at similar depth as glycerol/carbonyl groups (Fig. 2).Open in a separate windowFigure 2Melittin partitioned into the polar headgroup region of the lipid bilayer. (A) Snapshot of the simulation cell showing two melittin molecules (MLT1 and MLT2, in yellow) at the lipid-water interface. (B) Density cross-section of the simulation cell extracted from the 17-μs simulation. The peptides are typically located below the lipid phosphate (PO4) groups, in a similar depth as the glycerol/carbonyl (G/C) groups.To describe the secondary structure for each residue, we defined helicity by backbone dihedral angles (φ, ψ) within 30° from the ideal α-helical values (–57°, –47°). The per-residue helicity in the CHARMM simulation displays excellent agreement with amide exchange rates from NMR measurements that show a proline residue to separate two helical segments, which are unfolded below Ala⁵ and above Arg²² (15) (Fig. 3 A). In contrast, the OPLS simulation failed to reproduce the per-residue helicity except for a short central segment (see Fig. S3).Open in a separate windowFigure 3Helicity and conformational distribution of melittin as determined via MD simulation. (A) Helicity per residue for MLT1 and MLT2. (B) Corresponding evolution of the helicity. (C) Conformational distributions over the entire 17-μs simulation.Circular dichroism experiments typically report an average helicity of ∼70% for melittin at membrane interfaces (3,6,16,17), but other methods yield average helicities as high as 85% (15,18). Our CHARMM simulations are generally consistent with the experimental results, especially amide-exchange measurements (15); melittin helicity averaged to 78% for MLT1, whereas MLT2 transitioned from 75% to 89% helicity at t ≈ 8 μs, with an overall average helicity of 82% (Fig. 3 B). However, in the OPLS simulation, melittin steadily unfolds over the first 1.3 μs, after which the peptide remains only partly folded, with an average helicity of 33% (see Fig. S3). Similar force-field-related differences in peptide helicity were recently reported, albeit at shorter timescales (19). Although suitable NMR data are not presently available, we have computed NMR quadrupolar splittings for future reference (see Fig. S4).To answer the question asked in this article—whether the conformational space of folded melittin in the membrane interface can be described by a narrow distribution—the helicity distributions for the equilibrated trajectories are shown in Fig. 3 C. Whereas MLT1 in the CHARMM simulation produces a single, narrow distribution of the helicity, MLT2 has a bimodal distribution as a consequence of the folding event at t ≈ 8 μs (Fig. 3 C). We note that CHARMM force fields have a propensity for helix-formation and this transition might therefore be an artifact. We performed a cluster analysis to describe the structure of the peptide in the membrane interface. The four most populated conformations in the CHARMM simulation are shown in Fig. 4. The dominant conformation for both peptides was a helix kinked at G12 and unfolded at the last 5–6 residues of the C-terminus. The folding transition of MLT2 into a complete helix is visible by the 48% occupancy of a fully folded helix.Open in a separate windowFigure 4Conformational clusters of the two melittin peptides (MLT1 and MLT2) from the 17-μs CHARMM simulation in DOPC. Clustering is based on Cα-RMSD with a cutoff criterion of 2 Å.We conclude that the general assumption when calculating folding energetics holds: Folded melittin partitioned into membrane interfaces can be described by a narrow distribution of conformations. Furthermore, extended (several microsecond) simulations are needed to differentiate force-field effects. Although the CHARMM and OPLS simulations would seem to agree for the first few hundred nanoseconds, the structural conclusions differ drastically with longer trajectories, with CHARMM parameters being more consistent with experiments. However, as implied by the difference in substate distributions between MLT1 and MLT2, 17 μs might not be sufficient to observe the fully equilibrated partitioning process. The abrupt change in MLT2 might indicate that the helicity will increase to greater than experimentally observed in a sufficiently long simulation. On the other hand, it could be nothing more than a transient fluctuation. Increased sampling will provide further indicators of convergence of the helix partitioning process.     

3D-SIM Super-resolution of FtsZ and Its Membrane Tethers in Escherichia coli Cells

FtsZ, a bacterial homolog of eukaryotic tubulin, assembles into the Z ring required for cytokinesis. In Escherichia coli, FtsZ interacts directly with FtsA and ZipA, which tether the Z ring to the membrane. We used three-dimensional structured illumination microscopy to compare the localization patterns of FtsZ, FtsA, and ZipA at high resolution in Escherichia coli cells. We found that FtsZ localizes in patches within a ring structure, similar to the pattern observed in other species, and discovered that FtsA and ZipA mostly colocalize in similar patches. Finally, we observed similar punctate and short polymeric structures of FtsZ distributed throughout the cell after Z rings were disassembled, either as a consequence of normal cytokinesis or upon induction of an endogenous cell division inhibitor.The assembly of the bacterial tubulin FtsZ has been well studied in vitro, but the fine structure of the cytokinetic Z ring it forms in vivo is not well defined. Super-resolution microscopy methods including photoactivated localization microscopy (PALM) and three-dimensional-structured illumination microscopy (3D-SIM) have recently provided a more detailed view of Z-ring structures. Two-dimensional PALM showed that Z rings in Escherichia coli are likely composed of loosely-bundled dynamic protofilaments (1,2). Three-dimensional PALM studies of Caulobacter crescentus initially showed that Z rings were comprised of loosely bundled protofilaments forming a continuous but dynamic ring (1–3). However, a more recent high-throughput study showed that the Z rings of this bacterium are patchy or discontinuous (4), similar to Z rings of Bacillus subtilis and Staphylococcus aureus using 3D-SIM (5). Strauss et al. (5) also demonstrated that the patches in B. subtilis Z rings are highly dynamic.Assembly of the Z ring is modulated by several proteins that interact directly with FtsZ and enhance assembly or disassembly (6). For example, FtsA and ZipA promote ring assembly in E. coli by tethering it to the cytoplasmic membrane (7,8). SulA is an inhibitor of FtsZ assembly, induced only after DNA damage, which sequesters monomers of FtsZ to prevent its assembly into a Z ring (9). Our initial goals were to visualize Z rings in E. coli using 3D-SIM, and then examine whether any FtsZ polymeric structures remain after SulA induction. We also asked whether FtsA and ZipA localized in patchy patterns similar to those of FtsZ.We used a DeltaVision OMX V4 Blaze microscope (Applied Precision, GE Healthcare, Issaquah, WA) to view the high-resolution localization patterns of FtsZ in E. coli cells producing FtsZ-GFP (Fig. 1). Three-dimensional views were reconstructed using softWoRx software (Applied Precision). To rule out GFP artifacts, we also visualized native FtsZ from a wild-type strain (WM1074) by immunofluorescence (IF).Open in a separate windowFigure 1Localization of FtsZ in E. coli. (A) Cell with a Z ring labeled with FtsZ-GFP. (B) Rotated view of Z ring in panel A. (C) Cell with a Z ring labeled with DyLight 550 (Thermo Fisher Scientific, Waltham, MA). (D) Rotated view of Z ring in panel C. (B1 and D1) Three-dimensional surface intensity plots of Z rings in panels B and D, respectively. (E) A dividing cell producing FtsZ-GFP. The cell outline is shown in the schematic. (Asterisk) Focus of FtsZ localization; (open dashed ovals) filamentous structures of FtsZ. Three-dimensional surface intensity plots were created using the software ImageJ (19). Scale bars, 1 μm.Both FtsZ-GFP (Fig. 1, A, B, and B1) and IF staining for FtsZ (Fig. 1, C, D, and D1) consistently localized to patches around the ring circumference, similar to the B. subtilis and C. crescentus FtsZ patterns (4,5). Analysis of fluorescence intensities (see Fig. S1, A and B, in the Supporting Material) revealed that the majority of Z rings contain one or more gaps in which intensity decreases to background levels (82% for FtsZ-GFP and 69% for IF). Most rings had 3–5 areas of lower intensity, but only a small percentage of these areas had fluorescence below background intensity (34% for FtsZ-GFP and 21% for IF), indicating that the majority of areas with lower intensity contain at least some FtsZ.To elucidate how FtsZ transitions from a disassembled ring to a new ring, we imaged a few dividing daughter cells before they were able to form new Z rings (Fig. 1 E). Previous conventional microscopy had revealed dynamic FtsZ helical structures (10), but the resolution had been insufficient to see further details. Here, FtsZ visualized in dividing cells by 3D-SIM localized throughout as a mixture of patches and randomly-oriented short filaments (asterisk and dashed oval in Fig. 1, respectively). These structures may represent oligomeric precursors of Z ring assembly.To visualize FtsZ after Z-ring disassembly another way, we overproduced SulA, a protein that blocks FtsZ assembly. We examined E. coli cells producing FtsZ-GFP after induction of sulA expression from a pBAD33-sulA plasmid (pWM1736) with 0.2% arabinose. After 30 min of sulA induction, Z rings remained intact in most cells (Fig. 2 A and data not shown). The proportion of cellular FtsZ-GFP in the ring before and after induction of sulA was consistent with previous data (data not shown) (1,11).Open in a separate windowFigure 2Localization of FtsZ after overproduction of SulA. (A) Cell producing FtsZ-GFP after 0.2% arabinose induction of SulA for 30 min. (B) After 45 min. (B1) Magnified cell shown in panel B. (C) Cell producing native FtsZ labeled with AlexaFluor 488 (Life Technologies, Carlsbad, CA) 30 min after induction; (D) 45 min after induction. (D1) Magnified cell shown in panel D. Scale bars, 1 μm. (Asterisk) Focus of FtsZ localization; (open dashed ovals) filamentous structures of FtsZ.Notably, after 45 min of sulA induction, Z rings were gone (Fig. 2, B and B1), replaced by numerous patches and randomly-oriented short filaments (asterisk and dashed ovals in Fig. 2), similar to those observed in a dividing cell. FtsZ normally rapidly recycles from free monomers to ring-bound polymers (11), but a critical concentration of SulA reduces the pool of available FtsZ monomers, resulting in breakdown of the Z ring (9). The observed FtsZ-GFP patches and filaments are likely FtsZ polymers that disassemble before they can organize into a ring.We confirmed this result by overproducing SulA in wild-type cells and detecting FtsZ localization by IF (Fig. 2, C, D, and D1). The overall fluorescence patterns in cells producing FtsZ-GFP versus cells producing only native FtsZ were similar (Fig. 2, B1 and D1), although we observed fewer filaments with IF, perhaps because FtsZ-GFP confers slight resistance to SulA, or because the increased amount of FtsZ in FtsZ-GFP producing cells might titrate the SulA more effectively.Additionally, we wanted to observe the localization patterns of the membrane tethers FtsA and ZipA. Inasmuch as both proteins bind to the same C-terminal conserved tail of FtsZ (12–14), they would be expected to colocalize with the circumferential FtsZ patches in the Z ring. We visualized FtsA using protein fusions to mCherry and GFP (data not shown) as well as IF using a wild-type strain (WM1074) (Fig. 3 A). We found that the patchy ring pattern of FtsA localization was similar to the FtsZ pattern. ZipA also displayed a similar patchy localization in WM1074 by IF (Fig. 3 B).Open in a separate windowFigure 3Localization of FtsA (A) and ZipA (B) by IF using AlexaFluor 488. (C) FtsA-GFP ring. (D) Same cell shown in panel C with ZipA labeled with DyLight 550. (C1 and D1) Three-dimensional surface intensity plots of FtsA ring from panel C or ZipA ring from panel D, respectively. (E) Merged image of FtsA (green) and ZipA (red) from the ring shown in panels C and D. (F) Intensity plot of FtsA (green) and ZipA (red) of ring shown in panel E. The plot represents intensity across a line drawn counterclockwise from the top of the ring around the circumference, then into its lumen. Red/green intensity plot and three-dimensional surface intensity plots were created using the software ImageJ (19). Scale bar, 1 μm.To determine whether FtsA and ZipA colocalized to these patches, we used a strain producing FtsA-GFP (WM4679) for IF staining of ZipA using a red secondary antibody. FtsA-GFP (Fig. 3 C) and ZipA (Fig. 3 D) had similar patterns of fluorescence, although the three-dimensional intensity profiles (Fig. 3, C1 and D1) reveal slight differences in intensity that are also visible in a merged image (Fig. 3 E). Quantitation of fluorescence intensities around the circumference of the rings revealed that FtsA and ZipA colocalized almost completely in approximately half of the rings analyzed (Fig. 3 F, and see Fig. S2 A), whereas in the other rings there were significant differences in localization in one or more areas (see Fig. S2 B). FtsA and ZipA bind to the same C-terminal peptide of FtsZ and may compete for binding. Cooperative self-assembly of FtsA or ZipA might result in large-scale differential localization visible by 3D-SIM.In conclusion, our 3D-SIM analysis shows that the patchy localization of FtsZ is conserved in E. coli and suggests that it may be widespread among bacteria. After disassembly of the Z ring either in dividing cells or by excess levels of the cell division inhibitor SulA, FtsZ persisted as patches and short filamentous structures. This is consistent with a highly dynamic population of FtsZ monomers and oligomers outside the ring, originally observed as mobile helices in E. coli by conventional fluorescence microscopy (10) and by photoactivation single-molecule tracking (15). FtsA and ZipA, which bind to the same segment of FtsZ and tether it to the cytoplasmic membrane, usually display a similar localization pattern to FtsZ and each other, although in addition to the differences we detect by 3D-SIM, there are also likely differences that are beyond its ∼100-nm resolution limit in the X,Y plane.As proposed previously (16), gaps between FtsZ patches may be needed to accommodate a switch from a sparse Z ring to a more condensed ring, which would provide force to drive ring constriction (17). If this model is correct, the gaps should close upon ring constriction, although this may be beyond the resolution of 3D-SIM in constricted rings. Another role for patches could be to force molecular crowding of low-abundance septum synthesis proteins such as FtsI, which depend on FtsZ/FtsA/ZipA for their recruitment, into a few mobile supercomplexes.How are FtsZ polymers organized within the Z-ring patches? Recent polarized fluorescence data suggest that FtsZ polymers are oriented both axially and circumferentially within the Z ring in E. coli (18). The seemingly random orientation of the non-ring FtsZ polymeric structures we observe here supports the idea that there is no strong constraint requiring FtsZ oligomers to follow a circumferential path around the cell cylinder. The patches of FtsZ in the unperturbed E. coli Z ring likely represent randomly oriented clusters of FtsZ filaments that are associated with ZipA, FtsA, and essential septum synthesis proteins. New super-resolution microscopy methods should continue to shed light on the in vivo organization of these protein assemblies.     

Mitochondrial Protein Synthesis,Import, and Assembly

The mitochondrion is arguably the most complex organelle in the budding yeast cell cytoplasm. It is essential for viability as well as respiratory growth. Its innermost aqueous compartment, the matrix, is bounded by the highly structured inner membrane, which in turn is bounded by the intermembrane space and the outer membrane. Approximately 1000 proteins are present in these organelles, of which eight major constituents are coded and synthesized in the matrix. The import of mitochondrial proteins synthesized in the cytoplasm, and their direction to the correct soluble compartments, correct membranes, and correct membrane surfaces/topologies, involves multiple pathways and macromolecular machines. The targeting of some, but not all, cytoplasmically synthesized mitochondrial proteins begins with translation of messenger RNAs localized to the organelle. Most proteins then pass through the translocase of the outer membrane to the intermembrane space, where divergent pathways sort them to the outer membrane, inner membrane, and matrix or trap them in the intermembrane space. Roughly 25% of mitochondrial proteins participate in maintenance or expression of the organellar genome at the inner surface of the inner membrane, providing 7 membrane proteins whose synthesis nucleates the assembly of three respiratory complexes.TO think about how mitochondrial proteins are synthesized, imported, and assembled, it is useful to have a clear picture of the organellar structures that they, along with membrane lipids, compose and the functions that they carry out. As almost every schoolchild learns, mitochondria carry out oxidative phosphorylation, the controlled burning of nutrients coupled to ATP synthesis. Since Saccharomyces cerevisiae prefers to ferment sugars, respiration is a dispensable function and nonrespiring mutants are viable [although they cannot undergo meiosis (Jambhekar and Amon 2008)]. However, mitochondria themselves are not dispensable. A substantial fraction of intermediary metabolism occurs in mitochondria (Strathern et al. 1982), and at least one of these pathways, iron–sulfur cluster assembly, is essential for growth (Kispal et al. 2005). Thus, any mutation that prevents the biogenesis of mitochondria by, for example, preventing the import of protein constituents from the cytoplasm, is lethal (Baker and Schatz 1991).The mitochondria of S. cerevisiae are tubular structures at the cell cortex. While the number of distinct compartments can range from 1 to ∼50 depending upon conditions (Stevens 1981; Pon and Schatz 1991), continual fusion and fission events among them effectively form a single dynamic network (Nunnari et al. 1997). The outer membrane surrounds the tubules. The inner membrane has a boundary domain closely juxtaposed beneath the outer membrane and cristae domains that project internally from the boundary into the matrix (Figure 1A). The matrix is the aqueous compartment surrounded by the inner membrane. The aqueous intermembrane space lies between the membranes and is continuous with the space within cristae.Open in a separate windowFigure 1 Overview of mitochondrial structure in yeast. (A) Schematic of compartments comprising mitochondrial tubules. The outer membrane surrounds the organelle. The inner membrane surrounds the matrix and consists of two domains, the inner boundary membrane and the cristae membranes, which are joined at cristae junctions. The intermembrane space lies between the outer membrane and inner membrane. (B) Electron tomograph image of a highly contracted yeast mitochondrion observed en face (a) with the outer membrane (red) and (b) without the outer membrane. Reprinted by permission from John Wiley & Sons from Mannella et al. (2001).Inner membrane cristae are often depicted as baffles emanating from the boundary domain. However, electron tomography of mitochondria from several species, including yeast, shows that cristae actually emanate from the boundary membrane as narrow tubular structures at sites termed “crista junctions” and expand as they project into the matrix (Frey and Mannella 2000; Mannella et al. 2001) (Figure 1B). It seems clear that the boundary and cristae domains of the inner membrane have distinct compositions with respect to the respiratory complexes that are embedded preferentially in the cristae membrane domains, as well as other components (Vogel et al. 2006; Wurm and Jakobs 2006; Rabl et al. 2009; Suppanz et al. 2009; Zick et al. 2009; Davies et al. 2011).The outer and inner boundary membranes are connected at multiple contact sites, at least some of which are involved in protein translocation and may be transient (Pon and Schatz 1991). In addition, there appear to be firm contact sites, not directly involved with protein translocation, preferentially colocalized with crista junctions (Harner et al. 2011a).Overall, there appear to be ∼1000 distinct proteins in yeast mitochondria (Premsler et al. 2009). One series of proteomic studies on highly purified organelles identified 851 proteins thought to represent 85% of the total number of species (Sickmann et al. 2003; Reinders et al. 2006; Zahedi et al. 2006). Another study identified an additional 209 candidates (Prokisch et al. 2004). A computationally driven search for candidates involved in yeast mitochondrial function, coupled with experiments to assay respiratory function and maintenance of mitochondrial DNA (mtDNA), identified 109 novel candidates, although many of these may not be mitochondrial per se (Hess et al. 2009). Taking the boundary and cristae domains together, the inner membrane is the most protein-rich mitochondrial compartment, followed by the matrix (Daum et al. 1982).Only eight of the yeast mitochondrial proteins detected in proteomic studies are encoded by mtDNA and synthesized within the organelle. They are hydrophobic subunits of respiratory complexes III (bc₁ complex or ubiquinol-cytochrome c reductase), IV (cytochrome c oxidase), and V (ATP synthase), as well as a hydrophilic mitochondrial small subunit ribosomal protein. The remaining ∼99% of yeast mitochondrial proteins are encoded by nuclear genes, synthesized in cytoplasmic ribosomes, and imported into the organelle.An overview of known nuclearly encoded mitochondrial protein functions (Figure 2) reveals that ∼25% of them are involved directly in genome maintenance and expression of the eight major mitochondrial genes (Schmidt et al. 2010). The functions of ∼20% of the proteins are not known. Fifteen percent are involved in the well-known processes of energy metabolism. Protein translocation, folding, and turnover functions occupy ∼10% of mitochondrial proteins.Open in a separate windowFigure 2 Classification of identified mitochondrial proteins according to function. Reprinted by permission from Nature Publishing Group from Schmidt et al. (2010).The following discussion reviews our understanding of the biogenesis of mitochondria starting on the outside, the cytoplasm, and working inward through the mitochondrial compartments.     

Novel Growth Regime of MDCK II Model Tissues on Soft Substrates

It is well established that MDCK II cells grow in circular colonies that densify until contact inhibition takes place. Here, we show that this behavior is only typical for colonies developing on hard substrates and report a new growth phase of MDCK II cells on soft gels. At the onset, the new phase is characterized by small, three-dimensional droplets of cells attached to the substrate. When the contact area between the agglomerate and the substrate becomes sufficiently large, a very dense monolayer nucleates in the center of the colony. This monolayer, surrounded by a belt of three-dimensionally packed cells, has a well-defined structure, independent of time and cluster size, as well as a density that is twice the steady-state density found on hard substrates. To release stress in such dense packing, extrusions of viable cells take place several days after seeding. The extruded cells create second-generation clusters, as evidenced by an archipelago of aggregates found in a vicinity of mother colonies, which points to a mechanically regulated migratory behavior.Studying the growth of cell colonies is an important step in the understanding of processes involving coordinated cell behavior such as tissue development, wound healing, and cancer progression. Apart from extremely challenging in vivo studies, artificial tissue models are proven to be very useful in determining the main physical factors that affect the cooperativity of cells, simply because the conditions of growth can be very well controlled. One of the most established cell types in this field of research is the Madin-Darby canine kidney epithelial cell (MDCK), originating from the kidney distal tube (1). A great advantage of this polarized epithelial cell line is that it retained the ability for contact inhibition (2), which makes it a perfect model system for studies of epithelial morphogenesis.Organization of MDCK cells in colonies have been studied in a number of circumstances. For example, it was shown that in three-dimensional soft Matrigel, MDCK cells form a spherical enclosure of a lumen that is enfolded by one layer of polarized cells with an apical membrane exposed to the lumen side (3). These structures can be altered by introducing the hepatocyte growth factor, which induces the formation of linear tubes (4). However, the best-studied regime of growth is performed on two-dimensional surfaces where MDCK II cells form sheets and exhibit contact inhibition. Consequently, the obtained monolayers are well characterized in context of development (5), mechanical properties (6), and obstructed cell migration (7–9).Surprisingly, in the context of mechanics, several studies of monolayer formation showed that different rigidities of polydimethylsiloxane gels (5) and polyacrylamide (PA) gels (9) do not influence the nature of monolayer formation nor the attainable steady-state density. This is supposedly due to long-range forces between cells transmitted by the underlying elastic substrate (9). These results were found to agree well with earlier works on bovine aortic endothelial cells (10) and vascular smooth muscle cells (11), both reporting a lack of sensitivity of monolayers to substrate elasticity. Yet, these results are in stark contrast with single-cell experiments (12–15) that show a clear response of cell morphology, focal adhesions, and cytoskeleton organization to substrate elasticity. Furthermore, sensitivity to the presence of growth factors that are dependent on the elasticity of the substrate in two (16) and three dimensions (4) makes this result even more astonishing. Therefore, we readdress the issue of sensitivity of tissues to the elasticity of the underlying substrate and show that sufficiently soft gels induce a clearly different tissue organization.We plated MDCK II cells on soft PA gels (Young’s modulus E = 0.6 ± 0.2 kPa), harder PA gels (E = 5, 11, 20, 34 kPa), and glass, all coated with Collagen-I. Gels were prepared following the procedure described in Rehfeldt et al. (17); rigidity and homogeneity of the gels was confirmed by bulk and microrheology (see the Supporting Material for comparison). Seeding of MDCK II cells involved a highly concentrated solution dropped in the middle of a hydrated gel or glass sample. For single-cell experiments, cells were dispersed over the entire dish. Samples were periodically fixed up to Day 12, stained for nuclei and actin, and imaged with an epifluorescence microscope. Details are described in the Supporting Material.On hard substrates and glass it was found previously that the area of small clusters expands exponentially until the movement of the edge cannot keep up with the proliferation in the bulk (5). Consequently, the bulk density increases toward the steady state, whereas the density of the edge remains low. At the same time, the colony size grows subexponentially (5). This is what we denote “the classical regime of growth”. Our experiments support these observations for substrates with E ≥ 5 kPa. Specifically, on glass, colonies start as small clusters of very low density of 700 ± 200 cells/mm² (Fig. 1, A and B), typically surrounded by a strong actin cable (Fig. 1, B and C). Interestingly, the spreading area of single cells (Fig. 1 A) on glass was found to be significantly larger, i.e., (2.0 ± 0.9) × 10⁻³ mm². After Day 4 (corresponding cluster area of 600 ± 100 mm²), the density in the center of the colony reached the steady state with 6,800 ± 500 cells/mm², whereas the mean density of the edge profile grew to 4,000 ± 500 cells/mm². This density was retained until Day 12 (cluster area 1800 ± 100 mm²), which is in agreement with previous work (9).Open in a separate windowFigure 1Early phase of cluster growth on hard substrates. (A) Well-spread single cells, and small clusters with a visible actin cable 6 h after seeding. (B) Within one day, clusters densify and merge, making small colonies. (C) Edge of clusters from panel B.In colonies grown on 0.6 kPa gels, however, we encounter a very different growth scenario. The average spreading area of single cells is (0.34 ± 0.3) × 10⁻³ mm², which is six times smaller than on glass substrates (Fig. 2 A). Clusters of only few cells show that cells have a preference for cell-cell contacts (a well-established flat contact zone can be seen at the cell-cell interface in Fig. 2 A) rather than for cell-substrate contacts (contact zone is diffusive and the shape of the cells appears curved). The same conclusion emerges from the fact that dropletlike agglomerates, resting on the substrate, form spontaneously (Fig. 2 A), and that attempts to seed one single cluster of 90,000 cells fail, resulting in a number of three-dimensional colonies (Fig. 2 A). When the contact area with the substrate exceeds 4.7 × 10⁻³ mm², a monolayer appears in the center of such colonies (Fig. 2 B). The colonies can merge, and if individual colonies are small, the collapse into a single domain is associated with the formation of transient irregular structures (Fig. 2 B). Ultimately, large elliptical colonies (average major/minor axis of e = 1.8 ± 0.6) with a smooth edge are formed (Fig. 2 C), unlike on hard substrates where circular clusters (e = 1.06 ± 0.06) with a ragged edge comprise the characteristic phenotype.Open in a separate windowFigure 2Early phase of cluster growth on soft substrates. (A) Twelve hours after seeding, single cells remain mostly round and small. They are found as individual, or within small, three-dimensional structures (top). The latter nucleate a monolayer in their center (bottom), if the contact area with the substrate exceeds ∼5 × 10⁻³ mm². (B) Irregularly-shaped clusters appear due to merging of smaller droplets. A stable monolayer surrounded by a three-dimensional belt of densely packed cells is clearly visible, even in larger structures. (C) All colonies are recorded on Day 4.Irrespective of cluster size, in the new regime of growth, the internal structure is built of two compartments (Fig. 2 B):

• 1.The first is the edge (0.019 ± 0.05-mm wide), a three-dimensional structure of densely packed cells. This belt is a signature of the new regime because on hard substrates the edge is strictly two-dimensional (Fig. 1 C).
• 2.The other is the centrally placed monolayer with a spatially constant density that is very weakly dependent on cluster size and age (Fig. 3). The mean monolayer density is 13,000 ± 2,000 cells/mm², which is an average over 130 clusters that are up to 12 days old and have a size in the range of 10⁻³ to 10 mm², each shown by a data point in Fig. 3. This density is twice the steady-state density of the bulk tissue in the classical regime of growth.Open in a separate windowFigure 3Monolayer densities in colonies grown on 0.6 kPa substrates, as a function of the cluster size and age. Each cluster is represented by a single data point signifying its mean monolayer density. (Black lines) Bulk and (red dashed lines) edge of steady-state densities from monolayers grown on glass substrates. Error bars are omitted for clarity, but are discussed in the Supporting Material.

Until Day 4, the monolayer is very homogeneous, showing a nearly hexagonal arrangement of cells. From Day 4, however, defects start to appear in the form of small holes (typical size of (0.3 ± 0.1) × 10⁻³ mm²). These could be attributed to the extrusions of viable cells, from either the belt or areas of increased local density in the monolayer (inset in Fig. 4). This suggests that extrusions serve to release stress built in the tissue, and, as a consequence, the overall density is decreased.Open in a separate windowFigure 4Cell nuclei within the mother colony and in the neighboring archipelago of second-generation clusters grown on 0.6 kPa gels at Day 12. (Inset; scale bar = 10 μm) Scar in the tissue, a result of a cell-extrusion event. (Main image; scale bar = 100 μm) From the image of cell nuclei (left), it is clear that there are no cells within the scar, whereas the image of actin (right) shows that the cytoplasm of the cells at the edge has closed the hole.Previous reports suggest that isolated MDCK cells undergo anoikis 8 h after losing contact with their neighbors (18). However, in this case, it appears that instead of dying, the extruded cells create new colonies, which can be seen as an archipelago surrounding the mother cluster (Fig. 4). The viability of off-cast cells is further evidenced by the appearance of single cells and second-generation colonies with sizes varying over five orders of magnitude, from Day 4 until the end of the experiment, Day 12. Importantly, no morphological differences were found in the first- and second-generation colonies.In conclusion, we show what we believe to be a novel phase of growth of MDCK model tissue on soft PA gels (E = 0.6 kPa) that, to our knowledge, despite previous similar efforts (9), has not been observed before. This finding is especially interesting in the context of elasticity of real kidneys, for which a Young’s modulus has been found to be between 0.05 and 5 kPa (19,20). This coincides with the elasticity of substrates studied herein, and opens the possibility that the newly found phase of growth has a particular biological relevance. Likewise, the ability to extrude viable cells may point to a new migratory pathway regulated mechanically by the stresses in the tissue, the implication of which we hope to investigate in the future.     

Cell biology in development: Cell division and the maintenance of epithelial order

Epithelia are polarized layers of adherent cells that are the building blocks for organ and appendage structures throughout animals. To preserve tissue architecture and barrier function during both homeostasis and rapid growth, individual epithelial cells divide in a highly constrained manner. Building on decades of research focused on single cells, recent work is probing the mechanisms by which the dynamic process of mitosis is reconciled with the global maintenance of epithelial order during development. These studies reveal how symmetrically dividing cells both exploit and conform to tissue organization to orient their mitotic spindles during division and establish new adhesive junctions during cytokinesis.The association of large numbers of cells in tightly organized epithelial layers is a unique and defining feature of Metazoa. Although classical studies of development once labeled distinct embryonic regions as territories, fields, layers, placodes, and primordia, we now know many of these structures to be primarily constructed from epithelial sheets. Epithelial structure and function are critically dependent on cell polarization, which is coupled to the targeted assembly of adhesive junctions along the apicolateral membranes of adjacent cells (Tepass et al., 2001; Cavey and Lecuit, 2009). In brief, the plasma membrane of epithelial cells is polarized into apical and basolateral domains, each enriched with distinct lipid and protein components (Fig. 1; Rodriguez-Boulan et al., 2005; St Johnston and Ahringer, 2010). At the molecular level, E-cadherins are the major class of adhesion proteins that establish cell–cell connections through homophilic interaction across cell membranes (Takeichi, 1991, 2011; Halbleib and Nelson, 2006; Harris and Tepass, 2010). Whereas E-cadherin is apically enriched in invertebrate epithelia, it is localized along the lateral domain of vertebrate epithelial cells. In both cases, E-cadherin interacts with cytoplasmic actin filaments via the catenin class of adaptor proteins, thus coupling intercellular adhesive contacts to the cytoskeleton (Cavey and Lecuit, 2009; Harris and Tepass, 2010; Gomez et al., 2011). Within this framework, the maintenance of both polarity and cell–cell adhesion are essential for epithelial barrier function and tissue architecture during growth and morphogenesis (Papusheva and Heisenberg, 2010; Guillot and Lecuit, 2013b).Open in a separate windowFigure 1.Architectural implications of orthogonal and planar spindle orientations during epithelial cell division. (A) Programmed orthogonal orientation of the mitotic spindle can promote epithelial stratification, although the remodeling of adhesion and polarity complexes during this process remains an important area for further study. (B) Planar spindle orientation is coordinated with the overall cell polarity machinery and thus facilitates conservation of monolayer organization during rapid cell proliferation.During development, epithelia expand by the combined effects of cell growth (increase in cell size) and cell division (increase in cell numbers). Division events are typically oriented either parallel or orthogonal to the plane of the layer and less frequently at oblique angles (Gillies and Cabernard, 2011). When cells divide orthogonally (perpendicular to the plane of the epithelium), the two daughters will be at least initially nonequivalent with respect to position within the cell layer (Fig. 1 A). Under normal conditions, such programmed orthogonal divisions can be used to effect asymmetric segregation of cell fates or to establish distinct cell types, such as in the developing cortex (Fietz et al., 2010; Hansen et al., 2010) or during morphogenesis of stratified epithelia (Lechler and Fuchs, 2005; Williams et al., 2011). Conversely, when cells divide parallel to the plane of the epithelium (planar orientation; Fig. 1 B), both daughter cells are equivalent with respect to mother cell polarity and tightly integrated in the growing monolayer (Morin and Bellaïche, 2011).During planar division, epithelial cells typically round up, constrict in the middle to form the cytokinetic furrow, and divide symmetrically with respect to the apicobasal axis to produce two equal daughter cells. These daughters construct new cell–cell junctions at their nascent interface, thus integrating into the monolayer (Fig. 2, A–G). Although the intricate relationship between cell polarity and cell division has been explored for many years in the context of asymmetric cell division (Rhyu and Knoblich, 1995; Siller and Doe, 2009; Williams and Fuchs, 2013), recent studies have also begun to explore how epithelia maintain their morphology, integrity, and barrier function during continuous rounds of planar cell division and junction assembly. In this review, we highlight recent findings that provide new insights into the problem of symmetric planar cell division in diverse polarized epithelia, with a focus on two crucial mitotic events: (1) the orientation of cell division and (2) the formation of new cell junctions.Open in a separate windowFigure 2.Progression of planar cell division in an epithelial monolayer. Apical cross section (xy, top row) and longitudinal (xz, bottom row) view of a dividing cell (red). (A) At the level of apical junctions, cells are packed in a polygonal cell arrangement during interphase. (B) In prophase, the dividing nucleus begins to translocate apically as the cell rounds up and maintains a thin basal projection enriched with nonmuscle myosin II and actin (light blue). Notably, this type of nuclear migration is typically observed in pseudostratified columnar epithelia and does not occur in cuboidal and squamous epithelial tissues. (C) Localized molecular landmarks (apical complexes marked as gray bars on cell sides) direct orientation of the mitotic spindle to the plane of the epithelium (arrows). (D) Within the plane of the cell layer, the spindle can be further oriented (arrows) in response to molecular cues, global tissue tension, and local cell geometry. (E and F) After chromosome segregation during anaphase, the cell constricts in the middle and cleaves orthogonal to the plane of the monolayer. (G) After cytokinesis, daughter nuclei move basally and daughter cells form new junctions at their nascent interface (white) while elongating along the apicobasal axis.

Mitotic spindle position and orientation in epithelial cells

Planar orientation of epithelial cell division requires coordinated interaction between the cell polarization machinery and the mitotic spindle itself (Morin and Bellaïche, 2011). In animal cells, the spindle is organized by two symmetrically positioned poles or centrosomes, which nucleate three forms of microtubules (Tanaka, 2010): kinetochore microtubules that attach to the chromosomes, polar microtubules that overlap in an antiparallel fashion over the midplane, and astral microtubules that extend to the cell cortex, which is the actin-rich layer beneath the cell membrane (Lancaster and Baum, 2014). Work in Drosophila melanogaster and vertebrates reveals that at least three factors influence the orientation of this spindle machinery with respect to polarized epithelial architecture: cytoskeletal forces, localized cortical cues, and tissue tension.

Cytoskeletal forces position mitotic nuclei near the apical cell membrane.

In columnar and pseudostratified epithelia where cells elongate along their apicobasal axes, mitotic events are typically restricted to the apical domain of the epithelium (corresponding to the apical membrane of each cell; Fig. 2, C–F). How does the mitotic nucleus achieve the correct apical position? In Drosophila wing discs and zebrafish neuroepithelia, mitotic nuclei and the bulk of the cell cytoplasm are driven apically by actomyosin-dependent cortical contractility at prophase entry (Norden et al., 2009; Leung et al., 2011; Meyer et al., 2011). These events are fundamentally similar to mitotic cell rounding in tissue culture cells (Kunda and Baum, 2009; Lancaster et al., 2013; Lancaster and Baum, 2014). In many epithelia, as the cell rounds up and the nucleus translocates apically, a thin actin-rich projection maintains contact with the basal lamina (Fig. 2, B and C). It remains poorly understood how this structure behaves during cleavage and whether this basal process plays any role in the correct reintegration of the postmitotic daughter cells into the monolayer. Although actomyosin may be the primary driver of apical rounding in many cases, evidence also supports a role for microtubule-based mechanisms in the positioning of premitotic nuclei. In chicken neural tube and mouse cerebral cortex, nuclei migrate apically on microtubules before actomyosin-dependent rounding (Spear and Erickson, 2012a). Centrosomes provide directionality to the microtubules on which the nucleus migrates and organize the spindle once the mitotic chromatin reaches the apical domain (Peyre et al., 2011; Spear and Erickson, 2012a; Nakajima et al., 2013). Collectively, current evidence suggests that both actomyosin- and microtubule-dependent forces conspire to effect mitotic nuclear translocation in a highly context- and species-specific manner. One possibility is that the varying physical dimensions of epithelial cells require varying mechanisms for apical nuclear translocation. For example, highly elongated radial glial cells require active transport of the nucleus on microtubules before mitotic rounding, whereas cortical actomyosin contractility may be sufficient in less elongated cells (Spear and Erickson, 2012b). A major outstanding problem is how cortical contractility triggers cell rounding that is polarized along the apicobasal axis of the cell. Whereas centrosomes function as an apical landmark for nuclei moving on microtubules, it remains unclear what provides the directionality for the basal-to-apical actomyosin contraction. One hypothesis is that certain proteins can restrict the localization of nonmuscle myosin II at the basal domain of epithelial cells. The microcephaly protein Asp interacts with myosin II and regulates its polarized localization along the apicobasal axis in the fly optic lobe neuroepithelium. In asp mutant flies, myosin II is enriched apically instead of basally. Many dividing nuclei fail to reach the apical domain and are thus broadly distributed along the apicobasal axis of the epithelium, leading to a disorganized tissue (Rujano et al., 2013). Interestingly, Asp also interacts with microtubules, associates with spindle poles, and is essential for positioning the spindle in fly and vertebrate epithelia (Saunders et al., 1997; do Carmo Avides and Glover, 1999; Wakefield et al., 2001; Fish et al., 2006). Elucidating the function of proteins such as Asp at the interface of microtubules and actomyosin will be essential to our understanding of how the cytoskeleton drives apical mitotic rounding.

Localized molecular landmarks direct planar spindle orientation.

In most animal cells, the mitotic spindle is anchored to the cell cortex by astral microtubules (Fig. 2, C–E; Théry and Bornens, 2006). Translocation of the dynein–dynactin motor toward the astral microtubule minus ends provides a pulling force on centrosomes and is essential for spindle orientation and pole separation during cell division (Dujardin and Vallee, 2002; Kotak et al., 2012). Molecular cues embedded in the cortex can thus determine spindle orientation by anchoring the dynein–dynactin complex in restricted domains. In cultured MDCK and chick neuroepithelia cells, the Gαi–LGN–nuclear mitotic apparatus (NuMA) complex serves this function (Busson et al., 1998; Hao et al., 2010; Zheng et al., 2010; Peyre et al., 2011). Knockdown or mislocalization of these factors leads to spindle orientation defects that ultimately lead to removal of cell progenitors from the monolayer (Peyre et al., 2011). LGN (Pins in Drosophila) localizes to the lateral cell cortex by binding to the membrane-bound Gαi and enforces spindle orientation by recruiting NuMA (Mud in Drosophila), which binds directly to the dynein–dynactin motor. In certain epithelia, including MDCK cells and Drosophila wing discs, LGN is excluded from the apical domain by atypical PKC (aPKC) phosphorylation, thus restricting it at the lateral cell cortex (Konno et al., 2008; Hao et al., 2010; Zheng et al., 2010; Guilgur et al., 2012). In chick neuroepithelia, however, LGN is restricted at the lateral cortex independently of aPKC, suggesting that other cues control its localization (Morin et al., 2007; Peyre et al., 2011). In the mouse embryonic neocortex, the actin–membrane linkers ERM (ezrin/radixin/moesin) promote the association of LGN with NuMA (Machicoane et al., 2014), indicating that organized cortical actin is critical for correct LGN localization.Cell–cell junctions have been implicated in planar cell division in mammalian epithelia, suggesting a possible direct link between the polarity apparatus and the spindle machinery (Reinsch and Karsenti, 1994; den Elzen et al., 2009). Interfering with E-cadherin function or reducing E-cadherin levels abolishes junctional localization of APC (adenomatous polyposis coli), a microtubule-interacting protein that is required for planar spindle orientation and chromosome alignment (Green et al., 2005; den Elzen et al., 2009). However, spindle orientation may not directly depend on E-cadherin or adherens junctions (AJs) in all cases. In Drosophila follicle cells and imaginal discs as well as Xenopus laevis embryonic epithelia, mitotic spindles exhibit planar orientation but do not align with the AJs (Woolner and Papalopulu, 2012; Bergstralh et al., 2013; Nakajima et al., 2013). Moreover, disruption of AJs in Drosophila follicle cells does not affect spindle position (Bergstralh et al., 2013).In Drosophila wing discs, the spindle poles localize in close proximity to septate junctions, which are positioned immediately basal to AJs (Nakajima et al., 2013). Septate junctions are enriched with many proteins, including the neoplastic tumor suppressors SCRIB (Scribbled) and DLG1 (Discs large 1; Bilder and Perrimon, 2000; Bilder et al., 2000). In asymmetrically dividing cells, such as Drosophila sensory organ precursors and neuroblasts, DLG1 interacts with LGN at the cortex and is required for proper spindle orientation (Bellaïche et al., 2001; Siegrist and Doe, 2005; Johnston et al., 2009). Recent findings indicate that DLG1 is also essential for planar spindle orientation in the symmetric division of epithelial cells. In wing discs, knockdown of scrib or dlg1 leads to randomized spindle orientations. scrib knockdown wing discs exhibit diffuse DLG1 localization but no obvious apicobasal polarity defect, suggesting that epithelial disorganization could be a consequence of aberrant spindle orientation (Nakajima et al., 2013). However, it is not clear whether the septate junctions themselves are important. In Drosophila follicle epithelial cells where septate junctions do not form until relatively late in development (Oshima and Fehon, 2011), DLG1 is localized at the lateral cell cortex and is essential for planar spindle orientation (Bergstralh et al., 2013). Interestingly, dlg1 mutant follicle cells display misoriented divisions yet normal epithelial polarity and tissue organization. In this case, planar spindle orientation appears to be independent of junctions per se but still depends on a DLG1–LGN–NuMA complex, similar to asymmetrically dividing cells (Bergstralh et al., 2013).

Global stress and local cell geometry influence mitotic spindle orientation within the plane of the epithelium.

During planar divisions, the mitotic spindle aligns to the plane of the epithelium (xz; Fig. 2 C) and also within the plane of the cell layer (xy; Fig. 2 D). Studies in gastrulating zebrafish embryos revealed a role for the Wnt–Frizzled–planar cell polarity signaling pathway in orienting cell divisions (Concha and Adams, 1998; Gong et al., 2004). Similarly, the atypical cadherins Fat and Dachsous are involved in orienting cell divisions in the Drosophila wing and in developing mouse kidneys (Baena-López et al., 2005; Saburi et al., 2008). Although both of these pathways have been reviewed elsewhere (Morin and Bellaïche, 2011), recent studies also point to at least two other mechanisms that may independently influence spindle orientation within the plane of the monolayer: (1) global tissue stress and (2) local epithelial cell geometry.Epithelial cell shape and spindle orientation are modulated by global stress that accumulates during tissue growth. In Drosophila wing discs, cells in the center of the wing blade primordium proliferate at a faster rate than in the periphery. Consequently, cells in the periphery are mechanically stretched, and cells in the center are compressed. As a result of stretching, peripheral cells localize myosin II at their cortex and align their mitotic spindle with the stretch axis (LeGoff et al., 2013; Mao et al., 2013). Similarly, epithelial cells of the enveloping cell layer in gastrulating zebrafish embryos elongate and orient their spindle along the direction of tension generated by spreading during epiboly (Campinho et al., 2013). It is unclear whether myosin II directly conveys cell tension to the mitotic apparatus, and it will be necessary to dissect whether cell elongation alone or additional mechanosensing pathways signal cell tension to the mitotic spindle. Keratinocytes from the mammalian epidermis reorient their mitotic spindle in response to mechanical stretch in a NuMA-dependent manner. The mitotic spindle aligns with the cortical NuMA-localized crescent upon stretch and fails to orient when NuMA levels are reduced (Seldin et al., 2013). In summary, global tension generated by growth and cell spreading impact division orientation, suggesting that shape changes in proximity to dividing cells may also lead to a similar effect.Although variations certainly exist, the apical surfaces of proliferating epithelia tend to feature a consistent percentage of hexagonal, pentagonal, heptagonal, and octagonal cell shapes (Gibson et al., 2006; Farhadifar et al., 2007; Aegerter-Wilmsen et al., 2010). In Drosophila imaginal discs, these local patterns of cell packing may systematically influence spindle orientation, as mitotic cells are biased toward cleaving their common interfaces with subhexagonal neighbors (less than six sides) and avoid cleaving their interfaces with superhexagonal neighbors (more than six sides; Gibson et al., 2011). Although the mechanisms underlying the effect of local cell geometry remain elusive, cell packing influences mitotic cell shape and the distribution of adhesive cues, both of which could, in turn, bias spindle orientation. Indeed, dividing cells maintain contacts with their neighbors, which can influence the cell cortex and direct spindle orientation (Goldstein, 1995; Wang et al., 1997). The distribution of adhesions between epithelial cells may also alter the position or action of cortical force generators that interact with spindle microtubules in the mitotic cell. In support of this idea, when single cells are placed on micropatterned substrates, they orient their spindle relative to the geometry of their adhesion pattern and not their cell shape (Théry et al., 2005, 2007). Alternatively, neighbors of different polygonal shapes could stretch the mitotic cell, thus imposing a bias on its long axis. Indeed, sea urchin embryos orient their spindles to divide their longest axis (Hertwig, 1884) and can even sense complex cell geometries to orient their spindles accordingly (Minc et al., 2011). Still, precisely how the interphase morphology of epithelial cells might impinge on mitotic spindle orientation remains an open question.

Genesis of nascent junctions during epithelial cell division

After spindle orientation, the essential processes of cytokinesis and abscission are driven by the assembly and contraction of an actomyosin ring positioned in the cleavage plane (Fededa and Gerlich, 2012). In epithelia, ring contraction accompanied by membrane invagination ultimately gives rise to a new junctional interface between nascent daughter cells. Precisely how this new interface forms remains poorly understood. Recent studies in Drosophila epithelia reveal that, during cytokinesis, (a) E-cadherin levels are reduced at the interface between the cleavage furrow of dividing cells and their neighbors (Fig. 3), and (b) neighbor tension and midbody position guide establishment of new AJs in context with local epithelial geometry (Fig. 4).Open in a separate windowFigure 3.Cytokinetic membrane dynamics in epithelial cells. (A) Cytokinesis of a dividing epithelial cell (yellow) presents several unique structural considerations not addressed by the analysis of single cells. Recent studies (Founounou et al., 2013; Guillot and Lecuit, 2013a; Herszterg et al., 2013, 2014) report a local reduction of E-cadherin levels in proximity to the contractile ring in the dividing cell and its neighbor (red). How cytokinesis is resolved from there may vary in a context-dependent manner. (B) In Drosophila embryos, ring contraction leads to E-cadherin disengagement, and a gap forms between the mitotic cell and its neighbor (Guillot and Lecuit, 2013a). (C) In the Drosophila pupal notum, the contractile ring pulls the neighbor cell plasma membrane into the cleavage furrow, perhaps enabled by uncoupling of the membrane and the cortex in the neighbor (Herszterg et al., 2013, 2014).Open in a separate windowFigure 4.New AJ formation in dividing epithelial cells. Apical cross section (xy, top row) and longitudinal (xz, bottom row) view of a dividing epithelial cell (red). (A) Opposing forces (black vertical arrows) develop between the contractile ring in the dividing cell and the two neighboring cells (orange) in proximity to the cleavage furrow. E-cadherin clusters are reduced at the furrow/neighbor interface. (B) Myosin II and tension build up in the neighboring cell, causing the nascent daughter cells to juxtapose their plasma membranes at the presumptive site of junction assembly (black horizontal arrows). (C) Arp2/3 and Rac1 drive actin polymerization at the daughter cell interface around the midbody, stabilizing the nascent junction as the neighboring cell membrane withdraws. (D) The new junction is complete and of suitable length in context with the local epithelial geometry.

Mitotic cells remodel their adhesion junctions during cytokinesis.

Two kinds of forces are at work during cytokinesis: an active force in the dividing cell caused by ring contraction and a reactive force in contacting neighbors caused by their resistance to pulling to maintain their shape (Fig. 3 and Fig. 4 A; Founounou et al., 2013; Guillot and Lecuit, 2013a; Herszterg et al., 2013). Recent results indicate that these opposing forces can lead to a transient and partial reduction of cell adhesion after mitotic exit. In Drosophila epithelia, E-cadherin levels are reduced at the interface between the cleavage furrow of the dividing cell and its neighbors (Fig. 3, B and C; Founounou et al., 2013; Guillot and Lecuit, 2013a; Herszterg et al., 2013; Morais-de-Sá and Sunkel, 2013). Specifically in embryonic epithelia, the local reduction of E-cadherin facilitates membrane separation, and a gap appears between the dividing cell and its neighbors (Fig. 3 B; Guillot and Lecuit, 2013a). In the dorsal thorax, in contrast, the neighbor cell plasma membrane detaches from the cortex and is drawn into the cleavage furrow (Fig. 3 C; Herszterg et al., 2013, 2014). What triggers E-cadherin modulation in cells after mitotic exit? The loss of overall cell polarity is one possible mechanism. During mitosis in Drosophila, follicular epithelial cells lose cortical enrichment of some apical polarity proteins (aPKC, Crumbs, and Bazooka/Par3; Bergstralh et al., 2013; Morais-de-Sá and Sunkel, 2013), and embryonic cells lose localization of lethal giant larvae, a basolateral cortical protein (Huang et al., 2009). Contrasting with these observations, however, MDCK cells and Drosophila embryonic and dorsal thorax epithelial cells appear to maintain apicobasal polarity as they divide (Reinsch and Karsenti, 1994; Founounou et al., 2013; Guillot and Lecuit, 2013a; Herszterg et al., 2013). Furthermore, E-cadherin reduction is limited to the furrow/neighbor interface and is not observed in other areas of cell contact. Therefore, an alternative mechanism that explains local E-cadherin modulation is mechanical tension that arises precisely at the area between the contractile ring and the neighboring cell membrane (Founounou et al., 2013; Guillot and Lecuit, 2013a).Does E-cadherin modulation serve a functional role in mitotic cells? In Drosophila embryonic and dorsal thorax epithelia, E-cadherin decrease leads to a local adhesion disengagement proposed to facilitate the formation of new AJs between daughter cells (Founounou et al., 2013; Guillot and Lecuit, 2013a). It has been previously reported that cells maintain their AJs throughout division. For example, intercellular junctions are maintained in dividing cells of human colonic mucosal crypt cells and basal keratinocytes (Baker and Garrod, 1993). Similarly, mitotic MDCK cultured cells maintain tight junctions apically and E-cadherin basolaterally (Reinsch and Karsenti, 1994). The E-cadherin loss in certain Drosophila epithelia may be either a tissue-specific phenomenon or a highly dynamic process only observable with the temporal resolution of live-cell imaging. Moreover, dividing cells in the Drosophila dorsal thorax show decreased levels of E-cadherin yet maintain their cohesiveness (Herszterg et al., 2013). Interestingly, E-cadherin is internalized in mitotic MDCK cells (Bauer et al., 1998). It will therefore be important to investigate whether loss of E-cadherin leads to adhesion disengagement in other epithelial tissues and whether tension alone or in combination with biochemical pathways is responsible for E-cadherin modulation.

Epithelial neighbors exert tension on daughter cell membranes to facilitate new AJ formation.

How new junctional contacts form during mitosis is a poorly understood problem at the heart of epithelial cell biology. In Drosophila, new membrane interfaces between nascent daughter cells initially show only a weak level of E-cadherin clusters (Guillot and Lecuit, 2013a; Herszterg et al., 2013). Subsequently, the daughter cells assemble their AJs de novo. How is the length of these new junctions determined with respect to cell geometry? Recent evidence indicates that AJ length is a function of local cell packing within the epithelium. In dividing cells of the Drosophila dorsal thorax, the contractile ring triggers tension and accumulation of myosin II in neighbors at the furrow/neighbor interface (Fig. 4 B; Founounou et al., 2013; Herszterg et al., 2013). Myosin II in the neighboring cells in turn contracts and creates tension at the furrowing membrane of the nascent daughter cells, keeping them tightly pressed against each other (Fig. 4, B and C). This local membrane juxtaposition facilitates AJ formation. To allow expansion of the daughter cell interface and maintain AJ length, branched actin polymerization via Rac1 and Arp2/3 is oriented to the midbody, which serves as a positional landmark for new AJs (Fig. 4 C; Herszterg et al., 2013). The midbody is a narrow intercellular bridge that remains after the contracted cytokinetic ring has driven membrane invagination, and it recruits the abscission factors that will eventually separate the daughter cells (Fededa and Gerlich, 2012). Interestingly, the midbody is positioned apically as a result of the presence of AJs. In Drosophila follicular epithelia, the midbody also provides cues for the formation of the apical daughter cell interface, suggesting that it plays a role in both AJ and epithelial cell polarity establishment and maintenance in dividing epithelial cells (Morais-de-Sá and Sunkel, 2013). Thus, examples from Drosophila epithelia show that cohesion between dividing cells and their neighbors together with the apically positioned midbody provides a spatial template and polarized positional cue for de novo AJ assembly (Herszterg et al., 2013; Morais-de-Sá and Sunkel, 2013). Further work on other epithelial tissues may provide alternative mechanisms of junction biogenesis.

Growth and order in the epithelium: Thinking outside the cell

During development, epithelial monolayers have the remarkable capacity to maintain specialized morphologies and barrier functions during rapid cell proliferation. Mitotic cells remain adherent to their neighbors throughout cell division. Cell cohesion enables local geometry and global tissue tension to instruct mitotic cells where to position their cleavage plane and how to assemble their junctions. However, local tension may also lead to a transient disengagement of dividing cells from their neighbors after mitotic exit. How is global and local tension conveyed to protein complexes in mitotic cells so that different outcomes take place? Moreover, it is unclear whether and how tissue tension instructs synchronously dividing epithelial cells how to divide and reestablish their junctions after division. Clearly, this is a fundamental problem for the maintenance of epithelial order and may be linked to the origin of epithelial cancers, in which cells undergo rapid proliferation but fail to remain integrated into the monolayer.The selected studies discussed here hint at the remarkable level of coordination that occurs during epithelial cell division, recasting mitosis as a truly multicellular process. Looking ahead, understanding the interface between cells, proteins, and mechanical forces that each operate on different scales will require creative multidisciplinary approaches in diverse organismal systems. Indeed, epithelial organization is widespread in nature and is encountered among even the most basal animals, including sponges and cnidarians as well as the fruiting body of the nonmetazoan social amoeba Dictyostelium discoideum (Wood, 1959; Ereskovsky et al., 2009; Houliston et al., 2010; Dickinson et al., 2011; Meyer et al., 2011). Combined, future interdisciplinary studies and a fresh look at diverse animal models should yield new insight into epithelial cell division for many years to come.     

AUXIN RESPONSE FACTOR17 Is Essential for Pollen Wall Pattern Formation in Arabidopsis

In angiosperms, pollen wall pattern formation is determined by primexine deposition on the microspores. Here, we show that AUXIN RESPONSE FACTOR17 (ARF17) is essential for primexine formation and pollen development in Arabidopsis (Arabidopsis thaliana). The arf17 mutant exhibited a male-sterile phenotype with normal vegetative growth. ARF17 was expressed in microsporocytes and microgametophytes from meiosis to the bicellular microspore stage. Transmission electron microscopy analysis showed that primexine was absent in the arf17 mutant, which leads to pollen wall-patterning defects and pollen degradation. Callose deposition was also significantly reduced in the arf17 mutant, and the expression of CALLOSE SYNTHASE5 (CalS5), the major gene for callose biosynthesis, was approximately 10% that of the wild type. Chromatin immunoprecipitation and electrophoretic mobility shift assays showed that ARF17 can directly bind to the CalS5 promoter. As indicated by the expression of DR5-driven green fluorescent protein, which is an synthetic auxin response reporter, auxin signaling appeared to be specifically impaired in arf17 anthers. Taken together, our results suggest that ARF17 is essential for pollen wall patterning in Arabidopsis by modulating primexine formation at least partially through direct regulation of CalS5 gene expression.In angiosperms, the pollen wall is the most complex plant cell wall. It consists of the inner wall, the intine, and the outer wall, the exine. The exine is further divided into sexine and nexine layers. The sculptured sexine includes three major parts: baculum, tectum, and tryphine (Heslop-Harrison, 1971; Piffanelli et al., 1998; Ariizumi and Toriyama, 2011; Fig. 1A). Production of a functional pollen wall requires the precise spatial and temporal cooperation of gametophytic and sporophytic tissues and metabolic events (Blackmore et al., 2007). The intine layer is controlled gametophytically, while the exine is regulated sporophytically. The sporophytic tapetum cells provide material for pollen wall formation, while primexine determines pollen wall patterning (Heslop-Harrison, 1968).Open in a separate windowFigure 1.Schematic representation of the pollen wall and primexine development. A, The innermost layer adjacent to the plasma membrane is the intine. The bacula (Ba), tectum (Te), and tryphine (T) make up the sexine layer. The nexine is located between the intine and the sexine layers. The exine includes the nexine and sexine layers. B, Primexine (Pr) appears between callose (Cl) and plasma membrane (Pm) at the early tetrad stage (left panel). Subsequently, the plasma membrane becomes undulated (middle panel) and sporopollenin deposits on the peak of the undulated plasma membrane to form bacula and tectum (right panel).After meiosis, four microspores were encased in callose to form a tetrad. Subsequently, the primexine develops between the callose layer and the microspore membrane (Fig. 1B), and the microspore plasma membrane becomes undulated (Fig. 1B; Fitzgerald and Knox, 1995; Southworth and Jernstedt, 1995). Sporopollenin precursors then accumulate on the peak of the undulated microspore membrane to form the bacula and tectum (Fig. 1B; Fitzgerald and Knox, 1995). After callose degradation, individual microspores are released from the tetrad, and the bacula and tectum continue to grow into exine with further sporopollenin deposition (Fitzgerald and Knox, 1995; Blackmore et al., 2007).The callose has been reported to affect primexine deposition and pollen wall pattern formation. The peripheral callose layer, secreted by the microsporocyte, acts as the mold for primexine (Waterkeyn and Bienfait, 1970; Heslop-Harrison, 1971). CALLOSE SYNTHASE5 (CalS5) is the major enzyme responsible for the biosynthesis of the callose peripheral of the tetrad (Dong et al., 2005; Nishikawa et al., 2005). Mutation of Cals5 and abnormal CalS5 pre-mRNA splicing resulted in defective peripheral callose deposition and primexine formation (Dong et al., 2005; Nishikawa et al., 2005; Huang et al., 2013). Besides CalS5, four membrane-associated proteins have also been reported to be involved in primexine formation: DEFECTIVE EXINE FORMATION1 (DEX1; Paxson-Sowders et al., 1997, 2001), NO EXINE FORMATION1 (NEF1; Ariizumi et al., 2004), RUPTURED POLLEN GRAIN1 (RPG1; Guan et al., 2008; Sun et al., 2013), and NO PRIMEXINE AND PLASMA MEMBRANE UNDULATION (NPU; Chang et al., 2012). Mutation of DEX1 results in delayed primexine formation (Paxson-Sowders et al., 2001). The primexine in nef1 is coarse compared with the wild type (Ariizumi et al., 2004). The loss-of-function rpg1 shows reduced primexine deposition (Guan et al., 2008; Sun et al., 2013), while the npu mutant does not deposit any primexine (Chang et al., 2012). Recently, it was reported that Arabidopsis (Arabidopsis thaliana) CYCLIN-DEPENDENT KINASE G1 (CDKG1) associates with the spliceosome to regulate the CalS5 pre-mRNA splicing for pollen wall formation (Huang et al., 2013). Clearly, disrupted primexine deposition leads to aberrant pollen wall patterning and ruptured pollen grains in these mutants.The plant hormone auxin has multiple roles in plant reproductive development (Aloni et al., 2006; Sundberg and Østergaard, 2009). Knocking out the two auxin biosynthesis genes, YUC2 and YUC6, caused an essentially sterile phenotype in Arabidopsis (Cheng et al., 2006). Auxin transport is essential for anther development; defects in auxin flow in anther filaments resulted in abnormal pollen mitosis and pollen development (Feng et al., 2006). Ding et al. (2012) showed that the endoplasmic reticulum-localized auxin transporter PIN8 regulates auxin homeostasis and male gametophyte development in Arabidopsis. Evidence for the localization, biosynthesis, and transport of auxin indicates that auxin regulates anther dehiscence, pollen maturation, and filament elongation during late anther development (Cecchetti et al., 2004, 2008). The role of auxin in pollen wall development has not been reported.The auxin signaling pathway requires the auxin response factor (ARF) family proteins (Quint and Gray, 2006; Guilfoyle and Hagen, 2007; Mockaitis and Estelle, 2008; Vanneste and Friml, 2009). ARF proteins can either activate or repress the expression of target genes by directly binding to auxin response elements (AuxRE; TGTCTC/GAGACA) in the promoters (Ulmasov et al., 1999; Tiwari et al., 2003). The Arabidopsis ARF family contains 23 members. A subgroup in the ARF family, ARF10, ARF16, and ARF17, are targets of miRNA160 (Okushima et al., 2005b; Wang et al., 2005). Plants expressing miR160-resistant ARF17 exhibited pleiotropic developmental defects, including abnormal stamen structure and reduced fertility (Mallory et al., 2005). This indicates a potential role for ARF17 in plant fertility, although the detailed function remains unknown. In addition, ARF17 was also proposed to negatively regulate adventitious root formation (Sorin et al., 2005; Gutierrez et al., 2009), although an ARF17 knockout mutant was not reported and its phenotype is unknown.In this work, we isolated and characterized a loss-of-function mutant of ARF17. Results from cytological observations suggest that ARF17 controls callose biosynthesis and primexine deposition. Consistent with this, the ARF17 protein is highly abundant in microsporocytes and tetrads. Furthermore, we demonstrate that the ARF17 protein is able to bind the promoter region of CalS5. Our results suggest that ARF17 regulates pollen wall pattern formation in Arabidopsis.     

Interplay between Velocity and Travel Distance of Kinesin-based Transport in the Presence of Tau

Although the disease-relevant microtubule-associated protein tau is known to severely inhibit kinesin-based transport in vitro, the potential mechanisms for reversing this detrimental effect to maintain healthy transport in cells remain unknown. Here we report the unambiguous upregulation of multiple-kinesin travel distance despite the presence of tau, via decreased single-kinesin velocity. Interestingly, the presence of tau also modestly reduced cargo velocity in multiple-kinesin transport, and our stochastic simulations indicate that the tau-mediated reduction in single-kinesin travel underlies this observation. Taken together, our observations highlight a nontrivial interplay between velocity and travel distance for kinesin transport, and suggest that single-kinesin velocity is a promising experimental handle for tuning the effect of tau on multiple-kinesin travel distance.Conventional kinesin is a major microtubule-based molecular motor that enables long-range transport in living cells. Although traditionally investigated in the context of single-motor experiments, two or more kinesin motors are often linked together to transport the same cargo in vivo (1–4). Understanding the control and regulation of the group function of multiple kinesins has important implications for reversing failure modes of transport in a variety of human diseases, particularly neurodegenerative diseases. Tau is a disease-relevant protein enriched in neurons (5,6). The decoration of microtubules with tau is known to strongly inhibit kinesin transport in vitro (7–9), but how kinesin-based transport is maintained in the presence of high levels of tau, particularly in healthy neurons, remains an important open question. To date, no mechanism has been directly demonstrated to reverse the inhibitory effect of tau on kinesin-based transport. Here we present a simple in vitro study that demonstrates the significant upregulation of multiple-kinesin travel distance with decreasing ATP concentration, despite the presence of tau.This investigation was motivated by our recent finding that single-kinesin velocity is a key controller for multiple-kinesin travel distance along bare microtubules (10). The active stepping of each kinesin motor is stimulated by ATP (11), and each kinesin motor remains strongly bound to the microtubule between successive steps (10,11). As demonstrated for bare microtubules (10), with decreasing ATP concentrations, each microtubule-bound kinesin experiences a decreased stepping rate per unit time and spends an increased fraction of time in the strongly bound state; additional unbound kinesins on the same cargo have more time to bind to the microtubule before cargo travel terminates. Thus, reductions in single-kinesin velocity increase the probability that at least one kinesin motor will remain bound to the microtubule per unit time, thereby increasing the travel distance of each cargo (10). Because this effect only pertains to the stepping rate of each individual kinesin and does not address the potential presence of roadblocks such as tau on the microtubules, we hypothesized in this study that single-kinesin velocity may be exploited to relieve the impact of tau on multiple-kinesin travel distance.We focused our in vitro investigation on human tau 23 (htau23, or 3RS tau), an isoform of tau that exhibits the strongest inhibitory effect on kinesin-based transport (7–9). Importantly, htau23 does not alter the stepping rate of individual kinesins (7,9), supporting our hypothesis and enabling us to decouple single-kinesin velocity from the potential effects of tau. We carried out multiple-kinesin motility experiments using polystyrene beads as in vitro cargos (8,10), ATP concentration as an in vitro handle to controllably tune single-kinesin velocity (10,11), and three input kinesin concentrations to test the generality of potential findings for multiple-kinesin transport. Combined with previous two-kinesin studies (10,12), our measurements of travel distance (Fig. 1 A) indicate that the lowest kinesin concentration employed (0.8 nM) corresponds to an average of ∼2–3 kinesins per cargo. Note that in the absence of tau, the observed decrease in bead velocity at the higher kinesin concentrations (Fig. 1 A) is consistent with a recent in vitro finding (13). At 1 mM ATP, htau23 reduced kinesin-based travel distance by a factor of two or more (Fig. 1, A and B). This observation is in good agreement with previous reports (7,8).Open in a separate windowFigure 1Distributions of multiple-kinesin travel distances measured at three experimental conditions, to verify the effect of tau (A and B) and to investigate the impact of single-kinesin velocity on the tau effect (B and C). Shaded bars at 8.7 μm indicate counts of travel exceeding the field of view. The mean travel distance (d; ± standard error of mean, SEM), sample size (n), and corresponding mean velocity (v; ± SEM) are indicated. MT denotes microtubule. Mean travel distance increased substantially at 20 μM ATP (C), despite the presence of htau23. This effect persisted across all three kinesin concentrations tested (left to right).Consistent with our hypothesis, reducing the available ATP concentration to 20 μM increased the multiple-kinesin travel distance by >1.4-fold for all three input kinesin concentrations (Fig. 1, B and C), despite the presence of htau23. The corresponding reduction in single-kinesin velocity with decreasing ATP concentration (10,11) is reflected in the ∼3.4-fold reduction in the measured bead velocities (Fig. 1, B and C). Therefore, the strong negative relationship between single-kinesin velocity and multiple-kinesin travel distance occurs not only for bare microtubules (10), but also for tau-decorated microtubules.What causes the observed increase in travel distance at the lower ATP concentration (Fig. 1, B and C)? In addition to the mechanism discussed above for the case of bare microtubules (10), an intriguing mechanism was suggested by recent studies of tau-microtubule interactions in which htau23 was observed to dynamically diffuse along microtubule lattices (14,15): reducing the stepping rate of a microtubule-bound kinesin may effectively increase the probability that a tau roadblock can diffuse away before the kinesin takes its next step.Perhaps surprisingly, although htau23 does not impact single-kinesin velocity (7,9), we observed a modest reduction in the average velocity of multiple-kinesin transport in experiments using tau-decorated microtubules (Fig. 1, A and B). This decreased velocity reflects a substantially larger variance in the instantaneous velocity for bead trajectories in the presence of htau23 (see Fig. S1 in the Supporting Material), as quantified by parsing each bead trajectory into a series of constant-velocity segments using a previously developed automatic software incorporating Bayesian statistics (16).To test the possibility that single-kinesin travel distance impacts multiple-kinesin velocity, we performed stochastic simulations (see the Supporting Material) that assumed N identical kinesin motors available for transport and included kinesin’s detachment kinetics (17). Previously, this model successfully captured multiple-dynein travel distances in vivo using single-dynein characteristics measured in vitro (18). In this study, we introduced one (and only one) free parameter to reflect the probability of each bound kinesin encountering tau at each step. When encountering tau, each kinesin has a 54% probability of detaching from the microtubule (interpolated from Fig. 2A of Dixit et al. (7)); the undetached kinesin is assumed to remain engaged in transport and completes its step along the microtubule despite the presence of tau.Remarkably, our simple simulation suggested that the tau-mediated reduction in single-kinesin travel is sufficient to reduce multiple-kinesin velocity (Fig. 2 A). The majority of the velocity decrease is predicted to occur at the transition from single-kinesin to two-kinesin transport (Fig. 2). Further decreases in cargo velocity with increasing motor number are predicted to be modest and largely independent of tau (Fig. 2 B). The results of our simulation remain qualitatively the same when evaluated at two bounds (40 and 65%) encompassing the interpolated 54% probability of kinesin detaching at tau (see Fig. S2).Open in a separate windowFigure 2Stochastic simulations predict a tau-dependent reduction in multiple-kinesin velocity, assuming that the only effect of tau protein is to prematurely detach kinesin from the microtubule (or, to reduce single-kinesin travel distance). (A) Average velocity of cargos carried by the indicated number of kinesins was evaluated at 1 mM ATP, and for four probabilities that a kinesin may encounter tau at each step. Mean velocity was evaluated using 600 simulated trajectories for all simulation conditions. Error bars indicate SEM. (B) Change in cargo velocity with each additional kinesin (ΔVel/kinesin) as a function of tau-encounter probability. These values were calculated from cargo velocities shown in panel A. Error bars indicate SEM.We note that our simple simulations do not consider the possibility that kinesin may pause in front of a tau roadblock, as previously reported in Dixit et al. (7). We omitted this consideration because the interaction strength between kinesin and the microtubule in such a paused state is unknown. In a multiple-motor geometry, could a paused kinesin be dragged along by the other motors bound to the same cargo? Could a tau roadblock be forcefully swept off the microtubule surface by the collective motion of the cargo-motor complex? Significant experimental innovations are necessary to specifically address these questions in future multiple-motor assays and to guide modeling efforts. Nonetheless, our simple simulation demonstrates that reducing single-kinesin travel distance is sufficient to decrease multiple-kinesin travel distance.Taken together, our observations highlight a nontrivial interplay between velocity and travel distance for kinesin-based transport in the presence of tau. We uncover a previously unexplored dual inhibition of tau on kinesin-transport: in addition to limiting cargo travel distance, the tau-mediated reduction in single-kinesin travel distance also leads to a modest reduction in multiple-kinesin velocity. We provide what we believe to be the first demonstration of the unambiguous upregulation of multiple-kinesin travel distance despite the presence of tau, via reducing single-kinesin velocity, suggesting a mechanism that could be harnessed for future therapeutic interventions in diseases that result from aberrant kinesin-based transport.     

Voltage Sensor Gating Charge Transfer in a hERG Potassium Channel?Model

Relaxation of a hERG K⁺ channel model during molecular-dynamics simulation in a hydrated POPC bilayer was accompanied by transitions of an arginine gating charge across a charge transfer center in two voltage sensor domains. Inspection of the passage of arginine side chains across the charge transfer center suggests that the unique hydration properties of the arginine guanidine cation facilitates charge transfer during voltage sensor responses to changes in membrane potential, and underlies the preference of Arg over Lys as a mobile charge carrier in voltage-sensitive ion channels.The response of voltage-sensitive ion channels to changes in membrane potential is mediated by voltage sensor domains (VSD) containing a transmembrane helical segment (S4) with a repeating motif of positively charged and hydrophobic amino acids (Fig. 1) (1,2). Changes in membrane potential drive the S4 helix through the membrane plane with the charged side chains (largely arginine) on S4 swapping Glu/Asp carboxylate partners that lie on less mobile elements of the VSD (2). Movement of S4 is coupled to the ion-conducting pore to transmit changes in membrane potential to channel gating (3).Open in a separate windowFigure 1Structures of the VSD of membrane domains before MD in a POPC bilayer. The S2 (pink) and S4 (blue) helices of the VSD of the hERG model (A) and Kv1.2/2.1 chimera structure (B) are highlighted. (C) Sequence alignment of S2 and S4 among homologous voltage-sensitive K⁺ channels.The VSD charge-pairing motif of K⁺ and Na⁺ channels is best represented in VSD states at zero membrane potential (S4 helix up) for which crystal structures exist for Kv1.2 (4), Kv1.2/2.1 chimera (5), and Na_v channels (6,7). In these states, positively charged residues on the intra- and extracellular sections of the S4 helix are separated by a hydrophobic charge-transfer center (CTC) (1) or plug (8) containing a highly conserved Phe residue (Fig. 1). This plug restricts water incursion across the VSD, focusing the electric field across a narrow region near the bilayer center. In voltage-driven transitions between S4 down- and up-states, positively charged S4 side chains move across the CTC.The ether-à-go-go (eag) and eag-related family of voltage-sensitive K⁺ channels likely share similar charge pairing interactions with VSDs in other channels (9,10). However, eag VSDs contain an extra negative charge on S2 (underlined in Fig. 1 C) so that in hERG, Asp residues (D460 and D466) lie approximately one helical turn above and below the conserved charge-transfer center Phe (F463) (Fig. 1). This eag-specific motif might be expected to facilitate transfer of Arg side chains through the CTC and to stabilize the voltage sensor (VS) in the up state. We recently described an open state (VS-up) hERG model built on the crystal structure template of the Kv1.2/2.1 chimera and molecular-dynamics (MD) simulation of this model in a hydrated POPC bilayer (11). We have inspected an extended version of this simulation and identified transitions of a gating charge into the CTC despite the absence of a membrane potential change. These transitions are absent in equivalent MD simulations of the chimera structure in a POPC bilayer.Fig. 1 shows a single VS from starting structures of the hERG model and the chimera structure in a hydrated POPC bilayer, after restrained MD to anneal the protein-lipid interface (see Methods in the Supporting Material). Because the hERG model is constructed on the chimera structure according to the alignment in Fig. 1 the pattern of pairing between S4 charges and acidic VS side chains is equivalent in the hERG model and chimera structure.The arrangement of charge-paired side chains remains constant during MD in all subunits of the chimera (e.g., Fig. 2 E and see Fig. S2 in the Supporting Material). However, in two subunits of the hERG model the R534 side chain moves toward the extracellular side of the bilayer, sliding into the CTC to form a charge interaction with the extra Asp residue (D460 in hERG) that lies just above F463 (Fig. 2, A–C). This transition is facilitated by changes in side-chain rotamers of R534 and F463 as the planar Arg guanidine group rotates past the F463 ring, and the availability of D460 as a counterion for the R534 guanidine (Fig. 2). Movement of an Arg guanidine past the Phe side chain of the CTC is similar to that described in steered MD of an isolated VS domain (12).Open in a separate windowFigure 2Movement of the R534 side chain across the CTC in chain a of the hERG model simulation (A). Similar transitions are observed in chains a and b (panels B and C), but not chains c (D) or d (not shown), where the R534 side chain remains close to D466. In all subunits of the Kv1.2/2.1 chimera simulation, charge pairing of the starting structure (Fig. 1B) was maintained throughout (e.g., panel E and see Fig. S2 in the Supporting Material). (Black and blue lines) Distances from the Arg CZ or Lys ε atom to the two O atoms, respectively, of Asp or Glu.Mason et al. (13) have shown, using neutron scattering, that the low charge density guanidine cation (Gdm⁺) corresponding to the Arg side chain is poorly hydrated above and below the molecular plane. This property may underlie the universal preference for Arg (over Lys) in voltage sensor charge transfer. Although the poorly-hydrated surfaces of Gdm⁺ interact favorably with nonpolar (especially planar) surfaces (14,15), Gdm⁺ retains in-plane hydrogen bonding (13). In the transition of R534 across the CTC, in-plane solvation of the guanidine side chain is provided initially by D466, D501, and water molecules below the CTC, and during and after the transition by D501 and D460 side chains and waters above the CTC (Fig. 3, A and B). Complete transfer of the R534 side chain across the CTC was not observed, but would be expected to involve movement of the guanidine group away from H-bonding distance with D501.Open in a separate windowFigure 3In-plane solvation of R534 guanidine in the charge transfer center during the hERG model MD (A). (Dotted lines) H-bond distances of <2.5 Å. The right-hand group consists of top-down (B) and end-on (C) views of the distribution of oxygen atoms around the side chain of hERG R534 at 20-ns intervals during MD (subunit a). (D) End-on view of equivalent atom distributions around the K302 side chain during the Kv1.2/2.1 chimera MD (subunit c). (Red spheres, water O; pink, Asp OD1 and OD2; purple:, Glu OE1 and OE2.)The atom distribution around the R534 side chain during MD (Fig. 3, B and C) conforms to the experimental Gdm⁺ hydration structure (13), with H-bonding to waters and side-chain Asp O atoms exclusively in the guanidine plane. The passage of Gdm⁺ through the CTC is facilitated by the hydrophobic nature of Gdm⁺ above and below the molecular plane (13), which allows interaction with the nonpolar groups (especially F463) in the CTC (Fig. 3 A and see Fig. S3). This contrasts with the solvation properties of the Lys amino group (e.g., K302 of the Kv1.2/2.1 chimera (Fig. 1), which has a spherical distribution of H-bonding and charge-neutralizing oxygen atoms (Fig. 3 D and see Fig. S4).To further test these interpretations, we ran additional MD simulations of the isolated hERG VS domain model and an R534K mutant in a hydrated POPC bilayer. Again, the R534 side chain entered the CTC in the wild-type model simulation whereas the K534 side chain did not (see Fig. S5). Inspection of the atom distributions in Fig. 3 D (and see Fig. S4) indicates that the pocket below the conserved Phe of the CTC is particularly favorable for a Lys side chain, with waters and acidic side chains that satisfy the spherical solvation requirements of the terminal amino group, and nonpolar side chains that interact with the aliphatic part of the side chain.The occurrence of transitions of the R534 side chain through the CTC in the hERG model, in the absence of a change in membrane potential, indicates a relaxation from a less-stable starting structure. However, the path of the R534 side chain provides useful molecular-level insight into the nature of charge transfer in voltage sensors. How do these observations accord with broader evidence of charge transfer in voltage-sensitive channels in general, and hERG in particular? Studies with fluorinated analogs of aromatic side chains equivalent to F463 of hERG or F233 of the chimera indicate the absence of a significant role for cation-π interactions involving the CTC aromatic group in K⁺ and Na_v channels, although a planar side chain is preferred in some cases (1,16). In hERG, F463 can be replaced by M, L, or V with small effects on channel gating (17), indicating that the hERG CTC requires only a bulky nonpolar side chain to seal the hydrophobic center of the VS and allow passage of the Arg side chain through the CTC. Both absence of requirement for cation-π interactions, and accommodation of nonplanar hydrophobic side chains in a functional hERG CTC, are broadly consistent with the interpretation that it is the poorly-hydrated nature of the Arg guanidine group above and below the molecular plane (together with its tenacious proton affinity (18)) that governs its role in carrying gating charge in voltage sensors.While the simulations suggest that R534 may interact with D460 in the open channel state, the possibility that the extra carboxylate side chain above the CTC might facilitate gating charge transfer is seemingly inconsistent with the slow activation of hERG, although hERG D460C does activate even more slowly than the WT channel (9). However, S4 movement in hERG occurs in advance of channel opening (19), and slow gating is partly mediated by interactions involving hERG cytoplasmic domains (20); thus, slow S4 movement may not be an inherent property of the hERG voltage sensor. Recent studies show that when hERG gating is studied at very low [Ca²⁺] (50 μM) and low [H⁺] (pH 8.0), the channel is strongly sensitized in the direction of the open state; this effect is reduced in hERG D460C (and hERG D509C) (10). These observations support a role for the extra hERG Asp residues in binding Ca²⁺ (and H⁺) (10), allowing the channel to be allosterically responsive to changes in pH and [Ca²⁺]. A true comparison of a hERG model with experimental channel gating might involve studies on a channel lacking cytoplasmic domains that modulate gating, and using conditions (high pH and low [Ca²⁺]) that leave the eag-specific Asp residues unoccupied. This could reveal the inherent current-voltage relationships and kinetics of the hERG voltage sensor.     

Identification and Characterization of the Missing Pyrimidine Reductase in the Plant Riboflavin Biosynthesis Pathway

Riboflavin (vitamin B₂) is the precursor of the flavin coenzymes flavin mononucleotide and flavin adenine dinucleotide. In Escherichia coli and other bacteria, sequential deamination and reduction steps in riboflavin biosynthesis are catalyzed by RibD, a bifunctional protein with distinct pyrimidine deaminase and reductase domains. Plants have two diverged RibD homologs, PyrD and PyrR; PyrR proteins have an extra carboxyl-terminal domain (COG3236) of unknown function. Arabidopsis (Arabidopsis thaliana) PyrD (encoded by At4g20960) is known to be a monofunctional pyrimidine deaminase, but no pyrimidine reductase has been identified. Bioinformatic analyses indicated that plant PyrR proteins have a catalytically competent reductase domain but lack essential zinc-binding residues in the deaminase domain, and that the Arabidopsis PyrR gene (At3g47390) is coexpressed with riboflavin synthesis genes. These observations imply that PyrR is a pyrimidine reductase without deaminase activity. Consistent with this inference, Arabidopsis or maize (Zea mays) PyrR (At3g47390 or GRMZM2G090068) restored riboflavin prototrophy to an E. coli ribD deletant strain when coexpressed with the corresponding PyrD protein (At4g20960 or GRMZM2G320099) but not when expressed alone; the COG3236 domain was unnecessary for complementing activity. Furthermore, recombinant maize PyrR mediated NAD(P)H-dependent pyrimidine reduction in vitro. Import assays with pea (Pisum sativum) chloroplasts showed that PyrR and PyrD are taken up and proteolytically processed. Ablation of the maize PyrR gene caused early seed lethality. These data argue that PyrR is the missing plant pyrimidine reductase, that it is plastid localized, and that it is essential. The role of the COG3236 domain remains mysterious; no evidence was obtained for the possibility that it catalyzes the dephosphorylation that follows pyrimidine reduction.Riboflavin is the substrate for biosynthesis of the essential flavocoenzymes FMN and FAD, which occur in all kingdoms of life and have roles in diverse redox reactions as well as in other processes such as DNA repair, light sensing, and bioluminescence (Fischer and Bacher, 2005). Plants and many microorganisms can synthesize riboflavin, but humans and other animals cannot, so they must obtain it from the diet (Powers, 2003). Plant foods are important sources of riboflavin for humans, and the riboflavin pathway is a target for engineering biofortified crops (Fitzpatrick et al., 2012).Riboflavin biosynthesis proceeds via the same pathway in bacteria and plants (Fischer and Bacher, 2005; Roje, 2007). This pathway starts from GTP, which is converted by GTP cyclohydrolase II (named RibA in Escherichia coli) to the pyrimidine derivative 2,5-diamino-6-ribosylamino-4(3H)-pyrimidinone 5′-P. Deamination of the pyrimidine ring then yields 5-amino-6-ribosylamino-2,4(1H,3H)-pyrimidinedione 5′-P, and subsequent reduction of the ribosyl moiety gives 5-amino-6-ribitylamino-2,4(1H,3H)-pyrimidinedione 5′-P. After dephosphorylation, this product is condensed with 3,4-dihydroxy-2-butanone 4-P to give 6,7-dimethyl-8-ribityllumazine, whose dismutation yields riboflavin. Figure 1 shows the first four steps of this pathway.Open in a separate windowFigure 1.The first four steps of the riboflavin biosynthesis pathway in bacteria and plants. The enzymes involved are GTP cyclohydrolase II (RibA), pyrimidine deaminase (Deam), pyrimidine reductase (Red), and a specific phosphatase (Pase). Enzymes for which the plant genes are not known are colored red. Intermediates are as follows: 1, 2,5-diamino-6-ribosylamino-4(3H)-pyrimidinone 5′-P; 2, 5-amino-6-ribosylamino-2,4(1H,3H)-pyrimidinedione 5′-P; 3, 5-amino-6-ribitylamino-2,4(1H,3H)-pyrimidinedione 5′-P; 4, 5-amino-6-ribitylamino-2,4(1H,3H)-pyrimidinedione.In E. coli, the deamination and reduction steps are catalyzed by a single bifunctional enzyme, RibD, which has N-terminal deaminase and C-terminal reductase domains that retain their respective activities when expressed separately (Richter et al., 1997; Magalhães et al., 2008). The situation in plants seems superficially similar but is in fact more complex (Gerdes et al., 2012). The bidomain bacterial RibD protein has two types of homologs in plants (Fischer et al., 2004; Chatwell et al., 2006; Chen et al., 2006), here called PyrD and PyrR, both with apparent deaminase and reductase domains (Fig. 2A). Only PyrD, represented by At4g20960, has been studied biochemically; it was found to have pyrimidine deaminase but not reductase activity (Fischer et al., 2004). The function of PyrR, represented by At3g47390, is unknown, although it has been inferred to have reductase activity (Chatwell et al., 2006; Chen et al., 2006; Ouyang et al., 2010) and perhaps to lack deaminase activity (Ouyang et al., 2010). Another mystery surrounding PyrR proteins is the presence of an extra C-terminal domain of unknown function (COG3236 in the Clusters of Orthologous Groups database; Fig. 2A); this domain occurs as a stand-alone protein in many bacteria. One possibility is that it catalyzes the dephosphorylation that follows the pyrimidine reduction step in the pathway (Fig. 1). The phosphatase involved is most likely substrate specific, but it has not been identified in plants or any other organism (Roje, 2007; Gerdes et al., 2012), and genes for enzymes in the same pathway, especially for successive steps, are quite commonly fused (Suhre, 2007). A mutation (phs1) that deleted the COG3236 domain from Arabidopsis (Arabidopsis thaliana) PyrR resulted in a photosensitive phenotype that could be rescued by supplied FAD (Ouyang et al., 2010).Open in a separate windowFigure 2.Structure and phylogeny of plant PyrD and PyrR proteins. A, Domain architectures. The examples shown are Arabidopsis At4g20960 and At3g47390; the predicted plastid targeting peptide (TP) varies in length between species. B, Phylogenetic tree of PyrD and PyrR proteins. Sequences were aligned with ClustalW; the tree was built by the neighbor-joining method with MEGA5. Bootstrap values (%) for 1,000 replicates are next to the nodes. Only the tree topology is shown. Note that the PyrD proteins of green algae (underlined) lack a reductase domain. C, Alignments showing the conservation of the zinc-binding residues (arrowheads) in the deaminase domain of PyrD but not PyrR proteins and the conservation of the predicted substrate-binding residues (asterisks) in the reductase domain of PyrR but not PyrD proteins. The deaminase sequences correspond to residues 45 to 85 of B. subtilis RibD (synonym RibG); the reductase sequences correspond to residues 150 to 210 and (separated by dots) 288 to 292 of B. subtilis RibD. Identical zinc- or substrate-binding residues are black, and conservative replacements are gray. Dashes indicate gaps that maximize the alignment.The plant riboflavin synthesis pathway is considered to be plastidial (Roje, 2007), but this location is based almost solely on bioinformatics and high-throughput proteome analyses (Gerdes et al., 2012). In only one case is there more definitive experimental support: in vitro chloroplast import data for the pathway’s penultimate enzyme, 6,7-dimethyl-8-ribityllumazine synthase (Jordan et al., 1999). Similarly, clear genetic support for the function of most plant riboflavin synthesis enzymes is lacking (Gerdes et al., 2012), the exception being an Arabidopsis RibA homolog (Hedtke and Grimm, 2009).The work reported here established, using maize (Zea mays) and Arabidopsis, that PyrR is indeed the missing pyrimidine reductase, that it lacks deaminase activity, and that its COG3236 domain is not essential for pyrimidine reductase activity and most likely lacks phosphatase activity. We also demonstrated the import of PyrR and PyrD into chloroplasts in vitro and confirmed that the gene for PyrR is essential.     

Chlorophyll f-driven photosynthesis in a cavernous cyanobacterium

Chlorophyll (Chl) f is the most recently discovered chlorophyll and has only been found in cyanobacteria from wet environments. Although its structure and biophysical properties are resolved, the importance of Chl f as an accessory pigment in photosynthesis remains unresolved. We found Chl f in a cyanobacterium enriched from a cavernous environment and report the first example of Chl f-supported oxygenic photosynthesis in cyanobacteria from such habitats. Pigment extraction, hyperspectral microscopy and transmission electron microscopy demonstrated the presence of Chl a and f in unicellular cyanobacteria found in enrichment cultures. Amplicon sequencing indicated that all oxygenic phototrophs were related to KC1, a Chl f-containing cyanobacterium previously isolated from an aquatic environment. Microsensor measurements on aggregates demonstrated oxygenic photosynthesis at 742 nm and less efficient photosynthesis under 768- and 777-nm light probably because of diminished overlap with the absorption spectrum of Chl f and other far-red absorbing pigments. Our findings suggest the importance of Chl f-containing cyanobacteria in terrestrial habitats.The textbook concept that oxygenic phototrophs primarily use radiation in the visible range (400–700 nm) has been challenged by several findings of unique cyanobacteria and chlorophylls (Chl) over the past two decades (Miyashita et al., 1996; Chen et al., 2010; Croce and van Amerongen, 2014) Unicellular cyanobacteria in the genus Acaryochloris primarily employ Chl d for oxygenic photosynthesis at 700–720 nm (Miyashita et al., 1996) and thrive in shaded habitats with low levels of visible light but replete of near-infrared radiation (NIR, >700 nm, Kühl et al., 2005; Behrendt et al., 2011, 2012). Furthermore, Chl f was recently discovered in filamentous (Chen et al., 2010; Airs et al., 2014; Gan et al., 2014) and unicellular cyanobacteria (Miyashita et al., 2014), enabling light harvesting even further into the NIR region up to ∼740 nm, often aided by employing additional far-red light-absorbing pigments such as Chl d and phycobiliproteins (Gan et al., 2014). Whereas the biochemical structure (Willows et al., 2013) and biophysical properties (Li et al., 2013; Tomo et al., 2014) of Chl f have been studied in detail, the actual importance of this new chlorophyll for photosynthesis is hardly explored (Li et al., 2014).Chlorophyll f has been found in cyanobacteria originating from aquatic/wet environments: the filamentous Halomicronema hongdechloris from stromatolites in Australia (Chen et al., 2012), a unicellar morphotype (Strain KC1) from Lake Biwa in Japan (Akutsu et al., 2011; Miyashita et al., 2014) and a filamentous Leptolyngbya sp. strain (JSC-1, Gan et al., 2014) from a hot-spring and in a unicellular Chlorogloeopsis fritschii strain from rice paddies (Airs et al., 2014). In this study, we report on a unicellular Chl f-containing cyanobacterium originating from a wet cavernous habitat and demonstrate its capability of NIR-driven oxygenic photosynthesis. Enrichments of the new cyanobacterium were obtained from a dense dark green-blackish biofilm dominated by globular morphotypes of Nostocaceae growing on moist limestone outside Jenolan Caves, NSW, Australia. The sampling site was heavily shaded even during mid-day with low irradiance levels of 400- to 700-nm light varying from 0.5 to 5 μmol photons m⁻² s⁻¹. Biofilms were carefully scraped off the substratum and kept in shaded zip-lock bags in a moist atmosphere until further processing. Samples were then incubated at 28 °C in a f/2 medium under NIR at 720 nm (∼10 μmol photons m⁻² s⁻¹) yielding conspicuous green cell aggregates after several months of incubation. Repeated transfer of the aggregates into fresh medium resulted in a culture predominated by green cell clusters (Figure 1a), exhibiting orange-red fluorescence upon excitation with blue light (Figure 1b). Transmission electron microscopy revealed that the green clusters consisted of slightly elongated unicellular cyanobacteria (∼1- to 2-μm wide and ∼2- to 3-μm long), with stacked thylakoids and embedded in a joint polymer matrix (Figure 1c). Hyperspectral microscopy (Kühl and Polerecky, 2008) of the clusters revealed distinct troughs in the transmission spectra at absorption maxima indicative of Chl a (675–680 nm) and Chl f (∼720 nm; Figure 1d, red line). In situ spectral irradiance measurements at the sampling site showed strong depletion of visible wavelengths in the 480- to 710-nm range (Figure 1d, gray line), whereas highest light levels were found in the near-infrared region of the solar spectrum at 710–900 nm. The presence of Chl a and f was further confirmed in enrichment cultures using high-performance liquid chromatography-based pigment analysis (Figure 1e, Supplementary Figure S1), while no Chl d was detected. In addition, weak spectral signatures of carotenoids and phycobilins, with absorption occurring at ∼495 and 665 nm, were evident in the hyperspectral data. Cyanobacteria, including those producing Chl d/f, are known to actively remodel their pigment content in response to the available light spectrum (Stomp et al., 2007; Chen and Scheer, 2013; Gan et al., 2014) and Chl d/f has almost exclusively been found in cyanobacteria grown under far-red light and not under visible light (Kühl et al., 2005; Chen et al., 2010; Airs et al., 2014; Gan et al., 2014; Li et al., 2014; Miyashita et al., 2014). Recent work describes this acclimation response as ‘Far-Red Light photoacclimation'' (FaRLiP), which, in strain JSC-1, comprises a global change in gene expression and structural remodeling of the PSII/PSI core proteins and phycobilisome constituents (Gan et al., 2014). The extent to which this arrangement results in optimized photosynthetic performance is only known for the NIR (=710 nm)-acclimated strain JSC-1, where exposure to wavelengths >695 nm resulted in 40% higher O₂ evolution rates as compared with cells that were previously adapted to red light (645 nm; Gan et al., 2014). Yet the discrimination of actinic wavelengths and their relative effect on gross photosynthesis in Chl f-containing cells needs further investigation. Using an O₂ microsensor and the light–dark shift method (Revsbech et al., 1983) on embedded Chl f-containing aggregates, we found maximal gross photosynthesis rates (∼1.06 μmol O₂ cm⁻³ s⁻¹) to occur at irradiances of ∼250 μmol photons m⁻² s⁻¹ of 742 nm (half-bandwidth, HBW, 25 nm, Figures 2a and b) with light saturation to occur very early at ∼35 μmol photons m⁻² s⁻¹. Further red-shifted actinic light, that is, 768 nm (HBW 28 nm) and 777 nm (HBW 30 nm), yielded lower O₂ evolution rates, which, in all likelihood, are an effect of the diminished overlap with far-red light-absorbing pigments, including Chl f (Figures 2a and b). As O₂ evolution rates were measured on non-axenic cell aggregates, 16S rDNA amplicon sequencing was employed to determine the microbial diversity found within the enrichment culture. This revealed the presence of a variety of bacterial types, including anoxygenic phototrophs, yet all sequences for known oxygenic phototrophs in the data set (∼9.3% of all reads on the order level, Supplementary Figure S2) formed a single operational taxonomic unit (OTU) closely affiliated with the Chl f-containing strain KC1 (Miyashita et al., 2014, Figure 2c).Open in a separate windowFigure 1Imaging and pigment analysis of Chl f-containing cyanobacteria isolated from a cavernous low-light environment. (a) Representative bright field microscope image of cultured cells grown under 720 nm NIR. (b) Fluorescence image of the same cells as in a, excited at 450–490 nm, with emission being detected at >510 nm. (c) Transmission electron microscopy of a Chl f-containing cyanobacterium with densely stacked thylakoid membranes. (d) Transmittance spectrum of cell aggregate determined by hyperspectral imaging (red line). Ambient light conditions at the site of isolation (gray line), as measured by a spectroradiometer. Note the Chl f-specific in vivo absorption at ∼720 nm in the transmittance spectrum (dotted line). Small insert picture denotes the cells and area of interest (black arrow) from which the spectrum was taken. (e) In vitro absorption spectrum of Chl f extracted from enrichment cultures and analyzed via high-performance liquid chromatography. The two Chl f-specific absorption peaks (404 and 704 nm in acetone:MeOH solvent) are indicated.Open in a separate windowFigure 2Taxonomic affiliation and O₂ evolution of Chl f-containing cells as determined by O₂ microelectrode measurements and 16 S rDNA amplicon sequencing. (a) Emission spectra of narrow-band light-emitting diodes (LEDs) used in this study, with peak emissions at 742, 768 and 777 nm indicated by a–c, respectively. (b) Gross photosynthesis measured via an O₂ microsensor placed in a clump of agarose-embedded Chl f-containing cells. Different NIR irradiance was administered by the LEDs in a and by altering the distance of the LEDs to the embedded cells. (c) Phylogenetic affiliation of known Chl f and/or Chl d-containing cyanobacteria (highlighted in gray) and their respective habitat/place of isolation. Taxonomy was determined by clustering all known oxygenic phototrophs found in enrichment cultures from this study (at order level) into a single OTU (=292 bp length, see Supplementary Materials for details). Phylogeny was inferred using Maximum-likelihood in conjunction with the GTR +I +G nucleotide substitution model, tree stability was tested using bootstrapping with 100 replicates. The analysis involved 39 nucleotide sequences each truncated to a length of 292 bp. Here, the green-sulphur bacterium Chlorobium tepidum TLS was chosen as the outgroup.This advocates that cells from our enrichment culture are related to KC1 cells and supports, in conjunction with further morphological-, physiological- and ultrastructural evidence, that Chl f is extending the usable light spectrum for oxygenic photosynthesis in a cavernous low-light environment. Given the lifestyle and known habitats of recognized Chl d/f-producing cyanobacteria (Figure 2c), we propose that many, if not all, surface-associated cyanobacteria are intrinsically capable of producing far-red light-absorbing pigments and to actively employ them in oxygenic photosynthesis as a result of FaRLiP or similar, yet unknown, mechanisms.     

Pharmacological targeting of guanosine monophosphate synthase suppresses melanoma cell invasion and tumorigenicity

Malignant melanoma possesses one of the highest metastatic potentials among human cancers. Acquisition of invasive phenotypes is a prerequisite for melanoma metastases. Elucidation of the molecular mechanisms underlying melanoma invasion will greatly enhance the design of novel agents for melanoma therapeutic intervention. Here, we report that guanosine monophosphate synthase (GMPS), an enzyme required for the de novo biosynthesis of GMP, has a major role in invasion and tumorigenicity of cells derived from either BRAF^V600E or NRAS^Q61R human metastatic melanomas. Moreover, GMPS levels are increased in metastatic human melanoma specimens compared with primary melanomas arguing that GMPS is an attractive candidate for anti-melanoma therapy. Accordingly, for the first time we demonstrate that angustmycin A, a nucleoside-analog inhibitor of GMPS produced by Streptomyces hygroscopius efficiently suppresses melanoma cell invasion in vitro and tumorigenicity in immunocompromised mice. Our data identify GMPS as a powerful driver of melanoma cell invasion and warrant further investigation of angustmycin A as a novel anti-melanoma agent.Malignant melanoma is one of the most aggressive types of human cancers. Its ability to metastasize in combination with resistance to conventional anticancer chemotherapy makes melanoma extremely difficult to cure, and the median survival of patients with metastatic melanoma is 8.5 months.^{1, 2, 3} A better understanding of the biology behind melanoma aggressiveness is imperative to facilitate the development of novel anti-melanoma strategies.Melanoma and other cancers cells have been shown to strongly rely on de novo nucleotide biosynthesis^{4, 5} and often overexpress several biosynthetic enzymes involved in these pathways.^{6, 7, 8, 9} Recently, we have identified a fundamental connection between melanoma invasion and biosynthesis of guanylates,⁸ suggesting that distortion of the guanylate metabolism facilitates melanoma progression.Guanosine monophosphate reductase (GMPR) reduces GMP to one of its precursors, inosine monophosphate (IMP), and depletes intracellular GTP pools (Figure 1a). We have recently demonstrated that GMPR suppresses melanoma cell invasion and growth of human melanoma cell xenografts. These findings tightly linked guanylate production to the invasive potential of melanoma cells.⁸Open in a separate windowFigure 1GMPS contributes to the invasive capability of melanoma cells. (a) Simplified schematic of the metabolic pathway for guanylates production. (b) SK-Mel-103 and SK-Mel-28 cells were transduced with a control vector or two independent shRNAs to GMPS and tested for invasion through Matrigel (left panel). Where indicated, cells were incubated with 100 μM guanosine for 24 h before the assay and guanosine supplementation was maintained throughout the experimental procedure. The data represent the average ± S.E.M. of at least two independent experiments. GMPS suppression was verified by immunoblotting (right panel). (c) Cells transduced as in (a) were plated on coverslips coated with FITC-conjugated gelatin. After 16 h cells were fixed with 4% PFA and stained for actin (rhodamine-conjugated phallodin) and nuclei (Hoechst). Where indicated, cells were incubated with 100 μM guanosine for 24 h before the assay and guanosine supplementation was maintained throughout the experimental procedure. At least 25 cells/sample were imaged to assess the number of cells with gelatin degradation. The data represent the average ± S.E.M. of two independent experiments. *P<0.05, **P<0.001 compared with control; ^#P<0.05, ^##P<0.001 compared with untreated cells. Statistics performed with Student''s t-Test. See also Supplementary Figure S1Of the several enzymes involved in guanylate biosynthesis, inositol monophosphate dehydrogenases 1 and 2 (IMPDH-1, -2), functional antagonists of GMPR (Figure 1a), have been targeted clinically with several drugs including the most specific one, mycophenolic acid (MPA) and its salt, mycophenolate mofetil (MMF).^{10, 11, 12, 13} Nonetheless, prior studies demonstrated that MPA possesses poor anti-tumor activity,^{14, 15} and it is primarily used as an immunosuppressing agent in organ transplantation.^{10, 11, 12}GMP synthase (GMPS) is the other functional antagonist of GMPR. GMPS catalyzes the amination of xanitol monophosphate (XMP) to GMP to promote GTP synthesis (Figure 1a).^{16, 17} Most of the studies on GMPS have been performed in bacteria, yeast, and insects, where GMPS has been shown to have a key role in sporulation, pathogenicity, and axon guidance, respectively.^{18, 19, 20} Mammalian GMPS has been the subject of several studies addressing its unconventional (GMP-unrelated) roles in the regulation of activity of ubiquitin-specific protease 7 (USP7).^{21, 22, 23, 24} However, because of the newly revealed importance of guanylate metabolism in control of melanoma cell invasion and tumorigenicity,⁸ GMPS emerges as an attractive target for anti-cancer therapy.In the late 1950s, a specific inhibitor of bacterial GMPS, angustmycin A (also known as decoyinine), has been isolated from Streptomyces hygroscopius as a potential antibiotic with sporulation-inducing activity in Bacillus subtilis.^{25, 26, 27, 28, 29} Its anti-tumor activity has never been experimentally explored. In the current study, we investigated the role of GMPS in regulation of melanoma invasion and tumorigenicity, and explored the possibility of targeting GMPS by angustmycin A as a novel anti-melanoma strategy.     

Molecular Flow Quantified beyond the Diffraction Limit by Spatiotemporal Image Correlation of Structured Illumination Microscopy Data

We combine total internal reflection fluorescence structured illumination microscopy with spatiotemporal image correlation spectroscopy to quantify the flow velocities and directionality of filamentous-actin at the T cell immunological synapse. These techniques demonstrate it is possible to image retrograde flow of filamentous-actin at superresolution and provide flow quantification in the form of velocity histograms and flow vector maps. The flow was found to be retrograde and radially directed throughout the periphery of T-cells during synapse formation.Many biological processes are now being visualized with the use of superresolution fluorescence microscopy techniques. However, localization-based techniques primarily rely on fixed or slow moving samples to permit the collection of structural information. The 10-fold gains in resolution afforded by these superresolution techniques are usually possible through sacrificing the factors that originally made microscopy such a powerful tool: the ability to image live cells. In the case of stimulated emission depletion imaging, the scanning approach associated with this technique may fail to detect faster molecular events when imaging whole cellular regions.Structured illumination microscopy (SIM) is an alternative to these methods (1). It increases the resolution of conventional fluorescence microscopy twofold; it has the advantage of using a wide-field system, providing fast acquisition speeds of whole cells with relatively low laser powers; and it is compatible with standard fluorophores. By using a physical grating to produce interference patterns from a laser, periodic illumination is created. This patterned illumination causes information from higher spatial frequencies to be downmodulated (i.e., shifted) into the optical transfer function (support region) of the lens, resulting in higher-resolution spatial information being captured than is ordinarily obtainable.To quantify the directional motion of intracellular molecules, spatiotemporal image correlation spectroscopy (STICS (2)) was applied. Using spatial image correlation in time, STICS measures the similarity of pixels with those surrounding in lagging frames via a correlation function. The correlation function provides information on both flow velocities and directionality, while discounting static structures through the immobile object filter, achieved by subtracting a moving average of pixel intensities.The formation of an immunological synapse between T cells and antigen-presenting cells is a process requiring many dynamic (3) and subdiffraction-limited clustering events (4–6) to take place. The polymerization of actin is important for the spreading of cells over their target antigen-presenting cells (7), as well as cell mobility and migration (8). Retrograde flow of densely meshed cortical actin is observed at the basal membrane of synapse-forming T cells, where it may have a role in the corralling and clustering of signaling molecules at the plasma membrane (9), as well as at the leading edge of migrating cells (10). Filamentous actin is an extremely dynamic (7), densely packed, and thin (7-nm) structure (11,12).Here, we perform STICS on SIM data acquired on a total internal reflection fluorescence (TIRF) microscope system, which generated an evanescent field of 75-nm depth for excitation. To our knowledge, this is the first demonstration of an image correlation approach to quantify molecular dynamics on subresolution length scales using wide-field microscopy. To demonstrate the technique, we analyze two-dimensional actin flows in CD4⁺ T cells during immunological synapse formation, performed after cross-linking of antigen T cell receptors on a coverslip coated with specific antibodies.Fig. 1 a shows a schematic of the TIRF SIM setup. Excitation light (488 nm) passes through a polarizing module and then a phase-grating block, producing diffracted beams. These are then passed through a diffraction filter module to isolate the −1 and +1 order laser beams. These first-order laser beams are angled through the objective to produce total internal reflection conditions at the glass-water interface. The two evanescent waves interfere at the sample, producing structured illumination. The setup then produces lateral and rotational shifts through three orientations, producing nine raw images containing higher spatial frequencies than can normally be acquired by an objective using standard light microscopy. Fig. 1 b demonstrates the increased resolution obtained from TIRF SIM. Shown are the collected Fourier frequencies compared to those of a conventional microscope (dotted red line). Resolution was also measured using sparse 100-nm diameter fluorescent beads. Fig. 1 c shows a magnified image of these beads from which a line profile was obtained (yellow arrow). The full width at half-maximum of this profile (Fig. 1 d) gives a lateral resolution for the system of 120 nm.Open in a separate windowFigure 1(a) Schematic of the TIRF SIM setup. (b) Demonstration of the doubling of spatial resolution of collected frequencies through a Fourier transform (superimposed red circle demonstrating regular spatial frequency limits). (c) SIM reconstructed image of 100-nm bead (scale bar 0.5 μm). (d) (Plotted line) Bead showing full width at half-maximum of 120 nm.We then applied STICS analysis to quantify actin flow in T cell synapses acquired using TIRF SIM (Fig. 2). Fig. 2 a shows a schematic of the STICS analysis. From the raw data, immobile objects are first filtered by subtracting a moving average of the pixel values. Vector maps were obtained from correlation analysis of the time-series as previously published in Hebert et al. (2) and Brown et al. (13). Fig. 2 b shows a reconstructed TIRF SIM image of a mature T cell immunological synapse, representative of a time-point derived from the time series acquired at 1.28 fps (see Movie S1 in the Supporting Material). From this reconstructed image, two representative regions have been selected. In these regions, pseudo-colored actin flow vectors are overlaid onto the fluorescence intensity image. These range in magnitude from 0.01 μm/min (blue) to 5.61 μm/min (red). It can be observed that all flow vectors are directed radially toward the synapse center. A histogram of this flow is shown in Fig. 2 c. The histogram shows a peak retrograde flow velocity of 1.91 ± 1.27 μm/min. These data are representative of n = 7 T-cell synapses imaged by TIRF SIM.Open in a separate windowFigure 2(a) STICS analysis, performed by isolating mobile from immobile structures through a moving average filter (i) and binning a subset of pixels into blocks of superpixels (ii); the STICS software correlates spatial fluorescence fluctuations through time (iii). The code then outputs vector maps showing directionality and flow velocities. (b) TIRF SIM image of actin flow in a T cell 5 min after contact with a stimulatory coverslip. (Zoomed regions) Retrograde actin flow at the synapse periphery. (c) Histograms showing flow speed statistics of vectors from T-cell synapses (n = 7).