(C) PLOS One [1]. This unaltered content originally appeared in journals.plosone.org.
Licensed under Creative Commons Attribution (CC BY) license.
url:https://journals.plos.org/plosone/s/licenses-and-copyright

------------



Multiple molecular events underlie stochastic switching between 2 heritable cell states in fungi

['Naomi Ziv', 'Department Of Microbiology', 'Immunology', 'University Of California', 'San Francisco', 'California', 'United States Of America', 'Lucas R. Brenes', 'Alexander Johnson']

Date: 2022-05

Abstract Eukaryotic transcriptional networks are often large and contain several levels of feedback regulation. Many of these networks have the ability to generate and maintain several distinct transcriptional states across multiple cell divisions and to switch between them. In certain instances, switching between cell states is stochastic, occurring in a small subset of cells of an isogenic population in a seemingly homogenous environment. Given the scarcity and unpredictability of switching in these cases, investigating the determining molecular events is challenging. White-opaque switching in the fungal species Candida albicans is an example of stably inherited cell states that are determined by a complex transcriptional network and can serve as an experimentally accessible model system to study characteristics important for stochastic cell fate switching in eukaryotes. In standard lab media, genetically identical cells maintain their cellular identity (either “white” or “opaque”) through thousands of cell divisions, and switching between the states is rare and stochastic. By isolating populations of white or opaque cells, previous studies have elucidated the many differences between the 2 stable cell states and identified a set of transcriptional regulators needed for cell type switching and maintenance of the 2 cell types. Yet, little is known about the molecular events that determine the rare, stochastic switching events that occur in single cells. We use microfluidics combined with fluorescent reporters to directly observe rare switching events between the white and opaque states. We investigate the stochastic nature of switching by beginning with white cells and monitoring the activation of Wor1, a master regulator and marker for the opaque state, in single cells and throughout cell pedigrees. Our results indicate that switching requires 2 stochastic steps; first an event occurs that predisposes a lineage of cells to switch. In the second step, some, but not all, of those predisposed cells rapidly express high levels of Wor1 and commit to the opaque state. To further understand the rapid rise in Wor1, we used a synthetic inducible system in Saccharomyces cerevisiae into which a controllable C. albicans Wor1 and a reporter for its transcriptional control region have been introduced. We document that Wor1 positive autoregulation is highly cooperative (Hill coefficient > 3), leading to rapid activation and producing an “all or none” rather than a graded response. Taken together, our results suggest that reaching a threshold level of a master regulator is sufficient to drive cell type switching in single cells and that an earlier molecular event increases the probability of reaching that threshold in certain small lineages of cells. Quantitative molecular analysis of the white-opaque circuit can serve as a model for the general understanding of complex circuits.

Citation: Ziv N, Brenes LR, Johnson A (2022) Multiple molecular events underlie stochastic switching between 2 heritable cell states in fungi. PLoS Biol 20(5): e3001657. https://doi.org/10.1371/journal.pbio.3001657 Academic Editor: Joseph Heitman, Duke University Medical Center, UNITED STATES Received: December 7, 2021; Accepted: May 4, 2022; Published: May 20, 2022 Copyright: © 2022 Ziv et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Data Availability: All data and scripts used for analysis are available through the Open Science Framework, DOI: 10.17605/OSF.IO/6E9VZ (https://osf.io/6e9vz/). Funding: This study was supported by NIH grant R01AI049187 (A.D.J.) and aided by a grant from The Jane Coffin Childs Memorial Fund for Medical Research (to N.Z.).The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing interests: The authors have declared that no competing interests exist

