A Modeling Study on How Cell Division Affects Properties of Epithelial Tissues Under Isotropic Growth

Cell proliferation affects both cellular geometry and topology in a growing tissue, and hence rules for cell division are key to understanding multicellular development. Epithelial cell layers have for long times been used to investigate how cell proliferation leads to tissue-scale properties, including organism-independent distributions of cell areas and number of neighbors. We use a cell-based two-dimensional tissue growth model including mechanics to investigate how different cell division rules result in different statistical properties of the cells at the tissue level. We focus on isotropic growth and division rules suggested for plant cells, and compare the models with data from the Arabidopsis shoot. We find that several division rules can lead to the correct distribution of number of neighbors, as seen in recent studies. In addition we find that when also geometrical properties are taken into account other constraints on the cell division rules result. We find that division rules acting in favor of equally sized and symmetrically shaped daughter cells can best describe the statistical tissue properties.     

Effects of Growth and Mutation on Pattern Formation in Tissues

In many developing tissues, neighboring cells enter different developmental pathways, resulting in a fine-grained pattern of different cell states. The most common mechanism that generates such patterns is lateral inhibition, for example through Delta-Notch coupling. In this work, we simulate growth of tissues consisting of a hexagonal arrangement of cells laterally inhibiting their neighbors. We find that tissue growth by cell division and cell migration tends to produce ordered patterns, whereas lateral growth leads to disordered, patchy patterns. Ordered patterns are very robust to mutations (gene silencing or activation) in single cells. In contrast, mutation in a cell of a disordered tissue can produce a larger and more widespread perturbation of the pattern. In tissues where ordered and disordered patches coexist, the perturbations spread mostly at boundaries between patches. If cell division occurs on time scales faster than the degradation time, disordered patches will appear. Our work suggests that careful experimental characterization of the disorder in tissues could pinpoint where and how the tissue is susceptible to large-scale damage even from single cell mutations.     

Shear and extensional rheology of cellulose/ionic liquid solutions

In this study, we characterize the shear and extensional rheology of dilute to semidilute solutions of cellulose in the ionic liquid 1-ethyl-3-methylimidazolium acetate (EMIAc). In steady shear flow, the semidilute solutions exhibit shear thinning, and the high-frequency complex modulus measured in small amplitude oscillatory shear flow exhibits the characteristic scaling expected for solutions of semiflexible chains. Flow curves of the steady shear viscosity plotted against shear rate closely follow the frequency dependence of the complex viscosity acquired using oscillatory shear, thus satisfying the empirical Cox-Merz rule. We use capillary thinning rheometry (CaBER) to characterize the relaxation times and apparent extensional viscosities of the semidilute cellulose solutions in a uniaxial extensional flow that mimics the dynamics encountered in the spin-line during fiber spinning processes. The apparent extensional viscosity and characteristic relaxation times of the semidilute cellulose/EMIAc solutions increase dramatically as the solutions enter the entangled concentration regime at which fiber spinning becomes viable.     

Viscoelastic properties of bovine orbital connective tissue and fat: constitutive models

Reported mechanical properties of orbital connective tissue and fat have been too sparse to model strain–stress relationships underlying biomechanical interactions in strabismus. We performed rheological tests to develop a multi-mode upper convected Maxwell (UCM) model of these tissues under shear loading. From 20 fresh bovine orbits, 30 samples of connective tissue were taken from rectus pulley regions and 30 samples of fatty tissues from the posterior orbit. Additional samples were defatted to determine connective tissue weight proportion, which was verified histologically. Mechanical testing in shear employed a triborheometer to perform: strain sweeps at 0.5–2.0 Hz; shear stress relaxation with 1% strain; viscometry at 0.01−0.5 s⁻¹ strain rate; and shear oscillation at 1% strain. Average connective tissue weight proportion was 98% for predominantly connective tissue and 76% for fatty tissue. Connective tissue specimens reached a long-term relaxation modulus of 668 Pa after 1,500 s, while corresponding values for fatty tissue specimens were 290 Pa and 1,100 s. Shear stress magnitude for connective tissue exceeded that of fatty tissue by five-fold. Based on these data, we developed a multi-mode UCM model with variable viscosities and time constants, and a damped hyperelastic response that accurately described measured properties of both connective and fatty tissues. Model parameters differed significantly between the two tissues. Viscoelastic properties of predominantly connective orbital tissues under shear loading differ markedly from properties of orbital fat, but both are accurately reflected using UCM models. These viscoelastic models will facilitate realistic global modeling of EOM behavior in binocular alignment and strabismus.     

