Predictability of evolutionary trajectories in fitness landscapes

Experimental studies on enzyme evolution show that only a small fraction of all possible mutation trajectories are accessible to evolution. However, these experiments deal with individual enzymes and explore a tiny part of the fitness landscape. We report an exhaustive analysis of fitness landscapes constructed with an off-lattice model of protein folding where fitness is equated with robustness to misfolding. This model mimics the essential features of the interactions between amino acids, is consistent with the key paradigms of protein folding and reproduces the universal distribution of evolutionary rates among orthologous proteins. We introduce mean path divergence as a quantitative measure of the degree to which the starting and ending points determine the path of evolution in fitness landscapes. Global measures of landscape roughness are good predictors of path divergence in all studied landscapes: the mean path divergence is greater in smooth landscapes than in rough ones. The model-derived and experimental landscapes are significantly smoother than random landscapes and resemble additive landscapes perturbed with moderate amounts of noise; thus, these landscapes are substantially robust to mutation. The model landscapes show a deficit of suboptimal peaks even compared with noisy additive landscapes with similar overall roughness. We suggest that smoothness and the substantial deficit of peaks in the fitness landscapes of protein evolution are fundamental consequences of the physics of protein folding.Comment: 14 pages, 7 figure

Reciprocal Sign Epistasis between Frequently Experimentally Evolved Adaptive Mutations Causes a Rugged Fitness Landscape

The fitness landscape captures the relationship between genotype and evolutionary fitness and is a pervasive metaphor used to describe the possible evolutionary trajectories of adaptation. However, little is known about the actual shape of fitness landscapes, including whether valleys of low fitness create local fitness optima, acting as barriers to adaptive change. Here we provide evidence of a rugged molecular fitness landscape arising during an evolution experiment in an asexual population of Saccharomyces cerevisiae. We identify the mutations that arose during the evolution using whole-genome sequencing and use competitive fitness assays to describe the mutations individually responsible for adaptation. In addition, we find that a fitness valley between two adaptive mutations in the genes MTH1 and HXT6/HXT7 is caused by reciprocal sign epistasis, where the fitness cost of the double mutant prohibits the two mutations from being selected in the same genetic background. The constraint enforced by reciprocal sign epistasis causes the mutations to remain mutually exclusive during the experiment, even though adaptive mutations in these two genes occur several times in independent lineages during the experiment. Our results show that epistasis plays a key role during adaptation and that inter-genic interactions can act as barriers between adaptive solutions. These results also provide a new interpretation on the classic Dobzhansky-Muller model of reproductive isolation and display some surprising parallels with mutations in genes often associated with tumors

EQUIP: Implementing chronic care principles and applying formative evaluation methods to improve care for schizophrenia: QUERI Series

<p>Abstract</p> <p>Background</p> <p>This paper presents a case study that demonstrates the evolution of a project entitled "Enhancing QUality-of-care In Psychosis" (EQUIP) that began approximately when the U.S. Department of Veterans Affairs' Quality Enhancement Research Initiative (QUERI), and implementation science were emerging. EQUIP developed methods and tools to implement chronic illness care principles in the treatment of schizophrenia, and evaluated this implementation using a small-scale controlled trial. The next iteration of the project, EQUIP-2, was further informed by implementation science and the use of QUERI tools.</p> <p>Methods</p> <p>This paper reports the background, development, results and implications of EQUIP, and also describes ongoing work in the second phase of the project (EQUIP-2). The EQUIP intervention uses implementation strategies and tools to increase the adoption and implementation of chronic illness care principles. In EQUIP-2, these strategies and tools are conceptually grounded in a stages-of-change model, and include clinical and delivery system interventions and adoption/implementation tools. Formative evaluation occurs in conjunction with the intervention, and includes developmental, progress-focused, implementation-focused, and interpretive evaluation.</p> <p>Results</p> <p>Evaluation of EQUIP provided an understanding of quality gaps <it>and </it>how to address related problems in schizophrenia. EQUIP showed that solutions to quality problems in schizophrenia differ by treatment domain and are exacerbated by a lack of awareness of evidence-based practices. EQUIP also showed that improving care requires creating resources for physicians to help them easily implement practice changes, plus intensive education as well as product champions who help physicians use these resources. Organizational changes, such as the addition of care managers and informatics systems, were shown to help physicians with identifying problems, making referrals, and monitoring follow-up. In EQUIP-2, which is currently in progress, these initial findings were used to develop a more comprehensive approach to implementing and evaluating the chronic illness care model.</p> <p>Discussion</p> <p>In QUERI, small-scale projects contribute to the development and enhancement of hands-on, action-oriented service-directed projects that are grounded in current implementation science. This project supports the concept that QUERI tools can be useful in implementing complex care models oriented toward evidence-based improvement of clinical care.</p

