As animal cells develop, they pass through different states to mature into specific cell types. For some cells, this development depends on the cell’s ability to switch between two stable states, a property called bistability.

Many bistable systems operate during development and often feature proteins called transcription factors that regulate a few cell states in specific tissues. These proteins activate or repress a range of target genes, and cells can adjust the levels of the transcription factors to bring about transitions between states. In fruit flies, two transcription factors, called Yan and Pnt, regulate numerous processes throughout development.

In the developing embryo of a fruit fly, Yan and Pnt are regulated by signals induced by a receptor called EGFR. When EGFR is activated, Pnt is produced and Yan is degraded. When this receptor is not activated, Yan is produced and represses Pnt. Mathematical modelling of these interactions has concluded that this is a bistable system: that is, cells should either have high levels of Yan and low levels Pnt, or low Yan levels and high Pnt levels. However, larval eye cells first activate, and then turn off, both proteins together. This argues against bistability and raises questions about how these proteins regulate cell fate transitions in the eye, and perhaps in other organs.

To investigate this question, Peláez et al. tagged Yan with a fluorescent marker to track its activity in the eyes of fruit fly larvae as they develop. A combination of fluorescence-based microscopy and an automated imaging analysis system were then used to score fluorescence in thousands of individual eye cells and assess the changes in the levels of Yan over time. This approach revealed some unexpected results.

Yan levels were seen to vary in both immature and maturing cells. Thus eye cells transition between states as Yan levels increase rapidly, suggesting the need for a mechanism that is distinct from bistability. Peláez et al. suggest that larger changes in the seemingly random fluctuations in the levels of the transcription factor (also known as “expression noise”) might play a role in this mechanism. In particular, Yan expression noise briefly peaks as cells transition to a more mature state, and the transient nature of this ‘noisy’ response requires the activation of EGFR.

One possible explanation for these observations is that Yan’s effect on cell states depends on this variability in its levels, which might prime cells to change states when they receive another signal. These findings also raise many questions for future studies to explore, including how this increase in the noise level of Yan is triggered as cells begin their transitions towards specific cell types.

Cells within complex organisms exist in different states that confer specific functionalities to each cell. Cellular states can be organized into cyclic cascades such as the G1-S-G2-M cell cycle, or into linear cascades as observed in cell differentiation. A common feature of cells that undergo state transitions is the apparent irreversibility of the transitions even when such transitions are triggered by transient stimuli. Modeling of these transitions often assumes system bistability, in which cells can be resting in one of two stable states.

Animal cells frequently utilize transcription factors to enforce a given state, and transitions to another state are driven by increasing or decreasing the levels of these transcription factors (Yao et al., 2008; Kueh et al., 2013; Laslo et al., 2006; Park et al., 2012). The Rb-E2F pathway generates bistability in E2F expression, which dictates the transition from G1 to S phase (Yao et al., 2008). Expression of the transcription factor PU.1 determines lymphoid versus myeloid hematopoietic cell lineages (Kueh et al., 2013; Laslo et al., 2006). Adipocyte differentiation is controlled by differential expression of C/EBP and PPARγ proteins (Park et al., 2012). Regulation by positive feedback is a hallmark of bistable systems, and in all of the above cases, the transcription factors act in one or more positive feedback circuits. In some systems, positive feedback is generated by two transcription factors that mutually repress each other’s expression. In these scenarios, cells in one state continually express a transcription factor that represses a second transcription factor, and when these cells transit to another state, they continually express the second transcription factor, which represses its antagonist. Fate restriction in hematopoietic, neural, pancreatic, and muscle cell lineages is regulated by such double-negative feedback circuits (Briscoe et al., 2000; Laslo et al., 2006; Olguin et al., 2007; Schaffer et al., 2010).

