Feed-forward regulation adaptively evolves via dynamics rather than topology when there is intrinsic noise.

In transcriptional regulatory networks (TRNs), a canonical 3-node feed-forward loop (FFL) is hypothesized to evolve to filter out short spurious signals. We test this adaptive hypothesis against a novel null evolutionary model. Our mutational model captures the intrinsically high prevalence of weak affinity transcription factor binding sites. We also capture stochasticity and delays in gene expression that distort external signals and intrinsically generate noise. Functional FFLs evolve readily under selection for the hypothesized function but not in negative controls. Interestingly, a 4-node "diamond" motif also emerges as a short spurious signal filter. The diamond uses expression dynamics rather than path length to provide fast and slow pathways. When there is no idealized external spurious signal to filter out, but only internally generated noise, only the diamond and not the FFL evolves. While our results support the adaptive hypothesis, we also show that non-adaptive factors, including the intrinsic expression dynamics, matter.

Inhibiting mitochondrial phosphate transport as an unexploited antifungal strategy.

The development of effective antifungal therapeutics remains a formidable challenge because of the close evolutionary relationship between humans and fungi. Mitochondrial function may present an exploitable vulnerability because of its differential utilization in fungi and its pivotal roles in fungal morphogenesis, virulence, and drug resistance already demonstrated by others. We now report mechanistic characterization of ML316, a thiohydantoin that kills drug-resistant Candida species at nanomolar concentrations through fungal-selective inhibition of the mitochondrial phosphate carrier Mir1. Using genetic, biochemical, and metabolomic approaches, we established ML316 as the first Mir1 inhibitor. Inhibition of Mir1 by ML316 in respiring yeast diminished mitochondrial oxygen consumption, resulting in an unusual metabolic catastrophe marked by citrate accumulation and death. In a mouse model of azole-resistant oropharyngeal candidiasis, ML316 reduced fungal burden and enhanced azole activity. Targeting Mir1 could provide a new, much-needed therapeutic strategy to address the rapidly rising burden of drug-resistant fungal infection.

A Fungal-Selective Cytochrome bc1 Inhibitor Impairs Virulence and Prevents the Evolution of Drug Resistance.

To cause disease, a microbial pathogen must adapt to the challenges of its host environment. The leading fungal pathogen Candida albicans colonizes nutrient-poor bodily niches, withstands attack from the immune system, and tolerates treatment with azole antifungals, often evolving resistance. To discover agents that block these adaptive strategies, we screened 300,000 compounds for inhibition of azole tolerance in a drug-resistant Candida isolate. We identified a novel indazole derivative that converts azoles from fungistatic to fungicidal drugs by selective inhibition of mitochondrial cytochrome bc1. We synthesized 103 analogs to optimize potency (half maximal inhibitory concentration 0.4 ?M) and fungal selectivity (28-fold over human). In addition to reducing azole resistance, targeting cytochrome bc1 prevents C. albicans from adapting to the nutrient-deprived macrophage phagosome and greatly curtails its virulence in mice. Inhibiting mitochondrial respiration and restricting metabolic flexibility with this synthetically tractable chemotype provides an attractive therapeutic strategy to limit both fungal virulence and drug resistance.

Scalable and cost-effective NGS genotyping in the cloud.

BACKGROUND: While next-generation sequencing (NGS) costs have plummeted in recent years, cost and complexity of computation remain substantial barriers to the use of NGS in routine clinical care. The clinical potential of NGS will not be realized until robust and routine whole genome sequencing data can be accurately rendered to medically actionable reports within a time window of hours and at scales of economy in the 10's of dollars. RESULTS: We take a step towards addressing this challenge, by using COSMOS, a cloud-enabled workflow management system, to develop GenomeKey, an NGS whole genome analysis workflow. COSMOS implements complex workflows making optimal use of high-performance compute clusters. Here we show that the Amazon Web Service (AWS) implementation of GenomeKey via COSMOS provides a fast, scalable, and cost-effective analysis of both public benchmarking and large-scale heterogeneous clinical NGS datasets. CONCLUSIONS: Our systematic benchmarking reveals important new insights and considerations to produce clinical turn-around of whole genome analysis optimization and workflow management including strategic batching of individual genomes and efficient cluster resource configuration.

