Condensation of Ede1 promotes the initiation of endocytosis

  1. Geneva
    Switzerland

Abstract

Clathrin-mediated endocytosis is initiated by a network of weakly interacting proteins through a poorly understood mechanism. Ede1, the yeast homolog of mammalian Eps15, is an early-arriving endocytic protein and a key initiation factor. In the absence of Ede1, most other early endocytic proteins lose their punctate localization and endocytic uptake is decreased. We show that in yeast cells, cytosolic concentration of Ede1 is buffered at a critical level. Excess amounts of Ede1 form large condensates which recruit other endocytic proteins and exhibit properties of phase-separated liquid droplets. We demonstrate that the central region of Ede1, containing a coiled-coil and a prion-like region, is essential for both the condensate formation and the function of Ede1 in endocytosis. The functionality of Ede1 mutants lacking the central region can be partially rescued by an insertion of heterologous prion-like domains. Conversely, fusion of a heterologous lipid-binding domain with the central region of Ede1 can promote clustering into stable plasma membrane domains. We propose that the ability of Ede1 to form condensed networks supports the clustering of early endocytic proteins and promotes the initiation of endocytosis.

Introduction

Clathrin-mediated endocytosis is a process which eukaryotic cells use to produce small transport vesicles from their plasma membrane. These vesicles deliver cargo molecules to the endolysosomal trafficking network where they are sorted for recycling or degradation 24Grant and Donaldson2009. This process is the primary route of internalization of extracellular and surface molecules in eukaryotic cells 31Kirchhausen et al.2014.

Clathrin-mediated endocytosis requires a complex protein machinery to assemble on the plasma membrane in a specific sequence 28Kaksonen et al.200571Sirotkin et al.201078Taylor et al.2011. Endocytosis starts with the arrival of pioneer proteins which select the site and initiate the assembly of the endocytic coat 29Kaksonen and Roux2018. Different pioneer proteins 26Henne et al.201013Cocucci et al.201245Ma et al.2016, lipids 2Antonescu et al.2011 and cargo molecules 41Liu et al.201038Layton et al.2011 have been shown to promote the initiation step (reviewed by 23Godlee and Kaksonen2013). However, the exact mechanism by which the endocytic sites are initiated remains poorly understood.

The molecular mechanisms of clathrin-mediated endocytosis have been studied extensively in yeasts, such as the budding yeast Saccharomyces cerevisiae. The pioneer module in yeast includes highly conserved adaptor proteins which bind membrane and cargo, such as the adaptor protein complex 2 (AP-2 complex), Syp1 (FCHo1/2 in mammals) and Yap1801/2 (AP180). Two conserved scaffold proteins, clathrin and Ede1 (Eps15), also arrive during the early phase. In contrast to the remarkably well-ordered assembly of the membrane-bending phase 63Picco et al.2015, the early proteins lack a specific sequence of recruitment 10Carroll et al.201259Pedersen et al.2020. Curiously, all the genes coding for the earliest-arriving proteins can be deleted without completely blocking endocytosis 9Brach et al.2014. However, the frequency of endocytic events is decreased in such cells, and the ability to regulate cargo recruitment is drastically compromised.

Ede1, a homologue of mammalian Eps15, is one of the key early proteins. It is among the earliest to appear at the nascent endocytic site 10Carroll et al.2012. The deletion of the EDE1 gene reduces the overall membrane uptake by 35 % 22Gagny et al.2000 and decreases the number of productive ndocytic events by 50 % 10Carroll et al.201228Kaksonen et al.2005. Ede1 interacts via its three Eps15 homology (EH) domains and other interaction motifs with proteins such as the AP-2 complex, Hrr25, epsins, Sla2, Syp1, and Yap1801/2 46Maldonado-Báez et al.200864Reider et al.200960Peng et al.2015. Several of these proteins depend on interaction with Ede1 to become enriched at the endocytic sites 10Carroll et al.2012. Ede1 also oligomerizes via its coiled-coil domain, which is required for it to properly localise and function 8Boeke et al.201443Lu and Drubin2017.

In recent years, liquid-liquid phase separation of biomolecules has garnered much attention as a mechanism for assembly of cellular organelles which lack surrounding membranes, such as P granules, nucleoli and stress granules 5Banani et al.2017. Such membrane-less compartments can accelerate reactions, sequester molecules from the cytoplasm or establish spatial organisation 44Lyon et al.2021. The phase separation framework can also apply to sub-micrometre compartments that balance the needs of concentrating select components with allowing a dynamic exchange and rearrangement of molecules. Examples of such compartments include transcriptional superenhancers 67Sabari et al.2018 and a wide array of membrane receptor clusters 11Case et al.2019.

In this work, we show that Ede1 has a propensity to form cellular condensates. We demonstrate that the cytosolic concentration of Ede1 is buffered at a critical concentration. We identify the molecular features driving the condensation of Ede1 and show that they are essential for the normal function of Ede1 during endocytosis. Our findings suggest that Ede1 has a ability form molecular condensates, which promote the initiation and maturation of endocytic sites. The Ede1 condensates exhibit many of the hallmark properties of phase separated liquids.

Results

Ede1 can form dynamic protein condensates

In normal yeast cells, fluorescently tagged Ede1 localizes to endocytic sites at the plasma membrane 34Kukulski et al.2012. However, we discovered previously that under certain experimental conditions Ede1 can also assemble into large condensates 8Boeke et al.2014. These condensates were seen in cells that either overexpressed Ede1, or expressed Ede1 at normal levels, but lacked three early endocytic adaptors. Although these condensates are abnormal structures that have not been observed in wild-type cells, we reasoned that studying them in more detail might provide insights into the mechanism by which Ede1 promotes the assembly of the early endocytic proteins.

To visualize the condensates, we expressed Ede1 fused to enhanced green fluorescent protein (EGFP) from its endogenous locus in haploid wild-type cells, or in cells from which three endocytosis-related genes were deleted (yap1801Δ yap1802Δ apl3Δ, called 3×ΔEA for short). These genes code for early-arriving endocytic adaptor proteins Yap1801, Yap1802, and the α-subunit of the AP-2 complex. Alternatively, we overexpressed EGFP-Ede1 from its endogenous locus under the control of a strong heterologous promoter.

We observed that part of the cellular Ede1-EGFP in the mutant strains localized into condensates that were much brighter than the normal endocytic sites (Figure 1A). The condensates in the overexpression strain were larger and brighter than those in the 3×ΔEA cells (Figure 1B). The condensates usually associated with the plasma membrane, but were also observed away from it (Figure 1A). In contrast, normal endocytic sites are always associated with the plasma membrane. The condensates in the Ede1 overexpression strain were often large enough that their shape was resolvable (Figure 1B). They appeared circular in surface view, and as dome-like structures limited by the plasma membrane in the side view. The Ede1 condensates were remarkably long lived and we have observed individual condensates for up to 1 hour (Figure 1—video 1). This stands in contrast with normal endocytic sites, where Ede1-EGFP typically persists for 1-2 minutes 75Stimpson et al.2009. Despite their stability, some of the condensates appeared to undergo dynamic fission and fusion events, suggesting that they are not solid aggregates (Figure 1C).

Excess cytosolic Ede1 assembles into condensates in vivo.

(A) Representative images of yeast cells expressing Ede1-EGFP in wild-type (wt) and 3×ΔEA genetic backgrounds, or overexpressing EGFP-Ede1 under the control of the ADH1 promoter. Mutant cell micrographs are shown using the same display range as the wt (top) or their full display range (bottom). Two cells are shown for 3×ΔEA background to display the membrane-associated and cytoplasmic localizations of Ede1 condensates. Scale bars: 2 μm. (B) Representative images of Ede1-EGFP at endocytic sites in wt background, Ede1-EGFP condensates in 3×ΔEA cells, and EGFP-Ede1 overexpression-induced condensates (OE). OE condensates are shown in two different orientations. Each frame is 1.5 μm × 1.5 μm; dotted white line represents the approximate position of the plasma membrane. (C) Two time series of Ede1-EGFP condensates undergoing apparent fusion (top) and fission (bottom) events. Scale bar: 1 μm. (D) Representative images of diploid cells homozygous for the 3×ΔEA background, each expressing Ede1-EGFP and differing in the second Ede1 locus: EDE1-EGFP, EDE1, or ede1Δ. Scale bars: 2 μm. (E) The fraction of cells containing condensates in each strain from panel D. Bars and whiskers show mean ± SD of three independent experiments. A range of 40–70 cells were analyzed per data point.

#' @width 28
#' @height 20
library(tidyverse)
library(Hmisc) # mean_sdl and binomial confs
library(broom)
library(ggbeeswarm)
library(ggsignif)
library(lmerTest)
library(emmeans)
library(rstatix)
library(multcompView)

strains <- read_csv("data/fig1/raw/strain_lut.csv", 
                    col_types = cols(second_ede1  = col_factor(NULL))) # order factor levels in csv file
#knitr::kable(strains)

# initialize a data frame with counts and get information about experiment from file names
files <- list.files(path = 'data/fig1/raw/', pattern = 'dataset',
                    full.names = TRUE)

counts <- map(files, read_csv,
              col_names = c('cell', 'condensates', 'dataset'),
              cols(dataset = col_factor(NULL))) %>%
          bind_rows() %>%
          separate(cell, c('date', 'strain', 'stack', 'cell',
                   'channel'), '_') %>%
          left_join(strains, by = 'strain') %>%
          select(strain, stack, second_ede1, condensates, dataset)

proc <- counts %>%
  group_by(second_ede1, dataset) %>%
  summarise(x = sum(condensates), n = n(),
            p = binconf(x, n, method = 'wilson')[1],
            p_min = binconf(x, n, method = 'wilson')[2],
            p_max = binconf(x, n, method = 'wilson')[3])

proc_stats <- proc %>%
  select(second_ede1, dataset, p) %>%
  group_by(second_ede1) %>%
  summarise(mean = mean(p), sd = sd(p), n = n())

#knitr::kable(proc)

proc_stats <- proc %>%
  select(second_ede1, dataset, p) %>%
  group_by(second_ede1) %>%
  summarise(mean = mean(p), sd = sd(p), n = n())

#knitr::kable(proc_stats)

theme_frap <- function(base_size = 11, base_family = "",
                       base_line_size = base_size / 22,
                       base_rect_size = base_size / 22) {
  # minimal theme with border
  # based on theme_linedraw without the grid lines
  # also trying to remove all backgrounds but plot_background won't set to blank
  theme_linedraw(
    base_size = base_size,
    base_family = base_family,
    base_line_size = base_line_size,
    base_rect_size = base_rect_size
  ) %+replace%
    theme(
      # no grid and no backgrounds if I can help it
      plot.background = element_blank(),
      panel.background = element_blank(),
      panel.grid = element_blank(),
      plot.margin = margin(0, 0, 0, 0),
      complete = TRUE
    )
}

myplot <- proc %>% ggplot(aes(x = second_ede1, y = p)) + 
  geom_quasirandom(aes(shape = dataset), fill = 'gray80',
              size = 2, width = 0.3, show.legend = FALSE)+
  stat_summary(fun.y = mean, geom = "errorbar", 
               aes(ymax = ..y.., ymin = ..y..),
               width = 0.5, size = 1)+
  stat_summary(fun.data = 'mean_sdl', geom = 'errorbar',
               fun.args = list(mult = 1), width = 0.2)+
  #scale_color_brewer(palette = 'Set2')+
  scale_shape_manual(values = c(21:25))+
  theme_frap(base_size = 18, base_family = "Helvetica")+
  scale_y_continuous(labels = scales::label_percent(suffix = ''))+
  theme(legend.position = c(0.95, 0.95),
        legend.justification = c("right", "top"),
        legend.background = element_rect(fill = NULL, colour = 'black',
                                         size = 0.5))+
  labs(title = NULL,
       y = 'Cells with condensates (%)',
       x = 'Second Ede1 allele')
myplot

Legend

Representative 3×ΔEA cells expressing Ede1 tagged with mNeonGreen, msGFP2, or mCherry.

All three adaptors absent from 3×ΔEA cells interact with Ede1, as well as membrane lipids and protein cargo. As Ede1 does not have known membrane-binding activity, it is likely that the cytosolic pool of Ede1 is increased both in the overexpression and the 3×ΔEA backgrounds, and the excess protein assembles into the condensates. To test whether the formation of condensates in 3×ΔEA background depends on Ede1 concentration, we generated diploid cells homozygous for the three adaptor deletions. We expressed Ede1-EGFP in these cells either from both EDE1 alleles (EDE1-EGFP/EDE1-EGFP), or from one allele in combination with untagged EDE1 (EDE1-EGFP/EDE1) or a deletion of EDE1 (EDE1-EGFP/ede1Δ). The condensates formed in both strains expressing two alleles of EDE1, but not in the strain where only one EDE1 allele was present (Figure 1D). This result suggests that the condensate assembly depends on Ede1 concentration.

The EGFP tag used in these experiments has a weak tendency to dimerize and can induce protein clustering in some in vivo contexts 14Costantini et al.2012. We expressed Ede1 tagged with monomeric fluorescent proteins mNeonGreen, msGFP2, and mCherry in 3×ΔEA cells and observed the same phenotype as with EGFP-tagged Ede1 (Figure 1—figure supplement 1). The observed phenotype is therefore unlikely to have been caused by the choice of the fluorescent tag.

