Towards robust network based complex systems : from evolutionary cellular automata to biological models

Abstract Sitting between your past and your future doesn't mean you are in the present. Dakota Skye Complex systems science is an interdisciplinary field grouping under the same umbrella dynamical phenomena from social, natural or mathematical sciences. The emergence of a higher order organization or behavior, transcending that expected of the linear addition of the parts, is a key factor shared by all these systems. Most complex systems can be modeled as networks that represent the interactions amongst the system's components. In addition to the actual nature of the part's interactions, the intrinsic topological structure of underlying network is believed to play a crucial role in the remarkable emergent behaviors exhibited by the systems. Moreover, the topology is also a key a factor to explain the extraordinary flexibility and resilience to perturbations when applied to transmission and diffusion phenomena. In this work, we study the effect of different network structures on the performance and on the fault tolerance of systems in two different contexts. In the first part, we study cellular automata, which are a simple paradigm for distributed computation. Cellular automata are made of basic Boolean computational units, the cells; relying on simple rules and information from- the surrounding cells to perform a global task. The limited visibility of the cells can be modeled as a network, where interactions amongst cells are governed by an underlying structure, usually a regular one. In order to increase the performance of cellular automata, we chose to change its topology. We applied computational principles inspired by Darwinian evolution, called evolutionary algorithms, to alter the system's topological structure starting from either a regular or a random one. The outcome is remarkable, as the resulting topologies find themselves sharing properties of both regular and random network, and display similitudes Watts-Strogtz's small-world network found in social systems. Moreover, the performance and tolerance to probabilistic faults of our small-world like cellular automata surpasses that of regular ones. In the second part, we use the context of biological genetic regulatory networks and, in particular, Kauffman's random Boolean networks model. In some ways, this model is close to cellular automata, although is not expected to perform any task. Instead, it simulates the time-evolution of genetic regulation within living organisms under strict conditions. The original model, though very attractive by it's simplicity, suffered from important shortcomings unveiled by the recent advances in genetics and biology. We propose to use these new discoveries to improve the original model. Firstly, we have used artificial topologies believed to be closer to that of gene regulatory networks. We have also studied actual biological organisms, and used parts of their genetic regulatory networks in our models. Secondly, we have addressed the improbable full synchronicity of the event taking place on. Boolean networks and proposed a more biologically plausible cascading scheme. Finally, we tackled the actual Boolean functions of the model, i.e. the specifics of how genes activate according to the activity of upstream genes, and presented a new update function that takes into account the actual promoting and repressing effects of one gene on another. Our improved models demonstrate the expected, biologically sound, behavior of previous GRN model, yet with superior resistance to perturbations. We believe they are one step closer to the biological reality.

Artificial neural network study of whole-cell bacterial bioreporter response determined using fluorescence flow cytometry.

Genetically engineered bioreporters are an excellent complement to traditional methods of chemical analysis. The application of fluorescence flow cytometry to detection of bioreporter response enables rapid and efficient characterization of bacterial bioreporter population response on a single-cell basis. In the present study, intrapopulation response variability was used to obtain higher analytical sensitivity and precision. We have analyzed flow cytometric data for an arsenic-sensitive bacterial bioreporter using an artificial neural network-based adaptive clustering approach (a single-layer perceptron model). Results for this approach are far superior to other methods that we have applied to this fluorescent bioreporter (e.g., the arsenic detection limit is 0.01 microM, substantially lower than for other detection methods/algorithms). The approach is highly efficient computationally and can be implemented on a real-time basis, thus having potential for future development of high-throughput screening applications.

Evolution des transferts interhospitaliers au départ d'un centre hospitalier universitaire suisse [Trends in interhospital transfers from a Swiss university hospital center].

