Literature DB >> 30111539

Control of the Restriction Point by Rb and p21.

Justin Moser1,2, Iain Miller1,2, Dylan Carter1,2, Sabrina L Spencer3,2.   

Abstract

The Restriction Point was originally defined as the moment that cells commit to the cell cycle and was later suggested to coincide with hyperphosphorylation of the retinoblastoma protein (Rb). Current cell cycle models posit that cells exit mitosis into a pre-Restriction Point state, where they have low cyclin-dependent kinase (CDK) activity and hypophosphorylated Rb; passage through the Restriction Point then occurs in late G1. Recent single-cell studies have challenged the current paradigm, raising questions about the location of the Restriction Point and the notion that cells exit mitosis into a pre-Restriction Point state. Here, we use a variety of single-cell techniques to show that both noncancer and cancer cells bifurcate into two subpopulations after anaphase, marked by increasing vs. low CDK2 activity and hyper- vs. hypophosphorylation of Rb. Notably, subpopulations with hyper- and hypophosphorylated Rb are present within minutes after anaphase, delineating one subpopulation that never "uncrosses" the Restriction Point and continues cycling and another subpopulation that exits mitosis into an uncommitted pre-Restriction Point state. We further show that the CDK inhibitor p21 begins rising in G2 in mother cells whose daughters exit mitosis into the pre-Restriction Point, CDK2low state. Furthermore, degradation of p21 coincides with escape from the CDK2low state and passage through the Restriction Point. Together, these data support a model in which only a subset of cells returns to a pre-Restriction Point state after mitosis and where the Restriction Point is sensitive to not only mitogens, but also inherited DNA replication stress via p21.
Copyright © 2018 the Author(s). Published by PNAS.

Entities:  

Keywords:  CDK2; G0; cell cycle; quiescence; restriction point

Mesh:

Substances:

Year:  2018        PMID: 30111539      PMCID: PMC6126733          DOI: 10.1073/pnas.1722446115

Source DB:  PubMed          Journal:  Proc Natl Acad Sci U S A        ISSN: 0027-8424            Impact factor:   11.205


The ability of cells to transition between proliferative and quiescent states is critical for organismal health, as this enables tissue development and maintenance while preventing cancer (1). To commit to the cell cycle and proliferate, cells must cross the Restriction Point, after which they will complete the current cell cycle, even if serum or mitogens are withdrawn. Early serum withdrawal–restimulation–withdrawal experiments in synchronized cells suggested that cells arrest at a single point between mitosis and S phase until serum and mitogen conditions become favorable to proliferation once again (2). It was later shown that pulsed, rather than constant, mitogen exposure was sufficient to cross the Restriction Point and commit cells to a round of proliferation (3, 4). Similarly, time-lapse microscopy of asynchronously cycling Swiss 3T3 cells suggested that cycling cells are sensitive to serum withdrawal for only the first 3–4 h after mitosis, putting the Restriction Point in this interval (5). Based on these data, a model emerged that both cycling cells and cells emerging from serum starvation were subject to a mid- to late-G1 Restriction Point. Cells pre-Restriction Point are uncommitted to the cell cycle and can arrest at the Restriction Point, whereas cells post-Restriction Point are no longer dependent on mitogens and will complete one round of division, even in the absence of mitogens. Molecular biological and biochemical investigations later uncovered the molecular and systems-level basis for the Restriction Point, with the retinoblastoma protein (Rb) and cyclin:cyclin-dependent kinase (CDK) complexes coming to the fore. Based on work in cells emerging from serum starvation, serum restimulation results in the buildup of cyclin D (6), which in complex with CDK4/6, initiates the process of Rb phosphorylation (7). In the canonical model, this liberates some E2F, which initiates transcription of cyclins E and A. These cyclins in complex with CDK2 help generate the positive feedback loop that triggers the switch from hypo- to hyperphosphorylated Rb (8), marking cell cycle commitment (9). This fully releases E2F and leads to the production of other genes critical for S-phase entry. Based on work in cells synchronized in mitosis by nocodazole and subsequently released, the switch from hyperphosphorylated to hypophosphorylated Rb was shown to begin in late anaphase and continue through early G1 (10). The activity of CDKs, in particular CDK2, then triggers the switch from hypo- to hyperphosphorylated Rb at the Restriction Point and underlies the bistability of this system (7, 11, 12). Thus, activation of CDK2 and hyperphosphorylation of Rb indicate passage through the Restriction Point. The synthesis of these observations led to a model of the cell cycle in which cells are born into an uncommitted state characterized by low CDK activity and hypophosphorylated Rb (Fig. 1) (13). On crossing the Restriction Point several hours after mitosis with the hyperphosphorylation of Rb, cells are committed to one round of the cell cycle, giving rise to two daughter cells again born into a state of low CDK activity and hypophosphorylated Rb.
Fig. 1.

Two models of the Restriction Point depicted using simulated data. (A) The canonical model of the cell cycle predicts that cells are born into a pre-Restriction Point state characterized by low CDK2 activity, hypophosphorylated Rb, mitogen sensitivity, and lack of cell cycle commitment. This state has been referred to variously as “early G1” and “G1pm” (5, 7). Then, 3–4 h after mitosis, cells reach the Restriction Point (R). Given optimal conditions, cells cross the Restriction Point and become committed to the cell cycle. This state has been referred to as “late G1” and “G1ps” (5, 7). (B) An alternative cell-cycle model that takes into account recent studies, in which cells integrate mitogen and stress signaling during a maternal window of signal integration that informs the G0–G1 decision that occurs after mitosis. One subpopulation is born with moderate CDK2 activity and hyperphosphorylated Rb, is insensitive to mitogen withdrawal, and is committed to the cell cycle, whereas the other subpopulation is born with low CDK2 activity and hypophosphorylated Rb and remains sensitive to mitogen withdrawal, similar to the cells diagrammed in A. The proliferative subpopulation thus never uncrosses the Restriction Point, whereas the CDK2low subpopulation is born into a pre-Restriction Point state that can be considered a transient G0. These cells emerge from the CDK2low state to reengage with the cell cycle by building up CDK2 activity and hyperphosphorylating Rb.