It must also be noted that there is a large difference in brightness between genuine endocytic sites and condensates in mutant cells. Therefore, we chose to saturate the display of those images which show both classes of objects simultaneously, in order for the cells and endocytic sites to remain visible. This difference is demonstrated in panel A of Figure 1.

Ede1 condensates exhibit liquid-like properties

Because of their spherical shapes, concentration dependence, and dynamic behaviors, we hypothesized that the Ede1 condensates might be phase-separated liquid droplets. To test this idea, we first performed fluorescence recovery after photobleaching (FRAP) experiments on Ede1-EGFP condensates in the 3×ΔEA background. After photobleaching, the condensates rapidly recovered most of their fluorescence (Figure 2A). The recovery half-time of a single-exponential FRAP model fitted to an average of 36 events was 22 s, and the mobile fraction was 63%.

Ede1 structures exchange molecules with the cytoplasm and respond to temperature changes.

(A, B) Fluorescence recovery after photobleaching (FRAP) of Ede1-EGFP condensates in 3×ΔEA cells (A) and endocytic sites in normal cells (B). Plots show mean fluorescence recovery ± SD; n = 36 across four independent experiments (panel A) and n = 14 across three independent experiments (panel B). Representative time series are shown below each plot. Each frame is 1 μm × 1 μm. (C) Time series of a partial bleaching of a condensate in a cell overexpressing EGFP-Ede1. A perceptually uniform color lookup table has been applied to highlight the changes in intensity. Each frame is 1.5 μm × 1.5 μm.(D) Average fluorescence recovery (n = 14) after partial bleaching of condensates in cells overexpressing EGFP-Ede1, as in panel (C). (E) Representative cells after 5-min treatment with indicated concentrations of 1,6-hexanediol. Maximum Z-projections. (F) Ede1-EGFP was imaged in 3×ΔEA cells at different temperatures. Cells were grown and imaged at 24°C. The temperature was raised to 37 and 42°C and returned to 24°C for the indicated amounts of time. Maximum Z-projections. All scale bars: 2 μm.

We then examined the recovery of Ede1 in normal endocytic sites, photobleaching them during total internal reflection fluorescence (TIRF) imaging. We found that Ede1-EGFP also turns over fast at endocytic sites (half-time of 7.8 s and mobile fraction of 91%). The turnover at the endocytic sites was faster and the mobile fraction higher than in the condensates.

We also performed experiments where we bleached the Ede1 condensates only partially, in order to visualize diffusion of fluorescent molecules within the condensates. For these experiments we used cells overexpressing Ede1-EGFP, in which the condensates are larger than in the 3×ΔEA cells. When a subregion of an endocytic condensate was bleached, the fluorescence in the bleached and unbleached regions equalized within several seconds (Figure 2C and D; Figure 2—video 1).

The formation of phase-separated condensates depends on protein concentration and can be affected by environmental factors such as temperature 50Molliex et al.201556Nott et al.201520Franzmann et al.2018. When the 3×ΔEA cells cultured at 24°C were incubated at 37°C for 5 min, the condensates dissolved while the endocytic sites persisted (Figure 2D). When the temperature was raised to 42°C, Ede1 signal became entirely diffuse. This effect was reversible for both the endocytic sites, which reformed after several minutes, and the condensates, which reappeared within 30 min after return to 24°C.

1,6-Hexanediol is an aliphatic alcohol that disrupts weak protein-protein interactions 57Patel et al.2007 and is used to distinguish between solid and liquid protein aggregates 33Kroschwald et al.2015. Both Ede1 condensates and endocytic patches rapidly disappeared in 3×ΔEA cells upon 1,6-hexanediol treatment (Figure 2E). Interestingly, endocytic patches only fully dissolved at higher 1,6-hexanediol concentrations than Ede1 condensates.

Taken together our results show that Ede1 both in the condensates and at the endocytic sites behaves in a highly dynamic, liquid-like manner.

Next, we sought to determine the concentration of Ede1 in the 3×ΔEA and Ede1 overexpression cells relative to the wild-type cells. We used spinning-disk confocal microscopy to limit the influence of out-of-focus condensates on the quantification of cytosolic intensity. We quantified the mean pixel intensity in entire cell volumes and in small cytosolic regions of cell cross-sections (Figure 3A). The total cellular intensity was on average 122 and 321% of the wild-type in 3×ΔEA and overexpression cells, respectively. The cytosolic intensity was nevertheless uniform across all three strains and we could detect no statistically significant differences (p=0.45 in F-test).

# Clean raw data
cyto_full <- dir('data/fig3_concentrations/cyto', pattern = 'dataset', full.names = T) %>%
  map(dir, pattern = '.csv', full.names = T) %>%
  map(
    ~ set_names(.x,
                nm = (basename(.x) %>% tools::file_path_sans_ext()))) %>%
  map(
    ~ map_df(.x, read_csv,
         col_types = cols(), # suppress the column spec message
         col_names = c('frame', 'intensity'), # set up names
         skip = 1, # skip the ImageJ header
         .id = "filename")) %>% # id observations by filename
  imap_dfr(~.x, .id = 'dataset') %>%
  select(-frame) %>%
  separate(filename, c('date', 'strain', 'image', 'channel', 'processing', 'window'), '_')

cell_full <- dir('data/fig3_concentrations/cell', pattern = 'dataset', full.names = T) %>%
  map(read_csv, col_types = cols()) %>%
  imap_dfr(~.x, .id = 'dataset') %>%
  separate(cell, c('date', 'strain', 'image', 'window', 'channel', 'processing'), '_') %>%
  rename(int_mean = mean, int_median = median, int_sd = sd)

cyto <- cyto_full %>%
  select(date, strain, image, window, intensity, dataset) %>%
  tibble(., localization = 'Cytosolic')

cell <- cell_full %>%
  rename(intensity = int_mean) %>%
  select(date, strain, image, window, intensity, dataset) %>%
  tibble(., localization = 'Cellular')

strains <- read_csv("data/fig3_concentrations/strain_lut.csv", 
                    col_types = cols(nickname = col_factor(NULL)))

intensities <- bind_rows(cyto, cell) %>%
  left_join(strains, by = 'strain') %>%
  mutate(strain = nickname,
         localization = recode_factor(localization,
                                      'Cellular' = 'Cellular',
                                      'Cytosolic' = 'Cytosolic'),
         .keep = 'unused') %>%
  select(-mCherry)

# Custom ggplot2 theme
# --------------------
theme_clean <- function(base_size = 11, base_family = "",
                        base_line_size = base_size / 22,
                        base_rect_size = base_size / 22) {
  theme_linedraw(
    base_size = base_size,
    base_family = base_family,
    base_line_size = base_line_size,
    base_rect_size = base_rect_size
  ) %+replace%
    theme(
      # no grid and no backgrounds if I can help it
      legend.background =  element_blank(),
      panel.background = element_blank(),
      panel.grid = element_blank(),
      plot.background = element_blank(),
      plot.margin = margin(0, 0, 0, 0),
      complete = TRUE
    )
}

theme_set(theme_clean(base_size = 14, base_family = "Helvetica"))

# generate normalized intensities (within each dataset)

# pull out scaling factor
wt_means <- intensities %>%
  filter(strain == 'wt') %>%
  group_by(dataset, localization) %>%
  summarise(wt_mean = mean(intensity), .groups = 'drop')

# join with the dataframe, generate new variable
intensities <- left_join(intensities, wt_means,
                         by = c('dataset', 'localization')) %>%
  mutate(normalized = intensity / wt_mean)

# generate summary statistics
intensity_stats <- intensities %>%
  group_by(localization, strain, dataset) %>%
  summarise(n = n(),
            across(c(intensity, normalized),
                   list(mean = mean, sd = sd, median = median,
                        se = ~ sd(.x) / sqrt(n())
                        )), .groups = 'drop')

#intensity_stats %>%
#  kable()

mean_ci <- intensity_stats %>%
  group_by(localization, strain)%>%
  summarise(round(mean_cl_normal(normalized_mean), 2)) %>%
  rename(mean = y, lower = ymin, upper = ymax)

#kable(mean_ci)

plot_base <- ggplot(intensity_stats,
                   aes(x = strain, y = intensity_mean)) +
  labs(title = NULL, x = 'Strain', y = 'Intensity (a. u.)')+
  scale_y_continuous(breaks = scales::breaks_extended(6))+
  scale_shape_manual(values = c(21:25)) +
  scale_color_brewer(palette = 'Set2') +
  scale_fill_brewer(palette = 'Set2') +
  facet_wrap(. ~localization,
             scales = 'free')


plot_scatter <- plot_base + 
  geom_quasirandom(inherit.aes = F, data = intensities,
                   aes(x = strain, y = intensity,
                       shape = dataset),
                   colour = 'gray75', #width = 0.3,
                   size = 0.8, #varwidth = TRUE,
                   show.legend = F
                   )

plot_violin <- plot_base + 
  geom_violin(inherit.aes = F, data = intensities,
              aes(x = strain, y = intensity),
              scale = 'width',
              colour = 'gray75', fill = 'transparent'
              )

plot_super <- plot_scatter +
  geom_quasirandom(aes(shape = dataset, fill = dataset),
                   show.legend = F, width = 0.3, size = 2) +
  stat_summary(fun = mean, geom = "crossbar", 
               width = 0.5, fatten = 1) +
  stat_summary(fun.data = 'mean_sdl', geom = 'errorbar',
               fun.args = list(mult = 1), width = 0.2)
plot_super

Condensate formation limits the concentration of Ede1 in the cytoplasm.

Fluorescence intensity was measured in wild-type (wt) and 3×ΔEA cells expressing Ede1-EGFP, or cells overexpressing (oe) EGFP-Ede1 from the ADH1 promoter. Gray points represent mean pixel intensities of entire cell volumes (Cellular), or small regions manually selected from cell cross-sections (Cytosolic). Large points: mean values from independent replicates; central line and whiskers: mean ± SD of replicate means. Pairwise comparisons based on a linear mixed model (n.s., not significant; ***, p<0.001).

Two yeast strains were imaged for 8 hr after change of carbon source from glucose to galactose (see Materials and methods).

Cells express EGFP-Ede1 under the control of GALS promoter (green channel only), or Ede1-EGFP and Sla1-mCherry expressed from the endogenous loci (green and magenta, respectively). Scale bars: 2 μm.

intensities3c <- list.files("data/fig3_gal-induction/imagej/", pattern = '.csv' , full.names = TRUE ) %>%
  set_names(nm = (basename(.) %>% tools::file_path_sans_ext())) %>% # names for the map output, wthout foder or extension
  map_df(read_csv,
         col_types = cols(), # this doesn't do anything but suppresses the column spec message
         col_names = c('frame', 'intensity'), # set up names
         skip = 1, # skip the ImageJ header
         .id = "filename") %>% # id observations by filename
  separate(filename, c('date', 'strain', 'movie', 'channel', 'processing', 'position', 'signal'), '_') %>%
  mutate(time = frame * 5 + 10) %>% # imaging started ~15 minutes after galactose addition
  select(signal, movie, position, time, intensity)

mean_intensities <- intensities3c %>%
  group_by(time, signal) %>%
  # signal variable combines strain and localization for now
  # generate statistics for each signal at each timepoint
  summarise(mean = mean(intensity), sd = sd(intensity), n = n(), se = sd(intensity) / sqrt(n())) %>%
  # generate wide format data: each statistic gets its own column for each signal
  pivot_wider(names_from = signal, values_from = c(mean, sd, se, n)) %>%
  # use the new columns to generate background-subtracted calues (x_corr)
  # and the ratio (corrected oe divided by wt for each localization)
  # and propagate the errors
  mutate(corr_wt.cyto = mean_wt.cyto - mean_bg,
         corr_wt.total = mean_wt.total - mean_bg,
         corr_oe.cyto = mean_oe.cyto - mean_bg,
         corr_oe.total = mean_oe.total - mean_bg,
         corr.sd_wt.cyto = sqrt(sd_wt.cyto ^ 2 + sd_bg ^ 2),
         corr.sd_wt.total = sqrt(sd_wt.total ^ 2 + sd_bg ^ 2),
         corr.sd_oe.cyto = sqrt(sd_oe.cyto ^ 2 + sd_bg ^ 2),
         corr.sd_oe.total = sqrt(sd_oe.total ^ 2 + sd_bg ^ 2),
         corr.se_wt.cyto = sqrt(se_wt.cyto ^ 2 + se_bg ^ 2),
         corr.se_wt.total = sqrt(se_wt.total ^ 2 + sd_bg ^ 2),
         corr.se_oe.cyto = sqrt(se_oe.cyto ^ 2 + se_bg ^ 2),
         corr.se_oe.total = sqrt(se_oe.total ^ 2 + sd_bg ^ 2),
         corr_ratio.cyto = corr_oe.cyto / corr_wt.cyto,
         corr.sd_ratio.cyto = corr_ratio.cyto* sqrt(
           (corr.sd_oe.cyto / corr_oe.cyto) ^ 2 + (corr.sd_wt.cyto / corr_wt.cyto) ^ 2),
         corr.se_ratio.cyto= corr_ratio.cyto* sqrt(
           (corr.se_oe.cyto / corr_oe.cyto) ^ 2 + (corr.se_wt.cyto / corr_wt.cyto) ^ 2), 
         corr_ratio.total = corr_oe.total / corr_wt.total,
         corr.sd_ratio.total = corr_ratio.total * sqrt(
           (corr.sd_oe.total / corr_oe.total) ^ 2 + (corr.sd_wt.total / corr_wt.total) ^ 2),
         corr.se_ratio.total = corr_ratio.total * sqrt(
           (corr.se_oe.total / corr_oe.total) ^ 2 + (corr.se_wt.total / corr_wt.total) ^ 2) 
  )

