UPLC-ESI-TOFMS-based metabolomics and gene expression dynamics inspector self-organizing metabolomic maps as tools for understanding the cellular response to ionizing radiation

Global transcriptomic and proteomic profiling platforms have yielded important insights into the complex response to ionizing radiation (IR). Nonetheless, little is known about the ways in which small cellular metabolite concentrations change in response to IR. Here, a metabolomics approach using ultraperformance liquid chromatography coupled with electrospray time-of-flight mass spectrometry was used to profile, over time, the hydrophilic metabolome of TK6 cells exposed to IR doses ranging from 0.5 to 8.0 Gy. Multivariate data analysis of the positive ions revealed dose- and time-dependent clustering of the irradiated cells and identified certain constituents of the water-soluble metabolome as being significantly depleted as early as 1 h after IR. Tandem mass spectrometry was used to confirm metabolite identity. Many of the depleted metabolites are associated with oxidative stress and DNA repair pathways. Included are reduced glutathione, adenosine monophosphate, nicotinamide adenine dinucleotide, and spermine. Similar measurements were performed with a transformed fibroblast cell line, BJ, and it was found that a subset of the identified TK6 metabolites were effective in IR dose discrimination. The GEDI (Gene Expression Dynamics Inspector) algorithm, which is based on self-organizing maps, was used to visualize dynamic global changes in the TK6 metabolome that resulted from IR. It revealed dose-dependent clustering of ions sharing the same trends in concentration change across radiation doses. "Radiation metabolomics," the application of metabolomic analysis to the field of radiobiology, promises to increase our understanding of cellular responses to stressors such as radiation.

Investigation of the relationship between protein, message and inducer concentrations in recombinant E-coli cells

Chloramphenicol acetyl transferase (CAT) protein and mRNA levels in E. coli were determined following induction of a tac::cat construct by isopropyl-beta-thiogalactopyranoside (IPTG). High cat mRNA levels did not directly reflect CAT protein levels, in either shakeflask experiments or fermentations. Furthermore, concentrations of IPTG resulting in the highest levels of expression of cat mRNA, were different to those resulting in highest levels of CAT protein. The data suggest that high transcriptional activities lead to limitations at the translational level.

Bristle pattern formation in the thorax formation of Drosophila: a systems level approach

A fundamental question in developmental biology is how tissues are patterned to give rise to differentiated body structures with distinct morphologies. The Drosophila wing disc offers an accessible model to understand epithelial spatial patterning. It has been studied extensively using genetic and molecular approaches. Bristle patterns on the thorax, which arise from the medial part of the wing disc, are a classical model of pattern formation, dependent on a pre-pattern of trans-activators and –repressors. Despite of decades of molecular studies, we still only know a subset of the factors that determine the pre-pattern. We are applying a novel and interdisciplinary approach to predict regulatory interactions in this system. It is based on the description of expression patterns by simple logical relations (addition, subtraction, intersection and union) between simple shapes (graphical primitives). Similarities and relations between primitives have been shown to be predictive of regulatory relationships between the corresponding regulatory factors in other Systems, such as the Drosophila egg. Furthermore, they provide the basis for dynamical models of the bristle-patterning network, which enable us to make even more detailed predictions on gene regulation and expression dynamics. We have obtained a data-set of wing disc expression patterns which we are now processing to obtain average expression patterns for each gene. Through triangulation of the images we can transform the expression patterns into vectors which can easily be analysed by Standard clustering methods. These analyses will allow us to identify primitives and regulatory interactions. We expect to identify new regulatory interactions and to understand the basic Dynamics of the regulatory network responsible for thorax patterning. These results will provide us with a better understanding of the rules governing gene regulatory networks in general, and provide the basis for future studies of the evolution of the thorax-patterning network in particular.

Scalable inference of heterogeneous reaction kinetics from pooled single-cell recordings.

Mathematical methods combined with measurements of single-cell dynamics provide a means to reconstruct intracellular processes that are only partly or indirectly accessible experimentally. To obtain reliable reconstructions, the pooling of measurements from several cells of a clonal population is mandatory. However, cell-to-cell variability originating from diverse sources poses computational challenges for such process reconstruction. We introduce a scalable Bayesian inference framework that properly accounts for population heterogeneity. The method allows inference of inaccessible molecular states and kinetic parameters; computation of Bayes factors for model selection; and dissection of intrinsic, extrinsic and technical noise. We show how additional single-cell readouts such as morphological features can be included in the analysis. We use the method to reconstruct the expression dynamics of a gene under an inducible promoter in yeast from time-lapse microscopy data.