Two models of the Restriction Point depicted using simulated data. (A) The canonical model of the cell cycle predicts that cells are born into a pre-Restriction Point state characterized by low CDK2 activity, hypophosphorylated Rb, mitogen sensitivity, and lack of cell cycle commitment. This state has been referred to variously as “early G1” and “G1pm” (5, 7). Then, 3–4 h after mitosis, cells reach the Restriction Point (R). Given optimal conditions, cells cross the Restriction Point and become committed to the cell cycle. This state has been referred to as “late G1” and “G1ps” (5, 7). (B) An alternative cell-cycle model that takes into account recent studies, in which cells integrate mitogen and stress signaling during a maternal window of signal integration that informs the G0–G1 decision that occurs after mitosis. One subpopulation is born with moderate CDK2 activity and hyperphosphorylated Rb, is insensitive to mitogen withdrawal, and is committed to the cell cycle, whereas the other subpopulation is born with low CDK2 activity and hypophosphorylated Rb and remains sensitive to mitogen withdrawal, similar to the cells diagrammed in A. The proliferative subpopulation thus never uncrosses the Restriction Point, whereas the CDK2low subpopulation is born into a pre-Restriction Point state that can be considered a transient G0. These cells emerge from the CDK2low state to reengage with the cell cycle by building up CDK2 activity and hyperphosphorylating Rb. Although this model of the cell cycle is widely accepted, recent studies in asynchronously cycling cells have challenged this paradigm. The development of a CDK2 activity sensor revealed two classes of cellular behavior after mitosis in the mammary epithelial cell line MCF10A (14). Some cells follow the expected paradigm, wherein after mitosis, CDK2 activity turns off, Rb becomes dephosphorylated, and cells are once again mitogen sensitive (“CDK2low cells”). CDK2low cells are not senescent, and cells can remain in this state anywhere from a few hours to multiple days before reentering the cell cycle (15). While the CDK2low state could be conceptualized as “early G1,” it also has many features consistent with a transient G0 or quiescence, including declining Ki67 protein levels (16), little to no CDK2 activity, hypophosphorylated Rb, and elevated p21 (14). In contrast, other cells show increasing CDK2 activity after mitosis, hyperphosphorylated Rb, and mitogen insensitivity, and they immediately commit to another cell cycle (“CDK2inc cells”) (14). Moreover, the CDK2 activity sensor can be used to visualize passage through the Restriction Point: cells in the CDK2low state are mitogen sensitive and become locked in the CDK2low state if mitogens are withdrawn, but after cells build up CDK2 activity above a threshold level, they become mitogen insensitive and continue progression through the cell cycle, even if mitogens are withdrawn or the MAPK pathway is inhibited (14, 17, 18). Recent papers suggest that entry into quiescence is facilitated during the maternal G2 (14, 19) and that spontaneous entry into the CDK2low state is caused in part by low-level endogenous DNA replication stress that is passed from mother cells to daughter cells (20–23). Furthermore, the CDK2 inhibitor, p21, was found to be up-regulated in these stressed cells and was required for entry into the CDK2low state after both endogenous and exogenous replication stress (20, 21). Thus, p21 seems to be a critical component of the Restriction Point decision, but its exact connection to the Restriction Point is unclear. Based on these recent findings, we propose here that healthy cells in optimal growth conditions cycle continuously without ever “uncrossing” the Restriction Point or losing Rb phosphorylation (Fig. 1), whereas the classic Restriction Point model predicts that cells are born into a pre-Restriction Point state with hypophosphorylated Rb and low CDK2 activity (Fig. 1). Distinguishing between these two models and establishing their generalizability are necessary for improving our understanding of how cells behave both in cell culture and in vivo. The canonical model of the cell cycle ignores the behavioral heterogeneity as well as its underlying causes that have been unearthed in recent years. It thus incompletely acknowledges the various mechanisms and timing thereof that regulate cell cycle commitment, exit, and reentry. In this manuscript, we show the applicability of the alternative model across five widely used transformed and nontransformed cell lines and one type of primary cell, all of which have functional Rb and p53, establishing its validity beyond MCF10A cells. We show that all of these cell lines display a bifurcation in CDK2 activity and Rb phosphorylation at mitotic exit. This result contrasts with the canonical cell cycle model, in which cells have hypophosphorylated Rb after anaphase. Furthermore, by tracking p21 under endogenous control in single cells over multiple cell cycles, we detect the up-regulation of p21 in G2 phase of mother cells whose daughters enter the CDK2low state after mitosis. We further show that these CDK2low cells reenter the cell cycle by degrading p21 at the Restriction Point. Together, these data reveal that only a fraction of cells in the cell types examined here display features of the canonical model and exit mitosis into a pre-Restriction Point state of hypophosphorylated Rb. Our results instead support an alternative model in which, under saturating mitogens, the proliferation–quiescence decision is influenced strongly by the G2/M levels of p21, which if high enough, result in the dephosphorylation of Rb in telophase, and the reset of those cells to a pre-Restriction Point state.

Results

A Bifurcation in CDK2 Activity and Rb Phosphorylation After Mitosis Is Present in both Transformed and Nontransformed Cell Types.

The observation of a bifurcation in CDK2 activity after mitosis in unperturbed cells was initially made in nontransformed MCF10A mammary epithelial cells (14). To determine the generalizability of this finding, we examined several noncancerous (MCF10A and retinal RPE-hTERT) and cancerous (mammary MCF7, osteosarcoma U2OS, and colorectal HCT116) cells as well as primary human lung fibroblasts (HLFs). HLFs were recently reported to follow the canonical model depicted in Fig. 1 rather than the alternative model in Fig. 1 (17). We transduced all six cell types with fluorescent histone 2B (H2B) as a nuclear marker as well as the CDK2 activity sensor and monitored CDK2 activity by time-lapse imaging and cell tracking in asynchronously cycling cells. Despite the wide variety of cell types examined, each showed a bifurcation in CDK2 activity at mitotic exit (Fig. 2 and Movies S1–S6). Notably, the cancer lines did not always show a lower fraction of CDK2low cells than MCF10A. Only HCT116 cells, with 1% CDK2low, were significantly more proliferative than MCF10A cells, with 26% CDK2low, whereas MCF7 and U2OS cells were intermediate with 18 and 17% CDK2low, respectively. HLF cells were comparable with these, with 23% of cells entering the CDK2low state, whereas RPE-hTERT cells were the least likely to immediately enter another cell cycle, with 31% entering the CDK2low state after mitosis under optimal growth conditions.
Fig. 2.