# go back to long format data: it's near unreadable, but easier to plot
mean_intensities_long <- mean_intensities %>%
  pivot_longer(!time, #process everything except for time
               names_to = c('.value', 'signal'), # new variables
               names_sep = '_') %>% # generate values for new variables
  separate(signal, c('strain', 'localization'), sep = '[.]') %>%
  mutate(localization = recode_factor(localization,
                               'total' = 'Cellular',
                               'cyto' = 'Cytosolic',
                               ))

raw_plot <- ggplot(intensities3c,
                   aes(x = time, y = intensity)) +
  scale_color_brewer(palette = 'Set2')+
  scale_fill_brewer(palette = 'Set2')+
  stat_summary(aes(colour = signal),
               fun.y = mean, geom = 'line', na.rm = T)+
  stat_summary(aes(fill = signal),
               fun.data = 'mean_cl_normal',
               #fun.args = list(mult = 1), # mult = how many deviations
               geom = 'ribbon', alpha = 0.3, na.rm = T)
#raw_plot


base_plot <- ggplot(mean_intensities)

corr_plot <- mean_intensities_long %>%
  filter(strain == 'wt' | strain == 'oe' ) %>%
  mutate(time = time / 60) %>%
    ggplot(aes(x = time)) +
    geom_line(aes(y = corr,
                  colour = strain))+
    geom_ribbon(aes(ymin = corr - 2 * corr.se,
                    ymax = corr + 2 * corr.se,
                    fill = strain),
                alpha = 0.3)+
    scale_color_brewer(palette = 'Set2')+
    scale_fill_brewer(palette = 'Set2')+
    facet_grid(rows = vars(localization),
               scales = 'free',
               )+
    scale_x_continuous(breaks = seq(0, 8, 2))+
    labs(x="Time (h)", y = "Mean fluorescence intensity (a.u.)")

corr_plot

Fluorescence intensity during the expression induction was measured in regions representing entire cells (2 and 4 in panel B) and their cytoplasm (1 and 3).

Mean intensity is shown for endogenously expressed (wt) or overexpressed (oe) Ede1 after background subtraction, ±2 × SEM (n = 40 cells for each strain).

We also induced overexpression of EGFP-tagged Ede1 using the weakened galactokinase promoter GALS 52Mumberg et al.1994. We followed the changes in fluorescence intensity over 8 hr after switching carbon source from glucose to galactose (Figure 3B and C; Figure 3—video 1). The cells, in which overexpression was induced after glucose repression, initially showed no Ede1 sites at the membrane. After several hours, Ede1 appeared as transient endocytic sites, and later formed large and stable condensates. The intensity of the cytosolic regions in the overexpression cells never surpassed the cytosolic intensity of wild-type cells, even as the total intensity of the mutant reached approximately 150% of wild-type intensity by the end of the experiment.

These experiments suggest that the cytosolic concentration of Ede1 is buffered by the formation of the condensates. Moreover, the cytoplasmic fluorescence intensity in wild-type cells is already at the limit observed during overexpression. This signifies that the total concentration in wild-type cells is above the critical concentration required for phase separation.

Ede1 condensates recruit other endocytic proteins

We then imaged double-tagged strains to test whether other endocytic proteins colocalize with Ede1 condensates in the 3×ΔEA background (Figure 4A and B). We found that the condensates contain multiple early (Syp1, Ent1, Sla2) and late (End3, Pan1, Sla1) coat proteins, as well as the actin nucleation-promoting factor Las17, known to physically interact with the End3/Pan1/Sla1 complex 77Sun et al.201519Feliciano and Di Pietro2012.

Endocytic condensates recruit many proteins.

(A) Images of representative cells expressing Ede1-EGFP and indicated endocytic proteins tagged with mCherry in 3×ΔEA background. (C) Montage from timelapse imaging of Ede1-EGFP and Abp1-mCherry during apparent fission of an Ede1 condensate.

strains <- read_csv("data/fig4/raw/strain_lut.csv", 
                    col_types = cols(target  = col_factor(NULL))) # order factor levels per csv file
#knitr::kable(strains)

files <- list.files(path = 'data/fig4/raw/', pattern = 'dataset',
                    full.names = TRUE)
counts <- map(files, read_csv,
              col_names = c('cell', 'colocalizes', 'dataset'),
              cols(dataset = col_factor(NULL))) %>%
          bind_rows() %>% # combine datasets
          separate(cell, c('date', 'strain', 'stack', 
                   'channel', 'processing', 'window'), '_') %>%
          left_join(strains, by = 'strain') %>%
          select(strain, stack, window, target, colocalizes, dataset)

proc <- counts %>%
  group_by(target, dataset) %>%
  summarise(x = sum(colocalizes), n = n(),
            p = binconf(x, n, method = 'wilson')[1],
            p_min = binconf(x, n, method = 'wilson')[2],
            p_max = binconf(x, n, method = 'wilson')[3])

#knitr::kable(proc)


proc_stats <- proc %>%
  select(target, dataset, p) %>%
  group_by(target) %>%
  summarise(mean = mean(p), sd = sd(p), n = n())

#knitr::kable(proc_stats)

theme_frap <- function(base_size = 11, base_family = "",
                       base_line_size = base_size / 22,
                       base_rect_size = base_size / 22) {
  # minimal theme with border
  # based on theme_linedraw without the grid lines
  # also trying to remove all backgrounds but plot_background won't set to blank
  theme_linedraw(
    base_size = base_size,
    base_family = base_family,
    base_line_size = base_line_size,
    base_rect_size = base_rect_size
  ) %+replace%
    theme(
      # no grid and no backgrounds if I can help it
      plot.background = element_blank(),
      panel.background = element_blank(),
      panel.grid = element_blank(),
      plot.margin = margin(0, 0, 0, 0),
      complete = TRUE
    )
}

fig4b <- proc %>% ggplot(aes(x = target, y = p)) + 
  geom_quasirandom(aes(shape = dataset), fill = 'gray80',
              size = 2, width = 0.3, show.legend = FALSE)+
  stat_summary(fun.y = mean, geom = "errorbar", 
               aes(ymax = ..y.., ymin = ..y..),
               width = 0.5, size = 1)+
  stat_summary(fun.data = 'mean_sdl', geom = 'errorbar',
               fun.args = list(mult = 1), width = 0.2)+
  #scale_color_brewer(palette = 'Set2')+
  scale_shape_manual(values = c(21:25))+ # points with borders
  theme_frap(base_size = 18, base_family = "Helvetica")+
  scale_y_continuous(labels = scales::label_percent(suffix = ''))+
  theme(legend.position = c(0.95, 0.95),
        legend.justification = c("right", "top"),
        legend.background = element_rect(fill = NULL, colour = 'black',
                                         size = 0.5))+
  labs(title = NULL,
       y = 'Fraction of colocalizing condensates (%)',
       x = 'Target protein')
fig4b

Fraction of Ede1-EGFP condensates that colocalized with mCherry puncta in each strain from panel A.

Bars and whiskers show mean ± SD of three independent experiments. A range of 36–98 cells were analyzed per data point.

Ede1 condensates do not colocalize with membranes stained by FM4-64.

Representative cells expressing Ede1-EGFP in the 3×ΔEA background, or overexpressing EGFP-Ede1 under the control of the ADH1 promoter after a 60-min staining with FM4-64. Scale bar 2 μm.

Three proteins—Myo5, Abp1, and Rvs167—whose arrival overlaps with actin polymerization at the endocytic sites 76Sun et al.2006 localized to a minority of the condensates (42, 20, and 10%, respectively). We further examined the interaction of condensates and Abp1 using timelapse imaging and found that Abp1-mCherry patches assembled on the condensates transiently (Figure 4C), which explains the partial colocalization. Abp1-mCherry was also present whenever the apparent fission of condensates occurred (Figure 4C, Figure 4—video 1). The appearance of Abp1 during condensate fission suggests that fission could be caused by a force exerted by actin filaments. Some of the proteins we see in the condensates, such as Las17, could potentially trigger actin polymerization.

Costaining with FM4-64 showed that endocytic protein condensates do not contain a significant membrane fraction (Figure 4—figure supplement 1).

Overall, the colocalization experiments showed that the endocytic condensates are complex, and that the proteins contained within them are at least partially functional as they can recruit their interaction partners and are associated with cycles of assembly and disassembly of actin.

Ede1 central region is necessary and sufficient for phase separation

Ede1 is a 1381 amino acid long, multidomain protein (Figure 5A). Its N-terminal region contains three Eps15-homology (EH) domains that interact with asparagine-proline-phenylalanine motifs found on endocytic adaptors such as Ent1/2, Yap1801/2, and Sla2 46Maldonado-Báez et al.2008. Such repeats of domains interacting with linear motifs are known to promote liquid-liquid phase separation 40Li et al.20126Banjade and Rosen2014. The EH domains are followed by a proline-rich region and a coiled-coil domain 64Reider et al.200943Lu and Drubin2017. The C-terminal half of Ede1 contains a Syp1-interacting region 64Reider et al.2009 and a ubiquitin-associated (UBA) domain.

The central region of Ede1 is necessary and sufficient for condensate formation.

(A) The domain structure of Ede1 and prediction of disordered and prion-like regions. EH, Eps15-homology domain; UBA, ubiquitin-associated domain. Domains are drawn to scale according to UniProt entry P34216, and numbers above mark domain boundaries used in our constructs. The top plot represents IUPred2a disorder prediction score, with the shaded areas predicted to be disordered by MobiDB-lite consensus method. Prion-like domain (PLD) prediction score was calculated using the PLAAC software. (B) Representative cells expressing full-length (FL) Ede1 and its truncation mutants in wild-type and 3×ΔEA backgrounds. All constructs are C-terminally tagged with EGFP. Maximum intensity projections of 3D volumes, scale bars: 2 μm.

load('data/fig5_ede1-truncations/truncation_data.RData')

proc <- counts %>%
  group_by(bg, Ede1, dataset) %>%
  summarise(x = sum(nuages), n = n(),
            p = binconf(x, n, method = 'wilson')[1],
            p_min = binconf(x, n, method = 'wilson')[2],
            p_max = binconf(x, n, method = 'wilson')[3])

#knitr::kable(proc)

proc_stats <- proc %>%
  select(bg, Ede1, dataset, p) %>%
  group_by(bg, Ede1) %>%
  summarise(mean = mean(p), sd = sd(p), n = n())

#knitr::kable(proc_stats)

fig5c <- ggplot(proc, aes(x = Ede1, y = p, group = bg, fill = bg)) + 

  stat_summary(fun.y = mean, geom = "bar",
               position = position_dodge(),
               width = 0.5, colour = 'black')+
  stat_summary(fun.data = 'mean_sdl', geom = 'errorbar',
               position = position_dodge(width = 0.5),
               fun.args = list(mult = 1), width = 0.3)+
  geom_point(aes(shape = dataset),
             position = position_jitterdodge(jitter.width = 0.4, 
                                             dodge.width = 0.5),
             size = 2, show.legend = FALSE)+
  scale_color_brewer(palette = 'Set2')+
  scale_fill_grey(start = 1, end = 0.5)+
  scale_shape_manual(values = c(21:25))+
  theme_frap(base_size = 16, base_family = "Helvetica")+
  scale_y_continuous(labels = scales::label_percent(suffix = ''))+
  theme(legend.position = c(0.5, 0.95),
        legend.justification = c("center", "top"),
        #legend.box.just = "right",
        legend.margin = margin(0, 6, 6, 6),
        legend.background = element_rect(fill = NULL, colour = 'black',
                                         size = 0.5),
        legend.title = element_blank())+
  labs(title = NULL,
       y = 'Cells with condensates (%)')

print(plot)

fig5c

The fraction of cells containing condensates in each strain from panel B.

Bars and whiskers show mean ± SD of three independent experiments. A range of 37–78 cells were analyzed per data point.

Levels of Ede1 mutants assessed by western blotting.

(a) Representative western blot of EGFP-tagged Ede1 truncation constructs expressed in wild-type (wt) and 3×ΔEA backgrounds.

bands <- read_csv('data/fig5_blots/20200814_quant.csv',
                  col_types = cols (dataset = col_factor(NULL),
                        Ede1 = col_factor(NULL),
                        background = col_factor(NULL)))

plot <- ggplot(bands, aes(x = Ede1, y = ratio_norm,
                          group = background, fill = background))+
    theme_frap(base_size = 14, base_family = "Helvetica")+
    scale_y_continuous(labels = scales::label_percent(suffix = ''))+
    labs(title = NULL,
       y = 'Normalized band intensity (% of wt)')

plot + 
  stat_summary(fun.y = mean, geom = "bar",
               position = position_dodge(),
               width = 0.5, colour = 'black')+
  stat_summary(fun.data = 'mean_sdl', geom = 'errorbar',
               position = position_dodge(width = 0.5),
               fun.args = list(mult = 1), width = 0.3)+
  geom_point(aes(shape = dataset),
             position = position_jitterdodge(jitter.width = 0.4, 
                                             dodge.width = 0.5),
             size = 2, show.legend = FALSE)+
  scale_fill_grey(start = 1, end = 0.5)+
  scale_shape_manual(values = c(21:25))+
  theme(legend.position = c(0.95, 0.95),
        legend.justification = c("right", "top"),
        #legend.box.just = "right",
        legend.margin = margin(0, 6, 6, 6),
        legend.background = element_rect(fill = NULL, colour = 'black',
                                         size = 0.5),
        legend.title = element_blank())