Molecular characterization of the contribution of autophagy to antigen presentation using quantitative proteomics

L’autophagie est une voie hautement conservée de dégradation lysosomale des constituants cellulaires qui est essentiel à l’homéostasie cellulaire et contribue à l’apprêtement et à la présentation des antigènes. Les rôles relativement récents de l'autophagie dans l'immunité innée et acquise sous-tendent de nouveaux paradigmes immunologiques pouvant faciliter le développement de nouvelles thérapies où la dérégulation de l’autophagie est associée à des maladies auto-immunes. Cependant, l'étude in vivo de la réponse autophagique est difficile en raison du nombre limité de méthodes d'analyse pouvant fournir une définition dynamique des protéines clés impliquées dans cette voie. En conséquence, nous avons développé un programme de recherche en protéomique intégrée afin d’identifier et de quantifier les proteines associées à l'autophagie et de déterminer les mécanismes moléculaires régissant les fonctions de l’autophagosome dans la présentation antigénique en utilisant une approche de biologie des systèmes. Pour étudier comment l'autophagie et la présentation antigénique sont activement régulés dans les macrophages, nous avons d'abord procédé à une étude protéomique à grande échelle sous différentes conditions connues pour stimuler l'autophagie, tels l’activation par les cytokines et l’infection virale. La cytokine tumor necrosis factor-alpha (TNF-alpha) est l'une des principales cytokines pro-inflammatoires qui intervient dans les réactions locales et systémiques afin de développer une réponse immune adaptative. La protéomique quantitative d'extraits membranaires de macrophages contrôles et stimulés avec le TNF-alpha a révélé que l'activation des macrophages a entrainé la dégradation de protéines mitochondriales et des changements d’abondance de plusieurs protéines impliquées dans le trafic vésiculaire et la réponse immunitaire. Nous avons constaté que la dégradation des protéines mitochondriales était sous le contrôle de la voie ATG5, et était spécifique au TNF-alpha. En outre, l’utilisation d’un nouveau système de présentation antigènique, nous a permi de constater que l'induction de la mitophagie par le TNF-alpha a entrainée l’apprêtement et la présentation d’antigènes mitochondriaux par des molécules du CMH de classe I, contribuant ainsi la variation du répertoire immunopeptidomique à la surface cellulaire. Ces résultats mettent en évidence un rôle insoupçonné du TNF-alpha dans la mitophagie et permet une meilleure compréhension des mécanismes responsables de la présentation d’auto-antigènes par les molécules du CMH de classe I. Une interaction complexe existe également entre infection virale et l'autophagie. Récemment, notre laboratoire a fourni une première preuve suggérant que la macroautophagie peut contribuer à la présentation de protéines virales par les molécules du CMH de classe I lors de l’infection virale par l'herpès simplex virus de type 1 (HSV-1). Le virus HSV1 fait parti des virus humains les plus complexes et les plus répandues. Bien que la composition des particules virales a été étudiée précédemment, on connaît moins bien l'expression de l'ensemble du protéome viral lors de l’infection des cellules hôtes. Afin de caractériser les changements dynamiques de l’expression des protéines virales lors de l’infection, nous avons analysé par LC-MS/MS le protéome du HSV1 dans les macrophages infectés. Ces analyses nous ont permis d’identifier un total de 67 protéines virales structurales et non structurales (82% du protéome HSV1) en utilisant le spectromètre de masse LTQ-Orbitrap. Nous avons également identifié 90 nouveaux sites de phosphorylation et de dix nouveaux sites d’ubiquitylation sur différentes protéines virales. Suite à l’ubiquitylation, les protéines virales peuvent se localiser au noyau ou participer à des événements de fusion avec la membrane nucléaire, suggérant ainsi que cette modification pourrait influer le trafic vésiculaire des protéines virales. Le traitement avec des inhibiteurs de la réplication de l'ADN induit des changements sur l'abondance et la modification des protéines virales, mettant en évidence l'interdépendance des protéines virales au cours du cycle de vie du virus. Compte tenu de l'importance de la dynamique d'expression, de l’ubiquitylation et la phosphorylation sur la fonction des proteines virales, ces résultats ouvriront la voie vers de nouvelles études sur la biologie des virus de l'herpès. Fait intéressant, l'infection HSV1 dans les macrophages déclenche une nouvelle forme d'autophagie qui diffère remarquablement de la macroautophagie. Ce processus, appelé autophagie associée à l’enveloppe nucléaire (nuclear envelope derived autophagy, NEDA), conduit à la formation de vésicules membranaires contenant 4 couches lipidiques provenant de l'enveloppe nucléaire où on retrouve une grande proportion de certaines protéines virales, telle la glycoprotéine B. Les mécanismes régissant NEDA et leur importance lors de l’infection virale sont encore méconnus. En utilisant un essai de présentation antigénique, nous avons pu montrer que la voie NEDA est indépendante d’ATG5 et participe à l’apprêtement et la présentation d’antigènes viraux par le CMH de classe I. Pour comprendre l'implication de NEDA dans la présentation des antigènes, il est essentiel de caractériser le protéome des autophagosomes isolés à partir de macrophages infectés par HSV1. Aussi, nous avons développé une nouvelle approche de fractionnement basé sur l’isolation de lysosomes chargés de billes de latex, nous permettant ainsi d’obtenir des extraits cellulaires enrichis en autophagosomes. Le transfert des antigènes HSV1 dans les autophagosomes a été determine par protéomique quantitative. Les protéines provenant de l’enveloppe nucléaire ont été préférentiellement transférées dans les autophagosome lors de l'infection des macrophages par le HSV1. Les analyses protéomiques d’autophagosomes impliquant NEDA ou la macroautophagie ont permis de decouvrir des mécanismes jouant un rôle clé dans l’immunodominance de la glycoprotéine B lors de l'infection HSV1. Ces analyses ont également révélées que diverses voies autophagiques peuvent être induites pour favoriser la capture sélective de protéines virales, façonnant de façon dynamique la nature de la réponse immunitaire lors d'une infection. En conclusion, l'application des méthodes de protéomique quantitative a joué un rôle clé dans l'identification et la quantification des protéines ayant des rôles importants dans la régulation de l'autophagie chez les macrophages, et nous a permis d'identifier les changements qui se produisent lors de la formation des autophagosomes lors de maladies inflammatoires ou d’infection virale. En outre, notre approche de biologie des systèmes, qui combine la protéomique quantitative basée sur la spectrométrie de masse avec des essais fonctionnels tels la présentation antigénique, nous a permis d’acquérir de nouvelles connaissances sur les mécanismes moléculaires régissant les fonctions de l'autophagie lors de la présentation antigénique. Une meilleure compréhension de ces mécanismes permettra de réduire les effets nuisibles de l'immunodominance suite à l'infection virale ou lors du développement du cancer en mettant en place une réponse immunitaire appropriée.