Introduction Transcription circuits, defined as transcription regulators and the DNA cis-regulatory sequences they bind, control the expression of genes and thereby define cellular identity. Cellular identity is not static; gene expression programs change in response to external and internal stimuli. Certain transcriptional networks can produce distinct cell states that can be maintained through many cycles of cell division. The best characterized of such circuits come from microbes, such as the cI-Cro circuit of phage lambda [1], the control of lactose utilization in Escherichia coli [2], competence in Bacillus subtilis [3], and galactose utilization in Saccharomyces cerevisiae [4]. Such “bistability” has also been engineered in both prokaryotic and eukaryotic synthetic circuits [5–7] and is thought to underlie many instances of cell differentiation in multicellular organisms [8–10]. The production of 2 stable transcriptional states from a single genome depends on positive feedback regulation and nonlinearity (such as cooperativity) within the feedback circuit, which can convert graded inputs into discontinuous switch-like outputs [11–13]. The output of a given circuit “diagram,” however, is entirely dependent on the physical parameters of the components, such as dissociation constants and protein concentrations; even slight variations in circuit architecture or parameters can produce distinct outputs [14]. Many questions remain concerning the functional properties and parameters of network motifs found in real biological systems, particularly those based on the large regulatory networks commonly found in eukaryotic organisms. A major challenge to experimentally investigating the functional role of the multiple feedback loops present in complex circuits is the limited ability to independently manipulate the different components. In order to understand the output of these large networks, it is necessary to integrate the analysis of individual motifs with an analysis of the network as a whole. In some cases, switching between states depends on external signals, but in other cases switching between stable states appears stochastic, occurring in a subset of cells of an isogenic population in a seemingly homogenous environment. In microbes, stochastic phenotypic switching is often thought of in terms of adaptation to fluctuating environments [15,16], where there is a predicted optimal relationship between the frequency of stochastic switching and the frequency of environmental change. In pathogens, stochastic switching can create distinct subpopulations with different cell features, protecting against targeted host defense systems [17]. Stochastic switching also occurs during development to create cellular diversity [18]. For example, photoreceptor patterning in the fly retina is based on a choice between cell states following the stochastic expression of a single transcription factor [19]. Understanding stochastic phenotype switching at a molecular level is challenging, in particular, identifying the early, seemingly random, initiating events that determine which cells will undergo switching. To study the mechanism of stochastic switching between 2 stable transcription states, we investigate white-opaque switching in Candida albicans, a common component of the human gut microbiota but also an opportunistic pathogen causing life-threatening bloodstream infections [20]. White and opaque cells differ in the expression of hundreds of genes resulting in drastic differences in cell morphology, metabolism, the ability to mate, and interactions with the immune system [21–25]. A complex transcriptional network controls white-opaque switching, as many (at least 7) transcription factors are involved and are known to regulate one another, in a series of nested feedback loops [26–28]. The master regulator of white-opaque switching is Wor1, a transcription factor whose deletion completely blocks switching (locking cells in the white state) and whose ectopic expression in white cells converts them en masse to opaque cells, even in the absence of other critical regulators [26,29–31]. Wor1 is thought to positively regulate itself, as ectopic Wor1 expression activates expression from the endogenous Wor1 locus, which is necessary for maintaining the opaque cell state when ectopic expression is removed [31]. Wor1 is differently regulated between the 2 cell states, with 40-fold higher expression in opaque cells. In common with the other regulators of white-opaque switching, the Wor1 regulatory sequences (7 kb upstream control region and 2 kb 5′ UTR) are some of the longest in the C. albicans genome [23]. Its complexity is reminiscent of enhancers from multicellular organisms, even exhibiting the phase transitions described for some mammalian enhancers [32]. Wor1 is positioned in the center of a network of multiple positive feedback loops, binding both its own promoter as well as those of both activators and repressors of the opaque cell state [26,27]. In this work, we investigate the process of white-opaque switching at a single-cell level. To do so, we developed 2 complementary quantitative approaches to study the molecular mechanisms underlying switching. The first is the use of microfluidics and fluorescent microscopy to follow switching in C. albicans, allowing us to quantify the activation of the master regulator Wor1 in pedigrees where both switching and nonswitching cells are dividing. The second is the use of a synthetic inducible system in S. cerevisiae to investigate the function of white-opaque regulators and characterize specific regulatory interactions, such as Wor1 autoregulation. Our results reveal features of white-opaque switching that were unanticipated from previous “bulk” culture analyses and provide key insights into the mechanism behind stochastic switching.