plot + 
  stat_summary(fun.y = mean, geom = "errorbar",
               aes(ymax = ..y.., ymin = ..y..),
               width = 0.5, size = 1)+
  stat_summary(fun.data = 'mean_sdl', geom = 'errorbar',
               fun.args = list(mult = 1), width = 0.2)+
  geom_point(aes(shape = dataset, fill = dataset),
             position = position_jitter(width = 0.1),
             size = 2, show.legend = FALSE)+
  scale_fill_grey(start = 1, end = 0.5)+
  scale_shape_manual(values = c(21:25))+
  theme(axis.text.x = element_text(angle = -45,
                                   vjust = 0.1, hjust = 0.1))+
  facet_grid(rows = vars(background))

Band intensity relative to full-length (FL) Ede1-EGFP in wt cells was calculated for three independent experiments.

Hog1 was used as a loading control. Bars and whiskers show mean ± SD.

We noticed that the proline-rich region contains a high number of asparagine and glutamine residues, a hallmark of prion-like domains (PLDs) which are proposed to regulate phase separation of many proteins 1Alberti et al.200920Franzmann et al.201821Franzmann and Alberti2019. We used PLAAC 36Lancaster et al.2014, a web-based version of the algorithm used by 1Alberti et al.2009 to detect prion candidates in yeast proteome, to analyze the Ede1 sequence (Figure 5A). The algorithm detected a 99-amino acid-long prion-like sequence between amino acids 374 and 472, suggesting that this region could also be involved in the phase separation of endocytic proteins. We also consulted IUPred2a 49Mészáros et al.2018 and MobiDB-lite 55Necci et al.2017 algorithms to predict intrinsically disordered regions (IDRs) in Ede1. About 36% of Ede1 is predicted to be disordered; the unstructured regions are contained within the proline- and glutamine-rich region, and between the coiled-coil and the UBA domain.

To understand which features of Ede1 mediate phase-separation, we expressed a series of truncations of Ede1 in both wild-type and 3×ΔEA backgrounds, and analyzed their localization to condensates and endocytic sites (Figure 5B and C). In our 3×ΔEA background, the N- and C-terminal regions were dispensable for condensate formation. Surprisingly, the Ede1 central region consisting of amino acids 366–900 (the unstructured PQ-rich region and the coiled-coil domain) localized to large condensates in both the 3×ΔEA and wild-type backgrounds. It was also the minimal construct to form these condensates, as constructs containing only the coiled-coil or the PQ-rich region showed diffuse cytoplasmic localization.

For the wild-type background, our results are largely consistent with previously published results about the localization of Ede1 mutants 8Boeke et al.201443Lu and Drubin2017. Namely, the coiled-coil domain was necessary for Ede1 to assemble into endocytic sites, while the N- and C-terminal parts of Ede1 were individually dispensable.

We analyzed the concentration of truncated variants in both backgrounds by quantitative western blotting (Figure 5—figure supplement 1). The truncated Ede1 variants were expressed at higher levels than the full-length protein. This observation suggests that the loss of condensates in strains with truncated Ede1 is indeed caused by missing motifs rather than lowered concentration. We also confirmed that the concentration of full-length Ede1 increased in the 3×ΔEA background. It is unclear whether this represents a compensatory genetic mechanism, or if condensation is initially caused by reduced plasma membrane recruitment and subsequently interferes with protein degradation.

The functional significance of the Ede1 central region

To test the role of the central region of Ede1 in endocytosis, we created EGFP-tagged Ede1 mutants with internal deletions of amino acids 366–590 (Ede1ΔPQ), 591–900 (Ede1ΔCC), and 366–900 (Ede1ΔPQCC). Ede1ΔPQCC failed to localize to endocytic sites, whereas the two single-domain Ede1ΔPQ and Ede1ΔCC deletion mutants were still punctate, but more diffuse compared to full-length Ede1-EGFP (Figure 6A).

Ede1 features necessary for phase separation are also crucial for its function.

(A) Representative cells expressing full-length Ede1 and three internal Ede1 deletion mutants: Ede1ΔPQ (Δ366-590), Ede1ΔCC (Δ591-900), and Ede1ΔPQCC (Δ366-900) tagged with EGFP. (B) Representative cells expressing Sla1-EGFP and indicated Ede1 mutants.

# Data cleanup

strains <- read_csv("data/fig6_sla1-initiation/strain_lut.csv",
                    col_types = cols(ede1 = col_factor(NULL)))

sla1_density <- dir('data/fig6_sla1-initiation/density', pattern = 'dataset', full.names = T) %>%
  #map(dir, pattern = '.csv', full.names = T) %>%
  map(read_csv) %>%
  imap_dfr(~.x, .id = 'dataset')

# set pixel area; Orca pixels are 6.45um / 100x magnification 
pixel_area <- 0.0645**2

# clean up dataframe
sla1_density <- sla1_density %>%
  separate(Cell, c('date', 'strain', 'movie', 'cell', 'channel'), '_') %>% # separate filename into variables
  select(-c(channel, Threshold))%>% # remove stuff
  rename(cross_area = Cross_Area, patches = Patches) %>%
  left_join(strains, by = 'strain') %>%
  mutate(area = cross_area * 4 * pixel_area, density = patches / area)

# The first dataset was acquired in 2018
# with a shorter exposure time than repeats
# Therefore it has to be processed separately
# Let's see how different it really is before we include/exclude it
# from the analysis
frame_length <- 0.23

lifetime_00 <- list.files('data/fig6_sla1-initiation/lifetime/old_exposure',
                          pattern = '.txt', full.names = T) %>%
  set_names(nm = (basename(.) %>% tools::file_path_sans_ext())) %>% # names for the map output, wthout folder or extension
  map_df(read_table2, comment = '%%', col_names = F,
           col_types = cols(), .id = 'filename') %>%
  group_by(filename) %>%
  summarise(max = max(X1), min = min(X1)) %>%
  mutate(lifetime = (max - (min - 1)) * frame_length,
         dataset = '0',
         .keep = 'unused') %>%
  separate(filename, c('date', 'strain', 'movie', 'channel', 'processing', 'cell', 'md', 'trajectory'), '_') %>%
  select(-c(channel, processing, md)) %>%
  left_join(strains, by = 'strain')

frame_length <- 0.5

sla1_lifetime <- dir('data/fig6_sla1-initiation/lifetime', pattern = 'dataset', full.names = T) %>%
  map(dir, pattern = '.txt', full.names = T) %>%
  map(
    ~ set_names(.x,
                nm = (basename(.x) %>% tools::file_path_sans_ext()))) %>%
  map(
    ~ map_df(.x, read_table2,
             comment = '%%',
             col_names = F,
             col_types = cols(),
             .id = 'filename')) %>% # id observations by filename
  imap_dfr(~.x, .id = 'dataset') %>%
  group_by(filename, dataset) %>%
  summarise(max = max(X1), min = min(X1),
            .groups = 'drop') %>%
  mutate(lifetime = (max - (min - 1)) * frame_length) %>%
  select(filename, lifetime, dataset)

# Final touches on the lifetime, separated from main calculation
sla1_lifetime <- sla1_lifetime %>%
  separate(filename, c('date', 'strain', 'movie', 'cell', 'channel', 'processing', 'md', 'trajectory'), '_') %>%
  select(-c(channel, processing, md)) %>%
  left_join(strains, by = 'strain')

# Bind dataset 0 to the rest
sla1_lifetime <- bind_rows(sla1_lifetime, lifetime_00)

#Sla1 density

# couple of small functions for extracting comparisons
# into a format understandable by ggplot

#' Add letter group labels generated by Tukey's HSD
#'
#' This function performs ANOVA and Tukey's HSD,
#' extracts the letter labels and attaches them
#' to the original data.
#'
#' @param x: df or tibble, the data (long format!)
#' @param yvar: chr, name of the dependent variable
#' @param xvar: chr, name of the independent variable
#' @param alpha: dbl, confidence level passed on to Tukey's test
#'
add_tukey_labels <- function(x, yvar, xvar, alpha = 0.95){
  aov_form <- formula(paste(yvar, '~', xvar))
  anova <- aov(formula = aov_form, data = x)
  tukey <- TukeyHSD(anova, which = xvar, conf.level = alpha)
  # Extract labels and factor levels from Tukey post-hoc 
  x_labels <- as_tibble(
    multcompLetters4(anova, tukey)[[xvar]]$Letters,
    rownames = xvar
    )
  
  if (is.factor(x[[xvar]])) {
    x_labels[[xvar]] <- factor(
      x_labels[[xvar]], levels = levels(x[[xvar]])
    )
  }
  
  x_labels <- rename(x_labels, tukey_group = value)
  x <- left_join(x, x_labels, by = xvar)
  
  return(x)
}

#' Extract comparisons from rstatix tidy tests
#' 
#' This function subsets the selected comparisons in a Tukey,
#' Dunn or similar test done by rstatix.
#' It converts groups to a list of vectors that can be passed to geom_signif
#'
#'
#' @param x: df or tibble, the comparison results
#' @param rows: integer vector, the rows with desired comparisons
extract_comparisons <- function(x, rows){
  x_subset <- x[rows,] %>%
    .[nrow(.):1,]
  
  x_comparisons <- x_subset %>%
    select(group1, group2) %>%
    t() %>%
    as.data.frame() %>%
    as.list
  x_annotations <- x_subset$p.adj.signif %>%
    as.vector()
  
  significance <- list(
    comparisons = x_comparisons,
    annotations = x_annotations
  )
  
  return(significance)
}

sla1_density_stats <- sla1_density %>%
  group_by(ede1, dataset) %>%
  summarise(n = n(),
            across(density,
                   list(mean = mean, sd = sd, 
                        se = ~ sd(.x) / sqrt(n),
                        median = median, mad = mad)),
            .groups = 'drop')

#sla1_density_stats %>% knitr::kable()

plot_blank <- ggplot(sla1_density_stats,
                     aes(x = ede1, y = density_mean)) +
                         #shape = dataset, fill = dataset))+
  labs(title = NULL, x = 'Ede1', y = expression( "Sla1 patches/µm"^2)) +
  scale_y_continuous(breaks = seq(0.0, 1.0, 0.1)) +
  scale_shape_manual(values = c(21:25)) +
  scale_color_brewer(palette = 'Set2') +
  scale_fill_brewer(palette = 'Set2')

plot_scatter <- plot_blank +
  geom_quasirandom(inherit.aes = F, data = sla1_density,
                   aes(x = ede1, y = density,
                       shape = dataset,# colour = dataset
                       ),
                   colour = 'grey75',# shape = 1,
                   show.legend = F, size = 0.8
                   )

plot_violin <- plot_blank + 
  geom_violin(inherit.aes = F,
              data = sla1_density, aes(x = ede1, y = density),
              colour = 'grey75', fill = 'transparent'
              )

plot_super <- plot_scatter +
  geom_quasirandom(aes(shape = dataset, fill = dataset),
                   show.legend = F,
                   width = 0.3, size = 2)+
  stat_summary(fun = mean, geom = 'crossbar',
               width = 0.5, fatten = 1)+
  stat_summary(fun.data = 'mean_sdl',
               fun.args = list(mult = 1), 
               geom = 'errorbar', width = 0.2)
  
#plot_super

plot_super_violin <- plot_violin +
  geom_quasirandom(aes(shape = dataset, fill = dataset),
                   show.legend = F,
                   width = 0.3, size = 2)+
  stat_summary(fun = mean, geom = 'crossbar',
               width = 0.5, fatten = 1)+
  stat_summary(fun.data = 'mean_sdl',
               fun.args = list(mult = 1),
               geom = 'errorbar', width = 0.2)

#plot_super_violin

tukey <- tukey_hsd(sla1_density_stats, density_mean ~ ede1,
                   ordered = TRUE)

significance <- extract_comparisons(tukey, c(1:4, 10))

plot_super_signif <- plot_super +
  geom_signif(comparisons = significance$comparisons,
              annotations = significance$annotations,
              step_increase = 0.03,
              tip_length = 0.01, vjust = 0.8,
              margin_top = -0.1)
plot_super_signif

Sla1 patch density in Ede1 mutants.

Large points represent mean measurements from independently repeated datasets. Central line and whiskers denote the mean ± SD calculated from dataset averages. Gray points show individual observations. Letters denote pairwise comparisons based on Tukey-Kramer test; groups which do not share any letters are significantly different at α = 0.05. Scale bars: 2 μm.

#sla1 lifetime
sla1_lifetime_stats <- sla1_lifetime %>%
  group_by(ede1, dataset) %>%
  summarise(n = n(),
            across(lifetime,
                   list(mean = mean, sd = sd, 
                        se = ~ sd(.x) / sqrt(n()),
                        median = median, mad = mad)),
            .groups = 'drop')

#sla1_lifetime_stats %>% knitr::kable()

plot_blank <- ggplot(sla1_lifetime_stats,
                     aes(x = ede1, y = lifetime_mean)) +
                         #shape = dataset, fill = dataset))+
  labs(title = NULL, x = 'Ede1', y = "Sla1 lifetime (s)") +
  scale_y_continuous(breaks = scales::breaks_extended(6))+
  scale_shape_manual(values = c(21:25)) +
  scale_color_brewer(palette = 'Set2') +
  scale_fill_brewer(palette = 'Set2')

