- 1Department of Biology, Stanford University, Stanford, CA, United States
- 2Departments of Pediatrics and Genetics, School of Medicine, Stanford University, Stanford, CA, United States
- 3Chan Zuckerberg Biohub, San Francisco, CA, United States
Every type of cell in an animal maintains a specific size, which likely contributes to its ability to perform its physiological functions. While some cell size control mechanisms are beginning to be elucidated through studies of cultured cells, it is unclear if and how such mechanisms control cell size in an animal. For example, it was recently shown that RB, the retinoblastoma protein, was diluted by cell growth in G1 to promote size-dependence of the G1/S transition. However, it remains unclear to what extent the RB-dilution mechanism controls cell size in an animal. We therefore examined the contribution of RB-dilution to cell size control in the mouse liver. Consistent with the RB-dilution model, genetic perturbations decreasing RB protein concentrations through inducible shRNA expression or through liver-specific Rb1 knockout reduced hepatocyte size, while perturbations increasing RB protein concentrations in an Fah−/− mouse model increased hepatocyte size. Moreover, RB concentration reflects cell size in G1 as it is lower in larger G1 hepatocytes. In contrast, concentrations of the cell cycle activators Cyclin D1 and E2f1 were relatively constant. Lastly, loss of Rb1 weakened cell size control, i.e., reduced the inverse correlation between how much cells grew in G1 and how large they were at birth. Taken together, our results show that an RB-dilution mechanism contributes to cell size control in the mouse liver by linking cell growth to the G1/S transition.
Introduction
Cell size is important for cell physiology because it determines the size of organelles including the nucleus, spindle, centrosome, mitochondrial network, and vacuole (Jorgensen et al., 2007; Neumann and Nurse, 2007; Decker et al., 2011; Rafelski et al., 2012; Good et al., 2013; Hazel et al., 2013; Chan and Marshall, 2014). Through this control of organelle size as well as ribosome number, cell size impacts biosynthesis. In addition to biosynthesis, cell size is important for the function of diverse cell types. Erythrocytes are tiny because they must migrate through tight spaces, while macrophages must be large enough to engulf pathogens. The importance of cell size for diverse cell types is further reflected in the fact that the size of cells of a given type is remarkably uniform and deviations from the typical cell size are often associated with disease states. For example, the coefficient of variation of the size distribution of red blood cells is used as a common diagnostic parameter in comprehensive blood tests, and many cancers, such as small-cell lung cancer, are characterized by altered cell size (Bell and Waizbard, 1986; Evans and Jehle, 1991; Savage et al., 2007; Salvagno et al., 2015; Nguyen et al., 2016; Zatulovskiy and Skotheim, 2020; Sandlin, 2022).
To control their size to be within a target range, proliferating cells can either regulate their growth such that larger cells grow slower, or they can link cell growth to progression through the cell division cycle (Tzur et al., 2009; Kafri et al., 2013; Ginzberg et al., 2015; Ginzberg et al., 2018; Liu et al., 2022). Studies in cultured cells have now begun to elucidate some of the molecular mechanisms that are responsible. One model proposes that cell growth triggers the G1/S transition by diluting the cell cycle inhibitor RB, the retinoblastoma protein (Zatulovskiy et al., 2020). The amount of dilution required, and the target cell size, can then be modulated by the activity of the cyclin D-Cdk4 kinase complex (Tan et al., 2021). A second model proposes that small cell size activates the p38 stress kinase, which, in turn, inhibits cell division (Liu et al., 2018). Taken together, this recent work in the field is starting to uncover a series of molecular mechanisms through which cell growth triggers the G1/S transition to control cell size. However, these cell size control mechanisms are all based on in vitro cell culture studies, and their in vivo relevance remains unclear.
Here, we aimed to test the in vivo relevance of a cell size control model, namely, the RB-dilution model, in the liver of genetically engineered mice. Consistent with the RB-dilution model, increasing or decreasing RB protein concentrations in hepatocytes in vivo caused an increase or decrease in their cell size, respectively. In addition, we found that RB protein concentration decreased with cell size in G1 phase, consistent with RB concentration being a cell size sensor. Finally, we found that RB deletion reduced the ability of primary hepatocytes to control their size when grown in culture. More specifically, we found that deletion of Rb1 reduced the inverse correlation between how large hepatocytes were at birth and how much they grew in G1. Taken together, our results suggest that the RB-dilution mechanism plays an important role in controlling hepatocyte cell size.
Results
Decreasing RB protein concentration in mouse liver reduces hepatocyte size
If the RB-dilution model applies to the mouse liver, reducing RB protein concentration would reduce hepatocyte size. To test this, we utilized two different mouse models to reduce RB protein concentration. The first model harbors a Doxycycline-inducible transgene broadly expressing an shRNA against Rb1 (Rosa-rtTA; TRE-shRb1 mouse (Doan et al., 2021)). We gave Doxycycline containing water (1 g/L) from birth to 4 weeks old to Rosa-rtTA (control) or Rosa-rtTA; TRE-shRb1 mice from the same litter (Figure 1A). Then, the hepatocytes were isolated from these mice to check for knockdown efficiency (Figure 1B). To ensure that there was no difference in cell cycle status, we waited for another 3 weeks to measure their cell size distributions (Figure 1C, Supplementary Figures S1A,D). Since the liver is a highly polyploid organ (Zhang et al., 2019), and polyploid cells are bigger than diploid cells, we sorted the isolated hepatocytes based on viability and on ploidy before size measurements were made in order to eliminate the confounding effect of ploidy on size (Figure 1C, Supplementary Figure S1A). We found that within the same ploidy, the hepatocyte size in Rb1 knockdown mice was smaller than that in the control mice, supporting the RB-dilution model.
FIGURE 1. Decreasing RB protein concentration in mouse liver reduces hepatocyte size. (A). Schematic of experimental timeline performed using Rosa-rtTA; TRE-shRb1 mice. (B). Immunoblot and qPCR analysis following shRb1 induction in mouse liver (n = 4). Mice were harvested after 4 weeks of Dox treatment (1 g/L) that started at birth. (C). Coulter counter cell size measurement of the hepatocytes isolated from Rosa-rtTA (Ctrl, n = 4) and Rosa-rtTA; TRE-shRb1 (shRb1, n = 4) mice at 3 weeks after Dox removal. (D) Schematic of experimental timeline performed using Rb1 flx/flx; Alb-CreERT2 mice. (E) Immunoblot shows RB protein concentrations in the livers of Rb1 +/+, Rb1 +/-, and Rb1 −/− mice. Mice were harvested at 6 weeks after Tamoxifen treatment. The right panel shows the quantification of protein concentration. RB protein concentration was normalized using actin. (F) Cell size measurements were made using a Coulter counter to examine hepatocytes isolated from Rb1 +/+ (n = 6 for tetraploid, n = 4 for octoploid), Rb1 +/− (n = 8 for tetraploid, n = 4 for octoploid), and Rb1 −/− (n = 6 for tetraploid, n = 4 for octoploid) mice at 6 weeks after Tamoxifen treatment. Each dot represents data from one mouse. The error bars indicate standard deviation. **p < 0.01.
To further test the RB-dilution model in the liver, we used a second mouse model that also reduced the RB protein concentration, but only in the liver. This mouse model contains an Rb1flx/flx allele and expresses the Cre recombinase only in the liver through an Alb-CreERT2 allele. By injecting Tamoxifen into these mice, the Rb1 gene can be knocked out in a liver-specific manner. We generated Rb1−/−, Rb1−/+, and Rb1+/+ mice from the same litter by administering Tamoxifen to these mice when they were 15 days old (Figure 1D). At 6 weeks after Tamoxifen treatment, we harvested the livers and examined the hepatocyte size of these mice. We note that the Rb1+/− mice did not have significantly reduced RB protein concentrations compared to the Rb1+/+ mice (Figure 1E). Consistent with this similar RB protein expression, we found that hepatocyte size in Rb1+/− mice was not significantly different from that in Rb1+/+ mice (Figure 1F, Supplementary Figure S1B). In contrast, hepatocyte size was significantly smaller in Rb1−/− mice compared to Rb1+/+ and Rb1+/− mice, which again supports a role for RB in controlling cell size in the mouse liver (Figure 1F, Supplementary Figure S1B).
Increasing RB concentration in mouse liver increases hepatocyte cell size
After finding that decreasing RB concentration decreases cell size, we sought to test if increasing RB protein concentration increases hepatocyte cell size. To overexpress RB, we utilized the Fah−/− mouse system (Azuma et al., 2007). Deletion of the Fah gene causes toxin accumulation in hepatocytes that will eventually lead to hepatocyte death. Toxin accumulation is prevented in Fah−/− mice by treating them with NTBC (2-(2-nitro-4-trifluoromethylbenzoyl)-1,3-cyclohexanedione) (Grompe, 2017). When NTBC is withdrawn, cells expressing exogenous Fah, introduced by injecting Fah + transposons, will clonally expand to repopulate the injured liver (Figure 2A) (Wuestefeld et al., 2013). Importantly, other genetic elements, such as transgenes, can be added to the Fah transposon so that they are co-integrated into some hepatocyte genomes. This system provides a versatile platform for genetic manipulations in the liver. To increase RB expression, we delivered an Fah-P2A-Rb1 transposon and the SB100 transposase into Fah−/− mice via hydrodynamic transfection through the tail vein. As a control, we injected the transposon with Fah-P2A-mCherry (Figure 2B). 6 weeks after injection, the liver was almost fully repopulated by the cells containing Fah transposons, and there was an appreciable level of RB protein overexpression (Figures 2C,D, Supplementary Figure S2C). Consistent with the RB-dilution model, overexpression of RB resulted in bigger hepatocytes (Figure 2E, Supplementary Figure S1C).
FIGURE 2. Increasing Rb1 concentration in mouse liver increases hepatocyte cell size. (A) Schematic of Fah−/− system as a platform for genetic manipulation of hepatocytes in the liver. (B) Schematic of experimental timeline. Fah−/− mice were injected with Fah transposons as indicated, and their livers were harvested 4 weeks after the injection. (C,D) Immunofluorescence staining (C) and immunoblots (D) were performed to measure Fah and RB expression. (E) Cell size measurement of the hepatocytes isolated from Fah−/− mice injected with Fah transposons were made with a Coulter counter. Each dot represents data from one mouse (n = 8 for tetraploid, n = 7 for octoploid). The error bars indicate standard deviation. *p < 0.05.
In addition to RB dosage affecting cell size in cultured cells, it was also reported that the activity of the upstream kinase cyclin D-Cdk4/6 affected cell size (Tan et al., 2021). Lower Cdk4/6 activity resulted in larger cells, while higher activity reduced cell size. This can be integrated with the RB-dilution model since cyclin D-Cdk4/6-dependent RB phosphorylation contributes to RB inactivation. This means that the higher the Cdk4/6 activity, the more RB is phosphorylated and inactivated, and the less it needs to be diluted (Rubin et al., 2020). To test for the effect of Cdk4/6 activity on hepatocyte cell size, we used Palbociclib to inhibit Cdk4/6 activity during liver regeneration following a partial hepatectomy (Supplementary Figure S2A). In this experiment, the partial hepatectomy induces massive hepatocyte growth and proliferation during a short time window that can be monitored experimentally. More specifically, following partial hepatectomy, mice were treated with Palbociclib (80 mg/kg) for five consecutive days (Supplementary Figure S2A). We isolated hepatocytes from these mice at different time points and measured cell size. Consistent with the observations in cultured cells that the Cdk4/6-RB pathway controls cell size, Palbociclib treatment increased hepatocyte size following a partial hepatectomy (Supplementary Figure S2B).
RB protein is diluted in G1 in primary hepatocytes
So far, our data indicate that manipulations of RB protein concentrations in vivo changed hepatocyte size as the RB-dilution model predicted. However, these deletion and over-expression experiments cannot determine if RB protein concentration functions as a cell size sensor. For RB to be a cell size sensor responsible for cell size control, not only must its concentration affect cell size, but its concentration must also change with cell size. In other words, RB must be diluted in G1 for the model to work.
To test if RB is diluted in G1, we sought to measure the amounts of RB protein and cell size in single cells. However, detecting RB protein in single proliferating cells from the mouse liver tissue is challenging for several reasons. First, immunofluorescence measurements are challenging due to RB’s low expression level and the high autofluorescence background of hepatocytes. Second, RB-dilution would control cell size by coordinating cell growth with the G1/S transition only in proliferating cells, who represent only a minority of adult hepatocytes (Magami et al., 2002; Wei et al., 2021). To circumvent these challenges, we sought to specifically measure RB protein concentrations in proliferating isolated primary hepatocytes. To do this, we used a 2D primary hepatocyte culture protocol (Peng et al., 2018; Jin et al., 2021) (Figure 3A). In these cultures, primary hepatocytes actively enter the cell cycle, and we can detect RB protein using immunfluorescence because the autofluorescence is significantly reduced in the isolated hepatocytes compared to those in the mouse liver.
FIGURE 3. RB protein is diluted in G1 in primary hepatocytes. (A) Schematic of experiment examining RB protein concentrations in single primary hepatocytes in 2D cultures. (B) Representative immunofluorescence images of primary hepatocytes isolated from Fucci2 mice and illustration of gating procedure used to isolate G1 cells. (C) Immunofluorescence staining of RB in Fucci2 primary hepatocytes. The Rb1−/− hepatocytes were used as a negative control for RB antibody staining. (D) Immunofluorescence staining of Rpb1, Hmgb1, CycD1, and E2f1 in Fucci2 primary hepatocytes grown in 2D culture. (E) Nuclear area of G1 hepatocytes plotted against the indicated protein amount (left panels) or concentration (right panels). The shaded area indicates 95% confidence interval.
To measure the RB protein concentrations in hepatocytes, we performed an immunofluorescence staining on 2D hepatocyte cultures isolated from Fucci2 mice expressing G1 and S/G2 cell cycle phase reporters (Abe et al., 2013) (Figure 3A,B). Based on their DNA content, Cdt1-mCherry signal, and Geminin-Venus signal, we determined which cells in the population were in G1 phase (Figure 3B). We then examined the relationship between target protein amount and cell size, which was estimated using nuclear area as a proxy for cell size (Figure 3C,D) (Berenson et al., 2019). We also measured the RNA Polymerase II subunit Rpb1, whose concentration was expected to stay constant, and Hmgb1, which has been found to be diluted by cell growth similarly to RB in other cell lines (Figure 3D) (Lanz et al., 2021). We found that unlike Rpb1, the total amount of RB protein did not increase in proportion to cell size, but behaved more similarly to the anticipated subscaling protein Hmgb1 (Figure 3E). We also examined amounts of two key activators of the G1/S transition, Cyclin D1 and E2f1 (Figure 3D). While RB protein amounts remained relatively unchanged as cells became larger in G1, Cyclin D1 and E2f1 protein amounts increased with cell size (Figure 3E), suggesting that the activator (Cyclin D1, E2f1) to inhibitor (RB) ratio is increasing as cells grow in G1. While immunofluorescence measurements to quantify protein amounts are not necessarily linear due to the particulars of specific antibody binding, our measurements of RB concentration as a function of nuclear area clearly demonstrate that RB is diluted as cells grow in G1. Taken together, our immunofluorescence measurements of RB concentration in proliferating hepatocytes are consistent with RB operating as a cell size sensor in G1.
Knocking out Rb1 weakens G1/S cell size control in primary hepatocytes
At its most fundamental level, cell size control requires cells that are born smaller to grow more during the cell division cycle than cells that are born larger. Size control thus results in more similarly sized cells following a division cycle. The degree to which cells control their size can therefore be quantified by the inverse correlation between cell size at birth and the amount of cell growth. Therefore, the strength of cell size control can be evaluated by the slope of the linear fit between the cell size at birth and the amount of cell growth over the cell cycle (Cadart et al., 2018; Zatulovskiy and Skotheim, 2020). Perfect size control would yield a slope value of -1, and if the amount of growth was independent of the birth size, the slope would be 0. If RB contributed to cell size control, its deletion should reduce the degree of inverse correlation between cell size at birth and the amount of cell growth in the subsequent cell cycle (Zatulovskiy et al., 2020).
To test the extent to which RB was responsible for cell size control in G1, we sought to measure cell size and cell cycle phase durations in proliferating hepatocytes that we tracked through cell divisions. To do this, we performed live cell imaging on the 2D cultured hepatocytes isolated from Rb1flx/flx; Rosa26-CreERT2; Fucci2 mice. Four days before imaging, cells were treated with 4-Hydroxytamoxifen (4-OHT) to delete the Rb1 gene (Figure 4A). We used the inflection point of the Geminin-Venus S/G2 reporter signal and the maximum point of the Cdt1-mCherry G1 reporter signal to determine the G1/S transition time at the midpoint between these two events (Figure 4B). We again used the nuclear area as a proxy for cell size. Consistent with RB contributing to cell size control, we found that the Rb1−/− hepatocytes exhibited smaller size and shorter G1 length (Figure 4C). Moreover, the correlation between birth nuclear area and the amount of growth in G1 was reduced in Rb1−/− hepatocytes compared to Rb1+/+ hepatocytes (Figures 4D,E, Supplementary Figures S2D,E). We note that in each experiment, we used hepatocytes from different mice, but that experimental and control hepatocytes for each replicate were derived from the same mouse. Taken together, our data show that knocking out Rb1 weakened G1/S cell size control in hepatocytes, consistent with the RB-dilution model.
FIGURE 4. Deletion of Rb1 weakens G1/S cell size control in primary hepatocytes. (A) Schematic of the experiment used to examine the role of Rb1 in hepatocyte cell size control. (B) Live cell imaging traces of the Cdt1-mCherry and Geminin-Venus cell cycle reporters in Fucci2 hepatocytes grown in 2D culture. The G1/S transition is determined using the inflection point of the Geminin-Venus signal and the point of maximum Cdt1-mCherry signal. (C) Cell size and cell cycle statistics of Rb1+/+ and Rb1−/− hepatocytes. Cell size is approximated using nuclear area. (D) The correlations between nuclear area at birth and the nuclear area changes during G1 in Rb1+/+ and Rb1−/− hepatocytes. The shaded area indicates the 95% confidence interval. (E) The slope values of the linear relationship in (D) are plotted from three independent biological replicates. Each experiment used hepatocytes from different mice, and Rb1+/+ and Rb1−/− hepatocytes for each replicate experiment were derived from the same mouse. Dashed lines link experiments performed using hepatocytes derived from the same mouse. (F) p107 mRNA measured by RT-qPCR from Rb1+/+ and Rb1−/− 2D cultured primary hepatocytes (left panel) and Rosa-rtTA; TRE-shRb1 mice (right panel). The error bars indicate standard deviation. *p < 0.05, **p < 0.01. (G) Model indicating multiple mechanisms through which RB regulates cell size at the G1/S transition.
While Rb1 deletion weakened G1/S size control, it did not remove it completely as the slope between nuclear area at birth and G1 growth amount was not 0 (Figure 4D). It was previously shown that Rb1 deletion resulted in the upregulation of the related protein p107, which can compensate for Rb1 loss (Sage et al., 2000). Therefore, the partial loss of size control in Rb1−/− cells might be explained by p107 upregulation because p107 protein was previously found to be similarly diluted as RB (Zatulovskiy et al., 2020). Moreover, the combined deletion of Rb1 and p107 genes in mouse embryonic fibroblasts resulted in much smaller cells than either of the single deletions (Dannenberg et al., 2000; Sage et al., 2000). To test if p107 was upregulated in Rb1−/− hepatocytes, we measured the mRNA levels of p107 in Rb1−/− 2D cultured hepatocytes, as well as in the Rosa-rtTA; TRE-shRb1 mouse livers. Indeed, p107 mRNA was increased in Rb1−/− and Rb1 knock down cells (Figure 4F), while another member of the pocket protein family, p130, remained same, supporting our hypothesis that p107 is upregulated to compensate for Rb1 loss and contributes to G1/S size control when Rb1 was deleted (Figure 4G).
Discussion
Cell size is fundamental to cell physiology, yet we are only now beginning to understand the molecular mechanisms that control size through linking cell growth to division. Cell growth was found to trigger cell division by diluting a cell cycle inhibitor in budding yeast and Arabidopsis plant cells (Lanz et al., 2021; Schmoller et al., 2015; Xie et al., 2022; Swaffer et al., 2021; D’Ario et al., 2021). This general principle of inhibitor-dilution was also found in human cells grown in culture, who diluted the cell cycle inhibitor RB in G1 to trigger the G1/S transition (Zatulovskiy et al., 2020). This RB-dilution mechanism can be modulated by cyclin D-Cdk4/6 activity, which inactivates RB (Tan et al., 2021). However, it was not clear whether the RB-dilution mechanism was physiologically relevant to cells growing in tissues in vivo. This is because of the large differences between the conditions cells experience when growing in a 2D in vitro culture compared to a 3D tissue environment. In tissues, cells send and receive many more paracrine and endocrine signals, have multiple cell-cell contacts, experience dynamic forces and pressures, and experience an environment with very different nutrients and growth factors. In addition to these differences in culture and tissue environments, size control models developed in culture are cell autonomous and do not take into account the possibility that tissue level processes could be used to sense and control the size of neighboring cells.
Due to the significant differences of in vitro and in vivo growth conditions, it is crucial to test molecular models developed using cell culture in animal contexts. We therefore sought to test the RB-dilution model in vivo in the mouse liver. The RB-dilution model has two main requirements. First, that changing the concentration of RB should have an effect on the average cell size, and second, that the concentration of endogenous RB should reflect cell size. Both of these requirements were satisfied in mouse hepatocytes. By using genetically engineered mouse models to manipulate RB concentrations, we found hepatocyte size was correlated with RB protein concentrations in mouse liver. Moreover, the endogenous RB protein was diluted in larger G1 cells. Finally, we examined cell size control in proliferating isolated primary hepatocytes transferred into culture. Loss of Rb1 weakened G1/S cell size control, i.e., reduced the correlation between how much cells grew in G1 and how large they were at birth. Taken together, these results support a role for the RB-dilution model in regulating hepatocyte cell size. As previous work showed similar mechanisms operate in yeast and plants, it is now clear that examples of the inhibitor-dilution class of mechanisms can be found across the tree of life.
In addition to being broadly applicable, inhibitor dilution models can easily be modified to link cell size to cell cycle progression in polyploid cells, which are common in mouse liver. Indeed, polyploidy is a common feature of mammalian livers, where ploidy is generally proportional to cell size (Mu et al., 2020). The RB-dilution model is compatible with hepatocyte polyploidization, which is predominantly mediated by cytokinesis failure (Celton-Morizur and Desdouets, 2010; Gentric et al., 2012). This is because the link between cell growth and cell cycle progression through the dilution of a cell cycle inhibitor does not actually require division. RB may be diluted in G1 so that S phase is triggered. The, RB accumulates through S/G2 and cells can return to a subsequent stable tetraploid G1 phase without division. Crucially, the RB-dilution model is based on protein concentrations, which will be the same in two diploid G1 cells (if division was successful) as in one tetraploid G1 cell (if division was unsuccessful). Inhibitor dilution thus potentially provides an elegant mechanism to explain the repeated growth requirements for DNA synthesis and the large size of polyploid eukaryotic cells.
While our data here support a role for RB-dilution in controlling cell size in the mouse liver, they also indicate that this is unlikely to be the only size control mechanism working in this context. The presence of compensatory mechanisms can be inferred from the fact that Rb1 deletion and over-expression produce modest changes in cell size (Figures 1, 2). Moreover, Rb1 deletion in primary hepatocytes grown in culture only partially eliminated the correlation between growth and progression through G1 (Figure 4). One possible explanation for these results is that other cell size control mechanisms inhibit the division of small cells in the absence of RB. One possibility is that signaling through the stress activated kinase p38 inhibits cell division in small cells to prevent their division and allow them more time to grow in G121,22. A second possibility is that RB loss is compensated for by the upregulation of another cell cycle inhibitor. Indeed, RB regulates the expression of the related cell cycle inhibitor p107 (Dannenberg et al., 2000; Sage et al., 2000; Burkhart et al., 2010; Wirt and Sage, 2010). Consistent with prior work, we found that deletion of Rb1 In hepatocytes strongly increased p107 expression. This regulation of p107 by RB could buffer cell size against changes in RB protein concentrations (Figure 4G). Further investigation into this mechanism is needed, which can be achieved by incorporating p107 knockout in the Rb1 knockout background.
In addition to p107 regulation buffering the effect of Rb1 deletion, there could be additional size-control mechanisms. In yeast, Whi5-dilution is one of several mechanisms that operate based on size-dependent gene expression (Keifenheim et al., 2017; Chen et al., 2020). In this model, a series of cell cycle activators generally increase in concentration as cells grow, while cell cycle inhibitors decrease in concentration. The balance of activators and inhibitors then sets cell size. Removal of any one component would serve to shift the mean cell size, but not remove the ability of the cell to control size around the new setpoint. Supporting the applicability of this multiple size-scaling mechanism model to mammalian cells, recent proteomics analyses found that cell size-dependent concentration changes were widespread across the proteome of human cells grown in culture (Cheng et al., 2021; Lanz et al., 2021). Moreover, beyond such size-scaling concentration-based mechanisms, there may be additional mechanisms operating by other mechanistic principles yet to be discovered. Much more work remains to be done to identify the mechanistic links between cell growth and division. As a first step, we anticipate such work will need to be done in yeast and cultured cells. While this work will need to be validated in vivo, our results here, indicating that the RB-dilution model contributes to controlling cell size in hepatocytes, should serve as encouragement for future studies in vitro and their application in vivo.
Materials and methods
Mice
All mice were handled in accordance with the guidelines of the Institutional Animal Care and Use Committee at Stanford University. The Rosa-rtTA; TRE-shRb1 mice were described before (Doan et al., 2021). Rosa-rtTA or Rosa-rtTA; TRE-shRb1 mice were given Doxycyclin water (RPI, 1 g/L) from birth to when they were 4 weeks old. At 6 weeks old, the mice were sacrificed to isolate primary hepatocytes. Alb-CreERT2 mice were a kind gift of Dr. Pierre Chambon 51, and Rb1flx/flx; p107−/−; p130flx/flx; Alb-CreERT2 were described before 52. We bred Rb1flx/flx; p107−/−; p130flx/flx; Alb-CreERT2 mice to Rb1flx/flx mice to generate Rb1flx/flx; Alb-CreERT2, Rb1flx/+; Alb-CreERT2, and Rb1+/+; Alb-CreERT2 mice from the same cohort. Then, tamoxifen (Sigma, 50 mg/kg, five consecutive doses) was given to these mice when they were ∼15 days old through i. p. Injection. At 6 weeks old, the mice were sacrificed to isolate primary hepatocytes. The Fah−/− mice were kindly shared by Dr. Markus Grompe’s group at Oregon Health and Science University. At 6 weeks old, the Fah−/− mice were subjected to hydrodynamical transfection through their tail vein. The DNA plasmids we injected contain a transposon and transposase. The plasmid backbone was kindly provided by Dr. Hao Zhu’s lab at UT Southwestern Medical Center. We modified the transposons to express the Fah gene and mCherry or Rb1. The transposase was SB100.6 weeks after the injection, the mice were sacrificed to isolate primary hepatocytes. The Fucci2 mice (Abe et al., 2013) were kindly shared by Dr. Valentina Greco’s lab at Yale University. The Fucci2 mice contained the transgene expressing Cdt1-mCherry and Geminin-Venus from the ROSA26 locus.
Partial hepatectomy and palbociclib treatment
2/3rd partial hepatectomies were performed on 8-week-old male CD1 mice (JAX). The procedure followed an established protocol 53. At day 0 to day 5 after the surgery, mice were injected with either vehicle or Palbociclib HCl (80 mg/kg). Palbociclib HCl was dissolved in 50 mmol/L sodium lactate (pH 4) at a concentration of 10 μg/ml under a fume hood. Then, the drug was given to mice at a dosage of 80 mg/kg via oral gavage for five consecutive days. On day 0, 1, 2, 3, 4, 5 after surgery, one mouse from each treatment was sacrificed to isolate primary hepatocytes.
Primary hepatocyte isolation and size measurement
Primary hepatocytes were isolated by two-step collagenase perfusion 54. Liver perfusion medium (Thermo Fisher Scientific, 17701038), liver digest medium (Thermo Fisher Scientific, 17703034) and hepatocyte wash medium (Thermo Fisher Scientific, 17704024) were used. After isolation, cells were washed with PBS once and stained with Zombie Green (Biolegend, 423,111) at room temperature for 15 min to mark dead cells. After washing with cold PBS twice, cells were fixed in 4% PFA (Thermo Fisher Scientific, 28,906) at 37°C for 15 min. Then, cells were washed with PBS 3 times, subjected to DNA staining using Hoechst (Thermo Fisher Scientific, 62,249), and sorted on a Flow Cytometer (BD FACSAria II) based on their live/dead staining and their DNA staining. Only live (without Zombie Green signal) hepatocytes with same ploidy were sorted. 10,000–20000 cells were sorted for each population. After sorting, the cell size distributions for each population were measured using a Z2 Coulter counter (Beckman).
Since male and female mice have different hepatocyte size, they were analyzed separately, normalized, and then combined. Therefore, only when there were all genotypes within one gender, the mice of that gender in that litter would be used for the experiment. The mice genotypes in each litter used for this study can be found in Supplementary Figure S1.
Primary hepatocyte 2D culture
The protocol for 2D culture of primary hepatocytes was kindly shared by Yinhua Jin from Dr. Roeland Nusse’s lab at Stanford University (Peng et al., 2018; Jin et al., 2021). Briefly, primary hepatocytes from Rb1flx/flx; Rosa26-CreER2; Fucci2 mice were isolated by two-step collagenase perfusion. After isolation, cells were washed 3 times with hepatocyte wash medium (Thermo Fisher Scientific, 17704024). Cells were then plated in a 6-well plate precoated with collagen I (50 μg/ml) at a density of 200,000 per well. The culture medium contained 3 µM CHIR99021 (Peprotech), 25 ng/ml EGF (Peprotech), 50 ng/ml HGF (Peprotech), and 100 ng/ml TNFa (Peprotech) in Basal medium. The Basal medium contained William’s E medium (GIBCO), 1% Glutamax (GIBCO), 1% Non-Essential Amino Acids (GIBCO), 1% Penicillin/streptomycin (GIBCO), 0.2% normocin (Invitrogen), 2% B27 (GIBCO), 1% N2 supplement (GIBCO), 2% FBS (Corning), 10 mM nicotinamide (Sigma), 1.25 mM N-acetylcysteine (Sigma), 10 µM Y27632 (Peprotech) and 1 µM A83-01 (Tocris). The culture medium was refreshed every other day. Cells were passaged via trypsinization using TrypLE (Thermo-Fisher).
Live cell imaging and analysis
The cultured Fucci2 hepatocytes were seeded on collagen coated 35-mm glass-bottom dishes (MatTek) 1 day before imaging. Then, the cells were transferred to a Zeiss Axio Observer Z1 microscope equipped with an incubation chamber and imaged for 48 h at 37°C and 5% CO2. Brightfield and fluorescence images were collected from Rb1+/+ and Rb1−/− Fucci2 hepatocytes at multiple positions every 20 min using an automated stage controlled by the Micro-Manager software. We used a Zyla 5.5 sCMOS camera and an A-plan 10x/0.25NA Ph1 objective. For each cell, the transition from G1 to S phase was taken as the midpoint between when the inflection point of the Geminin-Venus signal and the point of maximum Cdt1-mCherry signal occurred. Cell size was approximated using the area of the nucleus (Berenson et al., 2019).
Immunoblots
Mouse liver tissues or isolated primary hepatocyte pellets were homogenized and lysed in RIPA buffer (Thermo Scientific) supplemented with protease inhibitors and phosphatase inhibitors. Lysates were separated on 8% SDS-PAGE gels and transferred to nitrocellulose membranes. Membranes were then blocked with SuperBlock™ (TBS) Blocking Buffer (Thermo Fisher Scientific) and incubated overnight at 4°C with primary antibodies in 3% BSA solution in PBS. The primary antibodies were detected using the fluorescently labeled secondary antibodies IRDye® 680LT Goat anti-Mouse IgG (LI-COR 926–68020) and IRDye® 800CW Goat anti-Rabbit IgG (LI-COR 926–32211). Membranes were imaged on a LI-COR Odyssey CLx and analyzed with LI-COR Image Studio software. The following antibodies were used: anti-β-Actin (Sigma, A2103, 1:2000), Anti-Rb (Santa cruz, sc-74570, 1:500).
Immunofluorescence staining
Fucci2 hepatocytes were seeded on a 35-mm collagen-coated glass-bottom dish (MatTek) 1 day before immunofluorescence staining and analysis. For the staining, cells were fixed with 4% paraformaldehyde for 20 min at room temperature, permeabilized with 0.5% Triton™ X-100 (Sigma-Aldrich) for 15 min at room temperature, and then blocked with 3% BSA in PBS. Then, the cells were incubated with primary antibodies overnight at 4°C. After three washes with PBS, the cells were incubated with Alexa Fluor 647 conjugated goat anti-mouse secondary antibodies (Invitrogen, A32728) at 1:1,000 for 1 h at room temperature. After three washes with PBS, cells were incubated with 500 nM DAPI for 30 min at room temperature before imaging. The primary antibodies used for immunofluorescence were Anti-Rb (Santa cruz, sc-74570, 1:100), Anti-Rpb1-CTD (Abcam, ab252854, 1:250), Anti-HMGB1 (Abcam, ab79823, 1:250), Anti-E2F1 (Santa cruz, sc-251, 1:100), and Anti-Ccnd1 (Abcam, ab16663, 1:250). The cells were imaged using a Zeiss Axio Observer Z1 microscope with an A-plan 10x/0.25NA objective.
Image analysis
For microscopy data quantification, cell nuclei were segmented using either the mCherry channel for cells in G1 that express Cdt1-mCherry or the GFP channel for cells in S/G2 that express Geminin-Venus. Segmentation was performed using the Fiji plugin StarDist, which is a deep-learning tool for segmenting nuclei in images that are difficult to segment using thresholding-based methods. The total pixel intensities within the segmented masks in each channel were recorded, and each object’s background was subtracted based on the median intensity of the image. G1 cells were selected based on their Cdt1-mCherry and Geminin-Venus signals. The cells that had high Cdt1-mCherry and low Geminin-Venus were selected as G1 cells. When analyzing the scaling behavior of immunostained proteins, nuclear area was used as a proxy for cell size 7.
RNA extraction and qRT-PCR
Total RNA was isolated using Direct-zol RNA Miniprep kit (Zymo Research). For qRT-PCR, cDNA synthesis was performed with 1 μg of total RNA using an iScript Reverse Transcription Kit (Biorad). qPCR reactions were made with the 2x SYBR Green Master Mix (Biorad). Gene expression levels were measured using the ΔΔCt method.
Statistical analysis
The data in most figure panels reflect multiple experiments performed on different days using mice derived from different litters. For comparison between groups, we conducted the hypothesis test using the two-tailed student’s t-test. Statistical significance is displayed as p < 0.05 (∗) or p < 0.01 (∗∗) unless specified otherwise.
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Ethics statement
The animal study was reviewed and approved by the Administrative panel on laboratory animal care of Stanford University.
Author contributions
SZ and JMS designed the project, performed the experiments, and wrote the manuscript. JS and JA shared the genetically engineered mouse models, and assisted with manuscript preparation. EZ assisted with project design and manuscript preparation.
Funding
This work was supported by the National Institutes of Health (R01 DK128578, P01 CA254867).
Acknowledgments
We thank Roel Nusse and Jinghua Li for assistance in establishing the protocol to culture primary hepatocytes in 2D, Thuyen Nguyen for maintaining the mouse colonies in the Sage lab, Dr. Marcus Grompe for sharing the Fah−/− mice, and the members of the Skotheim laboratory for discussions and feedback on the work.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2022.965595/full#supplementary-material
References
Abe, T., Sakaue-Sawano, A., Kiyonari, H., Shioi, G., Inoue, K. i., Horiuchi, T., et al. (2013). Visualization of cell cycle in mouse embryos with Fucci2 reporter directed by Rosa26 promoter. Development 140, 237–246. doi:10.1242/dev.084111
Azuma, H., Paulk, N., Ranade, A., Dorrell, C., Al-Dhalimy, M., Ellis, E., et al. (2007). Robust expansion of human hepatocytes in Fah-/-/Rag2-/-/Il2rg-/- mice.-/-/Il2rg-/- mice. Nat. Biotechnol. 25, 903–910. doi:10.1038/nbt1326
Bell, C. D., and Waizbard, E. (1986). Variability of cell size in primary and metastatic human breast carcinoma. Invasion Metastasis 6, 11–20.
Berenson, D. F., Zatulovskiy, E., Xie, S., and Skotheim, J. M. (2019). Constitutive expression of a fluorescent protein reports the size of live human cells. Mol. Biol. Cell 30, 2985–2995. doi:10.1091/mbc.E19-03-0171
Burkhart, D. L., Wirt, S. E., Zmoos, A.-F., Kareta, M. S., and Sage, J. (2010). Tandem E2F binding sites in the promoter of the p107 cell cycle regulator control p107 expression and its cellular functions. PLoS Genet. 6, e1001003. doi:10.1371/journal.pgen.1001003
Cadart, C., Monnier, S., Grilli, J., Saez, P. J., Srivastava, N., Attia, R., et al. (2018). Size control in mammalian cells involves modulation of both growth rate and cell cycle duration. Nat. Commun. 9, 3275. doi:10.1038/s41467-018-05393-0
Celton-Morizur, S., and Desdouets, C. (2010). Polyploidization of liver cells. Adv. Exp. Med. Biol. 676, 123–135. doi:10.1007/978-1-4419-6199-0_8
Chan, Y.-H. M., and Marshall, W. F. (2014). Organelle size scaling of the budding yeast vacuole is tuned by membrane trafficking rates. Biophys. J. 106, 1986–1996. doi:10.1016/j.bpj.2014.03.014
Chen, Y., Zhao, G., Zahumensky, J., Honey, S., and Futcher, B. (2020). Differential scaling of gene expression with cell size may explain size control in budding yeast. Mol. Cell 78, 359.e6–370.e6. doi:10.1016/j.molcel.2020.03.012
Cheng, L., Chen, J., Kong, Y., Tan, C., Kafri, R., and Björklund, M. (2021). Size-scaling promotes senescence-like changes in proteome and organelle content. 2021.08.05.455193. doi:10.1101/2021.08.05.455193
Dannenberg, J. H., van Rossum, A., Schuijff, L., and te Riele, H. (2000). Ablation of the retinoblastoma gene family deregulates G(1) control causing immortalization and increased cell turnover under growth-restricting conditions. Genes Dev. 14, 3051–3064. doi:10.1101/gad.847700
D’Ario, M., Tavares, R., Schiessl, K., Desvoyes, B., Gutierrez, C., Howard, M., et al. (2021). Cell size controlled in plants using DNA content as an internal scale. Science 372, 1176–1181. doi:10.1126/science.abb4348
Decker, M., Jaensch, S., Pozniakovsky, A., Zinke, A., O'Connell, K. F., Zachariae, W., et al. (2011). Limiting amounts of centrosome material set centrosome size in C. elegans embryos. Curr. Biol. 21, 1259–1267. doi:10.1016/j.cub.2011.06.002
Doan, A., Arand, J., Gong, D., Drainas, A. P., Shue, Y. T., Lee, M. C., et al. (2021). RB depletion is required for the continuous growth of tumors initiated by loss of RB. PLoS Genet. 17, e1009941. doi:10.1371/journal.pgen.1009941
Evans, T. C., and Jehle, D. (1991). The red blood cell distribution width. J. Emerg. Med. 9 (1), 71–74. doi:10.1016/0736-4679(91)90592-4
Gentric, G., Desdouets, C., and Celton-Morizur, S. (2012). Hepatocytes polyploidization and cell cycle control in liver physiopathology. Int. J. Hepatol. 2012, 282430. doi:10.1155/2012/282430
Ginzberg, M. B., Chang, N., D'Souza, H., Patel, N., Kafri, R., and Kirschner, M. W. (2018). Cell size sensing in animal cells coordinates anabolic growth rates and cell cycle progression to maintain cell size uniformity. eLife 7, e26957. doi:10.7554/eLife.26957
Ginzberg, M. B., Kafri, R., and Kirschner, M. (2015). Cell biology. On being the right (cell) size. Science 348, 1245075. doi:10.1126/science.1245075
Good, M. C., Vahey, M. D., Skandarajah, A., Fletcher, D. A., and Heald, R. (2013). Cytoplasmic volume modulates spindle size during embryogenesis. Science 342, 856–860. doi:10.1126/science.1243147
Grompe, M. (2017). Fah knockout animals as models for therapeutic liver repopulation. Adv. Exp. Med. Biol. 959, 215–230. doi:10.1007/978-3-319-55780-9_20
Hazel, J., Krutkramelis, K., Mooney, P., Tomschik, M., Gerow, K., Oakey, J., et al. (2013). Changes in cytoplasmic volume are sufficient to drive spindle scaling. Science 342, 853–856. doi:10.1126/science.1243110
Jin, Y., Anbarchian, T., Wu, P., Sarkar, A., Fish, M., and Nusse, R. (2021). Hepatocyte cell cycle progression depends on a transcriptional repressor cascade downstream of wnt signaling. bioRxiv
Jorgensen, P., Edgington, N. P., Schneider, B. L., Rupes, I., Tyers, M., and Futcher, B. (2007). The size of the nucleus increases as yeast cells grow. Mol. Biol. Cell 18, 3523–3532. doi:10.1091/mbc.e06-10-0973
Kafri, R., Levy, J., Ginzberg, M. B., Oh, S., Lahav, G., and Kirschner, M. W. (2013). Dynamics extracted from fixed cells reveal feedback linking cell growth to cell cycle. Nature 494, 480–483. doi:10.1038/nature11897
Keifenheim, D., Sun, X. M., D'Souza, E., Ohira, M. J., Magner, M., Mayhew, M. B., et al. (2017). Size-dependent expression of the mitotic activator Cdc25 suggests a mechanism of size control in fission yeast. Curr. Biol. 27, 1491e4–1497.e4. doi:10.1016/j.cub.2017.04.016
Lanz, M. C., Zatulovskiy, E., Swaffer, M. P., Zhang, L., Ilerten, I., Zhang, S., et al. (2021). Increasing cell size remodels the proteome and promotes senescence. 2021.07.29.454227. doi:10.1101/2021.07.29.454227
Liu, S., Ginzberg, M. B., Patel, N., Hild, M., Leung, B., Li, Z., et al. (2018). Size uniformity of animal cells is actively maintained by a p38 MAPK-dependent regulation of G1-length. eLife 7, e26947. doi:10.7554/eLife.26947
Liu, X., Yan, J., and Kirschner, M. W. (2022). Beyond G1/S regulation: How cell size homeostasis is tightly controlled throughout the cell cycle? 2022.02.03.478996. doi:10.1101/2022.02.03.478996
Magami, Y., Azuma, T., Inokuchi, H., Kokuno, S., Moriyasu, F., Kawai, K., et al. (2002). Cell proliferation and renewal of normal hepatocytes and bile duct cells in adult mouse liver. Liver 22, 419–425. doi:10.1034/j.1600-0676.2002.01702.x
Mu, L., Kang, J. H., Olcum, S., Payer, K. R., Calistri, N. L., Kimmerling, R. J., et al. (2020). Mass measurements during lymphocytic leukemia cell polyploidization decouple cell cycle- and cell size-dependent growth. Proc. Natl. Acad. Sci. U. S. A. 117, 15659–15665. doi:10.1073/pnas.1922197117
Neumann, F. R., and Nurse, P. (2007). Nuclear size control in fission yeast. J. Cell Biol. 179, 593–600. doi:10.1083/jcb.200708054
Nguyen, A., Yoshida, M., Goodarzi, H., and Tavazoie, S. F. (2016). Highly variable cancer subpopulations that exhibit enhanced transcriptome variability and metastatic fitness. Nat. Commun. 7, 11246. doi:10.1038/ncomms11246
Peng, W. C., Logan, C. Y., Fish, M., Anbarchian, T., Aguisanda, F., Alvarez-Varela, A., et al. (2018). Inflammatory cytokine TNFα promotes the long-term expansion of primary hepatocytes in 3D culture. Cell 175, 1607–1619.e15. doi:10.1016/j.cell.2018.11.012
Rafelski, S. M., Viana, M. P., Zhang, Y., Chan, Y. H. M., Thorn, K. S., Yam, P., et al. (2012). Mitochondrial network size scaling in budding yeast. Science 338, 822–824. doi:10.1126/science.1225720
Rubin, S. M., Sage, J., and Skotheim, J. M. (2020). Integrating old and new paradigms of G1/S control. Mol. Cell 80, 183–192. doi:10.1016/j.molcel.2020.08.020
Sage, J., Mulligan, G. J., Attardi, L. D., Miller, A., Chen, S., Williams, B., et al. (2000). Targeted disruption of the three Rb-related genes leads to loss of G1 control and immortalization. Genes Dev. 14, 3037–3050. doi:10.1101/gad.843200
Salvagno, G. L., Sanchis-Gomar, F., Picanza, A., and Lippi, G. (2015). Red blood cell distribution width: A simple parameter with multiple clinical applications. Crit. Rev. Clin. Lab. Sci. 52, 86–105. doi:10.3109/10408363.2014.992064
Sandlin, C. (2022). 3D characterization of cell size dysregulation in human lung adenocarcinoma reveals a network of fine processes connecting alveolar type 2 cells.
Savage, V. M., Allen, A. P., Brown, J. H., Gillooly, J. F., Herman, A. B., Woodruff, W. H., et al. (2007). Scaling of number, size, and metabolic rate of cells with body size in mammals. Proc. Natl. Acad. Sci. U. S. A. 104, 4718–4723. doi:10.1073/pnas.0611235104
Schmoller, K. M., Turner, J. J., Kõivomägi, M., and Skotheim, J. M. (2015). Dilution of the cell cycle inhibitor Whi5 controls budding-yeast cell size. Nature 526, 268–272. doi:10.1038/nature14908
Swaffer, M. P., Kim, J., Chandler-Brown, D., Langhinrichs, M., Marinov, G. K., Greenleaf, W. J., et al. (2021). Transcriptional and chromatin-based partitioning mechanisms uncouple protein scaling from cell size. Mol. Cell 81, 4861e7–4875.e7. doi:10.1016/j.molcel.2021.10.007
Tan, C., Ginzberg, M. B., Webster, R., Iyengar, S., Liu, S., Papadopoli, D., et al. (2021). Cell size homeostasis is maintained by CDK4-dependent activation of p38 MAPK. Dev. Cell 56, 1756e7–1769.e7. doi:10.1016/j.devcel.2021.04.030
Tzur, A., Kafri, R., LeBleu, V. S., Lahav, G., and Kirschner, M. W. (2009). Cell growth and size homeostasis in proliferating animal cells. Science 325, 167–171. doi:10.1126/science.1174294
Wei, Y., Wang, Y. G., Jia, Y., Li, L., Yoon, J., Zhang, S., et al. (2021). Liver homeostasis is maintained by midlobular zone 2 hepatocytes. Science 371, eabb1625. doi:10.1126/science.abb1625
Wirt, S. E., and Sage, J. (2010). p107 in the public eye: an Rb understudy and more. Cell Div. 5, 9. doi:10.1186/1747-1028-5-9
Wuestefeld, T., Pesic, M., Rudalska, R., Dauch, D., Longerich, T., Kang, T. W., et al. (2013). A Direct in vivo RNAi screen identifies MKK4 as a key regulator of liver regeneration. Cell 153, 389–401. doi:10.1016/j.cell.2013.03.026
Xie, S., Swaffer, M., and Skotheim, J. M. (2022). Eukaryotic cell size control and its relation to biosynthesis and senescence. Annu. Rev. Cell Dev. Biol. 38. doi:10.1146/annurev-cellbio-120219-040142
Zatulovskiy, E., and Skotheim, J. M. (2020). On the molecular mechanisms regulating animal cell size homeostasis. Trends Genet. 36, 360–372. doi:10.1016/j.tig.2020.01.011
Zatulovskiy, E., Zhang, S., Berenson, D. F., Topacio, B. R., and Skotheim, J. M. (2020). Cell growth dilutes the cell cycle inhibitor Rb to trigger cell division. Science 369, 466–471. doi:10.1126/science.aaz6213
Keywords: retinablastoma protein, cell size, mouse liver, G1/S, inhibitor dilution
Citation: Zhang S, Zatulovskiy E, Arand J, Sage J and Skotheim JM (2022) The cell cycle inhibitor RB is diluted in G1 and contributes to controlling cell size in the mouse liver. Front. Cell Dev. Biol. 10:965595. doi: 10.3389/fcell.2022.965595
Received: 09 June 2022; Accepted: 27 July 2022;
Published: 25 August 2022.
Edited by:
Eiman Aleem, University of Arizona, United StatesReviewed by:
Ariz Mohammad, Washington University in St. Louis, United StatesGabriel Neurohr, ETH Zürich, Switzerland
Ran Kafri, University of Toronto, Canada
Copyright © 2022 Zhang, Zatulovskiy, Arand, Sage and Skotheim. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jan M. Skotheim, skotheim@stanford.edu