OBJECTIVE: Research on interhospital transfers provides a basis for describing and quantifying patient flow and its evolution over time, offering an insight into hospital organization and management and hospital overcrowding. The purpose of this study was to conduct a qualitative and quantitative analysis of patient flow and to examine trends over an eight-year period. METHODS: A retrospective descriptive study of interhospital transfers was conducted between 2003 and 2011 based on an analysis of demographic, medical and operational characteristics. Ambulance transfers and transfers requiring physician assistance were analyzed separately. RESULTS: The number of interhospital transfers increased significantly over the study period,from 4,026 in 2003 to 6,481 in 2011 (+60.9%). The number of ambulance transfers increased by almost 300% (616 in 2003 compared to 2,460 in 2011). Most of the transfers (98%) were to hospitals located less than 75 km from the university hospital (median: 24 km, 5-44). In 2011, 24% of all transfers were to psychiatric institutions. 26% of all transfer cases were direct transfers from the emergency department. An increasing number of transfers required physician assistance. 18% of these patients required ventilatory support, whole 9.8% required vasoactive drugs. 11.6% of these transfers were due to hospital overcrowding. Conclusion: The study shows that there has been a significant increase in interhospital transfers. This increase is related to hospital overcrowding and to the network-based systems governing patient care strategies.

Connectivity differences in brain networks.

The scenario considered here is one where brain connectivity is represented as a network and an experimenter wishes to assess the evidence for an experimental effect at each of the typically thousands of connections comprising the network. To do this, a univariate model is independently fitted to each connection. It would be unwise to declare significance based on an uncorrected threshold of α=0.05, since the expected number of false positives for a network comprising N=90 nodes and N(N-1)/2=4005 connections would be 200. Control of Type I errors over all connections is therefore necessary. The network-based statistic (NBS) and spatial pairwise clustering (SPC) are two distinct methods that have been used to control family-wise errors when assessing the evidence for an experimental effect with mass univariate testing. The basic principle of the NBS and SPC is the same as supra-threshold voxel clustering. Unlike voxel clustering, where the definition of a voxel cluster is unambiguous, 'clusters' formed among supra-threshold connections can be defined in different ways. The NBS defines clusters using the graph theoretical concept of connected components. SPC on the other hand uses a more stringent pairwise clustering concept. The purpose of this article is to compare the pros and cons of the NBS and SPC, provide some guidelines on their practical use and demonstrate their utility using a case study involving neuroimaging data.

Transition to parenthood: The role of social interaction and endogenous networks

Empirical studies indicate that the transition to parenthood is influenced by an individual's peer group. To study the mechanisms creating interdepen- dencies across individuals' transition to parenthood and its timing we apply an agent-based simulation model. We build a one-sex model and provide agents with three different characteristics regarding age, intended education and parity. Agents endogenously form their network based on social closeness. Network members then may influence the agents' transition to higher parity levels. Our numerical simulations indicate that accounting for social inter- actions can explain the shift of first-birth probabilities in Austria over the period 1984 to 2004. Moreover, we apply our model to forecast age-specific fertility rates up to 2016.

The HIV care cascade in Switzerland: reaching the UNAIDS/WHO targets for patients diagnosed with HIV.

OBJECTIVES: To describe the HIV care cascade for Switzerland in the year 2012. DESIGN/METHODS: Six levels were defined: (i) HIV-infected, (ii) HIV-diagnosed, (iii) linked to care, (iv) retained in care, (v) on antiretroviral treatment (ART), and (vi) with suppressed viral load. We used data from the Swiss HIV Cohort Study (SHCS) complemented by a nationwide survey among SHCS physicians to estimate the number of HIV-patients not registered in the cohort. We also used Swiss ART sales data to estimate the number of patients treated outside the SHCS network. Based on the number of patients retained in care, we inferred the estimates for levels (i) to (iii) from previously published data. RESULTS: We estimate that (i) 15 200 HIV-infected individuals lived in Switzerland in 2012 (margins of uncertainty, 13 400-19 300). Of those, (ii) 12 300 (81%) were diagnosed, (iii) 12 200 (80%) linked, and (iv) 11 900 (79%) retained in care. Broadly based on SHCS network data, (v) 10 800 (71%) patients were receiving ART, and (vi) 10 400 (68%) had suppressed (<200 copies/ml) viral loads. The vast majority (95%) of patients retained in care were followed within the SHCS network, with 76% registered in the cohort. CONCLUSION: Our estimate for HIV-infected individuals in Switzerland is substantially lower than previously reported, halving previous national HIV prevalence estimates to 0.2%. In Switzerland in 2012, 91% of patients in care were receiving ART, and 96% of patients on ART had suppressed viral load, meeting recent UNAIDS/WHO targets.