Density-Profile Processes Describing Biological Signaling Networks: Almost Sure Convergence to Deterministic Trajectories

We introduce jump processes in R(k), called density-profile processes, to model biological signaling networks. Our modeling setup describes the macroscopic evolution of a finite-size spin-flip model with k types of spins with arbitrary number of internal states interacting through a non-reversible stochastic dynamics. We are mostly interested on the multi-dimensional empirical-magnetization vector in the thermodynamic limit, and prove that, within arbitrary finite time-intervals, its path converges almost surely to a deterministic trajectory determined by a first-order (non-linear) differential equation with explicit bounds on the distance between the stochastic and deterministic trajectories. As parameters of the spin-flip dynamics change, the associated dynamical system may go through bifurcations, associated to phase transitions in the statistical mechanical setting. We present a simple example of spin-flip stochastic model, associated to a synthetic biology model known as repressilator, which leads to a dynamical system with Hopf and pitchfork bifurcations. Depending on the parameter values, the magnetization random path can either converge to a unique stable fixed point, converge to one of a pair of stable fixed points, or asymptotically evolve close to a deterministic orbit in Rk. We also discuss a simple signaling pathway related to cancer research, called p53 module.

A hormonal misunderstanding in Acca sellowiana embryogenesis: levels of zygotic embryogenesis do not match those of somatic embryogenesis