Can tissue surface tension drive somite formation?

The prevailing model of somitogenesis supposes that the presomitic mesoderm is segmented into somites by a clock and wavefront mechanism. During segmentation, mesenchymal cells undergo compaction, followed by a detachment of the presumptive somite from the rest of the presomitic mesoderm and the subsequent morphological changes leading to rounded somites. We investigate the possibility that minimization of tissue surface tension drives the somite sculpting processes. Given the time in which somite formation occurs and the high bulk viscosities of tissues, we find that only small changes in shape and form of tissue typically occur through cell movement driven by tissue surface tension. This is particularly true for somitogenesis in the zebrafish. Hence it is unlikely that such processes are the sole and major driving force behind somite formation. We propose a simple chemotactic mechanism that together with heightened adhesion can account for the morphological changes in the time allotted for somite formation.     

Interpreting CFSE Obtained Division Histories of B Cells in Vitro with Smith–Martin and Cyton Type Models

The fluorescent dye carboxyfluorescin diacetate succinimidyl ester (CFSE) classifies proliferating cell populations into groups according to the number of divisions each cell has undergone (i.e., its division class). The pulse labeling of cells with radioactive thymidine provides a means to determine the distribution of times of entry into the first cell division. We derive in analytic form the number of cells in each division class as a function of time using the cyton approach that utilizes independent stochastic distributions for the time to divide and the time to die. We confirm that our analytic form for the number of cells in each division class is consistent with the numerical solution of a set of delay differential equations representing the generalized Smith–Martin model with cell death rates depending on the division class. Choosing the distribution of time to the first division to fit thymidine labeling data for B cells stimulated in vitro with lipopolysaccharide (LPS) and either with or without interleukin-4 (IL-4), we fit CFSE data to determine the dependence of B cell kinetic parameters on the presence of IL-4. We find when IL-4 is present, a greater proportion of cells are recruited into division with a longer average time to first division. The most profound effect of the presence of IL-4 was decreased death rates for smaller division classes, which supports a role of IL-4 in the protection of B cells from apoptosis.     

Highly protective alkalinization by ammonia vapor diffusion in viscosimetric DNA damage assessment

A method for the measurement of viscosities correlated to DNA alterations in alkaline homogenate suspensions is described. The alkaline pH shift to afford cell lysis, DNA unfolding, and denaturation is attained by gaseous ammonia diffusion, thus avoiding shear stress from mechanical mixing. At the same time a stabilizing density gradient is established. This solution is run through a plastic measuring tube that is wide enough to minimize the influence of uneven swelling of the lysing DNA-containing components. Flow times under a carefully controlled water head are registered, and their ratios to control solutions are evaluated. The relative viscosities show a strong and irreversible dependence on shear and on DNase treatment and therefore are considered as essentially DNA derived. The time dependence of the lysate viscosities with and without the DNA-damaging agent bleomycin is given and the dose:activity curves of this agent with sponge homogenates from two orders of Porifera are given with their 50% effective concentration values. The dose:activity curve of an extract from a polluted marine point source is demonstrated. The concentration changes in sponges exposed at differently polluted marine sites are shown. The idea of alkalinization through gaseous diffusion in conjunction with a simple measuring device has already proven a sensitive, reliable, and specific tool in the assessment of DNA damage produced under both laboratory and field conditions.     

