The Evolution of Combinatorial Gene Regulation in Fungi
[['The recent genome sequencing and annotation of the major model organisms established that organismal complexity does not scale in a simple way with gene count .', 'This discordance is consistent with earlier proposals that “tinkering” with gene regulation may be a particularly powerful mode of evolution [1–3] .', 'In principle , changes in when and where , and thereby in what combinations , genes are expressed can help to explain changes in organismal complexity over longer time scales .', 'Over shorter time scales , the contributions of changes in gene regulation to phenotypic variation have been clearly demonstrated [4 , 5] .', 'For example , small changes in gene regulation underlie the gain and loss of wing spots in Drosophila species [6] and armor in stickleback fish [7] .', 'The plasticity of gene regulatory networks is of interest because it presumably relates directly to the ability of these networks to generate phenotypic novelty [8] .', 'The potential for rapid turnover ( gains and losses ) of transcription factor binding sites was predicted on theoretical grounds [9–11] and was supported by comparisons of cis regulatory sequence both within and between species [12–15] .', 'Recently , experimental localization of four transcription factors across the mouse and human genomes revealed that binding sites have diverged appreciably between these two species [16] .', 'Analogous experiments performed on two transcription factors from closely related yeasts led to similar conclusions [17] , although in this case , it was not clear how the differences in binding related to gains and losses of cis-acting sequences .', 'The ascomycete lineage , which includes the model yeast S . cerevisiae , serves as a powerful framework for investigating the general impact of regulatory evolution , because several of its members are particularly easy to study experimentally .', 'These include the model yeast S . cerevisiae , the dairy yeast K . lactis , and the human pathogen C . albicans .', 'S . cerevisiae and K . lactis diverged more recently than either did from C . albicans; the divergence of S . cerevisiae and C . albicans is thought to have occurred on the order of 300 million years ago [18] .', 'To date , only a handful of comparative gene regulation studies have been carried out in fungi .', 'These include a few large-scale analyses of changes in gene expression [19] and cis regulatory motifs [20 , 21] as well as some smaller-scale studies [22–24] focusing on sets of co-regulated genes .', 'Whereas the whole-network studies have generally uncovered an abundance of divergence , the smaller-scale studies have characterized this divergence in greater detail or provided mechanistic insight into transcriptional rewiring .', 'Here we take an approach intermediate in scale and attempt to characterize the evolution of a large combinatorial circuit comprised of the MADS-box transcriptional regulator Mcm1 and each of its cofactors .', 'Mcm1 has been intensively studied in S . cerevisiae and , in most cases , it is found as a homodimer that binds DNA cooperatively with other sequence-specific DNA binding cofactors to regulate sets of genes , which we refer to here as regulons .', 'Five regulons have been identified in S . cerevisiae where Mcm1 acts in combination with a second transcriptional regulator .', 'Mcm1 joins with the following: ( 1 ) MATα2 to turn off the a-specific genes ( asgs ) , ( 2 ) MATα1 to turn on the α-specific genes ( αsgs ) , ( 3 ) Fkh2 and Ndd1 to activate G2/M-specific genes , ( 4 ) Yox1 to repress the M/G1-specific genes , and ( 5 ) Arg80 and Arg81 to either repress or activate the arginine metabolic genes [25] .', 'Because Mcm1 itself is not generally regulated , it is typically the regulation of its cofactors that produces the effect of differential gene regulation at each of the Mcm1–cofactor regulons [25] .', 'For example , it is the regulated binding of Mcm1′s cofactor Yox1 that leads to the M/G1-specific expression of genes in the Mcm1-Yox1 regulon [26] .', 'At these Mcm1–cofactor regulons , Mcm1 is thought to increase specificity through added protein–DNA and protein–protein interactions [27] .', 'Previously we showed Mcm1 to be at the center of a rewiring event that led to the replacement of one cofactor ( MATa2 ) with another ( MATα2 ) [24] .', 'In principle , the free energy gain contributed by the interaction between Mcm1 and its flanking cofactor could catalyze evolutionary change by expanding the space of cis-regulatory sequences that yield appropriate gene regulation .', 'For instance , mutations that strengthen Mcm1′s interaction with its cofactor or with DNA can compensate for mutations to the cofactor–DNA interaction , thereby expanding the possibilities for cross-reaction with a new DNA binding protein .', 'This idea bears at least a formal similarity to the neutral networks in RNA sequence space studied by Fontana and colleagues [28] .', 'Because Mcm1 participates in many combinatorial interactions in S . cerevisiae , and because it regulates a large number of genes , we felt that Mcm1 provided a particularly strong entry point to study the evolution of combinatorial networks .', 'To study this problem we performed ChIP-Chip ( chromatin immunoprecipitation , analyzed genome-wide using microarrays ) on Mcm1 in three species ( S . cerevisiae , K . lactis , and C . albicans ) and combined this data with informatics analyses across 32 fungal species .', 'We found that all five Mcm1–cofactor regulons currently characterized in S . cerevisiae are present at least in limited form in K . lactis and C . albicans , suggesting an ancient origin of these regulons .', 'Although the Mcm1–cofactor interaction is typically conserved and a small set of core target genes remains part of the regulon in each species , most of these regulons have undergone substantial divergence through gain and loss of cis-acting sequences .', 'On the global level , substantial gain and loss of Mcm1 binding sites is also evident .', 'Although some of this , as discussed above , is due to target genes moving in and out of existing regulons , much of it is due to the evolution of entirely new Mcm1–cofactor regulons .', 'We highlight two specific instances in which combinatorial regulation by Mcm1 and a cofactor is gained; in one case , we observe the large-scale convergent evolution of regulation at the ribosomal genes and in the other , we describe a very recent gain of regulation that was likely shaped by the selective pressures of the human host .', 'The picture that emerges from this study is one of massive transcriptional rewiring in species that span approximately the same range of protein sequence divergence as human , fish , and sea squirt [29 , 30] .', 'This rewiring consists of both rapid turnover of cis-acting sequence and the formation of new combinations of regulatory proteins .'], ['Mcm1 was chromatin immunoprecipitated ( ChIP-ed ) from S . cerevisiae , K . lactis , and C . albicans a cells using peptide antibodies custom designed for the Mcm1 ortholog of each species .', 'To maximize the detection of Mcm1 binding , each strain was grown under two different conditions known to stimulate binding of Mcm1: yeast extract peptone dextrose ( YEPD ) medium and pheromone-inducing medium with α pheromone ( details in the Text S1 ) .', 'Immunoprecipitate and whole-cell extract samples were competitively hybridized to custom-designed Agilent microarrays that tile 60mer probes at a median spacing of 66 , 59 , and 79 base pairs ( bp ) across the genomes of S . cerevisiae , K . lactis , and C . albicans , respectively ( Figure S1 ) .', 'For each species/condition , the ChIP-Chip was performed twice and the two biological replicates were combined in downstream data processing .', 'Data were processed by a variety of methods , and it was determined empirically that the Joint Binding Deconvolution ( JBD ) algorithm [31] provides the best combination of consistency across species and accuracy on a test set of previously characterized S . cerevisiae binding sites ( see Text S1 ) .', 'Complete ChIP profiles for all experiments can be viewed at: http://genome . ucsf . edu/mcm1_evolution/ .', 'The majority of regions that JBD called as bound by Mcm1 contained at least one instance of the well-characterized Mcm1 binding motif [32 , 33] .', 'We therefore decided to incorporate motif information into our final criteria for Mcm1-bound segments .', 'De novo motif finding by MEME [34] on a set of high-confidence bound regions from JBD yields Mcm1 binding site motifs that are roughly the same in each species ( Figure 1 ) ; the motif deduced de novo from S . cerevisiae closely resembles previously described Mcm1 recognition sequences .', 'In C . albicans , there was a large subset of bound regions without a canonical Mcm1 motif .', 'These regions are largely explained by the appearance of a noncanonical motif ( Figure 1 ) , discussed later .', 'Parameter cutoffs for JBD statistics and the motif p-value were chosen that correctly call 85% ( 28 of 33 genes ) of our S . cerevisiae test set as positives while also calling an additional 219 of 5 , 769 genes as bound .', 'Details regarding test set selection are provided in the Text S1 along with a discussion of false-positive rates and receiver operator characteristic plots ( Figure S8 ) evaluating a variety of parameter value choices .', 'These same cutoffs used for the S . cerevisiae data yield 626 of 5 , 327 genes bound in K . lactis and 761 of 6 , 090 genes bound in C . albicans ( gene lists in Table S1 ) .', 'For these and all subsequent calculations , Mcm1 targets from the two growth conditions examined have been pooled .', 'After defining Mcm1 targets in each species , we sought to evaluate the overlap of these targets between species .', 'We mapped orthologs using an existing algorithm [24] , which was modified to reduce directional bias ( Text S1 , section titled: “Mapping orthologous gene sets” ) , on an updated database of open reading frame ( ORF ) sequences from 32 fully sequenced genomes ( Table S2 ) .', 'Genes bound by Mcm1 in each species A were then “mapped to” one of the other two species B via our ortholog map .', 'The number of genes “mapped from” A and also found to be in the Mcm1 bound gene set of B was counted and is displayed as a fraction of the total genes bound in species A that can be mapped to species B ( Figure 2A ) .', 'Note the lack of symmetry; comparing the Mcm1 bound gene set of species A to that of B is not the same as comparing the bound gene set of species B to that of A because of the different total number of Mcm1-bound genes in the different species .', 'Overlap p-values were also calculated for each species pair by using the hypergeometric distribution ( Figure 2B ) .', 'There is significant overlap in Mcm1-targeted genes between each pair of species ( p < 10−3 ) .', 'As might be expected , conservation is strongest between the two more closely related species , S . cerevisiae and K . lactis , with 42% of mapped S . cerevisiae Mcm1 targets also bound by K . lactis Mcm1 .', 'However , as the lower frequency ( 16% ) of K . lactis Mcm1-mapped targets bound by S . cerevisiae Mcm1 indicates , the K . lactis Mcm1 target set is much larger .', 'Interestingly , the C . albicans Mcm1 target set overlaps more significantly with the K . lactis Mcm1 set than it does with the S . cerevisiae Mcm1 set , indicating that a sizeable fraction of the extra genes bound by K . lactis Mcm1 are shared with C . albicans Mcm1 and are therefore likely to have been lost as Mcm1 targets on the branch leading to S . cerevisiae ( see next section ) .', 'For simplicity , we have focused here on only those genes that can be mapped in a 1:1 fashion between species .', 'However , similar results are obtained when genes with more complex interspecies mappings ( e . g . , 2:1 ) are included .', 'To rule out the possibility that our results were biased by the exact parameters chosen , we repeated the analysis with a variety of parameter choices and obtained similar results ( Figure S9 ) .', 'To assess the prevalence of gain and loss of Mcm1 binding sites across the three-species phylogeny , we constructed a model with nine parameters: four gain rates and four loss rates , corresponding to each of the four branches of the rooted tree , and a single parameter representing the probability of an Mcm1 binding site at the root of the tree ( Figure 2C ) .', 'We take as our dataset the Mcm1-binding occurrence patterns at each of the 2 , 766 genes that can be mapped between S . cerevisiae , K . lactis , and C . albicans in a 1:1:1 fashion via our ortholog mapping .', 'There are eight such patterns , e . g . , the pattern “101” for hypothetical gene X indicates an Mcm1 binding site is present upstream of gene X in S . cerevisiae and C . albicans , but not in K . lactis .', 'We devised a modified maximum-likelihood algorithm to infer the gain and loss rates on each branch of the three-species phylogeny .', 'A more thorough description of this procedure is given in the Text S1 .', 'The results show a high degree of binding site turnover on all branches of the tree .', 'For example , we estimate that the last common ancestor of S . cerevisiae and K . lactis had Mcm1 binding sites at 156 genes .', 'Since divergence , Mcm1 binding sites were gained at 44 genes and lost at 109 genes in the S . cerevisiae lineage .', 'Likewise , Mcm1 binding sites were gained at 128 genes and lost at 38 genes in the K . lactis lineage .', 'Thus , present day S . cerevisiae and K . lactis have only 38 Mcm1-targeted genes in common .', 'We do not believe that this analysis is biased by any systematic failures to detect Mcm1 binding sites in our ChIP experiments—either through experimental biases or because the growth conditions chosen did not promote Mcm1 binding .', 'In cases where Mcm1 is bound upstream of a gene in one species but not in the other two species , the Mcm1 motif is generally not present in those other two species as well ( Text S1 , section titled: “Mcm1 DNA motifs are not present at genes that are not bound” ) .', 'In the sections that follow , we will further dissect the changes in combinatorial regulation that give rise to the conservation and divergence summarized in Figure 2C .', 'If we consider just the subset of genes that has a 1:1:1 mapping in our ortholog table , only 12 genes ( ∼13% of the genes bound in S . cerevisiae ) are part of the Mcm1 circuit in all three species ( Figure 3 ) .', 'If gene duplications are allowed , the number of genes in S . cerevisiae with at least one “ortholog” bound in K . lactis and C . albicans is 45 ( ∼18% of the genes bound in S . cerevisiae ) .', 'Presumably this conserved set of target genes reflects a conserved role played by Mcm1 in the common ancestor as well as in the three modern species .', 'The set of ancestral Mcm1-bound genes is clearly enriched for genes regulated by the cell cycle ( Figure 3 , shaded orange ) and mating type ( Figure 3 , shaded blue ) .', 'The latter is confirmation of results from our previous study [24] describing the conservation of membership within the a-specific gene regulon despite the dramatic switch from positive regulation by MATa2 to negative regulation by MATα2 .', 'In S . cerevisiae , the cell cycle genes listed are regulated by the Mcm1 cofactors Fkh2/Ndd1 and Yox1 .', 'The conservation of these genes as targets of Mcm1 prompted us to inquire whether combinatorial control by Mcm1 and each of its S . cerevisiae cofactors was also conserved since the time when S . cerevisiae , K . lactis , and C . albicans diverged from a common ancestor .', 'The Mcm1–cofactor regulons of S . cerevisiae were mapped to Mcm1-bound regions in K . lactis or C . albicans , and motif finding was performed to identify cis-regulatory elements controlling the orthologous regulons ( details presented in Text S1 ) .', 'The results of this analysis ( Figure 4A and 4B ) demonstrate that most known Mcm1–cofactor interactions from S . cerevisiae are present in K . lactis and C . albicans and are therefore likely of ancient origin .', 'In the description that follows , we first compare the cis-regulatory motifs of the more closely related S . cerevisiae and K . lactis and then compare these to the motifs of the more divergent C . albicans .', 'Here we use the term “interaction” to refer to both demonstrated protein–protein interactions as well as those inferred from the co-occurrence in cis of two or more regulatory motifs .', 'One caveat of this approach is that co-occurrence of motifs can arise from cooperative as well as competitive binding of two transcription factors .', 'However , we think the latter is unlikely for most cases documented in this work , because the spacing of the motifs tends to be highly constrained and nonoverlapping , a feature typically observed for cooperative binding with Mcm1 .', 'In general the cis-regulatory elements of the K . lactis and S . cerevisiae Mcm1–cofactor regulons are similar , suggesting that the corresponding Mcm1–cofactor interactions have changed little since these two lineages split .', 'Notable exceptions are the changes seen at asgs discussed previously [24] and the apparent added Fkh2 specificity flanking several of the Yox1-Mcm1 sites in K . lactis ( Figure 4B ) .', 'Although the latter is seen in at least a few genes in S . cerevisiae [26] , this Yox1-Mcm1-Fkh2 architecture appears much more prominent in K . lactis .', 'Comparison of K . lactis and S . cerevisiae to the more divergent C . albicans reveals that a number of changes to cis-regulatory motifs have occurred over longer time scales .', 'At the Fkh2-Mcm1 regulon , there is a shift in the placement of the Fkh2 site relative to the Mcm1 site by 1 base pair , which occurs across the entire regulon .', 'We note that species with the tighter Fkh2-Mcm1 spacing have clear orthologs to Ndd1 , a protein which in S . cerevisiae binds the Fkh2-Mcm1 complex periodically , thereby driving G2/M-specific expression [35] , while those with the lengthened spacing do not .', 'It is not known how the Fkh2-Mcm1 complex of C . albicans would function to drive G2/M-specific gene expression without an Ndd1 ortholog , although this altered spacing may provide a clue .', 'It is also noteworthy that Fkh2 is related to another protein , Fkh1 , which is derived from the yeast whole-genome duplication event [36] , meaning that these two genes found in S . cerevisiae map to a single gene in K . lactis and C . albicans .', 'It is known that Fkh2 binds DNA cooperatively with Mcm1 , but that Fkh1 does not [37] .', 'Given the evidence for the Fkh2-Mcm1 motif in K . lactis and C . albicans , we infer that this interaction is ancestral to the species under study and that after duplication , only Fkh2 retained the ability to bind cooperatively with Mcm1 .', 'The cis-regulatory motif at the MATα1-Mcm1 regulon has clearly changed as well , indicating that MATα1 , despite its obvious conservation , recognizes distinct DNA motifs in different species .', 'However , the altered MATα1 motif observed in C . albicans is not necessarily incompatible with the S . cerevisiae protein , a surmise based on previous mutagenesis studies [38] .', 'Despite this change in motif , experimental evidence indicates that MATα1′s function as an activator of αsgs is the same in S . cerevisiae and C . albicans [39] .', 'Once it was determined that most Mcm1–cofactor pairings are conserved across the species we examined , we then determined to what extent the set of genes in their corresponding regulons was also conserved .', 'The motif matrices for each of the Mcm1–cofactor pairs were used to score the entire set of Mcm1-bound sequences in each species and thus to define the members of each Mcm1–cofactor regulon in each species ( Text S1 ) .', 'We found that the number of targets in each regulon is roughly the same across the three species , but the precise set of members is not .', 'However , within each regulon , there is a small , core set of conserved genes ( Figure 4C ) .', 'For example , the Fkh2-Mcm1 regulon consists of roughly 20 genes in each species , but only three genes are part of the regulon in all three species .', 'Previously we showed that for the asg regulon at least , this core is conserved throughout the yeasts spanning the lineage of S . cerevisiae and C . albicans [24] .', 'A similar promoter sequence analysis with the Mcm1-Fkh2 matrices supports a conserved core within this regulon as well ( unpublished data ) .', 'For example , the promoters of BUD4 and CDC20 have strong matches to the Fkh2-Mcm1 matrix in most species within the lineage spanning S . cerevisiae and C . albicans .', 'Thus , turnover within these regulons is not a purely stochastic process , but rather is constrained in some respects by purifying selection .', 'As summarized in Figure 2C , this study revealed many specific instances of gains and losses of Mcm1 regulation across the ascomycete lineage .', 'The large number of changes seen at the global level , however , can not be fully accounted for by binding site turnover within the ancestral Mcm1–cofactor regulons alone ( Figure 4C ) .', 'In the following section , we highlight three examples of large-scale rewiring events , chosen for their particular clarity and their relevance to well-developed systems .'], ['In this work , we have tracked the evolution of combinatorial gene regulation by the highly conserved transcriptional regulator Mcm1 and each of its known cofactors across the ascomycete fungal lineage .', 'Our analysis shows that the genes regulated by Mcm1 have changed considerably over the evolutionary time scales represented by this lineage; our results reveal many more differences than similarities in the Mcm1 circuitry .', 'Regulation by Mcm1 is more pervasive in K . lactis and C . albicans , where 12% of all genes are bound , than in S . cerevisiae , where 4% of genes are bound .', 'The fraction of genes shared as targets between all three species is very low ( 13%–18% ) , and we have demonstrated that this is due to both substantial gain and loss of Mcm1 binding sites along each branch of this phylogeny ( Figure 2B ) .', 'The extensive amount of gain and loss observed is consistent with recent studies in mammals [16] and closely related yeasts [17] and suggests the following three possibilities: ( 1 ) there is a richness of selective advantages offered in the dynamic rewiring of gene regulatory networks , ( 2 ) there are a large number of neutral alternatives to gene regulation by Mcm1 , or ( 3 ) selection on gene expression is weak .', 'The latter possibility seems at odds with other observations such as the large fraction of genes devoted to transcriptional regulation in S . cerevisiae ( ∼3% ) , the greater-than-expected number of transcriptional regulators retained after the whole genome duplication ( ∼6% versus ∼3% ) , and the considerable conservation found in many S . cerevisiae promoters [54 , 55] .', 'Additionally , the fact that many of the Mcm1 sites are enriched at functionally related genes and often found in tandem with cofactor motifs argues strongly against the hypothesis that a large number of these sites are fortuitous and nonfunctional .', 'Gauging the relative contributions of selection versus neutral drift on the gene regulatory networks will be an exciting challenge for future research [56] .', 'Despite the highly dynamic nature of evolution of Mcm1 regulation , we find evidence that most Mcm1–cofactor interactions characterized in S . cerevisiae are also present in K . lactis and C . albicans ( Figure 4B ) .', 'Although the Mcm1–cofactor pairings are conserved , the set of genes that each regulates has diverged considerably across species .', 'Nonetheless , each Mcm1–cofactor pair targets a small core of genes conserved as part of the regulon .', 'These regulon cores are enriched for genes functioning in the cell cycle and mating .', 'Thus it would seem that Mcm1′s role in these processes evolved prior to the split of the species we have chosen to study .', 'Nevertheless , even at these conserved regulons , there are many species-specific differences .', 'For example , across an entire regulon , the spacing between Fkh2 and Mcm1 binding sites has changed in S . cerevisiae and K . lactis relative to C . albicans , as have the DNA recognition sequences of MATα1 .', 'This latter observation is particularly interesting because it suggests that the specificity of MATα1 has evolved without an accompanying gene duplication .', 'In addition to the conservation of Mcm1–cofactor interactions associated with cell cycle and mating , we see the evolution of new Mcm1–cofactor regulons .', 'For example , Mcm1 binding sites are gained at the majority of ribosomal genes in K . lactis in close proximity to binding motifs for another transcription factor , Rap1 ( Figure 5C and 5D ) .', 'The evolution of ribosomal gene regulation has been studied previously [21] , but a role for Mcm1 was not discussed .', 'Our new results support the idea , first proposed by Tanay et al . [21] , that while the protein sequence of this critical macromolecular machine has remained nearly constant , its regulation has undergone substantial diversification in yeasts .', 'What is perhaps most surprising is our finding that the set of species that contain Mcm1 binding motifs upstream of ribosomal genes ( Figure 5A and 5B; C . glabrata , K . lactis , Y . lipolytica , and the A . nidulans lineage ) do not cluster phylogenetically .', 'From this we inferred that Mcm1 binding at ribosomal genes likely evolved on four separate occasions .', 'If further genome sequencing continues to support this result , this will serve as the largest instance of convergent regulatory evolution yet reported .', 'The relatively sudden appearance of Mcm1 binding sites in close proximity to Rap1 sites at roughly 70 ribosomal genes in K . lactis raises another important question: Can the commonly accepted mutational processes , such as point mutation and recombination , support this scale of concomitant changes—or must some alternative mechanism for moving promoters around the genome be invoked [57 , 58] ?', 'One can argue that , without a redundant mechanism in place , loss or gain of Mcm1 regulation of even a single gene means losing precise control over one component of a macromolecular complex that is thought to need tight stochiometric control [40] .', 'With further sequencing and characterization of Mcm1′s functional role at the ribosomal genes , it may become clear how such a massive regulatory change can take place at a set of genes encoding such highly conserved , tightly regulated and essential proteins .', 'In C . albicans , we identified the presence of Mcm1 at a noncanonical motif upstream of roughly 110 genes .', 'The noncanonical motif differs significantly from the canonical Mcm1 motif ( Figure 1 ) , although in both cases GC-rich regions flank an AT-rich core .', 'To our knowledge no MADS-box domain has ever been shown to bind a sequence this far diverged from the canonical Mcm1 motif .', 'Even so , we find that noncanonical motifs tend to be centered with respect to peaks of ChIP-Chip enrichment and thus conclude that Mcm1 either binds this motif directly with some unknown cofactor or some unknown transcriptional regulator recognizes this motif and interacts strongly with Mcm1 .', 'The set of genes at which Mcm1 binds the noncanonical motif is enriched for processes such as adhesion and contains three of four known regulators of the white-opaque phenotypic switch [45] .', 'The white-opaque switch is of considerable interest because the white and opaque states are heritable and because the two states are thought to allow adaptation to different niches within a human host [49 , 50] .', 'In this vein , the evolution of regulation associated with the switch deserves special attention too , because the changes seen here represent , to our knowledge , the first gene regulatory changes to be associated with a heritable biological process and one of only a few instances implicated to play an adaptive role in fungal biology [59] .', 'The results of our comparative analysis of 32 yeast species demonstrate that Mcm1 binding at the noncanonical motif is found only in two very closely related species , C . albicans and C . dubliniensis , and thus likely arose only very recently ( Figure 7D ) .', 'Moreover white-opaque switching has been described only in these two species [60] , which are both pathogens of humans .', 'Thus , the evidence so far suggests that the white-opaque switch may have arisen just prior to the divergence of C . albicans and C . dubliniensis and that the emergence of the noncanonical Mcm1 motif at white-opaque regulators was crucial to this development .', 'Alternatively , the white-opaque switch may have arisen earlier , and the addition of Mcm1 regulation may have refined it in some way , affecting heritability , for example .', 'The picture that emerges from this study is one of massive transcriptional rewiring in species that span approximately the same range of divergence as human , fish , and sea squirt [29 , 30] .', 'Mcm1 regulates hundreds of genes in S . cerevisiae , K . lactis , and C . albicans , but less than 20% of Mcm1–target gene connections are preserved in all three species .', 'The differences arise from target genes moving in and out of ancient Mcm1–cofactor regulons , but also from the formation of new Mcm1–cofactor interactions and the loss of ancient ones .', 'Taken together with our previous work [24] , we have now provided evidence for the gain of three interactions: Mcm1 with MATα2 , Mcm1 with Rap1 , and Mcm1 with Wor1 .', 'We have also described loss of an interaction between Mcm1 and MATa2 and the loss of an interaction between Mcm1 and Arg81 that was preserved in an Mcm1 duplicate .', 'In attempting to judge the relative contributions of combinatorial control per se to the evolution of transcriptional circuits , we acknowledge that the ideal “control” datasets do not exist .', 'For example , data collected from a large noncombinatorial circuit ( should one even exist ) over several species would allow an objective assessment of the special contribution of combinatorial control to circuit evolution .', 'Nonetheless , our results provide experimental and informatic support for the idea that combinatorial networks are highly evolvable [61–64] , and perhaps more importantly , they document specific mechanisms by which one large combinatorial circuit has evolved .'], ['Detailed methods can be found in Text S1 , Figure S1–S10 , and Table S1–S3 .', 'The information can also be found in one complete file , Protocol S1 .']]
['Introduction', 'Results', 'Discussion', 'Methods']
['It is widely suspected that gene regulatory networks are highly plastic .', 'The rapid turnover of transcription factor binding sites has been predicted on theoretical grounds and has been experimentally demonstrated in closely related species .', 'We combined experimental approaches with comparative genomics to focus on the role of combinatorial control in the evolution of a large transcriptional circuit in the fungal lineage .', 'Our study centers on Mcm1 , a transcriptional regulator that , in combination with five cofactors , binds roughly 4% of the genes in Saccharomyces cerevisiae and regulates processes ranging from the cell-cycle to mating .', 'In Kluyveromyces lactis and Candida albicans , two other hemiascomycetes , we find that the Mcm1 combinatorial circuits are substantially different .', 'This massive rewiring of the Mcm1 circuitry has involved both substantial gain and loss of targets in ancient combinatorial circuits as well as the formation of new combinatorial interactions .', 'We have dissected the gains and losses on the global level into subsets of functionally and temporally related changes .', 'One particularly dramatic change is the acquisition of Mcm1 binding sites in close proximity to Rap1 binding sites at 70 ribosomal protein genes in the K . lactis lineage .', 'Another intriguing and very recent gain occurs in the C . albicans lineage , where Mcm1 is found to bind in combination with the regulator Wor1 at many genes that function in processes associated with adaptation to the human host , including the white-opaque epigenetic switch .', 'The large turnover of Mcm1 binding sites and the evolution of new Mcm1–cofactor interactions illuminate in sharp detail the rapid evolution of combinatorial transcription networks .']
['In explaining the diversity of organisms on Earth , it is increasingly evident that evolutionary changes in when and where genes are expressed provide a crucial source of variation .', 'By using genome-wide transcription factor localization experiments in S . cerevisiae , K . lactis , and C . albicans , combined with comparative genomics across many more yeast species , we examined how a large combinatorial transcription circuit evolves over the course of hundreds of millions of years .', 'Combinatorial regulation is pervasive in eukaryotic organisms and is thought to allow for increased specificity and integration of multiple signals in the control of gene expression .', 'Our studies focused on one prolific combinatorial regulator , Mcm1 , which , in combination with five cofactors , binds and regulates genes functioning in a diverse range of cellular processes in S . cerevisiae .', 'We found evidence of massive network rewiring , including high rates of gain and loss of Mcm1 binding sites and the formation of new Mcm1–cofactor combinations and the breaking of old ones .', 'We propose that the multiple protein–protein and protein–DNA interactions that specify transcription in combinatorial circuits allow for a richness of compensatory mutations and thereby provide ample opportunity for both adaptive and neutral evolution .']
['biochemistry', 'computational biology', 'evolutionary biology', 'genetics and genomics']
Experimental approaches combined with comparative genomics show how a large combinatorial transcription circuit has rewired over evolutionary time scales.