Accelerated search for biomolecular network models to interpret high-throughput experimental data

<p>Abstract</p> <p>Background</p> <p>The functions of human cells are carried out by biomolecular networks, which include proteins, genes, and regulatory sites within DNA that encode and control protein expression. Models of biomolecular network structure and dynamics can be inferred from high-throughput measurements of gene and protein expression. We build on our previously developed fuzzy logic method for bridging quantitative and qualitative biological data to address the challenges of noisy, low resolution high-throughput measurements, i.e., from gene expression microarrays. We employ an evolutionary search algorithm to accelerate the search for hypothetical fuzzy biomolecular network models consistent with a biological data set. We also develop a method to estimate the probability of a potential network model fitting a set of data by chance. The resulting metric provides an estimate of both model quality and dataset quality, identifying data that are too noisy to identify meaningful correlations between the measured variables.</p> <p>Results</p> <p>Optimal parameters for the evolutionary search were identified based on artificial data, and the algorithm showed scalable and consistent performance for as many as 150 variables. The method was tested on previously published human cell cycle gene expression microarray data sets. The evolutionary search method was found to converge to the results of exhaustive search. The randomized evolutionary search was able to converge on a set of similar best-fitting network models on different training data sets after 30 generations running 30 models per generation. Consistent results were found regardless of which of the published data sets were used to train or verify the quantitative predictions of the best-fitting models for cell cycle gene dynamics.</p> <p>Conclusion</p> <p>Our results demonstrate the capability of scalable evolutionary search for fuzzy network models to address the problem of inferring models based on complex, noisy biomolecular data sets. This approach yields multiple alternative models that are consistent with the data, yielding a constrained set of hypotheses that can be used to optimally design subsequent experiments.</p

Tempo and mode of performance evolution across multiple independent origins of adhesive toe pads in lizards

Understanding macroevolutionary dynamics of trait evolution is an important endeavor in evolutionary biology. Ecological opportunity can liberate a trait as it diversifies through trait space, while genetic and selective constraints can limit diversification. While many studies have examined the dynamics of morphological traits, diverse morphological traits may yield the same or similar performance and as performance is often more proximately the target of selection, examining only morphology may give an incomplete understanding of evolutionary dynamics. Here, we ask whether convergent evolution of pad‐bearing lizards has followed similar evolutionary dynamics, or whether independent origins are accompanied by unique constraints and selective pressures over macroevolutionary time. We hypothesized that geckos and anoles each have unique evolutionary tempos and modes. Using performance data from 59 species, we modified Brownian motion (BM) and Ornstein–Uhlenbeck (OU) models to account for repeated origins estimated using Bayesian ancestral state reconstructions. We discovered that adhesive performance in geckos evolved in a fashion consistent with Brownian motion with a trend, whereas anoles evolved in bounded performance space consistent with more constrained evolution (an Ornstein–Uhlenbeck model). Our results suggest that convergent phenotypes can have quite distinctive evolutionary patterns, likely as a result of idiosyncratic constraints or ecological opportunities

Non-equilibrium dynamics and floral trait interactions shape extant angiosperm diversity.

Why are some traits and trait combinations exceptionally common across the tree of life, whereas others are vanishingly rare? The distribution of trait diversity across a clade at any time depends on the ancestral state of the clade, the rate at which new phenotypes evolve, the differences in speciation and extinction rates across lineages, and whether an equilibrium has been reached. Here we examine the role of transition rates, differential diversification (speciation minus extinction) and non-equilibrium dynamics on the evolutionary history of angiosperms, a clade well known for the abundance of some trait combinations and the rarity of others. Our analysis reveals that three character states (corolla present, bilateral symmetry, reduced stamen number) act synergistically as a key innovation, doubling diversification rates for lineages in which this combination occurs. However, this combination is currently less common than predicted at equilibrium because the individual characters evolve infrequently. Simulations suggest that angiosperms will remain far from the equilibrium frequencies of character states well into the future. Such non-equilibrium dynamics may be common when major innovations evolve rarely, allowing lineages with ancestral forms to persist, and even outnumber those with diversification-enhancing states, for tens of millions of years

Evolutionary history of four binary blue stragglers from the globular clusters \omega Cen, M55, 47 Tuc and NGC 6752

Context. Origin and evolution of blue stragglers in globular clusters is still a matter of debate. Aims. The aim of the present investigation is to reproduce the evolutionary history of four binary blue stragglers in four different clusters, for which precise values of global parameters are known. Methods. Using the model for cool close binary evolution, developed by one of us (KS), progenitors of all investigated binaries were found and their parameters evolved into the presently observed values. Results. The results show that the progenitors of the binary blue stragglers are cool close binaries with period of a few days, which transform into stragglers by rejuvenation of the initially less massive component by mass transfer from its more massive companion overflowing the inner critical Roche surface. The parameters of V209 from \omega Cen indicate that the binary is substantially enriched in helium. This is an independent and strong evidence for the existence of the helium rich subpopulation in this cluster.Comment: Accepted by Astronomy and Astrophysic