SimuCell3D: three-dimensional simulation of tissue mechanics with cell polarization.

The three-dimensional (3D) organization of cells determines tissue function and integrity, and changes markedly in development and disease. Cell-based simulations have long been used to define the underlying mechanical principles. However, high computational costs have so far limited simulations to either simplified cell geometries or small tissue patches. Here, we present SimuCell3D, an efficient open-source program to simulate large tissues in three dimensions with subcellular resolution, growth, proliferation, extracellular matrix, fluid cavities, nuclei and non-uniform mechanical properties, as found in polarized epithelia. Spheroids, vesicles, sheets, tubes and other tissue geometries can readily be imported from microscopy images and simulated to infer biomechanical parameters. Doing so, we show that 3D cell shapes in layered and pseudostratified epithelia are largely governed by a competition between surface tension and intercellular adhesion. SimuCell3D enables the large-scale in silico study of 3D tissue organization in development and disease at a great level of detail.

Simulating 3D Cell Shape with the Cellular Potts Model.

Computer simulations have become a widely used method for the field of mechanobiology. An important question is whether one can predict the shape and forces of cells as a function of the extracellular environment. Different types of models have been described before to simulate cell and tissue shapes in structured environments. In this chapter, we give a brief overview of commonly used models and then describe the Cellular Potts Model, a lattice-based modelling framework, in more detail. We provide a hands-on guide on how to build a model that simulates the shape of a single cell on a micropattern in three dimensions in different open source software packages using the Cellular Potts framework. A simulation is set up with an initial configuration of generalized cells that change shape and position due to an energy function that incorporates cellular volume and surface area constraints as well as interaction energies between the generalized cells.

Actomyosin fibers DApPLE epithelial apical junctions.