Endogenous levels of IAA, ABA and four types of CKs were analyzed in zygotic and indirect (ISE) and direct somatic embryogenesis of Acca sellowiana. Zygotic and somatic embryos at different developmental stages were sampled for morphological and hormonal analysis. Both embryo types showed substantial asymmetry in hormone levels. Zygotic embryos displayed a conspicuous peak of IAA in early developmental stages. The results outlined the hormonal variations occurring during zygotic and somatic embryogenesis regarding the timing, nature and hormonal status involved in both processes. The short transient pulse of IAA observed on the 3rd day in culture was suggested to be involved with the signaling for the induction of somatic embryogenesis. Fertilized ovule development was associated with increased IAA levels 21-24 days after pollination, followed by a sharp decrease in the cotyledonary stage, both in zygotic and somatic embryos. There was a prominent increase in ABA levels in cultures which generated ISE 24-30 days after pollination, a period that corresponds to the heart and torpedo stages. The levels of total CKs (Z, [9R]Z, iP and [9R]iP) were also always higher in zygotic than in somatic embryogenesis. While zygotic embryogenesis was dominated by the presence of zeatin, the somatic process, contrarily, was characterized by a large variation of the other cytokinin forms and amounts studied. The above results, when taken together, could be related to the previously observed high frequency formation of anomalous somatic embryos formed in A. sellowiana, as well as to their low germination ability.

Radiation metabolomics. 1. Identification of minimally invasive urine biomarkers for gamma-radiation exposure in mice

Gamma-radiation exposure has both short- and long-term adverse health effects. The threat of modern terrorism places human populations at risk for radiological exposures, yet current medical countermeasures to radiation exposure are limited. Here we describe metabolomics for gamma-radiation biodosimetry in a mouse model. Mice were gamma-irradiated at doses of 0, 3 and 8 Gy (2.57 Gy/min), and urine samples collected over the first 24 h after exposure were analyzed by ultra-performance liquid chromatography-time-of-flight mass spectrometry (UPLC-TOFMS). Multivariate data were analyzed by orthogonal partial least squares (OPLS). Both 3- and 8-Gy exposures yielded distinct urine metabolomic phenotypes. The top 22 ions for 3 and 8 Gy were analyzed further, including tandem mass spectrometric comparison with authentic standards, revealing that N-hexanoylglycine and beta-thymidine are urinary biomarkers of exposure to 3 and 8 Gy, 3-hydroxy-2-methylbenzoic acid 3-O-sulfate is elevated in urine of mice exposed to 3 but not 8 Gy, and taurine is elevated after 8 but not 3 Gy. Gene Expression Dynamics Inspector (GEDI) self-organizing maps showed clear dose-response relationships for subsets of the urine metabolome. This approach is useful for identifying mice exposed to gamma radiation and for developing metabolomic strategies for noninvasive radiation biodosimetry in humans.

Phase transitions in gene networks evolved under different selection rules

Mathematical models of gene regulation are a powerful tool for understanding the complex features of genetic control. While various modeling efforts have been successful at explaining gene expression dynamics, much less is known about how evolution shapes the structure of these networks. An important feature of gene regulatory networks is their stability in response to environmental perturbations. Regulatory systems are thought to have evolved to exist near the transition between stability and instability, in order to have the required stability to environmental fluctuations while also being able to achieve a wide variety of functions (corresponding to different dynamical patterns). We study a simplified model of gene network evolution in which links are added via different selection rules. These growth models are inspired by recent work on `explosive' percolation which shows that when network links are added through competitive rather than random processes, the connectivity phase transition can be significantly delayed, and when it is reached, it appears to be first order (discontinuous, e.g., going from no failure at all to large expected failure) instead of second order (continuous, e.g., going from no failure at all to very small expected failure). We find that by modifying the traditional framework for networks grown via competitive link addition to capture how gene networks evolve to avoid damage propagation, we also see significant delays in the transition that depend on the selection rules, but the transitions always appear continuous rather than `explosive'.

Reprint of ""Dynamics of the expression of intermediate filaments vimentin and desmin during myofibroblast differentiation after corneal injury""