Emergence of terpene cyclization in Artemisia annua

The emergence of terpene cyclization was critical to the evolutionary expansion of chemical diversity yet remains unexplored. Here we report the first discovery of an epistatic network of residues that controls the onset of terpene cyclization in Artemisia annua. We begin with amorpha-4,11-diene synthase (ADS) and (E)-b-farnesene synthase (BFS), a pair of terpene synthases that produce cyclic or linear terpenes, respectively. A library of B27,000 enzymes is generated by breeding combinations of natural amino-acid substitutions from the cyclic into the linear producer. We discover one dominant mutation is sufficient to activate cyclization, and together with two additional residues comprise a network of strongly epistatic interactions that activate, suppress or reactivate cyclization. Remarkably, this epistatic network of equivalent residues also controls cyclization in a BFS homologue from Citrus junos. Fitness landscape analysis of mutational trajectories provides quantitative insights into a major epoch in specialized metabolism

Heterogeneous Adaptive Trajectories of Small Populations on Complex Fitness Landscapes

Background Small populations are thought to be adaptively handicapped, not only because they suffer more from deleterious mutations but also because they have limited access to new beneficial mutations, particularly those conferring large benefits. Methodology/Principal Findings Here, we test this widely held conjecture using both simulations and experiments with small and large bacterial populations evolving in either a simple or a complex nutrient environment. Consistent with expectations, we find that small populations are adaptively constrained in the simple environment; however, in the complex environment small populations not only follow more heterogeneous adaptive trajectories, but can also attain higher fitness than the large populations. Large populations are constrained to near deterministic fixation of rare large-benefit mutations. While such determinism speeds adaptation on the smooth adaptive landscape represented by the simple environment, it can limit the ability of large populations from effectively exploring the underlying topography of rugged adaptive landscapes characterized by complex environments. Conclusions Our results show that adaptive constraints often faced by small populations can be circumvented during evolution on rugged adaptive landscapes

Systematic Dissection and Trajectory-Scanning Mutagenesis of the Molecular Interface That Ensures Specificity of Two-Component Signaling Pathways

Two-component signal transduction systems enable bacteria to sense and respond to a wide range of environmental stimuli. Sensor histidine kinases transmit signals to their cognate response regulators via phosphorylation. The faithful transmission of information through two-component pathways and the avoidance of unwanted cross-talk require exquisite specificity of histidine kinase-response regulator interactions to ensure that cells mount the appropriate response to external signals. To identify putative specificity-determining residues, we have analyzed amino acid coevolution in two-component proteins and identified a set of residues that can be used to rationally rewire a model signaling pathway, EnvZ-OmpR. To explore how a relatively small set of residues can dictate partner selectivity, we combined alanine-scanning mutagenesis with an approach we call trajectory-scanning mutagenesis, in which all mutational intermediates between the specificity residues of EnvZ and another kinase, RstB, were systematically examined for phosphotransfer specificity. The same approach was used for the response regulators OmpR and RstA. Collectively, the results begin to reveal the molecular mechanism by which a small set of amino acids enables an individual kinase to discriminate amongst a large set of highly-related response regulators and vice versa. Our results also suggest that the mutational trajectories taken by two-component signaling proteins following gene or pathway duplication may be constrained and subject to differential selective pressures. Only some trajectories allow both the maintenance of phosphotransfer and the avoidance of unwanted cross-talk

Exploring the Complexity of the HIV-1 Fitness Landscape

Although fitness landscapes are central to evolutionary theory, so far no biologically realistic examples for large-scale fitness landscapes have been described. Most currently available biological examples are restricted to very few loci or alleles and therefore do not capture the high dimensionality characteristic of real fitness landscapes. Here we analyze large-scale fitness landscapes that are based on predictive models for in vitro replicative fitness of HIV-1. We find that these landscapes are characterized by large correlation lengths, considerable neutrality, and high ruggedness and that these properties depend only weakly on whether fitness is measured in the absence or presence of different antiretrovirals. Accordingly, adaptive processes on these landscapes depend sensitively on the initial conditions. While the relative extent to which mutations affect fitness on their own (main effects) or in combination with other mutations (epistasis) is a strong determinant of these properties, the fitness landscape of HIV-1 is considerably less rugged, less neutral, and more correlated than expected from the distribution of main effects and epistatic interactions alone. Overall this study confirms theoretical conjectures about the complexity of biological fitness landscapes and the importance of the high dimensionality of the genetic space in which adaptation takes place