An evolutionarily conserved prion-like element converts wild fungi from metabolic specialists to generalists.

[GAR(+)] is a protein-based element of inheritance that allows yeast (Saccharomyces cerevisiae) to circumvent a hallmark of their biology: extreme metabolic specialization for glucose fermentation. When glucose is present, yeast will not use other carbon sources. [GAR(+)] allows cells to circumvent this "glucose repression." [GAR(+)] is induced in yeast by a factor secreted by bacteria inhabiting their environment. We report that de novo rates of [GAR(+)] appearance correlate with the yeast's ecological niche. Evolutionarily distant fungi possess similar epigenetic elements that are also induced by bacteria. As expected for a mechanism whose adaptive value originates from the selective pressures of life in biological communities, the ability of bacteria to induce [GAR(+)] and the ability of yeast to respond to bacterial signals have been extinguished repeatedly during the extended monoculture of domestication. Thus, [GAR(+)] is a broadly conserved adaptive strategy that links environmental and social cues to heritable changes in metabolism.

Cross-kingdom chemical communication drives a heritable, mutually beneficial prion-based transformation of metabolism.

In experimental science, organisms are usually studied in isolation, but in the wild, they compete and cooperate in complex communities. We report a system for cross-kingdom communication by which bacteria heritably transform yeast metabolism. An ancient biological circuit blocks yeast from using other carbon sources in the presence of glucose. [GAR(+)], a protein-based epigenetic element, allows yeast to circumvent this "glucose repression" and use multiple carbon sources in the presence of glucose. Some bacteria secrete a chemical factor that induces [GAR(+)]. [GAR(+)] is advantageous to bacteria because yeast cells make less ethanol and is advantageous to yeast because their growth and long-term viability is improved in complex carbon sources. This cross-kingdom communication is broadly conserved, providing a compelling argument for its adaptive value. By heritably transforming growth and survival strategies in response to the selective pressures of life in a biological community, [GAR(+)] presents a unique example of Lamarckian inheritance.

COSMOS: Python library for massively parallel workflows.

SUMMARY: Efficient workflows to shepherd clinically generated genomic data through the multiple stages of a next-generation sequencing pipeline are of critical importance in translational biomedical science. Here we present COSMOS, a Python library for workflow management that allows formal description of pipelines and partitioning of jobs. In addition, it includes a user interface for tracking the progress of jobs, abstraction of the queuing system and fine-grained control over the workflow. Workflows can be created on traditional computing clusters as well as cloud-based services. AVAILABILITY AND IMPLEMENTATION: Source code is available for academic non-commercial research purposes. Links to code and documentation are provided at http://lpm.hms.harvard.edu and http://wall-lab.stanford.edu. CONTACT: dpwall@stanford.edu or peter_tonellato@hms.harvard.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

PLAAC: a web and command-line application to identify proteins with prion-like amino acid composition.

UNLABELLED: Prions are self-templating protein aggregates that stably perpetuate distinct biological states and are of keen interest to researchers in both evolutionary and biomedical science. The best understood prions are from yeast and have a prion-forming domain with strongly biased amino acid composition, most notably enriched for Q or N. PLAAC is a web application that scans protein sequences for domains with P: rion- L: ike A: mino A: cid C: omposition. Users can upload sequence files, or paste sequences directly into a textbox. PLAAC ranks the input sequences by several summary scores and allows scores along sequences to be visualized. Text output files can be downloaded for further analyses, and visualizations saved in PDF and PNG formats. AVAILABILITY AND IMPLEMENTATION: http://plaac.wi.mit.edu/. The Ruby-based web framework and the command-line software (implemented in Java, with visualization routines in R) are available at http://github.com/whitehead/plaac under the MIT license. All software can be run under OS X, Windows and Unix.

Yeast reveal a "druggable" Rsp5/Nedd4 network that ameliorates α-synuclein toxicity in neurons.

