Asymmetric clustering of centrosomes defines the early evolution of tetraploid cells
Abstract
Tetraploidy has long been of interest to both cell and cancer biologists, partly because of its documented role in tumorigenesis. A common model proposes that the extra centrosomes that are typically acquired during tetraploidization are responsible for driving tumorigenesis. However, tetraploid cells evolved in culture have been shown to lack extra centrosomes. This observation raises questions about how tetraploid cells evolve and more specifically about the mechanisms(s) underlying centrosome loss. Here, using a combination of fixed cell analysis, live cell imaging, and mathematical modeling, we show that populations of newly formed tetraploid cells rapidly evolve in vitro to retain a near-tetraploid chromosome number while losing the extra centrosomes gained at the time of tetraploidization. This appears to happen through a process of natural selection in which tetraploid cells that inherit a single centrosome during a bipolar division with asymmetric centrosome clustering are favored for long-term survival.
Introduction
Organismal polyploidy is confined to certain taxa, but many species across the tree of life are thought to have had polyploid ancestors at some point in their evolutionary history (Comai, 2005; Smith and Keinath, 2015; Smith et al., 2013; Wolfe, 2001) and polyploidy is thought to contribute to speciation and evolution (Bennett, 2004; Otto and Whitton, 2000).
In vertebrates, organismal polyploidy is rare, and among mammals it has only been described in a single species (Gallardo et al., 1999). However, within individual diploid mammals, some tissues physiologically develop to have a higher ploidy than the majority of somatic cells (Nagl, 1990; Orr-Weaver, 2015; Øvrebø and Edgar, 2018; Schoenfelder and Fox, 2015). Polyploidy can also occur outside of the context of normal development, and is linked with both pathology (particularly cancer [Ganem et al., 2007]) and aging (Tanaka et al., 2015). Tetraploid cells are commonly found in premalignant lesions and tumors at different stages (Davoli and de Lange, 2011; Galipeau et al., 1996; Olaharski et al., 2006). Furthermore, meta-analysis of catalogued tumor genomes has provided evidence that close to 40% of all cancers – even those that were not tetraploid at the time of sampling – had a tetraploid intermediate stage at some point during tumor evolution (Zack et al., 2013). Consistent with this, several studies have shown a direct, causative link between tetraploidy and tumorigenesis (Fujiwara et al., 2005; Nguyen et al., 2009).
In proliferating cells, tetraploidy can arise via abnormal cell cycle events, including cytokinesis failure, cell fusion, endoreduplication, and mitotic slippage (Davoli and de Lange, 2011; Dikovskaya et al., 2007; Edgar and Orr-Weaver, 2001; Larsson et al., 2008; Rieder, 2011). Most of these events result in the concomitant acquisition of extra centrosomes along with genome duplication. Importantly, both tetraploidy and extra centrosomes have been shown to trigger a p53-dependent arrest in some experimental systems (Andreassen et al., 2001; Fava et al., 2017). In experimental systems in which such an arrest does not occur, extra centrosomes have been shown to promote chromosomal instability (Ganem et al., 2009; Silkworth et al., 2009) and invasive/migratory behavior (Godinho et al., 2014). Recent studies have also shown that extra centrosomes promote and in some cases are sufficient to drive tumorigenesis in vivo (Levine et al., 2017; Serçin et al., 2016).
Based on these studies, it has been speculated that the extra centrosomes emerging as a result of tetraploidization may drive chromosomal instability and, in turn, tumorigenesis (Storchova and Pellman, 2004). However, it was previously reported that cytokinesis failure does not result in stable centrosome amplification in a cell population (Krzywicka-Racka and Sluder, 2011). Moreover, anecdotal reports (Ganem et al., 2009; Godinho et al., 2014; Kuznetsova et al., 2015; Potapova et al., 2016) have indicated that clones of tetraploid or near-tetraploid cells displayed normal centrosome numbers. This suggests that our understanding of the evolution of tetraploid cells is incomplete and how centrosome and chromosome numbers evolve after tetraploidization needs to be revisited. To address this problem, we studied the time period immediately following cytokinesis failure and investigated how centrosome and chromosome numbers change in newly formed tetraploid cells. Following the observation that the number of centrosomes, but not chromosomes, rapidly returns to normal, we combined computational and experimental approaches to identify a specific cellular mechanism that underlies the loss of extra centrosomes.
Results
To investigate the early consequences of tetraploidy and the evolution of newly formed tetraploid cells, we induced cytokinesis failure by dihydrocytochalasin B (DCB) treatment for 20 hr (Figure 1A) in both DLD-1 (pseudodiploid colorectal cancer cells) and p53-/- hTERT-immortalized RPE-1 cells (Izquierdo et al., 2014) (hereafter referred to as RPE-1 p53-/-; p53-null RPE-1 cells were used because the parental, p53-positive, cells display a G1 arrest after cytokinesis failure, as shown in Ganem et al., 2014). Cells generated by this method are referred to, throughout the paper, as ‘newly formed tetraploid cells’ (text) or ‘4N new’ (figures).
Newly formed tetraploid cells undergo diverse fates in their first mitotic division
To determine the fate of the first tetraploid mitosis, we performed live-cell phase contrast microscopy for 24 hr following DCB washout (Figure 1A; note, newly formed tetraploid cells can easily be identified by the presence of two nuclei). We found that multipolar divisions were frequent in both cell types (Figure 1B,C), consistent with the acquisition of extra centrosomes upon cytokinesis failure and with the ability of extra centrosomes to promote formation of multipolar mitotic spindles. We confirmed the high rates of multipolar divisions by analyzing ana-/telophase cells immunostained for α–tubulin and centrin (Figure 1—figure supplement 1A,B). This fixed-cell analysis and analysis of live RPE-1 p53-/- cells with GFP-tagged centrin (Figure 1—figure supplement 1C,D) also confirmed that spindle poles always contained centrosomes (i.e., two centrin dots), as acentrosomal poles were never observed in multipolar ana-/telophase cells. Furthermore, all of the live binucleate cells that we observed contained supernumerary centrosomes (Figure 1—figure supplement 1D). These centrosomes were duplicated prior to mitosis and were never lost/extruded during mitosis (Figure 1—figure supplement 1D). The observation that only ~20–30% of newly formed tetraploid cells underwent bipolar anaphase indicates that centrosome clustering is not prevalent in newly formed tetraploid cells. However, when we followed these multipolar mitoses through cytokinesis, we found that tripolar or tetrapolar anaphases did not always generate three or four daughter cells, respectively (Figure 1D–F). Instead, the DNA corresponding to two or more anaphase poles was often enclosed in a single daughter cell, giving rise to binucleated or, rarely, trinucleated daughter cells (Figure 1D–F), consistent with previous observations (Chen et al., 2016; Wheatley and Wang, 1996).
Highly aneuploid cells form early in the evolution of tetraploid cells, but quickly disappear from the population
We next investigated how these early cell divisions after tetraploidization may impact chromosome numbers in the proliferating cell population. To this end, we carried out a time-course experiment in which we performed chromosome counting in the cell population after the 20 hr DCB treatment and every two days thereafter for a 12 day period (Figure 2A–B). Immediately after drug washout, we observed a tetraploid fraction corresponding to approximately 80% and 60% of the population for DLD-1 and RPE-1 p53-/- cells, respectively (Figure 2C–E,G, day 0). Two days after DCB washout, we observed high frequencies of cells with chromosome counts in the hypotetraploid/hyperdiploid range (Figure 2B, middle panel; Figure 2C–E,G). It is conceivable that these highly aneuploid cells may originate from the multipolar mitoses we observed in our live-cell imaging experiments. Indeed, by quantifying DNA fluorescence of separated chromosome masses in fixed ana-/telophase cells (Figure 2—figure supplement 1A), we observed that while chromosome distribution to the daughter cells was balanced in bipolar divisions (Figure 2—figure supplement 1B–D), it greatly deviated from an equal distribution in multipolar mitoses (Figure 2—figure supplement 1B–D), indicating that the daughter cells in multipolar divisions inherit variable proportions of the genome. However, the fraction of cells with highly abnormal chromosome numbers rapidly decreased over the course of the 12 day experiment in both cell types and highly aneuploid cells were virtually eliminated from the DLD-1 population by day 12, leaving sub-populations of near-diploid cells (presumably derived from cells that did not respond to the initial DCB treatment) and near-tetraploid cells (Figure 2B,C–E,G). The appearance and loss of highly aneuploid cells from the RPE-1 p53-/- cell population was delayed compared to the DLD-1 cell population (Figure 2G), possibly suggesting that newly formed tetraploid RPE-1 p53-/- cells display lower proliferation and death rates. Indeed, the appearance and disappearance of highly aneuploid cells corresponded to an increase followed by a decline in the fraction of dead cells between day 2 and 6 in both cell types (Figure 2F,H). However, cell death rates were lower in RPE-1 p53-/- compared to DLD-1 cells (Figure 2F vs. Figure 2H), possibly explaining the delay in both elimination of highly aneuploid cells and proliferation of ~4N cells in the RPE-1 p53-/- cell population (Figure 2G).
To investigate a possible cause for the disappearance of cells with highly aneuploid chromosome counts, we built a probabilistic model to evaluate the karyotypic outcomes of multipolar divisions (see Materials and Methods and Figure 3A). The model predicted that daughter nuclei emerging from multipolar divisions in tetraploid cells were very likely to bear a monosomy or nullisomy for at least one chromosome (Figure 3B,C). Because nullisomic cells and cells with certain monosomies are expected to be unable to proliferate further, the daughters of multipolar divisions would be expected to display lower proliferation rates than the daughters of bipolar divisions. This was confirmed in long-term time lapse microscopy experiments in newly formed tetraploid DLD-1 cells (Figure 3D), which showed that daughter cells produced by multipolar divisions were more likely to die or arrest compared to cells produced by bipolar divisions (Figure 3E, Figure 3—figure supplement 1). Previous studies have shown that diploid and near-triploid cells undergoing multipolar divisions produce daughter cells that arrest or die (Ganem et al., 2009; Gisselsson et al., 2010). Our data indicate that this is also the case for tetraploid cells undergoing multipolar divisions.
Supernumerary centrosomes that arise through tetraploidization quickly disappear from the population
The observation that highly aneuploid cells make up only a very small fraction of the population by 12 days after tetraploidization suggests that the rate of multipolar divisions (which generate these cells) decreases over time. This could be due to either an increased ability of the extra centrosomes to cluster in the tetraploid cells or to elimination of the extra centrosomes. To explore this, we investigated if and how centrosome number varies over the same 12 day evolution period by analyzing cells immunostained for centrin immediately following cytokinesis failure and every two days thereafter (Figure 4A).
We performed this analysis in both mitotic and G1 cells and obtained similar results at all time points, except immediately following DCB washout (‘Day 0’ in Figure 4C,D). This discrepancy at day 0 could be explained by a delay in mitotic entry of newly formed 4N cells, particularly for the RPE-1 p53-/- cells. Despite this difference at day 0, the trend was clear: the fraction of the G1 cell population containing supernumerary centrioles after a 20 hr cytokinesis block was 90% and 87.3% in DLD-1 and RPE-1 p53-/- cells, respectively. However, this fraction rapidly diminished over the 12 day observation period, reaching frequencies that are close to the frequencies of cells with supernumerary centrioles in the parental populations. Moreover, the fraction of cells with supernumerary centrioles at day 12 (Figure 4C–D) was substantially smaller than the fraction of cells with ~4N chromosome number (16.3% and 13.3% vs. 63% and 33%, respectively, in DLD-1 and RPE-1 p53-/-; compare ‘Day 12' data from Figure 4C,D and Figure 2E,G). Indeed, statistical analysis showed a highly significant difference between the number of cells with ~4N chromosome number and cells with extra centrosomes at day 12 for both DLD-1 and RPE-1 p53-/- cells (two-sided Chi square test, p<0.0001 for both cell lines), indicating that a large fraction of the tetraploid cells that are present 12 days post-cytokinesis failure have lost their extra centrosomes.
Tetraploid cells can inherit a normal centrosome number through asymmetric centrosome clustering during cell division
We reasoned that one way in which tetraploid cells could regain a normal centrosome number while maintaining tetraploid chromosome numbers would be by asymmetrically clustering the centrosomes during formation of a bipolar spindle. As a result, one daughter cell would receive three centrosomes and the other daughter would receive one centrosome, but both would receive a chromosome number ~4N.
To explore this possibility, we analyzed bipolar mitotic DLD-1 and RPE-1 p53-/- cells fixed and immunostained for centrin immediately following washout of DCB (Figure 5A–B). We found nearly equal numbers of bipolar DLD-1 cells with symmetric vs. asymmetric centrosome clustering in late mitosis (metaphase, anaphase, or telophase), while bipolar RPE-1 p53-/- showed a slight bias towards symmetric clustering of centrosomes (Figure 5C).
For asymmetric centrosome clustering to explain evolution of a tetraploid cell population with normal centrosome number, one would also have to assume that the daughter cell inheriting a single centrosome has a selective advantage over the daughter cell inheriting extra centrosomes (e.g., due to the likelihood of multipolar division in cells with extra centrosomes). To test this, we built a mathematical model based on this assumption (for model details, see Materials and Methods, Figure 5—figure supplements 1–4, and Table 1 ). We started with a simple model (Model I, Figure 5—figure supplement 1A–C) in which, initially, 87–90% of the cells have two centrosomes in G1 (four in S/G2/M), corresponding to the experimentally observed frequencies after DCB treatment. Cells in the model can divide in a multipolar or bipolar fashion, and bipolar divisions can occur with either symmetric or asymmetric centrosome clustering (Figure 5—figure supplement 1B–C) – all with probabilities that reflect those observed experimentally (see Table 1 for details on which experimental data motivated various model parameters). The daughter cells from multipolar divisions have significantly reduced viability and are expected to be quickly eliminated by selection; based on this, in the model these cells were, for simplicity, assumed to die (Figure 5—figure supplement 1A–C). Cells inheriting a single centrosome were assumed to become stable cells that undergo bipolar divisions with high viability (Figure 5—figure supplement 1B–C). Cells inheriting two centrosomes would display the same fate as newly formed tetraploid cells, and cells inheriting three centrosomes were assumed to undergo multipolar division and consequently produce non-viable progeny (Figure 5—figure supplement 1B–C). Although this model (Model I) captured centrosome loss, it predicted centrosome loss over a much shorter time scale than was observed experimentally for either DLD-1 or RPE-1 p53-/- cells (Figure 5D–E, orange line). The final fraction of cells with extra centrosomes predicted from the model was also substantially lower than what was experimentally observed (Figure 5D–E, orange line). Parameter optimization within a reasonable range could not solve this discrepancy (Figure 5—figure supplement 3A). In particular, the final steady-state fraction of cells with extra centrosomes was strongly dependent on the probability of cytokinesis failure in cells with normal centrosome number (which generates new cells with extra centrosomes) (Figure 5—figure supplements 1D and 3A). For the experimentally quantified (Nicholson et al., 2015) probability (~2.5%) of spontaneous cytokinesis failure in DLD-1 cells, the steady-state fraction of cells with extra centrosomes cannot match the observed value.
We next considered the possibility that a sub-fraction of newly formed tetraploid cells may cluster their centrosomes more efficiently than other cells (Model II, Figure 5—figure supplement 2), which we herein dub ‘super-clustering’ (SC) cells. When such SC cells were included in the model (Model II, Figure 5—figure supplement 2), the model output was no longer constrained by the probability of cytokinesis failure (Figure 5—figure supplement 3B) and the final fraction of cells with extra centrosomes could match the experimentally observed values (Figure 5D–E, blue line). Moreover, our model results showed that these SC cells would persist in the population and therefore dominate the final population of cells with supernumerary centrosomes (Figure 5F–G).
Altogether, our modeling results show that asymmetric centrosome clustering, along with a selective advantage of cells that inherit a single centrosome, is sufficient to explain the loss of extra centrosomes in newly formed tetraploid cells, leading to the evolution of cell populations with tetraploid chromosome numbers but normal centrosome numbers (i.e., 1 centrosome, 2 centrioles in G1 tetraploid cells).
Long-term live-cell imaging confirms that centrosome elimination and stable tetraploid cells arise via asymmetric centrosome clustering and natural selection
To directly observe the process of centrosome loss and test the model assumption that cells inheriting a single centrosome from a bipolar division are the most likely to keep proliferating, we performed live cell imaging experiments in DLD-1 and RPE-1 p53-/- cells expressing GFP-tagged centrin.
Because previous observations (Ganem et al., 2009; Gisselsson et al., 2010) and our own data (Figure 3C,E) indicated that the progeny of multipolar divisions display reduced viability, and since only bipolar or near-bipolar divisions are likely to generate the evolved (day 12) near-tetraploid cell population observed in our time-course experiment, we focused on fates of daughter cells arising from bipolar divisions. We imaged newly generated tetraploid (binucleate) cells by phase contrast microscopy for 24 hr, after which we determined the number of GFP-centrin dots present in the daughter cells arising from bipolar divisions. These cells were then imaged for an additional 48 hr by phase contrast microscopy (Figure 6A) to determine their fates in relation to the number of centrosomes they inherited. We found that cells that divided in a bipolar manner showed a mix of symmetric and asymmetric centrosome clustering without a strong preference for one mode (Figure 6—figure supplement 1), consistent with our fixed-cell data (Figure 5A–C; two-sided Fisher’s exact test, p=0.8224 and p=0.2243 for fixed vs. live cell data in DLD-1 and RPE-1 p53-/-, respectively). Cells that inherited a normal centrosome number (1 centrosome/2 centrioles) were significantly more likely than cells that inherited supernumerary centrosomes to divide in a bipolar manner in both DLD-1 and RPE-1 p53-/- cells (Figure 6B–C). In contrast, cells that inherited too many centrosomes went through a mix of fates, dominated by multipolar divisions, arrest, and cell death (Figure 6—figure supplement 1). These data, together with our mathematical modeling, strongly suggest that populations of stably dividing tetraploid cells containing a normal number of centrosomes can arise via asymmetric clustering of centrosomes (3:1) in bipolar mitoses and selective pressure against cells that inherit extra centrosomes.
In generating our Model II, we included the assumption that a fraction of cells with extra centrosomes had a very high efficiency of centrosome clustering (‘SC cells’). This assumption was required to reproduce the observed evolution dynamics and final fraction of cells with extra centrioles, given the observed rate of cytokinesis failure. Based on this assumption, the model predicted that the fraction of SC cells rapidly increased over the first few days and that SC cells would make up about 90% of the remaining cell population with extra centrosomes at the end of the 12 day evolution period (Figure 5F,G). To test this model prediction, we analyzed fixed DLD-1 cells in ana-/telophase to determine the fractions of cells with extra centrioles that displayed bipolar vs. multipolar configurations. As predicted by our model, we found that approximately 90% of cells with extra centrioles displayed a bipolar configuration at day 12 compared with just 28% in newly formed tetraploid cells at Day 0 (Figure 6D). These results suggest that when extra centrosomes arise, they may only be retained when cells can cluster them efficiently, whereas cells that cannot cluster their centrosomes efficiently may disappear from the population.
Discussion
Newly formed tetraploid cells rapidly lose the extra centrosomes while maintaining the extra chromosomes
Here, we show, in two different cell lines, that populations of newly formed tetraploid cells rapidly evolve in vitro to retain a near-tetraploid chromosome number while losing the extra centrosomes gained at the time of tetraploidization. By combining fixed cell analysis, live cell imaging, and mathematical modeling, we show that this happens by a process of natural selection (Figure 7). Specifically, tetraploid cells that inherit a single centrosome during a bipolar division with asymmetric centrosome clustering are favored for long-term survival. Conversely, the majority of cells with extra centrosomes are eventually eliminated because of their high probability of undergoing multipolar division, which has a high likelihood of producing daughters with insufficient amounts of genetic material (Figure 7).
Our findings can explain previous anecdotal reports (Ganem et al., 2009; Godinho et al., 2014; Kuznetsova et al., 2015; Potapova et al., 2016) that clones isolated after experimental inhibition of cytokinesis consisted of tetraploid cells with a ‘normal’ number (i.e., same number as in diploid cells) of centrosomes. Our study also shows that this pattern of centrosome number evolution after tetraploidization is common to both cancer (DLD-1) and non-cancer (RPE-1) cells. Above all, our work reveals the mechanism (Figure 7) by which tetraploid cells containing a normal number of centrosomes emerge. Finally, our mathematical model successfully captures tetraploid cell evolution and may be used in the future to better understand how tetraploidy contributes to tumor initiation and progression in situ.
Tetraploidization and tumorigenesis: the case for extra chromosomes, extra centrosomes, or both as driving factors
The link between tetraploidization and tumorigenesis is supported by strong experimental evidence. Cancer genome sequencing data indicated that tetraploidization occurs at some point during the progression of a large fraction of tumors (Zack et al., 2013). Moreover, tetraploid mouse epithelial cells were shown to be more tumorigenic than their non-tetraploid counterparts when injected in nude mice (Fujiwara et al., 2005; Nguyen et al., 2009; Davoli and de Lange, 2012). A popular model for how tetraploidy may promote tumorigenesis is that the extra centrosomes (which arise concomitantly with tetraploidization) contribute to cancer phenotypes (Storchova and Pellman, 2004). This idea is supported by the following observations: centrosome amplification is frequently observed in the pre-malignant stages of certain cancers (Chan, 2011; Lopes et al., 2018) and is observed in a large fraction of human tumors (D'Assoro et al., 2002; Gustafson et al., 2000; Lingle et al., 1998; Pihan et al., 1998; Sato et al., 1999), in which it correlates with poor prognosis/advanced disease stage (Lopes et al., 2018; Godinho and Pellman, 2014); extra centrosomes can promote tumorigenesis in mouse (Levine et al., 2017; Serçin et al., 2016) and enhance the invasive behavior of mammary epithelial cells grown in 3D cultures (Godinho et al., 2014); finally, supernumerary centrosomes promote chromosome mis-attachment and mis-segregation (Ganem et al., 2009; Silkworth et al., 2009), leading to chromosomal instability, a hallmark of cancer believed to drive tumor evolution (Targa and Rancati, 2018). Together, these observations indicate that extra centrosomes are likely to contribute to tumor initiation and/or progression.
Our observation that extra centrosomes gained through tetraploidization are quickly lost raises the possibility that tetraploidy may drive tumorigenesis by means other than the acquisition of extra centrosomes. Indeed, a number of studies and observations suggest that tetraploidy per se may promote the emergence of cancer phenotypes. For instance, tetraploidy was shown to increase tolerance for genomic changes, leading to the rapid evolution of complex genomes (Dewhurst et al., 2014), as seen in cancer, and many tetraploid cells show increased chromosomal instability compared to diploid cells, even when no extra centrosomes are present (Kuznetsova et al., 2015). Additionally, polyploid cells were shown to be more resistant than their diploid counterparts to oxidative stress, genotoxic insult, irradiation, and certain chemotherapeutic drugs (Kuznetsova et al., 2015; Ianzini et al., 2009; Illidge et al., 2000). Lastly, in cancer patients, genome-doubling in early stage tumors was shown to correlate with poor relapse-free survival (Dewhurst et al., 2014), although centrosome number was not examined in these patients.
In light of our findings, one could imagine that in certain instances, cells experiencing a genome doubling event may initially lose their extra centrosomes and then re-acquire them at a later time, depending on additional factors. At least one example in the literature provides evidence for such a series of events. In Barrett’s esophagus, a pre-malignant condition that predisposes to esophageal cancer (Cameron et al., 1985; Hameeteman et al., 1989; Hvid-Jensen et al., 2011), accumulation of 4N cells has been shown to occur as the tissue transitions to metaplasia (Galipeau et al., 1996). A study on centrosome status in Barrett’s esophagus reported centrosome amplification prior to the transition to metaplasia (Segat et al., 2010), corresponding to the time when tetraploid cells accumulate (Galipeau et al., 1996), but also noted that the frequency of supernumerary centrosomes decreased with progression to metaplasia and neoplasia (Segat et al., 2010). Similarly, another study found an increase in centrosome amplification followed by a decrease during the progression from Barrett’s Esophagus to adenocarcinoma (Lopes et al., 2018). These results closely mirror the dynamics of evolution seen in our study and illustrate that extra centrosomes can be present early in tumor development (around the time when tetraploidy appears) but subsequently be lost. Therefore, while tetraploidy and supernumerary centrosomes are both individually linked with tumorigenesis, the link between tetraploidy, extra centrosomes, and disease progression may be less direct than conventionally thought.
Tetraploidization is intimately linked with the birth of extra centrosomes; however, tetraploidization may not lead to stable acquisition of supernumerary centrosomes unless (i) specific cellular/genetic changes have occurred to allow the cell to maintain its extra centrosomes and/or (ii) certain conditions in the tissue microenvironment exist that favor or necessitate the presence of extra centrosomes. Indeed, there is evidence that clustering of extra centrosomes into a bipolar configuration can be influenced by a number of cell intrinsic and extrinsic factors, and such factors may be important for determining the fraction of cells that retain extra centrosomes. Among cellular factors, the nonessential motor protein KIFC1 (also known as HSET), the epithelial cell protein E-cadherin, spindle assembly checkpoint components, the chromosome passenger complex, the NDC80 complex, and the augmin complex have been shown to affect (either positively or negatively) centrosome clustering efficiency (Kwon et al., 2008; Leber et al., 2010; Quintyne et al., 2005; Rhys et al., 2018; Sabino et al., 2015). Extracellular factors, such as geometric constraints imposed by the environment, have also been shown to alter centrosome clustering (Kwon et al., 2008). Alternatively, if tetraploidization occurs under circumstances that do not favor retention of extra centrosomes, the tetraploid cells may initially lose their extra centrosomes and then re-acquire them at a later time, as a result of genome instability, which may lead to the non-stoichiometric production of proteins involved in centrosome duplication. Thus, the evolutionary pattern that newly formed tetraploid cells will follow may vary depending on many factors, including genetic background, functional requirements in a given tissue/organ, or a variety of extracellular physical and physiological factors. All these potential factors could explain the high rates of extra centrosomes in certain tumors and animal models (Levine et al., 2017; Serçin et al., 2016).
Materials and methods
Experimental approaches
Cell lines and culture conditions
Request a detailed protocolDLD-1 cells (ATCC CCL-221) were purchased from the American Type Culture Collection (ATCC, Manassas, VA). The hTERT immortalized RPE-1 p53-/- cell line (Izquierdo et al., 2014) (referenced throughout the manuscript as RPE-1 p53-/-) was a gift from Dr. Meng-Fu Bryan Tsou (Memorial Sloan Kettering Cancer Center). Both the DLD-1 and the original hTERT RPE-1 cell lines originated from ATCC. The company provides certification. Potential mycoplasma infection was monitored regularly (no less than once every three weeks) by DNA staining (DAPI) and any cell batch with suspected mycoplasma infection was discarded. DLD-1 cells were cultured according to ATCC recommendations in RPMI 1640 medium with ATCC modification (Thermo Fisher Scientific – Gibco, CA, USA) supplemented with 10% fetal bovine serum (FBS; Thermo Fisher Scientific – Gibco, CA, USA) and 1% antibiotic-antimycotic (Thermo Fisher Scientific – Gibco, CA, USA). RPE-1 p53-/- cells were cultured according to the ATCC recommendations for hTERT-immortalized RPE-1 cells in 1:1 mixture of DMEM/F-12 with HEPES (Thermo Fisher Scientific – Gibco, CA, USA) also supplemented with 10% FBS and 1% antibiotic-antimycotic. All cells were grown on tissue culture polystyrene flasks (Corning, Tewksbury, MA) and were maintained in a humidified incubator at 37°C and 5% CO2.
Tetraploid DLD-1 and RPE-1 p53-/- cells were generated by treating diploid cell cultures with 1.5 µg/mL dihydrocytochalasin B (DCB; Sigma Aldrich, Saint Louis, MO) for 20 hr. For live cell experiments, cells were washed out (4 times with 1X PBS) into imaging medium and immediately taken to the microscope for imaging following.
Generating virally transduced cell lines
Request a detailed protocolThe GFP-Centrin 2 gene was PCR amplified from a modified pLL3.7 plasmid with unknown selection (gift of Tim Stearns, Stanford University), using forward and reverse primers designed to match the two termini of the fusion protein. The forward and reverse primers used (including restriction sites for NOTI and SALI and terminal non-sense nucleotides) were (with start and stop codons underlined):
(F)CAATAAAGCGGCCGCATGGTGAGCAAGGGCGAGGAGCTGT and
(R)GGACTGGTGGTCTGCGTCGACTTAATAGAGGCTGGTCTTTTTCATG.
Cleaned PCR product was ligated into the pLNXC2 retroviral expression vector by directional cloning using NOTI and SALI restriction enzymes. The presence of GFP-centrin 2 gene in plasmid DNA was confirmed by restriction digests visualized on DNA gels and via transient transfection into RPE-1 p53-/- cells to confirm centrosomal localization. GFP-Centrin expressing DLD-1 and RPE-1 p53-/- cells were generated by transduction with retroviral particles. GP-293 cells containing retroviral gag and pol genes (ClonTech Laboratories Inc, Mountain View, CA #631458) were co-transfected with the expression vector and the pVSV-G plasmid (Addgene, Cambridge, MA). 48 hr after transfection, supernatant was collected, filtered through a 0.45 µm pore (GD/X sterile 0.45 µm CA filter, GE Whatman PLC, Pittsburgh, PA), mixed with polybrene (Sigma-Aldrich, Saint Louis, MO) at a final concentration of 10 µg/ml, and added directly to the cells. After 24 hr, cell medium was replaced with fresh culture media. Starting 72 hr after viral transduction, transduced cells were selected with with G418 at a concentration of 500 µg/ml until negative control cells (untransduced cells treated with the same concentration of antibiotic) were dead, or approximately two weeks.
Cells co-expressing RFP-H2B were generated by further transducing GFP-Centrin 2 expressing cells, via the protocol described previously, using a pBABE retroviral plasmid containing RFP-H2B and a puromycin selection gene (gift from Neil Ganem, Boston University). Transduced cells were selected with puromycin at a concentration of 5 µg/ml (RPE-1 p53-/-) or 3.8 µg/ml (DLD-1).
Phase contrast live cell microscopy
Request a detailed protocolFor live-cell experiments, all cells were grown on MatTek glass bottom dishes with No. 1.5 glass (MatTek Corporation, Ashland, MA). At the time of imaging, cell medium was replaced with L-15 medium supplemented with 4.5 g/l glucose (high glucose). All live cell experiments were performed on a Nikon Eclipse Ti inverted microscope (Nikon instruments Inc, NY, USA) equipped with phase-contrast trans-illumination, transmitted light shutter, ProScan automated stage (Prior Scientific, Cambridge, UK), CoolSNAP HQ2 CCD camera (Photometrics, AZ, USA), Lumen200PRO light source (Prior Scientific, Cambridge, UK), and a temperature and humidity controlled incubator (Tokai Hit, Japan). For 24 hr and 72 hr live cell phase contrast videos, images were acquired every 6 min through a 20X/0.3 NA A Plan corrected phase contrast objective for the duration of the experiment. Time-lapse videos were analyzed using NIS Elements AR software (Nikon Instruments Inc, NY, USA) to determine the nature of division (bipolar, tripolar, tetrapolar) at anaphase and the subsequent number of daughter cells formed after cytokinesis.
Time course experimental procedure
Request a detailed protocolTime course (12 day) experiments were performed by seeding all cells needed for the first two time points (day 0 and day 2) along with a flask designated for propagating the experiment. For each replicate for DLD-1 cells, this included T-25 flasks seeded with 1 × 106 (day 0 metaphase spreads) and 5 × 105 (day 2 metaphase spreads), a T-75 flask with 1 × 106 cells, and acid-washed coverslips inside 35 mm Petri dishes with 2.5 × 105 (day 0) and 1 × 105 (day 2) cells for combined centrin/geminin immunostaining. On day 2, the T-75 flask was used to seed cells for the next two time points as follows: two T-25 flasks (metaphase spreads), one T-75 flask (propagating), and coverslips (centrin/geminin immunostaining). This was repeated for the entire 12 day period. The experiment was designed in the same way for RPE-1 p53-/- cells, but cell counts were as follows: T-25 flasks seeded at 1 × 106 cells (earlier time point, e.g. day 0) and 5 × 105 (later time point, e.g. day 2); T-75 seeded at 1.5 × 106 cells; coverslips seeded at 1.25 × 105 (earlier time point) and 8.5 × 104 cells (later time point).
Chromosome spread preparation and analysis
Request a detailed protocolCell cultures were grown in T-25 flasks as described in the previous section. For chromosome spread preparation, cells were incubated in their respective medium containing 50 ng/ml colcemid (Invitrogen – Karyomax, Waltham, MA) at 37°C for 5 hr to enrich for mitotically arrested cells. The cells were then collected by trypsinization and centrifuged at 1000 rpm for 5 min. Pre-warmed (37°C) hypotonic solution (0.075 M KCl) was added drop-wise to the cell pellet and incubated for 18 (DLD-1 cells) or 15 (RPE-1 p53-/- cells) minutes at 37°C. Several drops of freshly prepared fixative (3:1 methanol:glacial acetic acid) were added before centrifugation at 1000 rpm for 5 min. Supernatant was aspirated, fixative was added dropwise, and the cell suspension was again centrifuged at 1000 rpm for 5 min. The fixation step was repeated two more times and fixed cells were finally dropped on microscope slides. For RPE-1 p53-/- cells, a homemade humidity chamber constructed from PVC piping, plastic sheeting, and a nebulizer was used when drying slides to improve chromosome spread quality (effect of humidity on chromosome spread quality was described previously Deng et al., 2003). Chromosome spreads were then stained with 300 nM DAPI (Thermo Fisher Scientific – Invitrogen, Waltham, MA) for 10 min. DAPI-stained slides were mounted with an antifade solution containing 90% glycerol and 0.5% N-propyl gallate and sealed under a 22 × 50 mm coverslip (Corning Incorporated, Corning, NY) with nail polish. For chromosome counting, images of DAPI-stained chromosome spreads were acquired with the Nikon Eclipse Ti inverted microscope setup described in an earlier section. Images were acquired using either a 60X/1.4 NA or a 100X/1.4 NA Plan-Apochromatic phase contrast objective. After image acquisition, chromosome spreads were visualized in gray scale and chromosomes were counted using the count function in NIS elements.
Cell death assays
Request a detailed protocolTo measure cell death, 5 × 104 cells were plated in each of three wells of a 6-well plate and 1 × 106 cells were plated in a T-25 flask. The following day, cells were treated with 1.5 µg/ml DCB for 20 hr. After 20 hr, day 0 cells’ supernatant was collected, while the adherent cells were washed (3 times using PBS) and harvested in trypsin. The supernatant, all the washes, and the cell suspension were collected in the same tube, centrifuged, and re-suspended in 400 µl PBS for counting. Viable cells were differentiated from dead cells by trypan blue exclusion. The numbers of living and dead cells were counted and the fraction of dead cells out of the total number of cells was calculated. Cell counting was performed on days 0, 1, 2, and every 2 days for the remainder of the 12 day period (with new wells being seeded from T-25 flasks on day 2). Cell culture medium was changed 24 hr before counting each day in order to only count cells that died within a 24 hr period.
Immunofluorescence staining, image acquisition and data analysis
Request a detailed protocolFor centrin and geminin immunostaining, cells were grown on sterilized acid-washed glass coverslips inside 35 mm Petri dishes. Cells were fixed in freshly prepared 4% paraformaldehyde in PHEM buffer (60 mM Pipes, 25 mM HEPES, 10 mM EGTA, 2 mM MgSO4, pH 7.0) for 20 min at room temperature and then permeabilized for 10 min at room temperature in PHEM buffer containing 0.1% Triton-X 100. Following fixation and permeabilization, cells were washed three times with PBS and then blocked with 20% boiled goat serum (BGS) for 1 hr at room temperature. Cells were then incubated at 4°C overnight with primary antibodies diluted in 10% BGS. Next, cells were washed in PBS-T (PBS with 0.05% Tween 20) three times, and incubated at room temperature for 45 min with secondary antibodies diluted in 10% BGS. Cells were then washed four times with PBS-T, stained with DAPI (300 nM, Thermo Fisher Scientific – Invitrogen, Waltham, MA) for 5 min, and coverslips were mounted on microscope slides in an antifade solution containing 90% glycerol and 0.5% N-propyl gallate. For centrin/α-tubulin immunostaining, cells were washed in 1X PBS three times and fixed/permeabilized in 100% methanol for 10 min. After permeabilization, fixed cells were treated as described above for centrin/geminin staining. Primary antibodies were diluted as follows: rabbit anti-geminin (Abcam, Cambridge, MA), 1:100; mouse anti-centrin (Abnova, Zhongli, Taiwan), 1:100; rabbit anti-α-tubulin (Abcam, Cambridge, MA), 1:250. Secondary antibodies were diluted as follows: Rhodamine Red-X goat anti-rabbit (Jackson ImmunoResearch Laboratories, Inc, PA, USA), 1:100; Alexa 488 goat anti-mouse (Molecular Probes, Life Technologies, CA, USA), 1:200.
Centrin-stained samples were analyzed on a Nikon Eclipse TE2000 inverted microscope equipped with a 100X/1.4 NA Plan-Apochromatic phase contrast objective lens, motorized ProScan stage (Prior Scientific, Cambridge, UK), appropriate filter sets, and an XCITE 120Q light source (Excelitas Technologies, Waltham, MA, USA). Analysis was performed visually in both interphase cells and mitotic cells. The number of centrin dots was counted in cells that were determined to be in mitosis by DAPI staining. Mitotic cells with four centrin dots (i.e., two dots corresponding to each centrosome of a bipolar spindle) were categorized as normal; mitotic cells with greater than four dots were categorized as possessing supernumerary centrosomes. Interphase analysis was performed in G1/G0 cells, as determined by absence of nuclear geminin staining (McGarry and Kirschner, 1998). G0/G1 cells with two adjacent centrin dots (corresponding to a single centrosome) were classified as normal, whereas cells with greater than two centrin dots were classified as possessing supernumerary centrosomes. For centrosome clustering analysis (Figure 5A–C), bipolar metaphase, anaphase, or telophase cells were analyzed for the number of centrin dots present at respective spindle poles. For analysis of the fraction of fixed cells undergoing bipolar vs. multipolar division (Figure 1—figure supplement 1, Figure 6D), ana-/telophase cells stained with centrin and α-tubulin were analyzed for polarity (α-tubulin staining) and the presence or absence of centrioles at each spindle pole. To ensure that DCB treatment did not alter the polarity of mitotic cells, we co-treated cells with DCB and the CDK1 inhibitor RO-3306 (which causes a robust G2 arrest) for 20 hr, then washed both drugs out, waited for 1 hr for cells to proceed into mitosis, then fixed, stained and analyzed the relative proportion of bipolar and multipolar ana-/telophases in those cells naturally harboring supernumerary centrosomes. Representative z-stack image examples were acquired on the Nikon Eclipse Ti inverted microscope setup described in an earlier section. Images were acquired using either a 60X/1.4 NA or a 100X/1.4 NA Plan-Apochromatic phase contrast objective and appropriate filters.
For analysis of genome distribution in bipolar and multipolar divisions, images of ana-/telophase cells were acquired with a swept field confocal system (Prairie Technologies, WI, USA) on the same Nikon Eclipse TE2000-U inverted microscope described previously (Nikon Instruments Inc, NY, USA). The microscope was equipped with a CoolSNAP HQ2 CCD camera (Photometrics, AZ, USA), a multiband pass filter set (illumination at 405, 488, 561, and 640 nm), and an Agilent monolithic laser combiner (MLC400) controlled by a four channel acousto-optic tunable filter. Images were obtained by acquiring Z-stacks with 0.6 µm steps (Nyquist matched) so that the entire 3-D volume of the DNA was captured. Images were shading corrected using the NIS Elements shading correction function. Z-stacks were summed using the FIJI (Schindelin et al., 2012) sum slices function. The freehand selection tool was used to trace the signal area corresponding to an ana-/telophase chromosome cluster and the percentage of the signal intensity relative to total DNA for an ana-/telophase cell was determined. To calculate the symmetry score, the ratio between the actual fluorescence intensity percentage and the expected signal intensity percentage for an even distribution to 2 (50%), 3 (33.3%) or 4 (25%) poles (depending on the polarity of the division) was first calculated for each chromosome cluster. Then, the standard deviation of all measurements for a cell was calculated as a ‘symmetry score’ (ss). If a division was perfectly symmetrical, ss = 0 and any ss >0 indicates proportional increases in the asymmetry of DNA distribution to the poles.
Live cell imaging of fluorescently labeled cells
Request a detailed protocolFor live cell imaging of GFP-Centrin expressing cells, imaging was performed with a 60X/1.4 NA Plan-Apochromatic phase contrast objective lens (for RPE-1 p53-/- cells) or a 100X/1.4 NA Plan-Apochromatic phase contrast objective lens (for DLD-1 cells) controlled by Nikon Perfect Focus (Nikon Instruments Inc, NY, USA). In preparation for short-term live imaging of binucleate cells expressing GFP-centrin and RFP-H2B, the cells were washed out of DCB into medium containing 9 µM of the CDK1 inhibitor RO3306 to arrest cells at the G2/M transition. After 4 hr, the cells were again washed out into high glucose L-15 medium lacking phenol red. Imaging was performed by identifying individual binucleate cells in prophase or early prometaphase using RFP-H2B signal. Cells were imaged at the home Z-position in phase contrast every 4 min and the FITC channel every 4 min with asymmetrical Z-stacks defined by the home position and a range of −2.4 µm and +5.8 µm in 0.6 µm steps (14 steps). Chromosomes were imaged by phase contrast instead of fluorescence (RFP) to keep illumination, and hence photodamage, to a minimum, given that phase contrast imaging required lower exposure times than fluorescence imaging. Cells were imaged for a total of 3–4 hr. The videos were then analyzed to determine the number of centrin dots (centrioles) in the early mitotic cells and again in the resulting daughter cells after division.
For long-term cell fate experiments (Figure 6), GFP-Centrin expressing cells were used. Binucleate cells were imaged at 10 min intervals for 24 hr via phase contrast microscopy under a 60X/1.4 NA or 100X/1.4 NA Plan-Apochromatic phase contrast objective lens. Following this period, a number of daughter cells were selected and the number of centrioles was quickly counted for each by eye. A phase contrast image was obtained, along with asymmetric Z-stack images in the FITC channel, defined by the home position and a range of −2.4 µm and +5.8 µm in 0.6 µm steps. These daughter cells were then tracked via phase contrast microscopy at 10 min intervals for an additional 48 hr period to determine their subsequent fate.
Modeling approaches
Probabilistic model for karyotypic outcomes of multipolar divisions
Request a detailed protocolWe built the following model to evaluate the probabilities of nullisomy and/or monosomy in a cell division with p poles in a k-ploid mother cell, that is, a cell with k sets of M nonhomologous chromosomes (e.g., k = 2, M = 23 for normal, diploid human cells). For simplicity, we made the following assumptions:
The possibility of chromosome missegregation is ignored. Sister chromatids from each chromosome are partitioned to different spindle poles and end up in different daughter cells. The chromosome partitioning is otherwise random.
All chromosomes are partitioned in the same way as above and independent of one another.
Due to the second assumption, the probability of an event (e.g., nullisomy, monosomy, or nullisomy/monosomy) for at least one chromosome in a daughter cell reads as Equation 1. Because all chromosomes are equivalent in partitioning, the probability can be expressed in terms of the probability for Chr 1 without loss of generality.
Next, we need to determine the probability of each event of interest for Chr 1, and use Equation 1 to deduce the corresponding probability in the cell.
Probability of nullisomy
Request a detailed protocolBecause sister chromatids have to be partitioned to different poles, the total number of equal ways to partition one pair of sister chromatids to p poles reads as:
where the bracketed expression represents the binomial coefficient.
Because sister chromatids from each chromosome are independent of each other in the partitioning, the total number of equal ways to partition k pairs of sister chromatids to p poles reads as:
If any given pole receives 0 chromatids (i.e., nullisomy), then the total number of equal ways to partition k pairs of sister chromatids to the remaining p-1 poles reads as:
Hence, the probability that any given pole and the corresponding daughter cell bears a nullisomy for Chr 1 reads as:
Note that the probability in Equation 5 is not multiplied by another factor p for the number of possible poles/daughter cells, because we are looking for the probability of nullisomy of Chr 1 in a given daughter cell rather than in a given cell division.
Plugging Equation 5 into Equation 1 yields the probability of nullisomy in a cell.
Plugging M = 23, k = 4, p = 3 or 4 into Equation 6 yields the results presented in Figure 3C (white bars). Because all chromosomes are independent of each other, the number of nullisomies in a cell follows a binomial distribution , where according to Equation 5. The corresponding probability distribution for M = 23, k = 4, p = 3 or 4 is plotted in Figure 3B (top).
Probability of monosomy
Request a detailed protocolIf any given pole receives 1 chromatid (i.e., monosomy), then the total number of equal ways to partition the chromosomes reads as:
Hence, the probability that any given pole and the corresponding daughter cell bears a monosomy for Chr 1 reads as:
Plugging Equation 8 into Equation 1 yields the probability of monosomy in a cell.
Plugging M = 23, k = 4, p = 3 or 4 into Equation 9 yields the results presented in Figure 3C (grey bars). The number of monosomies in a cell follows a binomial distribution , where according to Equation 8. The corresponding probability distribution for M = 23, k = 4, p = 3 or 4 is plotted in Figure 3B (bottom).
Probability of nullisomy or monosomy
Request a detailed protocolBecause nullisomy and monosomy are mutually exclusive events for a given chromosome, e.g., Chr 1, the probability that any given pole and the corresponding daughter cell bears either nullisomy or monosomy for Chr 1 reads as:
Plugging Equation 10 into Equation 1 yields the probability of nullisomy or monosomy in a cell.
Plugging M = 23, k = 4, p=3 or 4 into Equation 11 yields the results presented in Figure 3C (black bars).
Model for centrosome number evolution in a cell population
Model I
Request a detailed protocolModel I was constructed based on the following minimal assumptions about how centrosome numbers vary during cell divisions (Figure 5—figure supplement 1A–C). The subscripts refer to the number of centrosomes in a cell during mitosis.
A cell with normal centrosome number (C2) undergoes normal division with probability q and cytokinesis failure (→ C4) with probability 1–q;
A cell with double centrosome number (C4) undergoes bipolar division with probability p and multipolar division with probability 1–p;
A bipolar division occurs in a symmetric fashion (2 C4) with probability r and in an asymmetric fashion (C2+C6) with probability 1–r.
A multipolar division of a C4 cell goes by 2 C2+C4 with probability s and 4 C2 with probability 1–s;
A multipolar division of a C4 cell in the fashion of 4 C2 is fatal;
A multipolar division of a C4 cell in the fashion of 2 C2+C4 only has C4 viable (equivalent to a normal C4) with probability f.
In addition,
C2 cells divide with rate bC2, and die with rate dC2;
C4 cells divide with rate bC4, and die with rate dC4;
C6 cells divide in multipolar fashion and die (there might be a small probability of viable division, which is neglected).
Based on the cell fate depicted in Figure 5—figure supplement 1B–C, the population dynamics are governed by the following ODEs:
with initial condition . The initial condition reflects the experimental observation that the newly induced 4N cell populations usually contain a certain fraction () of C2 (2N) cells.
Parameter sensitivity analysis (Figure 5—figure supplement 1D) indicated that the final fraction of cells with extra centrosomes strongly depends on q, the probability of cytokinesis failure in cells with normal centrosome number. In fact, the range of possible values for this final fraction is strongly constrained by the value of q, even if choice of all parameters could span a wide range (Figure 5—figure supplement 3A). This is because cytokinesis failure causes formation of new cells with extra centrosomes, and hence a large probability of cytokinesis failure is needed to maintain a higher fraction of these cells in the population.
Model II (with SC cells)
Request a detailed protocolIn the updated model (Figure 5—figure supplement 2A-C), we added SC cells, which are C4 cells that can effectively cluster extra centrosomes, and regularly undergo bipolar division. For this new cell type, we made the following assumptions.
Cytokinesis failure in cells with normal centrosome number results in SC cells with probability, v.
SC cells divide symmetrically (SC+SC) with a probability, rS. Otherwise, they divide asymmetrically (C2+C6).
SC cells have the same division and death rates as cells with normal centrosome number, because they are supposedly stable.
Based on the cell fate depicted in Figure 5—figure supplementv 2B-C, the population dynamics are governed by the following ODEs:
Parameter sensitivity analysis (Figure 5—figure supplement 2D) indicated that, based on Model II, the final fraction of cells with extra centrosomes is most sensitive to rS, the probability of symmetric division in SC cells, followed by q, the probability of cytokinesis failure in C2 cells, and v, the probability of getting SC cells upon cytokinesis failure. While Model I showed a strong constraint on q (Figure 5—figure supplement 3A), the strength of this constraint is relaxed in Model II (Figure 5—figure supplement 3B). In Model II, the major constraint is shifted to rS (Figure 5—figure supplement 3C), because asymmetric division (with probability 1- rS) controls the conversion of SC cells back to C2 cells. Nevertheless, ~90% probability of symmetric division is sufficient to maintain 20% cells with extra centrosomes in the steady state population.
Steady state of cell fractions
Request a detailed protocolWhen the cell division rate is sufficiently large compared to cell death rate in the models, the number of cells in each type will increase infinitely (Figure 5—figure supplement 4, left column). This case does reflect the experiments, in which the cell cultures were regularly sampled and re-populated on fresh medium, effectively creating a finite sample of the unlimited population growth. Although the total population grows infinitely, the fractions of each cell type approach fixed steady state values (Figure 5—figure supplement 4, right column). In fact, the steady state fraction of each cell type can be analytically solved as shown below.
Systems of homogenous linear ODE equations like Equations 12-14 and Equations 15-18 can be written in a vector form as
where is the list of variables.
The coefficient matrix, P, has the rate constants as entries. For Model I governed by Equations 12-14,
Likewise, for Model II governed by Equations 15-18,
If , then Equation 19 only has the trivial steady state where all variables equal zero. This trivial steady state is unstable if the overall proliferation rate is larger than the overall death rate. In other words, the whole cell population is expected to increase infinitely. Although the total population grows infinitely, the fraction of each cell type in the population could reach a steady state. To address this question via modeling, one can rewrite Equation 19 in terms of the fraction of each cell, that is,
Hence,
Equation 23 can be rewritten in vector format as
where and .
Because is a scalar function of time, at the steady state of Equation 24, approaches a constant, that is, . In other words, the steady state of Equation 24 is found when
Hence, the steady state solution of Equation 24 is an eigenvector of the coefficient matrix, P, normalized by the constraint, . equals the corresponding eigenvalue of P. We show in the following that is in fact the largest eigenvalue of P.
Theorem 1: The steady state solution of Equation 24 is given by the normalized eigenvector associated with the largest eigenvalue of the coefficient matrix, P, with the normalization condition, .
Heuristic proof:
At , the solution to Equation 24 approaches the solution to Equation 26.
The solution of Equation 26 reads
where ’s are eigenvalues of the matrix , and is the identity matrix.
At , Equation 27 is dominated by the term with the largest eigenvalue, that is,
A nonzero steady state solution requires . Note that the eigenvalues of P have one-to-one correspondence with the eigenvalues of Q. For each eigenvalue of Q, , is an eigenvalue of P. Because the largest eigenvalue of Q is 0, the largest eigenvalue of P is . The normalization constraint follows from the definition of fractions in Equation 22.
Based on Theorem 1, the steady state fractions of each cell type in the model can be obtained by computing the normalized eigenvector associated with the largest eigenvalue of the coefficient matrix, P, which can be easily done using a computation software, for example, MATLAB.
Data availability
All data generated during the study are provided in clearly labeled source data files in excel format.
References
-
Tetraploid state induces p53-dependent arrest of nontransformed mammalian cells in G1Molecular Biology of the Cell 12:1315–1328.https://doi.org/10.1091/mbc.12.5.1315
-
Perspectives on polyploidy in plants - ancient and neoBiological Journal of the Linnean Society 82:411–423.https://doi.org/10.1111/j.1095-8312.2004.00328.x
-
The incidence of adenocarcinoma in columnar-lined (Barrett's) esophagusNew England Journal of Medicine 313:857–859.https://doi.org/10.1056/NEJM198510033131404
-
A clinical overview of centrosome amplification in human cancersInternational Journal of Biological Sciences 7:1122–1144.https://doi.org/10.7150/ijbs.7.1122
-
Transient endoreplication down-regulates the kinesin-14 HSET and contributes to genomic instabilityMolecular Biology of the Cell 27:2911–2923.https://doi.org/10.1091/mbc.E16-03-0159
-
The advantages and disadvantages of being polyploidNature Reviews Genetics 6:836–846.https://doi.org/10.1038/nrg1711
-
The causes and consequences of polyploidy in normal development and CancerAnnual Review of Cell and Developmental Biology 27:585–610.https://doi.org/10.1146/annurev-cellbio-092910-154234
-
Loss of APC induces polyploidy as a result of a combination of defects in mitosis and apoptosisJournal of Cell Biology 176:183–195.https://doi.org/10.1083/jcb.200610099
-
The PIDDosome activates p53 in response to supernumerary centrosomesGenes & Development 31:34–45.https://doi.org/10.1101/gad.289728.116
-
Tetraploidy, aneuploidy and cancerCurrent Opinion in Genetics & Development 17:157–162.https://doi.org/10.1016/j.gde.2007.02.011
-
Causes and consequences of centrosome abnormalities in CancerPhilosophical Transactions of the Royal Society B: Biological Sciences 369:20130467.https://doi.org/10.1098/rstb.2013.0467
-
Barrett's esophagus: development of dysplasia and adenocarcinomaGastroenterology 96:1249–1256.https://doi.org/10.1016/S0016-5085(89)80011-3
-
Long-term complications to reflux disease in community practice A 17-year cohort study of 4706 patientsScandinavian Journal of Gastroenterology 46:1179–1186.https://doi.org/10.3109/00365521.2011.594084
-
Polyploid giant cells provide a survival mechanism for p53 mutant cells after DNA damageCell Biology International 24:621–633.https://doi.org/10.1006/cbir.2000.0557
-
Repeated cleavage failure does not establish centrosome amplification in untransformed human cellsThe Journal of Cell Biology 194:199–207.https://doi.org/10.1083/jcb.201101073
-
Mechanisms to suppress multipolar divisions in Cancer cells with extra centrosomesGenes & Development 22:2189–2203.https://doi.org/10.1101/gad.1700908
-
Cell fusions in mammalsHistochemistry and Cell Biology 129:551–561.https://doi.org/10.1007/s00418-008-0411-1
-
Proteins required for centrosome clustering in Cancer cellsScience Translational Medicine 2:33ra38.https://doi.org/10.1126/scitranslmed.3000915
-
Centrosome amplification arises before neoplasia and increases upon p53 loss in tumorigenesisJournal of Cell Biology 217:2353–2363.https://doi.org/10.1083/jcb.201711191
-
Polyploidy in differentiation and evolutionThe International Journal of Cell Cloning 8:216–223.https://doi.org/10.1002/stem.5530080404
-
Deregulated Aurora-B induced tetraploidy promotes tumorigenesisThe FASEB Journal 23:2741–2748.https://doi.org/10.1096/fj.09-130963
-
When bigger is better: the role of polyploidy in organogenesisTrends in Genetics 31:307–315.https://doi.org/10.1016/j.tig.2015.03.011
-
Polyploid incidence and evolutionAnnual Review of Genetics 34:401–437.https://doi.org/10.1146/annurev.genet.34.1.401
-
Polyploidy in tissue homeostasis and regenerationDevelopment 145:dev156034.https://doi.org/10.1242/dev.156034
-
Centrosome defects and genetic instability in malignant tumorsCancer Research 58:3974–3985.
-
Loss of E-cadherin provides tolerance to centrosome amplification in epithelial Cancer cellsJournal of Cell Biology 217:195–209.https://doi.org/10.1083/jcb.201704102
-
Centrosome abnormalities in pancreatic ductal carcinomaClinical Cancer Research 5:963–970.
-
Fiji: an open-source platform for biological-image analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019
-
The expanding implications of polyploidyJournal of Cell Biology 209:485–491.https://doi.org/10.1083/jcb.201502016
-
Pericentriolar material analyses in normal esophageal mucosa, Barrett's metaplasia and adenocarcinomaHistology and Histopathology 25:551–560.https://doi.org/10.14670/HH-25.551
-
Transient PLK4 overexpression accelerates tumorigenesis in p53-deficient epidermisNature Cell Biology 18:100–110.https://doi.org/10.1038/ncb3270
-
From polyploidy to Aneuploidy, genome instability and CancerNature Reviews Molecular Cell Biology 5:45–54.https://doi.org/10.1038/nrm1276
-
Cytokinetic Failure-induced tetraploidy develops into aneuploidy, triggering skin aging in Phosphovimentin-deficient miceJournal of Biological Chemistry 290:12984–12998.https://doi.org/10.1074/jbc.M114.633891
-
Cancer: a CINful evolutionCurrent Opinion in Cell Biology 52:136–144.https://doi.org/10.1016/j.ceb.2018.03.007
-
Midzone microtubule bundles are continuously required for cytokinesis in cultured epithelial cellsThe Journal of Cell Biology 135:981–989.https://doi.org/10.1083/jcb.135.4.981
-
Yesterday's polyploids and the mystery of diploidizationNature Reviews Genetics 2:333–341.https://doi.org/10.1038/35072009
-
Pan-cancer patterns of somatic copy number alterationNature Genetics 45:1134–1140.https://doi.org/10.1038/ng.2760
Article and author information
Author details
Funding
Virginia Tech (College of Science Dean's Discovery Fund)
- Daniela Cimini
Virginia Tech (Fralin Life Sciences Institute Discretionary funds)
- Daniela Cimini
Virginia Tech (ICTAS Center for Engineered Health Seed funding)
- Daniela Cimini
National Science Foundation (MCB-1517506)
- Daniela Cimini
Virginia Tech (BIOTRANS IGEP)
- Nicolaas C Baudoin
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We would like to thank the labs of Neil Ganem and Tim Stearns for reagents and Meng-Fu Bryan Tsou for the RPE-1 p53-/- cell line. We further acknowledge members of the Cimini, Chen, and Hauf labs for helpful discussion. Finally, we would like to thank Dr. Silke Hauf for critical reading of the manuscript.
Copyright
© 2020, Baudoin et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 3,021
- views
-
- 500
- downloads
-
- 29
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cell Biology
- Chromosomes and Gene Expression
During oncogene-induced senescence there are striking changes in the organisation of heterochromatin in the nucleus. This is accompanied by activation of a pro-inflammatory gene expression programme – the senescence-associated secretory phenotype (SASP) – driven by transcription factors such as NF-κB. The relationship between heterochromatin re-organisation and the SASP has been unclear. Here, we show that TPR, a protein of the nuclear pore complex basket required for heterochromatin re-organisation during senescence, is also required for the very early activation of NF-κB signalling during the stress-response phase of oncogene-induced senescence. This is prior to activation of the SASP and occurs without affecting NF-κB nuclear import. We show that TPR is required for the activation of innate immune signalling at these early stages of senescence and we link this to the formation of heterochromatin-enriched cytoplasmic chromatin fragments thought to bleb off from the nuclear periphery. We show that HMGA1 is also required for cytoplasmic chromatin fragment formation. Together these data suggest that re-organisation of heterochromatin is involved in altered structural integrity of the nuclear periphery during senescence, and that this can lead to activation of cytoplasmic nucleic acid sensing, NF-κB signalling, and activation of the SASP.
-
- Cell Biology
Protein aggregation increases during aging and is a pathological hallmark of many age-related diseases. Protein homeostasis (proteostasis) depends on a core network of factors directly influencing protein production, folding, trafficking, and degradation. Cellular proteostasis also depends on the overall composition of the proteome and numerous environmental variables. Modulating this cellular proteostasis state can influence the stability of multiple endogenous proteins, yet the factors contributing to this state remain incompletely characterized. Here, we performed genome-wide CRISPRi screens to elucidate the modulators of proteostasis state in mammalian cells, using a fluorescent dye to monitor endogenous protein aggregation. These screens identified known components of the proteostasis network and uncovered a novel link between protein and lipid homeostasis. Increasing lipid uptake and/or disrupting lipid metabolism promotes the accumulation of sphingomyelins and cholesterol esters and drives the formation of detergent-insoluble protein aggregates at the lysosome. Proteome profiling of lysosomes revealed ESCRT accumulation, suggesting disruption of ESCRT disassembly, lysosomal membrane repair, and microautophagy. Lipid dysregulation leads to lysosomal membrane permeabilization but does not otherwise impact fundamental aspects of lysosomal and proteasomal functions. Together, these results demonstrate that lipid dysregulation disrupts ESCRT function and impairs proteostasis.