Permalink
Cannot retrieve contributors at this time
Fetching contributors…
Cannot retrieve contributors at this time
| <?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.1d1 20130915//EN" "JATS-archivearticle1.dtd"><article article-type="research-article" dtd-version="1.1d1" xmlns:xlink="http://www.w3.org/1999/xlink"><front><journal-meta><journal-id journal-id-type="nlm-ta">eLife</journal-id><journal-id journal-id-type="hwp">elife</journal-id><journal-id journal-id-type="publisher-id">eLife</journal-id><journal-title-group><journal-title>eLife</journal-title></journal-title-group><issn publication-format="electronic">2050-084X</issn><publisher><publisher-name>eLife Sciences Publications, Ltd</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="publisher-id">00367</article-id><article-id pub-id-type="doi">10.7554/eLife.00367</article-id><article-categories><subj-group subj-group-type="display-channel"><subject>Research article</subject></subj-group><subj-group subj-group-type="heading"><subject>Genomics and evolutionary biology</subject></subj-group></article-categories><title-group><article-title>Improved use of a public good selects for the evolution of undifferentiated multicellularity</article-title></title-group><contrib-group><contrib contrib-type="author" corresp="yes" id="author-2995"><name><surname>Koschwanez</surname><given-names>John H</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="corresp" rid="cor1">*</xref><xref ref-type="other" rid="par-1"/><xref ref-type="other" rid="par-2"/><xref ref-type="fn" rid="con1"/><xref ref-type="fn" rid="conf1"/><xref ref-type="other" rid="dataro1"/><xref ref-type="other" rid="dataro2"/></contrib><contrib contrib-type="author" id="author-2624"><name><surname>Foster</surname><given-names>Kevin R</given-names></name><xref ref-type="aff" rid="aff2"/><xref ref-type="other" rid="par-2"/><xref ref-type="fn" rid="con2"/><xref ref-type="fn" rid="conf1"/><xref ref-type="other" rid="dataro1"/><xref ref-type="other" rid="dataro2"/></contrib><contrib contrib-type="author" id="author-3077"><name><surname>Murray</surname><given-names>Andrew W</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="other" rid="par-2"/><xref ref-type="fn" rid="con3"/><xref ref-type="fn" rid="conf1"/><xref ref-type="other" rid="dataro1"/><xref ref-type="other" rid="dataro2"/></contrib><aff id="aff1"><institution content-type="dept">FAS Center for Systems Biology and Department of Molecular and Cellular Biology</institution>, <institution>Harvard University</institution>, <addr-line><named-content content-type="city">Cambridge</named-content></addr-line>, <country>United States</country></aff><aff id="aff2"><institution content-type="dept">Department of Zoology</institution>, <institution>University of Oxford</institution>, <addr-line><named-content content-type="city">Oxford</named-content></addr-line>, <country>United Kingdom</country></aff></contrib-group><contrib-group content-type="section"><contrib contrib-type="editor"><name><surname>Tautz</surname><given-names>Diethard</given-names></name><role>Reviewing editor</role></contrib><aff><institution>Max Planck Institute for Evolutionary Biology</institution>, <country>Germany</country></aff></contrib-group><author-notes><corresp id="cor1"><label>*</label>For correspondence: <email>jkoschwanez@cgr.harvard.edu</email></corresp></author-notes><pub-date date-type="pub" publication-format="electronic"><day>02</day><month>04</month><year>2013</year></pub-date><pub-date pub-type="collection"><year>2013</year></pub-date><volume>2</volume><elocation-id>e00367</elocation-id><history><date date-type="received"><day>01</day><month>11</month><year>2012</year></date><date date-type="accepted"><day>14</day><month>02</month><year>2013</year></date></history><permissions><copyright-statement>© 2013, Koschwanez et al</copyright-statement><copyright-year>2013</copyright-year><copyright-holder>Koschwanez et al</copyright-holder><license xlink:href="http://creativecommons.org/licenses/by/3.0/"><license-p>This article is distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/3.0/">Creative Commons Attribution License</ext-link>, which permits unrestricted use and redistribution provided that the original author and source are credited.</license-p></license></permissions><self-uri content-type="pdf" xlink:href="elife00367.pdf"/><related-article ext-link-type="doi" id="ra1" related-article-type="commentary" xlink:href="10.7554/eLife.00655"/><abstract><object-id pub-id-type="doi">10.7554/eLife.00367.001</object-id><p>We do not know how or why multicellularity evolved. We used the budding yeast, <italic>Saccharomyces cerevisiae</italic>, to ask whether nutrients that must be digested extracellularly select for the evolution of undifferentiated multicellularity. Because yeast use invertase to hydrolyze sucrose extracellularly and import the resulting monosaccharides, single cells cannot grow at low cell and sucrose concentrations. Three engineered strategies overcame this problem: forming multicellular clumps, importing sucrose before hydrolysis, and increasing invertase expression. We evolved populations in low sucrose to ask which strategy they would adopt. Of 12 successful clones, 11 formed multicellular clumps through incomplete cell separation, 10 increased invertase expression, none imported sucrose, and 11 increased hexose transporter expression, a strategy we had not engineered. Identifying causal mutations revealed genes and pathways, which frequently contributed to the evolved phenotype. Our study shows that combining rational design with experimental evolution can help evaluate hypotheses about evolutionary strategies.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.001">http://dx.doi.org/10.7554/eLife.00367.001</ext-link></p></abstract><abstract abstract-type="executive-summary"><object-id pub-id-type="doi">10.7554/eLife.00367.002</object-id><title>eLife digest</title><p>Life first appeared on Earth more than 3 billion years ago in the form of single-celled microorganisms. The diverse array of complex life forms that we see today evolved from these humble beginnings, but it is not clear what triggered the evolution of multicellular organisms from single cells.</p><p>One of the simplest multicellular eukaryotes is the yeast, <italic>Saccharomyces cerevisiae</italic>—a fungus that has been used for centuries in baking and brewing and, more recently, as a model organism in molecular biology. Yeast cells feed on sugar (sucrose), but are unable to absorb it directly from their surroundings. Instead they secrete an enzyme called invertase, which breaks down the sucrose into simpler components that cells can take up with the help of sugar transporters.</p><p>However, single yeast cells living in a low-sucrose environment face a problem: most of the simple sugars that they produce diffuse out of reach. To overcome this difficulty, the cells could form multicellular clumps, which would enable each cell to consume the sugars that drift away from its neighbours. Alternatively, the cells could increase their production of invertase, or they could begin to take up sucrose directly.</p><p>Using genetic engineering, Koschwanez et al. produced three strains of yeast, each with one of these traits, and confirmed that all three strategies do indeed help fungi to grow in low sucrose. But could any of these traits evolve spontaneously? To test this possibility, Koschwanez et al. introduced wild-type yeast cells into a low-sucrose environment and studied any populations of cells that managed to survive. Of 12 that did, 11 had acquired the ability to form multicellular clumps, while 10 had increased their expression of invertase. Surprisingly, none had evolved the ability to import sucrose. However, 11 of the populations that survived also displayed an adaptation that the researchers had not predicted beforehand: they all expressed higher levels of the sugar transporters that take up sucrose breakdown products.</p><p>The work of Koschwanez et al. suggests that the benefits of being able to share invertase and, therefore, simple sugars, may have driven the evolution of multicellularity in ancient organisms. Moreover, their use of rational design (engineered mutations) combined with experimental evolution (allowing colonies to grow under selection pressure and studying the strategies that they adopt) offers a new approach to studying evolution in the lab.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.002">http://dx.doi.org/10.7554/eLife.00367.002</ext-link></p></abstract><kwd-group kwd-group-type="author-keywords"><title>Author keywords</title><kwd>Multicellularity</kwd><kwd>Experimental evolution</kwd><kwd>Evolution of cooperation</kwd></kwd-group><kwd-group kwd-group-type="research-organism"><title>Research organism</title><kwd>S. cerevisiae</kwd></kwd-group><funding-group><award-group id="par-1"><funding-source><institution-wrap><institution>National Institute of General Medical Science</institution></institution-wrap></funding-source><award-id>K25GM085806</award-id><principal-award-recipient><name><surname>Koschwanez</surname><given-names>John H</given-names></name></principal-award-recipient></award-group><award-group id="par-2"><funding-source><institution-wrap><institution>National Institute of General Medical Science</institution></institution-wrap></funding-source><award-id>GM06873</award-id><principal-award-recipient><name><surname>Koschwanez</surname><given-names>John H</given-names></name><name><surname>Foster</surname><given-names>Kevin R</given-names></name><name><surname>Murray</surname><given-names>Andrew W</given-names></name></principal-award-recipient></award-group><funding-statement>The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.</funding-statement></funding-group><custom-meta-group><custom-meta specific-use="meta-only"><meta-name>Author impact statement</meta-name><meta-value>The ability to share resources for the benefit of all members of a group may have driven ancient organisms to evolve from a unicellular to a multicellular state.</meta-value></custom-meta></custom-meta-group></article-meta></front><body><sec id="s1" sec-type="intro"><title>Introduction</title><p>Multicellular organisms have evolved from a unicellular ancestor at least 25 times (<xref ref-type="bibr" rid="bib28">Grosberg and Strathmann, 2007</xref>), but we know little about what selected for the simplest form of multicellularity: an undifferentiated clump of cells produced by the repeated division of a single cell. Two driving forces have been proposed, protection from a variety of factors (including predation [<xref ref-type="bibr" rid="bib35">Kessin et al., 1996</xref>], environmental stress [<xref ref-type="bibr" rid="bib68">Smukalla et al., 2008</xref>], and phagocytosis [<xref ref-type="bibr" rid="bib9">Boraas et al., 1998</xref>]) and more efficient nutrient usage (<xref ref-type="bibr" rid="bib23">Dworkin, 1972</xref>; <xref ref-type="bibr" rid="bib61">Pfeiffer and Bonhoeffer, 2003</xref>; <xref ref-type="bibr" rid="bib37">Koschwanez et al., 2011</xref>; <xref ref-type="bibr" rid="bib3">Alegado et al., 2012</xref>).</p><p>In earlier work, we showed that sharing public goods favors clumps over isolated cells and proposed that sharing could have selected for simple multicellularity (<xref ref-type="bibr" rid="bib37">Koschwanez et al., 2011</xref>). The budding yeast, <italic>Saccharomyces cerevisiae</italic>, utilizes sucrose by secreting invertase (<xref ref-type="bibr" rid="bib20">Dodyk and Rothstein, 1964</xref>; <xref ref-type="bibr" rid="bib15">Carlson et al., 1981</xref>). Over 95% of this enzyme remains in the cell wall (<xref ref-type="bibr" rid="bib25">Esmon et al., 1987</xref>; <xref ref-type="bibr" rid="bib70">Tammi et al., 1987</xref>), where it hydrolyzes sucrose into glucose and fructose, which are imported into the cell by a variety of hexose transporters (<xref ref-type="bibr" rid="bib52">Meijer et al., 1996</xref>; <xref ref-type="bibr" rid="bib63">Reifenberger et al., 1997</xref>). Lab yeast strains cannot grow from low density in low concentrations of sucrose because of diffusion: each cell captures only a small fraction of the sugars that sucrose hydrolysis releases, and the molecules released by other, distant cells are at very low concentration. As a result, cells cannot capture enough glucose and fructose to grow. Forming multicellular clumps overcomes this failure; cells in a clump can capture glucose and fructose diffusing from their neighbors and grow in concentrations of sucrose where low concentrations of individual cells cannot.</p><p>Speculating on evolution based on experiments with engineered strains is problematic. How well is the ease of engineering a strategy correlated with its evolutionary accessibility? Are multiple mutations required? Do these mutations reduce fitness in other environmental conditions? Are other strategies more accessible? Do certain combinations of strategies outcompete single strategies? And finally, how many different strategies does a set of parallel cultures adopt? Experimentally evolving populations, characterizing their phenotypes, and finding the mutations responsible for evolution allows us to address these questions.</p><p>We compared engineered and evolved solutions to the problem of growing on low sucrose concentrations. We engineered and tested three strategies that allow yeast cells to grow on low sucrose: forming multicellular clumps, importing sucrose before hydrolyzing it, and increasing invertase expression. We also evolved unengineered, laboratory strains to grow in the same conditions. All but one of the 12 clones selected from 10 independent populations form multicellular clumps as a result of incomplete cell separation, showing that this selection efficiently selects for multicellularity. In addition, 10 of the clones elevated invertase expression and 11 elevated hexose transporter expression, a strategy that we failed to anticipate, but none showed evidence of sucrose import. We combined bulk segregant analysis (<xref ref-type="bibr" rid="bib53">Michelmore et al., 1991</xref>; <xref ref-type="bibr" rid="bib10">Brauer et al., 2006</xref>; <xref ref-type="bibr" rid="bib67">Segrè et al., 2006</xref>; <xref ref-type="bibr" rid="bib7">Birkeland et al., 2010</xref>; <xref ref-type="bibr" rid="bib48">Magwene et al., 2011</xref>) with whole genome sequencing to identify putative causal mutations. We recreated two of the evolved clones, one with five mutations and one with eight mutations, to show that mutations we had identified were indeed causal. Finally, we competed the evolved clones against their ancestor and found that adaptation in sucrose severely reduces fitness in high glucose.</p></sec><sec id="s2" sec-type="results"><title>Results</title><sec id="s2-1"><title>Three engineered strategies for growth in low sucrose</title><p>We began by asking whether undifferentiated multicellularity was the only strategy that allowed yeast cells to grow on low sucrose. We hypothesized that there were two alternative strategies: increasing invertase expression, and importing sucrose and then hydrolyzing it inside the cell (<xref ref-type="fig" rid="fig1">Figure 1</xref>). Having previously tested undifferentiated multicellularity (<xref ref-type="bibr" rid="bib37">Koschwanez et al., 2011</xref>), we engineered the other two strategies and tested whether they allow cells to grow from low densities in low sucrose concentrations. Population growth requires both cell growth and cell proliferation. For simplicity, we refer to the combination of these properties as growth.<fig id="fig1" position="float"><object-id pub-id-type="doi">10.7554/eLife.00367.003</object-id><label>Figure 1.</label><caption><title>Three engineered strategies for growth in low sucrose.</title><p>Strategy 1, form multicellular clumps, was previously verified (<xref ref-type="bibr" rid="bib37">Koschwanez et al., 2011</xref>). The results of testing strategy 2, make more invertase, and strategy 3, import sucrose, are shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>. All three strategies outcompete wild-type strains when the sole carbon source is 1 mM sucrose (<xref ref-type="table" rid="tbl1">Table 1</xref>).</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.003">http://dx.doi.org/10.7554/eLife.00367.003</ext-link></p></caption><graphic xlink:href="elife00367f001"/></fig></p><p>Increased invertase expression will increase the rate of sucrose hydrolysis at the cell wall. Although cells will still lose most of the fructose and glucose to the environment, the increased hydrolysis rate will increase the monosaccharide concentration at the plasma membrane and lead to a higher rate of sugar import. We increased invertase expression by replacing the promoter of the invertase gene (<italic>SUC2</italic>) with the <italic>GAL1</italic> promoter in a yeast strain that is unable to utilize galactose (see <xref ref-type="supplementary-material" rid="SD3-data">supplementary file 3</xref> for all strains used in this study; <xref ref-type="bibr" rid="bib30">Ingolia and Murray 2007</xref>). As a result, galactose serves as a gratuitous inducer: it induces Suc2 but cannot itself be metabolized. Single cells were placed in each microwell of a 96-well plate, and incubated with a range of sucrose and galactose concentrations. <xref ref-type="fig" rid="fig2">Figure 2A</xref> shows that increasing invertase expression allows growth from a single cell in low sucrose concentrations; increasing inducer concentrations leads to better growth.<fig id="fig2" position="float"><object-id pub-id-type="doi">10.7554/eLife.00367.004</object-id><label>Figure 2.</label><caption><title>Two strategies for growth from low sucrose concentrations.</title><p>(<bold>A</bold>) Strong expression of secreted invertase allows growth from a single cell at low sucrose concentrations. All <italic>GAL1</italic> promoter induction data is from the same yeast strain yJHK312 in which transcription of <italic>SUC2</italic> is driven by the <italic>GAL1</italic> promoter. Galactokinase (<italic>GAL1</italic>) is deleted from this strain so that galactose acts as an inducer and not as a carbon source, and the Gal regulon has been engineered to produce a graded rather than a bistable response to increased galactose concentrations by overexpressing <italic>GAL3</italic> from the <italic>ACT1</italic> promoter (<xref ref-type="bibr" rid="bib30">Ingolia and Murray, 2007</xref>). (<bold>B</bold>) Sucrose import allows growth from a single cell in low sucrose concentrations. The ‘<italic>SUC2</italic>, import' strain yJHK372 expresses <italic>SUC2</italic> from the <italic>SUC2</italic> promoter and <italic>MAL11</italic> from the <italic>ACT1</italic> promoter. The ‘<italic>SUC2</italic>, no import' strain yJHK222 expresses <italic>SUC2</italic> from the <italic>SUC2</italic> promoter. The ‘<italic>suc2-1cyt</italic>, import' strain yJHK373 expresses cytoplasmic invertase from the <italic>SUC2</italic> promoter and <italic>MAL11</italic> from the <italic>ACT1</italic> promoter. The ‘<italic>suc2</italic>Δ, import' strain yJHK374 has <italic>SUC2</italic> deleted and expresses <italic>MAL11</italic> from the <italic>ACT1</italic> promoter. For both (<bold>A</bold>) and (<bold>B</bold>), single cells were inoculated by fluorescence activated cell sorting (FACS) into 150 µl wells at the given sugar and galactose concentration and grown without shaking for 85 hr at 30°C and the results shown are totals of three experiments; each experiment used one plate per sugar concentration, and each plate used 24 wells per strain or galactose concentration. In both figures, 2 mM glucose + 2 mM fructose is used as a positive control, and error bars refer to 95% binomial confidence interval using the adjusted Wald method. FRU is fructose, GLC is glucose, and SUC is sucrose.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.004">http://dx.doi.org/10.7554/eLife.00367.004</ext-link></p></caption><graphic xlink:href="elife00367f002"/></fig></p><p>An alternative strategy is importing sucrose before hydrolyzing it since some invertase molecules are retained in the cytoplasm rather than being exported by protein secretion. Mal11 is a maltose importer that can also import sucrose (<xref ref-type="bibr" rid="bib69">Stambuk et al., 1999</xref>), but is not expressed in our yeast strains (<xref ref-type="bibr" rid="bib12">Brown et al., 2010</xref>). We expressed Mal11 from a strong, constitutive promoter (<italic>P</italic><sub><italic>ACT1</italic></sub>) in three different strains: a standard lab strain (<italic>SUC2)</italic>, a strain that does not secrete invertase but produces cytoplasmic invertase (<italic>suc2-1cyt</italic>) (<xref ref-type="bibr" rid="bib37">Koschwanez et al., 2011</xref>), and a strain that lacks invertase (<italic>suc2</italic>Δ). <xref ref-type="fig" rid="fig2">Figure 2B</xref> shows that both the <italic>SUC2</italic> and <italic>suc2-1cyt</italic> strains that make the maltose importer can grow from single cells on as little as 1 mM sucrose. The nearly identical growth of the <italic>SUC2</italic> and <italic>suc2-1cyt</italic> strains shows that the importer makes extracellular sucrose digestion dispensable, and the failure of the <italic>suc2Δ</italic> strains to grow shows that sucrose utilization after import still requires invertase. Thus three different strategies, clumping, increased invertase expression, and sucrose import, each allows yeast to grow from low cell density at low sucrose concentrations.</p><p>We competed each of the three strategies against an unmodified strain to ask whether they could invade an ancestral population. Derivatives of the two competing strains, each expressing a different fluorescent protein, were mixed together and passaged on 1 mM sucrose. In each passage, the cells were grown together for almost eight generations and then diluted 200-fold into fresh medium. The competition was assessed by the number of passages required to eliminate the less fit strain, or the ratio between the two strains at the end of the sixth passage. The three engineered strains all outcompeted the ancestral wild-type lab strain in 1 mM sucrose (<xref ref-type="table" rid="tbl1">Table 1</xref>), although the strain engineered to express increased invertase (EngHiInvertase) was a much worse competitor than the other two (EngClumpy, EngSucImport).<table-wrap id="tbl1" position="float"><object-id pub-id-type="doi">10.7554/eLife.00367.005</object-id><label>Table 1.</label><caption><p>Fitness of engineered strains, evolved clones, recreated strains, and reverted strains</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.005">http://dx.doi.org/10.7554/eLife.00367.005</ext-link></p></caption><table frame="hsides" rules="groups"><thead><tr><td>Strain 1</td><td>Strain 2</td><td>1 mM sucrose</td><td>1 mM glucose + 1 mM fructose</td><td>80 mM glucose</td></tr></thead><tbody><tr><td>EngClumpy</td><td>Wild-type lab</td><td>+++</td><td>0</td><td>0</td></tr><tr><td>EngHiInvertase</td><td>Wild-type lab</td><td>+</td><td>0</td><td>0</td></tr><tr><td>EngSucImport</td><td>Wild-type lab</td><td>+++</td><td>0</td><td>0</td></tr><tr><td>EvoClone1</td><td>Ancestor</td><td>++++</td><td>−−−</td><td>−−−−</td></tr><tr><td>EvoClone2</td><td>Ancestor</td><td>++++</td><td>0</td><td>−−</td></tr><tr><td>EvoClone3</td><td>Ancestor</td><td>++++</td><td>−−−</td><td>−−−−</td></tr><tr><td>EvoClone4</td><td>Ancestor</td><td>++++</td><td>+</td><td>−−−</td></tr><tr><td>EvoClone5</td><td>Ancestor</td><td>++++</td><td>+</td><td>−−−</td></tr><tr><td>EvoClone6</td><td>Ancestor</td><td>++++</td><td>+</td><td>−−</td></tr><tr><td>EvoClone7A</td><td>Ancestor</td><td>++++</td><td>−−</td><td>−−−</td></tr><tr><td>EvoClone7B</td><td>Ancestor</td><td>++++</td><td>+</td><td>−</td></tr><tr><td>EvoClone7C</td><td>Ancestor</td><td>++++</td><td>−−−</td><td>−−−−</td></tr><tr><td>EvoClone8</td><td>Ancestor</td><td>++++</td><td>−</td><td>−−−−</td></tr><tr><td>EvoClone9</td><td>Ancestor</td><td>++++</td><td>−</td><td>−−−</td></tr><tr><td>EvoClone10</td><td>Ancestor</td><td>++++</td><td>−</td><td>−−−−</td></tr><tr><td>Recreated2</td><td>Ancestor</td><td>++++</td><td>−</td><td>−</td></tr><tr><td>Recreated9</td><td>Ancestor</td><td>++++</td><td>0</td><td>−</td></tr><tr><td>EvoClone2</td><td>Recreated2</td><td>+</td><td>+</td><td>0</td></tr><tr><td>EvoClone9</td><td>Recreated9</td><td>+</td><td>−</td><td>−−−</td></tr><tr><td>EvoClone2</td><td>Reverted2</td><td>++++</td><td>+</td><td>−</td></tr><tr><td>EvoClone9</td><td>Reverted9</td><td>++++</td><td>+</td><td>+</td></tr><tr><td>EvoClone2</td><td>EngClumpy</td><td>++</td><td>0</td><td>−−</td></tr><tr><td><italic>ace2Δ</italic></td><td>Ancestor</td><td>+++</td><td>0</td><td>0</td></tr><tr><td><italic>gin4-W19* irc8-G57V mck1-G227Vfs249</italic></td><td>Ancestor</td><td>+++</td><td>−</td><td>−</td></tr></tbody></table><table-wrap-foot><fn><p>++++Strain 1 eliminates strain 2 in 1–2 growth cycles.</p></fn><fn><p>+++ Strain 1 eliminates strain 2 in 3–4 growth cycles.</p></fn><fn><p>++ Strain 1 eliminates strain 2 in 5–6 growth cycles.</p></fn><fn><p>+ Strain 1 > 75% of population after 6 growth cycles.</p></fn><fn><p>0 Neither strain is >75% of population after 6 growth cycles.</p></fn><fn><p>− Strain 2 > 75% of population after 6 growth cycles.</p></fn><fn><p>−− Strain 2 eliminates strain 1 in 5–6 growth cycles.</p></fn><fn><p>−−− Strain 2 eliminates strain 1 in 3–4 growth cycles.</p></fn><fn><p>−−−− Strain 2 eliminates strain 1 in 1–2 growth cycles.</p></fn><fn><p>Growth cycle numbers are averaged over three independent experiments.</p></fn></table-wrap-foot></table-wrap></p></sec><sec id="s2-2"><title>Experimental evolution produces multicellularity</title><p>Having shown that engineering can produce three different strategies for growth on low sucrose, we asked, ‘What would evolution do?' We experimentally evolved multiple, parallel cultures to grow well on low sucrose. Ten independent populations of budding yeast were serially diluted in minimal, 1 mM sucrose-containing medium. The starting strain for each population was a haploid, non-clumpy, prototrophic strain (<xref ref-type="fig" rid="fig3">Figure 3A</xref>) that constitutively expressed YFP and carried a DNA polymerase mutation (<italic>POL3-L523D</italic>) that elevated its mutation rate roughly 100-fold (<xref ref-type="bibr" rid="bib32">Jin et al., 2005</xref>). We used a mutator strain to increase the speed of adaptation. At each growth cycle, 5×10<sup>5</sup> cells were inoculated into 50 ml of medium (<xref ref-type="fig" rid="fig3">Figure 3B</xref>). Over 25–35 serial dilutions, the time it took the culture to become cloudy fell from 2 weeks to 3 days, and all 10 populations (named EvoPopulation1–10) displayed a clumpy phenotype (<xref ref-type="fig" rid="fig3">Figure 3C</xref>). There was no loss of constitutive YFP expression.<fig-group><fig id="fig3" position="float"><object-id pub-id-type="doi">10.7554/eLife.00367.006</object-id><label>Figure 3.</label><caption><title>Evolved populations show a clumpy phenotype.</title><p>(<bold>A</bold>) An ancestor derivative (yJHK111) after growth in 1 mM sucrose. (<bold>B</bold>) Schematic of experimental evolution. Cells were inoculated in 1 mM sucrose media, grown to high density, and then 10<sup>5</sup> cells were reinoculated into fresh media for a total of 25–35 cycles. A sample was frozen down at each serial dilution. (<bold>C</bold>) Samples taken from the last time point of the evolved populations. Representative DIC images were taken with a 40× objective in a glass-bottomed, 96-well plate. All scale bars are 50 µm.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.006">http://dx.doi.org/10.7554/eLife.00367.006</ext-link></p></caption><graphic xlink:href="elife00367f003"/></fig><fig id="fig3fs1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.00367.007</object-id><label>Figure 3—figure supplement 1.</label><caption><title>Eleven of twelve clones show a clumpy phenotype.</title><p>EvoPopulation7 had three morphologically distinct clones, named EvoClone7A, 7B, and 7C. For each of the remaining populations, we selected one clone from a group of eight morphologically indistinguishable clones. Representative DIC images were taken with a 40× objective in a glass-bottom 96-well plate. Scale bars are 50 µm.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.007">http://dx.doi.org/10.7554/eLife.00367.007</ext-link></p></caption><graphic xlink:href="elife00367fs001"/></fig><fig id="fig3s2" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.00367.008</object-id><label>Figure 3—figure supplement 2.</label><caption><title>Size distribution of evolved clones.</title><p>Sizes were measured on a Multisizer 3 Coulter counter. The thick line in the center of the box is the median size of the cells. The left and right side of the box correspond to the first and third quartiles. The lines extending from the boxes (the whiskers) correspond to the values within 1.5 times the inter-quartile range (IQR). Clumps or cells with sizes outside the whiskers are shown as dots. Cells were grown in the media shown for at least 12 hr, and sizes were measured during log-phase growth. Note that the x-axis is log-scale.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.008">http://dx.doi.org/10.7554/eLife.00367.008</ext-link></p></caption><graphic xlink:href="elife00367fs002"/></fig><fig id="fig3s3" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.00367.009</object-id><label>Figure 3—figure supplement 3.</label><caption><title>EvoClone9 morphology changes in different media.</title><p>Representative 40× DIC images of EvoClone9 cells that were grown to log phase in (<bold>A</bold>) 1 mM sucrose, (<bold>B</bold>) 80 mM glucose and (<bold>C</bold>) 1 mM glucose + 1 mM fructose. Scale bars are 50 µm.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.009">http://dx.doi.org/10.7554/eLife.00367.009</ext-link></p></caption><graphic xlink:href="elife00367fs003"/></fig></fig-group></p><p>We asked which strategies individual populations had adopted. Eight clones were selected from each population and their morphology and growth on sucrose were examined. For nine of the ten populations (all but EvoPopulation7), all the clones were morphologically identical; one clone was selected from each of these populations and used for further studies. EvoPopulation7 produced three different phenotypes: one had medium-sized clumps (EvoClone7A), one was not clumpy (EvoClone7B), and one had large clumps (EvoClone7C). Images of the clones are shown in <xref ref-type="fig" rid="fig3fs1">Figure 3—figure supplement 1</xref>, and size distributions of each clone are shown in <xref ref-type="fig" rid="fig3s2">Figure 3—figure supplement 2</xref>.</p><p>We tested the ability of the evolved clones to compete with their ancestor. Each of the evolved clones strongly outcompeted their ancestor, eliminating it from cultures within two passages on 1 mM sucrose (<xref ref-type="table" rid="tbl1">Table 1</xref>). We conclude that efficient use of public goods, liberated through extracellular hydrolysis, selects for the evolution of undifferentiated multicellularity.</p></sec><sec id="s2-3"><title>Clumps form by failure in cell separation</title><p>We asked how the evolved clones formed clumps. Wild yeast isolates can form multicellular clumps in two ways: flocculation (<xref ref-type="bibr" rid="bib29">Guo et al., 2000</xref>), in which separate cells stick to each other via cell wall-bound adhesins, or failure to separate daughters from their mothers because the cell wall that joins them is not digested after cytokinesis (<xref ref-type="bibr" rid="bib79">Yvert et al., 2003</xref>). To find which method the evolved clones used, we co-cultured two different colored versions of each evolved clone, using constitutive expression of different fluorescent proteins to mark the two versions. If the clumps formed by flocculation, many clumps would contain cells of both colors; if cells fail to separate, each clump would contain only one color since it would arise by the repeated division of a single cell.</p><p>All the multicellular clones are the result of incomplete separation; each clump contains cells of only one or the other color. Representative images from two evolved clones are shown in <xref ref-type="fig" rid="fig4">Figure 4</xref> (images from the remaining clones are shown in <xref ref-type="fig" rid="fig4s1">Figure 4—figure supplement 1</xref>). We checked our method by using two control strains: a flocculating strain that expresses a high level of Flo1, a known adhesin, (<xref ref-type="bibr" rid="bib68">Smukalla et al., 2008</xref>) and a lab strain that contains a wild isolate allele of <italic>AMN1</italic> (<italic>AMN1-RM11</italic>), which prevents cell separation after many of the cell divisions (<xref ref-type="bibr" rid="bib79">Yvert et al., 2003</xref>). <xref ref-type="fig" rid="fig4">Figure 4</xref> shows that the flocculating clumps from the strain expressing Flo1 contain both colors, whereas clumps from the <italic>AMN1-RM11</italic> strain contain only one color.<fig-group><fig id="fig4" position="float"><object-id pub-id-type="doi">10.7554/eLife.00367.010</object-id><label>Figure 4.</label><caption><title>Clumpiness is due to failure to separate and not flocculation.</title><p>Each image shows two genetically identical strains that are labeled with different fluorescent proteins, shown as magenta and green in the image. The strains were grown together from low density in 1 mM sucrose. (<bold>A</bold>) Lab strain with constitutively expressed <italic>FLO1</italic>. Flocculation is evident from the mix of colors in a single clump. (<bold>B</bold>) Lab strain with the RM11 allele of <italic>AMN1</italic>. (<bold>C</bold>) EvoClone2. (<bold>D</bold>) EvoClone9. The clumps in the <italic>AMN1-RM11</italic> strain and the evolved clones are uniform in color, showing that clumpiness is due to failure to separate after cell division. Representative fluorescent images were taken with a 20× objective in a glass-bottomed, 96-well plate. All scale bars are 50 µm.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.010">http://dx.doi.org/10.7554/eLife.00367.010</ext-link></p></caption><graphic xlink:href="elife00367f004"/></fig><fig id="fig4s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.00367.011</object-id><label>Figure 4—figure supplement 1.</label><caption><title>Clumpiness is due to failure to separate and not flocculation.</title><p>See <xref ref-type="fig" rid="fig4">Figure 4</xref> for a description of the experiment. Images of EvoClone2 and EvoClone9 in <xref ref-type="fig" rid="fig4">Figure 4</xref> are cropped versions of the images shown here for EvoClone2 and EvoClone9 (cropped area outlined in yellow). Scale bars are 50 µm.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.011">http://dx.doi.org/10.7554/eLife.00367.011</ext-link></p></caption><graphic xlink:href="elife00367fs004"/></fig></fig-group></p></sec><sec id="s2-4"><title>RNA sequencing shows that most clones elevate invertase expression and hexose transporter expression</title><p>We used RNA sequencing to examine the pattern of gene expression in the evolved clones. We isolated and sequenced RNA from the 12 evolved clones and two independent ancestor derivatives, all grown in 1 mM sucrose to log phase. In 10 of the 12 clones invertase (<italic>SUC2</italic>) expression was elevated between 3- and 21-fold above their ancestor (<xref ref-type="table" rid="tbl2">Table 2</xref>); EvoClone7C and EvoClone8 were the only clones with no significant change in <italic>SUC2</italic> expression.<table-wrap id="tbl2" position="float"><object-id pub-id-type="doi">10.7554/eLife.00367.012</object-id><label>Table 2.</label><caption><p>Significant changes in invertase and hexose transporter expression</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.012">http://dx.doi.org/10.7554/eLife.00367.012</ext-link></p></caption><table frame="hsides" rules="groups"><thead><tr><td>Strain name</td><td>Significant increases in invertase (<italic>SUC2</italic>) expression</td><td>Significant increases in hexose transporter expression</td></tr></thead><tbody><tr><td>EvoClone1</td><td>5X</td><td><italic>HXT1</italic> 42X</td></tr><tr><td/><td/><td><italic>HXT2</italic> 5X</td></tr><tr><td/><td/><td><italic>HXT3</italic> 108X</td></tr><tr><td/><td/><td><italic>HXT4</italic> 91X</td></tr><tr><td>EvoClone2</td><td>7X</td><td><italic>HXT4</italic> 9X</td></tr><tr><td>EvoClone3</td><td>7X</td><td><italic>HXT4</italic> 4X</td></tr><tr><td>EvoClone4</td><td>9X</td><td><italic>HXT1</italic> 4X</td></tr><tr><td/><td/><td><italic>HXT3</italic> 32X</td></tr><tr><td/><td/><td><italic>HXT4</italic> 14X</td></tr><tr><td/><td/><td><italic>HXT9</italic> 5X</td></tr><tr><td/><td/><td><italic>HXT11</italic> 4X</td></tr><tr><td>EvoClone5</td><td>7X</td><td><italic>HXT1</italic> 4X</td></tr><tr><td/><td/><td><italic>HXT2</italic> 8X</td></tr><tr><td/><td/><td><italic>HXT3</italic> 6X</td></tr><tr><td/><td/><td><italic>HXT4</italic> 21X</td></tr><tr><td/><td/><td><italic>HXT6</italic> 7X</td></tr><tr><td/><td/><td><italic>HXT7</italic> 6X</td></tr><tr><td>EvoClone6</td><td>9X</td><td><italic>HXT2</italic> 9X</td></tr><tr><td/><td/><td><italic>HXT4</italic> 12X</td></tr><tr><td>EvoClone7A</td><td>3X</td><td><italic>HXT2</italic> 7X</td></tr><tr><td/><td/><td><italic>HXT3</italic> 2X</td></tr><tr><td/><td/><td><italic>HXT4</italic> 92X</td></tr><tr><td>EvoClone7B</td><td>3X</td><td><italic>HXT2</italic> 10X</td></tr><tr><td/><td/><td><italic>HXT3</italic> 22X</td></tr><tr><td/><td/><td><italic>HXT4</italic> 103X</td></tr><tr><td>EvoClone7C</td><td>Not significant</td><td><italic>HXT1</italic> 3X</td></tr><tr><td/><td/><td><italic>HXT2</italic> 4X</td></tr><tr><td/><td/><td><italic>HXT3</italic> 10X</td></tr><tr><td/><td/><td><italic>HXT4</italic> 52X</td></tr><tr><td/><td/><td><italic>HXT6</italic> 4X</td></tr><tr><td>EvoClone8</td><td>Not significant</td><td><italic>HXT2</italic> 6X</td></tr><tr><td/><td/><td><italic>HXT3</italic> 5X</td></tr><tr><td/><td/><td><italic>HXT4</italic> 13X</td></tr><tr><td>EvoClone9</td><td>4X</td><td><italic>HXT2</italic> 5X</td></tr><tr><td>EvoClone10</td><td>21X</td><td>None</td></tr></tbody></table><table-wrap-foot><fn><p>See <xref ref-type="supplementary-material" rid="SD1-data">supplementary file 1C</xref> for list of all genes that were significantly changed and their change in each evolved clone.</p></fn></table-wrap-foot></table-wrap></p><p>RNA sequencing revealed another change that likely led to better growth in 1 mM sucrose: hexose transporter expression was elevated in 11 of 12 evolved clones. Yeast encodes at least 16 different hexose transporters, most encoded by members of the <italic>HXT</italic> gene family. The expression of <italic>HXT4</italic>, which encodes a high-affinity glucose transporter, was elevated in 10 of these evolved clones (<xref ref-type="table" rid="tbl2">Table 2</xref>); EvoClone10 was the only strain without increased expression of any hexose transporter.</p><p>The effect of evolution on genes with known roles in cell separation is shown in <xref ref-type="table" rid="tbl3">Table 3</xref>. Four genes involved in cell separation, <italic>AMN1</italic>, <italic>CTS1</italic>, <italic>DSE2</italic>, and <italic>SCW11</italic>, had significantly reduced levels in all 11 clumpy evolved clones, and three other genes implicated in cell separation showed reduced expression in 10 (<italic>DSE1, SUN4</italic>), or 9 (<italic>DSE4</italic>) of the clumpy clones. These seven genes are not decreased in expression in the non-clumpy evolved clone (EvoClone7B) or the control strain. This supports the argument that the clumps form by failure to separate. <xref ref-type="supplementary-material" rid="SD1-data">Supplementary file 1A</xref> lists genes whose expression was elevated or reduced in nine or more evolved clones (and unchanged in the control strain). We suspect many of these genes may show differences in expression because the evolved clones grow much more rapidly than their ancestors in low sucrose medium. <xref ref-type="supplementary-material" rid="SD1-data">Supplementary file 1B</xref> lists the genes whose expression was increased or decreased at least tenfold in each of the ten evolved clones we examined. All of the approximately 1500 genes with significant expression level changes are listed in <xref ref-type="supplementary-material" rid="SD1-data">supplementary file 1C</xref> along with their level of change in each evolved clone.<table-wrap id="tbl3" position="float"><object-id pub-id-type="doi">10.7554/eLife.00367.013</object-id><label>Table 3.</label><caption><p>Cell separation genes whose expression fell significantly</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.013">http://dx.doi.org/10.7554/eLife.00367.013</ext-link></p></caption><table frame="hsides" rules="groups"><thead><tr><td>Gene</td><td>Reduced in multicellular clones</td><td>Reduced in single cell clone</td><td>Function</td></tr></thead><tbody><tr><td><italic>AMN1</italic></td><td>11/11</td><td>0/1</td><td>Cell separation protein</td></tr><tr><td><italic>CTS1</italic></td><td>11/11</td><td>0/1</td><td>Cell separation, chitinase</td></tr><tr><td><italic>DSE2</italic></td><td>11/11</td><td>0/1</td><td>Cell separation, possible glucanase</td></tr><tr><td><italic>SCW11</italic></td><td>11/11</td><td>0/1</td><td>Cell separation, possible glucanase</td></tr><tr><td><italic>DSE1</italic></td><td>10/11</td><td>0/1</td><td>Cell separation, protein of unknown function</td></tr><tr><td><italic>SUN4</italic></td><td>10/11</td><td>0/1</td><td>Cell separation, possible glucanase</td></tr><tr><td><italic>DSE4</italic></td><td>9/11</td><td>0/1</td><td>Cell separation, possible glucanase</td></tr></tbody></table></table-wrap></p></sec><sec id="s2-5"><title>All evolved clones depend on secreted invertase</title><p>Forming multicellular clumps and increasing hexose transporter expression both suggest that evolved clones are still hydrolyzing sucrose extracellularly. But since we had shown that engineering sucrose import allowed growth on low sucrose, we needed to rule out the possibility that some of the evolved clones were using this strategy. If an evolved clone depended on sucrose import, keeping it from secreting invertase would have little effect on its growth in sucrose. To test for this possibility, we removed the signal sequence (<xref ref-type="bibr" rid="bib33">Kaiser and Botstein, 1986</xref>; <xref ref-type="bibr" rid="bib59">Perlman et al., 1986</xref>) that directs Suc2's secretion from each of the evolved clones (now named EvoCloneX-<italic>suc2-1cyt</italic>) and competed these derivatives against the corresponding, unmodified, evolved clone. In each case, the evolved clone quickly outcompeted the version that did not secrete invertase (<xref ref-type="table" rid="tbl4">Table 4</xref>), demonstrating that growth of all the evolved clones depends on secreted invertase. To verify that this growth defect was not due to reduction in invertase expression, we measured, using reverse transcription followed by quantitative PCR (RT-qPCR), the <italic>SUC2</italic> expression of the two evolved clones (EvoClone2 and EvoClone9) and their <italic>suc2-1cyt</italic> counterparts grown in 1 mM glucose. In both cases, expression of <italic>SUC2</italic> was slightly greater in the <italic>suc2-1cyt</italic> strain but statistically insignificant over three independent trials.<table-wrap id="tbl4" position="float"><object-id pub-id-type="doi">10.7554/eLife.00367.014</object-id><label>Table 4.</label><caption><p>Fitness of evolved clones after removal of <italic>SUC2</italic> signal sequence</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.014">http://dx.doi.org/10.7554/eLife.00367.014</ext-link></p></caption><table frame="hsides" rules="groups"><thead><tr><td>Strain 1</td><td>Strain 2</td><td>1 mM sucrose</td><td>1 mM glucose + 1 mM fructose</td></tr></thead><tbody><tr><td>EvoClone1</td><td>EvoClone1-<italic>suc2-1cyt</italic></td><td>++++</td><td>0</td></tr><tr><td>EvoClone2</td><td>EvoClone2-<italic>suc2-1cyt</italic></td><td>+++</td><td>0</td></tr><tr><td>EvoClone3</td><td>EvoClone3-<italic>suc2-1cyt</italic></td><td>++++</td><td>0</td></tr><tr><td>EvoClone4</td><td>EvoClone4-<italic>suc2-1cyt</italic></td><td>+++</td><td>0</td></tr><tr><td>EvoClone5</td><td>EvoClone5-<italic>suc2-1cyt</italic></td><td>++++</td><td>0</td></tr><tr><td>EvoClone6</td><td>EvoClone6-<italic>suc2-1cyt</italic></td><td>++++</td><td>+</td></tr><tr><td>EvoClone7A</td><td>EvoClone7A-<italic>suc2-1cyt</italic></td><td>+++</td><td>0</td></tr><tr><td>EvoClone7B</td><td>EvoClone7B-<italic>suc2-1cyt</italic></td><td>+++</td><td>0</td></tr><tr><td>EvoClone7C</td><td>EvoClone7C-<italic>suc2-1cyt</italic></td><td>++++</td><td>+</td></tr><tr><td>EvoClone8</td><td>EvoClone8-<italic>suc2-1cyt</italic></td><td>++++</td><td>+</td></tr><tr><td>EvoClone9</td><td>EvoClone9-<italic>suc2-1cyt</italic></td><td>++++</td><td>0</td></tr><tr><td>EvoClone10</td><td>EvoClone10-<italic>suc2-1cyt</italic></td><td>++++</td><td>+</td></tr><tr><td>wt</td><td>wt-<italic>suc2-1cyt</italic></td><td>+++</td><td>0</td></tr><tr><td>wt-<italic>suc2-1cyt</italic>-importer</td><td>wt-<italic>suc2-1cyt</italic></td><td>++++</td><td>0</td></tr></tbody></table><table-wrap-foot><fn><p>See <xref ref-type="table" rid="tbl1">Table 1</xref> for definition of fitness measurements. We used two control competitions: in the first, a standard lab strain outcompeted a lab strain with a missing <italic>SUC2</italic> secretion signal sequence (<italic>suc2-1cyt)</italic>; in the second, a <italic>suc2-1cyt</italic> strain with <italic>MAL11</italic> expressed from the <italic>ACT1</italic> promoter outcompeted the <italic>suc2-1cyt</italic> strain that did not express a sucrose importer.</p></fn></table-wrap-foot></table-wrap></p></sec><sec id="s2-6"><title>Evolved clones have poor fitness on other carbon sources</title><p>Did evolving to grow faster on low sucrose affect the ability of the evolved clones to grow on other carbon sources? Other experiments have resulted in antagonistic pleiotropy, where improved fitness in one environment often corresponds to reduced fitness in another (<xref ref-type="bibr" rid="bib41">Lenski, 1988</xref>; <xref ref-type="bibr" rid="bib76">Wenger et al., 2011</xref>). We therefore determined the fitness of the evolved clones in three environments: low (1 mM) sucrose, low monosaccharide (1 mM glucose plus 1 mM fructose, the hydrolysis products of 1 mM sucrose), and high (80 mM) glucose.</p><p>Each evolved clone quickly outcompeted the ancestor in low sucrose and lost to its ancestor in high glucose (<xref ref-type="table" rid="tbl1">Table 1</xref>). Four of the evolved clones also lost quickly on low monosaccharide; the remaining clones had approximately the same fitness as the ancestor, suggesting that most of the mutations acquired in low sucrose were selectively neutral in the equivalent monosaccharide concentration. The size distributions of each evolved clone was similar in all three conditions (<xref ref-type="fig" rid="fig3s2">Figure 3—figure supplement 2</xref>) except for EvoClone9, which formed much larger clumps in sucrose than it did in the other two media (<xref ref-type="fig" rid="fig3s2">Figure 3—figure supplements 2</xref> and <xref ref-type="fig" rid="fig3s3">3</xref>).</p></sec><sec id="s2-7"><title>Putative causal mutations found through bulk segregant analysis</title><p>We began the genetic characterization of the evolved clones by looking for the causal mutations. Because we used a strain that produced roughly one mutation per cell cycle, we predicted that neutral or slightly deleterious mutations would substantially outnumber the mutations that made cells grow faster in low sucrose concentrations. The neutral and nearly neutral mutations accumulate because they occur in lineages that were lucky enough to have strongly beneficial mutations; in the absence of sex, these mutations thus hitchhike during selection.</p><p>To identify the causal mutations, we used bulk segregant analysis, which uses sexual reproduction to separate causal from hitchhiking mutations. The evolved strain was crossed to its ancestor, put through meiosis, and a large population of haploid spores was selected for the evolved phenotype. Mutations that confer a strong advantage on low sucrose will be present in almost all the selected spores, whereas those that do not will be present in roughly half the spores. We monitored the allele frequency in the selected pool by preparing DNA and sequencing it to roughly 100-fold coverage (<xref ref-type="fig" rid="fig5">Figure 5</xref>). We found a total of 80 putative causal mutations in the twelve evolved clones out of 1521 mutations, confirming that most mutations are non-causal (<xref ref-type="table" rid="tbl5">Table 5</xref>). To track the spread of the putative causal mutations through the evolved populations, we used Sanger sequencing to measure the allele frequency (<xref ref-type="bibr" rid="bib27">Gresham et al., 2008</xref>) over time. <xref ref-type="fig" rid="fig6">Figure 6</xref> shows the spread of putative causal mutations through two of the populations, EvoClone2 and EvoClone9; <xref ref-type="fig" rid="fig6s1">Figure 6—figure supplements 1</xref> and <xref ref-type="fig" rid="fig6s2">2</xref> show the remaining populations.<fig-group><fig id="fig5" position="float"><object-id pub-id-type="doi">10.7554/eLife.00367.015</object-id><label>Figure 5.</label><caption><title>Schematic of bulk segregant analysis and evolved clone reconstruction.</title><p>(<bold>A</bold>) A clone is selected from the population and then backcrossed to a derivative of its ancestor. The resulting diploid is sporulated, allowing the mutant alleles to randomly segregate among the haploid progeny. When the haploid progeny are selected for growth in low sucrose, only those cells with causal alleles (red triangles) remain; non-causal alleles (blue diamond, square, and circle) segregate randomly and are present in about half of the spores. (<bold>B</bold>) The ancestor, evolved clone, and pool of selected progeny are sequenced. Comparing the genome sequences of the ancestor and evolved clone reveals mutations. The allele frequency in the selected spores can then be estimated from the frequency of the reads in the pool of selected progeny. We classified any mutant allele present in >90% of the reads as a putative causal mutation (<xref ref-type="table" rid="tbl5">Table 5</xref>). (<bold>C</bold>) The wild-type alleles in the ancestor were replaced with the putative causal mutations to recreate the evolved clone (<xref ref-type="fig" rid="fig7">Figure 7</xref>). (<bold>D</bold>) Growth of the recreated strain was tested in low sucrose (<xref ref-type="table" rid="tbl1">Table 1</xref>).</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.015">http://dx.doi.org/10.7554/eLife.00367.015</ext-link></p></caption><graphic xlink:href="elife00367f005"/></fig><fig id="fig5s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.00367.016</object-id><label>Figure 5—figure supplement 1.</label><caption><title>Protocol for replacing alleles in yeast.</title><p>(1) design the region of the gene to be included in the plasmid as follows: (<bold>A</bold>). Determine the ‘usable' gene, or the promoter plus the open reading frame (ORF) plus the terminator. (<bold>B</bold>) Figure out which end of the usable gene is closer to the mutation. This is the side that will not be truncated. (<bold>C</bold>) Find a cut site >300 bp away from the mutation. If the new mutation is needed after transformation and before loopout, place the cut site toward far end of gene. This is the end of the gene that will be truncated. If old function is needed, place the cut site toward near end of gene. The region between the mutation and the cut site is now the ‘plasmid' region. The distance between the cut site and the mutation is necessary because of recision. (<bold>D</bold>) Extend the plasmid region at least 300 bp onto the opposite side of the cut site. This distance ensures efficient transformation. (<bold>E</bold>) Extend the plasmid region to the nearest end of the usable gene. This will ensure the gene is intact after transformation. (<bold>F</bold>) Find the distance from the mutation to each end of the plasmid region. Extend the plasmid region so that the length that does not include the cut site is larger than the length that does include the cut site. This ensures that loopout is more likely to result in the new mutation. (2) PCR amplify the plasmid region of the gene, treat it with polynucleotide kinase, and blunt-end ligate it into a digested <italic>URA3</italic> plasmid backbone that has been treated with Antarctic phosphatase. Only pertinent regions of the plasmid are shown in the diagram. An additional yeast drug marker is useful in the plasmid backbone to verify that the insert has looped out. (3) Transform the yeast with the cut plasmid and select on—URA. This will integrate the linearized plasmid into the chromosome. PCR amplify and sequence the new allele to verify insertion; one primer should be outside the included region. (4) Grow the transformed yeast strain overnight in YPD to allow the insert to loop out through homologous recombination and plate on 5FOA. Cells will only grow on 5FOA if the <italic>URA3</italic> marker has looped out. PCR amplify and sequence the new allele with both primers inside the included region to verify that the new allele is the only copy of the gene remaining. Replica plate to YP 2% acetate to eliminate petite mutants.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.016">http://dx.doi.org/10.7554/eLife.00367.016</ext-link></p></caption><graphic xlink:href="elife00367fs005"/></fig></fig-group><table-wrap id="tbl5" position="float"><object-id pub-id-type="doi">10.7554/eLife.00367.017</object-id><label>Table 5.</label><caption><p>Putative causal mutations in each evolved clone</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.017">http://dx.doi.org/10.7554/eLife.00367.017</ext-link></p></caption><table frame="hsides" rules="groups"><thead><tr><td>Strain name</td><td>Nominal generations</td><td>Number of mutations</td><td>Non-synonymous and promoter mutations segregating at evolved allele frequency >90%</td><td>Nucleotide change</td><td>Amino acid change</td><td>Mutant allele reads/total reads</td></tr></thead><tbody><tr><td>EvoClone1</td><td>307</td><td>71</td><td><italic>ACE2</italic></td><td>703 G→T</td><td>E235*</td><td>77/79</td></tr><tr><td/><td/><td/><td><italic>IRA1</italic></td><td>8987C→G</td><td>T2996S</td><td>99/102</td></tr><tr><td/><td/><td/><td><italic>PHO87</italic></td><td>196 G→C</td><td>A66P</td><td>55/55</td></tr><tr><td/><td/><td/><td><italic>RGT1</italic></td><td>3157C→T</td><td>Q1053*</td><td>72/74</td></tr><tr><td/><td/><td/><td><italic>SAN1</italic></td><td>1707C→A</td><td>N569K</td><td>76/77</td></tr><tr><td/><td/><td/><td><italic>SIN4</italic></td><td>383 G→A</td><td>G128D</td><td>106/106</td></tr><tr><td/><td/><td/><td><italic>UBR1</italic></td><td>1916T→A</td><td>L639*</td><td>42/44</td></tr><tr><td>EvoClone2</td><td>273</td><td>115</td><td><italic>ACE2</italic></td><td>968T→A</td><td>L323*</td><td>110/110</td></tr><tr><td/><td/><td/><td><italic>CSE2</italic></td><td>100_100delT</td><td>S35Rfs54</td><td>102/107</td></tr><tr><td/><td/><td/><td><italic>IRA1</italic></td><td>7657A→T</td><td>S2553C</td><td>72/75</td></tr><tr><td/><td/><td/><td><italic>MTH1</italic></td><td>459_459delC</td><td>H154Tfs156</td><td>93/100</td></tr><tr><td/><td/><td/><td><italic>UBR1</italic></td><td>524 G→A</td><td>C175Y</td><td>104/112</td></tr><tr><td>EvoClone3</td><td>301</td><td>252</td><td><italic>IRA2</italic></td><td>8081C→A</td><td>S2694Y</td><td>75/81</td></tr><tr><td/><td/><td/><td><italic>IRC8</italic></td><td>365T→C</td><td>L122P</td><td>84/88</td></tr><tr><td/><td/><td/><td><italic>NAT1</italic></td><td>1782 G→A</td><td>W594*</td><td>84/87</td></tr><tr><td/><td/><td/><td><italic>SYP1</italic></td><td>1376C→T</td><td>T459I</td><td>90/96</td></tr><tr><td>EvoClone4</td><td>229</td><td>95</td><td><italic>ACE2</italic></td><td>670 G→T</td><td>E224*</td><td>93/93</td></tr><tr><td/><td/><td/><td><italic>RGT1</italic></td><td>2494_2495insT</td><td>L832Ffs834</td><td>135/139</td></tr><tr><td/><td/><td/><td><italic>SIN4</italic></td><td>382 G→A</td><td>G128S</td><td>104/105</td></tr><tr><td/><td/><td/><td><italic>UBR1</italic></td><td>1916T→A</td><td>L639*</td><td>56/61</td></tr><tr><td>EvoClone5</td><td>237</td><td>120</td><td><italic>ACE2</italic></td><td>565C→T</td><td>Q189*</td><td>149/152</td></tr><tr><td/><td/><td/><td><italic>ARO2</italic></td><td>371C→G</td><td>A124G</td><td>93/101</td></tr><tr><td/><td/><td/><td><italic>MCK1</italic></td><td>38_38delG</td><td>G14Dfs22</td><td>152/158</td></tr><tr><td/><td/><td/><td><italic>SNF2</italic></td><td>71 G→T</td><td>R24I</td><td>94/100</td></tr><tr><td/><td/><td/><td><italic>SNF3</italic></td><td>1235T→A</td><td>V412E</td><td>148/148</td></tr><tr><td>EvoClone6</td><td>232</td><td>110</td><td><italic>ACE2</italic></td><td>507_507delT</td><td>N169Kfs177</td><td>142/146</td></tr><tr><td/><td/><td/><td><italic>GCN2</italic></td><td>892A→G</td><td>N298D</td><td>115/117</td></tr><tr><td/><td/><td/><td><italic>GPB2</italic></td><td>235 G→T</td><td>E79*</td><td>144/145</td></tr><tr><td/><td/><td/><td><italic>MTH1</italic></td><td>152_152delG</td><td>S51Ifs56</td><td>110/112</td></tr><tr><td/><td/><td/><td><italic>NRG1</italic></td><td>371C→A</td><td>S124*</td><td>92/95</td></tr><tr><td/><td/><td/><td><italic>RAD6</italic></td><td>191C→A</td><td>P64H</td><td>127/135</td></tr><tr><td>EvoClone7A</td><td>242</td><td>105</td><td><italic>ACE2</italic></td><td>1901C→A</td><td>S634*</td><td>84/87</td></tr><tr><td/><td/><td/><td><italic>RAD61</italic></td><td>225T→A</td><td>N75K</td><td>97/97</td></tr><tr><td>EvoClone7B</td><td>242</td><td>94</td><td><italic>GCN3</italic></td><td>176C→A</td><td>S59Y</td><td>79/87</td></tr><tr><td/><td/><td/><td><italic>IRA2</italic></td><td>7049_7049delC</td><td>A2350Gfs2354</td><td>129/132</td></tr><tr><td/><td/><td/><td><italic>RAM1</italic></td><td>566T→A</td><td>L189Q</td><td>133/137</td></tr><tr><td/><td/><td/><td><italic>SAN1</italic></td><td>1464C→G</td><td>N488K</td><td>97/100</td></tr><tr><td/><td/><td/><td><italic>SNF3</italic></td><td>692C→A</td><td>A231D</td><td>87/90</td></tr><tr><td>EvoClone7C</td><td>242</td><td>115</td><td><italic>GCR2</italic></td><td>533T→A</td><td>L178Q</td><td>148/154</td></tr><tr><td/><td/><td/><td><italic>IRA2</italic></td><td>7049_7049delC</td><td>A2350Gfs2354</td><td>166/169</td></tr><tr><td/><td/><td/><td><italic>PDR1</italic></td><td>2527 G→A</td><td>D843N</td><td>145/159</td></tr><tr><td/><td/><td/><td><italic>PUF4</italic></td><td>1960C→T</td><td>Q654*</td><td>184/202</td></tr><tr><td>EvoClone8</td><td>253</td><td>122</td><td><italic>ACE2</italic></td><td>−379 G→A</td><td>Promoter</td><td>103/103</td></tr><tr><td/><td/><td/><td><italic>AXL2</italic></td><td>700T→C</td><td>S234P</td><td>87/87</td></tr><tr><td/><td/><td/><td><italic>ERG1</italic></td><td>427 G→A</td><td>E143K</td><td>160/163</td></tr><tr><td/><td/><td/><td><italic>HXK1</italic></td><td>93A→T</td><td>E31D</td><td>200/217</td></tr><tr><td/><td/><td/><td><italic>IFM1</italic></td><td>1724T→A</td><td>I575N</td><td>130/134</td></tr><tr><td/><td/><td/><td><italic>MIT1</italic></td><td>188 G→A</td><td>W63*</td><td>131/142</td></tr><tr><td/><td/><td/><td><italic>SKS1</italic></td><td>1311C→G</td><td>Y437*</td><td>152/152</td></tr><tr><td/><td/><td/><td><italic>SNF3</italic></td><td>1237 G→A</td><td>E413K</td><td>101/101</td></tr><tr><td/><td/><td/><td><italic>UBC5</italic></td><td>443A→G</td><td>D118G</td><td>137/141</td></tr><tr><td/><td/><td/><td><italic>UBR1</italic></td><td>3859_3859delG</td><td>G1287Dfs1345</td><td>170/179</td></tr><tr><td>EvoClone9</td><td>265</td><td>126</td><td><italic>ARE1</italic></td><td>−10 G→T</td><td>Promoter</td><td>127/127</td></tr><tr><td/><td/><td/><td><italic>GCN2</italic></td><td>4582A→C</td><td>I1528L</td><td>125/137</td></tr><tr><td/><td/><td/><td><italic>GIN4</italic></td><td>57 G→A</td><td>W19*</td><td>79/81</td></tr><tr><td/><td/><td/><td><italic>IRC8</italic></td><td>170 G→T</td><td>G57V</td><td>118/119</td></tr><tr><td/><td/><td/><td><italic>MCD1</italic></td><td>524C→T</td><td>S175L</td><td>115/115</td></tr><tr><td/><td/><td/><td><italic>MCK1</italic></td><td>675_675delG</td><td>G227Vfs249</td><td>114/116</td></tr><tr><td/><td/><td/><td><italic>MED1</italic></td><td>1009C→G</td><td>L337V</td><td>104/113</td></tr><tr><td/><td/><td/><td/><td>1465 G→T</td><td>E489*</td><td>118/125</td></tr><tr><td/><td/><td/><td><italic>UBR1</italic></td><td>3148_3148delC</td><td>L1050Yfs1063</td><td>141/145</td></tr><tr><td>EvoClone10</td><td>242</td><td>196</td><td><italic>ACE2</italic></td><td>1874A→T</td><td>Q625L</td><td>121/122</td></tr><tr><td/><td/><td/><td><italic>AXL2</italic></td><td>432_432delC</td><td>Y145Mfs154</td><td>161/164</td></tr><tr><td/><td/><td/><td><italic>BPH1</italic></td><td>2369C→A</td><td>S790Y</td><td>144/144</td></tr><tr><td/><td/><td/><td><italic>DNF2</italic></td><td>2351T→C</td><td>F784S</td><td>90/97</td></tr><tr><td/><td/><td/><td><italic>ECM5</italic></td><td>3466 G→A</td><td>D1156N</td><td>127/128</td></tr><tr><td/><td/><td/><td><italic>ENP2</italic></td><td>1129T→A</td><td>F377I</td><td>127/132</td></tr><tr><td/><td/><td/><td><italic>GAC1</italic></td><td>−7T→A</td><td>Promoter</td><td>106/117</td></tr><tr><td/><td/><td/><td><italic>HTZ1</italic></td><td>−369T→C</td><td>Promoter</td><td>69/73</td></tr><tr><td/><td/><td/><td><italic>KEM1</italic></td><td>2268 G→A</td><td>M756I</td><td>148/148</td></tr><tr><td/><td/><td/><td><italic>MCD1</italic></td><td>−28 G→T</td><td>Promoter</td><td>153/153</td></tr><tr><td/><td/><td/><td><italic>MPT5</italic></td><td>2409T→A</td><td>L590*</td><td>146/146</td></tr><tr><td/><td/><td/><td><italic>MRPS17</italic></td><td>325 G→A</td><td>D109N</td><td>141/147</td></tr><tr><td/><td/><td/><td><italic>NUT1</italic></td><td>2582C→A</td><td>S861*</td><td>166/169</td></tr><tr><td/><td/><td/><td><italic>PRC1</italic></td><td>−283 G→A</td><td>Promoter</td><td>138/148</td></tr><tr><td/><td/><td/><td><italic>RGT1</italic></td><td>2060G→T</td><td>G687V</td><td>91/91</td></tr><tr><td/><td/><td/><td><italic>SAC6</italic></td><td>1736A→T</td><td>K542M</td><td>125/137</td></tr><tr><td/><td/><td/><td><italic>TOP3</italic></td><td>1679T→C</td><td>V560A</td><td>100/103</td></tr><tr><td/><td/><td/><td><italic>UBR1</italic></td><td>56T→A</td><td>L19Q</td><td>130/140</td></tr><tr><td/><td/><td/><td><italic>WHI2</italic></td><td>187 G→T</td><td>E63*</td><td>138/139</td></tr><tr><td/><td/><td/><td><italic>WTM2</italic></td><td>−297T→A</td><td>Promoter</td><td>125/129</td></tr></tbody></table><table-wrap-foot><fn><p>Nomenclature based on (<xref ref-type="bibr" rid="bib18">den Dunnen and Antonarakis 2000</xref>): Mx→y: nucleotide change from x to y at base M, starting at base 1 (negative indicates promoter region). M_Ndelx: Deletion (Insertion:ins) of nucleotide x from base M to N. XNY: amino acid change from X to Y at codon N. * indicates stop codon. XNYfsN: as above, plus a frame shift mutation that results in stop codon at N.</p></fn><fn><p>The following mutations likely hitchhiked and were not included in this table: EvoClone3: <italic>atg4</italic> with <italic>sin4</italic>; EvoClone7A: <italic>slx4</italic> promoter with <italic>ace2</italic>; EvoClone7B: <italic>thi3</italic> with <italic>ram1</italic>, <italic>crt10</italic> with <italic>ira2</italic>, EvoClone8: <italic>crh1</italic> promoter with <italic>ubr1</italic>; EvoClone9: <italic>brr1</italic> with <italic>MED1</italic>, <italic>nsp1</italic> promoter with <italic>irc8</italic>; EvoClone10: <italic>pri2</italic> promoter with <italic>rgt1</italic>, <italic>ino80</italic> with <italic>nut1</italic>. This claim is based on the genetic linkage between the two alleles and the lower allele frequency of the mutation we argue is hitch-hiking relative to the putative causal mutation. The following mutations are not shown in the time courses in <xref ref-type="fig" rid="fig6">Figure 6</xref> and <xref ref-type="fig" rid="fig6s1">Figure 6—supplements 1</xref> and <xref ref-type="fig" rid="fig6s2">2</xref> because they were present at frequencies of less than 5% of the final population: EvoClone7C: <italic>puf4-Q654*</italic>; EvoClone8: <italic>mit1-W63*</italic>; EvoClone9 <italic>mcd1-S175L</italic>; EvoClone10: <italic>top3-V560A</italic> The following mutations were in the original, time zero strain and are not included in this table even through they segregated at > 90%: EvoClone2: <italic>ira1-F664I</italic>; EvoClone8: <italic>phm8-I97N</italic>, <italic>rpl37a</italic>(−52T→G), <italic>vta1-A247V</italic>, <italic>yor1-E393D</italic>; EvoClone9: <italic>vta1-A247V</italic>; EvoClone10: <italic>aim32-E241G</italic>, <italic>irs4-N257S</italic>, <italic>nnt1</italic>(−427T→C), <italic>prp9-N155S</italic>, <italic>yrb1-N120I</italic>.</p></fn></table-wrap-foot></table-wrap><fig-group><fig id="fig6" position="float"><object-id pub-id-type="doi">10.7554/eLife.00367.018</object-id><label>Figure 6.</label><caption><title>Putative causal mutation frequency at time points during the evolution.</title><p>The alleles at the indicated time points were sequenced using Sanger sequencing, and frequencies were estimated from peaks in the trace plots. See <xref ref-type="fig" rid="fig6s1">Figure 6—figure supplements 1</xref> and <xref ref-type="fig" rid="fig6s2">2</xref> for the other evolved populations. See <xref ref-type="table" rid="tbl5">Table 5</xref> for amino acid and nucleotide changes.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.018">http://dx.doi.org/10.7554/eLife.00367.018</ext-link></p></caption><graphic xlink:href="elife00367f006"/></fig><fig id="fig6s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.00367.019</object-id><label>Figure 6—figure supplement 1.</label><caption><title>Putative causal mutation frequency at time points during the evolution for EvoClone 1, 3, 4, 5, 6, 7A, 7B, and 7C.</title><p>See <xref ref-type="fig" rid="fig6">Figure 6</xref>.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.019">http://dx.doi.org/10.7554/eLife.00367.019</ext-link></p></caption><graphic xlink:href="elife00367fs006"/></fig><fig id="fig6s2" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.00367.020</object-id><label>Figure 6—figure supplement 2.</label><caption><title>Putative causal mutation frequency at time points during the evolution for EvoClone 8 and 10.</title><p>See <xref ref-type="fig" rid="fig6">Figure 6</xref>.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.020">http://dx.doi.org/10.7554/eLife.00367.020</ext-link></p></caption><graphic xlink:href="elife00367fs007"/></fig></fig-group></p><p>Over all 12 clones, the 80 putative causal mutations lie in or near 53 genes. Two genes, <italic>ACE2</italic> and <italic>UBR1,</italic> were mutated in at least half the clones, and groups of genes in three pathways, involved in glucose sensing (the <italic>RGT1</italic> group), growth regulation (the <italic>IRA1/IRA2</italic> group), and transcription (the Mediator group) were also mutated in many clones. We suspect that most of the putative causal mutations are loss of function mutations. Many alleles of the two most frequently mutated genes, <italic>ACE2</italic> and <italic>UBR1</italic> were nonsense mutations. In addition, amongst the 39 genes that were only mutated once, 12 (31%) of the mutations were nonsense mutations. This suggests that many of the mutations in the remaining 27 genes are loss of function mutations: for two genes that have been extensively studied, <italic>URA3</italic> and <italic>CAN1</italic>, 31% of strong loss of function mutations are nonsense mutations (<xref ref-type="bibr" rid="bib40">Lang and Murray, 2008</xref>). (Because loss-of-function mutations are typically recessive, we use lower case italic nomenclature to indicate the putative causal allele (e.g., <italic>irc8-G57V</italic>), and upper case italics to indicate the wild-type allele (e.g., <italic>IRC8</italic>). See the notes for <xref ref-type="table" rid="tbl5">Table 5</xref> for a description of the mutation nomenclature. Roman text (e.g., Irc8) indicates the protein). <xref ref-type="supplementary-material" rid="SD2-data">Supplementary file 2</xref> is a summary of the mutations in each gene pathway.</p><p>Ace2 activates the transcription of enzymes that degrade the septum that connects mother to daughter cell (<xref ref-type="bibr" rid="bib17">Colman-Lerner et al., 2001</xref>; <xref ref-type="bibr" rid="bib65">Sbia et al., 2008</xref>). <italic>ACE2</italic> was mutated in eight evolved clones; six of these mutations were nonsense mutations scattered through the open reading frame. All six of these clones had at least tenfold reductions in the expression of three genes known to aid in cell separation (<italic>CTS1</italic> [chitinase], <italic>DSE1</italic>, and <italic>DSE2</italic> [<xref ref-type="supplementary-material" rid="SD1-data">supplementary file 1B</xref>]) in agreement with the original finding that <italic>ace2Δ</italic> mutants are clumpy (<xref ref-type="bibr" rid="bib22">Dohrmann et al., 1992</xref>).</p><p>Ubr1 is an E3 ubiquitin ligase in the N-end rule pathway. Ubr1 targets proteins with certain N-terminal amino acids, which are exposed by proteolytic cleavage or removal of the N-terminal methionine, for degradation (<xref ref-type="bibr" rid="bib6">Bartel et al., 1990</xref>). <italic>UBR1</italic> was mutated in 6 of the 12 clones; four mutations were nonsense mutations, strongly suggesting that the remaining two are also a loss of function mutations. Rad6, the ubiquitin-conjugating enzyme that forms a heterodimer with Ubr1 (<xref ref-type="bibr" rid="bib21">Dohmen et al., 1991</xref>), had a missense mutation in a seventh evolved clone.</p><p>The Rgt1 pathway controls cells' response to external glucose concentrations. In high glucose, Rgt1, a transcription factor, represses the expression of high affinity hexose transporters such as <italic>HXT4</italic> (<xref ref-type="bibr" rid="bib57">Ozcan and Johnston, 1999</xref>). In low glucose, Snf3, a membrane-spanning glucose sensor, induces degradation of Mth1, a protein that is needed for Rgt1 to extert its repressive effects (<xref ref-type="bibr" rid="bib11">Broach, 2012</xref>). One of the genes in the Rgt1 pathway was mutated in 8 of the 12 clones: <italic>SNF3</italic> (3 missense), <italic>MTH1</italic> (2 nonsense), and <italic>RGT1</italic> (2 nonsense, 1 missense); no clones had two mutations in the pathway. Because Mth1 and Rgt1 work together to repress genes involved in growth on poor carbon sources, we expect that the mutations we isolated are loss of function mutations. <italic>MTH1</italic> loss of function mutations have been selected by others in glucose-limited chemostats (<xref ref-type="bibr" rid="bib34">Kao and Sherlock, 2008</xref>; <xref ref-type="bibr" rid="bib38">Kvitek and Sherlock, 2011</xref>). On the other hand, because Snf3's activity leads to the induction of glucose-repressed genes, we speculate that the <italic>SNF3</italic> mutations may result in gain of function. One of our mutations, <italic>snf3-A231D,</italic> lies near a known dominant mutation, <italic>SNF3</italic>-<italic>R229K</italic> (<xref ref-type="bibr" rid="bib56">Ozcan et al., 1996</xref>), that elevates hexose transporter expression.</p><p>The Ras-cAMP pathway regulates cell growth in response to nutrients (<xref ref-type="bibr" rid="bib11">Broach, 2012</xref>). Ras1 and Ras2 are small G proteins whose activity responds to nutrient sensing and who stimulate the activity of adenyl cyclase, which produces cAMP. Ira1 and Ira2 both encode GTPase activating proteins that inactivate Ras by stimulating its intrinsic GTPase activity (<xref ref-type="bibr" rid="bib71">Tanaka et al., 1989</xref>, <xref ref-type="bibr" rid="bib72">1990</xref>). <italic>IRA1</italic> or <italic>IRA2</italic> were mutated in four unrelated evolved clones (3 missense, 1 nonsense). Mutations that increase cAMP-dependent protein kinase activity are known to prevent starvation-induced cell cycle arrest (<xref ref-type="bibr" rid="bib50">Matsumoto et al., 1983</xref>). By increasing Ras activity, loss of function mutations in <italic>IRA1</italic> and <italic>IRA2</italic> will elevate cAMP, which may allow cells to grow at lower external sugar concentrations than their ancestors.</p><p>Mediator is a multiprotein, global regulator of eukaryotic transcription. We found mutations in four of the at least 20 proteins that make up Mediator (<xref ref-type="bibr" rid="bib55">Myers and Kornberg, 2000</xref>): <italic>SIN4</italic> (2 missense mutations), <italic>CSE2</italic> (1 nonsense mutation), <italic>MED1</italic> (1 nonsense mutation), and <italic>NUT1</italic> (1 nonsense mutation). Both mutations in <italic>SIN4</italic> occurred in residue 128 (G128D, G128S) and one, <italic>sin4-G128D</italic>, was also identified in a screen for suppressors of a deletion in <italic>SWI6.</italic> Swi6 associates with Swi4 to form SBF (SCB [Swi4-Swi6 cell cycle box] binding factor), a complex that regulates transcription early in the yeast cell cycle (<xref ref-type="bibr" rid="bib46">Li et al., 2005</xref>). The <italic>SIN4</italic> mutations thus may cause increased expression of a gene normally activated by SBF. The hexose transporter gene, <italic>HXT3</italic>, (<xref ref-type="bibr" rid="bib31">Iyer et al., 2001</xref>) is a known target of SBF and its was strongly elevated (108-fold and 32-fold) in the two clones (EvoClone1 and EvoClone4) carrying <italic>SIN4</italic> mutation (<xref ref-type="table" rid="tbl2 tbl5">Table 2 and 5</xref>).</p><p>We analyzed read depth across all evolved clones and saw three major duplication or deletion events: (1) chromosome 3 in EvoClone10 was duplicated in the region between 152 and 172 kb from the left end of the chromosome. This doubling of read depth was reduced to a 50% increase in the segregated pool, indicating that the duplication was not causal. (2) The number of P-Type ATPases at the <italic>ENO1/2/5</italic> locus was reduced from 3 to 2 in EvoClone6 and EvoClone10. This reduction was also present in the backcrossed strain, indicating that it may have been causal. (3) The number of hexose transporter gene repeats at the <italic>HXT6/7</italic> locus increased from two in the ancestral strain to three in EvoClone7C. This amplification was also present in the backcrossed strain, indicating that it also may have been causal. Amplification of the <italic>HXT6/7</italic> locus has been found in other experimental evolutions (<xref ref-type="bibr" rid="bib13">Brown et al., 1998</xref>; <xref ref-type="bibr" rid="bib27">Gresham et al., 2008</xref>; <xref ref-type="bibr" rid="bib34">Kao and Sherlock, 2008</xref>).</p></sec><sec id="s2-8"><title>Recreation of evolved clones verifies causal mutations</title><p>Are the putative causal mutations really responsible for the evolved phenotypes? We addressed this question by engineering the candidate mutations from two clones, EvoClone2 and EvoClone9, into their ancestor and asking if this manipulation reproduced the behavior of these clones. To recreate the sets of putative causal mutations, we replaced the five ancestral alleles in the ancestor with the five putative causal mutations in EvoClone2 to make Recreated2 and the eight putative causal mutations in EvoClone9 to make Recreated9. <xref ref-type="fig" rid="fig7">Figure 7</xref> and <xref ref-type="fig" rid="fig3s2">Figure 3—figure supplement 2</xref> shows that the morphology and the clump size distribution of the recreated strains are similar to that of the evolved clones.<fig id="fig7" position="float"><object-id pub-id-type="doi">10.7554/eLife.00367.021</object-id><label>Figure 7.</label><caption><title>Engineering in alleles can recreate the evolved and ancestral morphologies.</title><p>The ancestral strain was converted to the evolved morphology by converting ancestral alleles to those of the putative causal mutations and the evolved strains were converted to the ancestral morphology by converting the putative causal mutations to their ancestral alleles. The strains were grown separately in 1 mM sucrose and then mixed. The top row shows EvoClone2 strains and the bottom row shows EvoClone9 strains. The ancestor constitutively expresses mCherry and is shown in yellow; the evolved clone constitutively expresses YFP and is shown in green. The recreated strain (left) and the reverted strain (right) constitutively express CFP and are shown in magenta (recreated evolved) and cyan (reverted to ancestral). Representative confocal fluorescent (left) and brightfield (right) images were taken with a 60× objective in a glass-bottomed, 96-well plate. All scale bars are 50 µm.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.021">http://dx.doi.org/10.7554/eLife.00367.021</ext-link></p></caption><graphic xlink:href="elife00367f007"/></fig></p><p>To assess the fitness of the reconstructed strains, we competed them against their ancestor and the evolved clones. Both the evolved and recreated clones outcompeted their ancestor in 1–2 passages (<xref ref-type="table" rid="tbl1">Table 1</xref>) showing that they are much fitter on low sucrose. Both recreated strains were slightly less fit their evolved counterparts: there were fewer recreated than evolved cells after six passages (46 generations). This fitness defect has three possible causes: (1) the recreated clones are missing minor causal alleles, which we failed to find; (2) detrimental mutations were introduced during the multiple transformations required to make the recreated strains; or, (3) the evolved clone (which is a mutator, unlike the recreated strain) continued to evolve and adapt during the competition.</p><p>To confirm that the putative causal alleles accounted for the evolved phenotype, we reverted these alleles in the evolved clones to their ancestral state. We replaced each of the mutant alleles in EvoClone2 and EvoClone9 with the ancestral allele. Both the resulting strains, Reverted2 and Reverted9, grew poorly in 1 mM sucrose, mostly existed as single cells, and were quickly outcompeted by the evolved clones (<xref ref-type="table" rid="tbl1">Table 1</xref>, <xref ref-type="fig" rid="fig7">Figure 7</xref>, and <xref ref-type="fig" rid="fig3s2">Figure 3—figure supplement 2</xref>).</p><p>We used bulk segregant analysis to ask whether individual alleles were required for the evolved phenotype. We crossed the recreated strains to their ancestors, isolated haploid spores, selected this population for growth on 1 mM sucrose, and measured the frequency of the evolved alleles in the selected population. All five mutations in Recreated2 were strongly selected for, six of the eight mutations in Recreated9 were strongly selected for, and one mutation in Recreated9 was moderately selected for, confirming that all but one of the thirteen putative causal mutations were causal for improved growth on low sucrose (<xref ref-type="fig" rid="fig8">Figure 8</xref>). We suspect that the eighth mutation in Recreated9, <italic>gcn2-I1528L</italic>, may be a false positive. In the original spore selection, the mutant allele segregated at only 91%, the low end of our threshold, and the mutation changed a poorly conserved isoleucine to leucine, suggesting it is unlikely to have a major effect on the protein's activity. In sum, our recreation experiments show that the putative causal alleles, with a single exception, are both necessary and sufficient to produce the evolved phenotype.<fig-group><fig id="fig8" position="float"><object-id pub-id-type="doi">10.7554/eLife.00367.022</object-id><label>Figure 8.</label><caption><title>Bulk segregant analysis with the recreated strains verifies causal alleles and shows that alleles responsible for clumpiness are selected in low sucrose and not low monosaccharide.</title><p>The recreated strains were backcrossed, sporulated, and selected in three different media: 1 mM sucrose (low sucrose), 80 mM glucose (high glucose), and 1 mM glucose plus 1 mM fructose (low monosaccharide). The mutant allele fraction was estimated from Sanger sequencing across the allelic variants. The size of the data point (small, medium, or large) for each allele and media combination refers to one of three independently derived diploids. (<bold>A</bold>) Recreated2 allele segregation. (<bold>B</bold>) Recreated9 allele segregation. (<bold>C</bold>) Ancestor strain with <italic>ACE2</italic> deleted (<italic>ace2Δ</italic>) has a clumpy phenotype. The <italic>ACE2</italic> mutation in Recreated2, a likely loss of function mutation that caused the clumpiness in EvoClone2, was selected for in low sucrose and was not selected for in low monosaccharide. (<bold>D</bold>) Ancestor strain that has wild type alleles of <italic>IRC8</italic>, <italic>MCK1</italic>, and <italic>GIN4</italic> replaced with the EvoClone9 alleles has a clumpy phenotype. All three mutant alleles were selected for in low sucrose and were not selected for in low monosaccharide. Representative DIC images were taken with a 40× objective in a glass-bottomed, 96-well plate. Scale bars are 50 µm.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.022">http://dx.doi.org/10.7554/eLife.00367.022</ext-link></p></caption><graphic xlink:href="elife00367f008"/></fig><fig id="fig8s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.00367.023</object-id><label>Figure 8—figure supplement 1.</label><caption><title>Change in <italic>HXT4</italic> and <italic>SUC2</italic> expression in various allelic combinations of Recreated2 compared to the ancestor.</title><p>Change in RNA expression was measured using RT-qPCR using RNA isolated from a culture grown in 1 mM sucrose minimal media. Each point is the mean change in expression from three independent trials. The total length of the error bar is twice the size of the standard deviation. Ancestor alleles are shown in upper case and mutant alleles are shown in lower case. (<bold>A</bold>) Change in <italic>HXT4</italic> expression. Allelic combinations containing <italic>mth1</italic> have higher expression than those combinations with <italic>MTH1</italic>. (<bold>B</bold>) Change in <italic>SUC2</italic> expression. Note that the axes are different in the two plots.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.023">http://dx.doi.org/10.7554/eLife.00367.023</ext-link></p></caption><graphic xlink:href="elife00367fs008"/></fig></fig-group></p></sec><sec id="s2-9"><title>Some causal mutations show antagonistic pleiotropy</title><p>We studied the role of the causal mutations in different growth conditions. Do the mutations selected on low sucrose increase, decrease, or have no effect on fitness on other carbon sources? To ask this question, we selected spores from crosses between ancestral and recreated clones in low monosaccharide and high glucose. Several mutations were selected for in low sucrose but not in low monosaccharide: mutations in <italic>MTH1</italic>, <italic>CSE2</italic>, and <italic>ACE2</italic> in Recreated2; and in <italic>MCK1</italic>, <italic>IRC8</italic>, and <italic>GIN4</italic> in Recreated9. The primary selective force for four of these six mutations is likely to be their ability to produce multicellularity: an <italic>ace2</italic>Δ strain is morphologically similar to EvoClone2 and a <italic>mck1-G227Vfs249 irc8-G57V gin4-W19*</italic> strain is morphologically similar to EvoClone9 (<xref ref-type="fig" rid="fig8">Figure 8</xref>). Each of these two strains also outcompetes the ancestor strain in low sucrose but not in low monosaccharide (<xref ref-type="table" rid="tbl1">Table 1</xref>). This shows that multicellular clumps were specifically selected by growth in low sucrose, rather than being a general response to low sugar concentrations. But clumpiness alone was not sufficient to match the growth of the evolved strains in sucrose—the engineered clumpy strain (EngClumpy) was not as fit in 1 mM sucrose as one of the evolved clones that had similar clump size (EvoClone2) (<xref ref-type="table" rid="tbl1">Table 1</xref>).</p><p>Our evolved clones are less fit than their ancestor on high glucose (<xref ref-type="table" rid="tbl1">Table 1</xref>). There are two explanations for this observation. The first is a direct result of selection: some mutations that improve growth on low sucrose reduce fitness on high glucose. The second appeals to the large number of mutations that hitchhiked with the causal mutations: some of the hitchhikers are neutral on low sucrose, but slow cell growth on high glucose. We distinguished these two possibilities by competing the evolved and recreated strains against their ancestor in high glucose. For both EvoClone2 and EvoClone9, the recreated strains performed as well as or better than the evolved clones, but worse than the ancestor, suggesting that both causal alleles and hitchhiking alleles contribute to the evolved clones' reduced fitness in high glucose.</p><p>We identified the causal alleles that impair growth on high glucose through bulk segregant analysis. When the progeny of the cross between ancestor and recreated strains are grown on high glucose, one of the causal mutations (<italic>ubr1-C175Y</italic>) in Recreated2 is selected against, three are approximately neutral, and one (<italic>ira1-S2553C</italic>) is strongly selected for. For Recreated 9, one mutation (<italic>gin4-W19*</italic> ) is very strongly selected against, one (<italic>ubr1-L1050Yfs1063</italic>) is weakly selected against, three are roughly neutral, two are weakly selected for (<italic>mck1-G227Vfs249</italic> and <italic>are1(-10 G</italic>→<italic>T</italic>)), and one (<italic>mcd1</italic>-<italic>S175L</italic>) is very strongly selected for in high glucose (<xref ref-type="fig" rid="fig8">Figure 8</xref>). Note that the mutations in <italic>UBR1</italic> and <italic>GIN4</italic> have opposite effects in low monosaccharide: <italic>ubr1</italic> mutations are selected for in both low sucrose and low monosaccharide, whereas the <italic>gin4-W19*</italic> mutation is selected against in low monosaccharide.</p><p>To analyze the effect of individual mutations on expression of <italic>HXT4</italic> and <italic>SUC2</italic> in EvoClone2, we made seven strains, each with a different allelic combination of the five causal mutations. We then measured <italic>HXT4</italic> and <italic>SUC2</italic> expression using RT-qPCR in the seven new strains and in Recreated2. <xref ref-type="fig" rid="fig8s1">Figure 8—figure supplement 1</xref> shows the results. The increase in <italic>HXT4</italic> expression is clearly a result of the <italic>MTH1</italic> loss of function mutation (<italic>mth1-H154Tfs156</italic>): all combinations containing this allele have higher <italic>HXT4</italic> expression than those containing the ancestral <italic>MTH1</italic> allele. This result confirms the role of <italic>MTH1</italic> in hexose transporter repression. The gradient in <italic>SUC2</italic> expression across the allelic combinations indicates that <italic>SUC2</italic> expression is a complex phenotype, which is under the quantitative control alleles at several genes. We also checked the clumpiness of each strain by manual observation using a microscope: the strains with the wild type <italic>ACE2</italic> allele were not clumpy, while each strain with <italic>ace2-L323*</italic> was clumpy. This result confirms that the clumpiness in EvoClone2 is caused by the loss of function mutation in <italic>ACE2</italic>.</p></sec></sec><sec id="s3" sec-type="discussion"><title>Discussion</title><p>In earlier work, we showed that clumps of yeast could grow in concentrations of sucrose where single cells could not. We speculated that the sharing of public goods (in this case, the hydrolysis products of sucrose) could select for multicellularity. Here, we confirmed this hypothesis: 11 of 12 clones evolved in low sucrose formed multicellular clumps caused by cells failing to separate after cell division (<xref ref-type="fig" rid="fig3">Figure 3</xref> and <xref ref-type="fig" rid="fig3fs1 fig3s2 fig3s3">Figure 3—figure supplements 1–3</xref>). RNA sequencing revealed two additional strategies that most clones used: 10 elevated invertase expression, a strategy we had engineered, and 11 elevated hexose transporter expression, a strategy we had not engineered. We found no evidence of another strategy we had engineered: sucrose import. We cannot eliminate the possibility that the clones import some sucrose, but when we kept cells from secreting invertase, their fitness in low sucrose was severely reduced, showing that their growth still depends on invertase secretion. We identified the mutations that led to increased growth on low sucrose and showed that these were necessary and sufficient to recreate the evolved phenotype in two of the evolved clones. The genes that were mutated identified pathways that control cell growth and nutrient utilization.</p><p>The evolved cells clump because they fail to separate after cell division, rather than by flocculating (<xref ref-type="fig" rid="fig4">Figure 4</xref>). A loss of function mutation in <italic>ACE2</italic> was likely the primary contributor to the clumpy phenotype in at least 6 of the 11 multicellular clones. We were able to recreate the clumpy morphology of EvoClone9 with mutations in <italic>MCK1</italic>, <italic>IRC8</italic>, and <italic>GIN4</italic>; note that EvoClone3, one of the two other clumpy clones with wild type <italic>ACE2</italic>, also had a mutation in <italic>IRC8</italic>. Mutations in <italic>ACE2</italic> are likely to be frequent because inactivation of this gene simultaneously reduces the expression of multiple genes needed for cell separation; mutations in the genes that Ace2 regulates would be likely to have smaller effects. None of our clones became clumpy by restoring function to <italic>AMN1</italic>, the gene whose loss of function allele was selected for during laboratory domestication of budding yeast (<xref ref-type="bibr" rid="bib79">Yvert et al., 2003</xref>). This is not surprising since the target size for mutations that restore <italic>AMN1</italic> to its wild type function is much smaller than the target size for inactivating <italic>ACE2.</italic></p><p>Was the clumpy phenotype selected because it confers some other advantage, such as faster settling in unstirred cultures (<xref ref-type="bibr" rid="bib62">Ratcliff et al., 2012</xref>), rather than because it improves growth in low sucrose? Three lines of evidence show that this is not true. First, none of the evolved clones had a strong fitness advantage over the ancestor in 1 mM glucose plus 1 mM fructose, the monosaccharide equivalent of 1 mM sucrose (<xref ref-type="table" rid="tbl1">Table 1</xref>) showing that these clones were specifically adapted to sucrose. Second, a strain that was engineered to be clumpy (EngClumpy) outcompeted the ancestor in low sucrose, but not low monosaccharide (<xref ref-type="table" rid="tbl1">Table 1</xref>), demonstrating that low sucrose concentrations specifically selects for clumps, reinforcing our earlier findings on engineered multicellularity (<xref ref-type="bibr" rid="bib37">Koschwanez et al., 2011</xref>). Third, the alleles that recreate the clumpy morphology (<italic>ace2-L323*</italic> in EvoClone2 and the combination of <italic>mck1-G227Vfs249</italic>, <italic>gin4-W19*</italic>, and <italic>irc8-G57V</italic> in EvoClone9) are strongly selected for in low sucrose and are not selected for in low monosaccharide (<xref ref-type="fig" rid="fig8">Figure 8</xref> and <xref ref-type="table" rid="tbl1">Table 1</xref>). Because the starting pool of segregants and the selection protocol was identical in both media, this result shows that more efficient sucrose utilization, rather than some other phenotype, selected for the clumpy phenotype.</p><p>Analyzing gene expression and the role of individual mutations demonstrated that we had selected for more than multicellularity. Most (10 out of 12) of the evolved clones had increased invertase expression, even though the strain we had engineered to have high invertase expression had only a small advantage over a wild-type strain (<xref ref-type="table" rid="tbl1">Table 1</xref>). The advantage of increased invertase expression in a shaken culture of single cells is likely to be small since all cells share the monosaccharides that escape from the small, unstirred volume that surrounds each cell. But high invertase expression makes sense in conjunction with the other two strategies we observed, multicellularity and elevated monosaccharide import. In a clump, multiple cells share a larger unstirred volume and have preferential access to the monosaccharides their neighbors release (<xref ref-type="bibr" rid="bib37">Koschwanez et al., 2011</xref>). Increasing invertase expression will increase the rate of sucrose hydrolysis, and making more hexose transporters will increase the fraction of the monosaccharides that are imported rather than escaping to the bulk medium. Thus we expect that increasing invertase expression, making cell clumps, and making more hexose transporters will act together to increase fitness on low sucrose. Because low extracellular glucose induces the expression of both high affinity hexose transporters and invertase, it is likely that single mutations can increase levels of both types of proteins. In addition, clumpiness may indirectly elevate invertase expression: because secreted invertase expression peaks at roughly 0.5 mM external glucose (<xref ref-type="bibr" rid="bib37">Koschwanez et al., 2011</xref>), clumps may increase invertase expression because concentrating the cells in space collectively raises the local glucose concentration.</p><p>Many of the mutations we selected are likely to affect the linked processes of sugar harvesting and the control of cell growth. As the level of potential nutrients in their environment falls, cells must make choices: how much of their resources they invest in trying to harvest nutrients from their environment, and how much of the imported nutrients they store to meet future challenges and how much they use to maximize their current growth rate. As an example, wild strains may stop growing and dividing at low external glucose levels because this decision improves their ability to survive if glucose levels continued to fall, even though they could support cell growth at the current glucose concentration if they expressed very high levels of hexose transporters and stored no sugar as glycogen. Influences on this type of decision are likely to explain repeated mutations in one gene (<italic>UBR1</italic>) and three groups of genes, which control glucose's effect on transcription (the <italic>RGT1</italic> group), how nutrient availability controls cell growth (the <italic>IRA</italic> group), and general transcription (the Mediator group).</p><p><italic>UBR1</italic> was mutated in half of the evolved clones. Four of the six mutations are nonsense mutations, suggesting that we selected for inactivation of this protein. We suspect that prevention of N-end rule mediated protein degradation increases the expression of genes that promote growth in carbon limited cells: in both evolved clones that we recreated, the evolved allele of <italic>UBR1</italic> was selected for in sucrose and low monosaccharide and selected against in high glucose (<xref ref-type="fig" rid="fig8">Figure 8</xref>).</p><p>The <italic>RGT1</italic> group of genes gives external glucose concentrations the ability to control the expression of genes needed for growth in low glucose concentrations. Genes in this pathway were causally mutated in 8 of 12 evolved clones. The pattern of mutations we observe in these genes is consistent with selection to reduce the repression mediated by Rgt1 and Mth1. We also found, by using RT-qPCR on one of the reconstructed clones, that a loss-of-function mutation in <italic>MTH1</italic> was responsible for the increase in <italic>HXT4</italic> expression (<xref ref-type="fig" rid="fig8s1">Figure 8—figure supplement 1</xref>).</p><p>Mutations in one of the components of the Mediator complex, Sin4, had previously been shown to be involved in genes required for passage through Start, the transition in the yeast cell cycle that is most sensitive to nutritional signals. We recovered two mutations in this gene and found they were both associated with a strong increase in expression of Hxt3, a low affinity glucose transporter.</p><p>We saw mutations in <italic>IRA1</italic> and <italic>IRA2</italic> in four independent lineages. These mutations are likely to reduce the activity of Ira1 and Ira2 and thus increase the activity of Ras at low sugar concentrations, promoting growth. A mutant of <italic>IRA1</italic> in EvoClone2 was one of only two mutants of 13 in the recreated strains that was strongly selected for in all three conditions, and mutants of <italic>IRA1</italic> have been found in a yeast experimental evolution done in low glucose (<xref ref-type="bibr" rid="bib34">Kao and Sherlock, 2008</xref>) and rich media (Lang, personal communication, October 2012). We suspect that Ira1 and Ira2 reduce the activity of Ras and thus the growth rate of cells on suboptimal carbon sources (low sugar concentrations, ethanol). Because we used batch cultures, even the high glucose cultures went through an environmental shift after they had fermented all the glucose to ethanol. Disruption of either <italic>IRA1</italic> or <italic>IRA2</italic> increases sensitivity to heat shock (<xref ref-type="bibr" rid="bib71">Tanaka et al., 1989</xref>, <xref ref-type="bibr" rid="bib72">1990</xref>). We suggest that we have selected for mutants than ignore pathways designed to induce cells exposed to unfavorable conditions to reduce their growth rate and induce programs that will protect them if conditions worsen.</p><p>We took an integrated approach to understanding evolution. Experimental evolution has long been used to study nutrient adaptation in yeast (<xref ref-type="bibr" rid="bib58">Paquin and Adams, 1983</xref>; <xref ref-type="bibr" rid="bib26">Ferea et al., 1999</xref>) and bacteria (<xref ref-type="bibr" rid="bib43">Levin et al., 1977</xref>; <xref ref-type="bibr" rid="bib42">Lenski et al., 1991</xref>; <xref ref-type="bibr" rid="bib5">Barrick et al., 2009</xref>). Although evolved bacteriophage isolates were sequenced in the 1990s (<xref ref-type="bibr" rid="bib14">Bull et al., 1997</xref>), the cost of cellular genome sequencing has only recently dropped to a level where it is possible to find single nucleotide polymorphisms in an evolved microbe (<xref ref-type="bibr" rid="bib5">Barrick et al., 2009</xref>; <xref ref-type="bibr" rid="bib4">Araya et al., 2010</xref>). We combined engineering, experimental evolution, bulk segregant analysis, and genetic reconstruction to build a picture of how budding yeast evolve in response to a specific challenge: the need to grow on a nutrient that is released into the local environment by extracellular hydrolysis. Testing hypothetical strategies by engineering allowed us to look for their phenotypic signatures, such as the increased expression of hydrolytic enzymes, in evolved populations. Having a method to sift the minority of causal mutations from the majority of hitchhikers allowed us to use mutators and increase the speed of evolution. The ability to swap evolved and ancestral alleles made it possible to verify our methods for identifying causal mutations and show that we had identified mutations that were either generalists, promoting growth across all the environments we tested, or specialists showing antagonistic pleiotropy by being strongly beneficial in one environment and deleterious in another. Finally, performing multiple experiments in parallel allowed us to identify mutations in particular genes and pathways as contributing strongly to evolutionary success.</p><p>We speculate that unicellular organisms have evolved a combination of three strategies to access nutrients that are digested extracellularly: increasing secretion of an extracellular hydrolase, increasing the import of its products, and forming multicellular clumps. An increase in nutrient transport seems the least costly, but its effect as a sole strategy is limited if the local concentration of the hydrolysis products remains unchanged. An increase in enzyme secretion will increase the level of the hydrolysis products, but enzyme production is costly and allows unrelated cells from the same or different species to ‘cheat' and steal simple nutrients without incurring the cost. The formation of multicellular clumps allows the hydrolysis products to be shared among direct kin and reduces the risk of cheating (<xref ref-type="bibr" rid="bib37">Koschwanez et al., 2011</xref>). We have shown that yeast, when selected for growth in low sucrose media, used a mixture of all three strategies.</p><p>The combination of laboratory engineering and evolutionary analysis shown here and in our previous work supports the speculation that the sharing of public goods was the initial selection for undifferentiated multicellularity. For organisms with rigid cell walls like yeast, algae, and bacteria, cells must secrete enzymes to separate mother from daughter after cytokinesis. The target size for mutations in these enzymes and the control circuits that regulate their production will be large, making it easy to evolve undifferentiated multicellularity. Although lab strains of yeast were selected to be unicellular during their domestication (<xref ref-type="bibr" rid="bib54">Mortimer and Johnston, 1986</xref>), the level of clumpiness varies in wild isolates of yeast, and undifferentiated multicellularity may have been recurrently selected for its advantage in harvesting nutrients or for protection from stress (<xref ref-type="bibr" rid="bib68">Smukalla et al., 2008</xref>). Organisms without a cell wall, such as animal cells, would need another innovation, such as cell-adhering cadherins (<xref ref-type="bibr" rid="bib1">Abedin and King, 2008</xref>), in order to evolve multicellular clumps. But even in these organisms, cells must adhere to each other during sexual reproduction, implying the existence of adhesive proteins that could be modified to allow the evolution of multicellular clumps.</p><p>If undifferentiated multicellularity is so easy to evolve in eukaryotes, then why did it take so long for fossils of differentiated multicellular organisms to appear in the fossil record? Indeed, there is suggestive evidence of undifferentiated eukaryotic multicellularity 2 billion years ago (<xref ref-type="bibr" rid="bib24">El Albani et al. 2010</xref>), but the radiation of more complex multicellularity is not seen for over a billion years later. This disparity has received a large amount of attention in paleontology, particularly in the context of the seemingly rapid radiation of the animals at the base of the Cambrian. Possible explanations for this ‘Cambrian explosion' include novel selective forces due to the environmental changes that occurred shortly before the Cambrian or the possibility that evolutionary arms races among a few early species drove the radiation (<xref ref-type="bibr" rid="bib49">Marshall, 2006</xref>; <xref ref-type="bibr" rid="bib60">Peters and Gaines, 2012</xref>). Alternatively, the evolution of differentiated multicellularity may have been slow and gradual and we simply lack the fossil evidence. Whatever the case, it is interesting to ask why the first differentiation might have arisen. From our data, we speculate that the effectiveness of clumping as a strategy could be improved if a limited number of cells increase their enzyme secretion. Such a simple division of labor may have been the precursor to complex multicellularity in the lineages that lead to animals and other complex multicellular species.</p></sec><sec id="s4" sec-type="materials|methods"><title>Materials and methods</title><sec id="s4-1"><title>Media</title><p>Unless otherwise noted, the experimental evolution and all experiments were conducted in minimal (no amino acids or nucleotides) synthetic media with the following two exceptions: (1) the auxotrophic strains expressing high levels of Flo1 required the addition of leucine, histidine, and lysine, and (2) strains undergoing galactose induction in <xref ref-type="fig" rid="fig2">Figure 2</xref> were pregrown in YP 2% glycerol (10 g/l yeast extract, 20 g/l peptone, and 2% [v/v] glycerol) plus the indicated concentration of galactose. The minimal synthetic media was made with refrigerated sugar stocks and refrigerated 10× yeast nitrogen base (YNB) that was based on Wickerman's recipe (<xref ref-type="bibr" rid="bib77">Wickerham, 1951</xref>) but without riboflavin, folic acid, and inositol (see [<xref ref-type="bibr" rid="bib37">Koschwanez et al., 2011</xref>] for the recipe). 200 units/ml Penicillin and 0.2 mg/ml streptavidin were added to the experimental evolution media to prevent bacterial contamination. Unless otherwise noted, all chemicals used in this research were purchased from Sigma-Aldrich (<ext-link ext-link-type="uri" xlink:href="http://www.sigmasldrich.com">www.sigmasldrich.com</ext-link>).</p></sec><sec id="s4-2"><title>Experimental evolution</title><p>Each of the 10 evolved cultures was derived from a haploid <italic>MAT</italic><bold>a</bold> W303 prototrophic strain that expressed a yellow fluorescent protein (YFP) variant from the constitutive <italic>ACT1</italic> promoter. See <xref ref-type="supplementary-material" rid="SD3-data">supplementary file 3</xref> for a description of the starting strains and their construction. To start the evolution, a single colony from each strain was inoculated into 80 mM glucose media and grown for 7–20 generations, and 5×10<sup>5</sup> cells were then inoculated into a 125 ml flask containing 50 ml of 1 mM sucrose media. The flasks were placed at 30°C on a shaker rotating at 120 rpm. When there was visible growth in a culture (>10<sup>6</sup> cells/ml), it was passaged as follows. First, the culture was spun down and resuspended in 1 ml of yeast nitrogen base (YNB) to ensure that non-sucrose carbon sources, created by invertase activity or cellular metabolism, were not inoculated into the new cultures. Second, a sample of resuspended cells were mixed with glycerol and frozen at −80°C. Third, the cell density was measured with a Coulter counter (<ext-link ext-link-type="uri" xlink:href="http://www.beckmancoultermedia.com">www.beckmancoultermedia.com</ext-link>). Finally, the resuspended cells were vortexed and 5×10<sup>5</sup> cells were inoculated into a new 125 ml flask containing 50 ml of warm 1 mM sucrose media. When the cells became too clumpy to count accurately in the Coulter counter, the cultures were instead diluted at a 500:1 ratio into fresh media.</p></sec><sec id="s4-3"><title>Clonal selection</title><p>To pick clones, eight replicates of each culture were serial diluted 2:1 twelve times in 1 mM sucrose media in 96-well plates such that the last dilution contained an estimated average of less than one cell per well for each replicate. For each population, there were thus eight parallel sets of dilutions. The plates were incubated on a plate shaker at 30°C for 3 days, and each well was examined using a microscope for growth in sucrose and for morphology. A well was assumed to contain a clonal population at the lowest starting cell density where growth was seen if there was a lower density where growth was not seen in any of the eight wells. Each of the selected clones was verified to grow from low density in 1 mM sucrose within 2 days.</p></sec><sec id="s4-4"><title>Bulk segregant analysis</title><p><italic>URA3</italic> was deleted from each YFP-labeled <italic>MAT</italic><bold>a</bold> evolved clone. The <italic>ura3Δ</italic> strain was mated with yJHK519, a <italic>MAT</italic>α, <italic>trp1Δ::kanMX4 ura3::P</italic><sub><italic>STE2</italic></sub><italic>-URA3</italic> derivative of the ancestor labeled with CFP. In this strain, the endogenous <italic>URA3</italic> promoter is replaced with the <italic>STE2</italic> promoter, which is only induced in <italic>MAT</italic><bold>a</bold> cells, making it possible to select for <italic>MAT</italic><bold>a</bold> spores after meiosis. Mating was performed by mixing cells from the two strains together on a YPD plate with a toothpick and growing overnight at 30°C. The mating mixtures were then plated onto G418, -TRP plates, and a diploid strain was selected from a colony on the plate.</p><p>To sporulate the diploid strains, cultures were grown to saturation in YPD, and then diluted 1:50 into YP 2% acetate. The cells were grown in acetate for 12–24 hr and then pelleted and resuspended in 2% acetate. After 4–5 days of incubation on a roller drum at 25°C, sporulation was verified by observing the cells in a microscope.</p><p>To digest ascii, 1 ml of the sporulated culture was pelleted and resuspended in 50 µl 10% Zymolyase (<ext-link ext-link-type="uri" xlink:href="http://www.zymoresearch.com">www.zymoresearch.com</ext-link>) for 1 hr at 30°C. 400 µl of water and 50 µl of 10% Triton X-100 were added, and the digested spores were sonicated for 5–10 s to separate the tetrads. Tetrad separation was verified by observing the cells in a microscope. The separated tetrads were then spun down slowly (6000 rpm) and resuspended in growth media.</p><p>To select haploid spores, the entire digested spore culture was added to 50 ml of 1 mM glucose + 1 mM fructose minimal media and growth to saturation. This media selected for <italic>TRP1,</italic> haploid <italic>MAT</italic><bold>a</bold> cells: neither haploid <italic>MAT</italic>α nor diploid <italic>MAT</italic><bold>a</bold><italic>/MAT</italic>α cells can express <italic>URA3</italic> from the <italic>STE2</italic> promoter. Each culture was selected for growth in 1 mM sucrose through four passages (inoculation to saturated growth) as follows: one 100:1 dilution, and three 500:1 dilutions (≈34 generations total). The cells were spun down and resuspended in YNB during each dilution in order to eliminate the non-sucrose carbon sources from the media. Genomic DNA was made from the final saturated culture.</p><p>Bulk segregant analysis of Recreated2 and Recreated9 followed the same procedure with one exception: the cells were not spun down in between dilutions. Instead, 100 µl of culture was added directly to 50 ml of fresh media for a 500:1 dilution.</p></sec><sec id="s4-5"><title>Genomic DNA preparation</title><p>To prepare genomic DNA, the culture was pelleted and resuspended in 50 µl of 1% Zymolyase in 0.1 M, pH 8.0 NaEDTA. The cells were incubated for 30 min at 37°C to digest the cell wall, and then the cells were lysed by adding 50 µl 0.2 M NaEDTA, 0.4 M pH 8.0 Tris, 2% SDS and incubated at 65°C for 30 min. 63 µl of 5 M potassium acetate was added, and the mixture was incubated for 30 min on ice. The insoluble residue was then pelleted, and 750 µl of ice-cold ethanol was added to 300 µl of the supernatant to precipitate the DNA. The DNA was pelleted, and the pellet was resuspended in 0.2 mg/ml RNAase A. After 1 hr of incubation at 37°C, 2 µl of 20 mg/ml Proteinase K was added, and the solution was incubated for an additional 2 hr at 37°C. The DNA was again precipitated by adding 130 µl isopropanol. The DNA was pelleted, briefly washed with 70% ethanol, repelleted, and resuspended in 100 µl 10 mM Tris, pH 8.0.</p></sec><sec id="s4-6"><title>RNA isolation</title><p>To isolate RNA, cells in log-phase growth were fixed by adding 6 ml of culture to 9 ml of ice-cold methanol and then incubating at −20°C for 10 min. Cells were pelleted at 4°C, resuspended in 1 ml of RNAase-free ice-cold water in a 2 ml cryogenic storage vial, and then repelleted at 4°C. RNA was then isolated with acidic phenol using a published protocol (<xref ref-type="bibr" rid="bib16">Collart and Oliviero, 2001</xref>).</p></sec><sec id="s4-7"><title>RT-qPCR</title><p>Isolated RNA was treated with DNase I (Thermo Scientific EN0525, <ext-link ext-link-type="uri" xlink:href="http://www.thermofisher.com">www.thermofisher.com</ext-link>) according to manufacturer's instructions. Intact RNA was verified by observing two sharp rRNA bands using agarose gel electrophoresis. cDNA was made from 100–200 ng of RNA using Thermo Scientific Maxima First Strand cDNA Synthesis Kit for RT-qPCR (K1641). The 20 µl reaction was diluted tenfold and 11 µl of the diluted sample was used for real time qPCR using Thermo Scientific Maxima SYBR Green/ROX qPCR Master Mix (K0222). Amplification efficiency of the primers was verified by generating standard curves with four serial dilutions. The following primers were used: <italic>SUC2</italic> (GCTTTCCTTTTC-CTTTTGGCTGG, TCATTCATCCAGCCCTTGTTGG); <italic>HXT4</italic> (TTGGGTTACTGTACAAACTACG, TGTCA-TACCACCAATCATAAAC); <italic>ALG9</italic> (GTTTAATCCGGGCTGGTTCCAT, TAGACCCAGTGGACAGATAGCG). <italic>ALG9</italic> was used for normalization (<xref ref-type="bibr" rid="bib73">Teste et al., 2009</xref>). The ΔΔC<sub>T</sub> method was used to find change in RNA expression (<xref ref-type="bibr" rid="bib47">Livak and Schmittgen, 2001</xref>). Three independent trials were performed for each reported RT-qPCR result and the data reported is the mean difference in expression between two strains.</p></sec><sec id="s4-8"><title>Sequencing</title><p>DNA and RNA libraries were prepared for sequencing using the Illumina TruSEQ kit (<ext-link ext-link-type="uri" xlink:href="http://www.illumina.com">www.illumina.com</ext-link>) and were sequenced on an Illumina HiSeq 2000. Mean coverage across the genome was as follows: ancestor DNA—56×, evolved clone DNA—10× minimum, sporulated pools—90× minimum, RNA—40× minimum. Single end, 50 bp reads were used for the ancestor DNA, evolved clone DNA, and RNA. Paired end, 100 bp reads were used for the sporulated pool DNA. Sequencing data is deposited in the Sequence Read Archive.</p></sec><sec id="s4-9"><title>Sequence analysis</title><p>DNA sequences were aligned to the S288C reference genome r64 (downloaded from the Saccharomyces Genome Database, <ext-link ext-link-type="uri" xlink:href="http://www.yeastgenome.org">www.yeastgenome.org</ext-link>) using the Burrows-Wheeler Aligner (<ext-link ext-link-type="uri" xlink:href="http://bio-bwa.sourceforge.net">bio-bwa.sourceforge.net</ext-link>) (<xref ref-type="bibr" rid="bib44">Li and Durbin, 2009</xref>). The resulting SAM (Sequence Alignment/Map) file was converted to a BAM (binary SAM) file, sorted, indexed, and made into a pileup format file using the samtools software (<ext-link ext-link-type="uri" xlink:href="http://samtools.sourceforge.net">samtools.sourceforge.net</ext-link>) (<xref ref-type="bibr" rid="bib45">Li et al., 2009</xref>). Indels were realigned locally using GATK (<ext-link ext-link-type="uri" xlink:href="http://www.broadinstitute.org/gatk/">www.broadinstitute.org/gatk/</ext-link>) (<xref ref-type="bibr" rid="bib51">McKenna et al., 2010</xref>), and variants were called from the pileup file using the Varscan software (<ext-link ext-link-type="uri" xlink:href="http://varscan.sourceforge.net">varscan.sourceforge.net</ext-link>) (<xref ref-type="bibr" rid="bib36">Koboldt et al., 2012</xref>).</p><p>To perform the segregation analysis, we wrote a custom sequencing pipeline in Python (<ext-link ext-link-type="uri" xlink:href="http://www.python.org">www.python.org</ext-link>), using the Biopython (<ext-link ext-link-type="uri" xlink:href="http://biopython.org">biopython.org</ext-link>) and pysam (<ext-link ext-link-type="uri" xlink:href="http://code.google.com/p/pysam/">code.google.com/p/pysam/</ext-link>) modules, that finds sequence variants between the ancestor and clone, classifies each variant as a nonsynonymous coding region, synonymous coding region, or promoter mutation, and ranks each mutation by its segregation frequency. All software written for this analysis is publicly available at <ext-link ext-link-type="uri" xlink:href="https://github.com/koschwanez">https://github.com/koschwanez</ext-link>. We used the following criteria to select putative causal mutations: (1) the evolved clone was mutated relative to the strain used to inoculate the time zero culture, (2) the mutation caused a nonsynonymous substitution in a coding region or changed the promoter sequence ≤500 bp upstream of the coding region start site, and (3) the mutation appeared in over 90% of the reads in backcrossed pool. Putative causal mutations were manually verified by looking at aligned reads.</p><p>RNA sequences were aligned to the S288C reference genome using TopHat (<xref ref-type="bibr" rid="bib74">Trapnell et al., 2009</xref>), and significant differences in expression between the ancestor and the evolved clone were called using the default setting in Cufflinks (<xref ref-type="bibr" rid="bib75">Trapnell et al., 2010</xref>). The Cufflinks package uses the log of the ratio of expression in two conditions together with an estimate of the gene's variance to generate a t-value that it uses in a Students t-test. We used one ancestor derivative, yJHK110, as the reference strain to compare expression with the evolved clones, and the other ancestor derivative, yJHK111, as a control strain in the comparison.</p></sec><sec id="s4-10"><title>Data analysis and figures</title><p>Unless otherwise noted, data analysis was performed in the R programming language (<ext-link ext-link-type="uri" xlink:href="http://www.r-project.org">www.r-project.org</ext-link>) and plots were generated using the R library ggplot2 (<xref ref-type="bibr" rid="bib78">Wickham, 2009</xref>). The Adjusted Wald method of calculating 95% binomial confidence intervals (<xref ref-type="bibr" rid="bib2">Agresti and Coull, 1998</xref>) was used in <xref ref-type="fig" rid="fig2">Figure 2</xref> because a low number of samples were used to generate a binomial mean. Cartoons were made using Adobe Illustrator CS5 (<ext-link ext-link-type="uri" xlink:href="http://www.adobe.com">www.adobe.com</ext-link>).</p></sec><sec id="s4-11"><title>Sanger trace plot analysis</title><p>Commercial Sanger sequencing returns trace plots, chromatograms that indicate the relative frequency of each base at each position in the sequenced DNA. Trace plots were used in two analyses. First, we estimated the spread of mutations through a population by Sanger sequencing time points from the evolution. Second, we estimated the segregation of alleles in the backcrossed, recreated strains by Sanger sequencing the putative causal alleles in the final, selected population. The fraction of mutant alleles in the population was assumed to be the height of the mutant allele peak divided by the height of the mutant allele peak plus the ancestor allele peak. In the time course analysis, values below 5% (the approximate background level) were assumed to be zero, and values above 95% were assumed to be 100%. In the segregation analysis, actual values, not corrected for background, are shown.</p><p>We note that many of the mutations appear to rise with very similar time courses to each other. This could reflect the modest time resolution of our measurements and the insensitivity of using Sanger sequencing to estimate allele frequencies, but we suspect it reflects the strong advantages that accrue to those rare lineages where two or three beneficial mutations occur in quick succession (<xref ref-type="bibr" rid="bib39">Lang et al., 2011</xref>).</p></sec><sec id="s4-12"><title>Competitions</title><p>To compare the fitness of the recreated and evolved strains, we developed a fitness assay that was suitable for strains with non-uniform morphologies. We could not use standard, quantitative methods of measuring fitness, such as FACS-based competitions (<xref ref-type="bibr" rid="bib19">Desai et al., 2007</xref>), because we were not able to accurately count the number of cells per clump in a large population. We therefore used a qualitative microscopy-based fitness assay to compare growth of a YFP-labeled population and a CFP-labeled population. We started the competitions by growing each population separately to saturation in 1 mM glucose + 1 mM fructose. Equal volumes of each strain were mixed together, and then checked under a microscope to verify that both strains were equally present. 50 µl of the mixed culture was then inoculated into 10 ml of the test media. When the culture reached visible density (>10<sup>6</sup> cells/ml), three blindly chosen images of a 20 µl sample were taken on the microscope with a 20× objective. Extermination was defined as zero cells in one strain and more than 30 total cells of the other strain. If both strains were still present, 50 µl of the culture was inoculated into 10 ml of fresh media. After six growth cycles, a winner was declared if there was a clear majority (>75%) of one strain. Otherwise, the competition was declared a tie. The number of cycles until extermination was averaged across three independent experiments, and qualitative values (++++ through −−−−) were assigned as described in <xref ref-type="table" rid="tbl1">Table 1</xref>. + or − was assigned if one of the strains won two out of the three independent competitions that lasted through all six growth cycles.</p></sec><sec id="s4-13"><title>Check for flocculation or cell separation</title><p>The evolved clones already constitutively expressed YFP, and a constitutively expressed CFP version of each clone was made by deleting the P<sub><italic>ACT1</italic></sub>-ymCitrine-<italic>URA3</italic> construct from the <italic>URA3</italic> locus and then transforming with a <italic>ura3</italic>Δ::P<sub><italic>ACT1</italic></sub>-yCerulean-<italic>URA3</italic> construct. The two versions of each strain were grown separately in 1 mM glucose + 1 mM fructose media to saturation, and then the two versions were vortexed or sonicated (both gave identical results) for 10 s and then combined at low density in fresh 1 mM sucrose media. The mixed cultures were grown to visible density at 30°C and then pipetted into a well of a glass-bottomed 96-well plate and examined by microscope. In each of the three independent experiments, at least 50 clumps were checked and three representative images were taken for each strain.</p></sec><sec id="s4-14"><title>Microscopy</title><p>Images except those in <xref ref-type="fig" rid="fig7">Figure 7</xref> were taken in a 96-well glass-bottomed plate (greiner bio-one, <ext-link ext-link-type="uri" xlink:href="http://www.gbo.com">www.gbo.com</ext-link>) on a Nikon Ti inverted microscope (<ext-link ext-link-type="uri" xlink:href="http://www.nikoninstruments.com">www.nikoninstruments.com</ext-link>) with a Photometrics CoolSnap HQ camera (<ext-link ext-link-type="uri" xlink:href="http://www.photometrics.com">www.photometrics.com</ext-link>) and MetaMorph software (<ext-link ext-link-type="uri" xlink:href="http://www.metamorph.com">www.metamorph.com</ext-link>). <xref ref-type="fig" rid="fig7">Figure 7</xref> images were taken on a Nikon Ti inverted microscope with a Yokogawa spinning disc confocal unit (<ext-link ext-link-type="uri" xlink:href="http://www.yokogawa.com">www.yokogawa.com</ext-link>), 447, 515, and 594 nm lasers (<ext-link ext-link-type="uri" xlink:href="http://www.spectraloptics.com">www.spectraloptics.com</ext-link>), a Hamamatsu Orca camera (<ext-link ext-link-type="uri" xlink:href="http://www.hamamatsu.com">www.hamamatsu.com</ext-link>), and MetaMorph software. Images were converted to 8-bit, projected, adjusted for contrast, and annotated with scale bars using the Fiji distribution of ImageJ (<xref ref-type="bibr" rid="bib66">Schindelin et al., 2012</xref>). Contrast was changed for visibility only; our results do not depend on illumination levels and are not affected by the change in contrast.</p></sec><sec id="s4-15"><title>Strain reconstruction and reversion</title><p>Alleles were replaced in both the evolved and ancestor strains by transforming with a <italic>URA3</italic> plasmid that contained a portion of the targeted gene with the desired mutation and was digested at a single cut site within the target gene (<xref ref-type="bibr" rid="bib64">Rothstein, 1991</xref>). See <xref ref-type="fig" rid="fig5s1">Figure 5—figure supplement 1</xref> for the complete protocol. Transformants were selected on—URA plates, colonies that grew were streaked out on—URA plates, and the target sequence was amplified and sequenced to verify insertion. The colony was grown overnight in YPD (10 g/l yeast extract, 20 g/l peptone, and 2% [w/v] dextrose) and selected on 5FOA plates to identify cells in which the <italic>URA3</italic> construct had looped out (<xref ref-type="bibr" rid="bib8">Boeke et al., 1984</xref>). Colonies that grew were streaked out on 5FOA, replica plated to YP 2% Acetate (20 g/l agar, 10 g/l yeast extract, 20 g/l peptone, and 2% [w/v] potassium acetate) to eliminate petites, and the target region was amplified and sequenced to ensure that only the desired allele was present. Restriction enzymes, DNA polymerase, polynucleotide kinase, and ligation enzymes were purchased from New England Biolabs (<ext-link ext-link-type="uri" xlink:href="http://www.neb.com">www.neb.com</ext-link>).</p></sec></sec></body><back><ack id="ack"><title>Acknowledgements</title><p>The authors thank Erin O'Shea, Vlad Denic, Mike Laub, Greg Lang, Jim Bull, Rama Ranganathan, and members of the Murray Lab for reviewing the manuscript.</p></ack><sec sec-type="additional-information"><title>Additional information</title><fn-group content-type="competing-interest"><title>Competing interests</title><fn fn-type="conflict" id="conf1"><p>The authors declare that no competing interests exist.</p></fn></fn-group><fn-group content-type="author-contribution"><title>Author contributions</title><fn fn-type="con" id="con1"><p>JHK, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con2"><p>KRF, Conception and design, Drafting or revising the article</p></fn><fn fn-type="con" id="con3"><p>AWM, Conception and design, Analysis and interpretation of data, Drafting or revising the article</p></fn></fn-group></sec><sec sec-type="supplementary-material"><title>Additional files</title><supplementary-material id="SD1-data"><object-id pub-id-type="doi">10.7554/eLife.00367.024</object-id><label>Supplementary file 1.</label><caption><p>(<bold>A</bold>) Genes whose expression changed significantly in nine or more evolved clones. (<bold>B</bold>) Genes whose expression was increased or reduced at least tenfold. (<bold>C</bold>) All genes whose expression changed significantly in any clone or the control strain. Values are log<sub>2</sub> change from ancestor to sample strain. Genes whose value did not change significantly are shown as zero.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.024">http://dx.doi.org/10.7554/eLife.00367.024</ext-link></p></caption><media mime-subtype="xlsx" mimetype="application" xlink:href="elife00367s001.xlsx"/></supplementary-material><supplementary-material id="SD2-data"><object-id pub-id-type="doi">10.7554/eLife.00367.025</object-id><label>Supplementary file 2.</label><caption><p>Putative causal mutation pathway summary.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.025">http://dx.doi.org/10.7554/eLife.00367.025</ext-link></p></caption><media mime-subtype="doc" mimetype="application" xlink:href="elife00367s002.doc"/></supplementary-material><supplementary-material id="SD3-data"><object-id pub-id-type="doi">10.7554/eLife.00367.026</object-id><label>Supplementary file 3.</label><caption><p>Strain list.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.00367.026">http://dx.doi.org/10.7554/eLife.00367.026</ext-link></p></caption><media mime-subtype="doc" mimetype="application" xlink:href="elife00367s003.doc"/></supplementary-material></sec><sec sec-type="supplementary-material"><sec sec-type="datasets"><title>Major datasets</title><p>The following dataset was generated:</p><p><related-object content-type="generated-dataset" document-id="Dataset ID and/or url" document-id-type="dataset" document-type="data" id="dataro1"><name><surname>Foster</surname><given-names>KR</given-names></name>, <name><surname>Murray</surname><given-names>AW</given-names></name>, <name><surname>Koschwanez</surname><given-names>JH</given-names></name>, <year>2013</year><x>, </x><source>Sequences for: Improved use of a public good selects for the evolution of undifferentiated multicellularity</source><x>, </x><object-id pub-id-type="art-access-id">PRJNA192270</object-id><x>; </x><ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/bioproject/192270">http://www.ncbi.nlm.nih.gov/bioproject/192270</ext-link><x>, </x><comment>Publicly available at the Sequence Read Archive (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/sra">http://www.ncbi.nlm.nih.gov/sra</ext-link>).</comment></related-object></p><p>The following previously published dataset was used:</p><p><related-object content-type="generated-dataset" document-id="Dataset ID and/or url 2" document-id-type="dataset" document-type="data" id="dataro2"><name><surname>Cherry</surname><given-names>JM</given-names></name>, <name><surname>Hong</surname><given-names>EL</given-names></name>, <name><surname>Amundsen</surname><given-names>C</given-names></name>, <name><surname>Balakrishnan</surname><given-names>R</given-names></name>, <name><surname>Binkley</surname><given-names>G</given-names></name>, <name><surname>Chan</surname><given-names>ET</given-names></name>, <name><surname>Christie</surname><given-names>KR</given-names></name>, <name><surname>Costanzo</surname><given-names>MC</given-names></name>, <name><surname>Dwight</surname><given-names>SS</given-names></name>, <name><surname>Engel</surname><given-names>SR</given-names></name>, <name><surname>Fisk</surname><given-names>DG</given-names></name>, <name><surname>Hirschman</surname><given-names>JE</given-names></name>, <name><surname>Hitz</surname><given-names>BC</given-names></name>, <name><surname>Karra</surname><given-names>K</given-names></name>, <name><surname>Krieger</surname><given-names>CJ</given-names></name>, <name><surname>Miyasato</surname><given-names>SR</given-names></name>, <name><surname>Nash</surname><given-names>RS</given-names></name>, <name><surname>Park</surname><given-names>J</given-names></name>, <name><surname>Skrzypek</surname><given-names>MS</given-names></name>, <name><surname>Simison</surname><given-names>M</given-names></name>, <name><surname>Weng</surname><given-names>S</given-names></name>, <name><surname>Wong</surname><given-names>ED</given-names></name>, <year>2012</year><x>, </x><source>Saccharomyces Genome Database: the genomics resource of budding yeast</source><x>, </x><object-id pub-id-type="art-access-id">R64-1-1</object-id><x>; </x><ext-link ext-link-type="uri" xlink:href="http://downloads.yeastgenome.org/sequence/S288C_reference/genome_releases/">http://downloads.yeastgenome.org/sequence/S288C_reference/genome_releases/</ext-link><x>, </x><comment>Publicly available at <ext-link ext-link-type="uri" xlink:href="http://www.yeastgenome.org/">http://www.yeastgenome.org/</ext-link>.</comment></related-object></p></sec></sec><ref-list><title>References</title><ref id="bib1"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Abedin</surname><given-names>M</given-names></name><name><surname>King</surname><given-names>N</given-names></name></person-group><year>2008</year><article-title>The premetazoan ancestry of cadherins</article-title><source>Science</source><volume>319</volume><fpage>946</fpage><lpage>8</lpage><pub-id pub-id-type="doi">10.1126/science.1151084</pub-id></element-citation></ref><ref id="bib2"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Agresti</surname><given-names>A</given-names></name><name><surname>Coull</surname><given-names>B</given-names></name></person-group><year>1998</year><article-title>Approximate is better than ‘exact' for interval estimation of binomial proportions</article-title><source>Am Stat</source><volume>52</volume><fpage>119</fpage><lpage>26</lpage></element-citation></ref><ref id="bib3"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Alegado</surname><given-names>RA</given-names></name><name><surname>Brown</surname><given-names>LW</given-names></name><name><surname>Cao</surname><given-names>S</given-names></name><name><surname>Dermenjian</surname><given-names>RK</given-names></name><name><surname>Zuzow</surname><given-names>R</given-names></name><name><surname>Fairclough</surname><given-names>SR</given-names></name><etal/></person-group><year>2012</year><article-title>A bacterial sulfonolipid triggers multicellular development in the closest living relatives of animals</article-title><source>elife</source><volume>1</volume><fpage>e00013</fpage><pub-id pub-id-type="doi">10.7554/eLife.00013</pub-id></element-citation></ref><ref id="bib4"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Araya</surname><given-names>CL</given-names></name><name><surname>Payen</surname><given-names>C</given-names></name><name><surname>Dunham</surname><given-names>MJ</given-names></name><name><surname>Fields</surname><given-names>S</given-names></name></person-group><year>2010</year><article-title>Whole-genome sequencing of a laboratory-evolved yeast strain</article-title><source>BMC Genomics</source><volume>11</volume><fpage>88</fpage><pub-id pub-id-type="doi">10.1186/1471-2164-11-88</pub-id></element-citation></ref><ref id="bib5"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Barrick</surname><given-names>J</given-names></name><name><surname>Yu</surname><given-names>D</given-names></name><name><surname>Yoon</surname><given-names>S</given-names></name><name><surname>Jeong</surname><given-names>H</given-names></name><name><surname>Oh</surname><given-names>T</given-names></name><name><surname>Schneider</surname><given-names>D</given-names></name><etal/></person-group><year>2009</year><article-title>Genome evolution and adaptation in a long-term experiment with Escherichia coli</article-title><source>Nature</source><volume>461</volume><fpage>1243</fpage><lpage>7</lpage><pub-id pub-id-type="doi">10.1038/nature08480</pub-id></element-citation></ref><ref id="bib6"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bartel</surname><given-names>B</given-names></name><name><surname>Wünning</surname><given-names>I</given-names></name><name><surname>Varshavsky</surname><given-names>A</given-names></name></person-group><year>1990</year><article-title>The recognition component of the N-end rule pathway</article-title><source>EMBO J</source><volume>9</volume><fpage>3179</fpage><lpage>89</lpage></element-citation></ref><ref id="bib7"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Birkeland</surname><given-names>SR</given-names></name><name><surname>Jin</surname><given-names>N</given-names></name><name><surname>Ozdemir</surname><given-names>AC</given-names></name><name><surname>Lyons</surname><given-names>RH</given-names></name><name><surname>Weisman</surname><given-names>LS</given-names></name><name><surname>Wilson</surname><given-names>TE</given-names></name></person-group><year>2010</year><article-title>Discovery of mutations in <italic>Saccharomyces cerevisiae</italic> by pooled linkage analysis and whole genome sequencing</article-title><source>Genetics</source><pub-id pub-id-type="doi">10.1534/genetics.110.123232</pub-id></element-citation></ref><ref id="bib8"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Boeke</surname><given-names>JD</given-names></name><name><surname>LaCroute</surname><given-names>F</given-names></name><name><surname>Fink</surname><given-names>GR</given-names></name></person-group><year>1984</year><article-title>A positive selection for mutants lacking orotidine-5′-phosphate decarboxylase activity in yeast: 5-fluoro-orotic acid resistance</article-title><source>Mol Gen Genet</source><volume>197</volume><fpage>345</fpage><lpage>6</lpage><pub-id pub-id-type="doi">10.1007/BF00330984</pub-id></element-citation></ref><ref id="bib9"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Boraas</surname><given-names>M</given-names></name><name><surname>Seale</surname><given-names>D</given-names></name><name><surname>Boxhorn</surname><given-names>J</given-names></name></person-group><year>1998</year><article-title>Phagotrophy by a flagellate selects for colonial prey: a possible origin of multicellularity</article-title><source>Evol Ecol</source><volume>12</volume><fpage>153</fpage><lpage>64</lpage><pub-id pub-id-type="doi">10.1023/A:1006527528063</pub-id></element-citation></ref><ref id="bib10"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Brauer</surname><given-names>MJ</given-names></name><name><surname>Christianson</surname><given-names>CM</given-names></name><name><surname>Pai</surname><given-names>DA</given-names></name><name><surname>Dunham</surname><given-names>MJ</given-names></name></person-group><year>2006</year><article-title>Mapping novel traits by array-assisted bulk segregant analysis in <italic>Saccharomyces cerevisiae</italic></article-title><source>Genetics</source><volume>173</volume><fpage>1813</fpage><lpage>6</lpage><pub-id pub-id-type="doi">10.1534/genetics.106.057927</pub-id></element-citation></ref><ref id="bib11"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Broach</surname><given-names>JR</given-names></name></person-group><year>2012</year><article-title>Nutritional control of growth and development in yeast</article-title><source>Genetics</source><volume>192</volume><fpage>73</fpage><lpage>105</lpage><pub-id pub-id-type="doi">10.1534/genetics.111.135731</pub-id></element-citation></ref><ref id="bib12"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Brown</surname><given-names>CA</given-names></name><name><surname>Murray</surname><given-names>AW</given-names></name><name><surname>Verstrepen</surname><given-names>KJ</given-names></name></person-group><year>2010</year><article-title>Rapid expansion and functional divergence of subtelomeric gene families in yeasts</article-title><source>Curr Biol</source><volume>20</volume><fpage>895</fpage><lpage>903</lpage><pub-id pub-id-type="doi">10.1016/j.cub.2010.04.027</pub-id></element-citation></ref><ref id="bib13"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Brown</surname><given-names>CJ</given-names></name><name><surname>Todd</surname><given-names>KM</given-names></name><name><surname>Rosenzweig</surname><given-names>RF</given-names></name></person-group><year>1998</year><article-title>Multiple duplications of yeast hexose transport genes in response to selection in a glucose-limited environment</article-title><source>Mol Biol Evol</source><volume>15</volume><fpage>931</fpage><lpage>42</lpage><pub-id pub-id-type="doi">10.1093/oxfordjournals.molbev.a026009</pub-id></element-citation></ref><ref id="bib14"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bull</surname><given-names>JJ</given-names></name><name><surname>Badgett</surname><given-names>MR</given-names></name><name><surname>Wichman</surname><given-names>HA</given-names></name><name><surname>Huelsenbeck</surname><given-names>JP</given-names></name><name><surname>Hillis</surname><given-names>DM</given-names></name><name><surname>Gulati</surname><given-names>A</given-names></name><etal/></person-group><year>1997</year><article-title>Exceptional convergent evolution in a virus</article-title><source>Genetics</source><volume>147</volume><fpage>1497</fpage><lpage>507</lpage></element-citation></ref><ref id="bib15"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Carlson</surname><given-names>M</given-names></name><name><surname>Osmond</surname><given-names>BC</given-names></name><name><surname>Botstein</surname><given-names>D</given-names></name></person-group><year>1981</year><article-title>Mutants of yeast defective in sucrose utilization</article-title><source>Genetics</source><volume>98</volume><fpage>25</fpage><lpage>40</lpage></element-citation></ref><ref id="bib16"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Collart</surname><given-names>MA</given-names></name><name><surname>Oliviero</surname><given-names>S</given-names></name></person-group><year>2001</year><article-title>Preparation of Yeast RNA</article-title><source>Curr Protoc Mol Biol</source><volume>23</volume><fpage>13.12.1</fpage><lpage>13.12.5</lpage><pub-id pub-id-type="doi">10.1002/0471142727.mb1312s23</pub-id></element-citation></ref><ref id="bib17"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Colman-Lerner</surname><given-names>A</given-names></name><name><surname>Chin</surname><given-names>TE</given-names></name><name><surname>Brent</surname><given-names>R</given-names></name></person-group><year>2001</year><article-title>Yeast Cbk1 and Mob2 activate daughter-specific genetic programs to induce asymmetric cell fates</article-title><source>Cell</source><volume>107</volume><fpage>739</fpage><lpage>50</lpage><pub-id pub-id-type="doi">10.1016/S0092-8674(01)00596-7</pub-id></element-citation></ref><ref id="bib18"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>den Dunnen</surname><given-names>JT</given-names></name><name><surname>Antonarakis</surname><given-names>SE</given-names></name></person-group><year>2000</year><article-title>Mutation nomenclature extensions and suggestions to describe complex mutations: a discussion</article-title><source>Hum Mutat</source><volume>15</volume><fpage>7</fpage><lpage>12</lpage><pub-id pub-id-type="doi">10.1002/(SICI)1098-1004(200001)15:1<7::AID-HUMU4>3.0.CO;2-N</pub-id></element-citation></ref><ref id="bib19"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Desai</surname><given-names>MM</given-names></name><name><surname>Fisher</surname><given-names>DS</given-names></name><name><surname>Murray</surname><given-names>AW</given-names></name></person-group><year>2007</year><article-title>The speed of evolution and maintenance of variation in asexual populations</article-title><source>Curr Biol</source><volume>17</volume><fpage>385</fpage><lpage>94</lpage><pub-id pub-id-type="doi">10.1016/j.cub.2007.01.072</pub-id></element-citation></ref><ref id="bib20"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dodyk</surname><given-names>F</given-names></name><name><surname>Rothstein</surname><given-names>A</given-names></name></person-group><year>1964</year><article-title>Factors influencing the appearance of invertase in <italic>Saccharomyces cerevisiae</italic></article-title><source>Arch Biochem Biophys</source><volume>104</volume><fpage>478</fpage><lpage>86</lpage><pub-id pub-id-type="doi">10.1016/0003-9861(64)90492-8</pub-id></element-citation></ref><ref id="bib21"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dohmen</surname><given-names>RJ</given-names></name><name><surname>Madura</surname><given-names>K</given-names></name><name><surname>Bartel</surname><given-names>B</given-names></name><name><surname>Varshavsky</surname><given-names>A</given-names></name></person-group><year>1991</year><article-title>The N-end rule is mediated by the <italic>UBC2(RAD6)</italic> ubiquitin-conjugating enzyme</article-title><source>Proc Natl Acad Sci USA</source><volume>88</volume><fpage>7351</fpage><lpage>5</lpage><pub-id pub-id-type="doi">10.1073/pnas.88.16.7351</pub-id></element-citation></ref><ref id="bib22"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dohrmann</surname><given-names>PR</given-names></name><name><surname>Butler</surname><given-names>G</given-names></name><name><surname>Tamai</surname><given-names>K</given-names></name><name><surname>Dorland</surname><given-names>S</given-names></name><name><surname>Greene</surname><given-names>JR</given-names></name><name><surname>Thiele</surname><given-names>DJ</given-names></name><etal/></person-group><year>1992</year><article-title>Parallel pathways of gene regulation: homologous regulators <italic>SWI5</italic> and <italic>ACE2</italic> differentially control transcription of HO and chitinase</article-title><source>Genes Dev</source><volume>6</volume><fpage>93</fpage><lpage>104</lpage><pub-id pub-id-type="doi">10.1101/gad.6.1.93</pub-id></element-citation></ref><ref id="bib23"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dworkin</surname><given-names>M</given-names></name></person-group><year>1972</year><article-title>The myxobacteria: new directions in studies of procaryotic development</article-title><source>CRC Crit Rev Microbiol</source><volume>2</volume><fpage>435</fpage><lpage>52</lpage><pub-id pub-id-type="doi">10.3109/10408417209103873</pub-id></element-citation></ref><ref id="bib24"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>El Albani</surname><given-names>A</given-names></name><name><surname>Bengtson</surname><given-names>S</given-names></name><name><surname>Canfield</surname><given-names>DE</given-names></name><name><surname>Bekker</surname><given-names>A</given-names></name><name><surname>Macchiarelli</surname><given-names>R</given-names></name><name><surname>Mazurier</surname><given-names>A</given-names></name><etal/></person-group><year>2010</year><article-title>Large colonial organisms with coordinated growth in oxygenated environments 2.1 Gyr ago</article-title><source>Nature</source><volume>466</volume><fpage>100</fpage><lpage>4</lpage><pub-id pub-id-type="doi">10.1038/nature09166</pub-id></element-citation></ref><ref id="bib25"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Esmon</surname><given-names>PC</given-names></name><name><surname>Esmon</surname><given-names>BE</given-names></name><name><surname>Schauer</surname><given-names>IE</given-names></name><name><surname>Taylor</surname><given-names>A</given-names></name><name><surname>Schekman</surname><given-names>R</given-names></name></person-group><year>1987</year><article-title>Structure, assembly, and secretion of octameric invertase</article-title><source>J Biol Chem</source><volume>262</volume><fpage>4387</fpage><lpage>94</lpage></element-citation></ref><ref id="bib26"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ferea</surname><given-names>TL</given-names></name><name><surname>Botstein</surname><given-names>D</given-names></name><name><surname>Brown</surname><given-names>PO</given-names></name><name><surname>Rosenzweig</surname><given-names>RF</given-names></name></person-group><year>1999</year><article-title>Systematic changes in gene expression patterns following adaptive evolution in yeast</article-title><source>Proc Natl Acad Sci USA</source><volume>96</volume><fpage>9721</fpage><lpage>6</lpage><pub-id pub-id-type="doi">10.1073/pnas.96.17.9721</pub-id></element-citation></ref><ref id="bib27"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gresham</surname><given-names>D</given-names></name><name><surname>Desai</surname><given-names>MM</given-names></name><name><surname>Tucker</surname><given-names>CM</given-names></name><name><surname>Jenq</surname><given-names>HT</given-names></name><name><surname>Pai</surname><given-names>DA</given-names></name><name><surname>Ward</surname><given-names>A</given-names></name><etal/></person-group><year>2008</year><article-title>The repertoire and dynamics of evolutionary adaptations to controlled nutrient-limited environments in yeast</article-title><source>PLoS Genet</source><volume>4</volume><fpage>e1000303</fpage><pub-id pub-id-type="doi">10.1371/journal.pgen.1000303</pub-id></element-citation></ref><ref id="bib28"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Grosberg</surname><given-names>RK</given-names></name><name><surname>Strathmann</surname><given-names>RR</given-names></name></person-group><year>2007</year><article-title>The evolution of multicellularity: a minor major transition?</article-title><source>Annu Rev Ecol Evol Syst</source><volume>38</volume><fpage>621</fpage><lpage>54</lpage><pub-id pub-id-type="doi">10.1146/annurev.ecolsys.36.102403.114735</pub-id></element-citation></ref><ref id="bib29"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Guo</surname><given-names>B</given-names></name><name><surname>Styles</surname><given-names>CA</given-names></name><name><surname>Feng</surname><given-names>Q</given-names></name><name><surname>Fink</surname><given-names>GR</given-names></name></person-group><year>2000</year><article-title>A <italic>Saccharomyces</italic> gene family involved in invasive growth, cell-cell adhesion, and mating</article-title><source>Proc Natl Acad Sci USA</source><volume>97</volume><fpage>12158</fpage><lpage>63</lpage><pub-id pub-id-type="doi">10.1073/pnas.220420397</pub-id></element-citation></ref><ref id="bib30"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ingolia</surname><given-names>NT</given-names></name><name><surname>Murray</surname><given-names>AW</given-names></name></person-group><year>2007</year><article-title>Positive-feedback loops as a flexible biological module</article-title><source>Curr Biol</source><volume>17</volume><fpage>668</fpage><lpage>77</lpage><pub-id pub-id-type="doi">10.1016/j.cub.2007.03.016</pub-id></element-citation></ref><ref id="bib31"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Iyer</surname><given-names>VR</given-names></name><name><surname>Horak</surname><given-names>CE</given-names></name><name><surname>Scafe</surname><given-names>CS</given-names></name><name><surname>Botstein</surname><given-names>D</given-names></name><name><surname>Snyder</surname><given-names>M</given-names></name><name><surname>Brown</surname><given-names>PO</given-names></name></person-group><year>2001</year><article-title>Genomic binding sites of the yeast cell-cycle transcription factors SBF and MBF</article-title><source>Nature</source><volume>409</volume><fpage>533</fpage><lpage>8</lpage><pub-id pub-id-type="doi">10.1038/35054095</pub-id></element-citation></ref><ref id="bib32"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jin</surname><given-names>YH</given-names></name><name><surname>Garg</surname><given-names>P</given-names></name><name><surname>Stith</surname><given-names>CMW</given-names></name><name><surname>Al-Refai</surname><given-names>H</given-names></name><name><surname>Sterling</surname><given-names>JF</given-names></name><name><surname>Murray</surname><given-names>LJW</given-names></name><etal/></person-group><year>2005</year><article-title>The multiple biological roles of the 3′→5′ exonuclease of Saccharomyces cerevisiae DNA polymerase delta require switching between the polymerase and exonuclease domains</article-title><source>Mol Cell Biol</source><volume>25</volume><fpage>461</fpage><lpage>71</lpage><pub-id pub-id-type="doi">10.1128/MCB.25.1.461-471.2005</pub-id></element-citation></ref><ref id="bib33"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kaiser</surname><given-names>CA</given-names></name><name><surname>Botstein</surname><given-names>D</given-names></name></person-group><year>1986</year><article-title>Secretion-defective mutations in the signal sequence for <italic>Saccharomyces cerevisiae</italic> invertase</article-title><source>Mol Cell Biol</source><volume>6</volume><fpage>2382</fpage><lpage>91</lpage></element-citation></ref><ref id="bib34"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kao</surname><given-names>KC</given-names></name><name><surname>Sherlock</surname><given-names>G</given-names></name></person-group><year>2008</year><article-title>Molecular characterization of clonal interference during adaptive evolution in asexual populations of <italic>Saccharomyces cerevisiae</italic></article-title><source>Nat Genet</source><volume>40</volume><fpage>1499</fpage><lpage>504</lpage><pub-id pub-id-type="doi">10.1038/ng.280</pub-id></element-citation></ref><ref id="bib35"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kessin</surname><given-names>RH</given-names></name><name><surname>Gundersen</surname><given-names>GG</given-names></name><name><surname>Zaydfudim</surname><given-names>V</given-names></name><name><surname>Grimson</surname><given-names>M</given-names></name><name><surname>Blanton</surname><given-names>RL</given-names></name></person-group><year>1996</year><article-title>How cellular slime molds evade nematodes</article-title><source>Proc Natl Acad Sci USA</source><volume>93</volume><fpage>4857</fpage><lpage>61</lpage><pub-id pub-id-type="doi">10.1073/pnas.93.10.4857</pub-id></element-citation></ref><ref id="bib36"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Koboldt</surname><given-names>DC</given-names></name><name><surname>Zhang</surname><given-names>Q</given-names></name><name><surname>Larson</surname><given-names>DE</given-names></name><name><surname>Shen</surname><given-names>D</given-names></name><name><surname>McLellan</surname><given-names>MD</given-names></name><name><surname>Lin</surname><given-names>L</given-names></name><etal/></person-group><year>2012</year><article-title>VarScan 2: somatic mutation and copy number alteration discovery in cancer by exome sequencing</article-title><source>Genome Res</source><volume>22</volume><fpage>568</fpage><lpage>76</lpage><pub-id pub-id-type="doi">10.1101/gr.129684.111</pub-id></element-citation></ref><ref id="bib37"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Koschwanez</surname><given-names>JH</given-names></name><name><surname>Foster</surname><given-names>KR</given-names></name><name><surname>Murray</surname><given-names>AW</given-names></name></person-group><year>2011</year><article-title>Sucrose utilization in budding yeast as a model for the origin of undifferentiated multicellularity</article-title><source>Plos Biol</source><volume>9</volume><fpage>e1001122</fpage><pub-id pub-id-type="doi">10.1371/journal.pbio.1001122</pub-id></element-citation></ref><ref id="bib38"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kvitek</surname><given-names>DJ</given-names></name><name><surname>Sherlock</surname><given-names>G</given-names></name></person-group><year>2011</year><article-title>Reciprocal sign epistasis between frequently experimentally evolved adaptive mutations causes a rugged fitness landscape</article-title><source>PLoS Genet</source><volume>7</volume><fpage>e1002056</fpage><pub-id pub-id-type="doi">10.1371/journal.pgen.1002056</pub-id></element-citation></ref><ref id="bib39"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lang</surname><given-names>GI</given-names></name><name><surname>Botstein</surname><given-names>D</given-names></name><name><surname>Desai</surname><given-names>MM</given-names></name></person-group><year>2011</year><article-title>Genetic variation and the fate of beneficial mutations in asexual populations</article-title><source>Genetics</source><volume>188</volume><fpage>647</fpage><lpage>61</lpage><pub-id pub-id-type="doi">10.1534/genetics.111.128942</pub-id></element-citation></ref><ref id="bib40"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lang</surname><given-names>GI</given-names></name><name><surname>Murray</surname><given-names>AW</given-names></name></person-group><year>2008</year><article-title>Estimating the per-base-pair mutation rate in the yeast <italic>Saccharomyces cerevisiae</italic></article-title><source>Genetics</source><volume>178</volume><fpage>67</fpage><lpage>82</lpage><pub-id pub-id-type="doi">10.1534/genetics.107.071506</pub-id></element-citation></ref><ref id="bib41"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lenski</surname><given-names>R</given-names></name></person-group><year>1988</year><article-title>Experimental studies of pleiotropy and epistasis in <italic>Escherichia coli</italic>. I. Variation in competitive fitness among mutants resistant to virus T4</article-title><source>Evolution</source><volume>42</volume><fpage>425</fpage><lpage>32</lpage><pub-id pub-id-type="doi">10.2307/2409028</pub-id></element-citation></ref><ref id="bib42"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lenski</surname><given-names>RE</given-names></name><name><surname>Rose</surname><given-names>MR</given-names></name><name><surname>Simpson</surname><given-names>SC</given-names></name><name><surname>Tadler</surname><given-names>SC</given-names></name></person-group><year>1991</year><article-title>Long-term experimental evolution in <italic>Escherichia coli</italic>. I. Adaptation and divergence during 2,000 generations</article-title><source>Am Nat</source><volume>138</volume><fpage>1315</fpage><lpage>41</lpage><pub-id pub-id-type="doi">10.1086/285289</pub-id></element-citation></ref><ref id="bib43"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Levin</surname><given-names>BR</given-names></name><name><surname>Stewart</surname><given-names>FM</given-names></name><name><surname>Chao</surname><given-names>L</given-names></name></person-group><year>1977</year><article-title>Resource-limited growth, competition, and predation: A model and experimental studies with bacteria and bacteriophage</article-title><source>Am Nat</source><volume>111</volume><fpage>3</fpage><lpage>24</lpage><pub-id pub-id-type="doi">10.1016/j.mbs.2007.02.005</pub-id></element-citation></ref><ref id="bib44"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname><given-names>H</given-names></name><name><surname>Durbin</surname><given-names>R</given-names></name></person-group><year>2009</year><article-title>Fast and accurate short read alignment with Burrows-Wheeler transform</article-title><source>Bioinformatics</source><volume>25</volume><fpage>1754</fpage><lpage>60</lpage><pub-id pub-id-type="doi">10.1093/bioinformatics/btp324</pub-id></element-citation></ref><ref id="bib45"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname><given-names>H</given-names></name><name><surname>Handsaker</surname><given-names>B</given-names></name><name><surname>Wysoker</surname><given-names>A</given-names></name><name><surname>Fennell</surname><given-names>T</given-names></name><name><surname>Ruan</surname><given-names>J</given-names></name><name><surname>Homer</surname><given-names>N</given-names></name><etal/></person-group><year>2009</year><article-title>The sequence alignment/map format and SAMtools</article-title><source>Bioinformatics</source><volume>25</volume><fpage>2078</fpage><lpage>9</lpage><pub-id pub-id-type="doi">10.1093/bioinformatics/btp352</pub-id></element-citation></ref><ref id="bib46"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname><given-names>L</given-names></name><name><surname>Quinton</surname><given-names>T</given-names></name><name><surname>Miles</surname><given-names>S</given-names></name><name><surname>Breeden</surname><given-names>LL</given-names></name></person-group><year>2005</year><article-title>Genetic interactions between mediator and the late G<sub>1</sub>-specific transcription factor Swi6 in <italic>Saccharomyces cerevisiae</italic></article-title><source>Genetics</source><volume>171</volume><fpage>477</fpage><lpage>88</lpage><pub-id pub-id-type="doi">10.1534/genetics.105.043893</pub-id></element-citation></ref><ref id="bib47"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Livak</surname><given-names>KJ</given-names></name><name><surname>Schmittgen</surname><given-names>TD</given-names></name></person-group><year>2001</year><article-title>Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method</article-title><source>Methods</source><volume>25</volume><fpage>402</fpage><lpage>8</lpage><pub-id pub-id-type="doi">10.1006/meth.2001.1262</pub-id></element-citation></ref><ref id="bib48"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Magwene</surname><given-names>PM</given-names></name><name><surname>Willis</surname><given-names>JH</given-names></name><name><surname>Kelly</surname><given-names>JK</given-names></name></person-group><year>2011</year><article-title>The statistics of bulk segregant analysis using next generation sequencing</article-title><source>PLoS Comput Biol</source><volume>7</volume><fpage>e1002255</fpage><pub-id pub-id-type="doi">10.1371/journal.pcbi.1002255</pub-id></element-citation></ref><ref id="bib49"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Marshall</surname><given-names>CR</given-names></name></person-group><year>2006</year><article-title>Explaining the Cambrian “Explosion” of animals</article-title><source>Annu Rev Earth Planet Sci</source><volume>34</volume><fpage>355</fpage><lpage>84</lpage><pub-id pub-id-type="doi">10.1146/ annurev.earth.33.031504.103001</pub-id></element-citation></ref><ref id="bib50"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Matsumoto</surname><given-names>K</given-names></name><name><surname>Uno</surname><given-names>I</given-names></name><name><surname>Ishikawa</surname><given-names>T</given-names></name></person-group><year>1983</year><article-title>Control of cell division in <italic>Saccharomyces cerevisiae</italic> mutants defective in adenylate cyclase and cAMP-dependent protein kinase</article-title><source>Exp Cell Res</source><volume>146</volume><fpage>151</fpage><lpage>61</lpage><pub-id pub-id-type="doi">10.1016/0014-4827(83)90333-6</pub-id></element-citation></ref><ref id="bib51"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>McKenna</surname><given-names>A</given-names></name><name><surname>Hanna</surname><given-names>M</given-names></name><name><surname>Banks</surname><given-names>E</given-names></name><name><surname>Sivachenko</surname><given-names>A</given-names></name><name><surname>Cibulskis</surname><given-names>K</given-names></name><name><surname>Kernytsky</surname><given-names>A</given-names></name><etal/></person-group><year>2010</year><article-title>The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data</article-title><source>Genome Res</source><volume>20</volume><fpage>1297</fpage><lpage>303</lpage><pub-id pub-id-type="doi">10.1101/gr.107524.110</pub-id></element-citation></ref><ref id="bib52"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Meijer</surname><given-names>MM</given-names></name><name><surname>Boonstra</surname><given-names>J</given-names></name><name><surname>Verkleij</surname><given-names>AJ</given-names></name><name><surname>Verrips</surname><given-names>CT</given-names></name></person-group><year>1996</year><article-title>Kinetic analysis of hexose uptake in <italic>Saccharomyces cerevisiae</italic> cultivated in continuous culture</article-title><source>Biochim Biophys Acta</source><volume>1277</volume><fpage>209</fpage><lpage>16</lpage><pub-id pub-id-type="doi">10.1016/S0005-2728(96)00098-9</pub-id></element-citation></ref><ref id="bib53"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Michelmore</surname><given-names>RW</given-names></name><name><surname>Paran</surname><given-names>I</given-names></name><name><surname>Kesseli</surname><given-names>RV</given-names></name></person-group><year>1991</year><article-title>Identification of markers linked to disease-resistance genes by bulked segregant analysis: a rapid method to detect markers in specific genomic regions by using segregating populations</article-title><source>Proc Natl Acad Sci USA</source><volume>88</volume><fpage>9828</fpage><lpage>32</lpage><pub-id pub-id-type="doi">10.1073/pnas.88.21.9828</pub-id></element-citation></ref><ref id="bib54"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Mortimer</surname><given-names>RK</given-names></name><name><surname>Johnston</surname><given-names>JR</given-names></name></person-group><year>1986</year><article-title>Genealogy of principal strains of the yeast genetic stock center</article-title><source>Genetics</source><volume>113</volume><fpage>35</fpage><lpage>43</lpage></element-citation></ref><ref id="bib55"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Myers</surname><given-names>LC</given-names></name><name><surname>Kornberg</surname><given-names>RD</given-names></name></person-group><year>2000</year><article-title>Mediator of transcriptional regulation</article-title><source>Annu Rev Biochem</source><volume>69</volume><fpage>729</fpage><lpage>49</lpage><pub-id pub-id-type="doi">10.1146/annurev.biochem.69.1.729</pub-id></element-citation></ref><ref id="bib56"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ozcan</surname><given-names>S</given-names></name><name><surname>Dover</surname><given-names>J</given-names></name><name><surname>Rosenwald</surname><given-names>AG</given-names></name><name><surname>Wölfl</surname><given-names>S</given-names></name><name><surname>Johnston</surname><given-names>M</given-names></name></person-group><year>1996</year><article-title>Two glucose transporters in Saccharomyces cerevisiae are glucose sensors that generate a signal for induction of gene expression</article-title><source>Proc Natl Acad Sci USA</source><volume>93</volume><fpage>12428</fpage><lpage>32</lpage><pub-id pub-id-type="doi">10.1073/pnas.93.22.12428</pub-id></element-citation></ref><ref id="bib57"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ozcan</surname><given-names>S</given-names></name><name><surname>Johnston</surname><given-names>M</given-names></name></person-group><year>1999</year><article-title>Function and regulation of yeast hexose transporters</article-title><source>Microbiol Mol Biol Rev</source><volume>63</volume><fpage>554</fpage><lpage>69</lpage></element-citation></ref><ref id="bib58"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Paquin</surname><given-names>C</given-names></name><name><surname>Adams</surname><given-names>J</given-names></name></person-group><year>1983</year><article-title>Frequency of fixation of adaptive mutations is higher in evolving diploid than haploid yeast populations</article-title><source>Nature</source><volume>302</volume><fpage>495</fpage><lpage>500</lpage><pub-id pub-id-type="doi">10.1038/302495a0</pub-id></element-citation></ref><ref id="bib59"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Perlman</surname><given-names>D</given-names></name><name><surname>Raney</surname><given-names>P</given-names></name><name><surname>Halvorson</surname><given-names>HO</given-names></name></person-group><year>1986</year><article-title>Mutations affecting the signal sequence alter synthesis and secretion of yeast invertase</article-title><source>Proc Natl Acad Sci USA</source><volume>83</volume><fpage>5033</fpage><lpage>7</lpage><pub-id pub-id-type="doi">10.1073/pnas.83.14.5033</pub-id></element-citation></ref><ref id="bib60"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Peters</surname><given-names>SE</given-names></name><name><surname>Gaines</surname><given-names>RR</given-names></name></person-group><year>2012</year><article-title>Formation of the ‘Great Unconformity' as a trigger for the Cambrian explosion</article-title><source>Nature</source><volume>484</volume><fpage>363</fpage><lpage>6</lpage><pub-id pub-id-type="doi">10.1038/nature10969</pub-id></element-citation></ref><ref id="bib61"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pfeiffer</surname><given-names>T</given-names></name><name><surname>Bonhoeffer</surname><given-names>S</given-names></name></person-group><year>2003</year><article-title>An evolutionary scenario for the transition to undifferentiated multicellularity</article-title><source>Proc Natl Acad Sci USA</source><volume>100</volume><fpage>1095</fpage><lpage>8</lpage><pub-id pub-id-type="doi">10.1073/pnas.0335420100</pub-id></element-citation></ref><ref id="bib62"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ratcliff</surname><given-names>WC</given-names></name><name><surname>Denison</surname><given-names>RF</given-names></name><name><surname>Borrello</surname><given-names>M</given-names></name><name><surname>Travisano</surname><given-names>M</given-names></name></person-group><year>2012</year><article-title>Experimental evolution of multicellularity</article-title><source>Proc Natl Acad Sci USA</source><volume>109</volume><fpage>1595</fpage><lpage>600</lpage><pub-id pub-id-type="doi">10.1073/pnas.1115323109</pub-id></element-citation></ref><ref id="bib63"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Reifenberger</surname><given-names>E</given-names></name><name><surname>Boles</surname><given-names>E</given-names></name><name><surname>Ciriacy</surname><given-names>M</given-names></name></person-group><year>1997</year><article-title>Kinetic characterization of individual hexose transporters of <italic>Saccharomyces cerevisiae</italic> and their relation to the triggering mechanisms of glucose repression</article-title><source>Eur J Biochem</source><volume>245</volume><fpage>324</fpage><lpage>33</lpage><pub-id pub-id-type="doi">10.1111/j.1432-1033.1997.00324.x</pub-id></element-citation></ref><ref id="bib64"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rothstein</surname><given-names>R</given-names></name></person-group><year>1991</year><article-title>Targeting, disruption, replacement, and allele rescue: integrative DNA transformation in yeast</article-title><source>Meth Enzymol</source><volume>194</volume><fpage>281</fpage><lpage>301</lpage><pub-id pub-id-type="doi">10.1016/0076-6879(91)94022-5</pub-id></element-citation></ref><ref id="bib65"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sbia</surname><given-names>M</given-names></name><name><surname>Parnell</surname><given-names>EJ</given-names></name><name><surname>Yu</surname><given-names>Y</given-names></name><name><surname>Olsen</surname><given-names>AE</given-names></name><name><surname>Kretschmann</surname><given-names>KL</given-names></name><name><surname>Voth</surname><given-names>WP</given-names></name><etal/></person-group><year>2008</year><article-title>Regulation of the yeast Ace2 transcription factor during the cell cycle</article-title><source>J Biol Chem</source><volume>283</volume><fpage>11135</fpage><lpage>45</lpage><pub-id pub-id-type="doi">10.1074/jbc.M800196200</pub-id></element-citation></ref><ref id="bib66"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schindelin</surname><given-names>J</given-names></name><name><surname>Arganda-Carreras</surname><given-names>I</given-names></name><name><surname>Frise</surname><given-names>E</given-names></name><name><surname>Kaynig</surname><given-names>V</given-names></name><name><surname>Longair</surname><given-names>M</given-names></name><name><surname>Pietzsch</surname><given-names>T</given-names></name><etal/></person-group><year>2012</year><article-title>Fiji: an open-source platform for biological-image analysis</article-title><source>Nat Meth</source><volume>9</volume><fpage>676</fpage><lpage>82</lpage><pub-id pub-id-type="doi">10.1038/nmeth.2019</pub-id></element-citation></ref><ref id="bib67"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Segrè</surname><given-names>AV</given-names></name><name><surname>Murray</surname><given-names>AW</given-names></name><name><surname>Leu</surname><given-names>J-Y</given-names></name></person-group><year>2006</year><article-title>High-resolution mutation mapping reveals parallel experimental evolution in yeast</article-title><source>Plos Biol</source><volume>4</volume><fpage>e256</fpage><pub-id pub-id-type="doi">10.1371/journal.pbio.0040256</pub-id></element-citation></ref><ref id="bib68"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Smukalla</surname><given-names>S</given-names></name><name><surname>Caldara</surname><given-names>M</given-names></name><name><surname>Pochet</surname><given-names>N</given-names></name><name><surname>Beauvais</surname><given-names>A</given-names></name><name><surname>Guadagnini</surname><given-names>S</given-names></name><name><surname>Yan</surname><given-names>C</given-names></name><etal/></person-group><year>2008</year><article-title>FLO1 is a variable green beard gene that drives biofilm-like cooperation in budding yeast</article-title><source>Cell</source><volume>135</volume><fpage>726</fpage><lpage>37</lpage><pub-id pub-id-type="doi">10.1016/j.cell.2008.09.037</pub-id></element-citation></ref><ref id="bib69"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Stambuk</surname><given-names>BU</given-names></name><name><surname>da Silva</surname><given-names>MA</given-names></name><name><surname>Panek</surname><given-names>AD</given-names></name><name><surname>de Araujo</surname><given-names>PS</given-names></name></person-group><year>1999</year><article-title>Active α-glucoside transport in <italic>Saccharomyces cerevisiae</italic></article-title><source>FEMS Microbiol Lett</source><volume>170</volume><fpage>105</fpage><lpage>10</lpage></element-citation></ref><ref id="bib70"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tammi</surname><given-names>M</given-names></name><name><surname>Ballou</surname><given-names>L</given-names></name><name><surname>Taylor</surname><given-names>A</given-names></name><name><surname>Ballou</surname><given-names>CE</given-names></name></person-group><year>1987</year><article-title>Effect of glycosylation on yeast invertase oligomer stability</article-title><source>J Biol Chem</source><volume>262</volume><fpage>4395</fpage><lpage>401</lpage></element-citation></ref><ref id="bib71"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tanaka</surname><given-names>K</given-names></name><name><surname>Matsumoto</surname><given-names>K</given-names></name><name><surname>Toh-E</surname><given-names>A</given-names></name></person-group><year>1989</year><article-title>IRA1, an inhibitory regulator of the RAS-cyclic AMP pathway in <italic>Saccharomyces cerevisiae</italic></article-title><source>Mol Cell Biol</source><volume>9</volume><fpage>757</fpage><lpage>68</lpage></element-citation></ref><ref id="bib72"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tanaka</surname><given-names>K</given-names></name><name><surname>Nakafuku</surname><given-names>M</given-names></name><name><surname>Tamanoi</surname><given-names>F</given-names></name><name><surname>Kaziro</surname><given-names>Y</given-names></name><name><surname>Matsumoto</surname><given-names>K</given-names></name><name><surname>Toh-e</surname><given-names>A</given-names></name></person-group><year>1990</year><article-title>IRA2, a second gene of <italic>Saccharomyces cerevisiae</italic> that encodes a protein with a domain homologous to mammalian <italic>ras</italic> GTPase-activating protein</article-title><source>Mol Cell Biol</source><volume>10</volume><fpage>4303</fpage><lpage>13</lpage></element-citation></ref><ref id="bib73"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Teste</surname><given-names>M-A</given-names></name><name><surname>Duquenne</surname><given-names>M</given-names></name><name><surname>François</surname><given-names>JM</given-names></name><name><surname>Parrou</surname><given-names>J-L</given-names></name></person-group><year>2009</year><article-title>Validation of reference genes for quantitative expression analysis by real-time RT-PCR in <italic>Saccharomyces cerevisiae</italic></article-title><source>BMC Mol Biol</source><volume>10</volume><fpage>99</fpage><pub-id pub-id-type="doi">10.1186/1471-2199-10-99</pub-id></element-citation></ref><ref id="bib74"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Trapnell</surname><given-names>C</given-names></name><name><surname>Pachter</surname><given-names>L</given-names></name><name><surname>Salzberg</surname><given-names>SL</given-names></name></person-group><year>2009</year><article-title>TopHat: discovering splice junctions with RNA-Seq</article-title><source>Bioinformatics</source><volume>25</volume><fpage>1105</fpage><lpage>11</lpage><pub-id pub-id-type="doi">10.1093/bioinformatics/btp120</pub-id></element-citation></ref><ref id="bib75"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Trapnell</surname><given-names>C</given-names></name><name><surname>Williams</surname><given-names>BA</given-names></name><name><surname>Pertea</surname><given-names>G</given-names></name><name><surname>Mortazavi</surname><given-names>A</given-names></name><name><surname>Kwan</surname><given-names>G</given-names></name><name><surname>van Baren</surname><given-names>MJ</given-names></name><etal/></person-group><year>2010</year><article-title>Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation</article-title><source>Nat Biotechnol</source><volume>28</volume><fpage>511</fpage><lpage>5</lpage><pub-id pub-id-type="doi">10.1038/nbt.1621</pub-id></element-citation></ref><ref id="bib76"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wenger</surname><given-names>JW</given-names></name><name><surname>Piotrowski</surname><given-names>J</given-names></name><name><surname>Nagarajan</surname><given-names>S</given-names></name><name><surname>Chiotti</surname><given-names>K</given-names></name><name><surname>Sherlock</surname><given-names>G</given-names></name><name><surname>Rosenzweig</surname><given-names>F</given-names></name></person-group><year>2011</year><article-title>Hunger artists: yeast adapted to carbon limitation show trade-offs under carbon sufficiency</article-title><source>PLoS Genet</source><volume>7</volume><fpage>e1002202</fpage><pub-id pub-id-type="doi">10.1371/journal.pgen.1002202</pub-id></element-citation></ref><ref id="bib77"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wickerham</surname><given-names>LJ</given-names></name></person-group><year>1951</year><article-title>Taxonomy of yeasts</article-title><source>USDA Technical Bulletin</source><volume>1029</volume><fpage>1</fpage><lpage>56</lpage></element-citation></ref><ref id="bib78"><element-citation publication-type="book"><person-group person-group-type="author"><name><surname>Wickham</surname><given-names>H</given-names></name></person-group><year>2009</year><source>ggplot2: elegant graphics for data analysis</source><publisher-loc>New York</publisher-loc><publisher-name>Springer</publisher-name></element-citation></ref><ref id="bib79"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yvert</surname><given-names>G</given-names></name><name><surname>Brem</surname><given-names>RB</given-names></name><name><surname>Whittle</surname><given-names>J</given-names></name><name><surname>Akey</surname><given-names>JM</given-names></name><name><surname>Foss</surname><given-names>E</given-names></name><name><surname>Smith</surname><given-names>EN</given-names></name><etal/></person-group><year>2003</year><article-title>Trans-acting regulatory variation in Saccharomyces cerevisiae and the role of transcription factors</article-title><source>Nat Genet</source><volume>35</volume><fpage>57</fpage><lpage>64</lpage><pub-id pub-id-type="doi">10.1038/ng1222</pub-id></element-citation></ref></ref-list></back><sub-article article-type="article-commentary" id="SA1"><front-stub><article-id pub-id-type="doi">10.7554/eLife.00367.027</article-id><title-group><article-title>Decision letter</article-title></title-group><contrib-group content-type="section"><contrib contrib-type="editor"><name><surname>Tautz</surname><given-names>Diethard</given-names></name><role>Reviewing editor</role></contrib><aff><institution>Max Planck Institute for Evolutionary Biology</institution>, <country>Germany</country></aff></contrib-group></front-stub><body><boxed-text><p>eLife posts the editorial decision letter and author response on a selection of the published articles (subject to the approval of the authors). An edited version of the letter sent to the authors after peer review is shown, indicating the substantive concerns or comments; minor concerns are not usually shown. Reviewers have the opportunity to discuss the decision before the letter is sent (see <ext-link ext-link-type="uri" xlink:href="http://www.elifesciences.org/the-journal/review-process">review process</ext-link>). Similarly, the author response typically shows only responses to the major concerns raised by the reviewers.</p></boxed-text><p>Thank you for choosing to send your work entitled “Improved use of a public good selects for the evolution of undifferentiated multicellularity” for consideration at <italic>eLife</italic>. Your article has been favorably evaluated by a Senior editor and 2 reviewers, one of whom is a member of our Board of Reviewing Editors.</p><p>The Reviewing editor and the other reviewer discussed their comments before we reached this decision, and the Reviewing editor has assembled the following comments to help you prepare a revised submission.</p><p>In this paper the authors have set out to test the hypothesis that multicellularity can emerge in otherwise single-celled microbes by selection for improved growth in low-sucrose environments. Sucrose utilization in budding yeast is unusual as it must be hydrolyzed extracellularly by invertase, resulting in glucose and fructose, which can then be imported into the cell. In low sucrose conditions, cells can't grow when they exist as single cells as hydrolyzed glucose and fructose is lost through diffusion. However, the authors have previously shown that if cells are engineered to clump together they can grow in this condition. This is presumably because the local concentration of the products of sucrose hydrolysis (glucose and fructose) is sufficiently high for cells to grow when they clump but not when they are present as individual cells. In this paper the authors have performed a long-term selection in which wild type, non-clumpy cells were propagated by serial dilution in an environment containing low sucrose. They find that this condition does indeed result in selection for clumpy cells. In addition, most of the clones that they analyze have increased expression of the invertase gene (SUC2) and, surprisingly, increased expression of a large number of hexose transporter genes. The authors used a mutator strain and thus mutants with increased fitness have acquired a large number of mutations. To identify causative alleles the authors performed bulk segregant analysis, which allowed them to identify the likely alleles that confer increased fitness. They prove that a fitness increase is attributable to some of these mutations by engineering strains containing only these mutations in addition to reverting acquired mutations in the evolved strains. They speculate that clumping of unicellular organisms may have led to the development of multicellularity.</p><p>This is an excellent dataset that should definitively be published in a major journal. In fact, the experiments, as well as the analysis and the presentation, are mostly very clear. Still, there are some minor comments that need attention in a revised version.</p><p>1) There is little discussion about the general implications of the findings. Hence, the reader is left with the impression that it would have been rather easy to evolve the first steps to multicellularity, but would ask him or herself why it took such a long time until the first multicellular organisms emerged. Maybe the system studied is actually an oversimplification that has only an indirect implication for the evolution of multicellularity. We are of course aware that any such discussion would be rather speculative, but it might nonetheless be worthwhile to address at least some of the pertinent issues.</p><p>2) Regrettably, the authors have neglected to include an essential control: a population of cells propagated in parallel in a non-sucrose environment. Ideally, this would be a low glucose/fructose environment using identical methods. Using long term-selections it is very easy to select for clumping cells either intentionally (Ratcliff et al 2012) or unintentionally (Lang et al., 2011). As the authors set out to demonstrate that it is a benefit of sharing secreted invertase that drives selection for clumpiness, we would have expected this simple control. Alternatively, the same selection using a cytoplasmically expressed invertase and a <italic>MAL11</italic> overexpression allele would also serve as a negative control. This would provide the proof that the observed outcome is unique to the sucrose selection. Nevertheless, we do not consider the lack of this control as a fatal flaw as the authors present multiple post hoc analyses to demonstrate that the fitness increase is unique to a low sucrose environment. The most compelling being the demonstration that the evolved clones have no fitness difference (or slightly reduced fitness) compared with wild type in a low glucose/fructose environment.</p><p>3) The authors use of bulk segregant analysis to identify beneficial mutations and engineering to recreate strains is impressive. However, some simple tetrad dissection would enable the authors to make more definitive statements about the effect of ace2 mutations on clumpiness and Rgt1 pathway mutations' HXT expression (assayed using qPCR). In addition they could test whether the Ubr1 mutations segregate with any of the phenotypes. Segregation analysis of <italic>SUC2</italic> expression in a cross would test whether it is always associated with clumpiness.</p><p>4) The fact that the “recreated” strains are less fit than the evolved clones suggests that something is being missed. One obvious explanation is copy number variation. Studies from the Rosenzweig, Botstein, and Sherlock labs have shown that long-term propagation of yeast strains in nutrient poor environments using chemostats results in selection for transporter amplifications. As the low sucrose environment is also a nutrient poor one, we would expect amplification alleles to be beneficial. However, the authors' data do not allow testing for this as they have performed RNASeq analysis of clones. Were any individual clones sequenced that can be analyzed for copy number variation? The authors could also check in the bulk segregant sequencing data for evidence of copy number variation at the <italic>SUC2</italic> locus or any of the HXT loci that appear to be massively increased in expression.</p><p>5) Many wild <italic>S. cerevisiae</italic> strains are clumpy and loss of clumpiness is likely a lab-adapted trait. The authors should make it clear that loss and gain of clumpiness is likely to have occurred many times in yeast populations and provide other explanations for why clumpiness in yeast cells might be recurrently selected.</p></body></sub-article><sub-article article-type="reply" id="SA2"><front-stub><article-id pub-id-type="doi">10.7554/eLife.00367.028</article-id><title-group><article-title>Author response</article-title></title-group></front-stub><body><p><italic>1) There is little discussion about the general implications of the findings. Hence, the reader is left with the impression that it would have been rather easy to evolve the first steps to multicellularity, but would ask him or herself why it took such a long time until the first multicellular organisms emerged. Maybe the system studied is actually an oversimplification that has only an indirect implication for the evolution of multicellularity. We are of course aware that any such discussion would be rather speculative, but it might nonetheless be worthwhile to address at least some of the pertinent issues</italic>.</p><p>We added the following as the last paragraph of the Discussion section:</p><p>“If undifferentiated multicellularity is so easy to evolve in eukaryotes, then why did it take so long for fossils of differentiated multicellular organisms to appear in the fossil record? Indeed, there is suggestive evidence of undifferentiated eukaryotic multicellularity 2 billion years ago (El Albani et al. 2010), but the radiation of more complex multicellularity is not seen for over a billion years later. This disparity has received a large amount of attention in paleontology, particularly in the context of the seemingly rapid radiation of the animals at the base of the Cambrian. Possible explanations for this “Cambrian explosion” include novel selective forces due to the environmental changes that occurred shortly before the Cambrian or the possibility that evolutionary arms races among a few early species drove the radiation (Marshall 2006; Peters and Gaines 2012). Alternatively, the evolution of differentiated multicellularity may have been slow and gradual and we simply lack the fossil evidence. Whatever the case, it is interesting to ask why the first differentiation might have arisen. From our data, we speculate that the effectiveness of clumping as a strategy could be improved if a limited number of cells increase their enzyme secretion. Such a simple division of labor may have been the precursor to complex multicellularity in the lineages that lead to animals and other complex multicellular species.”</p><p><italic>2) Regrettably, the authors have neglected to include an essential control: a population of cells propagated in parallel in a non-sucrose environment. Ideally, this would be a low glucose/fructose environment using identical methods. Using long term-selections it is very easy to select for clumping cells either intentionally (Ratcliff et al 2012) or unintentionally (Lang et al., 2011). As the authors set out to demonstrate that it is a benefit of sharing secreted invertase that drives selection for clumpiness, we would have expected this simple control. Alternatively, the same selection using a cytoplasmically expressed invertase and a</italic> MAL11 <italic>overexpression allele would also serve as a negative control. This would provide the proof that the observed outcome is unique to the sucrose selection. Nevertheless, we do not consider the lack of this control as a fatal flaw as the authors present multiple post hoc analyses to demonstrate that the fitness increase is unique to a low sucrose environment. The most compelling being the demonstration that the evolved clones have no fitness difference (or slightly reduced fitness) compared with wild type in a low glucose/fructose environment</italic>.</p><p>We are frequently asked this question and the simplest answer would indeed be to be able to state that we had done the experiment suggested by the reviewer and found no multicellularity in strains evolved on low monosaccharide concentrations. There are three strong lines of evidence that argue that the evolution of multicellularity was a specific response to low sucrose concentrations: (1) the failure of the evolved clones to show increased fitness on low monosaccharide concentrations; (2) the fitness advantage of engineered clumpy strains in sucrose and not in low monosaccharide; (3) the failure to select for the mutations that cause multicellularity (<italic>ace2</italic> in Recreated 2 and <italic>mck1</italic>, <italic>irc8</italic>, and <italic>gin8</italic>) on low monosaccharide concentrations when we analyze the effect of individual mutations as they segregate in a cross between the ancestral and evolved strains.</p><p><italic>3) The authors use of bulk segregant analysis to identify beneficial mutations and engineering to recreate strains is impressive. However, some simple tetrad dissection would enable the authors to make more definitive statements about the effect of ace2 mutations on clumpiness and Rgt1 pathway mutations' HXT expression (assayed using qPCR). In addition they could test whether the Ubr1 mutations segregate with any of the phenotypes. Segregation analysis of</italic> SUC2 <italic>expression in a cross would test whether it is always associated with clumpiness</italic>.</p><p>We made 9 allelic combinations of strain Recreated2 and tested them for clumpiness, <italic>HXT4</italic> expression, and <italic>SUC2</italic> expression using qPCR. We concluded that the <italic>ACE2</italic> mutation caused clumpiness and the <italic>MTH1</italic> mutation caused increase in <italic>HXT4</italic> expression (as the reviewers suspected). The <italic>SUC2</italic> expression was more complex and we were not able to make any conclusions from these data. The <italic>UBR1</italic> mutation did not seem to influence any of the phenotypes. These additional experiments are detailed in the Results and the Materials and methods sections. The data from the experiment is in Figure 8–figure supplement 1.</p><p><italic>4) The fact that the “recreated” strains are less fit than the evolved clones suggests that something is being missed. One obvious explanation is copy number variation. Studies from the Rosenzweig, Botstein, and Sherlock labs have shown that long-term propagation of yeast strains in nutrient poor environments using chemostats results in selection for transporter amplifications. As the low sucrose environment is also a nutrient poor one, we would expect amplification alleles to be beneficial. However, the authors' data do not allow testing for this as they have performed RNASeq analysis of clones. Were any individual clones sequenced that</italic> can <italic>be analyzed for copy number variation? The authors could also check in the bulk segregant sequencing data for evidence of copy number variation at the</italic> SUC2 <italic>locus or any of the HXT loci that appear to be massively increased in expression</italic>.</p><p>We had looked for evidence of copy number variation across all genomes and found a noncausal major duplication event, but we failed to explain this as clearly as we should have done. At the reviewers' suggestion, we looked more carefully at read depth, particularly for amplification of the hexose transporters and <italic>SUC2</italic>. We found one case of amplification at these loci: the <italic>HXT6/HXT7</italic> locus was duplicated in EvoClone7C. This amplification was also present in the backcrossed strain, indicating that it may have been causal. We now have the following paragraph in the Results section describing the three total events we found:</p><p>“We analyzed read depth across all evolved clones and saw three major duplication or deletion events: (1) Chromosome 3 in EvoClone10 was duplicated in the region between 152 and 172 kb from the left end of the chromosome. This doubling of read depth was reduced to a 50% increase in the segregated pool, indicating that the duplication was not causal. (2) The number of P-Type ATPases at the <italic>ENO1/2/5</italic> locus was reduced from 3 to 2 in EvoClone6 and EvoClone10. This reduction was also present in the backcrossed strain, indicating that it may have been causal. (3) The number of hexose transporter gene repeats at the <italic>HXT6/7</italic> locus increased from two in the ancestral strain to three in EvoClone7C. This amplification was also present in the backcrossed strain, indicating that it also may have been causal. Amplification of the <italic>HXT6/7</italic> locus has been found in other experimental evolutions (Brown et al. 1998; Gresham et al. 2008; Kao and Sherlock 2008).”</p><p><italic>5) Many wild</italic> S. cerevisiae <italic>strains are clumpy and loss of clumpiness is likely a lab-adapted trait. The authors should make it clear that loss and gain of clumpiness is likely to have occurred many times in yeast populations and provide other explanations for why clumpiness in yeast cells might be recurrently selected</italic>.</p><p>We agree and we have added the following text to the discussion:</p><p>“Although lab strains of yeast were selected to be unicellular during their domestication (Mortimer and Johnston 1986), the level of clumpiness varies in wild isolates of yeast, and undifferentiated multicellularity may have been recurrently selected for its advantage in harvesting nutrients or for protection from stress (Smukalla et al. 2008).”</p></body></sub-article></article> |