Evolution of biological cooperation: an algorithmic approach.

This manuscript presents an algorithmic approach to cooperation in biological systems, drawing on fundamental ideas from statistical mechanics and probability theory. Fisher's geometric model of adaptation suggests that the evolution of organisms well adapted to multiple constraints comes at a significant complexity cost. By utilizing combinatorial models of fitness, we demonstrate that the probability of adapting to all constraints decreases exponentially with the number of constraints, thereby generalizing Fisher's result.

Global repression by tailless during segmentation.

The orphan nuclear receptor Tailless (Tll) exhibits conserved roles in brain formation and maintenance that are shared, for example, with vertebrate orthologous forms (Tlx). However, the early expression of tll in two gap domains in the segmentation cascade of Drosophila is unusual even for most other insects. Here we investigate tll regulation on pair-rule stripes.

The effect of mutational robustness on the evolvability of multicellular organisms and eukaryotic cells.

Canalization involves mutational robustness, the lack of phenotypic change as a result of genetic mutations. Given the large divergence in phenotype across species, understanding the relationship between high robustness and evolvability has been of interest to both theorists and experimentalists. Although canalization was originally proposed in the context of multicellular organisms, the effect of multicellularity and other classes of hierarchical organization on evolvability has not been considered by theoreticians.

Fully interpretable deep learning model of transcriptional control.

Motivation: The universal expressibility assumption of Deep Neural Networks (DNNs) is the key motivation behind recent worksin the systems biology community to employDNNs to solve important problems in functional genomics and moleculargenetics. Typically, such investigations have taken a 'black box' approach in which the internal structure of themodel used is set purely by machine learning considerations with little consideration of representing the internalstructure of the biological system by the mathematical structure of the DNN. DNNs have not yet been applied to thedetailed modeling of transcriptional control in which mRNA production is controlled by the binding of specific transcriptionfactors to DNA, in part because such models are in part formulated in terms of specific chemical equationsthat appear different in form from those used in neural networks.

An in silico analysis of robust but fragile gene regulation links enhancer length to robustness.

Organisms must ensure that expression of genes is directed to the appropriate tissues at the correct times, while simultaneously ensuring that these gene regulatory systems are robust to perturbation. This idea is captured by a mathematical concept called r-robustness, which says that a system is robust to a perturbation in up to r - 1 randomly chosen parameters. r-robustness implies that the biological system has a small number of sensitive parameters and that this number can be used as a robustness measure.

Adaptation, fitness landscape learning and fast evolution.

We consider evolution of a large population, where fitness of each organism is defined by many phenotypical traits. These traits result from expression of many genes. Under some assumptions on  fitness we prove that such model organisms  are capable, to some extent, to recognize the fitness landscape.

Physical implications of so(2, 1) symmetry in exact solutions for a self-repressing gene.

We chemically characterize the symmetries underlying the exact solutions of a stochastic negatively self-regulating gene. The breaking of symmetry at a low molecular number causes three effects. Two branches of the solution exist, having high and low switching rates, such that the low switching rate branch approaches deterministic behavior and the high switching rate branch exhibits sub-Fano behavior.

Correction: A sequence level model of an intact locus predicts the location and function of nonadditive enhancers.