Most examples of this type of bistable mechanism have transcription factors that act specifically within a handful of cell states limited to a single tissue or organ system. One remarkable exception to this rule is found in Drosophila. There, two ETS-domain transcription factors act in a wide assortment of cell types across the body and across the life cycle. Yan and Pointed (Pnt) act downstream of signals mediated by receptor tyrosine kinases (RTKs) that regulate cell differentiation, migration, and division in tissues ranging from ovarian follicular cells, dorsal and ventral neuroectoderm, embryonic mesoderm, the embryonic trachea, and the post-embryonic compound eye (Dumstrei et al., 1998; Gabay et al., 1996; Halfon et al., 2000; Jurgens et al., 1984; Morimoto et al., 1996; O'Neill et al., 1994; Ohshiro et al., 2002; Yao et al., 2008) It is thought that Yan and Pnt control such diverse cell states by acting in concert with tissue-specific transcription factors to regulate transcription of appropriate target genes. For instance, transcription of even-skipped occurs in mesoderm only if Yan/Pnt act in combination with Tinman and Twist proteins (Halfon et al., 2000), whereas transcription of prospero (pros) in the eye only occurs if Yan/Pnt act in combination with Lozenge, Sine Oculis, and Glass proteins (Hayashi et al., 2008; Xu et al., 2000).

Several tissues show mutually exclusive expression of Yan and Pnt, suggestive of cross-repression (Boisclair Lachance et al., 2014). The best characterized is the embryonic ventral ectoderm, where ventral-most cells express Pnt and more lateral cells express Yan (Gabay et al., 1996). This pattern is established by secretion of a ligand for the EGF Receptor (EGFR) from the ventral midline (Golembo et al., 1996). EGFR activation in nearby ventral cells induces the Ras-MAPK pathway to express Pnt and degrade Yan (Gabay et al., 1996; Melen et al., 2005). Cells with insufficient EGFR activation express Yan, which represses Pnt. Mathematical modeling has described this as a bistable system in which cells are either in a High Yan/Low Pnt state or a Low Yan/High Pnt state (Graham et al., 2010; Melen et al., 2005). Transition from one state to the other is ultrasensitive to the strength of EGFR activation (Melen et al., 2005).

Paradoxically, other Drosophila tissues show co-expression of Yan and Pnt (Boisclair Lachance et al., 2014). The larval eye is one such tissue. Retinal progenitor cells initiate expression of both proteins, and when they transit to differentiated photoreceptor fates, these cells reduce expression of both proteins. In contrast, when retinal progenitor cells transit to differentiated cone cell fates, they maintain their expression of both proteins. These observations are at odds with long-standing genetic studies that support a standard bistable mechanism acting in the eye (Lai and Rubin, 1992; O'Neill et al., 1994; Rebay and Rubin, 1995). Thus, new approaches to studying these transitions in the eye are needed.

Here, we have adopted a systems-level approach to study Yan dynamics in the larval eye. A yellow fluorescent protein (YFP) — based isoform of Yan was developed as a reporter for Yan protein levels. Fluorescence-based microscopic imaging of cells was coupled with automated high-throughput image analysis to score fluorescence in each cell and annotate the data in a quantitative and unbiased fashion. Yan exhibits monostability, both in progenitor and differentiating cells, with Yan levels varying in cells in either state. Cell state transitions occur independent of absolute Yan concentrations, suggesting that some other mechanism allows Yan to regulate transitions. One such mechanism might be the noise in Yan levels, which undergoes a transient spike as cells begin to transition to differentiated states. Loss of EGFR signaling, which prevents cells from differentiating, causes these cells to have prolonged noisy Yan expression, and suggests that Yan noise is key for cell state transitions in the eye.

The compound eye epithelium is established during embryogenesis as an internal disc of cells called the eye imaginal disc (Wolff and Ready, 1993). During the larval phase of the life cycle, the disc grows in size by asynchronous cell division. During the final 50 hr of the larval phase, a morphogenetic furrow (MF) moves across the eye disc from posterior to anterior (Figure 1A,B). All cells arrest in G1 phase within five cell diameters anterior to the furrow, and then as the furrow passes through them, periodic clusters of cells express the proneural gene atonal (Jarman et al., 1994). Atonal expression is subsequently restricted to one cell per cluster, which becomes the R8 photoreceptor. Each R8 cell then secretes an EGFR ligand that activates the receptor in neighboring cells and causes them to transit from multipotent progenitor to differentiated states (Figure 1C) (Freeman, 1996). Transitions occur in a sequence of symmetric pairs of multipotent progenitor cells that differentiate into R2/R5, R3/R4, and R1/R6 photoreceptors (Figure 1C) (Wolff and Ready, 1993). Thereafter, a single progenitor transits to a R7 photoreceptor fate followed by two pairs of cells, C1/C2 and C3/C4, that differentiate into cone cells. These cone cells are non-neuronal and form the simple lens that overlies each cluster of eight photoreceptors. The furrow induces the nearly simultaneous differentiation of a column of R8 cells, with repeated column inductions producing approximately 800 units or ommatidia as the furrow moves across the eye.

Figure 1 with 4 supplements see all

Download asset Open asset

A central tenet of the bistable model of cell differentiation in the eye posits that differentiation is marked by a transition from high Yan protein levels in multipotent progenitor cells to low Yan levels in differentiating cells (Graham et al., 2010). Formulation of this model stemmed from studies of R7 cell differentiation, the final photoreceptor recruited to each ommatidium. Reduced Yan causes inappropriate expression of the R7 determinant pros and ectopic R7 cells in yan hypomorphic mutants (Kauffmann et al., 1996; Lai and Rubin, 1992). Conversely, a Yan isoform that is resistant to MAPK-dependent degradation, blocks R7 differentiation and pros expression (Kauffmann et al., 1996; Rebay and Rubin, 1995).

Quantifying Yan dynamics

To quantitatively test the bistable model, we used BAC recombineering to insert fast-fold yellow fluorescent protein (YFP) in-frame at the carboxy-terminus of the yan coding sequence (Webber et al., 2013). The Yan-YFP transgene fully complemented null mutations in the endogenous yan gene and completely restored normal eye development (Figure 1—figure supplement 1), demonstrating that the YFP tag does not compromise Yan function and that all essential genomic regulatory sequences are included. The pattern of Yan-YFP protein expression was qualitatively similar to that of endogenous Yan (Figure 1—figure supplement 1).

We used histone His2Av-mRFP fluorescence in fixed specimens to mark all eye cell nuclei for automated segmentation (Figure 1D—figure supplement 2). Nuclei were manually classified into the different cell types of the eye, which is possible because every cell undergoing differentiation can be unambiguously identified by its position and nuclear morphology without the need of cell-specific markers (Ready et al., 1976; Tomlinson, 1985; Tomlinson and Ready, 1987; Wolff and Ready, 1993). To validate our identification of all cell types using this method, we cross-checked with specific cell-specific markers, and found that our classification was accurate over 98% of the time (Figure 1—figure supplement 3). Cells were scored for nuclear Yan-YFP concentration and their exact position within 3D coordinate space of each fixed eye sample (Figure 1—figure supplement 2). Yan-YFP protein levels were normalized to His2Av-mRFP, which provided some control over measurement variation (Figure 1—figure supplement 4). We then mapped each cell’s spatial position in the x-y plane of the eye disc onto time. Two reasons made this possible. First, the furrow moves at an approximately constant velocity forming one column of R8 cells every two hours (Basler and Hafen, 1989; Campos-Ortega and Hofbauer, 1977). Second, each column of R8 cells induces the other cell fates at a constant rate (Lebovitz and Ready, 1986). Thus even in a fixed specimen, the temporal dynamics of cell state transitions are visible in the repetitive physical organization of ommatidia relative to the furrow (Figure 1C). Together these features allow the estimation of time based on a cell’s position relative to the furrow (Figure 2—figure supplement 1). This can be applied repeatedly to hundreds of cells in a single sample, creating a composite picture of the dynamics (Figure 2B–F). Although the developmental progression of an individual cell cannot be measured by this approach, it nevertheless provides a dynamic view of hundreds of cells across a developing epithelium as they respond to signaling processes. From this information, individual cell behaviors can be reconstructed and modeled.

Figure 2 with 1 supplement see all

Download asset Open asset

Yan-YFP expression in multipotent progentior cells was biphasic (Figure 2B). Cells anterior to the furrow expressed a basal level of Yan-YFP, but this level dramatically increased in cells immediately anterior to the furrow, starting eight hours before the first R8 cells were identifiable. Approximately eight hours after R8 definition, Yan-YFP levels peaked, and thereafter gradually decayed until Yan-YFP approached its basal level again. The results are surprising in two ways. First, Yan-YFP is not maintained at a stable steady state within progenitor cells, which would have been predicted by the bistable model. Rather, its dynamics are reminiscent of monostable responses to transient stimuli, with a single basal steady state. Second, at the peak of Yan-YFP expression, there is remarkably large heterogeneity in Yan-YFP levels across cells.

We also followed Yan-YFP dynamics in cells as they transited into a differentiated state and thereafter. Again, the results did not follow the expectations predicted by the bistable model. First, progenitors transited to a differentiated state at levels of Yan-YFP that varied, depending upon the type of differentiated state being adopted (Figure 2C–H). Cells entering the R3/R4 and R1/R6 states began with Yan-YFP levels that were approximately two-fold greater than cells entering the R2/R5 states. Cells entering the R7 state were intermediate between these two extremes. Despite these differences at the transition point, Yan-YFP levels decayed to a similar basal steady state irrespective of the photoreceptor type, and this basal level was at least three-fold lower than that which the progenitor cells achieved (Figure 2G). Thus, rather than the single high Yan progenitor state previously modeled, our results suggest a dynamic range of high Yan states from which different cell fates are specified according to the spatio-temporal sequence of differentiation.

We noted that for most photoreceptors, it took approximately 20 hr for Yan-YFP to decay to the basal steady state (Figure 2G), significantly longer than had been previously modeled (Graham et al., 2010). Since expression of several neural-specific genes is detected 2–8 hr after the transition (Tomlinson and Ready, 1987; Van Vactor et al., 1988), the slower than anticipated Yan decay indicates that early differentiation does not require cells to have assumed a basal steady-state level of Yan-YFP.

The last group of progenitors to differentiate form the non-neuronal cone cells. We also measured Yan-YFP in those cells. Yan-YFP dynamics in cone cells were more similar to progenitor cells over the same time period (Figure 2F,H). This behavior was in contrast to photoreceptors, which exhibited different decay dynamics from progenitor cells. Thus, accelerated degradation of Yan-YFP is not essential for cells to transition to all retinal cell states.

EGFR-ras signaling regulates Yan-YFP dynamics

The bistable model posits that the switch from one state to another is triggered by a signal that progenitor cells receive though the EGFR protein. Given the unanticipated Yan-YFP dynamics, we wanted to ask whether and how they were influenced by EGFR signaling. EGFR null mutants are inviable; however, a temperature sensitive (ts) allele of EGFR enables controlled inactivation of the RTK’s activity (Kumar et al., 1998). We grew EGFR^ts mutant larvae at a restrictive temperature for 18 hr before analyzing effects on Yan-YFP. Surprisingly, progenitor cells exhibited biphasic expression of Yan-YFP over time, but the amplitude of the pulse in expression was significantly reduced (Figure 3A). This suggests that EGFR signaling contributes to the stimuli that induce the Yan-YFP peak. To further test this hypothesis, we misexpressed a constitutively active form of Ras protein in eye cells. The peak of Yan-YFP in progenitors was now prolonged (Figure 3B). Together, these results suggest that EGFR-Ras signaling stimulates the transient appearance of Yan-YFP in progenitor cells, and that the decline in Yan-YFP within older progenitor cells is linked to a loss of signal reception by these cells over time.

Figure 3 with 1 supplement see all

Download asset Open asset

We next examined the effects of EGFR and Ras on Yan-YFP dynamics in cells as they differentiate. The bistable model predicts that EGFR is required for the loss of Yan-YFP in photoreceptors. Indeed, EGFR^ts mutant R2/R5 cells delayed their initial decline in Yan-YFP levels (Figure 3C). Conversely, misexpression of constitutively active Ras caused a premature decline in Yan-YFP (Figure 3D). These results are consistent with EGFR-Ras stimulating the degradation of Yan-YFP as cells switch their states. However, Yan-YFP dropped to below-normal levels in EGFR^ts mutant R2/R5 cells (Figure 3C). These complex effects suggest a dual role for EGFR signaling in photoreceptors. In the first few hours as cells transit to a photoreceptor state, EGFR stimulates the accelerated decay of Yan-YFP. Thereafter, EGFR inhibits the decay of Yan-YFP in a manner that might be related to that role that EGFR plays in boosting Yan-YFP levels in progenitor cells.

The source of EGFR ligand originates from the R8 cell (Freeman, 1996). If this diffusive ligand is responsible for controlling Yan-YFP levels in other photoreceptors as they are recruited to an ommatidium, we would predict a correlation between Yan-YFP levels in differentiating cells and their distances from adjacent R8 cells. Indeed, at the times when R2/R5 cells differentiate (~0–15 hr), their Yan-YFP levels are positively correlated (p<0.01) with their physical distance to the nearest R8 cells (Figure 3—figure supplement 1). These correlations are absent in EGFR^ts mutants, providing evidence that R8 cells act through EGFR to control Yan-YFP dynamics in differentiating cells.

Cells are indifferent to absolute levels of Yan

The bistable model predicts that different cell states depend upon discrete absolute concentration of Yan present inside cells. To test this idea, we varied the number of Yan-YFP gene copies. In general, protein output is proportional to gene copy number in Drosophila (Lucchesi and Rawls, 1973). We increased Yan-YFP copy number from its normal diploid number to tetraploid, and monitored Yan-YFP in progenitors and differentiating cells. As expected, four copies caused a higher steady-state level of Yan-YFP in progenitor cells, though this increase was less than two-fold (Figure 4A). The amplitude of the Yan-YFP pulse was also increased as progenitor cells aged. Strikingly, as four-copy progenitor cells transited to a differentiated state, the onset of Yan-YFP decay occurred at the same time as it occurred for two-copy cells (Figure 4A–C). Yan-YFP levels were much greater in four-copy cells compared to two-copy cells making their transit into the identical cell states. To confirm that absolute Yan-YFP concentration had little effect on cell state transitions, we examined expression of a direct target of Yan in R7 cells: the pros gene (Kauffmann et al., 1996; Xu et al., 2000). Expression was monitored with an antibody specific for Pros protein. We observed at most a one hour delay in the onset of Pros expression in R7 cells containing four copies of Yan-YFP (Figure 4D), far less than the ten-hour delay predicted if absolute concentration of Yan-YFP dictated when Pros expression begins (Figure 4D).

Figure 4 with 2 supplements see all

Download asset Open asset

Possibly, absolute concentration of Yan is unimportant when a cell transits to a different state, but the switch to a constant basal Yan level is robustly maintained regardless of starting concentration. An examination of Yan-YFP decay in photoreceptor cells makes that possibility less likely; Yan-YFP decays to different basal levels in two- versus four-copy differentiated cells (Figure 4A–C). To further test this notion, we fit the data to several plausible functional forms. We found that exponentially decaying functions systematically best-fit to the data (Figure 4—figure supplement 1). Thus, for each cell state we fit an exponential decay function to its Yan-YFP temporal profile (Figure 4—figure supplement 2). From these fits, we derived the average decay rate constants and half-lives of Yan-YFP for cells carrying two, four, or six copies of yan. As expected, we found that Yan-YFP half-life was different between progenitors and differentiating photoreceptors (Figure 5—figure supplement 1). The half-life in photoreceptors was two-fold lower than in progenitors, accounting for the more rapid loss of Yan-YFP in the former cells. Strikingly, Yan-YFP half-life was not significantly affected by yan copy number within either progenitor or photoreceptor cells (Figure 5). Thus, Yan-YFP concentration only affected its decay rate as a first-order function, implying that there is no higher order mechanism to accelerate Yan-YFP decay when cells contain greater concentrations of Yan-YFP.

Figure 5 with 2 supplements see all

Download asset Open asset

As a final test of the effects of Yan-YFP levels on cell states, we allowed 4X and 6X yan animals to complete eye development and then examined the external patterning of the fully differentiated compound eye. The compound eyes were completely normal in appearance (Figure 5—figure supplement 2), implying that differentiation was unaffected by the absolute concentrations of Yan inside eye cells.

Yan expression noise spikes during cell state transitions

Our results indicate that Yan’s effects on retinal cell states are not dictated by uniform and stable concentrations of Yan protein. One potential explanation is that Yan’s effect on cell states actually depends on variability in Yan protein levels. A growing body of evidence is pointing to the importance of transient fluctuations in expression of factors to control cell states (Cahan and Daley, 2013; Torres-Padilla and Chambers, 2014). Key regulators of stem cells fluctuate in abundance over time, and this is correlated with stem cells existing in multiple connected microstates, with the overall structure of the population remaining in a stable pluripotent macrostate (MacArthur and Lemischka, 2013). Heterogeneity among cells is not simply due to independent noise in expression of individual genes but is due to fluctuating networks of regulatory genes (Chang et al., 2008; Kumar et al., 2014). Such fluctuations appear to be stochastic in nature, priming cells to transit into differentiated states with a certain probability that is guided by extrinsic signals.

Our data revealed considerable heterogeneity in the level of Yan-YFP among cells at similar developmental times (Figure 2B). To quantify the noise, we used developmental time to bin cells of similar age, and measured detrended fluctuations of Yan-YFP by calculating residuals to the fitted function and normalizing binned residuals to the function (Goldberger et al., 2002). Progenitor cells showed a spike in Yan-YFP noise as they began to induce Yan-YFP expression (Figure 6A). The noise spike was short-lived (approximately 17 hr), and noise thereafter returned to a basal level with secondary minor spikes. The major peak in noise magnitude coincided with the time at which R8 cells are formed.

Figure 6 with 4 supplements see all

Download asset Open asset

Photoreceptor cells showed a large spike in Yan-YFP noise as they began to differentiate (Figure 6B–E). The magnitude of each noise peak varied with the photoreceptor cell state; R3 and R4 cells exhibited the greatest amplitude in noise (Figure 6F). These noise spikes showed a distinct temporal relationship, with spikes coinciding with the times at which individual cell states were switched (Figure 6F). Thus, the noise spikes are not a simple consequence of a global stimulus synchronously affecting noise in all cells. Thereafter, all cells reduced Yan-YFP noise to a basal level that was comparable to basal noise in the progenitor cells. However, each cell type exhibited a secondary minor spike at 30–35 hr, which might reflect a synchronous stimulus.

Detrended fluctuation is one method to measure expression variability, but it can suffer from distortion if the model fitting is not adequately weighted. Therefore, we also calculated the coefficient of variation, that is, the standard deviation of Yan-YFP intensity within a sliding window divided by its mean. This method yielded noise profiles with transient spikes for each cell type that was consistent with calculations using detrended fluctuation (Figure 6—figure supplements 1 and 2). However, while the coefficient of variation yielded results that varied strongly with bin width, the detrended fluctuations yielded profiles that were generally robust against changes in bin width (Figure 6—figure supplement 1).

To rule out the possibility that these unexpectedly dynamic features of Yan-YFP were caused by its transgenic origins or fusion with YFP, we compared Yan-YFP dynamics to those of endogenous Yan protein that was bound with an anti-Yan antibody. The profiles of Yan-YFP protein levels and noise were highly similar to endogenous Yan protein levels and noise, in both multipotent and differentiating cells (Figure 6—figure supplement 3). Thus, transient spikes of expression heterogeneity are a fundamental feature of Yan protein.

Since Yan regulates Pros expression in R7 cells, it was possible that Pros showed a transient noise spike as a consequence. Therefore, we measured Pros protein heterogeneity and found that its dynamics did not resemble that of Yan (Figure 6—figure supplement 4). Instead, Pros noise was high starting at the onset of expression, and thereafter gradually declined as Pros protein levels increased in R7 cells. We conclude that noise spikes are not a general feature of gene expression in the developing eye but might reflect unique roles of Yan in mediating cell state transitions.

We wondered what might cause these spikes in Yan-YFP noise during cell state transitions. Because EGFR signaling is important for regulating Yan-YFP concentration during these transitions (Figure 3), we analyzed Yan-YFP noise when EGFR signaling was inhibited in EGFR^ts mutant animals raised at a non-permissive temperature. The noise spike in progenitor cells was not significantly affected by loss of EGFR signaling (data not shown). We also examined the effects of EGFR^ts on noise in differentiating photoreceptor cells. Interestingly, noise increased at the normal time of transition but the elevated noise did not quickly drop to basal levels (Figure 6G). Rather, high noise was extended for an additional 10 to 15 hr. Conversely, misexpressing constitutively active Ras within differentiating cells caused a premature dropdown in Yan-YFP noise (Figure 6H). These results indicate that EGFR/Ras signaling is required for the rapid drop in Yan-YFP noise after it has peaked, creating a transient spike.

This study relied upon a set of methods that enable systems-level analysis of Yan expression dynamics in a developing animal tissue. Transgene recombineering was used to insert YFP into a genomic rescue fragment of the yan gene, which fully replaced endogenous yan. Yan-YFP protein was quantified in thousands of individual cells by fluorescence confocal microscopy and automated cell segmentation analysis. Based on the unique features of Drosophila eye development, every analyzed cell was assigned an age, and composites of cells across a time-spectrum of ages were derived. This allowed us to reconstruct the temporal dynamics of Yan protein expression in cells as they transited from one state to another or were maintained in a given state. The fact that both Yan concentration and noise were easily measured using our approach indicates that it provides a powerful method for studying how other molecular determinants regulate cell states.

Contrary to what is currently believed, the expression of Yan in progenitor cells has many hallmarks of monostability. A stable basal state exists in cells anterior to the furrow, and when the furrow passes, Yan rises and falls to form a biphasic profile (Figure 7A). If cells transit towards differentiation, then the fall in Yan is accelerated but the fundamental biphasic profile is preserved. This monostable-like behavior is not like a behavior where progenitor cells exist in a high Yan stable-state and switch to a low Yan stable-state when they transit towards differentiation (Figure 7A). Other lines of evidence also point away from a bistable switch mechanism. Yan reaches its basal steady state many hours after cells have adopted their differentiated photoreceptor state and are executing specialized gene expression programs. Thus, Yan levels are variable at the time when cells actually become differentiated.

Figure 7

Download asset Open asset

Since cell states are indifferent to Yan stability, it would suggest that absolute Yan levels do not dictate Yan’s effects on cell behavior. This conclusion is bolstered by experiments in which yan gene dosage affects absolute Yan levels but not cell behavior. Several mechanisms could explain how Yan regulates cells in a concentration-independent manner. First, total Yan might vary but a specific modified form of Yan might remain constant. For example, MAPK phosphorylation of Yan could operate under Michaelis-Menten saturation to generate constant levels of phospho-Yan that depend upon MAPK activity. Second, Yan’s transcriptional activities might be independent of Yan concentration due to limiting levels of other transcription factors such as Lozenge, Glass, and Sine Oculis, which are known to act combinatorially with Yan to regulate gene transcription (Flores et al., 2000; Hayashi et al., 2008; Xu et al., 2000). Third, cells might sense relative changes in Yan, and respond to a constant fold-change in Yan levels using integral negative feedback or incoherent feedforward loops (Cohen-Saidon et al., 2009; Goentoro and Kirschner, 2009; Goentoro et al., 2009; Wartlick et al., 2011). Indeed, Yan is predicted to function in an incoherent feedforward loop with Pnt. Pnt directly activates transcription of target genes such as pros and mir-7, but it also activates expression of Yan (Figure 3G), which in turn, directly represses transcription of pros and mir-7 (Li et al., 2009; Xu et al., 2000).

The most intriguing hypothesis, however, is that variation in Yan concentration might actually be exploited by cells to guide their state transition. Indeed, it has recently been noted that undifferentiated stem cells express highly variable levels of transcription factors such as Nanog, Myc, Otx2, and Pax6 (Cahan and Daley, 2013; Kumar et al., 2014; Torres-Padilla and Chambers, 2014). Conditions that repress variability keep these cells trapped in an undifferentiated state (Kumar et al., 2014). Noisy expression is thought to place cells into a transition state where a sub-population of cells at any given time express a particular combination of factors that spontaneously trigger differentiation (Martinez Arias and Brickman, 2011). Such a fluctuating transition state renders the probability of cell differentiation to be greater than zero, and thus provides cells with the ability to respond to extrinsic signals.

Consistent with this hypothesis, Yan levels not only vary over long time-scales but also show heterogeneity within a small time interval. This expression noise is reminiscent of the heterogeneity in Nanog expression observed in pluripotent stem cells. As progenitor cells transit to differentiated states, there is a sharp spike in Yan expression noise that coincides with the early stages of the transition. By analogy to the situation in embryonic stem cells, the spike in Yan noise could reflect the existence of a metastable intermediate state between the undifferentiated and differentiated states (Figure 7B).

A limitation to our interpretation is that because our data do not track individual cells over time, we cannot distinguish whether the expression heterogeneity reflects asynchronous fluctuations within individual cells or a wide range of stable microstates. Real-time single cell measurements will be required to distinguish between these possibilities.

The purpose of the elevated noise could be to prime cells to change states upon reception of a differentiation signal. Without elevated noise, differentiation or de-differentiation would be blocked. A prediction of this hypothesis is that differentiation signals would act downstream of noise elevation and might even be important for noise reduction once cells enter an irreversible differentiated state. Indeed, EGFR/Ras signaling is not required for elevation of Yan noise but rather is important in promoting the rapid reduction of noise back to its baseline level after differentiation. In the absence of signaling, the cell may remain trapped in the intermediate (undifferentiated) state.

The model raises several questions. What triggers the increase in the noise level of Yan as cells enter this transition state? Do other determinants such as Pnt similarly fluctuate? Do individual cells fluctuate over time or does this represent a population of cells that adopt stable microstates? And finally, there is the question of the noise spike in progenitor cells at the time of R8 formation. Does this represent another transition state that early progenitor cells enter before moving into a different progenitor state? One intriguing possibility is that it represents the transition from G1-arrest to synchronized cell division that progenitor cells undergo in the second mitotic wave immediately posterior to the furrow (Wolff and Ready, 1993). Loss-of-function studies have suggested that Yan prevents progenitors from undergoing multiple rounds of cell division rather than again arresting at G1 phase after one round of division (Rogge et al., 1995). Further studies should help elucidate the causes and functions of this and the other noise peaks.

1. 1. Basler K
   2. Hafen E

   (1989)

   Dynamics of Drosophila eye development and temporal requirements of sevenless expression

   Development 107:723–731.

   • Google Scholar
2. 1. Boisclair Lachance JF
   2. Peláez N
   3. Cassidy JJ
   4. Webber JL
   5. Rebay I
   6. Carthew RW

   (2014) A comparative study of Pointed and Yan expression reveals new complexity to the transcriptional networks downstream of receptor tyrosine kinase signaling

   Developmental Biology 385:263–278.

   https://doi.org/10.1016/j.ydbio.2013.11.002

   • Google Scholar
3. 1. Briscoe J
   2. Pierani A
   3. Jessell TM
   4. Ericson J

   (2000) A homeodomain protein code specifies progenitor cell identity and neuronal fate in the ventral neural tube

   Cell 101:435–445.

   https://doi.org/10.1016/S0092-8674(00)80853-3

   • Google Scholar
4. 1. Brunner D
   2. Oellers N
   3. Szabad J
   4. Biggs WH
   5. Zipursky SL
   6. Hafen E

   (1994) A gain-of-function mutation in drosophila MAP kinase activates multiple receptor tyrosine kinase signaling pathways

   Cell 76:875–888.

   https://doi.org/10.1016/0092-8674(94)90362-X

   • Google Scholar
5. 1. Cahan P
   2. Daley GQ

   (2013) Origins and implications of pluripotent stem cell variability and heterogeneity

   Nature Reviews Molecular Cell Biology 14:357–368.

   https://doi.org/10.1038/nrm3584

   • Google Scholar
6. 1. Campos-Ortega JA
   2. Hofbauer A

   (1977) Cell clones and pattern formation: on the lineage of photoreceptor cells in the compound eye of Drosophila

   Wilhelm Roux's Archives of Developmental Biology 181:227–245.

   https://doi.org/10.1007/BF00848423

   • Google Scholar
7. 1. Chang HH
   2. Hemberg M
   3. Barahona M
   4. Ingber DE
   5. Huang S

   (2008) Transcriptome-wide noise controls lineage choice in mammalian progenitor cells

   Nature 453:544–547.

   https://doi.org/10.1038/nature06965

   • Google Scholar
8. 1. Clifford RJ
   2. Schüpbach T

   (1989)

   Coordinately and differentially mutable activities of torpedo, the Drosophila melanogaster homolog of the vertebrate EGF receptor gene

   Genetics 123:771–787.

   • Google Scholar
9. 1. Cohen-Saidon C
   2. Cohen AA
   3. Sigal A
   4. Liron Y
   5. Alon U

   (2009) Dynamics and variability of ERK2 response to EGF in individual living cells

   Molecular Cell 36:885–893.

   https://doi.org/10.1016/j.molcel.2009.11.025

   • Google Scholar
10. 1. Dominguez M
    2. Hafen E

    (1997) Hedgehog directly controls initiation and propagation of retinal differentiation in the Drosophila eye

    Genes & Development 11:3254–3264.

    https://doi.org/10.1101/gad.11.23.3254

    • Google Scholar
11. 1. Dumstrei K
    2. Nassif C
    3. Abboud G
    4. Aryai A
    5. Aryai A
    6. Hartenstein V

    (1998)

    EGFR signaling is required for the differentiation and maintenance of neural progenitors along the dorsal midline of the Drosophila embryonic head

    Development 125:3417–3426.

    • Google Scholar
12. 1. Flores GV
    2. Duan H
    3. Yan H
    4. Nagaraj R
    5. Fu W
    6. Zou Y
    7. Noll M
    8. Banerjee U

    (2000) Combinatorial signaling in the specification of unique cell fates

    Cell 103:75–85.

    https://doi.org/10.1016/S0092-8674(00)00106-9

    • Google Scholar
13. 1. Fortini ME
    2. Simon MA
    3. Rubin GM

    (1992) Signalling by the Sevenless protein tyrosine kinase is mimicked by Ras1 activation

    Nature 355:559–561.

    https://doi.org/10.1038/355559a0

    • Google Scholar
14. 1. Freeman M

    (1996) Reiterative use of the EGF receptor triggers differentiation of all cell types in the Drosophila eye

    Cell 87:651–660.

    https://doi.org/10.1016/S0092-8674(00)81385-9

    • Google Scholar
15. 1. Gabay L
    2. Scholz H
    3. Golembo M
    4. Klaes A
    5. Shilo BZ
    6. Klämbt C

    (1996)

    EGF receptor signaling induces pointed P1 transcription and inactivates Yan protein in the Drosophila embryonic ventral ectoderm

    Development 122:3355–3362.

    • Google Scholar
16. 1. Goentoro L
    2. Kirschner MW

    (2009) Evidence that fold-change, and not absolute level, of Beta-Catenin dictates Wnt signaling

    Molecular Cell 36:872–884.

    https://doi.org/10.1016/j.molcel.2009.11.017

    • Google Scholar
17. 1. Goentoro L
    2. Shoval O
    3. Kirschner MW
    4. Alon U

    (2009) The incoherent feedforward loop can provide fold-change detection in gene regulation

    Molecular Cell 36:894–899.

    https://doi.org/10.1016/j.molcel.2009.11.018

    • Google Scholar
18. 1. Goldberger AL
    2. Amaral LA
    3. Hausdorff JM
    4. Ivanov PCh
    5. Peng CK
    6. Stanley HE

    (2002) Fractal dynamics in physiology: alterations with disease and aging

    Proceedings of the National Academy of Sciences of the United States of America 99:2466–2472.

    https://doi.org/10.1073/pnas.012579499

    • Google Scholar
19. 1. Golembo M
    2. Raz E
    3. Shilo BZ

    (1996)

    The Drosophila embryonic midline is the site of Spitz processing, and induces activation of the EGF receptor in the ventral ectoderm

    Development 122:3363–3370.

    • Google Scholar
20. 1. Graham TG
    2. Tabei SM
    3. Dinner AR
    4. Rebay I

    (2010) Modeling bistable cell-fate choices in the Drosophila eye: qualitative and quantitative perspectives

    Development 137:2265–2278.

    https://doi.org/10.1242/dev.044826

    • Google Scholar
21. 1. Halfon MS
    2. Carmena A
    3. Gisselbrecht S
    4. Sackerson CM
    5. Jiménez F
    6. Baylies MK
    7. Michelson AM

    (2000) Ras pathway specificity is determined by the integration of multiple signal-activated and tissue-restricted transcription factors

    Cell 103:63–74.

    https://doi.org/10.1016/S0092-8674(00)00105-7

    • Google Scholar
22. 1. Hayashi T
    2. Xu C
    3. Carthew RW

    (2008) Cell-type-specific transcription of prospero is controlled by combinatorial signaling in the Drosophila eye

    Development 135:2787–2796.

    https://doi.org/10.1242/dev.006189

    • Google Scholar
23. 1. Jarman AP
    2. Grell EH
    3. Ackerman L
    4. Jan LY
    5. Jan YN

    (1994) Atonal is the proneural gene for Drosophila photoreceptors

    Nature 369:398–400.

    https://doi.org/10.1038/369398a0

    • Google Scholar
24. 1. Jurgens G
    2. Wieschaus E
    3. Nusslein-Volhard C
    4. Kluding H

    (1984) Mutations affecting the pattern of the larval cuticle in Drosophila melanogaster

    Wilhelm Roux's Archives of Developmental Biology 193:283–295.

    https://doi.org/10.1007/BF00848157

    • Google Scholar
25. 1. Kauffmann RC
    2. Li S
    3. Gallagher PA
    4. Zhang J
    5. Carthew RW

    (1996) Ras1 signaling and transcriptional competence in the R7 cell of Drosophila

    Genes & Development 10:2167–2178.

    https://doi.org/10.1101/gad.10.17.2167

    • Google Scholar
26. 1. Kueh HY
    2. Champhekar A
    3. Champhekhar A
    4. Nutt SL
    5. Elowitz MB
    6. Rothenberg EV

    (2013) Positive feedback between PU.1 and the cell cycle controls myeloid differentiation

    Science 341:670–673.

    https://doi.org/10.1126/science.1240831

    • Google Scholar
27. 1. Kumar JP
    2. Tio M
    3. Hsiung F
    4. Akopyan S
    5. Gabay L
    6. Seger R
    7. Shilo BZ
    8. Moses K

    (1998)

    Dissecting the roles of the Drosophila EGF receptor in eye development and MAP kinase activation

    Development 125:3875–3885.

    • Google Scholar
28. 1. Kumar RM
    2. Cahan P
    3. Shalek AK
    4. Satija R
    5. DaleyKeyser AJ
    6. Li H
    7. Zhang J
    8. Pardee K
    9. Gennert D
    10. Trombetta JJ
    11. Ferrante TC
    12. Regev A
    13. Daley GQ
    14. Collins JJ

    (2014) Deconstructing transcriptional heterogeneity in pluripotent stem cells

    Nature 516:56–61.

    https://doi.org/10.1038/nature13920

    • Google Scholar
29. 1. Lai Z
    2. Rubin GM

    (1992) Negative control of photoreceptor development in Drosophila by the product of the yan gene, an ETS domain protein

    Cell 70:609–620.

    https://doi.org/10.1016/0092-8674(92)90430-K

    • Google Scholar
30. 1. Laslo P
    2. Spooner CJ
    3. Warmflash A
    4. Lancki DW
    5. Lee HJ
    6. Sciammas R
    7. Gantner BN
    8. Dinner AR
    9. Singh H

    (2006) Multilineage transcriptional priming and determination of alternate hematopoietic cell fates

    Cell 126:755–766.

    https://doi.org/10.1016/j.cell.2006.06.052

    • Google Scholar
31. 1. Lebovitz RM
    2. Ready DF

    (1986) Ommatidial development in Drosophila eye disc fragments

    Developmental Biology 117:663–671.

    https://doi.org/10.1016/0012-1606(86)90335-0

    • Google Scholar
32. 1. Legent K
    2. Treisman JE

    (2008) Wingless signaling in Drosophila eye development

    Methods in Molecular Biology 469:141–161.

    https://doi.org/10.1007/978-1-60327-469-2_12

    • Google Scholar
33. 1. Li X
    2. Cassidy JJ
    3. Reinke CA
    4. Fischboeck S
    5. Carthew RW

    (2009) A microRNA imparts robustness against environmental fluctuation during development

    Cell 137:273–282.

    https://doi.org/10.1016/j.cell.2009.01.058

    • Google Scholar
34. 1. Lucchesi JC
    2. Rawls JM

    (1973) Regulation of gene function: a comparison of enzyme activity levels in relation to gene dosage in diploids and triploids of Drosophila melanogaster

    Biochemical Genetics 9:41–51.

    https://doi.org/10.1007/BF00485589

    • Google Scholar
35. 1. MacArthur BD
    2. Lemischka IR

    (2013) Statistical mechanics of pluripotency

    Cell 154:484–489.

    https://doi.org/10.1016/j.cell.2013.07.024

    • Google Scholar
36. 1. Martinez Arias A
    2. Brickman JM

    (2011) Gene expression heterogeneities in embryonic stem cell populations: origin and function

    Current Opinion in Cell Biology 23:650–656.

    https://doi.org/10.1016/j.ceb.2011.09.007

    • Google Scholar
37. 1. Melen GJ
    2. Levy S
    3. Barkai N
    4. Shilo BZ

    (2005) Threshold responses to morphogen gradients by zero-order ultrasensitivity

    Molecular Systems Biology 1:0028.

    https://doi.org/10.1038/msb4100036

    • Google Scholar
38. 1. Morimoto AM
    2. Jordan KC
    3. Tietze K
    4. Britton JS
    5. O'Neill EM
    6. Ruohola-Baker H

    (1996)

    Pointed, an ETS domain transcription factor, negatively regulates the EGF receptor pathway in Drosophila oogenesis

    Development 122:3745–3754.

    • Google Scholar
39. 1. O'Neill EM
    2. Rebay I
    3. Tjian R
    4. Rubin GM

    (1994) The activities of two ETS-related transcription factors required for Drosophila eye development are modulated by the Ras/MAPK pathway

    Cell 78:137–147.

    https://doi.org/10.1016/0092-8674(94)90580-0

    • Google Scholar
40. 1. Ohshiro T
    2. Emori Y
    3. Saigo K

    (2002) Ligand-dependent activation of breathless FGF receptor gene in Drosophila developing trachea

    Mechanisms of Development 114:3–11.

    https://doi.org/10.1016/S0925-4773(02)00042-4

    • Google Scholar
41. 1. Olguin HC
    2. Yang Z
    3. Tapscott SJ
    4. Olwin BB

    (2007) Reciprocal inhibition between Pax7 and muscle regulatory factors modulates myogenic cell fate determination

    The Journal of Cell Biology 177:769–779.

    https://doi.org/10.1083/jcb.200608122

    • Google Scholar
42. 1. Park BO
    2. Ahrends R
    3. Teruel MN

    (2012) Consecutive positive feedback loops create a bistable switch that controls preadipocyte-to-adipocyte conversion

    Cell Reports 2:976–990.

    https://doi.org/10.1016/j.celrep.2012.08.038

    • Google Scholar
43. 1. Qi J
    2. Wang B
    3. Pelaez N
    4. Rebay I
    5. Carthew R
    6. Katsaggelos A
    7. Amaral L

    (2013)

    Drosophila eye nuclei segmentation based on graph cut and convex shape prior. IEEE intl conf signal processing

    978-1-4799-2341-0/13/$ 31:00.

    • Google Scholar
44. 1. Ready DF
    2. Hanson TE
    3. Benzer S

    (1976) Development of the Drosophila retina, a neurocrystalline lattice

    Developmental Biology 53:217–240.

    https://doi.org/10.1016/0012-1606(76)90225-6

    • Google Scholar
45. 1. Rebay I
    2. Rubin GM

    (1995) Yan functions as a general inhibitor of differentiation and is negatively regulated by activation of the Ras1/MAPK pathway

    Cell 81:857–866.

    https://doi.org/10.1016/0092-8674(95)90006-3

    • Google Scholar
46. 1. Rogge R
    2. Green PJ
    3. Urano J
    4. Horn-Saban S
    5. Mlodzik M
    6. Shilo BZ
    7. Hartenstein V
    8. Banerjee U

    (1995)

    The role of yan in mediating the choice between cell division and differentiation

    Development 121:3947–3958.

    • Google Scholar
47. 1. Schaffer AE
    2. Freude KK
    3. Nelson SB
    4. Sander M

    (2010) Nkx6 transcription factors and Ptf1a function as antagonistic lineage determinants in multipotent pancreatic progenitors

    Developmental Cell 18:1022–1029.

    https://doi.org/10.1016/j.devcel.2010.05.015

    • Google Scholar
48. 1. Shilo BZ

    (2014) The regulation and functions of MAPK pathways in Drosophila

    Methods 68:151–159.

    https://doi.org/10.1016/j.ymeth.2014.01.020

    • Google Scholar
49. 1. Surkova S
    2. Kosman D
    3. Kozlov K
    4. Manu
    5. Myasnikova E
    6. Samsonova AA
    7. Spirov A
    8. Vanario-Alonso CE
    9. Samsonova M
    10. Reinitz J
    11. Spirov CE
    12. Vanario-Alonso M
    13. Reinitz J

    (2008) Characterization of the drosophila segment determination morphome

    Developmental Biology 313:844–862.

    https://doi.org/10.1016/j.ydbio.2007.10.037

    • Google Scholar
50. 1. Tomlinson A
    2. Ready DF

    (1987) Neuronal differentiation in the Drosophila ommatidium

    Developmental Biology 120:366–376.

    https://doi.org/10.1016/0012-1606(87)90239-9

    • Google Scholar
51. 1. Tomlinson A

    (1985)

    The cellular dynamics of pattern formation in the eye of Drosophila

    Journal of Embryology and Experimental Morphology 89:.

    • Google Scholar
52. 1. Torres-Padilla ME
    2. Chambers I

    (2014) Transcription factor heterogeneity in pluripotent stem cells: a stochastic advantage

    Development 141:2173–2181.

    https://doi.org/10.1242/dev.102624

    • Google Scholar
53. 1. Van Vactor D
    2. Krantz DE
    3. Reinke R
    4. Zipursky SL

    (1988) Analysis of mutants in chaoptin, a photoreceptor cell-specific glycoprotein in Drosophila, reveals its role in cellular morphogenesis

    Cell 52:281–290.

    https://doi.org/10.1016/0092-8674(88)90517-X

    • Google Scholar
54. 1. Wartlick O
    2. Mumcu P
    3. Kicheva A
    4. Bittig T
    5. Seum C
    6. Jülicher F
    7. González-Gaitán M

    (2011) Dynamics of DPP signaling and proliferation control

    Science 331:1154–1159.

    https://doi.org/10.1126/science.1200037

    • Google Scholar
55. 1. Webber JL
    2. Zhang J
    3. Cote L
    4. Vivekanand P
    5. Ni X
    6. Zhou J
    7. Nègre N
    8. Carthew RW
    9. White KP
    10. Rebay I

    (2013) The relationship between long-range chromatin occupancy and polymerization of the Drosophila ETS family transcriptional repressor Yan

    Genetics 193:633–649.

    https://doi.org/10.1534/genetics.112.146647

    • Google Scholar
56. 1. Wernet MF
    2. Labhart T
    3. Baumann F
    4. Mazzoni EO
    5. Pichaud F
    6. Desplan C

    (2003) Homothorax switches function of Drosophila photoreceptors from color to polarized light sensors

    Cell 115:267–279.

    https://doi.org/10.1016/S0092-8674(03)00848-1

    • Google Scholar
57. Book

    1. Wolff T
    2. Ready D

    (1993)

    Pattern formation in the Drosophila retina

    In: Bate M, Martinez-Arias A, editors. In the Development of Drosophila Melanogaster (Vol. 2).  Plainview, NY: Cold Spring Harbor Press. pp. 1277–1325.

    • Google Scholar
58. 1. Xu C
    2. Kauffmann RC
    3. Zhang J
    4. Kladny S
    5. Carthew RW

    (2000) Overlapping activators and repressors delimit transcriptional response to receptor tyrosine kinase signals in the drosophila eye

    Cell 103:87–97.

    https://doi.org/10.1016/S0092-8674(00)00107-0

    • Google Scholar
59. 1. Yao G
    2. Lee TJ
    3. Mori S
    4. Nevins JR
    5. You L

    (2008) A bistable Rb-E2F switch underlies the restriction point

    Nature Cell Biology 10:476–482.

    https://doi.org/10.1038/ncb1711

    • Google Scholar

Author details

1. Nicolás Peláez

   1. Department of Molecular Biosciences, Northwestern University, Evanston, United States
   2. Department of Chemical and Biological Engineering, Howard Hughes Medical Institute, University Northwestern, Evanston, United States

   Contribution

   NP, Conception and design; Acquisition of data; Analysis and interpretation of data; Drafting or revising the article

   Competing interests

   The authors declare that no competing interests exist.

2. Arnau Gavalda-Miralles

   Department of Chemical and Biological Engineering, Howard Hughes Medical Institute, University Northwestern, Evanston, United States

   Present address

   Department of Chemical Engineering, Universitat Rovira i Virgili, Tarragona, United States

   Contribution

   AGM, Analysis and interpretation of data; Drafting or revising the article

   Competing interests

   The authors declare that no competing interests exist.

3. Bao Wang

   Department Electrical Engineering and Computer Science, Northwestern University, Evanston, United States

   Contribution

   BW, Acquisition of data, Drafting or revising the article

   Competing interests

   The authors declare that no competing interests exist.

4. Heliodoro Tejedor Navarro

   Department of Chemical and Biological Engineering, Howard Hughes Medical Institute, University Northwestern, Evanston, United States

   Contribution

   HTN, Acquisition of data, Drafting or revising the article

   Competing interests

   The authors declare that no competing interests exist.

5. Herman Gudjonson

   James Franck Institute, University of Chicago, Chicago, United States

   Contribution

   HG, Analysis and interpretation of data, Drafting or revising the article

   Competing interests

   The authors declare that no competing interests exist.

6. Ilaria Rebay

   Ben May Department for Cancer Research, University of Chicago, Chicago, United States

   Contribution

   IR, Conception and design, Analysis and interpretation of data, Drafting or revising the article

   Competing interests

   The authors declare that no competing interests exist.

7. Aaron R Dinner

   James Franck Institute, University of Chicago, Chicago, United States

   Contribution

   ARD, Analysis and interpretation of data, Drafting or revising the article

   Competing interests

   The authors declare that no competing interests exist.

8. Aggelos K Katsaggelos

   Department Electrical Engineering and Computer Science, Northwestern University, Evanston, United States

   Contribution

   AKK, Acquisition of data, Drafting or revising the article

   Competing interests

   The authors declare that no competing interests exist.

9. Luís AN Amaral

   1. Department of Chemical and Biological Engineering, Howard Hughes Medical Institute, University Northwestern, Evanston, United States
   2. Department of Physics and Astronomy, Northwestern University, Evanston, United States

   Contribution

   LANA, Conception and design; Analysis and interpretation of data; Drafting or revising the article

   For correspondence

   amaral@northwestern.edu

   Competing interests

   The authors declare that no competing interests exist.

10. Richard W Carthew

    Department of Molecular Biosciences, Northwestern University, Evanston, United States

    Contribution

    RWC, Conception and design, Analysis and interpretation of data, Drafting or revising the article

    For correspondence

    r-carthew@northwestern.edu

    Competing interests

    The authors declare that no competing interests exist.

• 4,438

  views

• 778

  downloads

• 36

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.