Previous studies have suggested that abnormal corneal wound healing in patients after photorefractive keratectomy (PRK) is associated with the appearance of myofibroblasts in the stroma between two and four weeks after surgery. The purpose of this study was to examine potential myofibroblast progenitor cells that might express other filament markers prior to completion of the differentiation pathway that yields alpha-smooth muscle actin (SMA)-expressing myofibroblasts associated with haze localized beneath the epithelial basement membrane after PRK. Twenty-four female rabbits that had -9 diopter PRK were sacrificed at 1 week, 2 weeks, 3 weeks or 4 weeks after surgery. Corneal rims were collected, frozen at -80 degrees C, and analyzed by immunocytochemistry using anti-vimentin, anti-desmin, and anti-SMA antibodies. Double immunostaining was performed for the co-localization of SMA with vimentin or desmin with SMA. An increase in vimentin expression in stromal cells is noted as early as 1 week after PRK in the rabbit cornea. As the healing response continues at two or three weeks after surgery, many stromal cells expressing vimentin also begin to express desmin and SMA. By 4 weeks after the surgery most, if not all, myofibroblasts express vimentin, desmin and SMA. Generalized least squares regression analysis showed that there was strong evidence that each of the marker groups differed in expression over time compared to the other two (p < 0.01). Intermediate filaments - vimentin and desmin co-exist in myofibroblasts along with SMA and may play an important role in corneal remodeling after photorefractive keratectomy. The earliest precursors of myofibroblasts destined to express SMA and desmin are detectible by staining for vimentin at 1 week after surgery. (C) 2009 Elsevier Ltd. All rights reserved.

Dynamics of the expression of intermediate filaments vimentin and desmin during myofibroblast differentiation after corneal injury

Previous studies have suggested that abnormal corneal wound healing in patients after photorefractive keratectomy (PRK) is associated with the appearance of myofibroblasts in the stroma between two and four weeks after surgery. The purpose of this study was to examine potential myofibroblast progenitor cells that might express other filament markers prior to completion of the differentiation pathway that yields a-smooth muscle actin (SMA)-expressing myofibroblasts associated with haze localized beneath the epithelial basement membrane after PRK. Twenty-four female rabbits that had -9 diopter PRK were sacrificed at I week, 2 weeks, 3 weeks or 4 weeks after surgery. Corneal rims were collected, frozen at -80 degrees C, and analyzed by immunocytochemistry using anti-vimentin, anti-desmin, and anti-SMA antibodies. Double immunostaining was performed for the co-localization of SMA with vimentin or desmin with SMA. An increase in vimentin expression in stromal cells is noted as early as 1 week after PRK in the rabbit cornea. As the healing response continues at two or three weeks after surgery, many stromal cells expressing vimentin also begin to express desmin and SMA. By 4 weeks after the surgery most, if not all, myofibroblasts express vimentin, desmin and SMA. Generalized least squares regression analysis showed that there was strong evidence that each of the marker groups differed in expression over time compared to the other two (p < 0.01). Intermediate filaments - vimentin and desmin co-exist in myofibroblasts along with SMA and may play an important role in corneal remodeling after photorefractive keratectomy. The earliest precursors of myofibroblasts destined to express SMA and desmin are detectible by staining for vimentin at I week after surgery. (C) 2009 Elsevier Ltd. All rights reserved.

FABP4 Dynamics in Obesity: Discrepancies in Adipose Tissue and Liver Expression Regarding Circulating Plasma Levels.

BACKGROUND FABP4 is predominantly expressed in adipose tissue, and its circulating levels are linked with obesity and a poor atherogenic profile. OBJECTIVE In patients with a wide BMI range, we analyze FABP4 expression in adipose and hepatic tissues in the settings of obesity and insulin resistance. Associations between FABP4 expression in adipose tissue and the FABP4 plasma level as well as the main adipogenic and lipolytic genes expressed in adipose tissue were also analyzed. METHODS The expression of several lipogenic, lipolytic, PPAR family and FABP family genes was analyzed by real time PCR. FABP4 protein expression in total adipose tissues and its fractions were determined by western blot. RESULTS In obesity FABP4 expression was down-regulated (at both mRNA and protein levels), with its levels mainly predicted by ATGL and inversely by the HOMA-IR index. The BMI appeared as the only determinant of the FABP4 variation in both adipose tissue depots. FABP4 plasma levels showed a significant progressive increase according to BMI but no association was detected between FABP4 circulating levels and SAT or VAT FABP4 gene expression. The gene expression of FABP1, FABP4 and FABP5 in hepatic tissue was significantly higher in tissue from the obese IR patients compared to the non-IR group. CONCLUSION The inverse pattern in FABP4 expression between adipose and hepatic tissue observed in morbid obese patients, regarding the IR context, suggests that both tissues may act in a balanced manner. These differences may help us to understand the discrepancies between circulating plasma levels and adipose tissue expression in obesity.