Discussion The ability of cells to switch between different cellular states, each of which are faithfully inherited through many rounds of cell division, underlie a wide range of important phenomena including cell differentiation [40], antibiotic resistance [41,42], and cancer progression [43]. The fungal pathogen C. albicans and its close relatives Candida dubliniensis and Candida tropicalis all undergo a phenomenon known as white-opaque switching, whereby 2 distinct states are stable across many cell divisions, and switching between them is stochastic [44,45]. Although the presence of these heritable states has been preserved by selection in these species, the mechanism of switching and the complete range of physiological roles of the 2 cell types are topics of active investigation. It is known that opaque cells are the mating competent cell type [21] and are less susceptible to phagocytosis by macrophages, potentially promoting immune evasion [24,46]. Our recent work has also shown that, due to a specialized secreted protease program, a small number of opaque cells can promote growth of a large white cell population; thus, the 2 cell types can act synergistically [47]. Aside from its physiological roles, white-opaque switching has many advantages as a eukaryotic system to study stochastic phenotypic switching. It exemplifies 2 crucial phenomena, rare stochastic phenotypic switching and cell states that are stable through thousands of cell divisions. It shares many features with similar networks in multicellular eukaryotes, such as a large number of transcription factors and long, complex regulatory regions (enhancers). White-opaque regulators have also been shown to form phase-separated condensates [32], which are a feature of some mammalian enhancers [48]. A particular advantage of studying phenotypic switching and heritability in C. albicans is the relative ease of performing genetic manipulations. By combining genetic modifications with direct observation of single cells, our goal is to describe in quantitative terms the molecular mechanisms that underlie cell type switching. Positive feedback, encompassing either positive interactions, double negative interactions, or autocatalysis is necessary, but not sufficient for the creation of heritable cell states [11]. To facilitate switching, the response must incorporate nonlinearity, also known as ultrasensitivity, in order to create sharp transitions from gradual changes [49]. Mechanisms for creating or increasing ultrasensitivity in regulatory responses include cooperative binding of transcriptional regulators [50–52], protein sequestration of regulators [53], and the number of steps in transcriptional cascades [54]. We show here that Wor1 regulation of its own promoter is ultrasensitive even in the absence of other C. albicans regulators. While our measurements provide a baseline for the autoregulation of Wor1, other factors in C. albicans may change these parameters. While difficult to rigorously test in the native C. albicans system, we propose this ultrasensitive response is a crucial feature of white-opaque switching that can be incorporated into mathematical models describing switching [55]. Furthermore, the propensity to dimerize added by the GFP causes a small increase in the steepness of the activation curve but a large reduction in the activation threshold. It is this large reduction in activation threshold that likely leads to the large increase in white-opaque switching in C. albicans cells that harbor Wor1-GFP. This observation suggests switching depends on the equilibrium between Wor1 monomers and higher oligomers. While the natural threshold results in rare switching under standard lab conditions, our results show that switching frequency is easily adjusted by a relatively small change in the propensity of Wor1 to dimerize. Indeed, other white-opaque regulators may function in a similar manner to the GFP fusion; any increase in Wor1 multimer formation would produce a large effect on switching and could amplify a small amount of Wor1 expression. Although many questions remain concerning the molecular details of Wor1 autoregulation, our results quantify, for the first time, the ultrasensitivity of the WOR1 control region (enhancer) in response to carefully measured levels of the Wor1 protein. A particular advantage of our study is the use of microfluidics, which has allowed us to observe white-opaque switching at a single-cell level and at high temporal resolution. The ability to track single cells and their lineages through time allowed us to make an unexpected observation that switching can occur in multiple groups of related cells, even when these groups are separated by cell divisions that produce cells that do not switch. The presence of this pattern did not differ between strains containing Wor1-GFP versus Wor1-mGFP despite different switching frequencies. We can draw 2 conclusions from this switching pattern. First, at least 2 stochastic events are needed for cells to switch from white to opaque. The first event defines small pedigrees of cells with increased probability of switching and the second event, which only occurs in a subset of cells in these pedigrees, directly leads to Wor1 activation and cell switching. Second, when a complete switching event occurs, several closely related cells in the pedigree switch simultaneously, a property that may be adaptive. The occurrence of coordinated switching would be advantageous in situations where the rapid generation of the rare cell type could be beneficial. Although our experiments do not address the adaptive value of white-opaque switching, increasing the number of opaque cells could be valuable both for immune evasion or metabolic cooperation with white cells. A similar switching pattern was observed in a lab-modified galactose utilization pathway in S. cerevisiae [56]. In that study, removal of negative feedback created a system dominated by a single positive feedback loop and led to cells stochastically turning on the GAL regulon in the absence of galactose. In a very similar manner to our observations, the authors saw correlated activation in cell pedigrees, at times separated by cells not activating GAL expression [56]. This comparison suggests these coordinated switching patterns might be a general feature in eukaryotic systems. One hypothesis to account for our observed correlations among related cells relies on the stochastic and bursty nature of gene expression [57]. According to this hypothesis, the first stochastic event during white-opaque switching might be a rare burst of low-level Wor1 transcription. This state would be heritable over only a small number of cell divisions, defining the small pedigrees of cells predisposed to switch. The second stochastic event would then create variation in the amount of Wor1 between cells, perhaps because of variation in translation efficiency, uneven partitioning during division, or the stochastic expression of an inhibitor. Only certain cells in the pedigree would accumulate sufficient Wor1 protein to excite the Wor1 positive feedback loop, leading to high level Wor1 expression and subsequent cell type switching. This hypothesis is subject to explicit testing and will be the focus of future work. Our work has resolved several issues needed for a quantitative molecular understanding of white-opaque switching and lays the foundation for subsequent hypothesis tests. While white-opaque switching and the genetic network that controls it are highly complex, it represents 1 of only a few cases of a true epigenetic change in a eukaryotic cell that does not require external signals and is simple enough to be analyzed in detail and ultimately, understood in great depth.