[This corrects the article DOI: 10.1371/journal.pone.

Synthetic enhancer design by in silico compensatory evolution reveals flexibility and constraint in cis-regulation.

Background: Models that incorporate specific chemical mechanisms have been successful in describing the activity of Drosophila developmental enhancers as a function of underlying transcription factor binding motifs. Despite this, the minimum set of mechanisms required to reconstruct an enhancer from its constituent parts is not known. Synthetic biology offers the potential to test the sufficiency of known mechanisms to describe the activity of enhancers, as well as to uncover constraints on the number, order, and spacing of motifs.

A sequence level model of an intact locus predicts the location and function of nonadditive enhancers.

Metazoan gene expression is controlled through the action of long stretches of noncoding DNA that contain enhancers-shorter sequences responsible for controlling a single aspect of a gene's expression pattern. Models built on thermodynamics have shown how enhancers interpret protein concentration in order to determine specific levels of gene expression, but the emergent regulatory logic of a complete regulatory locus shows qualitative and quantitative differences from isolated enhancers. Such differences may arise from steric competition limiting the quantity of DNA that can simultaneously influence the transcription machinery.

DNA Occupancy of Polymerizing Transcription Factors: A Chemical Model of the ETS Family Factor Yan.

Transcription factors use both protein-DNA and protein-protein interactions to assemble appropriate complexes to regulate gene expression. Although most transcription factors operate as monomers or dimers, a few, including the E26 transformation-specific family repressors Drosophila melanogaster Yan and its human homolog TEL/ETV6, can polymerize. Although polymerization is required for both the normal and oncogenic function of Yan and TEL/ETV6, the mechanisms by which it influences the recruitment, organization, and stability of transcriptional complexes remain poorly understood.

The analysis of novel distal Cebpa enhancers and silencers using a transcriptional model reveals the complex regulatory logic of hematopoietic lineage specification.

C/EBPα plays an instructive role in the macrophage-neutrophil cell-fate decision and its expression is necessary for neutrophil development. How Cebpa itself is regulated in the myeloid lineage is not known. We decoded the cis-regulatory logic of Cebpa, and two other myeloid transcription factors, Egr1 and Egr2, using a combined experimental-computational approach.

Parallel Simulated Annealing Using an Adaptive Resampling Interval.

This paper presents a parallel simulated annealing algorithm that is able to achieve 90% parallel efficiency in iteration on up to 192 processors and up to 40% parallel efficiency in time when applied to a 5000-dimension Rastrigin function. Our algorithm breaks scalability barriers in the method of Chu (1999) by abandoning adaptive cooling based on variance. The resulting gains in parallel efficiency are much larger than the loss of serial efficiency from lack of adaptive cooling.

Noise-Driven Phenotypic Heterogeneity with Finite Correlation Time in Clonal Populations.

There has been increasing awareness in the wider biological community of the role of clonal phenotypic heterogeneity in playing key roles in phenomena such as cellular bet-hedging and decision making, as in the case of the phage-λ lysis/lysogeny and B. Subtilis competence/vegetative pathways. Here, we report on the effect of stochasticity in growth rate, cellular memory/intermittency, and its relation to phenotypic heterogeneity.

Natural variation of the expression pattern of the segmentation gene even-skipped in melanogaster.

The evolution of canalized traits is a central question in evolutionary biology. Natural variation in highly conserved traits can provide clues about their evolutionary potential. Here we investigate natural variation in a conserved trait-even-skipped (eve) expression at the cellular blastoderm stage of embryonic development in Drosophila melanogaster.

Gene regulation and noise reduction by coupling of stochastic processes.

Here we characterize the low-noise regime of a stochastic model for a negative self-regulating binary gene. The model has two stochastic variables, the protein number and the state of the gene. Each state of the gene behaves as a protein source governed by a Poisson process.

Ancestral resurrection of the Drosophila S2E enhancer reveals accessible evolutionary paths through compensatory change.

Upstream regulatory sequences that control gene expression evolve rapidly, yet the expression patterns and functions of most genes are typically conserved. To address this paradox, we have reconstructed computationally and resurrected in vivo the cis-regulatory regions of the ancestral Drosophila eve stripe 2 element and evaluated its evolution using a mathematical model of promoter function. Our feed-forward transcriptional model predicts gene expression patterns directly from enhancer sequence.

A prognostic gene signature for metastasis-free survival of triple negative breast cancer patients.

Although triple negative breast cancers (TNBC) are the most aggressive subtype of breast cancer, they currently lack targeted therapies. Because this classification still includes a heterogeneous collection of tumors, new tools to classify TNBCs are urgently required in order to improve our prognostic capability for high risk patients and predict response to therapy. We previously defined a gene expression signature, RKIP Pathway Metastasis Signature (RPMS), based upon a metastasis-suppressive signaling pathway initiated by Raf Kinase Inhibitory Protein (RKIP).

Quantitative imaging of gene expression in Drosophila embryos.

Quantitative measurements derived using sophisticated microscopy techniques are essential for understanding the basic principles that control the behavior of biological systems. Here we describe a data pipeline developed to extract quantitative data on segmentation gene expression from confocal images of gene expression patterns in Drosophila. The pipeline consists of image segmentation, background removal, temporal characterization of an embryo, data registration, and data averaging.

Preparation of Drosophila embryos for quantitative imaging of gene expression.

Quantitative measurements derived using sophisticated microscopy techniques are essential for understanding the basic principles that control the behavior of biological systems. We have developed a five-step data pipeline to extract quantitative data on segmentation gene expression from confocal images of gene expression patterns in Drosophila. This protocol describes the preparation of Drosophila embryos for imaging by confocal microscopy.

A synthetic biology approach to the development of transcriptional regulatory models and custom enhancer design.

Synthetic biology offers novel opportunities for elucidating transcriptional regulatory mechanisms and enhancer logic. Complex cis-regulatory sequences--like the ones driving expression of the Drosophila even-skipped gene--have proven difficult to design from existing knowledge, presumably due to the large number of protein-protein interactions needed to drive the correct expression patterns of genes in multicellular organisms. This work discusses two novel computational methods for the custom design of enhancers that employ a sophisticated, empirically validated transcriptional model, optimization algorithms, and synthetic biology.

Rearrangements of 2.5 kilobases of noncoding DNA from the Drosophila even-skipped locus define predictive rules of genomic cis-regulatory logic.

Rearrangements of about 2.5 kilobases of regulatory DNA located 5' of the transcription start site of the Drosophila even-skipped locus generate large-scale changes in the expression of even-skipped stripes 2, 3, and 7. The most radical effects are generated by juxtaposing the minimal stripe enhancers MSE2 and MSE3 for stripes 2 and 3 with and without small "spacer" segments less than 360 bp in length.

Quantitative dynamics and increased variability of segmentation gene expression in the Drosophila Krüppel and knirps mutants.

Here we characterize the response of the Drosophila segmentation system to mutations in two gap genes, Kr and kni, in the form of single or double homozygotes and single heterozygotes. Segmentation gene expression in these genotypes was quantitatively monitored with cellular resolution in space and 6.5 to 13min resolution in time.

Drosophila blastoderm patterning.

The Drosophila blastoderm embryo is a classic model for the study of the genetics of pattern formation. In recent years, quantitative empirical approaches have been employed extensively in the study of blastoderm pattern formation. This quantitative work has enabled the development of a number of data-driven computational models.

Modeling of gap gene expression in Drosophila Kruppel mutants.

The segmentation gene network in Drosophila embryo solves the fundamental problem of embryonic patterning: how to establish a periodic pattern of gene expression, which determines both the positions and the identities of body segments. The gap gene network constitutes the first zygotic regulatory tier in this process. Here we have applied the systems-level approach to investigate the regulatory effect of gap gene Kruppel (Kr) on segmentation gene expression.