1
|
Sawicki J, Berner R, Loos SAM, Anvari M, Bader R, Barfuss W, Botta N, Brede N, Franović I, Gauthier DJ, Goldt S, Hajizadeh A, Hövel P, Karin O, Lorenz-Spreen P, Miehl C, Mölter J, Olmi S, Schöll E, Seif A, Tass PA, Volpe G, Yanchuk S, Kurths J. Perspectives on adaptive dynamical systems. CHAOS (WOODBURY, N.Y.) 2023; 33:071501. [PMID: 37486668 DOI: 10.1063/5.0147231] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/20/2023] [Accepted: 05/24/2023] [Indexed: 07/25/2023]
Abstract
Adaptivity is a dynamical feature that is omnipresent in nature, socio-economics, and technology. For example, adaptive couplings appear in various real-world systems, such as the power grid, social, and neural networks, and they form the backbone of closed-loop control strategies and machine learning algorithms. In this article, we provide an interdisciplinary perspective on adaptive systems. We reflect on the notion and terminology of adaptivity in different disciplines and discuss which role adaptivity plays for various fields. We highlight common open challenges and give perspectives on future research directions, looking to inspire interdisciplinary approaches.
Collapse
Affiliation(s)
- Jakub Sawicki
- Potsdam Institute for Climate Impact Research, Telegrafenberg, 14473 Potsdam, Germany
- Akademie Basel, Fachhochschule Nordwestschweiz FHNW, Leonhardsstrasse 6, 4009 Basel, Switzerland
| | - Rico Berner
- Department of Physics, Humboldt-Universität zu Berlin, Newtonstraße 15, 12489 Berlin, Germany
| | - Sarah A M Loos
- DAMTP, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, United Kingdom
| | - Mehrnaz Anvari
- Potsdam Institute for Climate Impact Research, Telegrafenberg, 14473 Potsdam, Germany
- Fraunhofer Institute for Algorithms and Scientific Computing, Schloss Birlinghoven, 53757 Sankt-Augustin, Germany
| | - Rolf Bader
- Institute of Systematic Musicology, University of Hamburg, Hamburg, Germany
| | - Wolfram Barfuss
- Transdisciplinary Research Area: Sustainable Futures, University of Bonn, 53113 Bonn, Germany
- Center for Development Research (ZEF), University of Bonn, 53113 Bonn, Germany
| | - Nicola Botta
- Potsdam Institute for Climate Impact Research, Telegrafenberg, 14473 Potsdam, Germany
- Department of Computer Science and Engineering, Chalmers University of Technology, 412 96 Göteborg, Sweden
| | - Nuria Brede
- Potsdam Institute for Climate Impact Research, Telegrafenberg, 14473 Potsdam, Germany
- Department of Computer Science, University of Potsdam, An der Bahn 2, 14476 Potsdam, Germany
| | - Igor Franović
- Scientific Computing Laboratory, Center for the Study of Complex Systems, Institute of Physics Belgrade, University of Belgrade, Pregrevica 118, 11080 Belgrade, Serbia
| | - Daniel J Gauthier
- Potsdam Institute for Climate Impact Research, Telegrafenberg, 14473 Potsdam, Germany
| | - Sebastian Goldt
- Department of Physics, International School of Advanced Studies (SISSA), Trieste, Italy
| | - Aida Hajizadeh
- Research Group Comparative Neuroscience, Leibniz Institute for Neurobiology, Magdeburg, Germany
| | - Philipp Hövel
- Potsdam Institute for Climate Impact Research, Telegrafenberg, 14473 Potsdam, Germany
| | - Omer Karin
- Department of Mathematics, Imperial College London, London SW7 2AZ, United Kingdom
| | - Philipp Lorenz-Spreen
- Center for Adaptive Rationality, Max Planck Institute for Human Development, Lentzeallee 94, 14195 Berlin, Germany
| | - Christoph Miehl
- Akademie Basel, Fachhochschule Nordwestschweiz FHNW, Leonhardsstrasse 6, 4009 Basel, Switzerland
| | - Jan Mölter
- Department of Mathematics, School of Computation, Information and Technology, Technical University of Munich, Boltzmannstraße 3, 85748 Garching bei München, Germany
| | - Simona Olmi
- Akademie Basel, Fachhochschule Nordwestschweiz FHNW, Leonhardsstrasse 6, 4009 Basel, Switzerland
| | - Eckehard Schöll
- Potsdam Institute for Climate Impact Research, Telegrafenberg, 14473 Potsdam, Germany
- Akademie Basel, Fachhochschule Nordwestschweiz FHNW, Leonhardsstrasse 6, 4009 Basel, Switzerland
| | - Alireza Seif
- Pritzker School of Molecular Engineering, The University of Chicago, Chicago, Illinois 60637, USA
| | - Peter A Tass
- Department of Neurosurgery, Stanford University School of Medicine, Stanford, California 94304, USA
| | - Giovanni Volpe
- Department of Physics, University of Gothenburg, Gothenburg, Sweden
| | - Serhiy Yanchuk
- Potsdam Institute for Climate Impact Research, Telegrafenberg, 14473 Potsdam, Germany
- Department of Physics, Humboldt-Universität zu Berlin, Newtonstraße 15, 12489 Berlin, Germany
| | - Jürgen Kurths
- Potsdam Institute for Climate Impact Research, Telegrafenberg, 14473 Potsdam, Germany
- Department of Physics, Humboldt-Universität zu Berlin, Newtonstraße 15, 12489 Berlin, Germany
| |
Collapse
|
2
|
Ristič D, Gosak M. Interlayer Connectivity Affects the Coherence Resonance and Population Activity Patterns in Two-Layered Networks of Excitatory and Inhibitory Neurons. Front Comput Neurosci 2022; 16:885720. [PMID: 35521427 PMCID: PMC9062746 DOI: 10.3389/fncom.2022.885720] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2022] [Accepted: 03/24/2022] [Indexed: 11/13/2022] Open
Abstract
The firing patterns of neuronal populations often exhibit emergent collective oscillations, which can display substantial regularity even though the dynamics of individual elements is very stochastic. One of the many phenomena that is often studied in this context is coherence resonance, where additional noise leads to improved regularity of spiking activity in neurons. In this work, we investigate how the coherence resonance phenomenon manifests itself in populations of excitatory and inhibitory neurons. In our simulations, we use the coupled FitzHugh-Nagumo oscillators in the excitable regime and in the presence of neuronal noise. Formally, our model is based on the concept of a two-layered network, where one layer contains inhibitory neurons, the other excitatory neurons, and the interlayer connections represent heterotypic interactions. The neuronal activity is simulated in realistic coupling schemes in which neurons within each layer are connected with undirected connections, whereas neurons of different types are connected with directed interlayer connections. In this setting, we investigate how different neurophysiological determinants affect the coherence resonance. Specifically, we focus on the proportion of inhibitory neurons, the proportion of excitatory interlayer axons, and the architecture of interlayer connections between inhibitory and excitatory neurons. Our results reveal that the regularity of simulated neural activity can be increased by a stronger damping of the excitatory layer. This can be accomplished with a higher proportion of inhibitory neurons, a higher fraction of inhibitory interlayer axons, a stronger coupling between inhibitory axons, or by a heterogeneous configuration of interlayer connections. Our approach of modeling multilayered neuronal networks in combination with stochastic dynamics offers a novel perspective on how the neural architecture can affect neural information processing and provide possible applications in designing networks of artificial neural circuits to optimize their function via noise-induced phenomena.
Collapse
Affiliation(s)
- David Ristič
- Faculty of Natural Sciences and Mathematics, University of Maribor, Maribor, Slovenia
| | - Marko Gosak
- Faculty of Natural Sciences and Mathematics, University of Maribor, Maribor, Slovenia
- Faculty of Medicine, University of Maribor, Maribor, Slovenia
| |
Collapse
|
3
|
Franović I, Eydam S, Yanchuk S, Berner R. Collective Activity Bursting in a Population of Excitable Units Adaptively Coupled to a Pool of Resources. FRONTIERS IN NETWORK PHYSIOLOGY 2022; 2:841829. [PMID: 36926089 PMCID: PMC10013072 DOI: 10.3389/fnetp.2022.841829] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/22/2021] [Accepted: 02/16/2022] [Indexed: 06/18/2023]
Abstract
We study the collective dynamics in a population of excitable units (neurons) adaptively interacting with a pool of resources. The resource pool is influenced by the average activity of the population, whereas the feedback from the resources to the population is comprised of components acting homogeneously or inhomogeneously on individual units of the population. Moreover, the resource pool dynamics is assumed to be slow and has an oscillatory degree of freedom. We show that the feedback loop between the population and the resources can give rise to collective activity bursting in the population. To explain the mechanisms behind this emergent phenomenon, we combine the Ott-Antonsen reduction for the collective dynamics of the population and singular perturbation theory to obtain a reduced system describing the interaction between the population mean field and the resources.
Collapse
Affiliation(s)
- Igor Franović
- Scientific Computing Laboratory, Center for the Study of Complex Systems, Institute of Physics Belgrade, University of Belgrade, Belgrade, Serbia
| | - Sebastian Eydam
- Neural Circuits and Computations Unit, RIKEN Center for Brain Science, Wako, Japan
| | - Serhiy Yanchuk
- Institut für Mathematik, Technische Universität Berlin, Berlin, Germany
- Potsdam Institute for Climate Impact Research, Potsdam, Germany
- Institut für Mathematik, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Rico Berner
- Institut für Physik, Humboldt-Universität zu Berlin, Berlin, Germany
- Institut für Theoretische Physik, Technische Universität Berlin, Berlin, Germany
| |
Collapse
|
4
|
Stožer A, Šterk M, Paradiž Leitgeb E, Markovič R, Skelin Klemen M, Ellis CE, Križančić Bombek L, Dolenšek J, MacDonald PE, Gosak M. From Isles of Königsberg to Islets of Langerhans: Examining the Function of the Endocrine Pancreas Through Network Science. Front Endocrinol (Lausanne) 2022; 13:922640. [PMID: 35784543 PMCID: PMC9240343 DOI: 10.3389/fendo.2022.922640] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/18/2022] [Accepted: 05/16/2022] [Indexed: 12/12/2022] Open
Abstract
Islets of Langerhans are multicellular microorgans located in the pancreas that play a central role in whole-body energy homeostasis. Through secretion of insulin and other hormones they regulate postprandial storage and interprandial usage of energy-rich nutrients. In these clusters of hormone-secreting endocrine cells, intricate cell-cell communication is essential for proper function. Electrical coupling between the insulin-secreting beta cells through gap junctions composed of connexin36 is particularly important, as it provides the required, most important, basis for coordinated responses of the beta cell population. The increasing evidence that gap-junctional communication and its modulation are vital to well-regulated secretion of insulin has stimulated immense interest in how subpopulations of heterogeneous beta cells are functionally arranged throughout the islets and how they mediate intercellular signals. In the last decade, several novel techniques have been proposed to assess cooperation between cells in islets, including the prosperous combination of multicellular imaging and network science. In the present contribution, we review recent advances related to the application of complex network approaches to uncover the functional connectivity patterns among cells within the islets. We first provide an accessible introduction to the basic principles of network theory, enumerating the measures characterizing the intercellular interactions and quantifying the functional integration and segregation of a multicellular system. Then we describe methodological approaches to construct functional beta cell networks, point out possible pitfalls, and specify the functional implications of beta cell network examinations. We continue by highlighting the recent findings obtained through advanced multicellular imaging techniques supported by network-based analyses, giving special emphasis to the current developments in both mouse and human islets, as well as outlining challenges offered by the multilayer network formalism in exploring the collective activity of islet cell populations. Finally, we emphasize that the combination of these imaging techniques and network-based analyses does not only represent an innovative concept that can be used to describe and interpret the physiology of islets, but also provides fertile ground for delineating normal from pathological function and for quantifying the changes in islet communication networks associated with the development of diabetes mellitus.
Collapse
Affiliation(s)
- Andraž Stožer
- Institute of Physiology, Faculty of Medicine, University of Maribor, Maribor, Slovenia
| | - Marko Šterk
- Institute of Physiology, Faculty of Medicine, University of Maribor, Maribor, Slovenia
- Department of Physics, Faculty of Natural Sciences and Mathematics, University of Maribor, Maribor, Slovenia
| | - Eva Paradiž Leitgeb
- Institute of Physiology, Faculty of Medicine, University of Maribor, Maribor, Slovenia
| | - Rene Markovič
- Department of Physics, Faculty of Natural Sciences and Mathematics, University of Maribor, Maribor, Slovenia
- Institute of Mathematics and Physics, Faculty of Electrical Engineering and Computer Science, University of Maribor, Maribor, Slovenia
| | - Maša Skelin Klemen
- Institute of Physiology, Faculty of Medicine, University of Maribor, Maribor, Slovenia
| | - Cara E. Ellis
- Department of Pharmacology and Alberta Diabetes Institute, University of Alberta, Edmonton, AB, Canada
| | | | - Jurij Dolenšek
- Institute of Physiology, Faculty of Medicine, University of Maribor, Maribor, Slovenia
- Department of Physics, Faculty of Natural Sciences and Mathematics, University of Maribor, Maribor, Slovenia
| | - Patrick E. MacDonald
- Department of Pharmacology and Alberta Diabetes Institute, University of Alberta, Edmonton, AB, Canada
| | - Marko Gosak
- Institute of Physiology, Faculty of Medicine, University of Maribor, Maribor, Slovenia
- Department of Physics, Faculty of Natural Sciences and Mathematics, University of Maribor, Maribor, Slovenia
- *Correspondence: Marko Gosak,
| |
Collapse
|
5
|
Gast R, Schmidt H, Knösche TR. A Mean-Field Description of Bursting Dynamics in Spiking Neural Networks with Short-Term Adaptation. Neural Comput 2020; 32:1615-1634. [PMID: 32687770 DOI: 10.1162/neco_a_01300] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Bursting plays an important role in neural communication. At the population level, macroscopic bursting has been identified in populations of neurons that do not express intrinsic bursting mechanisms. For the analysis of phase transitions between bursting and non-bursting states, mean-field descriptions of macroscopic bursting behavior are a valuable tool. In this article, we derive mean-field descriptions of populations of spiking neurons and examine whether states of collective bursting behavior can arise from short-term adaptation mechanisms. Specifically, we consider synaptic depression and spike-frequency adaptation in networks of quadratic integrate-and-fire neurons. Analyzing the mean-field model via bifurcation analysis, we find that bursting behavior emerges for both types of short-term adaptation. This bursting behavior can coexist with steady-state behavior, providing a bistable regime that allows for transient switches between synchronized and nonsynchronized states of population dynamics. For all of these findings, we demonstrate a close correspondence between the spiking neural network and the mean-field model. Although the mean-field model has been derived under the assumptions of an infinite population size and all-to-all coupling inside the population, we show that this correspondence holds even for small, sparsely coupled networks. In summary, we provide mechanistic descriptions of phase transitions between bursting and steady-state population dynamics, which play important roles in both healthy neural communication and neurological disorders.
Collapse
Affiliation(s)
- Richard Gast
- Max Planck Institute for Human Cognitive and Brain Sciences, 04103 Leipzig, Germany
| | - Helmut Schmidt
- Max Planck Institute for Human Cognitive and Brain Sciences, 04103 Leipzig, Germany
| | - Thomas R Knösche
- Max Planck Institute for Human Cognitive and Brain Sciences, 04103 Leipzig, Germany, and Institute for Biomedical Engineering and Informatics, TU 98693 Ilmenau, Germany
| |
Collapse
|
6
|
Berner R, Sawicki J, Schöll E. Birth and Stabilization of Phase Clusters by Multiplexing of Adaptive Networks. PHYSICAL REVIEW LETTERS 2020; 124:088301. [PMID: 32167358 DOI: 10.1103/physrevlett.124.088301] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/27/2019] [Revised: 12/05/2019] [Accepted: 01/16/2020] [Indexed: 06/10/2023]
Abstract
We propose a concept to generate and stabilize diverse partial synchronization patterns (phase clusters) in adaptive networks which are widespread in neuroscience and social sciences, as well as biology, engineering, and other disciplines. We show by theoretical analysis and computer simulations that multiplexing in a multilayer network with symmetry can induce various stable phase cluster states in a situation where they are not stable or do not even exist in the single layer. Further, we develop a method for the analysis of Laplacian matrices of multiplex networks which allows for insight into the spectral structure of these networks enabling a reduction to the stability problem of single layers. We employ the multiplex decomposition to provide analytic results for the stability of the multilayer patterns. As local dynamics we use the paradigmatic Kuramoto phase oscillator, which is a simple generic model and has been successfully applied in the modeling of synchronization phenomena in a wide range of natural and technological systems.
Collapse
Affiliation(s)
- Rico Berner
- Institut für Theoretische Physik, Technische Universität Berlin, Hardenbergstrasse 36, 10623 Berlin, Germany
- Institut für Mathematik, Technische Universität Berlin, Hardenbergstrasse 36, 10623 Berlin, Germany
| | - Jakub Sawicki
- Institut für Theoretische Physik, Technische Universität Berlin, Hardenbergstrasse 36, 10623 Berlin, Germany
| | - Eckehard Schöll
- Institut für Theoretische Physik, Technische Universität Berlin, Hardenbergstrasse 36, 10623 Berlin, Germany
| |
Collapse
|
7
|
Virkar YS, Restrepo JG, Shew WL, Ott E. Dynamic regulation of resource transport induces criticality in interdependent networks of excitable units. Phys Rev E 2020; 101:022303. [PMID: 32168577 DOI: 10.1103/physreve.101.022303] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2018] [Accepted: 12/24/2019] [Indexed: 11/06/2022]
Abstract
Various functions of a network of excitable units can be enhanced if the network is in the "critical regime," where excitations are, on average, neither damped nor amplified. An important question is how can such networks self-organize to operate in the critical regime. Previously, it was shown that regulation via resource transport on a secondary network can robustly maintain the primary network dynamics in a balanced state where activity doesn't grow or decay. Here we show that this internetwork regulation process robustly produces a power-law distribution of activity avalanches, as observed in experiments, over ranges of model parameters spanning orders of magnitude. We also show that the resource transport over the secondary network protects the system against the destabilizing effect of local variations in parameters and heterogeneity in network structure. For homogeneous networks, we derive a reduced three-dimensional map which reproduces the behavior of the full system.
Collapse
Affiliation(s)
- Yogesh S Virkar
- Department of Computer Science, University of Colorado at Boulder, Boulder, Colorado 80309, USA
| | - Juan G Restrepo
- Department of Applied Mathematics, University of Colorado at Boulder, Boulder, Colorado 80309, USA
| | - Woodrow L Shew
- Department of Physics, University of Arkansas, Fayetteville, Arkansas 72701, USA
| | - Edward Ott
- Departments of Electrical and Computer Engineering and of Physics, University of Maryland, College Park, Maryland 20742, USA
| |
Collapse
|
8
|
Khoshkhou M, Montakhab A. Spike-Timing-Dependent Plasticity With Axonal Delay Tunes Networks of Izhikevich Neurons to the Edge of Synchronization Transition With Scale-Free Avalanches. Front Syst Neurosci 2019; 13:73. [PMID: 31866836 PMCID: PMC6904334 DOI: 10.3389/fnsys.2019.00073] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2019] [Accepted: 11/19/2019] [Indexed: 11/13/2022] Open
Abstract
Critical brain hypothesis has been intensively studied both in experimental and theoretical neuroscience over the past two decades. However, some important questions still remain: (i) What is the critical point the brain operates at? (ii) What is the regulatory mechanism that brings about and maintains such a critical state? (iii) The critical state is characterized by scale-invariant behavior which is seemingly at odds with definitive brain oscillations? In this work we consider a biologically motivated model of Izhikevich neuronal network with chemical synapses interacting via spike-timing-dependent plasticity (STDP) as well as axonal time delay. Under generic and physiologically relevant conditions we show that the system is organized and maintained around a synchronization transition point as opposed to an activity transition point associated with an absorbing state phase transition. However, such a state exhibits experimentally relevant signs of critical dynamics including scale-free avalanches with finite-size scaling as well as critical branching ratios. While the system displays stochastic oscillations with highly correlated fluctuations, it also displays dominant frequency modes seen as sharp peaks in the power spectrum. The role of STDP as well as time delay is crucial in achieving and maintaining such critical dynamics, while the role of inhibition is not as crucial. In this way we provide possible answers to all three questions posed above. We also show that one can achieve supercritical or subcritical dynamics if one changes the average time delay associated with axonal conduction.
Collapse
Affiliation(s)
- Mahsa Khoshkhou
- Department of Physics, College of Sciences, Shiraz University, Shiraz, Iran
| | - Afshin Montakhab
- Department of Physics, College of Sciences, Shiraz University, Shiraz, Iran
| |
Collapse
|
9
|
Zeng G, Huang X, Jiang T, Yu S. Short-term synaptic plasticity expands the operational range of long-term synaptic changes in neural networks. Neural Netw 2019; 118:140-147. [DOI: 10.1016/j.neunet.2019.06.002] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2018] [Revised: 04/23/2019] [Accepted: 06/03/2019] [Indexed: 02/04/2023]
|
10
|
Cocchi L, Gollo LL, Zalesky A, Breakspear M. Criticality in the brain: A synthesis of neurobiology, models and cognition. Prog Neurobiol 2017; 158:132-152. [PMID: 28734836 DOI: 10.1016/j.pneurobio.2017.07.002] [Citation(s) in RCA: 226] [Impact Index Per Article: 32.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2016] [Revised: 06/15/2017] [Accepted: 07/13/2017] [Indexed: 11/26/2022]
Abstract
Cognitive function requires the coordination of neural activity across many scales, from neurons and circuits to large-scale networks. As such, it is unlikely that an explanatory framework focused upon any single scale will yield a comprehensive theory of brain activity and cognitive function. Modelling and analysis methods for neuroscience should aim to accommodate multiscale phenomena. Emerging research now suggests that multi-scale processes in the brain arise from so-called critical phenomena that occur very broadly in the natural world. Criticality arises in complex systems perched between order and disorder, and is marked by fluctuations that do not have any privileged spatial or temporal scale. We review the core nature of criticality, the evidence supporting its role in neural systems and its explanatory potential in brain health and disease.
Collapse
Affiliation(s)
- Luca Cocchi
- QIMR Berghofer Medical Research Institute, Brisbane, Australia.
| | | | - Andrew Zalesky
- Melbourne Neuropsychiatry Centre, The University of Melbourne, Melbourne, Australia
| | - Michael Breakspear
- QIMR Berghofer Medical Research Institute, Brisbane, Australia; Metro North Mental Health Service, Brisbane, Australia
| |
Collapse
|