Project description:Sea stars and sea urchins evolved from a last common ancestor that lived at the end of the Cambrian, approximately half a billion years ago. In a previous comparative study of the gene regulatory networks (GRNs) that embody the genomic program for embryogenesis in these animals, we discovered an almost perfectly conserved five-gene network subcircuit required for endoderm specification. We show here that the GRN structure upstream and downstream of the conserved network kernel has, by contrast, diverged extensively. Mesoderm specification is accomplished quite differently; the Delta-Notch signaling system is used in radically distinct ways; and various regulatory genes have been coopted to different functions. The conservation of the conserved kernel is thus the more remarkable. The results indicate types of network linkage subject to evolutionary change. An emergent theme is that subcircuit design may be preserved even while the identity of genes performing given roles changes because of alteration in their cis-regulatory control systems.
Project description:The modern evolutionary synthesis seemingly fails to explain how a population can survive a large environmental change: the pre-existence of heritable variants adapted to the novel environment is too opportunistic, whereas the search for new adaptive mutations after the environmental change is so slow that the population may go extinct. Plasticity-led evolution, the initial environmental induction of a novel adaptive phenotype followed by genetic accommodation, has been proposed to solve this problem. However, the mechanism enabling plasticity-led evolution remains unclear. Here, we present computational models that exhibit behaviors compatible with plasticity-led evolution by extending the Wagner model of gene regulatory networks. The models show adaptive plastic response and the uncovering of cryptic mutations under large environmental changes, followed by genetic accommodation. Moreover, these behaviors are consistently observed over distinct novel environments. We further show that environmental cues, developmental processes, and hierarchical regulation cooperatively amplify the above behaviors and accelerate evolution. These observations suggest plasticity-led evolution is a universal property of complex developmental systems independent of particular mutations.
Project description:Convergent phenotypic evolution is often caused by recurrent changes at particular nodes in the underlying gene regulatory networks (GRNs). The genes at such evolutionary 'hotspots' are thought to maximally affect the phenotype with minimal pleiotropic consequences. This has led to the suggestion that if a GRN is understood in sufficient detail, the path of evolution may be predictable. The repeated evolutionary loss of larval trichomes among Drosophila species is caused by the loss of shavenbaby (svb) expression. svb is also required for development of leg trichomes, but the evolutionary gain of trichomes in the 'naked valley' on T2 femurs in Drosophila melanogaster is caused by reduced microRNA-92a (miR-92a) expression rather than changes in svb. We compared the expression and function of components between the larval and leg trichome GRNs to investigate why the genetic basis of trichome pattern evolution differs in these developmental contexts. We found key differences between the two networks in both the genes employed, and in the regulation and function of common genes. These differences in the GRNs reveal why mutations in svb are unlikely to contribute to leg trichome evolution and how instead miR-92a represents the key evolutionary switch in this context. Our work shows that variability in GRNs across different developmental contexts, as well as whether a morphological feature is lost versus gained, influence the nodes at which a GRN evolves to cause morphological change. Therefore, our findings have important implications for understanding the pathways and predictability of evolution.
Project description:Developmental plasticity enables the production of alternative phenotypes in response to different environmental conditions. While significant advances in understanding the ecological and evolutionary implications of plasticity have been made, understanding its genetic basis has lagged. However, a decade of genetic screens in the model nematode Pristionchus pacificus has culminated in 30 genes which affect mouth-form plasticity. We also recently reported the critical window of environmental sensitivity, and therefore have clear expectations for when differential gene expression should matter. Here, we collated previous data into a gene-regulatory network (GRN), and performed developmental transcriptomics across different environmental conditions, genetic backgrounds, and mouth-form mutants to assess the regulatory logic of plasticity. We found that only two genes in the GRN (eud-1 and seud-1/sult-1) are sensitive to the environment during the critical window. Interestingly, the time points of their sensitivity differ, suggesting that they act as sequential checkpoints. We also observed temporal constraint upon the transcriptional effects of mutating the GRN and revealed unexpected feedback between mouth-form genes. Surprisingly, expression of seud-1/sult-1, but not eud-1, correlated with mouth form biases across different strains and species. Finally, a comprehensive analysis of all samples identified metabolism as a shared pathway for regulating mouth-form plasticity. These data are presented in a Shiny app to facilitate gene-expression comparisons across development in up to 14 different conditions. Collectively, our results suggest that mouth-form plasticity evolved a constrained, two-tiered logic to integrate environmental information leading up to the final developmental decision.
Project description:BackgroundSingle-cell RNA sequencing (scRNA-seq) plays a pivotal role in our understanding of cellular heterogeneity. Current analytical workflows are driven by categorizing principles that consider cells as individual entities and classify them into complex taxonomies.ResultsWe devise a conceptually different computational framework based on a holistic view, where single-cell datasets are used to infer global, large-scale regulatory networks. We develop correlation metrics that are specifically tailored to single-cell data, and then generate, validate, and interpret single-cell-derived regulatory networks from organs and perturbed systems, such as diabetes and Alzheimer's disease. Using tools from graph theory, we compute an unbiased quantification of a gene's biological relevance and accurately pinpoint key players in organ function and drivers of diseases.ConclusionsOur approach detects multiple latent regulatory changes that are invisible to single-cell workflows based on clustering or differential expression analysis, significantly broadening the biological insights that can be obtained with this leading technology.
Project description:AimsIndividual prognostic assessment and disease evolution pathways are undefined in chronic heart failure (HF). The application of unsupervised learning methodologies could help to identify patient phenotypes and the progression in each phenotype as well as to assess adverse event risk.Methods and resultsFrom a bulk of 7948 HF patients included in the MECKI registry, we selected patients with a minimum 2-year follow-up. We implemented a topological data analysis (TDA), based on 43 variables derived from clinical, biochemical, cardiac ultrasound, and exercise evaluations, to identify several patients' clusters. Thereafter, we used the trajectory analysis to describe the evolution of HF states, which is able to identify bifurcation points, characterized by different follow-up paths, as well as specific end-stages conditions of the disease. Finally, we conducted a 5-year survival analysis (composite of cardiovascular death, left ventricular assist device, or urgent heart transplant). Findings were validated on internal (n = 527) and external (n = 777) populations. We analyzed 4876 patients (age = 63 [53-71], male gender n = 3973 (81.5%), NYHA class I-II n = 3576 (73.3%), III-IV n = 1300 (26.7%), LVEF = 33 [25.5-39.9], atrial fibrillation n = 791 (16.2%), peak VO2% pred = 54.8 [43.8-67.2]), with a minimum 2-year follow-up. Nineteen patient clusters were identified by TDA. Trajectory analysis revealed a path characterized by 3 bifurcation and 4 end-stage points. Clusters survival rate varied from 44% to 100% at 2 years and from 20% to 100% at 5 years, respectively. The event frequency at 5-year follow-up for each study cohort cluster was successfully compared with those in the validation cohorts (R = 0.94 and R = 0.84, P < 0.001, for internal and external cohort, respectively). Finally, we conducted a 5-year survival analysis (composite of cardiovascular death, left ventricular assist device, or urgent heart transplant observed in 22% of cases).ConclusionsEach HF phenotype has a specific disease progression and prognosis. These findings allow to individualize HF patient evolutions and to tailor assessment.
Project description:Though the bacterial transcription regulation apparatus is distinct in terms of several structural and functional features from its eukaryotic counterpart, the gross structure of the transcription regulatory network (TRN) is believed to be similar in both superkingdoms. Here, we explore the fine structure of the bacterial TRN and the underlying "co-regulatory network" (CRN) to show that despite the superficial similarities to the TRN of the eukaryotic model organism yeast, the bacterial networks display entirely different organizational principles. In particular unlike in eukaryotes, hubs of the bacterial networks are both global regulators and integrators of diverse disparate transcriptional responses. These and other organizational differences might correlate with the fundamental differences in gene and promoter organization in the two superkingdoms, especially the presence of operons and regulons in bacteria. Further we explored to find the interplay, if any, between network structures, mode of regulatory interactions and signal sensing of transcription factors (TFs) in shaping up the bacterial transcriptional regulatory responses. For this purpose, we first classified TFs according to their regulatory mode (activator, repressor or dual regulator) and sensory mechanism (one-component systems responding to internal or external signals, TFs from two-component systems and chromosomal structure modifying TFs) in the bacterial model organism Escherichia coli and then we studied the overall evolutionary optimization of network structures. The incorporation of TFs in different hierarchical elements of the TRN appears to involve on a multi-dimensional selection process depending on regulatory and sensory modes of TFs in motifs, co-regulatory associations between TFs of different functional classes and transcript half-lives. As a result it appears to have generated circuits that allow intricately regulated physiological state changes. We identified the biological significance of most of these optimizations, which can be further used as the basis to explore similar controls in other bacteria. We also show that, though on the larger evolutionary scale, unrelated TFs have evolved to become hubs, within lineages like gamma-proteobacteria there is strong tendency to retain hubs, as well as certain higher-order network modules that have emerged through lineage specific paralog duplications.
Project description:Understanding the complex regulatory networks underlying development and evolution of multi-cellular organisms is a major problem in biology. Computational models can be used as tools to extract the regulatory structure and dynamics of such networks from gene expression data. This approach is called reverse engineering. It has been successfully applied to many gene networks in various biological systems. However, to reconstitute the structure and non-linear dynamics of a developmental gene network in its spatial context remains a considerable challenge. Here, we address this challenge using a case study: the gap gene network involved in segment determination during early development of Drosophila melanogaster. A major problem for reverse-engineering pattern-forming networks is the significant amount of time and effort required to acquire and quantify spatial gene expression data. We have developed a simplified data processing pipeline that considerably increases the throughput of the method, but results in data of reduced accuracy compared to those previously used for gap gene network inference. We demonstrate that we can infer the correct network structure using our reduced data set, and investigate minimal data requirements for successful reverse engineering. Our results show that timing and position of expression domain boundaries are the crucial features for determining regulatory network structure from data, while it is less important to precisely measure expression levels. Based on this, we define minimal data requirements for gap gene network inference. Our results demonstrate the feasibility of reverse-engineering with much reduced experimental effort. This enables more widespread use of the method in different developmental contexts and organisms. Such systematic application of data-driven models to real-world networks has enormous potential. Only the quantitative investigation of a large number of developmental gene regulatory networks will allow us to discover whether there are rules or regularities governing development and evolution of complex multi-cellular organisms.
Project description:Recent developments in measurement techniques have enabled us to observe the time series of many components simultaneously. Thus, it is important to understand not only the dynamics of individual time series but also their interactions. Although there are many methods for analysing the interaction between two or more time series, there are very few methods that describe global changes of the interactions over time. Here, we propose an approach to visualise time evolution for the global changes of the interactions in complex systems. This approach consists of two steps. In the first step, we construct a meta-time series of networks. In the second step, we analyse and visualise this meta-time series by using distance and recurrence plots. Our two-step approach involving intermediate network representations elucidates the half-a-day periodicity of foreign exchange markets and a singular functional network in the brain related to perceptual alternations.
Project description:We analyze the effect of spike-timing-dependent plasticity (STDP) on a system of pulse-coupled class I neurons. Our research begins with a system of two mutually connected quadratic integrate-and-fire (QIF) neurons, which are canonical representatives of class I neurons. Along with various asymptotic modes previously observed in other neuronal models with plastic synapses, we found a stable synchronous mode characterized by unidirectional link from a slower neuron to a faster neuron. In this frequency-locked mode, the faster neuron emits multiple spikes per cycle of the slower neuron. We analytically obtain the Arnold tongues for this mode without STDP and with STDP. We also consider larger plastic networks of QIF neurons and show that the detected mode can manifest itself in such a way that slow neurons become pacemakers. As a result, slow and fast neurons can form large synchronous clusters that generate low-frequency oscillations. We demonstrate the generality of the results obtained with two connected QIF neurons using Wang-Buzsáki and Morris-Lecar biophysically plausible class I neuron models.