Evidence for the generalizability of the bifurcation in CDK2 activity, Rb phosphorylation, and p21 expression. (A) Single-cell traces of CDK2 activity for CDK2inc cells (blue) and CDK2low cells (red). For clarity, 250 manually verified traces are plotted for each cell line of thousands analyzed; the numbers of CDK2inc and CDK2low traces plotted here are proportional to the fractions of each at the population level. CDK2low: MCF10A, 26%; RPE-hTERT, 31%; HLF, 23%; MCF7, 18%; U2OS, 17%; and HCT116, 1%. (B) Single-cell measurement of phospho-Rb (S807/811) as a function of time-since-anaphase for all cells analyzed. Cells were imaged for 24 h, fixed, and stained for phospho-Rb, and the immunofluorescence image was registered to the final frame of the movie. MCF10A, n = 15,446 cells; RPE-hTERT, n = 11,936 cells; HLF, n = 547 cells; MCF7, n = 7,659 cells; U2OS, n = 5,348 cells; and HCT116, n = 8,157 cells. (C) Single-cell measurement of p21 as a function of time-since-anaphase as in B. MCF10A, n = 19,052 cells; RPE-hTERT, n = 2,809 cells; HLF, n = 1,019 cells; MCF7, n = 1,929 cells; U2OS, n = 4,958 cells; and HCT116, n = 825 cells.

Evidence for the generalizability of the bifurcation in CDK2 activity, Rb phosphorylation, and p21 expression. (A) Single-cell traces of CDK2 activity for CDK2inc cells (blue) and CDK2low cells (red). For clarity, 250 manually verified traces are plotted for each cell line of thousands analyzed; the numbers of CDK2inc and CDK2low traces plotted here are proportional to the fractions of each at the population level. CDK2low: MCF10A, 26%; RPE-hTERT, 31%; HLF, 23%; MCF7, 18%; U2OS, 17%; and HCT116, 1%. (B) Single-cell measurement of phospho-Rb (S807/811) as a function of time-since-anaphase for all cells analyzed. Cells were imaged for 24 h, fixed, and stained for phospho-Rb, and the immunofluorescence image was registered to the final frame of the movie. MCF10A, n = 15,446 cells; RPE-hTERT, n = 11,936 cells; HLF, n = 547 cells; MCF7, n = 7,659 cells; U2OS, n = 5,348 cells; and HCT116, n = 8,157 cells. (C) Single-cell measurement of p21 as a function of time-since-anaphase as in B. MCF10A, n = 19,052 cells; RPE-hTERT, n = 2,809 cells; HLF, n = 1,019 cells; MCF7, n = 1,929 cells; U2OS, n = 4,958 cells; and HCT116, n = 825 cells. We next examined the phosphorylation state of Rb in these six cell types. Not only is Rb a canonical CDK2 substrate, but the phosphorylation status of this protein is considered an indicator of a cell’s position relative to the Restriction Point (18, 24, 25). After 24 h of time-lapse imaging, we fixed and stained cells for phospho-Rb at Serine 807/811. Phospho-Rb is bimodally distributed (), and we consider the high mode to represent hyperphosphorylated Rb as was done previously (14, 23, 26). Through image registration, we linked the immunofluorescence data to the time-lapse data to determine how long ago each cell underwent anaphase (Fig. 2 and ). Phospho-Rb staining already shows clear bimodality one frame after anaphase (12 or 20 min depending on cell type) in all six cell types examined, indicating that the proliferative subpopulation maintains hyperphosphorylation of Rb from one cycle to the next as has previously been shown for MCF10A (26). Strikingly, the fraction of cells with hyperphosphorylated Rb immediately after anaphase is similar to the fraction of the population born into the CDK2inc state (). We then repeated this experiment, staining for p21 after 24 h of time-lapse imaging. A weakly bimodal or long-tailed distribution of p21 abundance is apparent within the first frame after anaphase (Fig. 2 and ). Both the hypophosphorylated Rb state and the high p21 population decrease as time-since-anaphase increases, highlighting the transient nature of the CDK2low state. Previous work by this group and others has shown that in MCF10A, RPE-hTERT, and Hs68 primary human fibroblasts, high p21 levels in G0/G1 cells are associated with 53BP1 nuclear bodies, a marker of DNA lesions (20, 21). We therefore asked if this relationship was conserved in the cancer cell lines used here. As was previously shown with MCF10A and RPE-hTERT, we find that G0/G1 cells with 53BP1 nuclear bodies have higher p21 levels than G0/G1 cells without 53BP1 nuclear bodies (). Thus, the three cancer cell lines tested here seem to face the CDK2inc vs. CDK2low decision in the same way as reported previously for nontransformed cells, although the fraction of cells born committed to each fate varies across cell lines.

p21 Up-Regulation and Rb Dephosphorylation Are Temporally Separated.

Given that asynchronously cycling cells already show differential Rb phosphorylation and p21 levels immediately after mitosis, we next asked when exactly this bifurcation in Rb status and p21 expression occurs. To this end, we fixed and stained asynchronously cycling, wild-type MCF10A cells for the mitotic marker phospho-histone H3 (pHH3) and either phospho-Rb or p21 and measured these three signals by flow cytometry. p21 was up-regulated in 30% of pHH3+ cells, mirroring the fraction of cells born into the CDK2low state (Fig. 3, upper right quadrant). In contrast, there was no discrete hypophosphorylated Rb population among pHH3+ cells. We therefore asked whether this unexpected result could be explained by the facts that phosphorylation of histone H3 is lost before anaphase and that phosphorylation of Rb is lost in CDK2low cells after anaphase.
Fig. 3.

A subset of cells matching the CDK2low population up-regulates p21 before mitosis but does not dephosphorylate Rb until chromatin decondensation begins. (A and B) p21 signal (A) or phospho-Rb (pRb; S807/811) signal (B) as a function of pHH3 signal in asynchronously cycling MCF10A cells as measured by flow cytometry. Thirty-three percent of mitotic (pHH3+) cells have high p21 levels; 12% of mitotic cells have falling pRb, but there is no discrete hypophosphorylated Rb population (n = 10,000 cells). Gates were set according to the saddle point for p21 and pRb. (C) Asynchronously growing MCF10A cells were stained for pHH3, pRb (S807/811), and Hoechst. Mitotic cells were identified based on Hoechst staining and displayed as mock film strips. Anaphase cells are ordered based on distance between chromosomes. Histone H3 is dephosphorylated starting in late anaphase (green stars), whereas Rb is not dephosphorylated until after anaphase is complete and chromatin decondensation begins. (D) pRb (S807/811) in MCF10A cells from Fig. 2. Shown are newly born cells fixed 12 min (Upper) or 24 min (Lower) after anaphase was detected. Cells are ordered based on chromatin decondensation and distance between sister cells. The numbers of cells with hyper- and hypophosphorylated Rb presented here are proportional to the total population. Green stars mark cells that are pHH3−; magenta stars mark cells with hypophosphorylated Rb.