α-Synuclein (α-syn) is a small lipid-binding protein implicated in several neurodegenerative diseases, including Parkinson's disease, whose pathobiology is conserved from yeast to man. There are no therapies targeting these underlying cellular pathologies, or indeed those of any major neurodegenerative disease. Using unbiased phenotypic screens as an alternative to target-based approaches, we discovered an N-aryl benzimidazole (NAB) that strongly and selectively protected diverse cell types from α-syn toxicity. Three chemical genetic screens in wild-type yeast cells established that NAB promoted endosomal transport events dependent on the E3 ubiquitin ligase Rsp5/Nedd4. These same steps were perturbed by α-syn itself. Thus, NAB identifies a druggable node in the biology of α-syn that can correct multiple aspects of its underlying pathology, including dysfunctional endosomal and endoplasmic reticulum-to-Golgi vesicle trafficking.

Heritable remodeling of yeast multicellularity by an environmentally responsive prion.

Prion proteins undergo self-sustaining conformational conversions that heritably alter their activities. Many of these proteins operate at pivotal positions in determining how genotype is translated into phenotype. But the breadth of prion influences on biology and their evolutionary significance are just beginning to be explored. We report that a prion formed by the Mot3 transcription factor, [MOT3(+)], governs the acquisition of facultative multicellularity in the budding yeast Saccharomyces cerevisiae. The traits governed by [MOT3(+)] involved both gains and losses of Mot3 regulatory activity. [MOT3(+)]-dependent expression of FLO11, a major determinant of cell-cell adhesion, produced diverse lineage-specific multicellular phenotypes in response to nutrient deprivation. The prions themselves were induced by ethanol and eliminated by hypoxia-conditions that occur sequentially in the natural respiro-fermentative cycles of yeast populations. These data demonstrate that prions can act as environmentally responsive molecular determinants of multicellularity and contribute to the natural morphological diversity of budding yeast.

Inhibiting GPI anchor biosynthesis in fungi stresses the endoplasmic reticulum and enhances immunogenicity.

In fungi, the anchoring of proteins to the plasma membrane via their covalent attachment to glycosylphosphatidylinositol (GPI) is essential and thus provides a valuable point of attack for the development of antifungal therapeutics. Unfortunately, studying the underlying biology of GPI-anchor synthesis is difficult, especially in medically relevant fungal pathogens because they are not genetically tractable. Compounding difficulties, many of the genes in this pathway are essential in Saccharomyces cerevisiae. Here, we report the discovery of a new small molecule christened gepinacin (for GPI acylation inhibitor) which selectively inhibits Gwt1, a critical acyltransferase required for the biosynthesis of fungal GPI anchors. After delineating the target specificity of gepinacin using genetic and biochemical techniques, we used it to probe key, therapeutically relevant consequences of disrupting GPI anchor metabolism in fungi. We found that, unlike all three major classes of antifungals in current use, the direct antimicrobial activity of this compound results predominantly from its ability to induce overwhelming stress to the endoplasmic reticulum. Gepinacin did not affect the viability of mammalian cells nor did it inhibit their orthologous acyltransferase. This enabled its use in co-culture experiments to examine Gwt1's effects on host-pathogen interactions. In isolates of Candida albicans, the most common fungal pathogen in humans, exposure to gepinacin at sublethal concentrations impaired filamentation and unmasked cell wall ß-glucan to stimulate a pro-inflammatory cytokine response in macrophages. Gwt1 is a promising antifungal drug target, and gepanacin is a useful probe for studying how disrupting GPI-anchor synthesis impairs viability and alters host-pathogen interactions in genetically intractable fungi.

Prions are a common mechanism for phenotypic inheritance in wild yeasts.

