1
|
Stone E, Coombs S, Landguth E. Archetypal analysis of COVID-19 in Montana, USA, March 13, 2020 to April 26, 2022. PLoS One 2024; 19:e0283265. [PMID: 38170725 PMCID: PMC10763954 DOI: 10.1371/journal.pone.0283265] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2023] [Accepted: 10/19/2023] [Indexed: 01/05/2024] Open
Abstract
Infectious disease data can often involve complex spatial patterns intermixed with temporal trends. Archetypal Analysis is a method to mine complex spatio-temporal data, and can be used to discover the dynamics of spatial patterns. The application of Archetypal Analysis to epidemiological data is relatively new, and here we present one of the first applications on COVID-19 data from March 13, 2020 to April 26, 2022, for the counties of Montana, USA. We present three views of the data set decomposed with Archetypal Analysis. First, we evaluate the entire 56 county data set. Second, we use a mutual information calculation to remove counties whose dynamics are mainly independent from the other counties, reducing the set to 17 counties. Finally, we analyze the top ten counties in terms of population size to focus on the dynamics in the large cities in the state. For each data set, we analyze four significant disease outbreaks across Montana. Archetypal Analysis uncovers distinct spatial patterns for each outbreak and demonstrates that each has a unique trajectory across the state.
Collapse
Affiliation(s)
- Emily Stone
- Dept. of Mathematical Science, University of Montana, Missoula, MT, United States of America
| | - Sebastian Coombs
- Dept. of Mathematical Science, University of Montana, Missoula, MT, United States of America
| | - Erin Landguth
- Center for Population Health Research, School of Public and Community Health Sciences, University of Montana, Missoula, MT, United States of America
| |
Collapse
|
2
|
Stone E, Coombs S, Landguth E. Archetypal analysis of COVID-19 in Montana, USA, March 13, 2020 to April 26, 2022. MEDRXIV : THE PREPRINT SERVER FOR HEALTH SCIENCES 2023:2023.03.06.23286886. [PMID: 36945386 PMCID: PMC10029034 DOI: 10.1101/2023.03.06.23286886] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/12/2023]
Abstract
Given the potential consequences of infectious diseases, it is important to understand how broad scale incidence variability influences the probability of localized outbreaks. Often, these infectious disease data can involve complex spatial patterns intermixed with temporal trends. Archetypal Analysis is a method to mine complex spatiotemporal epidemiological data, and can be used to discover the dynamics of spatial patterns. The application of Archetypal Analysis to epistemological data is relatively new, and here we present one of the first applications using COVID-19 data from March 13, 2020 to April 26, 2022, in the counties of Montana, USA. We present three views of the data set with Archetypal Analysis. First, we evaluate the entire 56 county data set. Second, we compute mutual information of the 56 counties' time series to remove counties whose dynamics are mainly independent from most of the other counties. We choose the top 17 counties ranked in terms of increasing total mutual information. Finally, to compare how population size might influence results, we conducted an analysis with 10 of the largest counties. Using the Archetypal Analysis results, we analyze the disease outbreaks across Montana, comparing and contrasting the three different cases and showing how certain counties can be found in distinct sets of archetypes. Using the reconstruction time series, we show how each outbreak had a unique trajectory across the state in terms of the archetypes.
Collapse
Affiliation(s)
- Emily Stone
- Dept. of Mathematical Science, University of Montana, 32 Campus Drive, Missoula, MT 59812, USA
| | - Sebastian Coombs
- Dept. of Mathematical Science, University of Montana, 32 Campus Drive, Missoula, MT 59812, USA
| | - Erin Landguth
- Center for Population Health Research, School of Public and Community Health Sciences, University of Montana, 32 Campus Drive, Missoula, MT 59812, USA
| |
Collapse
|
3
|
Olsen AS, Høegh RMT, Hinrich JL, Madsen KH, Mørup M. Combining electro- and magnetoencephalography data using directional archetypal analysis. Front Neurosci 2022; 16:911034. [PMID: 35968377 PMCID: PMC9374169 DOI: 10.3389/fnins.2022.911034] [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: 04/01/2022] [Accepted: 07/11/2022] [Indexed: 11/20/2022] Open
Abstract
Metastable microstates in electro- and magnetoencephalographic (EEG and MEG) measurements are usually determined using modified k-means accounting for polarity invariant states. However, hard state assignment approaches assume that the brain traverses microstates in a discrete rather than continuous fashion. We present multimodal, multisubject directional archetypal analysis as a scale and polarity invariant extension to archetypal analysis using a loss function based on the Watson distribution. With this method, EEG/MEG microstates are modeled using subject- and modality-specific archetypes that are representative, distinct topographic maps between which the brain continuously traverses. Archetypes are specified as convex combinations of unit norm input data based on a shared generator matrix, thus assuming that the timing of neural responses to stimuli is consistent across subjects and modalities. The input data is reconstructed as convex combinations of archetypes using a subject- and modality-specific continuous archetypal mixing matrix. We showcase the model on synthetic data and an openly available face perception event-related potential data set with concurrently recorded EEG and MEG. In synthetic and unimodal experiments, we compare our model to conventional Euclidean multisubject archetypal analysis. We also contrast our model to a directional clustering model with discrete state assignments to highlight the advantages of modeling state trajectories rather than hard assignments. We find that our approach successfully models scale and polarity invariant data, such as microstates, accounting for intersubject and intermodal variability. The model is readily extendable to other modalities ensuring component correspondence while elucidating spatiotemporal signal variability.
Collapse
Affiliation(s)
- Anders S. Olsen
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Lyngby, Denmark
| | - Rasmus M. T. Høegh
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Lyngby, Denmark
- WS Audiology, Lynge, Denmark
| | - Jesper L. Hinrich
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Lyngby, Denmark
| | - Kristoffer H. Madsen
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Lyngby, Denmark
- Danish Research Centre for Magnetic Resonance, Centre for Functional and Diagnostic Imaging and Research, Copenhagen University Hospital Amager and Hvidovre, Hvidovre, Denmark
| | - Morten Mørup
- Department of Applied Mathematics and Computer Science, Technical University of Denmark, Lyngby, Denmark
| |
Collapse
|
4
|
Halloran PF, Böhmig GA, Bromberg J, Einecke G, Eskandary FA, Gupta G, Myslak M, Viklicky O, Perkowska-Ptasinska A, Madill-Thomsen KS. Archetypal Analysis of Injury in Kidney Transplant Biopsies Identifies Two Classes of Early AKI. Front Med (Lausanne) 2022; 9:817324. [PMID: 35463013 PMCID: PMC9021747 DOI: 10.3389/fmed.2022.817324] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2021] [Accepted: 03/07/2022] [Indexed: 01/07/2023] Open
Abstract
All transplanted kidneys are subjected to some degree of injury as a result of the donation-implantation process and various post-transplant stresses such as rejection. Because transplants are frequently biopsied, they present an opportunity to explore the full spectrum of kidney response-to-wounding from all causes. Defining parenchymal damage in transplanted organs is important for clinical management because it determines function and survival. In this study, we classified the scenarios associated with parenchymal injury in genome-wide microarray results from 1,526 kidney transplant indication biopsies collected during the INTERCOMEX study. We defined injury groups by using archetypal analysis (AA) of scores for gene sets and classifiers previously identified in various injury states. Six groups and their characteristics were defined in this population: No injury, minor injury, two classes of acute kidney injury ("AKI," AKI1, and AKI2), chronic kidney disease (CKD), and CKD combined with AKI. We compared the two classes of AKI, namely, AKI1 and AKI2. AKI1 had a poor function and increased parenchymal dedifferentiation but minimal response-to-injury and inflammation, instead having increased expression of PARD3, a gene previously characterized as being related to epithelial polarity and adherens junctions. In contrast, AKI2 had a poor function and increased response-to-injury, significant inflammation, and increased macrophage activity. In random forest analysis, the most important predictors of function (estimated glomerular filtration rate) and graft loss were injury-based molecular scores, not rejection scores. AKI1 and AKI2 differed in 3-year graft survival, with better survival in the AKI2 group. Thus, injury archetype analysis of injury-induced gene expression shows new heterogeneity in kidney response-to-wounding, revealing AKI1, a class of early transplants with a poor function but minimal inflammation or response to injury, a deviant response characterized as PC3, and an increased risk of failure. Given the relationship between parenchymal injury and kidney survival, further characterization of the injury phenotypes in kidney transplants will be important for an improved understanding that could have implications for understanding native kidney diseases (ClinicalTrials.gov #NCT01299168).
Collapse
Affiliation(s)
- Philip F Halloran
- Alberta Transplant Applied Genomics Centre, Edmonton, AB, Canada.,Division of Nephrology and Transplant Immunology, Department of Medicine, University of Alberta, Edmonton, AB, Canada
| | - Georg A Böhmig
- Division of Nephrology and Dialysis, Department of Medicine III, Medical University of Vienna, Vienna, Austria
| | - Jonathan Bromberg
- Department of Surgery, University of Maryland, Baltimore, MD, United States
| | - Gunilla Einecke
- Department of Nephrology, Hannover Medical School, Hannover, Germany
| | - Farsad A Eskandary
- Division of Nephrology and Dialysis, Department of Medicine III, Medical University of Vienna, Vienna, Austria
| | - Gaurav Gupta
- Division of Nephrology, Virginia Commonwealth University, Richmond, VA, United States
| | - Marek Myslak
- Department of Clinical Interventions, Department of Nephrology and Kidney Transplantation Samodzielny Publiczny Wojewódzki Szpital Zespolony (SPWSZ) Hospital, Pomeranian Medical University, Szczecin, Poland
| | - Ondrej Viklicky
- Department of Nephrology and Transplant Center, Institute for Clinical and Experimental Medicine, Prague, Czechia
| | - Agnieszka Perkowska-Ptasinska
- Department of Transplantation Medicine, Nephrology and Internal Diseases, Medical University of Warsaw, Warsaw, Poland
| | | | | |
Collapse
|
5
|
|
6
|
Mokhtari EB, Landguth EL, Anderson S, Stone E. Decoding influenza outbreaks in a rural region of the USA with archetypal analysis. Spat Spatiotemporal Epidemiol 2021; 38:100437. [PMID: 34353529 DOI: 10.1016/j.sste.2021.100437] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 02/11/2021] [Revised: 06/23/2021] [Accepted: 06/23/2021] [Indexed: 10/21/2022]
Abstract
We present the first application of archetypal analysis for influenza data from 2010 to 2018 in Montana, USA. Using archetypes, we decompose the data into spatial and temporal components, allowing for a more informed analysis of spatial-temporal dynamic trends during an influenza season. Initially, we reduce the dimension of the set of counties by using a mutual information measure on the influenza time series to create a smaller, maximal mutual information network. Archetypal analysis then describes the relationship between influenza cases across counties and regions in Montana. Finally, we discuss the potential implications this analysis can have for infectious disease modeling, particularly where data is sparse and limited.
Collapse
Affiliation(s)
- Elham Bayat Mokhtari
- Vaccine Research Center, National Institute of Allergy and Infectious Disease, National Institutes of Health, Bethesda, MD 20892, USA
| | - Erin L Landguth
- Center for Population Health Research, School of Public and Community Health Sciences, University of Montana, 32 Campus Drive, Missoula, MT 59812, USA
| | - Stacey Anderson
- Communicable Disease Control and Prevention Bureau, Department of Health and Human Services, Helena, MT 59620, USA
| | - Emily Stone
- Department of Mathematical Sciences, University of Montana, 32 Campus Drive, Missoula, MT 59812, USA.
| |
Collapse
|
7
|
Cabero I, Epifanio I, Piérola A, Ballester A. Archetype analysis: A new subspace outlier detection approach. Knowl Based Syst 2021. [DOI: 10.1016/j.knosys.2021.106830] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
8
|
Mikami T, Iwasaki W. The flipping
t
‐ratio test: Phylogenetically informed assessment of the Pareto theory for phenotypic evolution. Methods Ecol Evol 2021. [DOI: 10.1111/2041-210x.13553] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Tomoyuki Mikami
- Department of Biological Sciences Graduate School of Science The University of Tokyo Tokyo Japan
| | - Wataru Iwasaki
- Department of Biological Sciences Graduate School of Science The University of Tokyo Tokyo Japan
- Department of Computational Biology and Medical Sciences Graduate School of Frontier Sciences The University of Tokyo Kashiwa Japan
- Atmosphere and Ocean Research Institute The University of Tokyo Kashiwa Japan
- Institute for Quantitative Biosciences The University of Tokyo Tokyo Japan
- Collaborative Research Institute for Innovative Microbiology The University of Tokyo Tokyo Japan
| |
Collapse
|
9
|
Cao Q, Yang N, Wang Y, Xu C, Zhang X, Fan K, Chen F, Liang H, Zhang Y, Deng X, Feng Y, Yang CG, Wu M, Bae T, Lan L. Mutation-induced remodeling of the BfmRS two-component system in Pseudomonas aeruginosa clinical isolates. Sci Signal 2020; 13:13/656/eaaz1529. [PMID: 33144518 DOI: 10.1126/scisignal.aaz1529] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Genetic mutations are a primary driving force behind the adaptive evolution of bacterial pathogens. Multiple clinical isolates of Pseudomonas aeruginosa, an important human pathogen, have naturally evolved one or more missense mutations in bfmS, which encodes the sensor histidine kinase of the BfmRS two-component system (TCS). A mutant BfmS protein containing both the L181P and E376Q substitutions increased the phosphorylation and thus the transcriptional regulatory activity of its cognate downstream response regulator, BfmR. This reduced acute virulence and enhanced biofilm formation, both of which are phenotypic changes associated with a chronic infection state. The increased phosphorylation of BfmR was due, at least in part, to the cross-phosphorylation of BfmR by GtrS, a noncognate sensor kinase. Other spontaneous missense mutations in bfmS, such as A42E/G347D, T242R, and R393H, also caused a similar remodeling of the BfmRS TCS in P. aeruginosa This study highlights the plasticity of TCSs mediated by spontaneous mutations and suggests that mutation-induced activation of BfmRS may contribute to host adaptation by P. aeruginosa during chronic infections.
Collapse
Affiliation(s)
- Qiao Cao
- College of Life Science, Northwest University, Xi'an 710127, China.,State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai 201203, China
| | - Nana Yang
- State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai 201203, China
| | - Yanhui Wang
- State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai 201203, China
| | - Chenchen Xu
- State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai 201203, China
| | - Xue Zhang
- State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai 201203, China
| | - Ke Fan
- State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai 201203, China
| | - Feifei Chen
- College of Life Science, Northwest University, Xi'an 710127, China.,State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai 201203, China
| | - Haihua Liang
- College of Life Science, Northwest University, Xi'an 710127, China
| | - Yingchao Zhang
- Department of Biomedical Sciences, City University of Hong Kong, Kowloon Tong, Hong Kong 999077, China
| | - Xin Deng
- Department of Biomedical Sciences, City University of Hong Kong, Kowloon Tong, Hong Kong 999077, China
| | - Youjun Feng
- School of Medicine, Zhejiang University, Hangzhou 310058, China
| | - Cai-Guang Yang
- State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai 201203, China.,School of Pharmaceutical Science and Technology, Hangzhou Institute for Advanced Study, UCAS, Hangzhou 310024, China
| | - Min Wu
- Department of Biomedical Sciences, University of North Dakota, Grand Forks, ND 58203-9037, USA
| | - Taeok Bae
- Department of Microbiology and Immunology, Indiana University School of Medicine-Northwest, Gary, IN 46408, USA
| | - Lefu Lan
- College of Life Science, Northwest University, Xi'an 710127, China. .,State Key Laboratory of Drug Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai 201203, China.,School of Pharmaceutical Science and Technology, Hangzhou Institute for Advanced Study, UCAS, Hangzhou 310024, China.,NMPA Key Laboratory for Testing Technology of Pharmaceutical Microbiology, Shanghai Institute for Food and Drug Control, Shanghai, China
| |
Collapse
|
10
|
Abstract
Tumours vary in gene expression programmes and genetic alterations. Understanding this diversity and its biological meaning requires a theoretical framework, which could in turn guide the development of more accurate prognosis and therapy. Here, we review the theory of multi-task evolution of cancer, which is based upon the premise that tumours evolve in the host and face selection trade-offs between multiple biological functions. This theory can help identify the major biological tasks that cancer cells perform and the trade-offs between these tasks. It introduces the concept of specialist tumours, which focus on one task, and generalist tumours, which perform several tasks. Specialist tumours are suggested to be sensitive to therapy targeting their main task. Driver mutations tune gene expression towards specific tasks in a tissue-dependent manner and thus help to determine whether a tumour is specialist or generalist. We discuss potential applications of the theory of multi-task evolution to interpret the spatial organization of tumours and intratumour heterogeneity.
Collapse
Affiliation(s)
- Jean Hausser
- Department of Cellular and Molecular Biology, Karolinska Institutet, Solna, Sweden.
- SciLifeLab, Solna, Sweden.
| | - Uri Alon
- Department of Molecular and Cell Biology, Weizmann Institute of Science, Rehovot, Israel.
| |
Collapse
|
11
|
Alcacer A, Epifanio I, Ibáñez MV, Simó A, Ballester A. A data-driven classification of 3D foot types by archetypal shapes based on landmarks. PLoS One 2020; 15:e0228016. [PMID: 31999749 PMCID: PMC6991988 DOI: 10.1371/journal.pone.0228016] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2019] [Accepted: 01/05/2020] [Indexed: 12/29/2022] Open
Abstract
The taxonomy of foot shapes or other parts of the body is important, especially for design purposes. We propose a methodology based on archetypoid analysis (ADA) that overcomes the weaknesses of previous methodologies used to establish typologies. ADA is an objective, data-driven methodology that seeks extreme patterns, the archetypal profiles in the data. ADA also explains the data as percentages of the archetypal patterns, which makes this technique understandable and accessible even for non-experts. Clustering techniques are usually considered for establishing taxonomies, but we will show that finding the purest or most extreme patterns is more appropriate than using the central points returned by clustering techniques. We apply the methodology to an anthropometric database of 775 3D right foot scans representing the Spanish adult female and male population for footwear design. Each foot is described by a 5626 × 3 configuration matrix of landmarks. No multivariate features are used for establishing the taxonomy, but all the information gathered from the 3D scanning is employed. We use ADA for shapes described by landmarks. Women’s and men’s feet are analyzed separately. We have analyzed 3 archetypal feet for both men and women. These archetypal feet could not have been recovered using multivariate techniques.
Collapse
Affiliation(s)
- Aleix Alcacer
- Departament de Matemàtiques, Universitat Jaume I, Castelló, Spain
| | - Irene Epifanio
- Departament de Matemàtiques, Universitat Jaume I, Castelló, Spain
- Institut de Matemàtiques i Aplicacions de Castelló, Universitat Jaume I, Castelló, Spain
- * E-mail:
| | - M. Victoria Ibáñez
- Departament de Matemàtiques, Universitat Jaume I, Castelló, Spain
- Institut de Matemàtiques i Aplicacions de Castelló, Universitat Jaume I, Castelló, Spain
| | - Amelia Simó
- Departament de Matemàtiques, Universitat Jaume I, Castelló, Spain
- Institut de Matemàtiques i Aplicacions de Castelló, Universitat Jaume I, Castelló, Spain
| | | |
Collapse
|
12
|
Epifanio I, Ibáñez MV, Simó A. Archetypal Analysis With Missing Data: See All Samples by Looking at a Few Based on Extreme Profiles. AM STAT 2019. [DOI: 10.1080/00031305.2018.1545700] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Affiliation(s)
- Irene Epifanio
- Departament de Matemàtiques-IMAC, Universitat Jaume I, Spain
| | | | - Amelia Simó
- Departament de Matemàtiques-IMAC, Universitat Jaume I, Spain
| |
Collapse
|
13
|
Sheftel H, Szekely P, Mayo A, Sella G, Alon U. Evolutionary trade-offs and the structure of polymorphisms. Philos Trans R Soc Lond B Biol Sci 2019; 373:rstb.2017.0105. [PMID: 29632259 DOI: 10.1098/rstb.2017.0105] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/30/2017] [Indexed: 12/15/2022] Open
Abstract
Populations of organisms show genetic differences called polymorphisms. Understanding the effects of polymorphisms is important for biology and medicine. Here, we ask which polymorphisms occur at high frequency when organisms evolve under trade-offs between multiple tasks. Multiple tasks present a problem, because it is not possible to be optimal at all tasks simultaneously and hence compromises are necessary. Recent work indicates that trade-offs lead to a simple geometry of phenotypes in the space of traits: phenotypes fall on the Pareto front, which is shaped as a polytope: a line, triangle, tetrahedron etc. The vertices of these polytopes are the optimal phenotypes for a single task. Up to now, work on this Pareto approach has not considered its genetic underpinnings. Here, we address this by asking how the polymorphism structure of a population is affected by evolution under trade-offs. We simulate a multi-task selection scenario, in which the population evolves to the Pareto front: the line segment between two archetypes or the triangle between three archetypes. We find that polymorphisms that become prevalent in the population have pleiotropic phenotypic effects that align with the Pareto front. Similarly, epistatic effects between prevalent polymorphisms are parallel to the front. Alignment with the front occurs also for asexual mating. Alignment is reduced when drift or linkage is strong, and is replaced by a more complex structure in which many perpendicular allele effects cancel out. Aligned polymorphism structure allows mating to produce offspring that stand a good chance of being optimal multi-taskers in at least one of the locales available to the species.This article is part of the theme issue 'Self-organization in cell biology'.
Collapse
Affiliation(s)
- Hila Sheftel
- Department Molecular Cell Biology, Weizmann Institute of Science, Rehovot 76100, Israel
| | - Pablo Szekely
- Department Molecular Cell Biology, Weizmann Institute of Science, Rehovot 76100, Israel
| | - Avi Mayo
- Department Molecular Cell Biology, Weizmann Institute of Science, Rehovot 76100, Israel
| | - Guy Sella
- Department of Biological Sciences, Columbia University, New York, NY 10027, USA
| | - Uri Alon
- Department Molecular Cell Biology, Weizmann Institute of Science, Rehovot 76100, Israel
| |
Collapse
|
14
|
Bartell JA, Sommer LM, Haagensen JAJ, Loch A, Espinosa R, Molin S, Johansen HK. Evolutionary highways to persistent bacterial infection. Nat Commun 2019; 10:629. [PMID: 30733448 PMCID: PMC6367392 DOI: 10.1038/s41467-019-08504-7] [Citation(s) in RCA: 70] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2018] [Accepted: 01/10/2019] [Indexed: 01/18/2023] Open
Abstract
Persistent infections require bacteria to evolve from their naïve colonization state by optimizing fitness in the host via simultaneous adaptation of multiple traits, which can obscure evolutionary trends and complicate infection management. Accordingly, here we screen 8 infection-relevant phenotypes of 443 longitudinal Pseudomonas aeruginosa isolates from 39 young cystic fibrosis patients over 10 years. Using statistical modeling, we map evolutionary trajectories and identify trait correlations accounting for patient-specific influences. By integrating previous genetic analyses of 474 isolates, we provide a window into early adaptation to the host, finding: (1) a 2–3 year timeline of rapid adaptation after colonization, (2) variant “naïve” and “adapted” states reflecting discordance between phenotypic and genetic adaptation, (3) adaptive trajectories leading to persistent infection via three distinct evolutionary modes, and (4) new associations between phenotypes and pathoadaptive mutations. Ultimately, we effectively deconvolute complex trait adaptation, offering a framework for evolutionary studies and precision medicine in clinical microbiology. The pathogen Pseudomonas aeruginosa undergoes complex trait adaptation within cystic fibrosis patients. Here, Bartell, Sommer, and colleagues use statistical modeling of longitudinal isolates to characterize the joint genetic and phenotypic evolutionary trajectories of P. aeruginosa within hosts.
Collapse
Affiliation(s)
- Jennifer A Bartell
- The Novo Nordisk Foundation Center for Biosustainability, Technical University of Denmark, 2800, Kgs. Lyngby, Denmark.
| | - Lea M Sommer
- Department of Clinical Microbiology, Rigshospitalet, 2100, Copenhagen Ø, Denmark.
| | - Janus A J Haagensen
- The Novo Nordisk Foundation Center for Biosustainability, Technical University of Denmark, 2800, Kgs. Lyngby, Denmark
| | - Anne Loch
- The Novo Nordisk Foundation Center for Biosustainability, Technical University of Denmark, 2800, Kgs. Lyngby, Denmark
| | - Rocio Espinosa
- The Novo Nordisk Foundation Center for Biosustainability, Technical University of Denmark, 2800, Kgs. Lyngby, Denmark
| | - Søren Molin
- The Novo Nordisk Foundation Center for Biosustainability, Technical University of Denmark, 2800, Kgs. Lyngby, Denmark
| | - Helle Krogh Johansen
- Department of Clinical Microbiology, Rigshospitalet, 2100, Copenhagen Ø, Denmark.,Department of Clinical Medicine, Faculty of Health and Medical Sciences, University of Copenhagen, 2200, Copenhagen N, Denmark
| |
Collapse
|
15
|
Cona G, Koçillari L, Palombit A, Bertoldo A, Maritan A, Corbetta M. Archetypes of human cognition defined by time preference for reward and their brain correlates: An evolutionary trade-off approach. Neuroimage 2018; 185:322-334. [PMID: 30355533 DOI: 10.1016/j.neuroimage.2018.10.050] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2018] [Revised: 10/05/2018] [Accepted: 10/18/2018] [Indexed: 01/24/2023] Open
Abstract
Biological systems carry out multiple tasks in their lifetime, which, in the course of evolution, may lead to trade-offs. In fact phenotypes (different species, individuals within a species, circuits, bacteria, proteins, etc.) cannot be optimal at all tasks, and, according to Pareto optimality theory, lay into a well-defined geometrical distribution (polygons and/or polyhedrons) in the space of traits. The vertices of this distribution contain archetypes, namely phenotypes that are specialists at one of the tasks, whereas phenotypes toward the center of the geometrical distribution show average performance across tasks. We applied this theory to the variability of cognitive and behavioral scores measured in 1206 individuals from the Human Connectome Project. Among all possible combinations of pairs of traits, we found the best fit to Pareto optimality when individuals were plotted in the trait-space of time preferences for reward, evaluated with the Delay Discounting Task (DDT). The DDT measures subjects' preference in choosing either immediate smaller rewards or delayed larger rewards. Time preference for reward was described by a triangular distribution in which each of the three vertices included individuals who used a particular strategy to discount reward. These archetypes accounted for variability on many cognitive, personality, and socioeconomic status variables, as well as differences in brain structure and functional connectivity, with only a weak influence of genetics. In summary, time preference for reward reflects a core variable that biases human phenotypes via natural and cultural selection.
Collapse
Affiliation(s)
- Giorgia Cona
- Department of General Psychology, University of Padua, Italy; Padova Neuroscience Center (PNC), University of Padua, Italy
| | - Loren Koçillari
- Department of Physics, University of Padua, Italy; Padova Neuroscience Center (PNC), University of Padua, Italy
| | - Alessandro Palombit
- Department of Information Engineering, University of Padua, Italy; Padova Neuroscience Center (PNC), University of Padua, Italy
| | - Alessandra Bertoldo
- Department of Information Engineering, University of Padua, Italy; Padova Neuroscience Center (PNC), University of Padua, Italy
| | - Amos Maritan
- Department of Physics, University of Padua, Italy; Padova Neuroscience Center (PNC), University of Padua, Italy
| | - Maurizio Corbetta
- Department of Neuroscience, University of Padua, Italy; Departments of Neurology, Radiology, Neuroscience, Washington University School of Medicine, Saint Louis, USA; Padova Neuroscience Center (PNC), University of Padua, Italy.
| |
Collapse
|
16
|
Swann E, Sun B, Cleland DM, Barnard AS. Representing molecular and materials data for unsupervised machine learning. MOLECULAR SIMULATION 2018. [DOI: 10.1080/08927022.2018.1450982] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
Affiliation(s)
- E. Swann
- Molecular and Materials Modelling, Data61 CSIRO , Docklands, Victoria, Australia
| | - B. Sun
- Molecular and Materials Modelling, Data61 CSIRO , Docklands, Victoria, Australia
| | - D. M. Cleland
- Molecular and Materials Modelling, Data61 CSIRO , Docklands, Victoria, Australia
| | - A. S. Barnard
- Molecular and Materials Modelling, Data61 CSIRO , Docklands, Victoria, Australia
| |
Collapse
|
17
|
Epifanio I, Ibáñez MV, Simó A. Archetypal shapes based on landmarks and extension to handle missing data. ADV DATA ANAL CLASSI 2017. [DOI: 10.1007/s11634-017-0297-7] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
18
|
Reeve J, Böhmig GA, Eskandary F, Einecke G, Lefaucheur C, Loupy A, Halloran PF. Assessing rejection-related disease in kidney transplant biopsies based on archetypal analysis of molecular phenotypes. JCI Insight 2017; 2:94197. [PMID: 28614805 DOI: 10.1172/jci.insight.94197] [Citation(s) in RCA: 104] [Impact Index Per Article: 14.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2017] [Accepted: 05/05/2017] [Indexed: 01/15/2023] Open
Abstract
Conventional histologic diagnosis of rejection in kidney transplants has limited repeatability due to its inherent requirement for subjective assessment of lesions, in a rule-based system that does not acknowledge diagnostic uncertainty. Molecular phenotyping affords opportunities for increased precision and improved disease classification to address the limitations of conventional histologic diagnostic systems and quantify levels of uncertainty. Microarray data from 1,208 kidney transplant biopsies were collected prospectively from 13 centers. Cross-validated classifier scores predicting the presence of antibody-mediated rejection (ABMR), T cell-mediated rejection (TCMR), and 5 related histologic lesions were generated using supervised machine learning methods. These scores were used as input for archetypal analysis, an unsupervised method similar to cluster analysis, to examine the distribution of molecular phenotypes related to rejection. Six archetypes were generated: no rejection, TCMR, 3 associated with ABMR (early-stage, fully developed, and late-stage), and mixed rejection (TCMR plus early-stage ABMR). Each biopsy was assigned 6 scores, one for each archetype, representing a probabilistic assessment of that biopsy based on its rejection-related molecular properties. Viewed as clusters, the archetypes were similar to existing histologic Banff categories, but there was 32% disagreement, much of it probably reflecting the "noise" in the current histologic assessment system. Graft survival was lowest for fully developed and late-stage ABMR, and it was better predicted by molecular archetype scores than histologic diagnoses. The results provide a system for precision molecular assessment of biopsies and a new standard for recalibrating conventional diagnostic systems.
Collapse
Affiliation(s)
- Jeff Reeve
- Alberta Transplant Applied Genomics Centre, Alberta, Canada.,Department of Laboratory Medicine and Pathology, University of Alberta, Edmonton, Alberta, Canada
| | - Georg A Böhmig
- Division of Nephrology and Dialysis, Department of Medicine III, Medical University of Vienna, Vienna, Austria
| | - Farsad Eskandary
- Division of Nephrology and Dialysis, Department of Medicine III, Medical University of Vienna, Vienna, Austria
| | - Gunilla Einecke
- Department of Nephrology, Medizinische Hochschule Hannover, Hannover, Germany
| | - Carmen Lefaucheur
- Paris Translational Research Center for Organ Transplantation, INSERM, UMR-S970, Paris, France.,Saint-Louis Hospital, Assistance Publique-Hôpitaux de Paris, Paris, France
| | - Alexandre Loupy
- Paris Translational Research Center for Organ Transplantation, INSERM, UMR-S970, Paris, France.,Necker Hospital, Assistance Publique-Hôpitaux de Paris, Paris, France
| | - Philip F Halloran
- Alberta Transplant Applied Genomics Centre, Alberta, Canada.,Department of Medicine, Division of Nephrology and Transplant Immunology, University of Alberta, Edmonton, Alberta, Canada
| | | |
Collapse
|
19
|
|
20
|
Sun B, Fernandez M, Barnard AS. Statistics, damned statistics and nanoscience - using data science to meet the challenge of nanomaterial complexity. NANOSCALE HORIZONS 2016; 1:89-95. [PMID: 32260631 DOI: 10.1039/c5nh00126a] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
For many years dealing with the complexity of nanoscale materials, the polydispersivity of individual samples, and the persistent imperfection of individual nanostructures has been secondary to our search for novel properties and promising applications. For our science to translate into technology, however, we will inevitably need to deal with the issue of structural diversity and integrate this feature into the next generation of more realistic structure/property predictions. This is challenging in the field of nanoscience where atomic level precision is typically inaccessible (experimentally), but properties can depend on structural variations at the atomic scale. Fortunately there exists a range of reliable statistical methods that are entirely applicable to nanoscale materials; ideal for navigating and analysing enormous amount of information required to accurately describe realistic samples. Combined with advances in automation and information technology the field of data science can assist us in dealing with our big data, characterising our uncertainties, and more rapidly identifying useful structure/property relationships. Taking greater advantage of data-driven methods involves thinking differently about our research, but applied appropriately these methods can accelerate the discovery of nanomaterials that are optimised to make the transition from science to technology.
Collapse
Affiliation(s)
- Baichuan Sun
- CSIRO Virtual Nanoscience Laboratory, Parkville, VIC 3052, Australia.
| | | | | |
Collapse
|
21
|
Fernandez M, Barnard AS. Identification of Nanoparticle Prototypes and Archetypes. ACS NANO 2015; 9:11980-11992. [PMID: 26575441 DOI: 10.1021/acsnano.5b05788] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
High-throughput (HT) computational characterization of nanomaterials is poised to accelerate novel material breakthroughs. The number of possible nanomaterials is increasing exponentially along with their complexity, and so statistical and information technology will play a fundamental role in rationalizing nanomaterials HT data. We demonstrate that multivariate statistical analysis of heterogeneous ensembles can identify the truly significant nanoparticles and their most relevant properties. Virtual samples of diamond nanoparticles and graphene nanoflakes are characterized using clustering and archetypal analysis, where we find that saturated particles are defined by their geometry, while nonsaturated nanoparticles are defined by their carbon chemistry. At the complex hull of the nanostructure spaces, a combination of complex archetypes can efficiency describe a large number of members of the ensembles, whereas the regular shapes that are typically assumed to be representative can only describe a small set of the most regular morphologies. This approach provides a route toward the characterization of computationally intractable virtual nanomaterial spaces, which can aid nanomaterials discovery in the foreseen big data scenario.
Collapse
Affiliation(s)
- Michael Fernandez
- CSIRO Virtual Nanoscience Laboratory , 343 Royal Parade, Parkville, Victoria 3052, Australia
| | - Amanda S Barnard
- CSIRO Virtual Nanoscience Laboratory , 343 Royal Parade, Parkville, Victoria 3052, Australia
| |
Collapse
|
22
|
Szekely P, Korem Y, Moran U, Mayo A, Alon U. The Mass-Longevity Triangle: Pareto Optimality and the Geometry of Life-History Trait Space. PLoS Comput Biol 2015; 11:e1004524. [PMID: 26465336 PMCID: PMC4605829 DOI: 10.1371/journal.pcbi.1004524] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2014] [Accepted: 08/26/2015] [Indexed: 12/14/2022] Open
Abstract
When organisms need to perform multiple tasks they face a fundamental tradeoff: no phenotype can be optimal at all tasks. This situation was recently analyzed using Pareto optimality, showing that tradeoffs between tasks lead to phenotypes distributed on low dimensional polygons in trait space. The vertices of these polygons are archetypes--phenotypes optimal at a single task. This theory was applied to examples from animal morphology and gene expression. Here we ask whether Pareto optimality theory can apply to life history traits, which include longevity, fecundity and mass. To comprehensively explore the geometry of life history trait space, we analyze a dataset of life history traits of 2105 endothermic species. We find that, to a first approximation, life history traits fall on a triangle in log-mass log-longevity space. The vertices of the triangle suggest three archetypal strategies, exemplified by bats, shrews and whales, with specialists near the vertices and generalists in the middle of the triangle. To a second approximation, the data lies in a tetrahedron, whose extra vertex above the mass-longevity triangle suggests a fourth strategy related to carnivory. Each animal species can thus be placed in a coordinate system according to its distance from the archetypes, which may be useful for genome-scale comparative studies of mammalian aging and other biological aspects. We further demonstrate that Pareto optimality can explain a range of previous studies which found animal and plant phenotypes which lie in triangles in trait space. This study demonstrates the applicability of multi-objective optimization principles to understand life history traits and to infer archetypal strategies that suggest why some mammalian species live much longer than others of similar mass.
Collapse
Affiliation(s)
- Pablo Szekely
- Department of Molecular Cell Biology, The Weizmann Institute of Science, Rehovot, Israel
| | - Yael Korem
- Department of Molecular Cell Biology, The Weizmann Institute of Science, Rehovot, Israel
| | - Uri Moran
- Department of Plant Science, The Weizmann Institute of Science, Rehovot, Israel
| | - Avi Mayo
- Department of Molecular Cell Biology, The Weizmann Institute of Science, Rehovot, Israel
| | - Uri Alon
- Department of Molecular Cell Biology, The Weizmann Institute of Science, Rehovot, Israel
- * E-mail:
| |
Collapse
|
23
|
Korem Y, Szekely P, Hart Y, Sheftel H, Hausser J, Mayo A, Rothenberg ME, Kalisky T, Alon U. Geometry of the Gene Expression Space of Individual Cells. PLoS Comput Biol 2015; 11:e1004224. [PMID: 26161936 PMCID: PMC4498931 DOI: 10.1371/journal.pcbi.1004224] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2014] [Accepted: 03/04/2015] [Indexed: 12/14/2022] Open
Abstract
There is a revolution in the ability to analyze gene expression of single cells in a tissue. To understand this data we must comprehend how cells are distributed in a high-dimensional gene expression space. One open question is whether cell types form discrete clusters or whether gene expression forms a continuum of states. If such a continuum exists, what is its geometry? Recent theory on evolutionary trade-offs suggests that cells that need to perform multiple tasks are arranged in a polygon or polyhedron (line, triangle, tetrahedron and so on, generally called polytopes) in gene expression space, whose vertices are the expression profiles optimal for each task. Here, we analyze single-cell data from human and mouse tissues profiled using a variety of single-cell technologies. We fit the data to shapes with different numbers of vertices, compute their statistical significance, and infer their tasks. We find cases in which single cells fill out a continuum of expression states within a polyhedron. This occurs in intestinal progenitor cells, which fill out a tetrahedron in gene expression space. The four vertices of this tetrahedron are each enriched with genes for a specific task related to stemness and early differentiation. A polyhedral continuum of states is also found in spleen dendritic cells, known to perform multiple immune tasks: cells fill out a tetrahedron whose vertices correspond to key tasks related to maturation, pathogen sensing and communication with lymphocytes. A mixture of continuum-like distributions and discrete clusters is found in other cell types, including bone marrow and differentiated intestinal crypt cells. This approach can be used to understand the geometry and biological tasks of a wide range of single-cell datasets. The present results suggest that the concept of cell type may be expanded. In addition to discreet clusters in gene-expression space, we suggest a new possibility: a continuum of states within a polyhedron, in which the vertices represent specialists at key tasks. In the past, biological experiments usually pooled together millions of cells, masking the differences between individual cells. Current technology takes a big step forward by measuring gene expression from individual cells. Interpreting this data is challenging because we need to understand how cells are arranged in a high dimensional gene expression space. Here we test recent theory that suggests that cells facing multiple tasks should be arranged in simple low dimensional polygons or polyhedra (generally called polytopes). The vertices of the polytopes are gene expression profiles optimal for each of the tasks. We find evidence for such simplicity in a variety of tissues—spleen, bone marrow, intestine—analyzed by different single-cell technologies. We find that cells are distributed inside polytopes, such as tetrahedrons or four-dimensional simplexes, with cells closest to each vertex responsible for a different key task. For example, intestinal progenitor cells that give rise to the other cell types show a continuous distribution in a tetrahedron whose vertices correspond to several key sub-tasks. Immune dendritic cells likewise are continuously distributed between key immune tasks. This approach of testing whether data falls in polytopes may be useful for interpreting a variety of single-cell datasets in terms of biological tasks.
Collapse
Affiliation(s)
- Yael Korem
- Department of Molecular Cell Biology, Weizmann Institute of Science, Rehovot, Israel
| | - Pablo Szekely
- Department of Molecular Cell Biology, Weizmann Institute of Science, Rehovot, Israel
| | - Yuval Hart
- Department of Molecular Cell Biology, Weizmann Institute of Science, Rehovot, Israel
| | - Hila Sheftel
- Department of Molecular Cell Biology, Weizmann Institute of Science, Rehovot, Israel
| | - Jean Hausser
- Department of Molecular Cell Biology, Weizmann Institute of Science, Rehovot, Israel
| | - Avi Mayo
- Department of Molecular Cell Biology, Weizmann Institute of Science, Rehovot, Israel
| | - Michael E. Rothenberg
- Department of Medicine, Division of Gastroenterology and Hepatology, and Institute for Stem Cell Biology and Regenerative Medicine, Stanford University, Stanford, California, United States of America
| | - Tomer Kalisky
- Faculty of Engineering, Bar Ilan University, Ramat Gan, Israel
| | - Uri Alon
- Department of Molecular Cell Biology, Weizmann Institute of Science, Rehovot, Israel
- * E-mail:
| |
Collapse
|
24
|
Vinué G, Epifanio I, Alemany S. Archetypoids: A new approach to define representative archetypal data. Comput Stat Data Anal 2015. [DOI: 10.1016/j.csda.2015.01.018] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
25
|
|
26
|
Tendler A, Mayo A, Alon U. Evolutionary tradeoffs, Pareto optimality and the morphology of ammonite shells. BMC SYSTEMS BIOLOGY 2015; 9:12. [PMID: 25884468 PMCID: PMC4404009 DOI: 10.1186/s12918-015-0149-z] [Citation(s) in RCA: 49] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/17/2014] [Accepted: 01/29/2015] [Indexed: 11/24/2022]
Abstract
Background Organisms that need to perform multiple tasks face a fundamental tradeoff: no design can be optimal at all tasks at once. Recent theory based on Pareto optimality showed that such tradeoffs lead to a highly defined range of phenotypes, which lie in low-dimensional polyhedra in the space of traits. The vertices of these polyhedra are called archetypes- the phenotypes that are optimal at a single task. To rigorously test this theory requires measurements of thousands of species over hundreds of millions of years of evolution. Ammonoid fossil shells provide an excellent model system for this purpose. Ammonoids have a well-defined geometry that can be parameterized using three dimensionless features of their logarithmic-spiral-shaped shells. Their evolutionary history includes repeated mass extinctions. Results We find that ammonoids fill out a pyramid in morphospace, suggesting five specific tasks - one for each vertex of the pyramid. After mass extinctions, surviving species evolve to refill essentially the same pyramid, suggesting that the tasks are unchanging. We infer putative tasks for each archetype, related to economy of shell material, rapid shell growth, hydrodynamics and compactness. Conclusions These results support Pareto optimality theory as an approach to study evolutionary tradeoffs, and demonstrate how this approach can be used to infer the putative tasks that may shape the natural selection of phenotypes. Electronic supplementary material The online version of this article (doi:10.1186/s12918-015-0149-z) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Avichai Tendler
- Department of Molecular cell biology, Weizmann Institute of Science, Rehovot, 76100, Israel.
| | - Avraham Mayo
- Department of Molecular cell biology, Weizmann Institute of Science, Rehovot, 76100, Israel.
| | - Uri Alon
- Department of Molecular cell biology, Weizmann Institute of Science, Rehovot, 76100, Israel.
| |
Collapse
|
27
|
Inferring biological tasks using Pareto analysis of high-dimensional data. Nat Methods 2015; 12:233-5, 3 p following 235. [PMID: 25622107 DOI: 10.1038/nmeth.3254] [Citation(s) in RCA: 95] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2014] [Accepted: 11/17/2014] [Indexed: 01/16/2023]
Abstract
We present the Pareto task inference method (ParTI; http://www.weizmann.ac.il/mcb/UriAlon/download/ParTI) for inferring biological tasks from high-dimensional biological data. Data are described as a polytope, and features maximally enriched closest to the vertices (or archetypes) allow identification of the tasks the vertices represent. We demonstrate that human breast tumors and mouse tissues are well described by tetrahedrons in gene expression space, with specific tumor types and biological functions enriched at each of the vertices, suggesting four key tasks.
Collapse
|
28
|
Hosseinkhan N, Zarrineh P, Masoudi-Nejad A. Analysis of Genome-scale Expression Network in Four Major Bacterial Residents of Cystic Fibrosis Lung. Curr Genomics 2014; 15:408-18. [PMID: 25435803 PMCID: PMC4245700 DOI: 10.2174/1389202915666140818205444] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2013] [Revised: 07/06/2014] [Accepted: 08/15/2014] [Indexed: 11/22/2022] Open
Abstract
In polymicrobial communities where several species co-exist in a certain niche and consequently the possibility of interactions among species is very high, gene expression data sources can give better insights in to underlying adaptation mechanisms assumed by bacteria. Furthermore, several possible synergistic or antagonistic interactions among species can be investigated through gene expression comparisons. Lung is one of the habitats harboring several distinct pathogens during severe pulmonary disorders such as chronic obstructive pulmonary disease (COPD) and cystic fibrosis (CF). Expression data analysis of these lung residents can help to gain a better understanding on how these species interact with each other within the host cells. The first part of this paper deals with introducing available data sources for the major bacteria responsible for causing lung diseases and their genomic relations. In the second part, the main focus is on the studies concerning gene expression analyses of these species.
Collapse
Affiliation(s)
- Nazanin Hosseinkhan
- Laboratory of Systems Biology and Bioinformatics (LBB), Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran
| | - Peyman Zarrineh
- Laboratory of Systems Biology and Bioinformatics (LBB), Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran
| | - Ali Masoudi-Nejad
- Laboratory of Systems Biology and Bioinformatics (LBB), Institute of Biochemistry and Biophysics, University of Tehran, Tehran, Iran
| |
Collapse
|
29
|
Environmental heterogeneity drives within-host diversification and evolution of Pseudomonas aeruginosa. mBio 2014; 5:e01592-14. [PMID: 25227464 PMCID: PMC4172072 DOI: 10.1128/mbio.01592-14] [Citation(s) in RCA: 123] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
Microbial population polymorphisms are commonly observed in natural environments, including long-term infected hosts. However, the underlying processes promoting and stabilizing diversity are difficult to unravel and are not well understood. Here, we use chronic infection of cystic fibrosis airways by the opportunistic pathogen Pseudomonas aeruginosa as a system for investigating bacterial diversification processes during the course of infection. We analyze clonal bacterial isolates sampled during a 32-year period and map temporal and spatial variations in population diversity to different infection sites within the infected host. We show that the ancestral infecting strain diverged into distinct sublineages, each with their own functional and genomic signatures and rates of adaptation, immediately after initial colonization. The sublineages coexisted in the host for decades, suggesting rapid evolution of stable population polymorphisms. Critically, the observed generation and maintenance of population diversity was the result of partitioning of the sublineages into physically separated niches in the CF airway. The results reveal a complex within-host population structure not previously realized and provide evidence that the heterogeneity of the highly structured and complex host environment promotes the evolution and long-term stability of pathogen population diversity during infection. Within-host pathogen evolution and diversification during the course of chronic infections is of importance in relation to therapeutic intervention strategies, yet our understanding of these processes is limited. Here, we investigate intraclonal population diversity in P. aeruginosa during chronic airway infections in cystic fibrosis patients. We show the evolution of a diverse population structure immediately after initial colonization, with divergence into multiple distinct sublineages that coexisted for decades and occupied distinct niches. Our results suggest that the spatial heterogeneity in CF airways plays a major role in relation to the generation and maintenance of population diversity and emphasize that a single isolate in sputum may not represent the entire pathogen population in the infected individual. A more complete understanding of the evolution of distinct clonal variants and their distribution in different niches could have positive implications for efficient therapy.
Collapse
|
30
|
de Lorenzo V. From theselfish genetoselfish metabolism: Revisiting the central dogma. Bioessays 2014; 36:226-35. [DOI: 10.1002/bies.201300153] [Citation(s) in RCA: 51] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Víctor de Lorenzo
- Systems & Synthetic Biology Program; Centro Nacional de Biotecnología CSIC Cantoblanco; Madrid Spain
| |
Collapse
|