A rheological study on plasma lipoproteins (author's transl)]

Lipoproteins of the types LpB and LP-X are studied in a microviscometer to measure intrinsic viscosity. Up to a shear rate of 1700s-1 no shear dependence of viscosity is observed. Intrinsic viscosities are 3.5 and 4.0 for LpB, and 8.5 ml/g for LP-X. Density measurements are used to calculate apparent specific volumes. The results are compatible with a spherical model. An upper limit of 0.45 is estimated for the hydration.     

Assessing the use of ellipsoidal microparticles for determining lipid membrane viscosity

The viscosity of lipid membranes sets the timescales of membrane-associated motions, whether driven or diffusive, and therefore influences the dynamics of a wide range of cellular processes. Techniques to measure membrane viscosity remain sparse, however, and reported measurements to date, even of similar systems, give viscosity values that span orders of magnitude. To address this, we improve a method based on measuring both the rotational and translational diffusion of membrane-anchored microparticles and apply this approach and one based on tracking the motion of phase-separated lipid domains to the same system of phase-separated giant vesicles. We find good agreement between the two methods, with inferred viscosities within a factor of two of each other. Our single-particle tracking technique uses ellipsoidal microparticles, and we show that the extraction of physically meaningful viscosity values from their motion requires consideration of their anisotropic shape. The validation of our method on phase-separated membranes makes possible its application to other systems, which we demonstrate by measuring the viscosity of bilayers composed of lipids with different chain lengths ranging from 14 to 20 carbon atoms, revealing a very weak dependence of two-dimensional viscosity on lipid size. The experimental and analysis methods described here should be generally applicable to a variety of membrane systems, both reconstituted and cellular.     

Cell death and tissue remodeling in planarian regeneration

Many long-lived organisms, including humans, can regenerate some adult tissues lost to physical injury or disease. Much of the previous research on mechanisms of regeneration has focused on adult stem cells, which give rise to new tissue necessary for the replacement of missing body parts. Here we report that apoptosis of differentiated cells complements stem cell division during regeneration in the planarian Schmidtea mediterranea. Specifically, we developed a whole-mount TUNEL assay that allowed us to document two dramatic increases in the rate of apoptosis following amputation—an initial localized response near the wound site and a subsequent systemic response that varies in magnitude depending on the type of fragment examined. The latter cell death response can be induced in uninjured organs, occurs in the absence of planarian stem cells, and can also be triggered by prolonged starvation. Taken together, our results implicate apoptosis in the restoration of proper anatomical scale and proportion through remodeling of existing tissues. We also report results from initial mechanistic studies of apoptosis in planarians, which revealed that a S. mediterranea homolog of the antiapoptotic gene BCL2 is required for cell survival in adult animals. We propose that apoptosis is a central mechanism working in concert with stem cell division to restore anatomical form and function during metazoan regeneration.     

Non‐Newtonian rheology in suspension cell cultures significantly impacts bioreactor shear stress quantification

The fields of regenerative medicine and tissue engineering require large‐scale manufacturing of stem cells for both therapy and recombinant protein production, which is often achieved by culturing cells in stirred suspension bioreactors. The rheology of cell suspensions cultured in stirred suspension bioreactors is critical to cell growth and protein production, as elevated exposure to shear stress has been linked to changes in growth kinetics and genetic expression for many common cell types. Currently, little is understood on the rheology of cell suspensions cultured in stirred suspension bioreactors. In this study, we present the impact of three common cell culture parameters, serum content, cell presence, and culture age, on the rheology of a model cell line cultured in stirred suspension bioreactors. The results reveal that cultures containing cells, serum, or combinations thereof are highly shear thinning, whereas conditioned and unconditioned culture medium without serum are both Newtonian. Non‐Newtonian viscosity was modeled using a Sisko model, which provided insight on structural mechanisms driving the rheological behavior of these cell suspensions. A comparison of shear stress estimated by using Newtonian and Sisko relationships demonstrated that assuming Newtonian viscosity underpredicts both mean and maximum shear stress in stirred suspension bioreactors. Non‐Newtonian viscosity models reported maximum shear stresses exceeding those required to induce changes in genetic expression in common cell types, whereas Newtonian models did not. These findings indicate that traditional shear stress quantification of cell or serum suspensions is inadequate and that shear stress quantification methods based on non‐Newtonian viscosity must be developed to accurately quantify shear stress.     

Modeling cell proliferation for simulating three-dimensional tissue morphogenesis based on a reversible network reconnection framework

Tissue morphogenesis in multicellular organisms is accompanied by proliferative cell behaviors: cell division (increase in cell number after each cell cycle) and cell growth (increase in cell volume during each cell cycle). These proliferative cell behaviors can be regulated by multicellular dynamics to achieve proper tissue sizes and shapes in three-dimensional (3D) space. To analyze multicellular dynamics, a reversible network reconnection (RNR) model has been suggested, in which each cell shape is expressed by a single polyhedron. In this study, to apply the RNR model to simulate tissue morphogenesis involving proliferative cell behaviors, we model cell proliferation based on a RNR model framework. In this model, cell division was expressed by dividing a polyhedron at a planar surface for which cell division behaviors were characterized by three quantities: timing, intracellular position, and normal direction of the dividing plane. In addition, cell growth was expressed by volume growth as a function of individual cell times within their respective cell cycles. Numerical simulations using the proposed model showed that tissues grew during successive cell divisions with several cell cycle times. During these processes, the cell number in tissues increased while maintaining individual cell size and shape. Furthermore, tissue morphology dramatically changed based on different regulations of cell division directions. Thus, the proposed model successfully provided a basis for expressing proliferative cell behaviors during morphogenesis based on a RNR model framework.     

Quantitative analysis of T cell homeostatic proliferation

T cell homeostatic proliferation occurs on transfer of T cells into lymphopenic recipients; transferred cells undergo several rounds of division in the absence of specific antigen stimulation. For a quantitative analysis of this phenomenon, we applied a mathematical method to describe proliferating T cells to match peak distributions from actual CFSE dilution data. For in vitro stimulation of T cells with anti-CD3/anti-CD28, our simulation confirmed a high proportion of cells entering cell cycle with a low proportion undergoing apoptosis. When applied to homeostatic proliferation, it described striking differences in CD4 and CD8 T cell proliferation rates, and accurately predicted that successive divisions were accompanied by higher rates of apoptosis, limiting the accumulation of proliferating cells. Thus, the presence of multiple CFSE dilution peaks cannot be considered equivalent to lymphocyte expansion. Finally, genetic effects were identified that may help explain links between homeostatic proliferation and autoimmunity.     

Reusing and composing models of cell fate regulation of human bone precursor cells

In order to treat osteoporosis and other bone mass disorders it is necessary to understand the regulatory processes that control the cell fate decisions responsible for going from bone precursor cells to bone tissue. Many processes interact to regulate cell division, differentiation and apoptosis. There are models for these basic processes, but not for their interactions. In this work we use the theory of switched systems, reuse and composition of validated models to describe the cell fate decisions leading to bone and fat formation. We describe the differentiation of osteo-adipo progenitor cells by composing its model with differentiation stimuli. We use the activation of the Wnt pathway as stimulus to osteoblast lineage, including regulation of cell division and apoptosis. This model is our first step to simulate physiological responses in silico to treatments for bone mass disorders.     

Co-Ordination of Cell Division and Tissue Expansion in Sunflower, Tobacco, and Pea Leaves: Dependence or Independence of Both Processes?

Abstract Temporal analyses of cell division and tissue expansion in pea, tobacco, and sunflower leaves reveal that both processes follow similar patterns during leaf development. Relative cell division and relative tissue expansion rates are maximal and constant during early leaf development, but they decline later. In contrast, relative cell expansion rate follows a bell-shaped curve during leaf growth. Cell division and tissue expansion have common responses to temperature, intercepted radiation, and water deficit. As a consequence, final leaf area and cell number remain highly correlated throughout a large range of environmental conditions for these different plant species, indicating that cell division and tissue expansion are co-ordinated during leaf development. This co-ordination between processes has long been explained by dependence between both processes. Most studies on dicotyledonous leaf development indicate that leaf expansion rate depends on the number of cells in the leaf. We tested this hypothesis with a large range of environmental conditions and different plant species. Accordingly, we found a strong correlation between both absolute leaf expansion rate and leaf cell number. However, we showed that this relationship is not necessarily causal because it can be simulated by the hypothesis of independence between cell division and tissue expansion according to Green's theory of growth (1976). Received 23 February 2000; accepted 3 March 2000     

JNK Controls the Onset of Mitosis in Planarian Stem Cells and Triggers Apoptotic Cell Death Required for Regeneration and Remodeling

Regeneration of lost tissues depends on the precise interpretation of molecular signals that control and coordinate the onset of proliferation, cellular differentiation and cell death. However, the nature of those molecular signals and the mechanisms that integrate the cellular responses remain largely unknown. The planarian flatworm is a unique model in which regeneration and tissue renewal can be comprehensively studied in vivo. The presence of a population of adult pluripotent stem cells combined with the ability to decode signaling after wounding enable planarians to regenerate a complete, correctly proportioned animal within a few days after any kind of amputation, and to adapt their size to nutritional changes without compromising functionality. Here, we demonstrate that the stress-activated c-jun–NH₂–kinase (JNK) links wound-induced apoptosis to the stem cell response during planarian regeneration. We show that JNK modulates the expression of wound-related genes, triggers apoptosis and attenuates the onset of mitosis in stem cells specifically after tissue loss. Furthermore, in pre-existing body regions, JNK activity is required to establish a positive balance between cell death and stem cell proliferation to enable tissue renewal, remodeling and the maintenance of proportionality. During homeostatic degrowth, JNK RNAi blocks apoptosis, resulting in impaired organ remodeling and rescaling. Our findings indicate that JNK-dependent apoptotic cell death is crucial to coordinate tissue renewal and remodeling required to regenerate and to maintain a correctly proportioned animal. Hence, JNK might act as a hub, translating wound signals into apoptotic cell death, controlled stem cell proliferation and differentiation, all of which are required to coordinate regeneration and tissue renewal.     

Mechanisms involved in tissue-specific apopotosis regulated by glucocorticoids

Physiological cell turnover is under the control of a sharp and dynamic balance of different homeostatic mechanisms such as the equilibrium between cell proliferation and cell death. These mechanisms play an important role in maintaining normal tissue function and architecture. It is well known that apoptosis is the prevalent mode of physiological cell loss in most tissues. Steroid hormones like glucocorticoids have been identified as key signals controlling cell turnover by modulating programmed cell death in a tissue- and cell-specific manner. In this sense, several reports have demonstrated that glucocorticoids are able to induce apoptosis in cells of the hematopoietic system such as monocytes, macrophages, and T lymphocytes. In contrast, they protect against apoptotic signals evoked by cytokines, cAMP, tumor suppressors, in glandular cells such as the mammary gland epithelia, endometrium, hepatocytes, ovarian follicular cells, and fibroblasts. Although several studies have provided significant information on hormone-dependent apoptosis in an specific tissue, a clearly defined pathway that mediates cell death in response to glucocorticoids in different cell types is still misunderstood. The scope of this review is held to those mechanisms by which glucocorticoids control apoptosis, emphasizing tissue-specific expression of genes that are involved in the apoptotic pathway.