Materials and methods C. albicans strain construction The genetic background for C. albicans is SC5314. SNY425 was used as a prototrophic a/α reference strain [58]. For a prototrophic a/Δ reference strain, the α copy of the MTL in SNY250 [59] was replaced with a copy of ARG4 using pJD1 [60]. Genetic modifications were created using a SAT1 marker-based CRISPR protocol targeting Candida maltosa LEU2 [61]. Homozygous Wor1-GFP and Wor1-mGFP fusion proteins were created at the endogenous WOR1 locus using a combination of 2 guides encompassing the stop codon (GGTACTTAGTTGAATTAATA and GTACTTAGTTGAATTAATAC). Donor DNA was created by restricting plasmids containing the GFP (or mGFP) sequence flanked by 500 bp of homology to the 3′ end of the WOR1 ORF and 3′ UTR [62]. The WOR1 motif deleted from the WOR1 promoter consists of 3 overlapping binding motifs (17 bp, AAAGTTTAAACTTTAAA), found 5,802 bp upstream of the translation start site. A combination of 2 guides were used for creating the deletion (GAAATTTATGAAAGACGGGT and ACCAGTAGTGATTCATAAAT), donor DNA consisted of 90 bp oligos containing 45 bp homology to flanking sequences with 2 single bp substitutions to disrupt PAM or guide sequences (AATAATTAGAGTTTTACAAGAAATTTATGAAAGACGGGTGcGAAATACTATAAAATAGAGCAATAAATGACAGAACCTAgCAGTAGTGAT). S. cerevisiae strain construction The genetic background for S. cerevisiae experiments is S288C. Starting strains were mating type a or α with multiple auxotrophies, created by mating FY23 (leu2Δ1 ura3-52 trp1Δ63) [63] and a mating type α spore of BY4743 (his3Δ1 leu2Δ0 lys2Δ0 met15Δ0 ura3Δ0) [64]. Subsequent transformations integrated various constructs in the following order: the chimeric transcription factor, deletion of the WOR1 homologs, the induced protein (GFP, Wor1-GFP, Wor1-mGFP, or Wor1), and the WOR1 transcriptional reporter. Plasmids were constructed using either restriction digest and subsequent ligation or by using homology-based cloning (In-Fusion, Takara 638911) and were sequence verified, all PCRs utilized a high fidelity polymerase (CloneAmp HiFi PCR Premix, Takara 639298). Plasmids contained components of the synthetic inducible system flanked by homology to the intended genomic integration locus and were restricted prior to transformation. The chimeric transcription factor used is “Z3PM” [39] that consists of the Zif268 DNA-binding domain, hPR ligand-binding domain and the yeast Msn2 activation domain. The transcription factor was expressed from the S. cerevisiae ADH1 promoter and terminated by the C. albicans ADH1 terminator. The transcription factor was integrated into leu2 proceeded by the Candida glabrata LEU2 sequence, including promoter and terminator. Induced proteins were expressed from a modified GAL1 promoter containing 3 dimeric Zif268 binding sites [39] and terminated by the C. albicans ADH1 terminator. The proteins were integrated into the trp1 locus proceeded by the C. glabrata TRP1 sequence, including promoter and terminator. The WOR1 homologs (MIT1 and ROF1) were knocked out using a plasmid-based CRISPR approach [65]. A plasmid with a 2 μ origin of replication was constructed that contained Cas9, a KanMX selection cassette and 2 guides (1 for each of the homologs). The guides were expressed from the SNR52 and SUP4 promoters, respectively, and both were terminated by a copy of the SNR52 terminator. Donor DNA consisted of 90 bp complementary annealed oligos [61] containing 44 bp of homology upstream and downstream of the ORF, flanking a GG sequence. After verifying both homologs were deleted, cells were grown overnight in the absence of G418, plated on YPD and replica plated to YPD+G418 to identify colonies that had lost the CAS9 containing plasmid. The main WOR1 transcriptional reporter consists of 6656 bp of the WOR1 promoter (ending at a presumed TATA box, 2048 bp before the translational start site), 115 bp containing the end of the CYC1 promoter and it is 5′ UTR (starting at the “‒52” TATA box [66]) and a fluorescent protein ORF (mCherry or GFP) followed by the C. albicans ACT1 terminator. An alternative reporter contained an additional 70 bp of the end of the WOR1 promoter and 1,980 bp encompassing the WOR1 5′ UTR instead of the CYC1 sequence. Transcriptional reporters were preceded by a hygromycin resistance cassette and integrated into the ura3 locus. Integration was achieved using 2 separate restricted plasmids simultaneously transformed. One plasmid contained homology upstream to the URA3 locus, the hygromycin cassette and the WOR1 promoter and beginning of 5′ UTR. The second plasmid contained the end of the WOR1 promoter, the CYC1 (or WOR1) 5′ UTR, the fluorescent protein and terminator and homology downstream to the URA3 locus. With 1 exception, combinations of constructs were represented by 2 completely independent strains based on different starting strains (of mating type a or α). The induced Wor1 (Δmit1/Δrof1) combination was represented by 2 mating type a strains. Microscopy and microfluidics C. albicans strains were streaked from frozen stocks on SD+aa+uri plates and incubated for 5 to 9 days at 25°C. Colonies were directly diluted to an OD of 0.25 and grown in liquid SD+aa+uri for 4 hours. Cultures were washed and diluted to approximately 5 × 106 cells/ml in PBS without calcium or magnesium, and 50 ul were loaded into microfluidic plates. For microfluidic experiments, we used a CellASIC ONIX microfluidic platform and corresponding Y04T-04 plates. Each plate contains 4 separate chambers, each chamber containing 104 individual traps. Images were captured using a Nikon Ti2-E microscope equipped with a Photometrics Prime 95B-25mm Camera and a CFI60 Plan Apochromat Lambda 100× Oil immersion objective lens. Cells were loaded into the plate chambers using 55.1 kPa of pressure for 5 seconds, and media (SD+aa+uri) flow was started immediately. The media flow program consisted of cycles of three 5-second bursts of perfusion at 10 kPa from each of 6 inlet channels followed by an hour of perfusion at 1.7 kPa from each of 6 inlet channels. Coordinates for fields containing single traps loaded with at least 1 cell were manually determined. Nikon NIS Elements software was used to drive stage movement and acquisition; the Nikon Perfect Focus System was utilized. Images (DIC and GFP channels) were captured for each field (40 to 70 per chamber) every 12 minutes for 24 hours. Temperature was kept at 25°C using an OKOLab Cage Incubator. Image and data analysis Custom Matlab scripts were used for image analysis; statistical analysis and data visualization utilized R [67], with extensive use of tidyverse [68]. Images were first adjusted to correct for stage movement. Each image was compared to the previous image in the time series, and the optimal geometric transformation was found to account for any (x,y) movement (Matlab function imregtform). Next, cells were segmented automatically in each image. Cell borders were identified using a combination of bottom hat filtering, top hat filtering, dilation, erosion, and selection based on object properties (area, perimeter, and eccentricity). Individual cells were tracked through time by manually determining their position in a series of images; each position was then associated with the corresponding automatically segmented cell. Pedigree information was also manually determined. For a metric of Wor1 expression per cell, the 300 highest intensity pixels (in the GFP channel) were averaged. Typical white cell size is 1,750 pixels; typical opaque cell size is 2,500 pixels. Cell length (Fig 2 and Fig B in S1 File) was determined using Fiji [69] by recording the length of a line drawn along the cells long axis. Cell length was recorded when cells were beginning to bud a new cell. A total of 100 references of white and opaque cells were recorded. Cell length was recorded for all pairs of switching cells across 3 experiments from different days, resulting in 248 mother-daughter pairs of cells containing Wor1-GFP and 50 pairs of cells containing Wor1-mGFP. For determining if mixed-fate pedigrees represent multiple independent switching events, we considered 168 traps loaded with cells containing Wor1-mGFP. The numbers of switching events were counted in each trap according to 3 different ways of calculating independent events: (1) every pair of mother-daughter cells are considered independent; (2) pairs of direct mother-daughter pairs are no longer considered independent; and (3) pairs in mixed-fate pedigrees are also no longer considered independent. We calculated the expected Poisson probability distribution based on the unambiguous zero-class of 134 traps (λ = ‒log (134/168)). We performed chi-squared goodness-of-fit tests with p-values computed by Monte Carlo simulations with 2,000 replicates. Flow cytometry, experiments, and data analysis S. cerevisiae strains were streaked from frozen stocks on YPD plates and incubated for 2 days at 30°C. Colonies were inoculated into 3 ml of SD+aa+uri media and grown overnight. Cells were washed, resuspended in PBS without calcium or magnesium, and diluted 1:16 into 96-well plates containing SD+aa+uri with different concentrations of progesterone. Progesterone was diluted from 2 mM stocks (95% ethanol, 5% DMSO); lower concentrations were supplemented with ethanol and DMSO to match the high concentration. Plates contained 8 experimental strains measured at 10 different concentrations of progesterone and 3 control strains measured at 4 concentrations of progesterone. For each set of strains, 2 identical plates were created and incubated, shaking at 500 rpm at 30°C. After a set incubation time, plates were washed 3 times and diluted in PBS without calcium or magnesium with 1 mM EDTA. One plate was read after 4.5 hours of incubation, and the second was read at 24 hours. Data from both plates were similar, and the 24-hour readings are presented. Measurements were taken using a BD FACSCelesta. GFP fluorescence was measured by blue (488 nm) excitation and a 530/30 band pass emission filter. mCherry fluorescence was measured by yellow-green (561 nm) excitation and a 610/20 band pass emission filter. Flow cytometry data were analyzed using R [67], with extensive use of tidyverse [68]. Data were imported into R using flowCore [70]. Data were filtered to remove potential cell debris, cell aggregates, and contaminants by constructing 2 gates, 1 based on forward and side scatter and a second based on GFP and a violet (405 nm) excitation, 670/30 emission channel. GFP and mCherry measurements were normalized by forward scatter. Presented data encompass 376 wells with a range of 127 to 8,905 events per well after filtering (median: 7,095.5; 95% of samples have over 4,000 events). For fitting Hill equations, we used the drm function in the package drc [71] using robust median estimation. mCherry measurements were log10 transformed and fitted to a 4-parameter log-logistic function with a parameterization converting the natural log of the half-max value into a parameter. For Wor1-GFP, only data from progesterone concentrations less than 2 nM were used for fitting.