A subset of cells matching the CDK2low population up-regulates p21 before mitosis but does not dephosphorylate Rb until chromatin decondensation begins. (A and B) p21 signal (A) or phospho-Rb (pRb; S807/811) signal (B) as a function of pHH3 signal in asynchronously cycling MCF10A cells as measured by flow cytometry. Thirty-three percent of mitotic (pHH3+) cells have high p21 levels; 12% of mitotic cells have falling pRb, but there is no discrete hypophosphorylated Rb population (n = 10,000 cells). Gates were set according to the saddle point for p21 and pRb. (C) Asynchronously growing MCF10A cells were stained for pHH3, pRb (S807/811), and Hoechst. Mitotic cells were identified based on Hoechst staining and displayed as mock film strips. Anaphase cells are ordered based on distance between chromosomes. Histone H3 is dephosphorylated starting in late anaphase (green stars), whereas Rb is not dephosphorylated until after anaphase is complete and chromatin decondensation begins. (D) pRb (S807/811) in MCF10A cells from Fig. 2. Shown are newly born cells fixed 12 min (Upper) or 24 min (Lower) after anaphase was detected. Cells are ordered based on chromatin decondensation and distance between sister cells. The numbers of cells with hyper- and hypophosphorylated Rb presented here are proportional to the total population. Green stars mark cells that are pHH3−; magenta stars mark cells with hypophosphorylated Rb. To more carefully examine the timing of dephosphorylation of Rb, we repeated the experiment using immunofluorescence microscopy, computationally isolated cells at each phase of mitosis based on Hoechst staining, and then evaluated the levels of pHH3 and phospho-Rb in the resulting mock film strips. Histone H3 is dephosphorylated in late anaphase, in agreement with previous literature (27–30), whereas Rb is not dephosphorylated until after late anaphase (Fig. 3 ), explaining why no discrete pHH3+/hypophosphorylated Rb population is detectable by flow cytometry. Synthesizing our flow cytometry, film strip, and time-lapse microscopy data (where we have immunofluorescence data for phospho-Rb beginning 12 min after anaphase), we conclude that Rb dephosphorylation in CDK2low cells occurs in telophase after chromosome decondensation has started. The differential timing of p21 up-regulation in G2 and Rb dephosphorylation in telophase suggests that p21 up-regulation (by inhibiting CDK activity) causes entry into the CDK2low state, whereas Rb dephosphorylation (due to mitotic phosphatase activity and a dearth of CDK activity as CDK2low cells exit mitosis) is a consequence.

p21 Up-Regulation Precedes Mitosis in Mother Cells Whose Daughters Enter the CDK2low State After Mitosis.

The similarity between the fraction of cells exiting the cell cycle after anaphase and the fraction of mitotic cells with high p21 prompted us to ask when exactly the up-regulation of p21 occurs. Given the importance of p21 in controlling the bifurcation in CDK2 activity (14, 20, 21, 31), we used CRISPR-Cas9 to engineer MCF10A that express mCitrine-tagged p21 from the endogenous CDKN1A locus (). Although some studies have highlighted the potential for N-terminal acylation or ubiquitylation as a regulatory mechanism for p21 (32, 33), we targeted the N terminus of p21 because similarly tagged constructs have been shown to maintain functional p21 (14, 31). p21 tagged at the N terminus localizes normally, degrades similarly to wild-type p21, can still interact with CDK2, and does not dramatically alter cell cycle dynamics (). We also examined three other cell lines in which p21 is tagged at the C terminus with a fluorescent protein at the endogenous locus [RPE-hTERT (21), MCF7, and HCT116 (34)] or expressed from an inducible promoter at endogenous levels [U2OS (35)]. Each cell line was transduced with a fluorescently tagged H2B nuclear marker and the CDK2 sensor and then imaged for 24 h (Movies S2–S6). Cells from each population were classified according to CDK2 activity after anaphase: CDK2inc, CDK2low, and cells that emerge from the CDK2low state 4–7 or 7–10 h after anaphase (CDK2emerge4–7 h and CDK2emerge7–10 h, respectively). The single-cell CDK2 and p21 traces were then averaged within these four groups and aligned to the time of anaphase (Fig. 4 ).
Fig. 4.

p21 up-regulation precedes mitosis in mother cells whose daughters enter the CDK2low state after mitosis, and p21 is degraded at the Restriction Point. (A–C) Cells expressing the CDK2 sensor and fluorescent p21 from the endogenous locus (with the exception of U2OS, which were induced to endogenous levels with doxycycline) were imaged for 24 h, tracked, classified according to CDK2 activity after anaphase, averaged, and aligned to either the time of anaphase (A and B) or the Restriction Point (C), which was defined as the time that CDK2 activity begins to rise. (A) CDK2 activity traces were classified as CDK2low if the CDK2 activity dropped below the cutoff indicated in and remained below this cutoff for the remainder of the imaging period. Traces were classified as CDK2emerge 4–7 h or CDK2emerge 7–10 h if the CDK2 activity dropped below the cutoff indicated in and rose above this cutoff 4–7 or 7–10 h after anaphase, respectively; otherwise, traces were classified as CDK2inc. A CDK2emerge population does not exist for HCT116 cells because few cells enter the CDK2low state. (B) Averaged p21 traces for the categories described in A. The purple line above each plot indicates the hypothesis testing result of a two-sample t test with P value < 0.05, comparing p21 in CDK2inc with CDK2low cells at each time point. (C) Averaged p21 and CDK2 activity traces for CDK2emerge7–10 h cells aligned to the Restriction Point. (D) Averaged p21 (orange) and CDK2 activity (black) traces from MCF10A CDK2emerge4–7 h cells aligned to the Restriction Point. Only those cells receiving DMSO (Left) or 1.4 M MLN4924 (Right) 3–4 h (gray shading) after CDK2 had begun to rise are included. (E) Averaged p21 (orange) and CDK2 activity (black) traces from MCF10A CDK2low cells that received DMSO (Left) or 1.4 M MLN4924 (Right) 3–5 h (gray shading) after anaphase. Error bars represent SEM. Data presented here are representative of two biological replicates. Cell counts for each group can be found in . mCit-p21, mCitrine-tagged p21 from the endogenous CDKN1A locus. p21 axis units are in log(10).