plot_scatter <- plot_blank +
  geom_quasirandom(inherit.aes = F, data = sla1_lifetime,
                   aes(x = ede1, y = lifetime,
                       shape = dataset,# colour = dataset
                       ),
                   colour = 'grey75',# shape = 1,
                   show.legend = F, size = 0.8
                   )

plot_violin <- plot_blank + 
  geom_violin(inherit.aes = F,
              data = sla1_lifetime, aes(x = ede1, y = lifetime),
              colour = 'grey75', fill = 'transparent'
              )

plot_super <- plot_scatter +
  geom_quasirandom(aes(shape = dataset, fill = dataset),
                   show.legend = F,
                   width = 0.3, size = 2)+
  stat_summary(fun = mean, geom = 'crossbar',
               width = 0.5, fatten = 1)+
  stat_summary(fun.data = 'mean_sdl',
               fun.args = list(mult = 1), 
               geom = 'errorbar', width = 0.2)
  
#plot_super

plot_super_violin <- plot_violin +
  geom_quasirandom(aes(shape = dataset, fill = dataset),
                   show.legend = F,
                   width = 0.3, size = 2)+
  stat_summary(fun = mean, geom = 'crossbar',
               width = 0.5, fatten = 1)+
  stat_summary(fun.data = 'mean_sdl',
               fun.args = list(mult = 1),
               geom = 'errorbar', width = 0.2)

#plot_super_violin

tukey <- tukey_hsd(sla1_lifetime_stats, lifetime_mean ~ ede1,
                   ordered = TRUE)

significance <- extract_comparisons(tukey, c(1:4, 10))

plot_super_signif <- plot_super +
  geom_signif(comparisons = significance$comparisons,
              annotations = significance$annotations,
              step_increase = 0.03,
              tip_length = 0.01, vjust = 0.8,
              margin_top = -0.1)
plot_super_signif

Sla1 patch lifetime in Ede1 mutants.

Large points represent mean measurements from independently repeated datasets. Central line and whiskers denote the mean ± SD calculated from dataset averages. Gray points show individual observations. Letters denote pairwise comparisons based on Tukey-Kramer test; groups which do not share any letters are significantly different at α = 0.05. Scale bars: 2 μm.

Recruitment of Ede1ΔPQCC to other proteins cannot rescue the endocytic defect.

Ede1ΔPQCC-mCherry-FKBP was recruited to Syp1-FRB or Sla2-FRB by addition of 10 μg ml−1 of rapamycin to the growth medium. (A) Ede1ΔPQCC-mCherry-FKBP signal in cells coexpressing Syp1-FRB and Sla1-EGFP. (C) Ede1ΔPQCC-mCherry-FKBP signal in cells coexpressing Sla2-FRB and Sla1-EGFP. (D) Density of Sla1-EGFP patches in the same cells cultured with or without rapamycin. P-value from Welch’s t-test. Scale bars: 2 μm.

strains <- read_csv("data/fig6_syp1-fusion/strain_lut.csv", 
                    col_types = cols(ede1 = col_factor(NULL))) # order factor levels in csv file

sla1_density <- dir('data/fig6_syp1-fusion', pattern = 'dataset', full.names = T) %>%
  #map(dir, pattern = '.csv', full.names = T) %>%
  map(read_csv) %>%
  imap_dfr(~.x, .id = 'dataset')

# set pixel area; Orca pixels are 6.45um / 100x magnification 
pixel_area <- 0.0645**2

sla1_density <- sla1_density %>%
  separate(Cell, c('date', 'strain', 'movie', 'cell', 'channel'), '_') %>% # separate filename into variables
  select(-c(channel, Threshold))%>% # remove stuff
  rename(cross_area = Cross_Area, patches = Patches) %>%
  mutate(area = cross_area * 4 * pixel_area, density = patches / area) %>%
  left_join(strains, by = 'strain')

sla1_density_stats <- sla1_density %>%
  group_by(ede1, dataset) %>%
  summarise(n = n(),
            across(density,
                   list(mean = mean, sd = sd, median = median, mad = mad,
                        shapiro.p = ~ tidy(shapiro.test(.x))$p.value)),
            .groups = 'drop')

#sla1_density_stats %>% knitr::kable()

anova <- anova_test(sla1_density_stats, density_mean ~ ede1)

tukey <- tukey_hsd(sla1_density_stats, density_mean ~ ede1,
                   ordered = TRUE)

plot_blank <- ggplot(sla1_density_stats,
                     aes(x = ede1, y = density_mean))+
  labs(title = NULL, x = NULL, y = expression("Sla1 patches/µm"^2))+
  scale_y_continuous(breaks = seq(0.0, 1.0, 0.1))+
  scale_shape_manual(values = c(21:25))+
  scale_color_brewer(palette = 'Pastel2')+
  scale_fill_brewer(palette = 'Set2')

plot_scatter <- plot_blank + 
  geom_quasirandom(data = sla1_density,
                   aes(x = ede1, y = density,
                     shape = dataset, #colour = dataset
                     ),
                 colour = 'grey75',
                 show.legend = F, size = 0.8
                 )

plot_violin <- plot_blank + 
  geom_violin(inherit.aes = F,
              data = sla1_density, aes(x = ede1, y = density),
              colour = 'grey70', fill = 'transparent'
              )

plot_super_violin <- plot_violin +
  geom_quasirandom(aes(shape = dataset, fill = dataset),
                   show.legend = F,
                   width = 0.3, size = 2)+
  stat_summary(fun = mean, geom = 'crossbar',
               width = 0.5, fatten = 1)+
  stat_summary(fun.data = 'mean_sdl',
               fun.args = list(mult = 1), 
               geom = 'errorbar',
               width = 0.2)

plot_super <- plot_scatter +
  geom_quasirandom(aes(shape = dataset, fill = dataset),
                   show.legend = F,
                   width = 0.3, size = 2)+
  stat_summary(fun = mean, geom = 'crossbar',
               width = 0.5, fatten = 1)+
  stat_summary(fun.data = 'mean_sdl',
               fun.args = list(mult = 1), 
               geom = 'errorbar',
               width = 0.2)
  
significance <- extract_comparisons(tukey, c(1:6))

plot_super_signif <- plot_super +
  geom_signif(comparisons = significance$comparisons,
              annotations = significance$annotations,
              step_increase = 0.03,
              tip_length = 0.01, vjust = 0.8,
              margin_top = -0.1)

plot_super_signif

Density of Sla1-EGFP patches in cells expressing wild-type Ede1, Ede1ΔPQCC-mCherry-FKBP, and Syp1-FRB cultured with or without rapamycin, or no Ede1.

Large points: mean measurements from independently repeated datasets. Central line and whiskers: mean ± SD calculated from dataset averages. Gray points: individual cells from all datasets. Statistical significance of pairwise comparisons was determined by Tukey-Kramer test; n.s., not significant; ***, p<0.001.

Next, we tested if these Ede1 mutants had endocytic defects by using Sla1 as a reporter of the late phase of endocytosis. We tagged Sla1 with EGFP in Ede1 mutant strains (Figure 6B) and quantified the density and lifetimes of endocytic sites (Figure 6C and D). In ede1Δ and ede1ΔPQCC cells, the mean number of endocytic events marked by Sla1-EGFP per μm2 was reduced by 46 and 43% of the wild-type, respectively. Consistent with their effects on Ede1 recruitment, the ede1ΔPQ and ede1ΔCC mutations caused intermediate reduction in patch density (by 24 and 26%, respectively). All differences from the wild type were statistically significant (p<0.001 in Tukey-Kramer test). The difference between ede1Δ and ede1ΔPQCC was not statistically significant (p=0.86).

The Sla1 lifetimes were likewise affected by the deletion of Ede1 central region. In ede1Δ, Sla1-EGFP lifetime was decreased by 29% and in ede1ΔPQCC, by 28%. The deletions of individual regions again showed intermediate defects (13 and 18% for Ede1ΔPQ and Ede1ΔCC, respectively). All differences from the wild type were statistically significant (p<0.01 in Tukey-Kramer test). The deletion of the entire EDE1 gene was not significantly different from the ede1ΔPQCC mutant (p = 0.99).

The Ede1ΔPQCC mutant does not localize to endocytic sites. We wanted to test if the N- and C-terminal domains alone could support endocytosis if a strong interaction with another endocytic protein was introduced. We generated cells coexpressing Ede1ΔPQCC-mCherry-FKBP and Syp1-FRB or Sla2-FRB in order to link the Ede1 mutant to another component of the early coat via the rapamycin-inducible FKBP-FRB dimerization system 25Haruki et al.2008. Recruitment to Syp1 caused the Ede1ΔPQCC-mCherry-FKBP signal to become more prominent around the bud necks, but did not rescue membrane patch formation, whereas recruitment to Sla2 partially rescued the patch localization of Ede1ΔPQCC (Figure 6—figure supplement 1). The average Sla1 density did not significantly change in either of these yeast strains upon rapamycin treatment (p=0.53 and p=0.59 between treated and untreated cells in Syp1-FRB and Sla2-FRB strains, respectively).

Taken together, our results show that the central region is essential for Ede1 to promote efficient endocytosis and to regulate the timing of coat maturation.

In ede1Δ cells, many of the early endocytic proteins fail to localize to endocytic sites 75Stimpson et al.200910Carroll et al.2012. We therefore visualized the localization of early proteins in Ede1 mutants lacking the central region.

Different proteins were affected by the Ede1 central deletions in different ways (Figure 7), consistent with the work done on ede1Δ mutants. The localization of Apl1 (β-subunit of the AP-2 complex) was the most severely disrupted. Apl1-EGFP patches were less defined in Ede1ΔPQ background, and undetectable in Ede1ΔCC or Ede1ΔPQCC cells. In these cells, Apl1-EGFP signal was dispersed in the cytoplasm, with a faint presence around the bud neck. Syp1 and Yap1801 remained localized to the membrane in all of the mutants, but the signal became more diffuse along the membrane and the patches less defined. This effect was the strongest for Ede1ΔPQCC, with intermediate effects in Ede1ΔPQ and Ede1ΔCC cells. However, Ent1 and Sla2 still assembled into endocytic patches in all of the Ede1 mutants. Taken together, our results indicate that the Ede1 central region is essential to concentrate early endocytic proteins.

Ede1 central deletion mutants are defective in early protein localization.

Maximum-intensity projections of 3D volumes are shown for representative cells with different early proteins tagged with EGFP. The strains express Ede1 species indicated at the top. Scale bars: 2 μm.

Ede1 central region can be replaced by other prion-like domains

We wanted to test whether the loss of function in ede1ΔPQCC cells in respect to Sla1 density could be rescued by heterologous protein sequences, such as IDRs, globular domains, or coiled-coils (Figure 8A).

Prion-like domains can partially replace Ede1 central region.

(A) Domain structure of the Ede1 central region replacement constructs. Amino acids 366–900 of Ede1 were replaced with prion-like intrinsically disordered region (IDR) sequences, monomeric (mCherry) or dimeric (dTomato) fluorescent proteins, and dimeric (Khc) or tetrameric (Eg5) coiled-coils. All mutants were expressed from the Ede1 locus under the control of the native promoter. (B) Representative cells expressing indicated Ede1 mutants tagged C-terminally with msGFP2.

#Data cleanup
strains <- read_csv("data/fig8_sla1-density/strain_lut.csv", 
                    col_types = cols(ede1 = col_factor(NULL))) # order factor levels in csv file

sla1_density <- dir('data/fig8_sla1-density', pattern = 'dataset', full.names = T) %>%
  #map(dir, pattern = '.csv', full.names = T) %>%
  map(read_csv) %>%
  imap_dfr(~.x, .id = 'dataset')

# set pixel area; Orca pixels are 6.45um / 100x magnification 
pixel_area <- 0.0645**2

sla1_density <- sla1_density %>%
  separate(Cell, c('date', 'strain', 'movie', 'cell', 'channel'), '_') %>% # separate filename into variables
  select(-c(channel, Threshold))%>% # remove stuff
  rename(cross_area = Cross_Area, patches = Patches) %>%
  mutate(area = cross_area * 4 * pixel_area, density = patches / area) %>%
  left_join(strains, by = 'strain')

#Sla1 patch density in Ede1 internal replacement mutants
#sla1_density %>%
#  group_by(ede1, dataset) #%>%
#  summarise(n = n(),
#            across(c(patches, area), list(mean = mean, sd = sd,
#                                          se = ~sd(.x) / sqrt(n()))))

sla1_density_stats <- sla1_density %>%
  group_by(ede1, dataset) %>%
  summarise(n = n(),
            across(density,
                   list(mean = mean, sd = sd, 
                        se = ~ sd(.x) / sqrt(n),
                        median = median, mad = mad)),
            .groups = 'drop')

plot_blank <- ggplot(sla1_density_stats,
                     aes(x = ede1, y = density_mean)) +
                         #shape = dataset, fill = dataset))+
  labs(title = NULL, x = 'Ede1', y = expression( "Sla1 patches/µm"^2)) +
  scale_y_continuous(breaks = seq(0.0, 1.0, 0.1)) +
  scale_shape_manual(values = c(21:25, 8)) +
  scale_color_brewer(palette = 'Set2') +
  scale_fill_brewer(palette = 'Set2')