EEG-based functional brain networks: does the network size matter?

Functional connectivity in human brain can be represented as a network using electroencephalography (EEG) signals. These networks--whose nodes can vary from tens to hundreds--are characterized by neurobiologically meaningful graph theory metrics. This study investigates the degree to which various graph metrics depend upon the network size. To this end, EEGs from 32 normal subjects were recorded and functional networks of three different sizes were extracted. A state-space based method was used to calculate cross-correlation matrices between different brain regions. These correlation matrices were used to construct binary adjacency connectomes, which were assessed with regards to a number of graph metrics such as clustering coefficient, modularity, efficiency, economic efficiency, and assortativity. We showed that the estimates of these metrics significantly differ depending on the network size. Larger networks had higher efficiency, higher assortativity and lower modularity compared to those with smaller size and the same density. These findings indicate that the network size should be considered in any comparison of networks across studies.

A protein interaction atlas for the nuclear receptors: properties and quality of a hub-based dimerisation network.

BACKGROUND: The nuclear receptors are a large family of eukaryotic transcription factors that constitute major pharmacological targets. They exert their combinatorial control through homotypic heterodimerisation. Elucidation of this dimerisation network is vital in order to understand the complex dynamics and potential cross-talk involved. RESULTS: Phylogeny, protein-protein interactions, protein-DNA interactions and gene expression data have been integrated to provide a comprehensive and up-to-date description of the topology and properties of the nuclear receptor interaction network in humans. We discriminate between DNA-binding and non-DNA-binding dimers, and provide a comprehensive interaction map, that identifies potential cross-talk between the various pathways of nuclear receptors. CONCLUSION: We infer that the topology of this network is hub-based, and much more connected than previously thought. The hub-based topology of the network and the wide tissue expression pattern of NRs create a highly competitive environment for the common heterodimerising partners. Furthermore, a significant number of negative feedback loops is present, with the hub protein SHP [NR0B2] playing a major role. We also compare the evolution, topology and properties of the nuclear receptor network with the hub-based dimerisation network of the bHLH transcription factors in order to identify both unique themes and ubiquitous properties in gene regulation. In terms of methodology, we conclude that such a comprehensive picture can only be assembled by semi-automated text-mining, manual curation and integration of data from various sources.

An interaction network of the mammalian COP9 signalosome identifies Dda1 as a core subunit of multiple Cul4-based E3 ligases.

The COP9 signalosome (CSN) is an evolutionarily conserved macromolecular complex that interacts with cullin-RING E3 ligases (CRLs) and regulates their activity by hydrolyzing cullin-Nedd8 conjugates. The CSN sequesters inactive CRL4(Ddb2), which rapidly dissociates from the CSN upon DNA damage. Here we systematically define the protein interaction network of the mammalian CSN through mass spectrometric interrogation of the CSN subunits Csn1, Csn3, Csn4, Csn5, Csn6 and Csn7a. Notably, we identified a subset of CRL complexes that stably interact with the CSN and thus might similarly be activated by dissociation from the CSN in response to specific cues. In addition, we detected several new proteins in the CRL-CSN interactome, including Dda1, which we characterized as a chromatin-associated core subunit of multiple CRL4 proteins. Cells depleted of Dda1 spontaneously accumulated double-stranded DNA breaks in a similar way to Cul4A-, Cul4B- or Wdr23-depleted cells, indicating that Dda1 interacts physically and functionally with CRL4 complexes. This analysis identifies new components of the CRL family of E3 ligases and elaborates new connections between the CRL and CSN complexes.

Uses of cancer registries for public health and clinical research in Europe : results of the European Network of Cancer Registries survey among 161 population-based cancer registries during 2010-2012