p21 up-regulation precedes mitosis in mother cells whose daughters enter the CDK2low state after mitosis, and p21 is degraded at the Restriction Point. (A–C) Cells expressing the CDK2 sensor and fluorescent p21 from the endogenous locus (with the exception of U2OS, which were induced to endogenous levels with doxycycline) were imaged for 24 h, tracked, classified according to CDK2 activity after anaphase, averaged, and aligned to either the time of anaphase (A and B) or the Restriction Point (C), which was defined as the time that CDK2 activity begins to rise. (A) CDK2 activity traces were classified as CDK2low if the CDK2 activity dropped below the cutoff indicated in and remained below this cutoff for the remainder of the imaging period. Traces were classified as CDK2emerge 4–7 h or CDK2emerge 7–10 h if the CDK2 activity dropped below the cutoff indicated in and rose above this cutoff 4–7 or 7–10 h after anaphase, respectively; otherwise, traces were classified as CDK2inc. A CDK2emerge population does not exist for HCT116 cells because few cells enter the CDK2low state. (B) Averaged p21 traces for the categories described in A. The purple line above each plot indicates the hypothesis testing result of a two-sample t test with P value < 0.05, comparing p21 in CDK2inc with CDK2low cells at each time point. (C) Averaged p21 and CDK2 activity traces for CDK2emerge7–10 h cells aligned to the Restriction Point. (D) Averaged p21 (orange) and CDK2 activity (black) traces from MCF10A CDK2emerge4–7 h cells aligned to the Restriction Point. Only those cells receiving DMSO (Left) or 1.4 M MLN4924 (Right) 3–4 h (gray shading) after CDK2 had begun to rise are included. (E) Averaged p21 (orange) and CDK2 activity (black) traces from MCF10A CDK2low cells that received DMSO (Left) or 1.4 M MLN4924 (Right) 3–5 h (gray shading) after anaphase. Error bars represent SEM. Data presented here are representative of two biological replicates. Cell counts for each group can be found in . mCit-p21, mCitrine-tagged p21 from the endogenous CDKN1A locus. p21 axis units are in log(10). Contrary to early models of cell cycle-dependent p21 expression (36, 37), we find that p21 up-regulation is not a general feature of G2. Instead, daughter cells that enter the CDK2inc state after mitosis maintain low levels of p21 in the previous G2 and M, while daughter cells that enter the CDK2low state after mitosis start up-regulating p21 5–10 h before anaphase, depending on the cell line (Fig. 4). These CDK2low daughter cells then continue to increase p21 levels after anaphase, sustaining the CDK2low state. In contrast, CDK2emerge cells that initially enter the CDK2low state and then reenter the cell cycle show a decline in p21 levels around the time of cell cycle reentry.

p21 Degradation Is Initiated at the Restriction Point.

To determine when CDK2emerge cells down-regulate p21 relative to when they reactivate CDK2 and cross the Restriction Point, we aligned and averaged single-cell CDK2 activity traces from CDK2emerge cells to the time of the rise in CDK2 activity (the Restriction Point) and monitored the levels of p21. While there is some cell-to-cell heterogeneity in the timing of p21 degradation relative to the timing of CDK2 activation (), our analysis revealed that, on average, p21 continues to accumulate in newly born CDK2low cells until the Restriction Point, at which point p21 levels fall dramatically (Fig. 4 and ). The near-simultaneous nature of rising CDK2 activity and falling p21 is consistent with the notion that CDK2 activity promotes the degradation of p21, an idea that has been predicted by analogy to p27 (38, 39) but cannot as of this writing be directly tested due to a lack of selective CDK2 inhibitors. To test instead the hypothesis that the degradation of p21 begins at the Restriction Point, we treated asynchronously cycling MCF10A cells with MLN4924, which inhibits Cullin-based E3 ligase complexes, and selected for analysis only CDK2emerge4–7 h cells that received the drug 3–4 h after the Restriction Point (Fig. 4). Whereas cells receiving vehicle continue to down-regulate p21 after crossing the Restriction Point, cells receiving MLN4924 rapidly reaccumulate p21. In contrast, p21 levels do not deviate from their increasing trajectory in CDK2low cells on treatment with MLN4924 (Fig. 4). We conclude that CDK2low cells do not actively degrade p21 and that degradation of p21 begins coincident with the rise in CDK2 activity at the Restriction Point.

Discussion and Conclusions