The self-templating conformations of yeast prion proteins act as epigenetic elements of inheritance. Yeast prions might provide a mechanism for generating heritable phenotypic diversity that promotes survival in fluctuating environments and the evolution of new traits. However, this hypothesis is highly controversial. Prions that create new traits have not been found in wild strains, leading to the perception that they are rare 'diseases' of laboratory cultivation. Here we biochemically test approximately 700 wild strains of Saccharomyces for [PSI(+)] or [MOT3(+)], and find these prions in many. They conferred diverse phenotypes that were frequently beneficial under selective conditions. Simple meiotic re-assortment of the variation harboured within a strain readily fixed one such trait, making it robust and prion-independent. Finally, we genetically screened for unknown prion elements. Fully one-third of wild strains harboured them. These, too, created diverse, often beneficial phenotypes. Thus, prions broadly govern heritable traits in nature, in a manner that could profoundly expand adaptive opportunities.

The spontaneous appearance rate of the yeast prion [PSI+] and its implications for the evolution of the evolvability properties of the [PSI+] system.

Epigenetically inherited aggregates of the yeast prion [PSI+] cause genomewide readthrough translation that sometimes increases evolvability in certain harsh environments. The effects of natural selection on modifiers of [PSI+] appearance have been the subject of much debate. It seems likely that [PSI+] would be at least mildly deleterious in most environments, but this may be counteracted by its evolvability properties on rare occasions. Indirect selection on modifiers of [PSI+] is predicted to depend primarily on the spontaneous [PSI+] appearance rate, but this critical parameter has not previously been adequately measured. Here we measure this epimutation rate accurately and precisely as 5.8 x 10(-7) per generation, using a fluctuation test. We also determine that genetic "mimics" of [PSI+] account for up to 80% of all phenotypes involving general nonsense suppression. Using previously developed mathematical models, we can now infer that even in the absence of opportunities for adaptation, modifiers of [PSI+] are only weakly deleterious relative to genetic drift. If we assume that the spontaneous [PSI+] appearance rate is at its evolutionary optimum, then opportunities for adaptation are inferred to be rare, such that the [PSI+] system is favored only very weakly overall. But when we account for the observed increase in the [PSI+] appearance rate in response to stress, we infer much higher overall selection in favor of [PSI+] modifiers, suggesting that [PSI+]-forming ability may be a consequence of selection for evolvability.

The evolution of reversible switches in the presence of irreversible mimics.

Reversible phenotypic switching can be caused by a number of different mechanisms including epigenetic inheritance systems and DNA-based contingency loci. Previous work has shown that reversible switching systems may be favored by natural selection. Many switches can be characterized as "on/off" where the "off" state constitutes a temporary and reversible loss of function. Loss-of-function phenotypes corresponding to the "off" state can be produced in many different ways, all yielding identical fitness in the short term. In the long term, however, a switch-induced loss of function can be reversed, whereas many loss-of-function mutations, especially deletions, cannot. We refer to these loss-of-function mutations as "irreversible mimics" of the reversible switch. Here, we develop a model in which a reversible switch evolves in the presence of both irreversible mimics and metapopulation structure. We calculate that when the rate of appearance of irreversible mimics exceeds the migration rate, the evolved reversible switching rate will exceed the bet-hedging rate predicted by panmictic models.

Balancing selection and heterogeneity across the classical human leukocyte antigen loci: a meta-analytic review of 497 population studies.

This paper presents a meta-analysis of high-resolution human leukocyte antigen (HLA) allele frequency data describing 497 population samples. Most of the datasets were compiled from studies published in eight journals from 1990 to 2007; additional datasets came from the International Histocompatibility Workshops and from the AlleleFrequencies.net database. In all, these data represent approximately 66,800 individuals from throughout the world, providing an opportunity to observe trends that may not have been evident at the time the data were originally analyzed, especially with regard to the relative importance of balancing selection among the HLA loci. Population genetic measures of allele frequency distributions were summarized across populations by locus and geographic region. A role for balancing selection maintaining much of HLA variation was confirmed. Further, the breadth of this meta-analysis allowed the ranking of the HLA loci, with DQA1 and HLA-C showing the strongest balancing selection and DPB1 being compatible with neutrality. Comparisons of the allelic spectra reported by studies since 1990 indicate that most of the HLA alleles identified since 2000 are very-low-frequency alleles. The literature-based allele-count data, as well as maps summarizing the geographic distributions for each allele, are available online.