plot_scatter <- plot_blank +
  geom_quasirandom(inherit.aes = F, data = sla1_density,
                   aes(x = ede1, y = density,
                       shape = dataset,# colour = dataset
                       ),
                   colour = 'grey75',# shape = 1,
                   show.legend = F, size = 0.8
                   )

plot_violin <- plot_blank + 
  geom_violin(inherit.aes = F,
              data = sla1_density, aes(x = ede1, y = density),
              colour = 'grey75', fill = 'transparent'
              )

plot_super <- plot_scatter +
  geom_quasirandom(aes(shape = dataset, fill = dataset),
                   show.legend = F,
                   width = 0.3, size = 2)+
  stat_summary(fun = mean, geom = 'crossbar',
               width = 0.5, fatten = 1)+
  stat_summary(fun.data = 'mean_sdl',
               fun.args = list(mult = 1), 
               geom = 'errorbar', width = 0.2)+
  guides(x = guide_axis(angle = -45))

plot_super_flip <- plot_super +
  coord_flip()+
  scale_x_discrete(limits = rev(levels(sla1_density$ede1)))+
  theme(panel.border = element_blank(),
        axis.line.y = element_blank(),
        axis.ticks.y = element_blank(),
        axis.line.x = element_line())

plot_super_violin <- plot_violin +
  geom_quasirandom(aes(shape = dataset, fill = dataset),
                   show.legend = F,
                   width = 0.3, size = 2)+
  stat_summary(fun = mean, geom = 'crossbar',
               width = 0.5, fatten = 1)+
  stat_summary(fun.data = 'mean_sdl',
               fun.args = list(mult = 1),
               geom = 'errorbar', width = 0.2)

sla1_density_stats <- sla1_density_stats %>%
  add_tukey_labels('density_mean', 'ede1')

plot_super_letters <- plot_super +
  geom_text(data = sla1_density_stats,
            aes(label = tukey_group), y = Inf, vjust = 1.2)
  #stat_summary(aes(label = tukey_group),
  #             fun.y = Inf,
  #             geom = 'text', na.rm = T, vjust = -0.5,
  #             #fun.args = list(err = errors, mult = error_range)
  #             )

plot_super_letters

Quantification of Sla1-EGFP patch density in strains expressing indicated Ede1 mutants.

Large points represent mean measurements from independently repeated datasets. Central line and whiskers denote the mean ± SD calculated from dataset averages. Gray points show individual cells from all datasets. At α = 0.05, all mutants are significantly different from wild type, and groups marked with an asterisk (*) are significantly different from Ede1ΔPQCC (Tukey-Kramer test; a complete table of pairwise comparisons and effect sizes can be found in Figure 8—source data 1). All scale bars: 2 μm.

Ede1FUS and Ede1Whi3 intrinsically disordered region replacement constructs.

Representative cells expressing Ede1FUS and Ede1Whi3 tagged C-terminally with msGFP2.

First, we replaced amino acids 366–900 of Ede1 with different IDRs. We considered three factors for choosing the replacement IDRs: known phase separation activity, Q/N content, and similarity to the Ede1 IDR sequence. We chose the region of Sup35 spanning the N and M regions based on its known phase separation activity in yeast cells 20Franzmann et al.2018. We chose the IDR of Snf5 because of its high score in the prion-like screen by 1Alberti et al.2009, reflecting a Gln-rich amino acid sequence, and the IDR of Whi3, because the method developed by 84Zarin et al.2019 indicated it as one of the sequences most similar to the PQ-rich region of Ede1. Finally, we also included the low-complexity domain of human FUS for its well-known tendency to phase separate and form hydrogels 58Patel et al.201530Kato et al.2012.

We also replaced the central region of Ede1 with several structured domains. We used the red fluorophores mCherry and dTomato as globular linkers assuming respectively monomeric and dimeric states. We also used the coiled-coil domains from two kinesin motors as oligomeric rod-shaped linkers: the kinesin-1 heavy chain (Khc amino acids 335–931) from Drosophila melanogaster and the human kinesin-5 Eg5 (amino acids 358–797), which form dimers and tetramers in their respective contexts 18Cuevas et al.199270Scholey et al.2014.

All of the Ede1 central region replacement constructs were expressed from the endogenous genomic locus under the control of the native promoter.

We assessed the localization of Ede1 constructs tagged with msGFP2 (Figure 8B, Figure 8—figure supplement 1). We found that the localization was partially rescued by the insertion of Whi3, Snf5, and Sup35 IDRs, as well as the oligomeric linkers (dTomato, Khc, Eg5), in place of the deleted central region. The localization was not rescued by the FUS low-complexity region or the monomeric globular linker mCherry.

We also assessed the density of Sla1-EGFP patches in cells expressing untagged Ede1mCherry, Ede1dTomato, Ede1Khc, Ede1Eg5, Ede1Snf5, or Ede1Sup35 (Figure 8C). All differences from wild type were statistically significant (p<0.001). No single replacement mutant was able to fully rescue the Sla1 patch density defect present in ede1ΔPQCC cells. Only ede1Sup35 was able to significantly rescue the Sla1 density in respect to the deletion of Ede1 central region. In ede1ΔPQCC and ede1Δ cells, the mean density of Sla1 patches was reduced by 40 and 48%, respectively. In ede1Sup35, the mean density was only reduced by 25% from wild type, respectively.

These results support the hypothesis that prion-like domains can aid clustering of the endocytic proteins.

Ede1 central region can cluster a heterologous lipid-binding protein

We hypothesized that the phase separation mediated by Ede1 central region is able to cluster membrane-associated proteins. To test our hypothesis, we fused the Ede1 central region to a diffusely membrane-bound protein. We created a GFP-Ede1366-900-2×PH(PLCδ) construct, based on a phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2) probe developed by 74Stefan et al.2002, by inserting the Ede1 central region between GFP and a tandem repeat of pleckstrin homology (PH) domain of phospholipase C δ1 (PLCδ). The original GFP-2×PH(PLCδ) construct is distributed homogeneously on the plasma membrane, while GFP-Ede1366-900 alone localized to bright intracellular condensates. In contrast, the fusion construct localized to the plasma membrane, forming puncta that resembled endocytic sites (Figure 9A,B).

Fusion with Ede1 central region changes the distribution of a PI(4,5)P2 probe.

(A) Maximum projections of cells expressing GFP-2 × PH and GFP-Ede1366-900 from a yeast centromeric plasmid under the control of TDH3 promoter. White dotted line shows cell outline. (B) Cells expressing GFP-Ede1366-900-2×PH from a multicopy plasmid. Examples of different structures classified as ‘D_iffuse’, ‘_P_atches’, or ‘_N_etworks’. (C) and (D) The same construct was expressed from four centromeric plasmids under different promoters. Individual cells were classified as in panel B. Plots show percentage of cells falling into each class per promoter (C), and mean cell pixel intensity per class (D). (E) Movies of cells coexpressing GFP-Ede1366-900-2×PH and Sla1-mCherry were acquired using TIRF microscopy. Single frame from a representative movie; points labeled ‘1’ and ‘2’ mark the top and bottom of the kymograph (F). Scale bars: 2 μm. All cells in this figure: _SLA1-mCherry::KANMX4, ede1Δ::natNT2.

Fluorescence recovery after photobleaching and hexanediol treatment of the fusion construct.

(A) Kymographs of representative photobleaching experiments performed on the GFP-2×PH or GFP-Ede1366-900-2×PH constructs. (B) Effects of 10% 1,6-hexanediol (HD) treatment on cells overexpressing GFP-2 × PH, GFP-Ede1366-900-2×PH, or GFP-Ede1366-900 from a single-copy plasmid under the control of the TDH3 promoter.

We also noticed subpopulations of cells with different localization patterns of the construct (Figure 9B). We speculated that the variable patterns were caused by heterogeneity in protein expression level due to plasmid copy number variation. To test that hypothesis, we expressed GFP-Ede1366-900-2×PH(PLCδ) from centromeric plasmids containing four different promoters of increasing strength 53Mumberg et al.1995. We classified the localization of the construct in these cells as ‘diffuse’, ‘punctate’, or ‘networked’. We found that the tendency to cluster into different patterns correlated with promoter strength and the expression level. Low expressing cells had more diffuse localization of the construct, and separated into puncta or well-separated regions as the concentration increased (Figure 9C,D).

The puncta formed by the GFP-Ede1366-900-2×PH(PLCδ) construct were stable over long imaging periods, but dynamically recruited the late coat protein Sla1 (Figure 9F, Figure 9—video 1, Figure 9—video 2). Sla1-mCherry persisted at these sites with similar lifetimes as during normal endocytosis indicating that the condensates can recruit endocytic coat components. Unlike full-length Ede1 at the endocytic sites, the chimeric construct does not undergo cycles of assembly and disassembly. When we photobleached the structures formed by highly expressed GFP-Ede1366-900-2×PH(PLCδ) we saw no recovery (Figure 9—figure supplement 1A), indicating that without the terminal domains, Ede1 central region might form solid, rather than liquid-like, structures. 10% 1,6-hexanediol also failed to dissolve the GFP-Ede1366-900-2×PH(PLCδ) structures (Figure 9—figure supplement 1B). Structures formed by strongly overexpressed GFP-Ede1366-900 were partially dissolved by 10% 1,6-hexanediol. This suggests that the stability of GFP-Ede1366-900-2×PH(PLCδ) might also be modulated by the PH domains interacting with the membrane. In addition, we noticed that 1,6-hexanediol treatment caused membrane deformations, or possible clustering of the PH domains expressed without the Ede1 fusion (Figure 9—figure supplement 1B). This observation is consistent with the previous reports of wide-ranging effects of 1,6-hexanediol 33Kroschwald et al.2015, and underscores that 1,6-hexanediol experiments need to be interpreted with caution.

These results show that directing the Ede1 central region to the plasma membrane is sufficient to create puncta on the membrane in a concentration-dependent manner. Surprisingly, these long-lived sites can repeatedly recruit endocytic coat components.

Discussion

Ede1 is the key organizer of the early phase of endocytosis in yeast 75Stimpson et al.20098Boeke et al.201443Lu and Drubin2017. Our results indicate that the large clusters of Ede1 observed previously in mutant cells 8Boeke et al.2014 are in fact phase-separated protein droplets. Moreover, we found that the cytosolic concentration of Ede1 in normal cells is at the same critical limit as in the mutant cells harboring Ede1 droplets. This suggests that liquid phase separation might be the mechanism through which Ede1 concentrates proteins at the early endocytic sites. We identified the central region of Ede1— containing a coiled-coil and a prion-like domain— as necessary for both the condensate formation, and for Ede1 to promote the initiation of endocytosis. We also found that heterologous prion-like domains can partially replace the Ede1 central region in endocytosis. We also demonstrated that the central region of Ede1 fused to a lipid-binding domain can condense on the plasma membrane. These findings suggest a potential link between endocytic assembly and the phenomenon of protein phase separation and raise questions about the material properties of the endocytic sites at different stages. They also highlight a possible novel role for disordered, prion-like regions found in numerous endocytic proteins 47Malinovska et al.2013.

Ede1 forms liquid protein droplets

Ede1 forms large condensates under conditions where the stoichiometry between Ede1 and the endocytic adaptor proteins is altered, such as overexpression of Ede1 or deletion of three early adaptors (8Boeke et al.2014, Figure 1).

We showed that these condensates are liquid, phase-separated droplets according to the following criteria: (a) observation of liquid-like behaviors, (b) molecule turnover, (c) dependence on a critical component concentration, (d) dependence on temperature, and (e) susceptibility to dissolution by 1,6-hexanediol.

We show that the Ede1 condensates undergo apparent fusion and fission events, the latter caused possibly by polymerizing actin filaments (Figure 4C, Figure 4—video 1). The Ede1 molecules exchange between the condensates and the cytosolic pool (Figure 2A) and, importantly, Ede1 molecules also rapidly diffuse within the condensates (Figure 2C and D). The formation of the Ede1 condensates depends on the cellular concentration of Ede1 (Figures 1D and !number(3)). The condensates dissolve rapidly and reversibly in response to temperature changes (Figure 2F), and are sensitive to treatment with 1,6-hexanediol (Figure 2E). In agreement with our findings in yeast, 17Day et al.2021 have shown that Eps15, the mammalian homolog of Ede1, can phase separate in vitro.

The Ede1 condensates are clearly distinct from endocytic sites by the virtue of size, brightness, and long-term stability. We assign them no function, other than as a tool used to study the properties of Ede1. However, 81Wilfling et al.2020 studied Ede1 condensates in the 3×ΔEA and Ede1 overexpression strains in parallel to our work. They describe a selective autophagy pathway mediated by Ede1, and propose that autophagy of phase-separated condensates might be a major route through which cells remove misfolded or unneeded endocytic proteins. The ability of Ede1 to cluster endocytic proteins could therefore play a dual role in endocytosis and autophagy.

Endocytic sites: solid or liquid?

The prevailing model of endocytosis focuses on the growing clathrin lattice as the driver of protein assembly 31Kirchhausen et al.201413Cocucci et al.2012. Indeed, clathrin is a major interaction hub in endocytosis, and a scaffold with a well-defined structure. Nevertheless, several lines of evidence suggest that this model of assembly might be incomplete.

In yeast, the early proteins can assemble in the absence of clathrin. In fact, many endocytic sites in clc1Δ cells stall during the early phase, but persist on the plasma membrane 10Carroll et al.2012. On the other hand, in the absence of Ede1, numerous early adaptors do not assemble at all, or assemble only for the duration of the membrane-bending phase 75Stimpson et al.200910Carroll et al.2012. These two observations demonstrate that Ede1 can sustain the early sites independently of the clathrin lattice. We hypothesize that Ede1 performs this function by undergoing liquid-liquid phase separation on the plasma membrane.