A long-standing model of the cell cycle suggests that cells are born into a pre-Restriction Point state in which they are uncommitted to proliferation. For the first few hours after anaphase, cells are thought to integrate environmental signals to determine if they can cross the Restriction Point. After they cross this point, they are committed to one round of the cell cycle, and the resulting daughter cells are again born into an uncommitted pre-Restriction Point state. The groundbreaking studies that established this model relied predominately on cell cycle synchronization and bulk population analysis, which perturb the cell cycle and mask heterogeneity in cell behavior. The rise of single-cell analysis has challenged aspects of this model, suggesting instead that, in actively cycling cells, the uncommitted CDK2low state is sampled only by a subset of cells (14) that experienced stress (20–23, 40) or blockade of MAPK signaling (14, 23, 26, 41) during the previous cell cycle. In line with this recent trend, this study uses a combination of single-cell time-lapse imaging and fixed-cell analysis to show, across a number of primary, immortalized but not transformed, and cancerous cell types, that only a subset of cells in a population enters the uncommitted CDK2low state after mitosis. Furthermore, independent of the CDK2 sensor, this heterogeneity is visible by immunofluorescence staining of Rb phosphorylation and p21, where a subset of cells exits mitosis with hyperphosphorylated Rb and low p21, while the remainder has hypophosphorylated Rb and high p21. The conclusion that a subset of cells is born committed to proliferation is further supported by the observation that, when subjected to serum withdrawal or acute Mek inhibition, CDK2inc cells finish the current cell cycle, even if they are so perturbed in early G1 (14). Thus, immediately after anaphase, CDK2inc cells are already in a post-Restriction Point state. In contrast, CDK2low cells remain sensitive to serum withdrawal and Mek inhibition as long as they are in the CDK2low state. These cells become insensitive to these perturbations after CDK2 activity rises and they cross the Restriction Point (14, 17). Our results, therefore, argue that only a subset of cells exits mitosis into a pre-Restriction Point state. Contrary to recent work (17), we find this can be true even for primary HLFs. Additionally, the similarity between the fraction of the population with hyperphosphorylated Rb and the fraction CDK2inc suggests that the size of the proliferative subpopulation can be estimated from simple fixed cell immunofluorescence of phospho-Rb. Several lines of evidence suggest that a cause of entry into the pre-Restriction Point CDK2low state is high p21, including that p21−/− cells rarely enter the CDK2low state and that acute overexpression of p21 in G2 is sufficient to send all cells into the CDK2low state after mitosis (14, 20, 31). Here, we show that endogenous p21 begins to increase during G2 in mother cells whose daughters enter the CDK2low state after mitosis, whereas mothers of CDK2inc daughters have significantly less p21 both before and after anaphase. Notably, this pattern is not dependent on whether p21 is fused to a fluorescent protein at the N terminus (MCF10A, U2OS) or the C terminus (RPE-hTERT, MCF7, HCT116) or if the fluorescent p21 is expressed at endogenous levels off an inducible promoter (U2OS) or from the endogenous CDKN1A locus, suggesting the existence of both p53 transcription-dependent and p53 transcription-independent regulation of p21 in the proliferation–quiescence decision. The observation that p21 dynamics are independent of whether p21 expression is off the endogenous promoter or via a doxycycline-inducible promoter suggests that posttranslational modifications play an essential role in regulating abundance of this protein and by extension, cell fate. We note that tagged p21 seems to be more highly expressed than the wild type (), possibly due to the loss of N terminus-driven degradation (32, 33). Nevertheless, both N- and C-terminal tagged p21 dynamics are similar. Our results reveal that the graded input of p21 levels in G2/M is converted to a binary output, resulting in proliferative CDK2inc and transiently quiescent CDK2low daughter cells after mitosis. If replication stress at the end of S phase leads to incompletely replicated genomic loci (20, 21, 42, 43), why does this not trigger a G2/M checkpoint? Do cells simply not detect the problem? Our data suggest that, in response to low levels of endogenous replication stress, a cellular warning signal is indeed triggered in G2 in the form of up-regulation of the CDK inhibitor p21 but that the levels of p21 attained before mitosis are insufficient to block progression through mitosis due to the very high levels of CDK2 and CDK1 activity in cells during G2 and M. Cells therefore proceed through mitosis and then arrest at the start of the new cell cycle when CDK activity is low again and when p21 levels are sufficient to block CDK-mediated cell cycle progression. In general, arrest of daughter cells in a CDK2low state after mitosis is likely easier to sustain long term than a G2/M arrest, although recent work in Drosophila has suggested the possibility of a G2 quiescence, at least in stem cells (44). Thus, our data suggest that G2 phase in mother cells represents a window, referred to previously as R1 (14) or as the maternal window of signal integration here, where cells sense both mitogens and stress and initiate a response, which is then converted into a bifurcation in CDK2 activity after mitosis. For CDK2inc cells committed to proliferation, this window is closed by the start of the new cell cycle, whereas the window of signal integration remains open for CDK2low cells, allowing them to continue integrating mitogen and stress signals until they cross the Restriction Point and commit to a new cell cycle. These p21high/CDK2low cells can reenter the cell cycle by degrading p21 at the Restriction Point. Thus, p21 degradation reflects the decision to resume proliferation from the CDK2low state. This result is consistent with recent observations that p21 degradation can begin before S phase (21, 34) and contrasts with other models of p21 degradation, which hold that p21 degradation begins at the start of S phase (35–37). While the presence of the proliferating cell nuclear antigen (PCNA)-interacting protein (PIP) degron in p21 (35) and the potential for PCNA-dependent degradation of this protein by CRL4Cdt2 (45) are important for active degradation of p21 in S phase, it is possible that PCNA-dependent degradation is a redundant system to guarantee that p21 will not be present at high enough levels to impede S-phase progression (31, 42) and that a PCNA-independent degradation mechanism initiates p21 destruction at the Restriction Point, several hours before the start of S phase. The relationship between p21 degradation and passage through the Restriction Point further suggests an intimate connection between the two, wherein the Restriction Point is not only sensitive to mitogens but also to DNA damage via p21. Whereas the role of mitogens in Restriction Point control has been extensively studied, when and how cells down-regulate p21 to promote passage through the Restriction Point is an open area of investigation. Why does the fraction of cells entering the CDK2low state vary from one cell line to the next? This study and earlier studies have focused on the roles that DNA damage and endogenous p21 play in the proliferation–quiescence decision (20, 21). If endogenous DNA damage arises in S phase, any variation in DNA replication error rates between cell lines will tune the fraction of the population exiting the cell cycle after mitosis. Furthermore, recent work has suggested that both DNA repair efficiency and basal ataxia telangiectasia mutated/ataxia telangiectasia and Rad3-related protein (ATM/ATR) activity inform p53 dynamics (46); given that p53 is the main transcription factor of p21 (47), how the interplay between these factors influences p53 dynamics may in turn inform the fraction of a population that exits the cell cycle after anaphase. Other molecules may also play a role in regulating the proliferation–quiescence decision, such as cyclins D and E (15, 48), p27, p57, Wee1, Cdc25 (49, 50), PP2A (19), and USP11 (51). Alternate routes out of the cell cycle may explain the fraction of transits through the CDK2low state that are not explained by inherited DNA damage from mother cells (20). The up- or down-regulation of any one of these proteins in a given cell line would help determine the fraction of cells entering the CDK2low state after mitosis. Indeed, HCT116 cells may avoid cell-cycle exit via the overexpression of low-molecular weight cyclin E (52). Furthermore, the availability of mitogens is a key factor—under cell culture conditions with saturating growth factors as used here, cells are likely to achieve their proliferative potential. In this regime, the alternative cell cycle model depicted in Fig. 1 is clearly evident, with a majority of cells committing immediately to another cell cycle, and therefore, cells will be limited primarily by stress. Conversely, under more physiological conditions, in which mitogens are limiting, proliferation will be constrained by both mitogens and stress, and a higher fraction of cells will enter the CDK2low state after mitosis. Our model is predicated on cells having intact Rb and p53; cell lines without functional p53 may not be able to use p21 to enter the CDK2low state. Unless these cells have alternative means of suppressing CDK2 activity, we expect cells lacking p53 to cycle continuously, such as was shown previously in cells where p21 or p53 was depleted or eliminated (14, 20, 21). Our data show that the classic model of the cell cycle, in which cells are born into a pre-Restriction Point state with hypophosphorylated Rb, is not universally applicable, and in fact, it is not well-suited to any of the cell types examined here, including primary HLFs. Instead, all six cell types that we examined are split such that some cells reset to a pre-Restriction Point state, while others (the majority in the case of the six cell types examined here) are born without ever having “uncrossed” the Restriction Point and are thus immediately committed to another cell cycle. Put another way, healthy cells in optimal mitogenic conditions do not reset to a pre-Restriction point state and therefore continue cycling. Although all cell types examined here have the Restriction Point machinery intact, only CDK2low cells reset to a pre-Restriction Point state. In other words, whereas all or most cells seem to have a maternal window in which they integrate mitogenic and stress signaling, only daughter cells entering the CDK2low state experience the canonical Restriction Point in the daughter cell cycle. This then implies that the Restriction Point must be located at the end of the CDK2low period. Thus, rather than occurring specifically 3–4 h after mitosis, the Restriction Point is better conceptualized as the variable time point at which the cell decides to reenter the cell cycle, degrade p21, and build up CDK2 activity. Positioning the Restriction Point outside the cell cycle as in Fig. 1, where it is only ever reset in cells that enter the CDK2low state, also explains the dramatic heterogeneity in cell cycle duration across a given population. Residence time in the pre-Restriction Point CDK2low state is variable, and cell cycle duration varies according to the length of time that it takes for those cells to cross the Restriction Point. In summary, our data support the alternative model in Fig. 1, where the decision to proliferate is informed by events occurring at the end of the previous cell cycle and only a subset of cells resets to a pre-Restriction Point state after mitosis.