Modeling Evolutionary Dynamics of Epigenetic Mutations in Hierarchically Organized Tumors

The cancer stem cell (CSC) concept is a highly debated topic in cancer research. While experimental evidence in favor of the cancer stem cell theory is apparently abundant, the results are often criticized as being difficult to interpret. An important reason for this is that most experimental data that support this model rely on transplantation studies. In this study we use a novel cellular Potts model to elucidate the dynamics of established malignancies that are driven by a small subset of CSCs. Our results demonstrate that epigenetic mutations that occur during mitosis display highly altered dynamics in CSC-driven malignancies compared to a classical, non-hierarchical model of growth. In particular, the heterogeneity observed in CSC-driven tumors is considerably higher. We speculate that this feature could be used in combination with epigenetic (methylation) sequencing studies of human malignancies to prove or refute the CSC hypothesis in established tumors without the need for transplantation. Moreover our tumor growth simulations indicate that CSC-driven tumors display evolutionary features that can be considered beneficial during tumor progression. Besides an increased heterogeneity they also exhibit properties that allow the escape of clones from local fitness peaks. This leads to more aggressive phenotypes in the long run and makes the neoplasm more adaptable to stringent selective forces such as cancer treatment. Indeed when therapy is applied the clone landscape of the regrown tumor is more aggressive with respect to the primary tumor, whereas the classical model demonstrated similar patterns before and after therapy. Understanding these often counter-intuitive fundamental properties of (non-)hierarchically organized malignancies is a crucial step in validating the CSC concept as well as providing insight into the therapeutical consequences of this model

Break dosage, cell cycle stage and DNA replication influence DNA double strand break response

DNA double strand breaks (DSBs) can be repaired by non-homologous end joining (NHEJ) or homology-directed repair (HR). HR requires nucleolytic degradation of 5′ DNA ends to generate tracts of single-stranded DNA (ssDNA), which are also important for the activation of DNA damage checkpoints. Here we describe a quantitative analysis of DSB processing in the budding yeast Saccharomyces cerevisiae. We show that resection of an HO endonuclease-induced DSB is less extensive than previously estimated and provide evidence for significant instability of the 3′ ssDNA tails. We show that both DSB resection and checkpoint activation are dose-dependent, especially during the G1 phase of the cell cycle. During G1, processing near the break is inhibited by competition with NHEJ, but extensive resection is regulated by an NHEJ-independent mechanism. DSB processing and checkpoint activation are more efficient in G2/M than in G1 phase, but are most efficient at breaks encountered by DNA replication forks during S phase. Our findings identify unexpected complexity of DSB processing and its regulation, and provide a framework for further mechanistic insights

Colonic epithelial ion transport is not affected in patients with diverticulosis

<p>Abstract</p> <p>Background</p> <p>Colonic diverticular disease is a bothersome condition with an unresolved pathogenesis. It is unknown whether a neuroepithelial dysfunction is present. The aim of the study was two-fold; (1) to investigate colonic epithelial ion transport in patients with diverticulosis and (2) to adapt a miniaturized Modified Ussing Air-Suction (MUAS) chamber for colonic endoscopic biopsies.</p> <p>Methods</p> <p>Biopsies were obtained from the sigmoid part of the colon. 86 patients were included. All patients were referred for colonoscopy on suspicion of neoplasia and they were without pathological findings at colonoscopy (controls) except for diverticulosis in 22 (D-patients). Biopsies were mounted in MUAS chambers with an exposed area of 5 mm². Electrical responses to various stimulators and inhibitors of ion transport were investigated together with histological examination. The MUAS chamber was easy to use and reproducible data were obtained.</p> <p>Results</p> <p>Median basal short circuit current (SCC) was 43.8 μA·cm^-2(0.8 – 199) for controls and 59.3 μA·cm^-2(3.0 – 177.2) for D-patients. Slope conductance was 77.0 mS·cm^-2(18.6 – 204.0) equal to 13 Ω·cm²for controls and 96.6 mS·cm^-2(8.4 – 191.4) equal to 10.3 Ω·cm²for D-patients. Stimulation with serotonin, theophylline, forskolin and carbachol induced increases in SCC in a range of 4.9 – 18.6 μA·cm^-2, while inhibition with indomethacin, bumetanide, ouabain and amiloride decreased SCC in a range of 6.5 – 27.4 μA·cm^-2, and all with no significant differences between controls and D-patients. Histological examinations showed intact epithelium and lamina propria before and after mounting for both types of patients.</p> <p>Conclusion</p> <p>We conclude that epithelial ion transport is not significantly altered in patients with diverticulosis and that the MUAS chamber can be adapted for studies of human colonic endoscopic biopsies.</p