Epithelial cell morphology is essential for cellular homeostasis, but the mechanisms by which cell shape is established remain unclear. In this study, Marivin et al. (2022. J. Cell Biol.https://doi.org/10.1083/jcb.202111002) identify DAPLE as a linker between polarity complexes and the actomyosin network at apical junctions. By recruiting CD2P and activating Gαßγ-mediated RhoA signaling, DAPLE ensures proper cell shape and function.

Metastasising Fibroblasts Show an HDAC6-Dependent Increase in Migration Speed and Loss of Directionality Linked to Major Changes in the Vimentin Interactome.

Metastasising cells express the intermediate filament protein vimentin, which is used to diagnose invasive tumours in the clinic. We aimed to clarify how vimentin regulates the motility of metastasising fibroblasts. STED super-resolution microscopy, live-cell imaging and quantitative proteomics revealed that oncogene-expressing and metastasising fibroblasts show a less-elongated cell shape, reduced cell spreading, increased cell migration speed, reduced directionality, and stronger coupling between these migration parameters compared to normal control cells. In total, we identified and compared 555 proteins in the vimentin interactome. In metastasising cells, the levels of keratin 18 and Rab5C were increased, while those of actin and collagen were decreased. Inhibition of HDAC6 reversed the shape, spreading and migration phenotypes of metastasising cells back to normal. Inhibition of HDAC6 also decreased the levels of talin 1, tropomyosin, Rab GDI ß, collagen and emilin 1 in the vimentin interactome, and partially reversed the nanoscale vimentin organisation in oncogene-expressing cells. These findings describe the changes in the vimentin interactome and nanoscale distribution that accompany the defective cell shape, spreading and migration of metastasising cells. These results support the hypothesis that oncogenes can act through HDAC6 to regulate the vimentin binding of the cytoskeletal and cell-extracellular matrix adhesion components that contribute to the defective motility of metastasising cells.

A deep generative model of 3D single-cell organization.

We introduce a framework for end-to-end integrative modeling of 3D single-cell multi-channel fluorescent image data of diverse subcellular structures. We employ stacked conditional ß-variational autoencoders to first learn a latent representation of cell morphology, and then learn a latent representation of subcellular structure localization which is conditioned on the learned cell morphology. Our model is flexible and can be trained on images of arbitrary subcellular structures and at varying degrees of sparsity and reconstruction fidelity. We train our full model on 3D cell image data and explore design trade-offs in the 2D setting. Once trained, our model can be used to predict plausible locations of structures in cells where these structures were not imaged. The trained model can also be used to quantify the variation in the location of subcellular structures by generating plausible instantiations of each structure in arbitrary cell geometries. We apply our trained model to a small drug perturbation screen to demonstrate its applicability to new data. We show how the latent representations of drugged cells differ from unperturbed cells as expected by on-target effects of the drugs.

APOE4 confers transcriptomic and functional alterations to primary mouse microglia.

Common genetic variants in more than forty loci modulate risk for Alzheimer's disease (AD). AD risk alleles are enriched within enhancers active in myeloid cells, suggesting that microglia, the brain-resident macrophages, may play a key role in the etiology of AD. A major genetic risk factor for AD is Apolipoprotein E (APOE) genotype, with the ε4/ε4 (E4) genotype increasing risk for AD by approximately 15 fold compared to the most common ε3/ε3 (E3) genotype. However, the impact of APOE genotype on microglial function has not been thoroughly investigated. To address this, we cultured primary microglia from mice in which both alleles of the mouse Apoe gene have been humanized to encode either human APOE ε3 or APOE ε4. Relative to E3 microglia, E4 microglia exhibit altered morphology, increased endolysosomal mass, increased cytokine/chemokine production, and increased lipid and lipid droplet accumulation at baseline. These changes were accompanied by decreased translation and increased phosphorylation of eIF2ɑ and eIF2ɑ-kinases that participate in the integrated stress response, suggesting that E4 genotype leads to elevated levels of cellular stress in microglia relative to E3 genotype. Using live-cell imaging and flow cytometry, we also show that E4 microglia exhibited increased phagocytic uptake of myelin and other substrates compared to E3 microglia. While transcriptomic profiling of myelin-challenged microglia revealed a largely overlapping response profile across genotypes, differential enrichment of genes in interferon signaling, extracellular matrix and translation-related pathways was identified in E4 versus E3 microglia both at baseline and following myelin challenge. Together, our results suggest E4 genotype confers several important functional alterations to microglia even prior to myelin challenge, providing insight into the molecular and cellular mechanisms by which APOE4 may increase risk for AD.

Primordial mimicry induces morphological change in Escherichia coli.

The morphology of primitive cells has been the subject of extensive research. A spherical form was commonly presumed in prebiotic studies but lacked experimental evidence in living cells. Whether and how the shape of living cells changed are unclear. Here we exposed the rod-shaped bacterium Escherichia coli to a resource utilization regime mimicking a primordial environment. Oleate was given as an easy-to-use model prebiotic nutrient, as fatty acid vesicles were likely present on the prebiotic Earth and might have been used as an energy resource. Six evolutionary lineages were generated under glucose-free but oleic acid vesicle (OAV)-rich conditions. Intriguingly, fitness increase was commonly associated with the morphological change from rod to sphere and the decreases in both the size and the area-to-volume ratio of the cell. The changed cell shape was conserved in either OAVs or glucose, regardless of the trade-offs in carbon utilization and protein abundance. Highly differentiated mutations present in the genome revealed two distinct strategies of adaption to OAV-rich conditions, i.e., either directly targeting the cell wall or not. The change in cell morphology of Escherichia coli for adapting to fatty acid availability supports the assumption of the primitive spherical form.

Viral adventitious agent detection using laser force cytology: Intrinsic cell property changes with infection and comparison to in vitro testing.

Adventitious agent testing in biomanufacturing requires assays of broad detection capability to screen for as many infectious agents as possible. The current gold standard for general infectious adventitious virus screening is the in vitro assay in which test articles are cultured onto a panel of different cell lines and observed for cytopathic effect (CPE). However, this assay is inherently subjective due to the nature of visual observation of cell morphology and labor and time intensive, requiring highly trained personnel to identify CPE. Laser force cytology (LFC) is an alternative, automated analytical method that uses a combination of optical and fluidic forces along with imaging to objectively and quantitatively assess CPE in cell culture. Importantly, because LFC uses no labels or antibodies, the assay is appropriate for general adventitious agent testing. Using LFC, changes in cellular features associated with virally infected cells were identified using principal component analysis. Using these features of infected cells, the sensitivity and earliness of detection with LFC was directly compared with the in vitro assay for a diverse panel of viruses incubated with chinese hamster ovary (CHO), Vero, and Medical Research Council cell strain 5 (MRC-5) cells. LFC detected viral infection with a sensitivity equal to the in vitro assay on average, but in certain virus and cell combinations including mouse minute virus (MMV) and reovirus 3 in CHO cells, detection was 4 days earlier and for MMV, the limit of detection was 10-fold lower. Overall, these results demonstrate the ability of LFC to serve as a biopharmaceutical adventitious agent testing methodology with sensitivity equivalent to the in vitro assay, but in an objective and automated manner.

Essential Role of Sptan1 in Cochlear Hair Cell Morphology and Function Via Focal Adhesion Signaling.

Hearing loss is the most common human sensory deficit. Hearing relies on stereocilia, inserted into the cuticular plate of hair cells (HCs), where they play an important role in the perception of sound and its transmission. Although numerous genes have been associated with hearing loss, the function of many hair cell genes has yet to be elucidated. Herein, we focused on nonerythroid spectrin αII (SPTAN1), abundant in the cuticular plate, surrounding the rootlets of stereocilia and along the plasma membrane. Interestingly, mice with HC-specific Sptan1 knockout exhibited rapid deafness, abnormal formation of stereocilia and cuticular plates, and loss of HCs from middle and apical turns of the cochlea during early postnatal stages. Additionally, Sptan1 deficiency led to the decreased spreading of House Ear Institute-Organ of Corti 1 cells, and induced abnormal formation of focal adhesions and integrin signaling in mouse HCs. Altogether, our findings highlight SPTAN1 as a critical molecule for HC stereocilia morphology and auditory function via regulation of focal adhesion signaling.

Combined effect of cell geometry and polarity domains determines the orientation of unequal division.

Cell division orientation is thought to result from a competition between cell geometry and polarity domains controlling the position of the mitotic spindle during mitosis. Depending on the level of cell shape anisotropy or the strength of the polarity domain, one dominates the other and determines the orientation of the spindle. Whether and how such competition is also at work to determine unequal cell division (UCD), producing daughter cells of different size, remains unclear. Here, we show that cell geometry and polarity domains cooperate, rather than compete, in positioning the cleavage plane during UCDs in early ascidian embryos. We found that the UCDs and their orientation at the ascidian third cleavage rely on the spindle tilting in an anisotropic cell shape, and cortical polarity domains exerting different effects on spindle astral microtubules. By systematically varying mitotic cell shape, we could modulate the effect of attractive and repulsive polarity domains and consequently generate predicted daughter cell size asymmetries and position. We therefore propose that the spindle position during UCD is set by the combined activities of cell geometry and polarity domains, where cell geometry modulates the effect of cortical polarity domain(s).

Factor quinolinone inhibitors alter cell morphology and motility by destabilizing interphase microtubules.

Factor quinolinone inhibitors are promising anti-cancer compounds, initially characterized as specific inhibitors of the oncogenic transcription factor LSF (TFCP2). These compounds exert anti-proliferative activity at least in part by disrupting mitotic spindles. Herein, we report additional interphase consequences of the initial lead compound, FQI1, in two telomerase immortalized cell lines. Within minutes of FQI1 addition, the microtubule network is disrupted, resulting in a substantial, although not complete, depletion of microtubules as evidenced both by microtubule sedimentation assays and microscopy. Surprisingly, this microtubule breakdown is quickly followed by an increase in tubulin acetylation in the remaining microtubules. The sudden breakdown and partial depolymerization of the microtubule network precedes FQI1-induced morphological changes. These involve rapid reduction of cell spreading of interphase fetal hepatocytes and increase in circularity of retinal pigment epithelial cells. Microtubule depolymerization gives rise to FH-B cell compaction, as pretreatment with taxol prevents this morphological change. Finally, FQI1 decreases the rate and range of locomotion of interphase cells, supporting an impact of FQI1-induced microtubule breakdown on cell motility. Taken together, our results show that FQI1 interferes with microtubule-associated functions in interphase, specifically cell morphology and motility.

Origami: Single-cell 3D shape dynamics oriented along the apico-basal axis of folding epithelia from fluorescence microscopy data.

A common feature of morphogenesis is the formation of three-dimensional structures from the folding of two-dimensional epithelial sheets, aided by cell shape changes at the cellular-level. Changes in cell shape must be studied in the context of cell-polarised biomechanical processes within the epithelial sheet. In epithelia with highly curved surfaces, finding single-cell alignment along a biological axis can be difficult to automate in silico. We present 'Origami', a MATLAB-based image analysis pipeline to compute direction-variant cell shape features along the epithelial apico-basal axis. Our automated method accurately computed direction vectors denoting the apico-basal axis in regions with opposing curvature in synthetic epithelia and fluorescence images of zebrafish embryos. As proof of concept, we identified different cell shape signatures in the developing zebrafish inner ear, where the epithelium deforms in opposite orientations to form different structures. Origami is designed to be user-friendly and is generally applicable to fluorescence images of curved epithelia.

Programming cell growth into different cluster shapes using diffusible signals.

Advances in genetic engineering technologies have allowed the construction of artificial genetic circuits, which have been used to generate spatial patterns of differential gene expression. However, the question of how cells can be programmed, and how complex the rules need to be, to achieve a desired tissue morphology has received less attention. Here, we address these questions by developing a mathematical model to study how cells can collectively grow into clusters with different structural morphologies by secreting diffusible signals that can influence cellular growth rates. We formulate how growth regulators can be used to control the formation of cellular protrusions and how the range of achievable structures scales with the number of distinct signals. We show that a single growth inhibitor is insufficient for the formation of multiple protrusions but may be achieved with multiple growth inhibitors, and that other types of signals can regulate the shape of protrusion tips. These examples illustrate how our approach could potentially be used to guide the design of regulatory circuits for achieving a desired target structure.

Auxin-dependent control of cytoskeleton and cell shape regulates division orientation in the Arabidopsis embryo.

Premitotic control of cell division orientation is critical for plant development, as cell walls prevent extensive cell remodeling or migration. While many divisions are proliferative and add cells to existing tissues, some divisions are formative and generate new tissue layers or growth axes. Such formative divisions are often asymmetric in nature, producing daughters with different fates. We have previously shown that, in the Arabidopsis thaliana embryo, developmental asymmetry is correlated with geometric asymmetry, creating daughter cells of unequal volume. Such divisions are generated by division planes that deviate from a default "minimal surface area" rule. Inhibition of auxin response leads to reversal to this default, yet the mechanisms underlying division plane choice in the embryo have been unclear. Here, we show that auxin-dependent division plane control involves alterations in cell geometry, but not in cell polarity axis or nuclear position. Through transcriptome profiling, we find that auxin regulates genes controlling cell wall and cytoskeleton properties. We confirm the involvement of microtubule (MT)-binding proteins in embryo division control. Organization of both MT and actin cytoskeleton depends on auxin response, and genetically controlled MT or actin depolymerization in embryos leads to disruption of asymmetric divisions, including reversion to the default. Our work shows how auxin-dependent control of MT and actin cytoskeleton properties interacts with cell geometry to generate asymmetric divisions during the earliest steps in plant development.

Understanding cytoskeletal avalanches using mechanical stability analysis.

Eukaryotic cells are mechanically supported by a polymer network called the cytoskeleton, which consumes chemical energy to dynamically remodel its structure. Recent experiments in vivo have revealed that this remodeling occasionally happens through anomalously large displacements, reminiscent of earthquakes or avalanches. These cytoskeletal avalanches might indicate that the cytoskeleton's structural response to a changing cellular environment is highly sensitive, and they are therefore of significant biological interest. However, the physics underlying "cytoquakes" is poorly understood. Here, we use agent-based simulations of cytoskeletal self-organization to study fluctuations in the network's mechanical energy. We robustly observe non-Gaussian statistics and asymmetrically large rates of energy release compared to accumulation in a minimal cytoskeletal model. The large events of energy release are found to correlate with large, collective displacements of the cytoskeletal filaments. We also find that the changes in the localization of tension and the projections of the network motion onto the vibrational normal modes are asymmetrically distributed for energy release and accumulation. These results imply an avalanche-like process of slow energy storage punctuated by fast, large events of energy release involving a collective network rearrangement. We further show that mechanical instability precedes cytoquake occurrence through a machine-learning model that dynamically forecasts cytoquakes using the vibrational spectrum as input. Our results provide a connection between the cytoquake phenomenon and the network's mechanical energy and can help guide future investigations of the cytoskeleton's structural susceptibility.

Impaired semen quality, an increase of sperm morphological defects and DNA fragmentation associated with environmental pollution in urban population of young men from Western Siberia, Russia.

Poor sperm morphology and an elevated DNA fragmentation level are considered to be related to spermiogenesis malfunctions as a result of genetic mutations and effects of environmental factors, including industrial pollution. Standardized cross-sectional population studies of sperm morphology defects and sperm DNA fragmentation, especially in regions with increased environmental pollution may be helpful to investigate an influence of industrial pollution and other population-related factors on spermiogenesis process. The aim of present study was to estimate an influence industrial pollution on sperm morphogenesis and sperm DNA fragmentation in men from the general population of the Western Siberia. The Novosibirsk and Kemerovo cities are located to same climatic conditions in Western Siberia but the Kemerovo city is characterized by increased environmental pollution especially by particulate matter (PM). The male volunteers living in Novosibirsk (n = 278) and Kemerovo (n = 258) were enrolled. Percentages of sperm morphological defects are counted after staining native ejaculate smears by Diff-Quick kits. DNA fragmentation was estimated by a SCSA technique. The residents of Kemerovo were characterized by lowered sperm count and sperm motility, elevated DNA fragmentation, poor sperm morphology and increased incidence of morphological effects of head (pyriform, elongated, round, abnormal acrosome and vacuolated chromatine), asymmetrical neck insertion and excess residual cytoplasm. Moreover, elevated DNA fragmentation was associated with lowered sperm count, sperm motility and increased percentages of several sperm morphology defects, with the place of residence affecting the relationships between conventional semen parameters, sperm morphology and DNA fragmentations. Our study suggests that excessive sperm head elongation and impaired acrosome formation can contribute to sperm morphology deterioration in men from polluted areas. Regional features in the relationships between sperm morphology, sperm count and DNA fragmentation were shown, suggesting an importance of studying sperm morphology pattern in men from different regions.

Symmetry at the 4-Cell Stage Is Associated with Embryo Aneuploidy.

The purpose of this study was to determine if morphometric parameters that can be measured quantitatively using a time-lapse embryo incubator are associated with aneuploidy. Embryos cultured in a time-lapse incubator and assessed with preimplantation genetic testing for aneuploidy (PGT-A) were analyzed retrospectively. Morphokinetic analysis included timing of cell divisions. Quantitative morphometric measurements included the distance between the second and first polar body, zona pellucida thickness at the pronuclear stage and at the 2-cell stage, and blastomere area at the 2- and 4-cell stages. Symmetry at the 2-cell stage was determined by percent difference between blastomeres; symmetry at the 4-cell stage was the percent difference between the smallest and largest blastomeres. Maternal age, blastocyst grade and day of biopsy were recorded. Euploid embryo characteristics were compared to aneuploid embryos. A receiver operating characteristic (ROC) curve was used to evaluate cell symmetry as a predictor of aneuploidy. Embryos (n = 182) from 21 patients (age 22-43; median = 34) were analyzed. Of the 182 embryos, 45% were euploid. Euploid and aneuploid embryos had similar morphokinetics and morphometry across many measures. As expected, age and blastocyst grade were associated with embryo ploidy. It was notable that, additionally, symmetry at the 4-cell stage (27% vs 31%, p = 0.01) was also associated with embryo ploidy. The optimized cutoff from the ROC curve to predict aneuploidy was determined to be 21%. Embryos with > 21% asymmetry at the 4-cell stage had high rates of aneuploidy while morphokinetic parameters were similar. In conclusion, this suggests that embryo selection models using time-lapse parameters would improve if they incorporate cleavage-stage morphometrics.

Tumor-initiating stem cell shapes its microenvironment into an immunosuppressive barrier and pro-tumorigenic niche.

Tumor-initiating stem cells (TSCs) are critical for drug resistance and immune escape. However, the mutual regulations between TSC and tumor microenvironment (TME) remain unclear. Using DNA-label retaining, single-cell RNA sequencing (scRNA-seq), and other approaches, we investigated intestinal adenoma in response to chemoradiotherapy (CRT), thus identifying therapy-resistant TSCs (TrTSCs). We find bidirectional crosstalk between TSCs and TME using CellPhoneDB analysis. An intriguing finding is that TSCs shape TME into a landscape that favors TSCs for immunosuppression and propagation. Using adenoma-organoid co-cultures, niche-cell depletion, and lineaging tracing, we characterize a functional role of cyclooxygenase-2 (Cox-2)-dependent signaling, predominantly occurring between tumor-associated monocytes and macrophages (TAMMs) and TrTSCs. We show that TAMMs promote TrTSC proliferation through prostaglandin E2 (PGE2)-PTGER4(EP4) signaling, which enhances ß-catenin activity via AKT phosphorylation. Thus, our study shows that the bidirectional crosstalk between TrTSC and TME results in a pro-tumorigenic and immunosuppressive contexture.

A PtdIns(3,4,5)P3 dispersal switch engages cell ratcheting at specific cell surfaces.

Force generation in epithelial tissues is often pulsatile, with actomyosin networks generating contractile forces before cyclically disassembling. This pulsed nature of cytoskeletal forces implies that there must be ratcheting mechanisms that drive processive transformations in cell shape. Previous work has shown that force generation is coordinated with endocytic remodeling; however, how ratcheting becomes engaged at specific cell surfaces remains unclear. Here, we report that PtdIns(3,4,5)P3 is a critical lipid-based cue for ratcheting engagement. The Sbf RabGEF binds to PIP3, and disruption of PIP3 reveals a dramatic switching behavior in which medial ratcheting is activated and epithelial cells begin globally constricting apical surfaces. PIP3 enrichments are developmentally regulated, with mesodermal cells having high apical PIP3 while germband cells have higher interfacial PIP3. Finally, we show that JAK/STAT signaling constitutes a second pathway that combinatorially regulates Sbf/Rab35 recruitment. Our results elucidate a complex lipid-dependent regulatory machinery that directs ratcheting engagement in epithelial tissues.

Viral Manipulation of a Mechanoresponsive Signaling Axis Disassembles Processing Bodies.

Processing bodies (PBs) are ribonucleoprotein granules important for cytokine mRNA decay that are targeted for disassembly by many viruses. Kaposi's sarcoma-associated herpesvirus is the etiological agent of the inflammatory endothelial cancer, Kaposi's sarcoma, and a PB-regulating virus. The virus encodes kaposin B (KapB), which induces actin stress fibers (SFs) and cell spindling as well as PB disassembly. We now show that KapB-mediated PB disassembly requires actin rearrangements, RhoA effectors, and the mechanoresponsive transcription activator, YAP. Moreover, ectopic expression of active YAP or exposure of ECs to mechanical forces caused PB disassembly in the absence of KapB. We propose that the viral protein KapB activates a mechanoresponsive signaling axis and links changes in cell shape and cytoskeletal structures to enhanced inflammatory molecule expression using PB disassembly. Our work implies that cytoskeletal changes in other pathologies may similarly impact the inflammatory environment.