Materials and Methods

Cell Culture and Maintenance.

MCF10A (ATCC CRL-10317) were maintained in DMEM/F12 supplemented with 5% horse serum, 100 ng/mL cholera toxin, 20 ng/mL EGF, 10 g/mL insulin, 0.5 g/mL hydrocortisone, and 100 g/mL each penicillin and streptomycin. RPE-hTERT [ATCC CRL-4000 and those tagged with p21-GFP (21)] were maintained in DMEM/F12 supplemented with 10% FBS, 1× Glutamax, and 100 g/mL each penicillin and streptomycin. MCF7 were maintained in RPMI supplemented with 10% FBS, 1× Glutamax, and 100 g/mL each penicillin and streptomycin. U2OS and HCT116 were maintained in McCoy’s 5A supplemented with 10% FBS and 100 g/mL each penicillin and streptomycin. Primary fetal HLFs (507-75f; Cell Applications) were maintained in Human Lung Fibroblast Growth Medium (516-500; Cell Applications), in contrast with the growth conditions used in ref. 17. For live-cell imaging, each cell line was maintained in a phenol red-free version of their growth media; HLFs were imaged in Human Lung Fibroblast Growth Medium. U2OS were seeded with 10 ng/L doxycycline to induce endogenous-level expression of fluorescent p21 (35) for at least 24 h before imaging. All cell lines were grown in a humidified incubator (5% CO2, 37 ). MCF7 and HCT116 cells with p21 tagged at the endogenous locus were provided by the laboratory of Galit Lahav, Department of Systems Biology, Harvard Medical School, Boston (34). RPE-hTERT with p21 tagged at the endogenous locus was provided by the laboratory of Chris Bakal, Division of Cancer Biology, Institute of Cancer Research, London (21). U2OS cells with inducible p21 were provided by the laboratory of Jeanette Cook, Department of Biochemistry and Biophysics, University of North Carolina at Chapel Hill, Chapel Hill, NC (35).

Cell Line Generation.

Low-passage wild-type MCF10A were transfected with a plasmid encoding CRISPR-Cas9, single-guide (sg)RNA targeting the 5′ end of CDKN1A, and a repair template encoding the mCitrine gene; the left and right homology arms extended for 1,000 bp upstream and downstream of the CDKN1A start codon. The sgRNA sequence used was GCGCCATGTCAGAACCGGCTGGG. Cells were transfected in antibiotic-free medium and supplemented with 1 M SCR7, an inhibitor of nonhomologous end joining (SML1546; Sigma). Cells were allowed to recover and expand before the YFP-positive cells were isolated via FACS. Clonal lines were then expanded and validated by PCR, genomic DNA sequencing, immunofluorescence, and Western blotting (). Cells used for time-lapse live-cell microscopy were transduced with a nuclear marker [H2B-mTurquoise or H2B-miFP (53)] and the CDK2 sensor DHB (14) using established lentivirus protocols, and double-positive cells were sorted by FACS. CDK2 activity was read out as the cytoplasmic to nuclear ratio of the DHB sensor. Primary HLFs were transduced with H2B-mTurquoise and DHB-mCherry at passage 2 and were imaged at passage 5.

Flow Cytometry.

MCF10A were harvested through trypsinization and resuspended in DMEM/F12 supplemented with 20% horse serum. The cell pellet was washed twice with PBS, and the cells were fixed and permeabilized with ice-cold methanol at −20 . The population was then split, and cells were stained for pHH3 (CST 9706) and either p21 (CST 2947S) or phospho-Rb (S807/811; CST 8516S) followed by secondary antibody staining. Fluorescence intensities for each signal were read on a MoFlo Cytomation and analyzed using custom MATLAB scripts.

Imaging and Image Processing.

Time-lapse imaging, immunofluorescence, image processing, and classification of populations were conducted as previously described (15, 20). Both fixed and time-lapse microscopy images were processed as previously described (26). The tracking code is available for download here: https://github.com/scappell/Cell_tracking. Additional detail is available in .
  51 in total

Review 1.  Control of the cell cycle and apoptosis.

Authors:  A S Lundberg; R A Weinberg
Journal:  Eur J Cancer       Date:  1999-12       Impact factor: 9.162

Review 2.  CDK inhibitors: cell cycle regulators and beyond.

Authors:  Arnaud Besson; Steven F Dowdy; James M Roberts
Journal:  Dev Cell       Date:  2008-02       Impact factor: 12.270

3.  PP2A-mediated regulation of Ras signaling in G2 is essential for stable quiescence and normal G1 length.