AIM: To provide insight into cancer registration coverage, data access and use in Europe. This contributes to data and infrastructure harmonisation and will foster a more prominent role of cancer registries (CRs) within public health, clinical policy and cancer research, whether within or outside the European Research Area. METHODS: During 2010-12 an extensive survey of cancer registration practices and data use was conducted among 161 population-based CRs across Europe. Responding registries (66%) operated in 33 countries, including 23 with national coverage. RESULTS: Population-based oncological surveillance started during the 1940-50s in the northwest of Europe and from the 1970s to 1990s in other regions. The European Union (EU) protection regulations affected data access, especially in Germany and France, but less in the Netherlands or Belgium. Regular reports were produced by CRs on incidence rates (95%), survival (60%) and stage for selected tumours (80%). Evaluation of cancer control and quality of care remained modest except in a few dedicated CRs. Variables evaluated were support of clinical audits, monitoring adherence to clinical guidelines, improvement of cancer care and evaluation of mass cancer screening. Evaluation of diagnostic imaging tools was only occasional. CONCLUSION: Most population-based CRs are well equipped for strengthening cancer surveillance across Europe. Data quality and intensity of use depend on the role the cancer registry plays in the politico, oncomedical and public health setting within the country. Standard registration methodology could therefore not be translated to equivalent advances in cancer prevention and mass screening, quality of care, translational research of prognosis and survivorship across Europe. Further European collaboration remains essential to ensure access to data and comparability of the results.

Establishment of a Representative Practice-based Research Network (PBRN) for the Monitoring of Primary Care in Switzerland.

Data are urgently needed to better understand processes of care in Swiss primary care (PC). A total of 2027 PC physicians, stratified by canton, were invited to participate in the Swiss Primary care Active Monitoring network, of whom 200 accepted to join. There were no significant differences between participants and a random sample drawn from the same physician databases based on sex, year of obtaining medical school diploma, or location. The Swiss Primary care Active Monitoring network represents the first large-scale, nationally representative practice-based research network in Switzerland and will provide a unique opportunity to better understand the functioning of Swiss PC.

Cell-based computational modeling of vascular morphogenesis using Tissue Simulation Toolkit.

Computational modeling has become a widely used tool for unraveling the mechanisms of higher level cooperative cell behavior during vascular morphogenesis. However, experimenting with published simulation models or adding new assumptions to those models can be daunting for novice and even for experienced computational scientists. Here, we present a step-by-step, practical tutorial for building cell-based simulations of vascular morphogenesis using the Tissue Simulation Toolkit (TST). The TST is a freely available, open-source C++ library for developing simulations with the two-dimensional cellular Potts model, a stochastic, agent-based framework to simulate collective cell behavior. We will show the basic use of the TST to simulate and experiment with published simulations of vascular network formation. Then, we will present step-by-step instructions and explanations for building a recent simulation model of tumor angiogenesis. Demonstrated mechanisms include cell-cell adhesion, chemotaxis, cell elongation, haptotaxis, and haptokinesis.

Hard-wired heterogeneity in blood stem cells revealed using a dynamic regulatory network model.

MOTIVATION: Combinatorial interactions of transcription factors with cis-regulatory elements control the dynamic progression through successive cellular states and thus underpin all metazoan development. The construction of network models of cis-regulatory elements, therefore, has the potential to generate fundamental insights into cellular fate and differentiation. Haematopoiesis has long served as a model system to study mammalian differentiation, yet modelling based on experimentally informed cis-regulatory interactions has so far been restricted to pairs of interacting factors. Here, we have generated a Boolean network model based on detailed cis-regulatory functional data connecting 11 haematopoietic stem/progenitor cell (HSPC) regulator genes. RESULTS: Despite its apparent simplicity, the model exhibits surprisingly complex behaviour that we charted using strongly connected components and shortest-path analysis in its Boolean state space. This analysis of our model predicts that HSPCs display heterogeneous expression patterns and possess many intermediate states that can act as 'stepping stones' for the HSPC to achieve a final differentiated state. Importantly, an external perturbation or 'trigger' is required to exit the stem cell state, with distinct triggers characterizing maturation into the various different lineages. By focusing on intermediate states occurring during erythrocyte differentiation, from our model we predicted a novel negative regulation of Fli1 by Gata1, which we confirmed experimentally thus validating our model. In conclusion, we demonstrate that an advanced mammalian regulatory network model based on experimentally validated cis-regulatory interactions has allowed us to make novel, experimentally testable hypotheses about transcriptional mechanisms that control differentiation of mammalian stem cells. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Paper 3: EUROCAT data quality indicators for population-based registries of congenital anomalies.

