1
|
Wang S, Chakraborty S, Fu Y, Lee MP, Liu J, Waldhaus J. 3D reconstruction of the mouse cochlea from scRNA-seq data suggests morphogen-based principles in apex-to-base specification. Dev Cell 2024; 59:1538-1552.e6. [PMID: 38593801 PMCID: PMC11187690 DOI: 10.1016/j.devcel.2024.03.028] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2022] [Revised: 04/03/2023] [Accepted: 03/13/2024] [Indexed: 04/11/2024]
Abstract
In the mammalian auditory system, frequency discrimination depends on numerous morphological and physiological properties of the organ of Corti, which gradually change along the apex-to-base (tonotopic) axis of the organ. For example, the basilar membrane stiffness changes tonotopically, thus affecting the tuning properties of individual hair cells. At the molecular level, those frequency-specific characteristics are mirrored by gene expression gradients; however, the molecular mechanisms controlling tonotopic gene expression in the mouse cochlea remain elusive. Through analyzing single-cell RNA sequencing (scRNA-seq) data from E12.5 and E14.5 time points, we predicted that morphogens, rather than a cell division-associated mechanism, confer spatial identity in the extending cochlea. Subsequently, we reconstructed the developing cochlea in 3D space from scRNA-seq data to investigate the molecular pathways mediating positional information. The retinoic acid (RA) and hedgehog pathways were found to form opposing apex-to-base gradients, and functional interrogation using mouse cochlear explants suggested that both pathways jointly specify the longitudinal axis.
Collapse
Affiliation(s)
- Shuze Wang
- Department of Otolaryngology-Head and Neck Surgery, Kresge Hearing Research Institute, University of Michigan, Ann Arbor, MI 48109, USA; Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI 48109, USA
| | - Saikat Chakraborty
- Department of Otolaryngology-Head and Neck Surgery, Kresge Hearing Research Institute, University of Michigan, Ann Arbor, MI 48109, USA
| | - Yujuan Fu
- Biomedical Informatics and Medical Education, University of Washington, Seattle, WA 98195, USA
| | - Mary P Lee
- Department of Otolaryngology-Head and Neck Surgery, Kresge Hearing Research Institute, University of Michigan, Ann Arbor, MI 48109, USA
| | - Jie Liu
- Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI 48109, USA
| | - Joerg Waldhaus
- Department of Otolaryngology-Head and Neck Surgery, Kresge Hearing Research Institute, University of Michigan, Ann Arbor, MI 48109, USA.
| |
Collapse
|
2
|
Manser CL, Perez-Carrasco R. A mathematical framework for measuring and tuning tempo in developmental gene regulatory networks. Development 2024; 151:dev202950. [PMID: 38780527 PMCID: PMC11234385 DOI: 10.1242/dev.202950] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2024] [Accepted: 05/13/2024] [Indexed: 05/25/2024]
Abstract
Embryo development is a dynamic process governed by the regulation of timing and sequences of gene expression, which control the proper growth of the organism. Although many genetic programmes coordinating these sequences are common across species, the timescales of gene expression can vary significantly among different organisms. Currently, substantial experimental efforts are focused on identifying molecular mechanisms that control these temporal aspects. In contrast, the capacity of established mathematical models to incorporate tempo control while maintaining the same dynamical landscape remains less understood. Here, we address this gap by developing a mathematical framework that links the functionality of developmental programmes to the corresponding gene expression orbits (or landscapes). This unlocks the ability to find tempo differences as perturbations in the dynamical system that preserve its orbits. We demonstrate that this framework allows for the prediction of molecular mechanisms governing tempo, through both numerical and analytical methods. Our exploration includes two case studies: a generic network featuring coupled production and degradation, with a particular application to neural progenitor differentiation; and the repressilator. In the latter, we illustrate how altering the dimerisation rates of transcription factors can decouple the tempo from the shape of the resulting orbits. We conclude by highlighting how the identification of orthogonal molecular mechanisms for tempo control can inform the design of circuits with specific orbits and tempos.
Collapse
Affiliation(s)
- Charlotte L. Manser
- Department of Life Sciences, Imperial College London, South Kensington Campus, Imperial College London, London SW7 2AZ, UK
| | - Ruben Perez-Carrasco
- Department of Life Sciences, Imperial College London, South Kensington Campus, Imperial College London, London SW7 2AZ, UK
| |
Collapse
|
3
|
Abstract
The temporal coordination of events at cellular and tissue scales is essential for the proper development of organisms, and involves cell-intrinsic processes that can be coupled by local cellular signalling and instructed by global signalling, thereby creating spatial patterns of cellular states that change over time. The timing and structure of these patterns determine how an organism develops. Traditional developmental genetic methods have revealed the complex molecular circuits regulating these processes but are limited in their ability to predict and understand the emergent spatio-temporal dynamics. Increasingly, approaches from physics are now being used to help capture the dynamics of the system by providing simplified, generic descriptions. Combined with advances in imaging and computational power, such approaches aim to provide insight into timing and patterning in developing systems.
Collapse
|
4
|
Jutras-Dubé L, El-Sherif E, François P. Geometric models for robust encoding of dynamical information into embryonic patterns. eLife 2020; 9:55778. [PMID: 32773041 PMCID: PMC7470844 DOI: 10.7554/elife.55778] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2020] [Accepted: 08/07/2020] [Indexed: 12/30/2022] Open
Abstract
During development, cells gradually assume specialized fates via changes of transcriptional dynamics, sometimes even within the same developmental stage. For anterior-posterior (AP) patterning in metazoans, it has been suggested that the gradual transition from a dynamic genetic regime to a static one is encoded by different transcriptional modules. In that case, the static regime has an essential role in pattern formation in addition to its maintenance function. In this work, we introduce a geometric approach to study such transition. We exhibit two types of genetic regime transitions arising through local or global bifurcations, respectively. We find that the global bifurcation type is more generic, more robust, and better preserves dynamical information. This could parsimoniously explain common features of metazoan segmentation, such as changes of periods leading to waves of gene expressions, ‘speed/frequency-gradient’ dynamics, and changes of wave patterns. Geometric approaches appear as possible alternatives to gene regulatory networks to understand development.
Collapse
Affiliation(s)
| | - Ezzat El-Sherif
- Division of Developmental Biology, Department of Biology, Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen, Germany
| | - Paul François
- Department of Physics, McGill University, Montreal, Canada
| |
Collapse
|
5
|
Gupta S, Fancher S, Korswagen HC, Mugler A. Temporal precision of molecular events with regulation and feedback. Phys Rev E 2020; 101:062420. [PMID: 32688616 DOI: 10.1103/physreve.101.062420] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2019] [Accepted: 06/08/2020] [Indexed: 11/06/2022]
Abstract
Cellular behaviors such as migration, division, and differentiation rely on precise timing, and yet the molecular events that govern these behaviors are highly stochastic. We investigate regulatory strategies that decrease the timing noise of molecular events. Autoregulatory feedback increases noise. Yet we find that in the presence of regulation by a second species, autoregulatory feedback decreases noise. To explain this finding, we develop a method to calculate the optimal regulation function that minimizes the timing noise. The method reveals that the combination of feedback and regulation minimizes noise by maximizing the number of molecular events that must happen in sequence before a threshold is crossed. We compute the optimal timing precision for all two-node networks with regulation and feedback, derive a generic lower bound on timing noise, and discuss our results in the context of neuroblast migration during Caenorhabditis elegans development.
Collapse
Affiliation(s)
- Shivam Gupta
- Department of Physics and Astronomy, Purdue University, West Lafayette, Indiana 47907, USA
| | - Sean Fancher
- Department of Physics and Astronomy, Purdue University, West Lafayette, Indiana 47907, USA.,Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Hendrik C Korswagen
- Hubrecht Institute, Royal Netherlands Academy of Arts and Sciences and University Medical Center Utrecht, 3584 CT Utrecht, Netherlands
| | - Andrew Mugler
- Department of Physics and Astronomy, Purdue University, West Lafayette, Indiana 47907, USA
| |
Collapse
|
6
|
Morozov A, Abbott K, Cuddington K, Francis T, Gellner G, Hastings A, Lai YC, Petrovskii S, Scranton K, Zeeman ML. Long transients in ecology: Theory and applications. Phys Life Rev 2019; 32:1-40. [PMID: 31982327 DOI: 10.1016/j.plrev.2019.09.004] [Citation(s) in RCA: 62] [Impact Index Per Article: 12.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2019] [Accepted: 09/09/2019] [Indexed: 11/15/2022]
Abstract
This paper discusses the recent progress in understanding the properties of transient dynamics in complex ecological systems. Predicting long-term trends as well as sudden changes and regime shifts in ecosystems dynamics is a major issue for ecology as such changes often result in population collapse and extinctions. Analysis of population dynamics has traditionally been focused on their long-term, asymptotic behavior whilst largely disregarding the effect of transients. However, there is a growing understanding that in ecosystems the asymptotic behavior is rarely seen. A big new challenge for theoretical and empirical ecology is to understand the implications of long transients. It is believed that the identification of the corresponding mechanisms along with the knowledge of scaling laws of the transient's lifetime should substantially improve the quality of long-term forecasting and crisis anticipation. Although transient dynamics have received considerable attention in physical literature, research into ecological transients is in its infancy and systematic studies are lacking. This text aims to partially bridge this gap and facilitate further progress in quantitative analysis of long transients in ecology. By revisiting and critically examining a broad variety of mathematical models used in ecological applications as well as empirical facts, we reveal several main mechanisms leading to the emergence of long transients and hence lays the basis for a unifying theory.
Collapse
Affiliation(s)
- Andrew Morozov
- Mathematics, University of Leicester, UK; Shirshov Institute of Oceanology, Moscow, Russia
| | | | | | - Tessa Francis
- Tacoma Puget Sound Institute, University of Washington, USA
| | | | - Alan Hastings
- Environmental Science and Policy, University of California, Davis, USA; Santa Fe Institute, Santa Fe, New Mexico, USA
| | - Ying-Cheng Lai
- Electrical, Computer and Energy Engineering, Arizona State University, Tempe, USA
| | - Sergei Petrovskii
- Mathematics, University of Leicester, UK; Peoples Friendship University of Russia (RUDN University), Moscow, Russia.
| | | | | |
Collapse
|
7
|
Rudolf H, Zellner C, El-Sherif E. Speeding up anterior-posterior patterning of insects by differential initialization of the gap gene cascade. Dev Biol 2019; 460:20-31. [PMID: 31075221 DOI: 10.1016/j.ydbio.2019.04.015] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2018] [Revised: 04/25/2019] [Accepted: 04/25/2019] [Indexed: 01/28/2023]
Abstract
Recently, it was shown that anterior-posterior patterning genes in the red flour beetle Tribolium castaneum are expressed sequentially in waves. However, in the fruit fly Drosophila melanogaster, an insect with a derived mode of embryogenesis compared to Tribolium, anterior-posterior patterning genes quickly and simultaneously arise as mature gene expression domains that, afterwards, undergo slight posterior-to-anterior shifts. This raises the question of how a fast and simultaneous mode of patterning, like that of Drosophila, could have evolved from a rather slow sequential mode of patterning, like that of Tribolium. In this paper, we propose a mechanism for this evolutionary transition based on a switch from a uniform to a gradient-mediated initialization of the gap gene cascade by maternal Hb. The model is supported by computational analyses and experiments.
Collapse
Affiliation(s)
- Heike Rudolf
- Department of Biology, Friedrich-Alexander-Universität Erlangen-Nürnberg, Staudtstrasse 5, Erlangen, 91058, Germany
| | - Christine Zellner
- Department of Biology, Friedrich-Alexander-Universität Erlangen-Nürnberg, Staudtstrasse 5, Erlangen, 91058, Germany
| | - Ezzat El-Sherif
- Department of Biology, Friedrich-Alexander-Universität Erlangen-Nürnberg, Staudtstrasse 5, Erlangen, 91058, Germany.
| |
Collapse
|
8
|
|
9
|
Gupta S, Varennes J, Korswagen HC, Mugler A. Temporal precision of regulated gene expression. PLoS Comput Biol 2018; 14:e1006201. [PMID: 29879102 PMCID: PMC5991653 DOI: 10.1371/journal.pcbi.1006201] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2017] [Accepted: 05/14/2018] [Indexed: 11/18/2022] Open
Abstract
Important cellular processes such as migration, differentiation, and development often rely on precise timing. Yet, the molecular machinery that regulates timing is inherently noisy. How do cells achieve precise timing with noisy components? We investigate this question using a first-passage-time approach, for an event triggered by a molecule that crosses an abundance threshold and that is regulated by either an accumulating activator or a diminishing repressor. We find that either activation or repression outperforms an unregulated strategy. The optimal regulation corresponds to a nonlinear increase in the amount of the target molecule over time, arises from a tradeoff between minimizing the timing noise of the regulator and that of the target molecule itself, and is robust to additional effects such as bursts and cell division. Our results are in quantitative agreement with the nonlinear increase and low noise of mig-1 gene expression in migrating neuroblast cells during Caenorhabditis elegans development. These findings suggest that dynamic regulation may be a simple and powerful strategy for precise cellular timing. Cells control important processes with precise timing, even though their underlying molecular machinery is inherently imprecise. In the case of Caenorhabditis elegans development, migrating neuroblast cells produce a molecule until a certain abundance is reached, at which time the cells stop moving. Precise timing of this event is critical to C. elegans development, and here we investigate how it can be achieved. Specifically, we investigate regulation of the molecule production by either an accumulating activator or a diminishing repressor. Our results are consistent with the nonlinear increase and low noise of gene expression observed in the C. elegans cells.
Collapse
Affiliation(s)
- Shivam Gupta
- Department of Physics and Astronomy, Purdue University, West Lafayette, Indiana, United States of America
| | - Julien Varennes
- Department of Physics and Astronomy, Purdue University, West Lafayette, Indiana, United States of America
| | - Hendrik C. Korswagen
- Hubrecht Institute, Royal Netherlands Academy of Arts and Sciences and University Medical Center Utrecht, Utrecht, the Netherlands
| | - Andrew Mugler
- Department of Physics and Astronomy, Purdue University, West Lafayette, Indiana, United States of America
- * E-mail:
| |
Collapse
|
10
|
A damped oscillator imposes temporal order on posterior gap gene expression in Drosophila. PLoS Biol 2018; 16:e2003174. [PMID: 29451884 PMCID: PMC5832388 DOI: 10.1371/journal.pbio.2003174] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2017] [Revised: 03/01/2018] [Accepted: 01/31/2018] [Indexed: 12/21/2022] Open
Abstract
Insects determine their body segments in two different ways. Short-germband insects, such as the flour beetle Tribolium castaneum, use a molecular clock to establish segments sequentially. In contrast, long-germband insects, such as the vinegar fly Drosophila melanogaster, determine all segments simultaneously through a hierarchical cascade of gene regulation. Gap genes constitute the first layer of the Drosophila segmentation gene hierarchy, downstream of maternal gradients such as that of Caudal (Cad). We use data-driven mathematical modelling and phase space analysis to show that shifting gap domains in the posterior half of the Drosophila embryo are an emergent property of a robust damped oscillator mechanism, suggesting that the regulatory dynamics underlying long- and short-germband segmentation are much more similar than previously thought. In Tribolium, Cad has been proposed to modulate the frequency of the segmentation oscillator. Surprisingly, our simulations and experiments show that the shift rate of posterior gap domains is independent of maternal Cad levels in Drosophila. Our results suggest a novel evolutionary scenario for the short- to long-germband transition and help explain why this transition occurred convergently multiple times during the radiation of the holometabolan insects. Different insect species exhibit one of two distinct modes of determining their body segments (known as segmentation) during development: they either use a molecular oscillator to position segments sequentially, or they generate segments simultaneously through a hierarchical gene-regulatory cascade. The sequential mode is ancestral, while the simultaneous mode has been derived from it independently several times during evolution. In this paper, we present evidence suggesting that simultaneous segmentation also involves an oscillator in the posterior end of the embryo of the vinegar fly, Drosophila melanogaster. This surprising result indicates that both modes of segment determination are much more similar than previously thought. Such similarity provides an important step towards our understanding of the frequent evolutionary transitions observed between sequential and simultaneous segmentation.
Collapse
|
11
|
Speed regulation of genetic cascades allows for evolvability in the body plan specification of insects. Proc Natl Acad Sci U S A 2017; 114:E8646-E8655. [PMID: 28973882 DOI: 10.1073/pnas.1702478114] [Citation(s) in RCA: 39] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
During the anterior-posterior fate specification of insects, anterior fates arise in a nonelongating tissue (called the "blastoderm"), and posterior fates arise in an elongating tissue (called the "germband"). However, insects differ widely in the extent to which anterior-posterior fates are specified in the blastoderm versus the germband. Here we present a model in which patterning in both the blastoderm and germband of the beetle Tribolium castaneum is based on the same flexible mechanism: a gradient that modulates the speed of a genetic cascade of gap genes, resulting in the induction of sequential kinematic waves of gap gene expression. The mechanism is flexible and capable of patterning both elongating and nonelongating tissues, and hence converting blastodermal to germband fates and vice versa. Using RNAi perturbations, we found that blastodermal fates could be shifted to the germband, and germband fates could be generated in a blastoderm-like morphology. We also suggest a molecular mechanism underlying our model, in which gradient levels regulate the switch between two enhancers: One enhancer is responsible for sequential gene activation, and the other is responsible for freezing temporal rhythms into spatial patterns. This model is consistent with findings in Drosophila melanogaster, where gap genes were found to be regulated by two nonredundant "shadow" enhancers.
Collapse
|
12
|
|
13
|
Perez-Carrasco R, Guerrero P, Briscoe J, Page KM. Intrinsic Noise Profoundly Alters the Dynamics and Steady State of Morphogen-Controlled Bistable Genetic Switches. PLoS Comput Biol 2016; 12:e1005154. [PMID: 27768683 PMCID: PMC5074595 DOI: 10.1371/journal.pcbi.1005154] [Citation(s) in RCA: 43] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2016] [Accepted: 09/19/2016] [Indexed: 01/08/2023] Open
Abstract
During tissue development, patterns of gene expression determine the spatial arrangement of cell types. In many cases, gradients of secreted signalling molecules—morphogens—guide this process by controlling downstream transcriptional networks. A mechanism commonly used in these networks to convert the continuous information provided by the gradient into discrete transitions between adjacent cell types is the genetic toggle switch, composed of cross-repressing transcriptional determinants. Previous analyses have emphasised the steady state output of these mechanisms. Here, we explore the dynamics of the toggle switch and use exact numerical simulations of the kinetic reactions, the corresponding Chemical Langevin Equation, and Minimum Action Path theory to establish a framework for studying the effect of gene expression noise on patterning time and boundary position. This provides insight into the time scale, gene expression trajectories and directionality of stochastic switching events between cell states. Taking gene expression noise into account predicts that the final boundary position of a morphogen-induced toggle switch, although robust to changes in the details of the noise, is distinct from that of the deterministic system. Moreover, the dramatic increase in patterning time close to the boundary predicted from the deterministic case is substantially reduced. The resulting stochastic switching introduces differences in patterning time along the morphogen gradient that result in a patterning wave propagating away from the morphogen source with a velocity determined by the intrinsic noise. The wave sharpens and slows as it advances and may never reach steady state in a biologically relevant time. This could explain experimentally observed dynamics of pattern formation. Together the analysis reveals the importance of dynamical transients for understanding morphogen-driven transcriptional networks and indicates that gene expression noise can qualitatively alter developmental patterning. The bistable switch, a common regulatory sub-network, is found in many biological processes. It consists of cross-repressing components that generate a switch-like transition between two possible states. In developing tissues, bistable switches, created by cross-repressing transcriptional determinants, are often controlled by gradients of secreted signalling molecules—morphogens. These provide a mechanism to convert a morphogen gradient into stripes of gene expression that determine the arrangement of distinct cell types. Here we use mathematical models to analyse the temporal response of such a system. We find that the behaviour is highly dependent on the intrinsic fluctuations that result from the stochastic nature of gene expression. This noise has a marked effect on both patterning time and the location of the stripe boundary. One of the techniques we use, Minimum Action Path theory, identifies key features of the switch without computationally expensive calculations. The results reveal a noise driven switching wave that propels the stripe boundary away from the morphogen source to eventually settle, at steady state, further from the morphogen source than in the deterministic description. Together the analysis highlights the importance dynamics in patterning and demonstrates a set of mathematical tools for studying this problem.
Collapse
Affiliation(s)
- Ruben Perez-Carrasco
- Department of Mathematics, University College London, Gower Street, London WC1E 6BT, UK
- * E-mail:
| | - Pilar Guerrero
- Department of Mathematics, University College London, Gower Street, London WC1E 6BT, UK
| | - James Briscoe
- The Francis Crick Institute, 1 Midland Road, London, NW1 1AT, UK
| | - Karen M. Page
- Department of Mathematics, University College London, Gower Street, London WC1E 6BT, UK
| |
Collapse
|
14
|
Beaupeux M, François P. Positional information from oscillatory phase shifts : insights from in silico evolution. Phys Biol 2016; 13:036009. [PMID: 27346171 DOI: 10.1088/1478-3975/13/3/036009] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Complex cellular decisions are based on temporal dynamics of pathways, including genetic oscillators. In development, recent works on vertebrae formation have suggested that relative phase of genetic oscillators encode positional information, including differentiation front defining vertebrae positions. Precise mechanisms for this are still unknown. Here, we use computational evolution to find gene network topologies that can compute the phase difference between oscillators and convert it into a decoder morphogen concentration. Two types of networks are discovered, based on symmetry properties of the decoder gene. So called asymmetric networks are studied, and two submodules are identified converting phase information into an amplitude variable. Those networks naturally display a 'shock' for a well defined phase difference, that can be used to define a wavefront of differentiation. We show how implementation of these ideas reproduce experimental features of vertebrate segmentation.
Collapse
Affiliation(s)
- M Beaupeux
- Ernest Rutherford Physics Building, McGill University, H3A2T8 Montreal QC, Canada
| | | |
Collapse
|