Dynamics of multiple lin gene expression in Sphingomonas paucimobilis B90A in response to different hexachlorocyclohexane isomers.

Sphingomonas paucimobilis B90A is able to degrade the alpha-, beta-, gamma-, and delta-isomers of hexachlorocyclohexane (HCH). It contains the genes linA, linB, linC, linD, linE, and linR, which have been implicated in HCH degradation. In this study, dynamic expression of the lin genes was measured in chemostat-grown S. paucimobilis B90A by RNA dot blot hybridization and real-time reverse transcriptase PCR upon exposure to a pulse of different HCH isomers. Irrespective of the addition of HCH, linA, linB, and linC were all expressed constitutively. In contrast, linD and linE were induced with alpha-HCH (2 mg/liter) and gamma-HCH (7 mg/liter). A sharp increase in mRNA levels for linD and linE was observed from 10 to 45 min after the addition of alpha- or gamma-HCH. Induction of linD and linE was not detectable upon the addition of 0.7 mg of gamma-HCH per liter, although the compound was degraded by the cells. The addition of beta-HCH (5 mg/liter) or delta-HCH (20 mg/liter) did not lead to linE and linD induction, despite the fact that 50% of the compounds were degraded. This suggests that degradation of beta- and delta-HCH proceeds by a different pathway than that of alpha- and gamma-HCH.

Dynamics of Gene Expression Revealed by Comparison of Serial Analysis of Gene Expression Transcript Profiles from Yeast Grown on Two Different Carbon SourcesD⃞

We describe a genome-wide characterization of mRNA transcript levels in yeast grown on the fatty acid oleate, determined using Serial Analysis of Gene Expression (SAGE). Comparison of this SAGE library with that reported for glucose grown cells revealed the dramatic adaptive response of yeast to a change in carbon source. A major fraction (>20%) of the 15,000 mRNA molecules in a yeast cell comprised differentially expressed transcripts, which were derived from only 2% of the total number of ∼6300 yeast genes. Most of the mRNAs that were differentially expressed code for enzymes or for other proteins participating in metabolism (e.g., metabolite transporters). In oleate-grown cells, this was exemplified by the huge increase of mRNAs encoding the peroxisomal β-oxidation enzymes required for degradation of fatty acids. The data provide evidence for the existence of redox shuttles across organellar membranes that involve peroxisomal, cytoplasmic, and mitochondrial enzymes. We also analyzed the mRNA profile of a mutant strain with deletions of the PIP2 and OAF1 genes, encoding transcription factors required for induction of genes encoding peroxisomal proteins. Induction of genes under the immediate control of these factors was abolished; other genes were up-regulated, indicating an adaptive response to the changed metabolism imposed by the genetic impairment. We describe a statistical method for analysis of data obtained by SAGE.

Dynamics of active lamellae in cultured epithelial cells: effects of expression of exogenous N-ras oncogene.

We examined the functional consequences of cellular transformation of rat IAR-2 epithelial cells, by a mutant N-ras oncogene, on the dynamics of active lamellae, structures that play an important role in cell motility, adhesion, and surface-receptor capping. Lamellar activity was assessed by measuring the rate of outer-edge pseudopodial activity and by analyzing the motility of Con A-coated beads placed on lamellar surfaces with optical tweezers. Although transformation dramatically affected the shape and size of active cellular lamellae, there was little detectable effect on either pseudopodial activity or bead movement. To investigate the potential relationship between functional lamellar activity and the microtubule cytoskeleton, lamellar activity was examined in nontransformed and transformed cells treated with the microtubule-disrupting drug nocodazole. In the absence of microtubules, transformed cells were less polarized and possessed decreased rates of pseudopodial and bead motility. On the basis of these observations, it is suggested that ras-induced transformation of epithelial cells consists of two cytoskeletal modifications: overall diminished actin cytoskeletal dynamics in lamellae and reorganization of the microtubule cytoskeleton that directs pseudopodial activity to smaller polarized lamellae.