FRAP shows that numerous proteins involved in endocytosis and the actin cytoskeleton continuously turn over at the endocytic sites 72Skruzny et al.201235Lacy et al.2019. Ede1 is one of such proteins (Figure 2B). This observation is consistent with a phase separation mechanism, although it does not prove it. Alternative mechanisms could explain fast turnover, such as ‘treadmilling’ of actin monomers or dynamic binding of adaptors to the clathrin lattice. Even clathrin itself shows fluorescence recovery as individual triskelia are replaced within the scaffold 83Wu et al.20013Avinoam et al.201512Chen et al.2019.

Several hallmark criteria frequently associated with liquid phase separation can also be explained by other mechanisms of compartmentalization. This fact has become intensely debated in the context of the nucleus, where bridging of multiple DNA sites could create an appearance of phase-separated compartments 48McSwiggen et al.201961Peng and Weber2019. Similarly, weak binding to the clathrin lattice could explain the enrichment of proteins at endocytic sites with high turnover and susceptible to dissolution by temperature and 1,6-hexanediol.

Our observations point to a mixed model, in which the structured lattice exists alongside a liquid phase formed by unstructured interactions. The existence of large Ede1 condensates is in itself one of the predictions generated by a phase separation model. A significant consequence of phase separation is that above a critical value, further increase of total component concentration leads to changes in relative volume, but not concentration, of the dense and the light phases. This is in contrast with the scaffold-binding model, where the size of the assemblies formed by endocytic proteins would be limited by the clathrin lattice. In addition, we have observed buffering of cytoplasmic Ede1 concentration during overexpression and in different genetic backgrounds. The cytoplasmic concentration in wild-type cells reached the limit observed in Ede1 overexpression, suggesting that the normal Ede1 concentration is sufficient for its phase separation on the plasma membrane. The small fraction of wild-type cells which contain larger Ede1 condensates (Figures 1E and 5C, 81Wilfling et al.2020) could very well be the consequence of natural variability in expression levels. The concentration buffering also suggests that the phase separation of Ede1 is driven primarily by homotypic interactions 65Riback et al.2020.

The qualitative criteria for liquid phase separation—sensitivity to temperature and 1,6-hexanediol—apply to Ede1 at the endocytic sites as well as to the large endocytic condensates. Curiously, the sites appear more stable against both of these treatments than the large condensates. This could be explained by the fact that adaptor binding confines Ede1 at the endocytic sites to the plane of the membrane. As such, the critical concentration could be lower than for the formation of 3D droplets, as shown previously in vitro for the Nephrin/Nck/N-WASP system 6Banjade and Rosen2014. A complementary explanation could be that genuine endocytic sites are stabilized by other interactions, such as those within the clathrin lattice, or the lattice formed by Sla2 and Ent1/2 in the presence of PI(4,5)P2.

The small size of the endocytic sites prevents the direct visualization of liquid-like shape changes. Super-resolution imaging of endocytic proteins in fixed cells revealed that Ede1 and Syp1 form larger and more amorphous structures than clathrin and its adaptors 54Mund et al.2018. The super-resolution experiments also suggest that even before it disassembles, Ede1 becomes progressively excluded from the center of the sites. Eps15 and FCHo1 behave in a similar fashion in mammalian cells 73Sochacki et al.2017. As the intermediate and late coat proteins form stable patches with low turnover 72Skruzny et al.2012, we propose that a liquid-like early module is displaced from the center of the invagination by the formation of a solid coat.

Ede1 is one of the most heavily phosphorylated proteins in yeast endocytosis 42Lu et al.2016, and it can be ubiquitylated as well as bind ubiquitin. Phosphorylation can regulate phase separation 51Monahan et al.201737Larson et al.2017, and the phosphorylation of Ede1 might regulate its state at the endocytic sites.

Our chimeric construct GFP-Ede1PQCC-2×PH combines a lipid-binding domain with the central region of Ede1. This construct forms bright structures on the surface of the plasma membrane, and the area covered by these structures appears larger in cells with higher expression levels (Figure 9). The puncta formed by GFP-Ede1PQCC-2×PH repeatedly recruit transient assemblies of the late coat protein Sla1. This suggests that the chimeric construct can initiate functional endocytic events. However, the GFP-Ede1PQCC-2×PH puncta persist over long imaging periods and do not disassemble after Sla1 internalization. Structures formed by GFP-Ede1PQCC-2×PH also do not recover fluorescence (Figure 9—figure supplement 1). This suggests that while the central region of Ede1 can drive condensation, the terminal regions are needed to maintain the liquid state.

It must be noted that we do not fully understand the nature of the microdomains formed by the GFP-Ede1PQCC-2×PH construct. For example, all known interaction motifs of Ede1 are located inside the terminal regions. It is thus unclear how the central region could recruit other endocytic proteins.

Prion-like domains are enriched in yeast endocytic proteins 1Alberti et al.200947Malinovska et al.2013, many of which contain polyglutamine tracts longer than that of Ede1. These disordered, low-complexity regions had been previously considered mere linkers 16Dafforn and Smith2004, but we show that the prion-like region of Ede1 is important for its condensation. Our results also show that unrelated prion-like domains can partially replace the function of Ede1 central region, even without the coiled-coil domain (Figure 8). How—or if—the endocytic proteins achieve specificity during recruitment of prion-like domains is an open question. It has also been proposed that phase separation of prion-like domains could provide force for membrane bending 7Bergeron-Sandoval et al.2021.

The coiled-coil domain of Ede1 is also critical for its phase separation, and even more important to the function of Ede1 than the prion-like domain (Figure 6; 43Lu and Drubin2017). The coiled-coil of Eps15 can form dimers and tetramers 79Tebar et al.199715Cupers et al.1997, and fluorescence correlation spectroscopy (FCS) data suggests that cytosolic Ede1 can form dimers and higher-order oligomers 8Boeke et al.2014. Multivalency is a critical characteristic of phase-separating proteins 39Li et al.20124Banani et al.2016, and oligomerization via the coiled-coil domain could promote phase separation by increasing the valency of other interactions. Less commonly, coiled-coils can also form phase-separating networks in absence of other interaction motifs, such as in the case of centrosome scaffold SPD-5 82Woodruff et al.2017. In our experiments, replacing the central region of Ede1 with heterologous coiled-coils partially rescued Ede1 localization, but not the late phase defect, of ede1ΔPQCC cells.

Materials and methods

Yeast strains and plasmids

The list of yeast strains and yeast plasmids used in this study is provided in Supplementary file 1. These materials are available upon request to the corresponding author.

Cells were maintained on rich medium at 24 or 30°C. C-terminally tagged or truncated mutants were generated via homologous recombination with PCR cassettes as described by 27Janke et al.2004. N-terminal truncation and internal domain deletion or replacement mutants of Ede1 were generated by first constructing the desired mutant gene in a pET-based plasmid using PCR mutagenesis or ligation-independent cloning 39Li et al.2012. The mutated Ede1 sequence was then amplified by PCR using primers containing 50 bp overlap with 5’ (forward primer) 3’ (reverse primer) UTR sequences of EDE1. The PCR product was transformed into ede1Δ::klURA3 cells. The transformants were selected for on plates containing 5-fluorootic acid, and confirmed by colony PCR and genomic sequencing.