The European Surveillance of Congenital Anomalies (EUROCAT) network of population-based congenital anomaly registries is an important source of epidemiologic information on congenital anomalies in Europe covering live births, fetal deaths from 20 weeks gestation, and terminations of pregnancy for fetal anomaly. EUROCAT's policy is to strive for high-quality data, while ensuring consistency and transparency across all member registries. A set of 30 data quality indicators (DQIs) was developed to assess five key elements of data quality: completeness of case ascertainment, accuracy of diagnosis, completeness of information on EUROCAT variables, timeliness of data transmission, and availability of population denominator information. This article describes each of the individual DQIs and presents the output for each registry as well as the EUROCAT (unweighted) average, for 29 full member registries for 2004-2008. This information is also available on the EUROCAT website for previous years. The EUROCAT DQIs allow registries to evaluate their performance in relation to other registries and allows appropriate interpretations to be made of the data collected. The DQIs provide direction for improving data collection and ascertainment, and they allow annual assessment for monitoring continuous improvement. The DQI are constantly reviewed and refined to best document registry procedures and processes regarding data collection, to ensure appropriateness of DQI, and to ensure transparency so that the data collected can make a substantial and useful contribution to epidemiologic research on congenital anomalies.

Atomic force and electron microscopic-based study of sarcolemmal surface of living cardiomyocytes unveils unexpected mitochondrial shift in heart failure.

Loss of T-tubules (TT), sarcolemmal invaginations of cardiomyocytes (CMs), was recently identified as a general heart failure (HF) hallmark. However, whether TT per se or the overall sarcolemma is altered during HF process is still unknown. In this study, we directly examined sarcolemmal surface topography and physical properties using Atomic Force Microscopy (AFM) in living CMs from healthy and failing mice hearts. We confirmed the presence of highly organized crests and hollows along myofilaments in isolated healthy CMs. Sarcolemma topography was tightly correlated with elasticity, with crests stiffer than hollows and related to the presence of few packed subsarcolemmal mitochondria (SSM) as evidenced by electron microscopy. Three days after myocardial infarction (MI), CMs already exhibit an overall sarcolemma disorganization with general loss of crests topography thus becoming smooth and correlating with a decreased elasticity while interfibrillar mitochondria (IFM), myofilaments alignment and TT network were unaltered. End-stage post-ischemic condition (15days post-MI) exacerbates overall sarcolemma disorganization with, in addition to general loss of crest/hollow periodicity, a significant increase of cell surface stiffness. Strikingly, electron microscopy revealed the total depletion of SSM while some IFM heaps could be visualized beneath the membrane. Accordingly, mitochondrial Ca(2+) studies showed a heterogeneous pattern between SSM and IFM in healthy CMs which disappeared in HF. In vitro, formamide-induced sarcolemmal stress on healthy CMs phenocopied post-ischemic kinetics abnormalities and revealed initial SSM death and crest/hollow disorganization followed by IFM later disarray which moved toward the cell surface and structured heaps correlating with TT loss. This study demonstrates that the loss of crest/hollow organization of CM surface in HF occurs early and precedes disruption of the TT network. It also highlights a general stiffness increased of the CM surface most likely related to atypical IFM heaps while SSM died during HF process. Overall, these results indicate that initial sarcolemmal stress leading to SSM death could underlie subsequent TT disarray and HF setting.