ABSTRACT
Nuclear envelope (NE) expansion must be controlled to maintain nuclear shape and function. The nuclear membrane expands massively during closed mitosis, enabling chromosome segregation within an intact NE. Phosphatidic acid (PA) and diacylglycerol (DG) can both serve as biosynthetic precursors for membrane lipid synthesis. How they are regulated in time and space and what the implications are of changes in their flux for mitotic fidelity are largely unknown. Using genetically encoded PA and DG probes, we show that DG is depleted from the inner nuclear membrane during mitosis in the fission yeast Schizosaccharomyces pombe, but PA does not accumulate, indicating that it is rerouted to membrane synthesis. We demonstrate that DG-to-PA conversion catalyzed by the diacylglycerol kinase Dgk1 (also known as Ptp4) and direct glycerophospholipid synthesis from DG by diacylglycerol cholinephosphotransferase/ethanolaminephosphotransferase Ept1 reinforce NE expansion. We conclude that DG consumption through both the de novo pathway and the Kennedy pathway fuels a spike in glycerophospholipid biosynthesis, controlling NE expansion and, ultimately, mitotic fidelity.
INTRODUCTION
The double-membrane nuclear envelope (NE) is a hallmark of eukaryotic cells. It consists of the outer nuclear membrane that is continuous with the endoplasmic reticulum (ER), and the inner nuclear membrane (INM) facing the nucleoplasm. The NE is perforated by nuclear pore complexes that regulate the transport and exchange of macromolecules between the cytoplasm and the nucleus. Tight control of nuclear shape and size is essential for proper cell function, and aberrant regulation of nuclear morphology has been implicated in aging and disease (Cantwell and Dey, 2022; Webster et al., 2009). The spherical shape of the interphase nucleus minimizes the amount of membrane material required to make the NE. Changes in lipid biosynthesis, nuclear trafficking and/or INM–chromatin attachments can all lead to deviation from this simplest of shapes (Arnone et al., 2013; Cantwell and Dey, 2022; Webster et al., 2009).
Closed mitosis of lower eukaryotes, in which the nucleocytoplasmic integrity is maintained throughout chromosome partitioning, represents a fascinating example of highly regulated NE expansion that is not scaled to nuclear volume increase (Zhang and Oliferenko, 2013). The model fission yeast Schizosaccharomyces pombe undergoes closed mitosis, whereas its relative, Schizosaccharomyces japonicus, breaks and reforms the NE (Yam et al., 2011). We have previously linked this divergence to differences in nuclear membrane management. As nuclear volume remains constant in closed mitosis, the surface area of the mother nucleus must increase to allow for the formation of two daughter nuclei. The nuclear membrane expands dramatically in S. pombe but not in S. japonicus, necessitating NE breakdown in S. japonicus (Yam et al., 2011). Massive expansion of the NE during mitosis in S. pombe is driven by the CDK1 (Cdc2)-dependent inactivation of an evolutionarily conserved regulator of phosphatidic acid (PA) flux, lipin (Ned1). This does not occur in S. japonicus, likely due to differences in trans-regulation of lipin phosphorylation (Makarova et al., 2016). Lipin inactivation upon mitotic entry presumably causes a block in the conversion of PA to diacylglycerol (DG) and a diversion of PA towards the production of glycerophospholipids (GPLs), thus triggering NE expansion required for closed nuclear division (Makarova et al., 2016).
The subcellular distribution of the lipin substrate PA and its product DG at the NE and the possible implications of changes in their distribution for NE expansion and mitotic fidelity remain unknown. Current models of how the regulation of PA flux controls NE growth are largely based on lipidomics analyses of extracts from cells with a constitutively (e.g. gene deletion or mutation) or slowly (e.g. overexpression) modified PA regulation network (Han et al., 2008; Santos-Rosa et al., 2005). Yet, in such cases, cells settle at metabolic steady states different from the that of the wild type, limiting the insight into dynamic events occurring in mitosis.
In the budding yeast, which also undergoes closed mitosis, lipin deficiency leads to PA accumulation (Han et al., 2008). It was proposed that such an increase in PA might lead to changes in the biophysical properties of the nuclear membrane (Han et al., 2008). This, together with PA-dependent transcriptional upregulation of GPL biosynthetic genes (Loewen et al., 2004), could lead to increased nuclear membrane biogenesis. Yet, this regulation circuitry is not conserved in fission yeast (Rhind et al., 2011; Rutherford et al., 2022), and it is also not clear whether lipin inactivation leads to increased PA levels in S. pombe. To gain an insight into the mechanisms underlying mitotic NE expansion, it is imperative to distinguish between different subcellular pools of PA and DG, and to probe their fates in normal mitosis and upon deregulation of the relevant biosynthetic pathways. Here, we address these questions by using a combination of imaging, genetic and cell physiological approaches in S. pombe and the related fission yeast S. japonicus.
RESULTS
Development of genetically encoded biosensors to probe subcellular distribution of PA and DG in S. pombe
We built a set of genetically encoded PA and DG biosensors, extending the designs used previously to visualize these lipid classes in Saccharomyces cerevisiae (Romanauska and Köhler, 2018). Our PA and DG sensors were expressed from the medium-strength S. pombe cdc15 promoter and tagged at the C-terminus with a GST tag and a fluorophore (Fig. 1A). Appending or omitting the SV40 nuclear localization sequence (NLS) around the lipid-binding module allowed us to probe lipid enrichment in the nucleus or the cytoplasm.
Lipid biosensors for probing subcellular distribution of PA and DG in fission yeast. (A) Diagram depicting the design of lipid sensors used in this study. The lipid sensors consist of a lipid-binding module (S. cerevisiae Opi1 PA-binding amphipathic helix and R. norvegicus PKCβ DG-binding domain) under the control of the S. pombe cdc15 promoter, with a C-terminal GST tag and a fluorophore. In the NLS versions of the sensors, the lipid-binding module is flanked by two SV40 NLS motifs. (B) Single-plane spinning-disk confocal images of S. pombe expressing the PA and NLS–PA sensors show their enrichment at LDs and the INM, respectively. LDs were stained with BODIPY 558/568. (C) Mutations in the PA-binding domain abolished normal localization of the PA sensors. Note that the presence of a weak bipartite NLS in the mutated PA-binding domain resulted in its nuclear re-localization in the absence of the SV40 NLS. (D) Single-plane spinning-disk confocal images of S. pombe cells expressing the DG and NLS–DG sensors show their enrichment at the cortex, LDs and the INM. (E) Mutations predicted to abolish DG-binding led to re-localization of the sensors. The DG-binding mutant sensor became largely cytoplasmic, whereas the NLS version was enriched in the nucleoplasm. (F) Overexpression of the DG kinase Dgk1 increased the intensity of the PA sensor and decreased the intensity of the DG sensor at the cortex. No significant changes were observed upon overexpression of the catalytically inactive D119A mutant of Dgk1. The wild-type and the mutant Dgk1 proteins were expressed under the control of the thiamine-repressible nmt1 promoter. (G,H) Scatter plots showing the quantitation of these phenotypes. n=3 biological repeats, 20 cells each. Each biological replicate is color coded and the medians are indicated. P-values were derived using non-parametric heteroscedastic two-tailed unpaired t-test for all data points. In panels B–F, the insets represent magnification of the areas indicated by the dotted lines. Scale bars: 5 µm.
Lipid biosensors for probing subcellular distribution of PA and DG in fission yeast. (A) Diagram depicting the design of lipid sensors used in this study. The lipid sensors consist of a lipid-binding module (S. cerevisiae Opi1 PA-binding amphipathic helix and R. norvegicus PKCβ DG-binding domain) under the control of the S. pombe cdc15 promoter, with a C-terminal GST tag and a fluorophore. In the NLS versions of the sensors, the lipid-binding module is flanked by two SV40 NLS motifs. (B) Single-plane spinning-disk confocal images of S. pombe expressing the PA and NLS–PA sensors show their enrichment at LDs and the INM, respectively. LDs were stained with BODIPY 558/568. (C) Mutations in the PA-binding domain abolished normal localization of the PA sensors. Note that the presence of a weak bipartite NLS in the mutated PA-binding domain resulted in its nuclear re-localization in the absence of the SV40 NLS. (D) Single-plane spinning-disk confocal images of S. pombe cells expressing the DG and NLS–DG sensors show their enrichment at the cortex, LDs and the INM. (E) Mutations predicted to abolish DG-binding led to re-localization of the sensors. The DG-binding mutant sensor became largely cytoplasmic, whereas the NLS version was enriched in the nucleoplasm. (F) Overexpression of the DG kinase Dgk1 increased the intensity of the PA sensor and decreased the intensity of the DG sensor at the cortex. No significant changes were observed upon overexpression of the catalytically inactive D119A mutant of Dgk1. The wild-type and the mutant Dgk1 proteins were expressed under the control of the thiamine-repressible nmt1 promoter. (G,H) Scatter plots showing the quantitation of these phenotypes. n=3 biological repeats, 20 cells each. Each biological replicate is color coded and the medians are indicated. P-values were derived using non-parametric heteroscedastic two-tailed unpaired t-test for all data points. In panels B–F, the insets represent magnification of the areas indicated by the dotted lines. Scale bars: 5 µm.
The PA recognition module was engineered from the Q2 domain of the S. cerevisiae transcriptional factor Opi1 (Loewen et al., 2004), with further optimization enhancing its PA-binding properties (Hofbauer et al., 2018). Specifically, the PA sensor used throughout this study was constructed by removing the endogenous NLS within the codon-optimized Q2 helix, introducing a G120W point mutation predicted to increase PA-binding selectivity, and duplicating the helix to stabilize it at the membrane (Fig. S1A,B; see Table S1 for the DNA sequence). In exponentially growing S. pombe, the cytoplasmic PA sensor was enriched mildly at the cellular cortex and strongly in distinct puncta, confirmed to be lipid droplets (LDs), as visualized by staining cells with the neutral lipid dye BODIPY 558/568 (Fig. 1B, upper panel). The NLS version of the sensor was enriched at the INM (Fig. 1B, lower panel). Confirming sensor specificity, mutation of amino acid residues known to mediate PA binding (Loewen et al., 2004) abrogated the specific localization patterns (Fig. 1C). Note that these mutations yield a weak bipartite NLS and result in nuclear import even in the case of the cytoplasmic construct.
The DG biosensor was developed based on the duplicated codon-optimized Rattus norvegicus protein kinase C (PKCβ C1a/b) DG-binding domain (Lucic et al., 2016) (Fig. 1D; Fig. S1C,D; see Table S1 for the DNA sequence). The DG sensor was strongly enriched in the INM and the cellular cortex, as well as in some but not all LDs (Fig. 1D). The DG-binding specificity was validated by using DG-binding-deficient versions of the sensor containing two point mutations (Q63E in the C1a domain and Q128E in the C1b domain) (Lucic et al., 2016). The loss of DG binding resulted in the redistribution of constructs to the nucleoplasm and/or cytoplasm, depending on the presence of the NLS (Fig. 1E). The cortical localization of the DG biosensor remained unaffected in the absence of ER–plasma membrane attachments in scs2Δ scs22Δ genetic background (Zhang et al., 2012), indicating that DG is likely enriched at the plasma membrane but not at the peripheral ER in S. pombe (Fig. S1E).
Of note, significant vacuolar autofluorescence that is common in fission yeast cells expressing low-to-moderate amounts of GFP-tagged proteins can be largely corrected by utilizing a recently developed Spectral Autofluorescence Image Correction By Regression (SAIBR) protocol (Rodrigues et al., 2022). Processing images of cells expressing PA and DG sensors by SAIBR highlighted specific subcellular enrichments of these constructs (Fig. S1F).
We validated both sensors further by analyzing their distribution under conditions in which the relative abundance of PA and DG was expected to deviate from that of the wild type. Overexpression of the DG kinase Dgk1 (also known as Ptp4), which catalyzes the reverse reaction to lipin by phosphorylating DG to produce PA, from the strong thiamine-repressible nmt1 promoter resulted in an increase in PA sensor signal at the cellular cortex (Fig. 1F,G), with a concomitant decrease in the cortical intensity of the DG sensor (Fig. 1F,H). These phenotypes were not observed in cells overexpressing the D119A catalytically inactive mutant of Dgk1 (Han et al., 2008).
The phosphatase Spo7–Nem1 positively regulates lipin activity, and the loss of either subunit of the phosphatase complex phenocopies the lack of lipin (Makarova et al., 2016; Santos-Rosa et al., 2005; Siniossoglou et al., 1998). Liquid chromatography electrospray ionization tandem mass spectrometry (LC-MS/MS) analyses of the total cellular extracts of S. pombe lipin pathway mutants revealed a slight decrease in PA levels in the spo7Δ mutant and no significant changes in the nem1Δ and lipin ned1Δ mutants (Fig. 2A). As expected, lipin pathway deficiency resulted in a significant decrease in cellular DG and triacylglycerol (TG) levels (Fig. 2A,B). To investigate these changes at the subcellular level, we analyzed the distribution of our sensors in cells lacking Spo7–Nem1 phosphatase activity. Despite the overall similar or lower cellular abundance, the PA sensor was enriched at the INM of both spo7Δ and nem1Δ mutant cells grown in the rich yeast extract with supplements (YES) medium (Fig. 2C,D). The DG sensor levels at the INM were correspondingly lower (Fig. 2E,F).
Lipid biosensors are useful to detect changes in PA and DG distribution at the subcellular level. (A,B) Relative abundance of PA and DG (A), and TG (B) in S. pombe lipin pathway mutants. n=3 biological replicates. The results are presented as mean±s.d. P-values were estimated using non-parametric heteroscedastic two-tailed unpaired t-test. (C–F) PA sensor levels at the INM were mildly increased (C,D) and DG sensor levels were decreased (E,F) when lipin was inactive in the absence of its activator phosphatase Spo7–Nem1. Note that whereas NLS–DG–GFP sensor enrichment at the INM was decreased overall, we detected bright patches at the NE of some nuclei. Nucleoporin Nup40–mCherry marks the nuclear boundary. The insets represent magnification of the areas indicated by the dotted lines. Scale bars: 5 µm. n=50 cells (D,F). Box-and-whisker plots were plotted according to the Tukey method. P-values were derived from two-tailed unpaired t-test. A.U., arbitrary units.
Lipid biosensors are useful to detect changes in PA and DG distribution at the subcellular level. (A,B) Relative abundance of PA and DG (A), and TG (B) in S. pombe lipin pathway mutants. n=3 biological replicates. The results are presented as mean±s.d. P-values were estimated using non-parametric heteroscedastic two-tailed unpaired t-test. (C–F) PA sensor levels at the INM were mildly increased (C,D) and DG sensor levels were decreased (E,F) when lipin was inactive in the absence of its activator phosphatase Spo7–Nem1. Note that whereas NLS–DG–GFP sensor enrichment at the INM was decreased overall, we detected bright patches at the NE of some nuclei. Nucleoporin Nup40–mCherry marks the nuclear boundary. The insets represent magnification of the areas indicated by the dotted lines. Scale bars: 5 µm. n=50 cells (D,F). Box-and-whisker plots were plotted according to the Tukey method. P-values were derived from two-tailed unpaired t-test. A.U., arbitrary units.
Of note, the NLS–DG sensor was occasionally enriched at the pronounced NE ‘flares’ formed in S. pombe spo7Δ and nem1Δ mutants (Fig. 2E, insets). The DG kinase Dgk1 was not excluded from the NLS–DG sensor-rich domains (Fig. S2A, insets). These patches were found throughout the nuclear periphery, i.e. they were not associated with the nucleolus (Fig. S2B) (Witkin et al., 2012). These domains could be also induced by the major reduction in cellular DG through the overexpression of Dgk1 (Fig. S2C), or by increasing the absolute amount of the NLS–DG sensor by expressing it under the strong tdh1 promoter (Fig. S2D). We reason that such local enrichments of the NLS–DG sensor might result from the imbalance between the sensor and its target lipid, and care should be taken when designing the experiments aimed at evaluating the suborganellar distribution of DG, particularly in the mutant backgrounds affecting its abundance.
Taken together, our results indicate that the newly developed genetically encoded sensors are useful as reporters to probe changes in PA and DG distribution in live fission yeast cells.
DG at the INM is depleted during closed mitosis of S. pombe
Time-lapse imaging showed that the fluorescence intensity of the NLS–PA sensor at the INM decreased mildly during mitosis (Fig. 3A–C). This suggested that lipin inactivation did not cause a spike in PA abundance at the NE, as previously proposed for budding yeast (Han et al., 2008). Presumably, upon lipin inactivation, all excess PA was diverted for GPL production, which is necessary for NE expansion (Makarova et al., 2016). The mild decrease could reflect the dilution effect due to the expansion in the surface area of the NE (reviewed by Zhang and Oliferenko, 2013).
The decay of the DG sensor at the inner nuclear membrane during closed mitosis in S. pombe does not result in higher PA sensor enrichment at the NE. (A) Time-lapse microscopy montage of mitotic S. pombe expressing the NLS–PA sensor. (B) Average PA sensor intensity at the INM of S. pombe mildly decreased over the course of mitosis compared to interphase cells. (C) Linear regression analysis showed that the rate of the NLS–PA sensor decay in mitotic cells was significantly different compared to interphase cells. (D) Time-lapse microscopy montage of mitotic S. pombe expressing the NLS–DG sensor. (E) The NLS–DG sensor intensity at the INM rapidly decayed as mitosis progressed. (F) The rate of decrease of the NLS–DG sensor at the INM of mitotic S. pombe was significantly higher than in interphase. (G) Time-lapse microscopy montage of the related fission yeast S. japonicus that does not expand the NE during mitosis. The NE breakage occurred between 5 and 6 min. (H,I) DG levels at the INM remained relatively constant throughout the semi-open mitosis of S. japonicus (H), quantified in I. For A,D,G, the maximum projections of three z-slices are shown. Scale bars: 5 µm. Nucleoporins Nup40–mCherry (A,D) and Nup85–mCherry (G) mark the nuclear boundary. For B,C,E,F,H,I, n=10 cells. Box-and-whisker plots were plotted according to the Tukey method. P-values were derived from two-tailed unpaired t-test.
The decay of the DG sensor at the inner nuclear membrane during closed mitosis in S. pombe does not result in higher PA sensor enrichment at the NE. (A) Time-lapse microscopy montage of mitotic S. pombe expressing the NLS–PA sensor. (B) Average PA sensor intensity at the INM of S. pombe mildly decreased over the course of mitosis compared to interphase cells. (C) Linear regression analysis showed that the rate of the NLS–PA sensor decay in mitotic cells was significantly different compared to interphase cells. (D) Time-lapse microscopy montage of mitotic S. pombe expressing the NLS–DG sensor. (E) The NLS–DG sensor intensity at the INM rapidly decayed as mitosis progressed. (F) The rate of decrease of the NLS–DG sensor at the INM of mitotic S. pombe was significantly higher than in interphase. (G) Time-lapse microscopy montage of the related fission yeast S. japonicus that does not expand the NE during mitosis. The NE breakage occurred between 5 and 6 min. (H,I) DG levels at the INM remained relatively constant throughout the semi-open mitosis of S. japonicus (H), quantified in I. For A,D,G, the maximum projections of three z-slices are shown. Scale bars: 5 µm. Nucleoporins Nup40–mCherry (A,D) and Nup85–mCherry (G) mark the nuclear boundary. For B,C,E,F,H,I, n=10 cells. Box-and-whisker plots were plotted according to the Tukey method. P-values were derived from two-tailed unpaired t-test.
Strikingly, quantitation of the NLS–DG sensor intensity over time revealed that it decreased at approximately twice the rate of the NLS–PA sensor intensity during S. pombe mitosis (Fig. 3D–F). To probe whether such a decay could be due to changes in sensor abundance, we analyzed NLS–DG sensor levels in S. pombe cultures undergoing synchronous mitosis by western blotting. We detected no significant differences in the sensor protein levels over the course of mitosis, indicating that a drop in its fluorescence intensity at the NE was likely reflective of a decrease in DG at that membrane (Fig. S3A–C).
Such a depletion in DG could be due, at least in part, to CDK1-dependent lipin inactivation. In fact, the NLS–DG sensor levels remained constant at the INM of mitotic S. japonicus (Fig. 3G–I), as lipin regulation is not tied to the cell cycle in this organism (Makarova et al., 2016). Interestingly, the subcellular enrichment of both PA and DG sensors in S. japonicus was different from that in S. pombe. The PA sensor was only very weakly enriched on the NE and on the cortex of S. japonicus (Fig. S3D), whereas the DG sensor was detected at the cortex and on both sides of the NE (Fig. 3G; Fig. S3E). DG levels at the outer nuclear membrane of S. japonicus remained stable during mitosis (Fig. S3E–G), similar to its behavior at the INM (Fig. 3G–I). Similar to S. pombe, strong autofluorescence from the vacuoles and mitochondria could be reduced in S. japonicus cells expressing the PA and DG sensors by SAIBR correction (Fig. S3H).
Given that S. japonicus does not expand the NE during mitosis, our results suggest that the depletion of DG at the INM of S. pombe reflects rerouting of the lipid biosynthetic flow towards GPL production, which is necessary for nuclear membrane expansion.
DG-to-PA conversion by Dgk1 contributes to NE expansion during closed mitosis
We next wondered whether the DG kinase Dgk1 contributed to the decrease in DG levels at the INM during mitosis of S. pombe. Dgk1 was previously identified in S. cerevisiae (Han et al., 2008) and the dgk1Δ mutant was recently shown to display a mild under-expanded ER phenotype (Papagiannidis et al., 2021), hinting at its role in controlling ER growth. The S. pombe ortholog is encoded by the dgk1 (also known as ptp4) gene (SPBC3D6.05). However, the role of Dgk1 during mitosis and its effect on the subcellular distribution of PA and DG in any system remained unknown.
Consistent with previous observations in S. cerevisiae (Han et al., 2008), the loss of Dgk1 in S. pombe was able to rescue the NE and ER expansion phenotypes of the lipin pathway mutants (Fig. S4A,B), indicating the overall circuitry conservation. The loss of Dgk1 did not rescue the large LD phenotype in the lipin-deficient mutants (Fig. S4C) but restored the number of LDs to wild-type levels (Fig. S4D).
While analyzing dgk1Δ S. pombe strains, we consistently noticed the presence of larger diploids in exponentially growing cultures of haploid cells. Indeed, flow cytometric analysis showed progressive accumulation of diploids in S. pombe dgk1Δ cultures (Fig. 4A). In dgk1Δ S. japonicus, diploidization was not detected (Fig. 4B). We reasoned that if the loss of Dgk1 resulted in a deficiency in NE expansion, the elongating spindle would experience compressive stress along its longitudinal direction, and would buckle and break, similarly to cells unable to produce fatty acids (Yam et al., 2011) or deficient in mitotic lipin inactivation (Makarova et al., 2016). The spindle collapse would result in the failure to divide the nucleus and diploidization. Time-lapse sequences of mitotic wild-type S. pombe co-expressing α-tubulin (mCherry–Atb2) and the NLS–DG sensor tagged with GFP (NLS–DG–GFP), growing in rich medium, showed that the elongating anaphase spindles remained largely straight and DG sensor intensity at the INM decayed (Fig. 4C). Interestingly, the dgk1Δ S. pombe cells under the same conditions exhibited a range of mitotic phenotypes. Severe buckling of the mitotic spindle was observed in 31% of cells, resulting in bow-shaped nuclear intermediates during anaphase (Fig. 4D). It took longer for these nuclei to divide, and they often formed daughter nuclei of unequal sizes (Fig. S4E). In 11% of dgk1Δ S. pombe cells, the NE failed to expand completely, and the mitotic spindle buckled and broke within the confines of the NE (Fig. 4E). Such cells could end up containing the nuclei with twice the genetic content and, therefore, were the likely source of a diploid cell population in dgk1Δ cultures. The remaining 57% of dgk1Δ S. pombe had no obvious mitotic phenotypes and divided normally (Fig. S4F,G).
The DG kinase Dgk1 is essential for maintaining mitotic fidelity in S. pombe. (A,B) Diploids accumulated upon prolonged culturing of dgk1Δ S. pombe (A), as shown by flow cytometric analysis. In S. japonicus (B), the loss of Dgk1 did not result in diploidization. Representative graphs from at least three biological repeats are shown. (C–E) Representative time-lapse spinning-disk microscopy montages showed a range of anaphase phenotypes in dgk1Δ S. pombe. In the wild type (C), the spindle remained relatively straight throughout mitosis, indicating that it was not constrained by the NE. In dgk1Δ S. pombe, 31% of cells exhibited anaphase spindle buckling with eventual relaxation (D), whereas 11% broke spindles following their compression within the confines of the NE (E). The maximum projections of three z-slices are shown, and the number of cells analyzed is indicated at the bottom right of each time-lapse sequence. α-tubulin mCherry-Atb2 marks the mitotic spindle. Scale bars: 5 µm. (F–H) The decrease in the NLS–DG sensor intensity in cells exhibiting spindle buckling (F) and breaking (G) phenotypes was compared against that of the wild-type. The rate of decrease of the NLS–DG sensor intensity at the INM is predictive of the severity of the mitotic phenotype in dgk1Δ S. pombe (H). n=8 cells. Box-and-whisker plots were plotted according to the Tukey method. P-values were derived from two-tailed unpaired t-test.
The DG kinase Dgk1 is essential for maintaining mitotic fidelity in S. pombe. (A,B) Diploids accumulated upon prolonged culturing of dgk1Δ S. pombe (A), as shown by flow cytometric analysis. In S. japonicus (B), the loss of Dgk1 did not result in diploidization. Representative graphs from at least three biological repeats are shown. (C–E) Representative time-lapse spinning-disk microscopy montages showed a range of anaphase phenotypes in dgk1Δ S. pombe. In the wild type (C), the spindle remained relatively straight throughout mitosis, indicating that it was not constrained by the NE. In dgk1Δ S. pombe, 31% of cells exhibited anaphase spindle buckling with eventual relaxation (D), whereas 11% broke spindles following their compression within the confines of the NE (E). The maximum projections of three z-slices are shown, and the number of cells analyzed is indicated at the bottom right of each time-lapse sequence. α-tubulin mCherry-Atb2 marks the mitotic spindle. Scale bars: 5 µm. (F–H) The decrease in the NLS–DG sensor intensity in cells exhibiting spindle buckling (F) and breaking (G) phenotypes was compared against that of the wild-type. The rate of decrease of the NLS–DG sensor intensity at the INM is predictive of the severity of the mitotic phenotype in dgk1Δ S. pombe (H). n=8 cells. Box-and-whisker plots were plotted according to the Tukey method. P-values were derived from two-tailed unpaired t-test.
Strikingly, we observed a correlation between the rates of DG sensor decay at the INM and the severity of the mitotic phenotypes in dgk1Δ S. pombe cells (Fig. 4F–H; Fig. S4G). Whereas DG sensor intensity at the INM decayed similarly to the wild type in the normally dividing dgk1Δ mutants (Fig. S4G), it decayed to a lesser degree in cells with buckling spindles (Fig. 4F) and remained constant in those cells that failed to expand the NE (Fig. 4G).
Thus, both the activity of Dgk1 as well as concurrent inactivation of lipin contribute to the decrease in DG levels at the NE during closed mitosis in S. pombe. Importantly, this indicates that Dgk1 activity promotes NE expansion required for closed mitosis. The fact that DG might decay at different rates in mitotic dgk1Δ cells, yielding distinct nuclear division phenotypes, suggests an underlying metabolic heterogeneity even within isogenic cellular populations.
The balance between DG and PA, but not the elevated levels of PA, is the hallmark of NE expansion in S. pombe
Given that the lack of Dgk1 counteracted excessive NE–ER expansion in the lipin pathway mutants (Fig. S4A,B), we wondered whether the loss of lipin activity resulting in a drop in DG and an increase in PA levels at the NE (Fig. 2D,F) might in turn alleviate insufficient mitotic NE expansion observed in dgk1Δ mutant cells. Intriguingly, the spo7Δ dgk1Δ double mutants displayed normal mitosis and a decay of DG sensor intensity at the NE similar to that of the wild type (Fig. 5A,B; see also Fig. S5A,B for nem1Δ dgk1Δ data). Moreover, we observed a significant reduction in diploidization in the cultures of the double mutants compared to cultures of the dgk1Δ strain (Fig. 5C; Fig. S5C). Inactivating the lipin pathway might constitutively channel PA towards GPL synthesis through the cytidine diphosphate (CDP)–DG pathway, rendering the Dgk1-dependent conversion of DG to PA during mitosis unnecessary.
Changes in the PA-to-DG sensor ratio at the INM are consistent with NE expansion in S. pombe. (A) Time-lapse microscopy of the spo7Δ dgk1Δ double mutant undergoing mitosis showed a decay of NLS–DG sensor intensity at the INM, similar to that of the wild type. The maximum projections of three z-slices are shown. Nucleoporin Nup40–mCherry marks the nuclear boundary. Scale bars: 5 µm. (B) Quantification of the rate of decay. n=10 cells. Box-and-whisker plots were plotted according to the Tukey method. P-values were derived from two-tailed unpaired t-test. (C) Flow cytometry showed that the spo7Δ dgk1Δ double mutant had a lower proportion of diploids in the population compared to dgk1Δ cells. (D,E) Quantification of the levels of DG (D) and PA (E) sensors at the INM in the indicated genotypes of interphase S. pombe are shown as scatter plots. n=3 biological repeats, 50 cells each. Each biological replicate is color coded and the medians indicated. P-values were derived using non-parametric heteroscedastic two-tailed unpaired t-test for all data points. A.U., arbitrary units.
Changes in the PA-to-DG sensor ratio at the INM are consistent with NE expansion in S. pombe. (A) Time-lapse microscopy of the spo7Δ dgk1Δ double mutant undergoing mitosis showed a decay of NLS–DG sensor intensity at the INM, similar to that of the wild type. The maximum projections of three z-slices are shown. Nucleoporin Nup40–mCherry marks the nuclear boundary. Scale bars: 5 µm. (B) Quantification of the rate of decay. n=10 cells. Box-and-whisker plots were plotted according to the Tukey method. P-values were derived from two-tailed unpaired t-test. (C) Flow cytometry showed that the spo7Δ dgk1Δ double mutant had a lower proportion of diploids in the population compared to dgk1Δ cells. (D,E) Quantification of the levels of DG (D) and PA (E) sensors at the INM in the indicated genotypes of interphase S. pombe are shown as scatter plots. n=3 biological repeats, 50 cells each. Each biological replicate is color coded and the medians indicated. P-values were derived using non-parametric heteroscedastic two-tailed unpaired t-test for all data points. A.U., arbitrary units.
As expected, Dgk1 deficiency resulted in a mild increase in DG sensor intensity at the NE (Fig. 5D). Perhaps less intuitively, the PA sensor intensity increased substantially in this genetic background, similarly to lipin pathway mutants (Fig. 5E). This suggests that increasing PA levels alone is not sufficient for NE expansion (Han et al., 2008). Instead, we observed changes in the relative abundance of PA and DG sensors at the NE. Indeed, the increase in PA sensor intensity at the INM in the lipin pathway mutants was accompanied by a drop in DG sensor intensity, whereas this was not the case in cells lacking Dgk1. The double mutants lacking both lipin and Dgk1 activities largely restored the relative abundances of PA and DG sensors at the INM (Fig. 5D,E), and showed normal NE morphology (Fig. S4A,B). We confirmed that both PA and DG biosensors were expressed at comparable levels in cells of all genotypes (Fig. S5D–G).
Taken together, our results suggest that the regulation of DG levels at the INM is important in controlling NE expansion. Yet, it was unclear how DG at the NE was removed in cells that lacked Dgk1 but did not display any mitotic phenotype.
The Kennedy pathway contributes strongly to GPL synthesis required for NE expansion during closed mitosis
In addition to Dgk1-dependent rerouting of DG towards PA for de novo GPL synthesis, DG might be also used directly by the precursor-dependent Kennedy pathway to produce GPLs (Carman and Han, 2009b). Our observations suggest that the contribution of the Kennedy pathway is critical for proper nuclear membrane expansion during mitosis. First, dgk1Δ mutant cultures grown in the chemically defined medium (Edinburgh minimal medium or EMM) that does not contain choline and ethanolamine (Cho/Etn) precursors required for the Kennedy pathway exhibited a high incidence of cells with unequally dividing nuclei and complete failure in nuclear division (‘cut’ phenotype) (Fig. 6A–C). Many cells in these cultures became non-viable (Fig. 6D). Second, whereas the PA sensor was highly abundant at the INM of dgk1Δ cells grown in the presence of the Kennedy pathway (the rich YES medium or EMM+Cho/Etn), we did not observe any enrichment in the absence of the Kennedy pathway precursors (Fig. 6E, compare with Fig. 5E). This suggests that rerouting of DG to GPL biosynthesis via the Kennedy pathway might decrease the consumption of PA via the CDP–DG and/or lipin pathways. Differential activity of the DG-consuming Kennedy pathway within a population of cells could be an underlying cause for the heterogeneity in cellular phenotypes associated with the loss of Dgk1.
The Kennedy pathway of GPL synthesis collaborates with Dgk1 in sustaining NE expansion required for closed mitosis in S. pombe. (A–C) In the absence of the Kennedy pathway precursors, dgk1Δ cells exhibited a substantial increase in uneven nuclear divisions (quantitation in B) and mitotic failure (‘cut’ phenotype; quantitation in C). The single-plane spinning-disk confocal images of S. pombe cells co-expressing the NLS–DG–GFP sensor and Nup40–mCherry are shown. n=4 biological repeats, ≥10 septated cells each. Scale bars: 5 µm. (D) High levels of cell death were detected by propidium iodide staining in S. pombe dgk1Δ cells cultured in the absence of the Kennedy pathway precursors. n=5 biological repeats, ≥50 cells each. Box-and-whisker plots were plotted according to the Tukey method. P-values were derived from two-tailed unpaired t-test. (E) Quantification of the levels of PA at the INM in wild-type and dgk1Δ cells grown in the indicated media. n=3 biological repeats, 50 cells each. Each biological replicate is color coded and the medians are indicated. P-values were derived using non-parametric heteroscedastic two-tailed unpaired t-test for all data points. (F) A representative single-plane spinning-disk confocal image shows Ept1–GFP localization at the NE of S. pombe. (G) Diploids accumulated upon culturing of the Kennedy pathway-deficient ept1Δ mutant strain of S. pombe, as shown by flow cytometric analysis. (H) Tetrad analysis revealed the semi-lethality of the ept1Δ dgk1Δ double mutant. A representative tetrad grown on rich YES agar on days 5 and 12 after dissection is shown. n=54 tetrads. (I) Relative growth rates of the indicated genotypes in rich YES medium. n=5 biological repeats. Error bars represent mean±s.d.
The Kennedy pathway of GPL synthesis collaborates with Dgk1 in sustaining NE expansion required for closed mitosis in S. pombe. (A–C) In the absence of the Kennedy pathway precursors, dgk1Δ cells exhibited a substantial increase in uneven nuclear divisions (quantitation in B) and mitotic failure (‘cut’ phenotype; quantitation in C). The single-plane spinning-disk confocal images of S. pombe cells co-expressing the NLS–DG–GFP sensor and Nup40–mCherry are shown. n=4 biological repeats, ≥10 septated cells each. Scale bars: 5 µm. (D) High levels of cell death were detected by propidium iodide staining in S. pombe dgk1Δ cells cultured in the absence of the Kennedy pathway precursors. n=5 biological repeats, ≥50 cells each. Box-and-whisker plots were plotted according to the Tukey method. P-values were derived from two-tailed unpaired t-test. (E) Quantification of the levels of PA at the INM in wild-type and dgk1Δ cells grown in the indicated media. n=3 biological repeats, 50 cells each. Each biological replicate is color coded and the medians are indicated. P-values were derived using non-parametric heteroscedastic two-tailed unpaired t-test for all data points. (F) A representative single-plane spinning-disk confocal image shows Ept1–GFP localization at the NE of S. pombe. (G) Diploids accumulated upon culturing of the Kennedy pathway-deficient ept1Δ mutant strain of S. pombe, as shown by flow cytometric analysis. (H) Tetrad analysis revealed the semi-lethality of the ept1Δ dgk1Δ double mutant. A representative tetrad grown on rich YES agar on days 5 and 12 after dissection is shown. n=54 tetrads. (I) Relative growth rates of the indicated genotypes in rich YES medium. n=5 biological repeats. Error bars represent mean±s.d.
Underscoring the notion that the de novo and the Kennedy pathways both contribute to GPL synthesis, we detected a decrease in NE–ER proliferation in EMM, compared to either EMM containing Cho/Etn precursors or rich YES medium (Fig. S6A,B). As DG synthesis is decreased in the lipin-deficient mutants, the Kennedy pathway presumably becomes sufficient for consumption of DG at the NE, explaining the mitotic DG decay observed in spo7Δ dgk1Δ double mutants (Fig. 5A,B; see also Fig. S5A,B for nem1Δ dgk1Δ data).
Ept1 (SPAC22A12.10) is the DG cholinephosphotransferase/ethanolaminephosphotransferase that catalyzes the final reaction of the Kennedy pathway in fission yeast. Ept1–GFP is enriched at the NE of S. pombe throughout the cell cycle, suggesting that DG-dependent GPL biosynthesis might occur at this location (Fig. 6F). In order to test directly whether the Kennedy pathway contributes to NE expansion, we generated a mutant strain lacking Ept1. The ept1Δ mutants exhibited striking diploidization with a majority of the population, having double the genetic content within 3 days of consecutive culture (Fig. 6G). The introduction of the tagged α-tubulin (mCherry-Atb2) triggered a much higher incidence of mitotic failure, with mutant cultures containing diploids even after a short period of culturing (Fig. S6C). Importantly, the ept1Δ dgk1Δ double mutant phenotype was sub-lethal, as observed through genetic crosses and phase-contrast microscopy (Fig. 6H; Fig. S6D). The rare ept1Δ dgk1Δ double mutants that we were able to recover exhibited severe delay in growth, compared to the wild type and the single ept1Δ and dgk1Δ mutants (Fig. 6I).
In summary, our results show that controlling the utilization of DG at the INM is important for the regulation of NE expansion in S. pombe. Both the Dgk1-dependent DG-to-PA conversion and the Kennedy pathway collaborate to channel DG-to-GPL production required for a mitotic spike in NE growth.
DISCUSSION
Our results show that NE expansion in S. pombe appears to require a pronounced drop in DG levels at the INM. This DG depletion is the result of several enzymatic reactions. Lipin, which synthesizes DG from PA, is inactivated at mitotic entry by CDK1 (Carman and Han, 2009a; Makarova et al., 2016), reducing DG inflow. At the same time, DG can be either phosphorylated by Dgk1 to produce PA (Han et al., 2008; Kwiatek et al., 2020) or directly used for GPL biosynthesis through the Kennedy pathway (Carman and Han, 2009b). As we did not observe PA sensor levels spiking at the NE at mitosis, it is likely that the Dgk1-produced PA was being used up for the synthesis of membrane lipids through the de novo CDP–DG route (Han et al., 2008; Kwiatek et al., 2020) (Fig. 7).
Hypothesis on how PA and DG flux during NE expansion is controlled in S. pombe. In wild-type S. pombe, the inactivation of lipin Ned1 by CDK1, together with Dgk1 activity and the Kennedy pathway, channel DG towards GPL biosynthesis. In the absence of Dgk1, insufficient depletion of DG results in reduced membrane availability and/or changes in the biophysical properties of the membrane. When Dgk1 is inactive, cells rely on Kennedy pathway-dependent DG consumption to produce GPLs. Cell-to-cell heterogeneity in the activity of the Kennedy pathway results in mitotic phenotypes, such as uneven nuclear division or chromosome segregation defects (‘cut’ phenotype). Many cells end up with twice the genetic content, forming diploids. G3P, glycerol-3-phosphate.
Hypothesis on how PA and DG flux during NE expansion is controlled in S. pombe. In wild-type S. pombe, the inactivation of lipin Ned1 by CDK1, together with Dgk1 activity and the Kennedy pathway, channel DG towards GPL biosynthesis. In the absence of Dgk1, insufficient depletion of DG results in reduced membrane availability and/or changes in the biophysical properties of the membrane. When Dgk1 is inactive, cells rely on Kennedy pathway-dependent DG consumption to produce GPLs. Cell-to-cell heterogeneity in the activity of the Kennedy pathway results in mitotic phenotypes, such as uneven nuclear division or chromosome segregation defects (‘cut’ phenotype). Many cells end up with twice the genetic content, forming diploids. G3P, glycerol-3-phosphate.
Lipidomics analysis of the S. cerevisiae lipin mutant showed an increase in cellular PA levels (Han et al., 2008). Such an increase in PA, a conically shaped lipid capable of generating negative membrane curvature (Zhukovsky et al., 2019), was inferred to modify nuclear membrane properties and lead to the transcriptional upregulation of GPL biosynthetic genes through the Opi1p–Ino2p–Ino4p regulation circuitry (Carman and Henry, 2007; Han et al., 2008), although the latter is not conserved in fission yeasts (Rhind et al., 2011; Rutherford et al., 2022). Of note, PA levels do not seem to increase on a global scale in the lipin-deficient S. pombe cells (Fig. 2A), nor at the NE during mitosis, when lipin function is inhibited by CDK1 (Fig. 3A,B) (Makarova et al., 2016). They also do not increase in the budding yeast spo7Δ and nem1Δ mutants, despite the block to lipin activity (Papagiannidis et al., 2021).
Similarly to PA, DG induces negative spontaneous curvature and was shown to promote the fusion of biological membranes, NE assembly and LD formation (Choudhary et al., 2018; Chung et al., 2018; Domart et al., 2012; Dumas et al., 2010; Miner et al., 2017). It is possible that, in addition to fueling membrane lipid synthesis, depletion of DG or changes in the PA-to-DG ratio might modify the biophysical properties of the NE, priming it for expansion. However, we favor the possibility that rather than having functional implications for NE remodeling, the PA-to-DG ratio at the NE simply reflects the utilization of these lipids by different biosynthetic pathways. Both lipids are precursors for the biosynthesis of other lipid species. GPLs are mainly synthesized through two pathways, the de novo CDP–DG-dependent route, which uses PA, and the precursor-dependent Kennedy pathway, which uses DG (Gibellini and Smith, 2010; McMaster and Bell, 1994). In addition, DG can be used for the production of the storage lipid TG (Carman and Han, 2009b; Carman and Henry, 2007; Holic et al., 2020). The differential contribution of these biosynthetic pathways might manifest as changes in the steady-state PA-to-DG ratio.
The range of mitotic phenotypes within isogenic populations of Dgk1-deficient cells grown in rich media suggests an underlying metabolic heterogeneity. This might arise due to a number of reasons, e.g., differences in the Kennedy pathway precursor uptake (Gibellini and Smith, 2010), the differential activities of biosynthetic pathways (Stewart-Ornstein et al., 2012) or stochasticity in gene expression (Kaern et al., 2005; Raj and van Oudenaarden, 2008). The phenotypic heterogeneity largely collapses when the Kennedy pathway is disabled either by withdrawal of precursors or the loss of Ept1 (Fig. 6), suggesting that the Kennedy pathway activity might indeed show cell-to-cell variability.
Both Dgk1 and the Kennedy pathway have roles in NE expansion that enables closed mitosis and, hence, the maintenance of mitotic fidelity in S. pombe (Figs 4–7). Interestingly, a functional genomics screen has linked the loss of Dgk1 to meiotic chromosome segregation defects (Blyth et al., 2018). These phenotypes were speculated to be a result of global deregulation in lipid synthesis (Holic et al., 2020). Our data suggest that one possible alternative explanation could be membrane limitation during meiotic remodeling of the NE.
The bulk of NE expansion during closed mitosis likely occurs in situ or at the ER domain close to the NE (Fig. 6F; Fig. S2A). Our results might have implications for mitosis in other organisms. Mitotic strategies in nature span from completely closed to completely open mitosis, in which the NE breaks down in prophase and reforms around the segregated chromosomes upon mitotic exit (Makarova and Oliferenko, 2016). Recent evidence highlights the importance of regulating GPL synthesis for NE reformation and other dynamic events at the NE in metazoans (Bahmanyar and Schlieker, 2020), and it would be of interest to address whether the Kennedy pathway or the control of DG-to-PA conversion by DG kinases is also critical for these processes.
Advancements in the system-level analyses of lipids have allowed greater scrutiny of their roles in various cellular processes (Han, 2016; Kofeler et al., 2021). Yet, these approaches struggle with analyzing lipid dynamics with fine spatiotemporal resolution. Complementing lipidomics and genetics approaches with a set of newly developed lipid sensors for S. pombe has allowed us to track changes in the spatiotemporal distribution of PA and DG with subcellular resolution and to formulate and test a set of hypotheses on the roles of PA-to-DG interconversion during mitotic NE remodeling. We believe that these sensors will become an effective and widely used tool for the study of dynamic membrane processes in fission yeasts.
MATERIALS AND METHODS
Strains, media and molecular biology methods
S. pombe and S. japonicus strains used in this study are listed in Table S2. Standard fission yeast media and methods were used (Aoki et al., 2010; Moreno et al., 1991). All experiments were performed using the rich non-defined YES medium with the exception of lipidomics, overexpression of Dgk1 from the thiamine-repressible nmt1 promoter (Maundrell, 1990) and the Kennedy pathway-related experiments shown in Fig. 6A–E and Fig. S6A,B, which were performed using chemically defined EMM. All strains were grown at 30°C unless otherwise specified, in temperature-controlled 200 rpm shaking incubators. Mating was performed on SPA medium and spores dissected on YES agar plates. Homozygous S. pombe diploids were generated by ade6-M210/ade6-M216 heteroallelic complementation (Ekwall and Thon, 2017). Choline and ethanolamine (Sigma-Aldrich) were added to EMM at final concentrations of 1 mM.
Molecular genetic manipulations were performed using either plasmid-based (Keeney and Boeke, 1994) or PCR-based (Bähler et al., 1998) homologous recombination. Lipid sensors were integrated at the ura4 locus.
Overexpression of Dgk1 and its catalytic mutant was carried out by integration of pREP1-Dgk1 (our laboratory collection) or its mutated version into the leu1 locus of S. pombe. For overexpression experiments, cells were pre-grown in EMM containing 5 µg ml−1 thiamine at 30°C for 18 h. These cells were then collected by centrifugation and washed thrice in EMM. Cells were then diluted into separate flasks of EMM with or without thiamine and grown for 20 h prior to imaging.
For the experiments in Fig. 6A–E and Fig. S6A,B, cells were precultured in YES medium at 30°C for 18 h. They were then diluted into separate flasks at an OD595 of 0.001 and grown at 30°C for 18 h in the appropriate medium (YES, EMM or EMM+Cho/Etn) prior to imaging the following day. Propidium iodide (Sigma-Aldrich) was used at a final concentration of 0.1 µg ml−1.
Growth rates were measured using a VICTOR Nivo Multimode Microplate Reader (PerkinElmer). Cells were initially precultured in YES medium at 30°C for 18 h. They were then diluted to an OD595 of 0.05 and 200 µl of each culture was seeded into five wells as technical repeats and grown at 30°C for 50 h without shaking. Absorbance readings were obtained every 10 min. Doubling time was calculated using the formula T=, where T is the doubling time, x is the final OD595 and x0 is the initial OD595. Results were expressed as doubling time of the experimental strain relative to that of the wild type.
Lipid sensor design
All lipid sensors were cloned into the pJK210-based plasmid backbone (our laboratory collection). DNA sequences of the lipid biosensors were synthesized using the GeneWiz FragmentGENE service. S. pombe promoters (promtdh1, promcdc15) were inserted between the KpnI and ApaI restriction sites, which were followed by the biosensor sequence between the XhoI and EcoRI restriction sites. This was immediately followed by the GST tag located between the PacI and BamHI restriction sites in the S. pombe sensors, followed by the fluorescent tag (GFP or mCherry) between the BamHI and NotI restriction sites, with a stop codon located before the NotI restriction site (Fig. 1A). The NLS version of the lipid biosensors contained two SV40 NLS motif appended to both ends of the lipid sensor module.
The PA lipid-binding module consisted of the S. cerevisiae Opi1 residues 111–189 to exclude the endogenous NLS signal located at residues 109–112, a duplication of the PA-binding amphipathic helices (residues 114–131) and a point mutation introduced at G120W (Fig. S1A; see Table S1 for the DNA sequence) (Hofbauer et al., 2018). The DG biosensor consisted of the duplicated DG-binding domain of R. norvegicus PKCβ (residues 31–158) (Fig. S1C; see Table S1 for the DNA sequence) (Lucic et al., 2016). The PA-binding mutant was constructed with the following mutations: L124R, Y127A, L129R, M131A, I133R, K136A, K137A and R138A (Loewen et al., 2004). The DG-binding mutant contained point mutations at Q63E and Q128E (Lucic et al., 2016). The sequences of all sensor variants were codon optimized for S. pombe.
Lipidomics
Fission yeast cultures for lipidomics were grown in the defined EMM medium with all supplements (adenine, uracil, histidine and leucine) at 30°C. Cells were collected by filtration and snap frozen in liquid nitrogen. Lipid extraction was performed using a two-phase chloroform-methanol extraction protocol (Ejsing et al., 2009). Cells were first disrupted in 200 µl of water using a Bead Ruptor 12 (OMNI International) with ceramic beads at 4.75 speed setting for six cycles of 30 s at 4°C. Then, 100 µl of the diluted lysate was transferred to fresh Eppendorf tubes and 900 µl of chloroform-methanol 17:1 (v/v) was added. The mixture was left shaking at 4°C for 2 h at 200 rpm in an Eppendorf thermomixer shaker. Following that, the mixture was centrifuged at 10,000 g at 4°C for 2 min. The lower phase was transferred to a fresh Eppendorf tube and dried in a Thermo Fisher Scientific Savant SpeedVac, whereas a second chloroform-methanol 2:1 (v/v) extraction was performed with the upper phase. The second extraction was combined with the first and dried in the SpeedVac before resuspension in 200 µl of chloroform-methanol 1:1 (v/v) containing internal standards (Table S3) and stored at −80°C until further use. Batch quality-control (QC) samples were made by pooling and mixing 10 µl of each sample lysate. These were aliquoted into eight Eppendorf tubes for lipid extraction as described earlier. Blanks were prepared by using 100 µl of water in place of lysates for lipid extraction. All lipid extracts from samples were mixed and aliquoted into nine vials as technical QC samples. The technical QC samples were diluted with chloroform-methanol 1:1 (v/v) to prepare 100, 50, 25, 12.5 and 6.25% diluted samples to assess instrument response linearity.
Initial profiling of S. pombe lipids with isotopic correction were performed by Lipid Data Analyzer (https://www.lipidmaps.org/resources/tools/10?task=4.5) using the wild-type strain SO2865. Multiple reaction monitoring lists for each lipid class of interest were constructed from the data obtained from the initial profiling experiments and used for targeted LC-MS lipidomics (Tables S4 and S5). All LC-MS/MS was performed on an Agilent 6495A QqQ mass spectrometer connected to a 1290 series chromatographic system with electrospray ionization for lipid ionization. All sample injection volumes were set at 2 µl. The spray voltage and nozzle voltage were set at 3500 V and 500 V, respectively. The drying gas and sheath gas temperatures were maintained at 200°C and 250°C, respectively, with flow rates both set at 12 l min−1. The nebulizer setting was 25 psi. Following instrument stabilization with 15 injections of a QC sample, instrument stability was monitored by injection of a QC sample and blank sample every five sample injections.
Chromatographic separation of the GPLs phosphatidylinositol, phosphatidylcholine, phosphatidylethanolamine and phosphatidylglycerol was performed using a gradient elution on an ethylene-bridged hybrid column with solvent A containing 1:1 (v/v) acetonitrile/H2O with 10 mM ammonium formate, and solvent B containing 19:1 (v/v) acetonitrile/H2O with 10 mM ammonium formate. The flow rate was set at 0.4 ml min−1, with solvent B set at 40% upon injection and increasing to 100% in 7 min. This value was retained for 2 min, decreased back to 40% in 1 min, and then retained there until the end of the gradient by 14 min. The eluent was directed to the electron-spray ionization source of the mass spectrometer and operated in positive ionization.
Chromatographic separation of PA and phosphatidylserine was performed using a different liquid phase (Triebl et al., 2014). Mobile phase A was deionized water containing 10 mM ammonium formate and 0.5% formic acid. Mobile phase B was 2-propanol/acetonitrile 5:2 (v/v) containing 10 mM ammonium formate and 0.5% formic acid. Gradient elution began at 5% solvent A with a linear increase to 50% over 12 min, 50% solvent A was then held for 3 min and, lastly, the column was re-equilibrated for 15 min.
Chromatographic separation of the neutral lipids DG and TG was performed using reversed-phase liquid chromatography on an Agilent rapid-resolution HD Zorbax Eclipse-C18 column. The mobile phases A consisted of 3:2 (v/v) H2O/acetonitrile with 10 mM ammonium formate and B consisted of 1:9 (v/v) acetonitrile/isopropanol with 10 mM ammonium formate. The flow rate was set at 0.4 ml min−1, with solvent B set at 20% upon injection and increasing to 60% in 2 min, then increasing to 100% over 12 min, and held there for a further 2 min. This was then decreased to 20% for the next 2 min.
Raw data were processed with MassHunter QqQ Quantitative software (version B.08, Agilent, https://www.agilent.com/en/product/software-informatics/mass-spectrometry-software/data-analysis/quantitative-analysis). Areas under the curve of the chromatogram peaks for each transition were measured and exported to Excel. Normalized peak areas were calculated by dividing the peak areas of the analyte with the corresponding internal standard. Relative abundances were obtained by multiplying the normalized peak areas with the molar concentration of the corresponding internal standard.
All lipidomics experiments were performed with three biological replicates and results are presented as mean±s.d. Statistical analysis was performed using non-parametric heteroscedastic two-tailed unpaired t-test. All solvents used were LC-MS grade purchased from Sigma-Aldrich or Thermo Fisher Scientific. Annotation of lipid classes and species was done according to the classification system previously described (Liebisch et al., 2013).
Fluorescence-activated cell sorting
S. pombe cells were grown at 30°C in YES medium and 107 cells were collected. Cultures were diluted every 24 h to an OD595 of 0.0005, and samples at an OD595 of 0.5 were typically collected over the course of 3 days. Cell pellets were resuspended in 1 ml of 70% ice-cold ethanol and kept at 4°C until required. For staining of cells for fluorescence-activated cell sorting (FACS) analysis, cells were first centrifuged at 1000 g for 2 min at 4°C and the supernatant removed. The pellet was washed once in 1 ml of 50 mM sodium citrate and resuspended in 500 µl of 50 mM sodium citrate containing 0.1 mg ml−1 RNaseA (Thermo Fisher Scientific). This was incubated at 37°C for 18 h with gentle shaking. Then, 500 µl of 50 mM sodium citrate containing either 2 µM SYTOX green (Thermo Fisher Scientific) for S. pombe samples or 8 µg ml−1 propidium iodide (Sigma-Aldrich) for S. japonicus samples was added, and the suspension sonicated briefly to remove cell doublets before FACS analysis. All FACS experiments were performed on a LSRFortessa Cell Analyzer (BD Biosciences) at the Crick Flow Cytometry Science Technology Platform, with at least three biological replicates.
Quantification and statistical analyses
Image processing and quantifications were performed in Fiji (Schindelin et al., 2012). Circularity of the nucleus was calculated using the formula θ= where θ is the circularity, A is the area and P is the perimeter.
Quantification of PA and DG sensor intensities at the INM or cortex of mitotic cells was performed using the segmented line function in Fiji with a line width of 5 pixels every 1 min of the time-lapse images. Lipid sensor decay at the NE of mitotic cells was normalized by taking the ratio of the sensor intensity at the NE to the cortex. The lines and shaded area represent the mean±s.d. The rate of sensor decay at the NE was calculated by taking the gradient via linear regression of the sensor decay at the NE in individual cells. For quantification of PA and DG at the INM of interphase cells (Fig. 1G,H, Fig. 2D,F, Fig. 5D,E, Fig. 6E), the segmented line function in Fiji with a line width of 5 pixels was used. To account for variation in our sensor expression levels, the average whole-cell fluorescence intensity was also measured using the polygon function in Fiji, and values above 1.5 interquartile range were discarded.
Scatter plots were used to present the data quantifying NLS–PA and NLS–DG sensor levels in cell populations (Fig. 1G,H, Fig. 5D,E, Fig. 6E). Each biological replicate was color coded and the medians indicated. P-values were estimated using non-parametric heteroscedastic two-tailed unpaired t-test for all data points. Box-and-whisker plots were created using the Tukey method in Prism 7 (GraphPad software). Lipid sensor protein levels (Fig. S3C, Fig. S5E,G) were compared using parametric one-way ANOVA with Brown–Forsythe post hoc test (n=3 biological repeats).
Image acquisition
The working concentration of all BODIPY-based dyes (Thermo Fisher Scientific) used was 1 µM. Time-lapse microscopy of cells undergoing mitosis was performed using cells grown on agar pads (Pemberton, 2014). All images were obtained using a Yokogawa CSU-X1 spinning disk confocal system mounted on the Eclipse Ti-E inverted microscope with a Nikon CFI Plan Apo Lambda 100× NA 1.45 oil objective, 600 series SS 488 nm, SS 561 nm lasers and Andor iXon Ultra U3-888-BV monochrome EMCCD camera controlled by Andor IQ3 or Fusion software. The temperature was maintained at 30°C, unless otherwise specified, using an Okolab Cage Incubator. Maximum-intensity projections of three z-slices of 0.5 µm step-size are shown in time-lapse montages. For autofluorescence correction, images of PA and DG biosensors in wild-type S. pombe and S. japonicus (Fig. S1F, Fig. S3H) were processed using the SAIBR Fiji plug-in (Rodrigues et al., 2022).
Cell cycle synchronization and western blotting
S. pombe cdc25-22 temperature-sensitive mutant cells expressing the NLS–DG sensor were grown for 18 h at the permissive temperature of 24°C in YES medium until the OD595 reached 0.2–0.3. They were next shifted to the restrictive temperature of 36°C for 3.5 h. Cells were released from the G2/M block by rapid cooling in an ice-water bath to 24°C and transferred to a shaking incubator set at 24°C. Five OD595 equivalents of cells were collected every 20 min and snap frozen in liquid nitrogen for western blotting. Additional 5 ml aliquots of cells were fixed in 70% ethanol, and cell cycle synchronization was verified by measuring the number of binucleated and septated cells by staining DNA with 2 µM SYTOX green (Fig. S3A–C; n=3 biological repeats). For comparing expression levels of PA and DG sensors in cells of different genetic backgrounds (Fig. S5D–G), five OD595 equivalents of cells were collected at an OD595 of 0.4–0.6.
After collection, cells were pelleted for 1 min at 1000 g. Cell pellets were snap frozen in liquid nitrogen and stored at −20°C until further use. For lysis of cells, cell pellets were first washed in 1 ml of ice-cold water and then resuspended in 1 ml of water containing 10% trichloroacetic acid for 1 h on ice. Cells were next pelleted for 10 min at 18,213 g at 4°C and the supernatant was removed. The cell pellets were then washed once with 1 ml of ice-cold acetone, pelleted again and dried in a SpeedVac for 2 min at room temperature. Cell pellets were then resuspended in 300 µl of lysis buffer (50 mM Tris-HCl pH 8.0, 1 mM EDTA and 1% SDS) and transferred into 2 ml Lysing Matrix Y tubes (MP Biomedicals). Lysis was performed using a MP Biomedicals FastPrep24 5G homogenizer for 15 s at 6.5 m s−1 twice in a 4°C cold room, with 150 s interval between cycles. Lysates were separated from the beads by puncturing the bottom of the Lysing Matrix Y tubes using a hot needle, followed by spinning down at 500 g for 3 min at 4°C in a 15 ml Falcon tube containing a 1 ml pipette tip. Cellular debris were cleared by further centrifugation in a clean Eppendorf tube for 1000 g for 5 min at 4°C. Lysates were boiled at 95°C for 5 min, followed by addition of 100 µl of 4× NuPAGE LDS sample buffer (Invitrogen) containing 10% β-mercaptoethanol, before heating for a further 10 min at 65°C. Then, 20 µl of lysate was loaded per lane on a NuPAGE 4–12% Bis-Tris gel (Invitrogen) and ran at 160 V for 1 h. Proteins were transferred at 100 V for 1 h using the wet-transfer method onto PVDF membranes. Membranes were blocked for 1 h with TBS intercept blocking buffer (LI-COR Biosciences), followed by incubation overnight at 4°C with mouse anti-GFP primary antibodies (1:10,000, Roche, 11814460001). Membranes were washed thrice with TBS containing 0.05% Tween 20, followed by incubation with IRDye 800CW IgG-conjugated goat anti-mouse secondary antibodies (1:10,000, LI-COR Biosciences, 926-32210) for 1 h. Proteins were detected using the Odyssey Infrared Imaging System (LI-COR Biosciences). GFP-tagged sensor levels were normalized using mouse anti-β-actin (1:10,000, Abcam, ab8224) and expressed relative to the wild type. Quantifications of band intensities were performed on the LI-COR Image Studio software.
Acknowledgements
We are grateful to the Oliferenko laboratory for discussions and Eugene Makeyev for suggestions on the manuscript. We thank Lydia Thompson and the Crick Flow Cytometry STP for help with FACS experiments, and Foo Juat Chin and the Singapore Lipidomics Incubator laboratory members for help with lipidomics experiments. Work in the S.O. laboratory was supported by the Francis Crick Institute, which receives its core funding from Cancer Research UK (CC0102), the UK Medical Research Council (CC0102) and the Wellcome Trust (CC0102).
Footnotes
Author contributions
Conceptualization: S.F., S.O.; Methodology: S.F., A.C.-G.; Validation: S.F.; Formal analysis: S.F., A.C.-G.; Investigation: S.F.; Resources: S.F.; Data curation: S.F.; Writing - original draft: S.F.; Writing - review & editing: S.F., A.C.-G., M.R.W., S.O.; Supervision: M.R.W., S.O.; Project administration: S.O.; Funding acquisition: M.R.W., S.O.
Funding
S.F. was in part supported by a King's College London–National University of Singapore joint PhD scholarship. Work in the M.R.W. laboratory is supported by grants from the National University of Singapore via the Life Sciences Institute, the National Research Foundation Singapore (NRFSBP-P4) and the National Research Foundation Singapore and Agency for Science, Technology and Research Industry Alignment Fund–Industry Collaboration Project (IAF-ICP) grant I1901E0040. Work in the S.O. laboratory was supported by the Francis Crick Institute, a Wellcome Trust Senior Investigator Award (103741/Z/14/Z), a Wellcome Trust Investigator Award in Science (220790/Z/20/Z) and a Biotechnology and Biological Sciences Research Council grant (BB/T000481/1) to S.O. Open access funding provided by The Francis Crick Institute. Deposited in PMC for immediate release.
Data availability
All relevant data can be found within the article and its supplementary information.
Peer review history
The peer review history is available online at https://journals.biologists.com/jcs/article-lookup/doi/10.1242/jcs.260568.reviewer-comments.pdf
References
Competing interests
The authors declare no competing or financial interests.