Sequences coding for IDRs replacing Ede1 core region in Figure 8 were obtained from several different sources. The coding sequences of Snf5 and Whi3 IDRs were amplified by PCR from the yeast genome and were confirmed identical to the sequences in the S288C reference genome (SGD:S000000493 and SGD:S000005141, respectively). The coding sequence of amino acids 1–253 of the Sup35NM3 mutant 20Franzmann et al.2018 was cloned from a plasmid kindly provided by Titus Franzmann. The coding sequence of human FUS low-complexity region (amino acids 2–214 from UniProt entry Q6IBQ5) was codon optimized for yeast expression, and synthesized with the rest of the Ede1 sequence by Synbio Technologies (New Jersey). The sequence coding for amino acids 335–931 of D. melanogaster kinesin-1 was cloned from Addgene plasmid K980 (#129762), a gift from William Hancock. The sequence coding for amino acids 358–797 of Homo sapiens kinesin-5 was cloned from Addgene plasmid mCherry-Kinesin11-N-18 (#55067), a gift from Michael Davidson.

Plasmids used in Figure 9 were based on pRS426-GFP-2×PH(PLCδ) 74Stefan et al.2002, a kind gift from Scott Emr. The Ede1366-900 coding sequence was inserted into this plasmid after the last GFP codon using ligation-independent cloning. GFP-2×PH(PLCδ) and GFP-Ede1366-900-2×PH(PLCδ) were then subcloned to a pRS416-based plasmid p416-GPD under the control of TDH3 promoter 53Mumberg et al.1995 using BamHI and SalI restriction sites. GFP-Ede1366-900-2×PH(PLCδ) was also subcloned into p416-CYC1, p416-ADH1, and p416-TEF1 53Mumberg et al.1995 using the same restriction sites.

The sequence of the msGFP2 fluorophore 80Valbuena et al.2020 was cloned into a PFA6a-based tagging plasmid from Addgene plasmid #135301, a kind gift from Benjamin Glick.

Live cell imaging

Yeast cells were grown to OD600 between 0.3 and 0.8 at 24°C in low-fluorescence synthetic drop-out medium lacking tryptophan, or tryptophan and uracil if required for plasmid maintenance. Cells were attached to cover slips coated with 1 mg ml−1 concanavalin A.

Widefield microscopy

Widefield micrographs were obtained on an Olympus IX81 widefield microscope equipped with a 100×/NA1.45 objective and an ORCA-ER CCD camera (Hamamatsu), using an X-CITE 120 PC (EXFO) metal halide lamp as the illumination source. The excitation and emission light when imaging EGFP- and mCherry-tagged proteins were filtered through the U-MGFPHQ and U-MRFPHQ filter sets (Olympus). The 3D stacks were acquired with 0.2 μm vertical spacing. The microscope was controlled using the MetaMorph software (Molecular Dynamics).

Total internal reflection fluorescence microscopy

All TIRF movies were recorded on an Olympus IX83 widefield microscope equipped with a 150×/NA1.45 objective and an ImageEM X2 EM-CCD camera (Hamamatsu) under the control of the VisiView software (Visitron Systems). The 488 nm and 561 nm laser lines were used for illumination of GFP- and mCherry-tagged proteins. Excitation and emission were filtered using a TRF89902 405/488/561/647 nm quad-band filter set (Chroma). Laser angles were controlled by iLas2 (Roper Scientific).

Fluorescence recovery after photobleaching

Bleaching of Ede1-EGFP in endocytic condensates (Figure 2A and B) was performed using a custom-built set-up that focuses a 488-nm laser beam at the sample plane, on the Olympus IX81 widefield microscope described above. The diameter of the bleach spot was approximately 0.5 μm.

Bleaching of unperturbed endocytic sites (Figure 2C) was performed with a 405 nm laser line controlled by the iLas2 targeting system during simultaneous excitation with 488 nm and 561 nm lasers in TIRF mode on the Olympus IX83 microscope described above. The emission light was collected through a Gemini beam splitter (Hamamatsu) equipped with a Di03-R488/561-t1 dichroic, and FF03-525/50-25 and FF01-630/92-25 emission filters (Semrock).

Spinning disk microscopy

Spinning disk confocal imaging (Figure 3) was performed in the Photonic Bioimaging Center at the University of Geneva using a Nikon Eclipse Ti1 microscope equipped with a CSU-W1 spinning disk (Yokogawa) using a 100×/NA1.49 objective, an sCMOS Prime 95B camera (Photometrics), and 488 nm and 561 nm lasers as the illumination source.

Induction of protein expression

For the induction of expression from GALS promoter during live-cell imaging (Figure 3), cells were thawed and grown for several days on Synthetic Complete medium agar plates with 2% galactose as the sole carbon source. The cells were then cultured overnight in a low-fluorescence synthetic drop-out liquid medium with no tryptophan and 2% raffinose as the sole carbon source. The cells were diluted in the morning into the same medium with 2% glucose as the carbon source. The cells were attached to cover slips as described above. Finally, the carbon source in the medium was switched to 2% galactose before the start of the imaging.

Image and data analysis

All code used in this study is available as a single repository at https://github.com/matkozak/KozakAndKaksonen2022, (copy archived at swh:1:rev:5441acf218619f2b03d90633613cccc373c6fe8a; 32Kozak2022). General image analysis was performed using the Fiji distribution of ImageJ 69Schindelin et al.201266Rueden et al.2017. All display images were corrected for background fluorescence using the rolling ball algorithm of ImageJ, and movies were corrected for photobleaching using a custo m ImageJ macro. Plots and statistical analyses were generated using R.

FRAP experiments

FRAP experiments performed on Ede1-EGFP condensates were analyzed according to 62Phair et al.2004. Mean fluorescence values were measured from regions of interest representing the background, the cell, and the condensate. A custom-written R script (available in the article repository) was used to subtract background fluorescence, correct for photobleaching and normalize the values between 0 (corrected fluorescence immediately after photobleaching) and 1 (mean corrected fluorescence of 5 s before photobleaching). The recovery curves of individual experiments were aligned to bleach time and averaged. Condensates that showed lateral or axial movement during the acquisition were manually excluded from the averaging. The average was fitted to a single exponential equation from which the mobile fraction and recovery half-time were calculated.

For FRAP experiments performed on native endocytic sites, the background fluorescence was first subtracted from the TIRF images using the ImageJ rolling ball algorithm. EGFP and mCherry fluorescence of single endocytic patches were measured within a circle with a radius of three pixels around the patch centroid position. A custom-written R script was used to calculate the fluorescence recovery much in the same way as for the FRAP of condensates, but no further corrections were made for background signal or imaging-induced photobleaching. To calculate average recovery, we manually selected only events in which Abp1-mCherry signal peaked at least 60 s after bleach time to exclude the effect of Ede1 disassembly at the end of endocytic events.

Patch numbers and lifetimes

For estimating the number of patches per membrane area, we analyzed single nonbudding cells. The patches were thresholded and counted using a custom Python script available in the article repository. We estimated the cell surface area by measuring the area of the cross-section from maximum intensity projection and multiplying it by 4, under the assumption that an unbudded yeast cell is approximately spherical. For estimating patch lifetimes, we tracked endocytic events using ParticleTracker from the MOSAIC suite 68Sbalzarini and Koumoutsakos2005 and multiplied trajectory length by the frame rate.

Cytosolic and total cellular intensity

To obtain cytosolic intensity of Ede1-EGFP, 5 × 5-pixel square regions away from the condensates and vacuoles were manually measured in ImageJ. To measure total cellular intensity, individual cells were cropped in ImageJ. A custom Python script was applied to the cropped cells to generate masks based on Rvs167-mCherry fluorescence, and subsequently measure the Ede1-EGFP signal intensity in the masked region.

Cell classification

To calculate the percentages of Ede1-EGFP condensates colocalizing with mCherry puncta in Figure 4, cells containing Ede1-EGFP condensates in Figure 5 and cells showing different GFP-EDE1366-900−2×PH localization patterns in Figure 9, single cells were cropped from imaging fields based on a neutral signal (GFP in the case of Figure 4 and brightfield image for Figures 5 and !number(9)). Next, an ImageJ macro was used to display random images from the dataset and the experimenter would assess the presence of the tested phenotype with no knowledge of which strain was being analyzed.

Western blotting

A 300 μl of ice-cold trichloroacetic acid was added to 5 ml of exponentially growing yeast cultures. The cells were pelleted by centrifugation, washed with cold acetone, and dried in a vacuum concentrator. The pellets were resuspended in 100 μl of urea buffer (25 mM Tris-HCl pH 6.8, 6 M urea, 1% SDS) and homogenized by shaking with 200 μl of glass beads. The samples were heated at 95°C for 5 min, mixed with 100 μl 2× SDS loading buffer and centrifuged at 16, 000× g for 5 min.

The samples were subjected to electrophoresis on 4–20% Precast Protein Gels (Bio-Rad) and transferred onto a nitrocellulose membrane using an iBlot2 device (ThermoFischer Scientific). The membranes were blocked for 30 min with 5% bovine serum albumin in PBS-Tween and incubated with primary antibodies overnight at 4°C. The membranes were washed in PBS-Tween, incubated with fluorescent secondary antibodies for 1 hr and washed in PBS-Tween. The fluorescence was measured on an Odyssey scanner (LI-COR Biosciences).

Antibodies

Ede1 constructs were detected using an anti-GFP mouse monoclonal antibody (ab291, Abcam) at 1/2000 dilution, and an anti-Hog1 rabbit polyclonal antibody (sc-9079, Santa Cruz Biotechnology) at 1/1000 dilution was used as a loading control. Donkey antimouse IRDye 680 and anti-rabbit IRDye 800 secondary antibodies (926–68072 and 926–32213 respectively, LI-COR Biosciences) were used at a 1/10,000 dilution.

References

    A systematic survey identifies prions and illuminates sequence features of prionogenic proteins137Cell146158
    Phosphatidylinositol-(4,5)-bisphosphate regulates clathrin-coated pit initiation, stabilization, and size22Molecular Biology of the Cell25882600
    Endocytic sites mature by continuous bending and remodeling of the clathrin coat348Science (New York, N.Y.)13691372
    Compositional Control of Phase-Separated Cellular Bodies166Cell651663
    Biomolecular condensates: organizers of cellular biochemistry18Nature Reviews. Molecular Cell Biology285298
    Phase transitions of multivalent proteins can promote clustering of membrane receptors3eLife
    Endocytic proteins with prion-like domains form viscoelastic condensates that enable membrane remodeling118PNAS
    Quantification of cytosolic interactions identifies Ede1 oligomers as key organizers of endocytosis10Molecular Systems Biology
    The initiation of clathrin-mediated endocytosis is mechanistically highly flexible24Current Biology548554
    Analysis of yeast endocytic site formation and maturation through a regulatory transition point23Molecular Biology of the Cell657668
    Regulation of Transmembrane Signaling by Phase Separation48Annual Review of Biophysics465494
    Dynamic instability of clathrin assembly provides proofreading control for endocytosis218The Journal of Cell Biology32003211
    The first five seconds in the life of a clathrin-coated pit150Cell495507
    Assessing the tendency of fluorescent proteins to oligomerize under physiologic conditions13Traffic (Copenhagen, Denmark)643649
    Parallel dimers and anti-parallel tetramers formed by epidermal growth factor receptor pathway substrate clone 15272The Journal of Biological Chemistry3343033434
    Natively unfolded domains in endocytosis: hooks, lines and linkers5EMBO Reports10461052
    Liquid-like protein interactions catalyse assembly of endocytic vesicles23Nature Cell Biology366376
    Evidence that the stalk of Drosophila kinesin heavy chain is an alpha-helical coiled coil116The Journal of Cell Biology957965
    SLAC, a complex between Sla1 and Las17, regulates actin polymerization during clathrin-mediated endocytosis23Molecular Biology of the Cell42564272
    Phase separation of a yeast prion protein promotes cellular fitness359Science (New York, N.Y.)
    Prion-like low-complexity sequences: Key regulators of protein solubility and phase behavior294The Journal of Biological Chemistry71287136
    A novel EH domain protein of Saccharomyces cerevisiae, Ede1p, involved in endocytosis113Journal of Cell Science33093319
    Review series: From uncertain beginnings: initiation mechanisms of clathrin-mediated endocytosis203The Journal of Cell Biology717725
    Pathways and mechanisms of endocytic recycling10Nature Reviews Molecular Cell Biology597608
    The anchor-away technique: rapid, conditional establishment of yeast mutant phenotypes31Molecular Cell925932
    FCHo proteins are nucleators of clathrin-mediated endocytosis328Science (New York, N.Y.)12811284
    A versatile toolbox for PCR-based tagging of yeast genes: new fluorescent proteins, more markers and promoter substitution cassettes21Yeast (Chichester, England)947962
    A modular design for the clathrin- and actin-mediated endocytosis machinery123Cell305320
    Mechanisms of clathrin-mediated endocytosis19Nature Reviews. Molecular Cell Biology313326
    Cell-free formation of RNA granules: low complexity sequence domains form dynamic fibers within hydrogels149Cell753767
    Molecular structure, function, and dynamics of clathrin-mediated membrane traffic6Cold Spring Harbor Perspectives in Biology
    https://archive.softwareheritage.org/swh:1:dir:2782495bab5f864ac750b4113d6c7cd87ac229d2;origin=https://github.com/matkozak/KozakAndKaksonen2022;visit=swh:1:snp:ca9e80e9e888a6be19acdd04741370e9e7e85441;anchor=swh:1:rev:5441acf218619f2b03d90633613cccc373c6fe8a
    Promiscuous interactions and protein disaggregases determine the material state of stress-inducible RNP granules4eLife
    Plasma membrane reshaping during endocytosis is revealed by time-resolved electron tomography150Cell508520
    Single-molecule turnover dynamics of actin and membrane coat proteins in clathrin-mediated endocytosis8eLife
    PLAAC: A web and command-line application to identify proteins with prion-like amino acid composition30Bioinformatics (Oxford, England)25012502
    Liquid droplet formation by HP1α suggests a role for phase separation in heterochromatin547Nature236240
    Modeling vesicle traffic reveals unexpected consequences for Cdc42p-mediated polarity establishment21Current Biology184194
    Phase transitions in the assembly of multivalent signalling proteins483Nature336340
    SLIC: A method for sequence- and ligation-independent cloningGene Synthesis: Methods and Protocols, Methods in Molecular Biology5159Humana Press
    Local clustering of transferrin receptors promotes clathrin-coated pit initiation191The Journal of Cell Biology13811393
    Clathrin-mediated endocytosis in budding yeast at a glance129Journal of Cell Science15311536
    Selection and stabilization of endocytic sites by Ede1, a yeast functional homologue of human Eps1528Molecular Biology of the Cell567575
    A framework for understanding the functions of biomolecular condensates across scales22Nature Reviews. Molecular Cell Biology215235
    Transient Fcho1/2Eps15/RAP-2 Nanoclusters Prime the AP-2 Clathrin Adaptor for Cargo Binding37Developmental Cell428443
    Interaction between Epsin/Yap180 adaptors and the scaffolds Ede1/Pan1 is required for endocytosis19Molecular Biology of the Cell29362948
    Protein disorder, prion propensities, and self-organizing macromolecular collectives1834Biochimica et Biophysica Acta918931
    Evaluating phase separation in live cells: diagnosis, caveats, and functional consequences33Genes & Development16191634
    IUPred2A: context-dependent prediction of protein disorder as a function of redox state and protein binding46Nucleic Acids Research
    Phase separation by low complexity domains promotes stress granule assembly and drives pathological fibrillization163Cell123133
    Phosphorylation of the FUS low-complexity domain disrupts phase separation, aggregation, and toxicity36The EMBO Journal29512967
    Regulatable promoters of Saccharomyces cerevisiae: comparison of transcriptional activity and their use for heterologous expression22Nucleic Acids Research57675768
    Yeast vectors for the controlled expression of heterologous proteins in different genetic backgrounds156Gene119122
    Systematic Nanoscale Analysis of Endocytosis Links Efficient Vesicle Formation to Patterned Actin Nucleation174Cell884896
    MobiDB-lite: fast and highly specific consensus prediction of intrinsic disorder in proteins33Bioinformatics (Oxford, England)14021404
    Phase transition of a disordered nuage protein generates environmentally responsive membraneless organelles57Molecular Cell936947
    Natively unfolded nucleoporins gate protein diffusion across the nuclear pore complex129Cell8396
    A Liquid-to-Solid Phase Transition of the ALS Protein FUS Accelerated by Disease Mutation162Cell10661077
    Spatial regulation of clathrin-mediated endocytosis through position-dependent site maturation219The Journal of Cell Biology
    Casein kinase 1 promotes initiation of clathrin-mediated endocytosis32Developmental Cell231240
    Evidence for and against Liquid-Liquid Phase Separation in the Nucleus5Non-Coding RNA
    Measurement of dynamic protein binding to chromatin in vivo, using photobleaching microscopy375Methods in Enzymology393414
    Visualizing the functional architecture of the endocytic machinery4eLife
    Syp1 is a conserved endocytic adaptor that contains domains involved in cargo selection and membrane tubulation28The EMBO Journal31033116
    Composition-dependent thermodynamics of intracellular phase separation581Nature209214
    ImageJ2: ImageJ for the next generation of scientific image data18BMC Bioinformatics
    Coactivator condensation at super-enhancers links phase separation and gene control361Science (New York, N.Y.)
    Feature point tracking and trajectory analysis for video imaging in cell biology151Journal of Structural Biology182195
    Fiji: an open-source platform for biological-image analysis9Nature Methods676682
    Structural basis for the assembly of the mitotic motor Kinesin-5 into bipolar tetramers3eLife
    Quantitative analysis of the mechanism of endocytic actin patch assembly and disassembly in fission yeast21Molecular Biology of the Cell28942904
    Molecular basis for coupling the plasma membrane to the actin cytoskeleton during clathrin-mediated endocytosis109PNAS
    Endocytic proteins are partitioned at the edge of the clathrin lattice in mammalian cells19Nature Cell Biology352361
    The yeast synaptojanin-like proteins control the cellular distribution of phosphatidylinositol (4,5)-bisphosphate13Molecular Biology of the Cell542557
    Early-Arriving Syp1p and Ede1p Function in Endocytic Site Placement and Formation in Budding Yeast20Molecular Biology of the Cell46404651
    Endocytic internalization in budding yeast requires coordinated actin nucleation and myosin motor activity11Developmental Cell3346
    A Pan1/End3/Sla1 complex links Arp2/3-mediated actin assembly to sites of clathrin-mediated endocytosis26Molecular Biology of the Cell38413856
    A high precision survey of the molecular dynamics of mammalian clathrin-mediated endocytosis9PLOS Biology
    Eps15 is constitutively oligomerized due to homophilic interaction of its coiled-coil region272The Journal of Biological Chemistry1541315418
    A photostable monomeric superfolder green fluorescent protein21Traffic (Copenhagen, Denmark)534544
    A Selective Autophagy Pathway for Phase-Separated Endocytic Protein Deposits80Molecular Cell764778
    The Centrosome Is a Selective Condensate that Nucleates Microtubules by Concentrating Tubulin169Cell10661077
    Clathrin exchange during clathrin-mediated endocytosis155The Journal of Cell Biology291300
    Proteome-wide signatures of function in highly diverged intrinsically disordered regions8eLife