Authors:  Nana Naetar; Velmurugan Soundarapandian; Larisa Litovchick; Kelsey L Goguen; Anna A Sablina; Christian Bowman-Colin; Piotr Sicinski; William C Hahn; James A DeCaprio; David M Livingston
Journal:  Mol Cell       Date:  2014-05-22       Impact factor: 17.970

4.  Cell cycle heterogeneity directs the timing of neural stem cell activation from quiescence.

Authors:  L Otsuki; A H Brand
Journal:  Science       Date:  2018-04-06       Impact factor: 47.728

5.  Basal p21 controls population heterogeneity in cycling and quiescent cell cycle states.

Authors:  K Wesley Overton; Sabrina L Spencer; William L Noderer; Tobias Meyer; Clifford L Wang
Journal:  Proc Natl Acad Sci U S A       Date:  2014-09-29       Impact factor: 11.205

6.  The proliferation-quiescence decision is controlled by a bifurcation in CDK2 activity at mitotic exit.

Authors:  Sabrina L Spencer; Steven D Cappell; Feng-Chiao Tsai; K Wesley Overton; Clifford L Wang; Tobias Meyer
Journal:  Cell       Date:  2013-09-26       Impact factor: 41.582

7.  Regulation of retinoblastoma protein functions by ectopic expression of human cyclins.

Authors:  P W Hinds; S Mittnacht; V Dulic; A Arnold; S I Reed; R A Weinberg
Journal:  Cell       Date:  1992-09-18       Impact factor: 41.582

8.  Cell cycle expression and p53 regulation of the cyclin-dependent kinase inhibitor p21.

Authors:  Y Li; C W Jenkins; M A Nichols; Y Xiong
Journal:  Oncogene       Date:  1994-08       Impact factor: 9.867

9.  CDK4/6 and autophagy inhibitors synergistically induce senescence in Rb positive cytoplasmic cyclin E negative cancers.

Authors:  Smruthi Vijayaraghavan; Cansu Karakas; Iman Doostan; Xian Chen; Tuyen Bui; Min Yi; Akshara S Raghavendra; Yang Zhao; Sami I Bashour; Nuhad K Ibrahim; Meghan Karuturi; Jing Wang; Jeffrey D Winkler; Ravi K Amaravadi; Kelly K Hunt; Debu Tripathy; Khandan Keyomarsi
Journal:  Nat Commun       Date:  2017-06-27       Impact factor: 14.919

10.  DNA damage during S-phase mediates the proliferation-quiescence decision in the subsequent G1 via p21 expression.

Authors:  Alexis R Barr; Samuel Cooper; Frank S Heldt; Francesca Butera; Henriette Stoy; Jörg Mansfeld; Béla Novák; Chris Bakal
Journal:  Nat Commun       Date:  2017-03-20       Impact factor: 14.919

View more
  27 in total

Review 1.  Bistable switches as integrators and actuators during cell cycle progression.

Authors:  Wayne Stallaert; Katarzyna M Kedziora; Hui Xiao Chao; Jeremy E Purvis
Journal:  FEBS Lett       Date:  2019-10-16       Impact factor: 4.124

2.  A tumor suppressor Retinoblastoma1 is essential for embryonic development in the sea urchin.

Authors:  Ana Fernandez-Nicolas; Derek Xu; Mamiko Yajima
Journal:  Dev Dyn       Date:  2019-10-01       Impact factor: 3.780

Review 3.  Putting the brakes on the cell cycle: mechanisms of cellular growth arrest.

Authors:  Lindsey R Pack; Leighton H Daigh; Tobias Meyer
Journal:  Curr Opin Cell Biol       Date:  2019-06-25       Impact factor: 8.382

4.  The structure of the human cell cycle.

Authors:  Wayne Stallaert; Katarzyna M Kedziora; Colin D Taylor; Tarek M Zikry; Jolene S Ranek; Holly K Sobon; Sovanny R Taylor; Catherine L Young; Jeanette G Cook; Jeremy E Purvis
Journal:  Cell Syst       Date:  2021-11-19       Impact factor: 10.304

5.  Senescence Evasion in Chemotherapy: A Sweet Spot for p21.

Authors:  Humza M Ashraf; Justin Moser; Sabrina L Spencer
Journal:  Cell       Date:  2019-07-11       Impact factor: 41.582

Review 6.  Slow-cycling (dormant) cancer cells in therapy resistance, cancer relapse and metastasis.

Authors:  Sukanya Basu; Yang Dong; Rahul Kumar; Collene Jeter; Dean G Tang
Journal:  Semin Cancer Biol       Date:  2021-05-09       Impact factor: 15.707

7.  Visualizing the metazoan proliferation-quiescence decision in vivo.

Authors:  Rebecca C Adikes; Abraham Q Kohrman; Michael A Q Martinez; Nicholas J Palmisano; Jayson J Smith; Taylor N Medwig-Kinney; Mingwei Min; Maria D Sallee; Ononnah B Ahmed; Nuri Kim; Simeiyun Liu; Robert D Morabito; Nicholas Weeks; Qinyun Zhao; Wan Zhang; Jessica L Feldman; Michalis Barkoulas; Ariel M Pani; Sabrina L Spencer; Benjamin L Martin; David Q Matus
Journal:  Elife       Date:  2020-12-22       Impact factor: 8.140

Review 8.  Novel insights in cell cycle dysregulation during prostate cancer progression.

Authors:  Salma Ben-Salem; Varadha Balaji Venkadakrishnan; Hannelore V Heemers
Journal:  Endocr Relat Cancer       Date:  2021-05-11       Impact factor: 5.900

Review 9.  Imaging developmental cell cycles.

Authors:  Abraham Q Kohrman; Rebecca P Kim-Yip; Eszter Posfai
Journal:  Biophys J       Date:  2021-05-06       Impact factor: 3.699

10.  p27 allosterically activates cyclin-dependent kinase 4 and antagonizes palbociclib inhibition.

Authors:  Keelan Z Guiley; Jack W Stevenson; Kevin Lou; Krister J Barkovich; Vishnu Kumarasamy; Tilini U Wijeratne; Katharine L Bunch; Sarvind Tripathi; Erik S Knudsen; Agnieszka K Witkiewicz; Kevan M Shokat; Seth M Rubin
Journal:  Science       Date:  2019-12-13       Impact factor: 47.728

View more

北京卡尤迪生物科技股份有限公司 © 2022-2023.