Supporting information S1 File. A PDF file containing Supporting information figures (A–J) and table (A). Each figure is followed by a legend and data underlying the figures can be found at (https://osf.io/6e9vz/). Table A. Switching frequencies. Switching frequencies as calculated by counting sectors on colonies grown on solid agar plates. Data come from 3 independent experiments containing multiple plates per strain. Numbers of colonies and sectors are pooled across plates. The percent of colonies with at least 1 sector is calculated. Standard error is the standard error of a sample proportion (√ [p (1-p) / n)]). Fig A. Opaque cell growth in microfluidic traps. Two examples of microfluidic traps with growing opaque cells. Strain contains Wor1-GFP fusion protein. Fig B. Cell length distributions. (A) Cell lengths of pairs of switching cells (mothers versus their daughters), including switching cells containing Wor1-mGFP. Typical white and opaque cell distributions are shown on margins. (B) Cell length distributions (violin plots) of different cell types. For the Wor1-GFP strain, data come from 3 independent experiments that were pooled in main text Fig 2A. Fig C. Nucleus loss during mitosis in switching cells. (A) Representative example of a cell division where the mother cell nucleus enters the elongated daughter cell but returns to the mother cell. Strain contains Wor1-GFP fusion protein. (B) Representative example of cell divisions where the mother cell nucleus enters the elongated daughter cell and remains in the daughter cell, creating a polyploid cell. There are 2 cell divisions with this pattern in this example. Strain contains Wor1-GFP fusion protein. Fig D. Semiautomated image analysis pipeline. (A) Representative original images of 2 time points for 1 field. (B) Same images as (A) after adjustment for (x,y) movement (note black regions at image edge) and automated segmentation of cells (magenta colored lines defining cell borders). (C) Same images as (A) and (B) with overlay of manually defined cell identities, pedigree information is explicitly encoded in cell names (e.g., cell “2.1” is the first daughter of cell “2,” which, in turn, is the second daughter of cell “0”). Fig E. Fluorescence (GFP/mGFP) levels in single cells. (A) Additional Wor1 expression level traces in single cells. Traces are of white cells that do not contain Wor1-GFP (gray) and white cells that do contain Wor1-GFP (pink/purple). X-axis represents time from cell birth. Cells that are not actively switching have similar fluorescence levels to typical white cells that do not contain GFP. (B) The same examples of Wor1 activation in single switching cells shown in main text Fig 2D, lined up by the approximate time cells reach high Wor1 levels. Although the time between the birth of a switching cell and its switching varies, reaching high levels of Wor1 from background fluorescence levels takes approximately 3 hours. Fig F. Mixed-fate pedigree example, Wor1-GFP. (A) Representative pedigree in which multiple groups of cells activated Wor1. Horizontal lines represent single cells; vertical lines represent budding of a daughter cell. Every horizontal line is made up of small tiles colored by Wor1-GFP fluorescence. Numbers within parentheses on the right of the pedigree represent the number of budded daughter cells per cell within the time period shown; not all daughters are depicted in the pedigree. An asterisk represents that the cell lost its nucleus to its daughter cell. Colored square tiles single out particular cells that are depicted in (B). (B) Subset of images representing data shown in (A). Particular cells are outlined in colors corresponding to (A). Cell outlines are based on automated image analysis. Fig G. Mixed-fate pedigrees. Schematics of all mixed-fate pedigrees observed for the Wor1-mGFP strain. Numbers represent the number of cell divisions separating switching cells from the establishment of the predisposed pedigree. Representation is for illustrating cell relationships and do not accurately reflect cell division times or timing/intensity of Wor1 activation. Fig H. Alternative Wor1 transcriptional reporters in S. cerevisiae. Distributions of Wor1 transcriptional reporters; all strains shown are inducing Wor1 protein with increasing progesterone concentrations and were independently constructed from the strains shown in the main text. “Full” reporter contains both the Wor1 7 kb control region and promoter and the Wor1 2 kb 5′ UTR. The “No-UTR” reporter contains the Wor1 control region and the Cyc1 core promoter and 5′ UTR as explained in the Methods section. Fig I. Distributions of Wor1 protein and reporter levels. Distributions of GFP (A) and mCherry (B) fluorescence; data are shown for all strains in a subset of 4 hormone concentrations (nM progesterone). Note bimodal distributions of GFP fluorescence when inducing Wor1-GFP at high hormone concentrations in (A). Note low-level expression of the Wor1 transcriptional reporter compared to autofluorescence in (B). Fig J. Distributions of Wor1 protein and reporter levels. (A) Same data as main text Fig 4C, separated by hormone concentration (nM progesterone). (B) 2D density plots of flow cytometry data. S. cerevisiae Wor1-GFP and Wor1-mGFP data (same as A) are shown alongside data for white and opaque C. albicans cells containing a Wor1-GFP fusion protein and a transcriptional reporter encompassing the 7 kb control region of Wor1 followed by mCherry (no 5′ UTR sequence). https://doi.org/10.1371/journal.pbio.3001657.s001 (PDF) S1 Video. Example of time-lapse movie of a switching event Movie corresponds to “Example 1” in Fig 1D. Movie is shown twice; in the second repetition, tracked cells are labeled. This strain contains the Wor1-GFP fusion protein. https://doi.org/10.1371/journal.pbio.3001657.s002 (MP4) S2 Video. Example of time-lapse movie of a switching event containing a mixed-fate pedigree Movie corresponds to example in Fig 3C. Movie is shown twice; in the second repetition, tracked cells are labeled. This strain contains the Wor1-mGFP fusion protein. https://doi.org/10.1371/journal.pbio.3001657.s003 (MP4)

Acknowledgments We thank members of the Johnson lab, in particular, Matt Lohse and Kyle Fowler for helpful discussions and Ananda Mendoza for technical support. We thank David Gresham, Kerry Geiler-Samerotte, Mariana Gomez Schiavon, and Ranen Aviner for advice and the lab of Hana El-Samad and Dmitri Petrov for providing reagents.

[END]

[1] Url: https://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.3001657

(C) Plos One. "Accelerating the publication of peer-reviewed science."
Licensed under Creative Commons Attribution (CC BY 4.0)
URL: https://creativecommons.org/licenses/by/4.0/


via Magical.Fish Gopher News Feeds:
gopher://magical.fish/1/feeds/news/plosone/