The role of elasticity on adhesion and clustering of neurons on soft surfaces

The question of whether material stiffness enhances cell adhesion and clustering is still open to debate. Results from the literature are seemingly contradictory, with some reports illustrating that adhesion increases with surface stiffness and others suggesting that the performance of a system of cells is curbed by high values of elasticity. To address the role of elasticity as a regulator in neuronal cell adhesion and clustering, we investigated the topological characteristics of networks of neurons on polydimethylsiloxane (PDMS) surfaces - with values of elasticity (E) varying in the 0.55-2.65 MPa range. Results illustrate that, as elasticity increases, the number of neurons adhering on the surface decreases. Notably, the small-world coefficient - a topological measure of networks - also decreases. Numerical simulations and functional multi-calcium imaging experiments further indicated that the activity of neuronal cells on soft surfaces improves for decreasing E. Experimental findings are supported by a mathematical model, that explains adhesion and clustering of cells on soft materials as a function of few parameters - including the Young's modulus and roughness of the material. Overall, results indicate that - in the considered elasticity interval - increasing the compliance of a material improves adhesion, improves clustering, and enhances communication of neurons.

From pixels to connections: exploring in vitro neuron reconstruction software for network graph generation

Digital reconstruction has been instrumental in deciphering how in vitro neuron architecture shapes information flow. Emerging approaches reconstruct neural systems as networks with the aim of understanding their organization through graph theory. Computational tools dedicated to this objective build models of nodes and edges based on key cellular features such as somata, axons, and dendrites. Fully automatic implementations of these tools are readily available, but they may also be purpose-built from specialized algorithms in the form of multi-step pipelines. Here we review software tools informing the construction of network models, spanning from noise reduction and segmentation to full network reconstruction. The scope and core specifications of each tool are explicitly defined to assist bench scientists in selecting the most suitable option for their microscopy dataset. Existing tools provide a foundation for complete network reconstruction, however more progress is needed in establishing morphological bases for directed/weighted connectivity and in software validation.

Human skeletal muscle aging atlas

Skeletal muscle aging is a key contributor to age-related frailty and sarcopenia with substantial implications for global health. Here we profiled 90,902 single cells and 92,259 single nuclei from 17 donors to map the aging process in the adult human intercostal muscle, identifying cellular changes in each muscle compartment. We found that distinct subsets of muscle stem cells exhibit decreased ribosome biogenesis genes and increased CCL2 expression, causing different aging phenotypes. Our atlas also highlights an expansion of nuclei associated with the neuromuscular junction, which may reflect re-innervation, and outlines how the loss of fast-twitch myofibers is mitigated through regeneration and upregulation of fast-type markers in slow-twitch myofibers with age. Furthermore, we document the function of aging muscle microenvironment in immune cell attraction. Overall, we present a comprehensive human skeletal muscle aging resource ( https://www.muscleageingcellatlas.org/ ) together with an in-house mouse muscle atlas to study common features of muscle aging across species.

Integrative analysis of yeast colony growth

Yeast colonies are routinely grown on agar plates in everyday experimental settings to understand basic molecular processes, produce novel drugs, improve health, and so on. Standardized conditions ensure these colonies grow in a reproducible fashion, while in nature microbes are under a constantly changing environment. Here we combine the power of computational simulations and laboratory experiments to investigate the impact of non-standard environmental factors on colony growth. We present the developement and parameterization of a quantitative agent-based model for yeast colony growth to reproduce measurements on colony size and cell number in a colony at non-standard environmental conditions. Specifically, we establish experimental conditions that mimic the effects of humidity changes and nutrient gradients. Our results show how colony growth is affected by moisture changes, nutrient availability, and initial colony inoculation conditions. We show that initial colony spread, not initial cell number have higher impact on the final size and cell number of colonies. Parameters of the model were identified by fitting these experiments and the fitted model gives guidance to establish conditions which enable unlimited growth of yeast colonies.

Spatial transcriptomics reveals discrete tumour microenvironments and autocrine loops within ovarian cancer subclones

High-grade serous ovarian carcinoma (HGSOC) is genetically unstable and characterised by the presence of subclones with distinct genotypes. Intratumoural heterogeneity is linked to recurrence, chemotherapy resistance, and poor prognosis. Here, we use spatial transcriptomics to identify HGSOC subclones and study their association with infiltrating cell populations. Visium spatial transcriptomics reveals multiple tumour subclones with different copy number alterations present within individual tumour sections. These subclones differentially express various ligands and receptors and are predicted to differentially associate with different stromal and immune cell populations. In one sample, CosMx single molecule imaging reveals subclones differentially associating with immune cell populations, fibroblasts, and endothelial cells. Cell-to-cell communication analysis identifies subclone-specific signalling to stromal and immune cells and multiple subclone-specific autocrine loops. Our study highlights the high degree of subclonal heterogeneity in HGSOC and suggests that subclone-specific ligand and receptor expression patterns likely modulate how HGSOC cells interact with their local microenvironment.

Engineering pluripotent stem cells with synthetic biology for regenerative medicine

Pluripotent stem cells (PSCs), characterized by self-renewal and capacity of differentiating into three germ layers, are the programmable building blocks of life. PSC-derived cells and multicellular systems, particularly organoids, exhibit great potential for regenerative medicine. However, this field is still in its infancy, partly due to limited strategies to robustly and precisely control stem cell behaviors, which are tightly regulated by inner gene regulatory networks in response to stimuli from the extracellular environment. Synthetic receptors and genetic circuits are powerful tools to customize the cellular sense-and-response process, suggesting their underlying roles in precise control of cell fate decision and function reconstruction. Herein, we review the progress and challenges needed to be overcome in the fields of PSC-based cell therapy and multicellular system generation, respectively. Furthermore, we summarize several well-established synthetic biology tools and their applications in PSC engineering. Finally, we highlight the challenges and perspectives of harnessing synthetic biology to PSC engineering for regenerative medicine.

Longevity interventions modulate mechanotransduction and extracellular matrix homeostasis in C. elegans

Dysfunctional extracellular matrices (ECM) contribute to aging and disease. Repairing dysfunctional ECM could potentially prevent age-related pathologies. Interventions promoting longevity also impact ECM gene expression. However, the role of ECM composition changes in healthy aging remains unclear. Here we perform proteomics and in-vivo monitoring to systematically investigate ECM composition (matreotype) during aging in C. elegans revealing three distinct collagen dynamics. Longevity interventions slow age-related collagen stiffening and prolong the expression of collagens that are turned over. These prolonged collagen dynamics are mediated by a mechanical feedback loop of hemidesmosome-containing structures that span from the exoskeletal ECM through the hypodermis, basement membrane ECM, to the muscles, coupling mechanical forces to adjust ECM gene expression and longevity via the transcriptional co-activator YAP-1 across tissues. Our results provide in-vivo evidence that coordinated ECM remodeling through mechanotransduction is required and sufficient to promote longevity, offering potential avenues for interventions targeting ECM dynamics.

Multimodal spatiotemporal phenotyping of human retinal organoid development

Organoids generated from human pluripotent stem cells provide experimental systems to study development and disease, but quantitative measurements across different spatial scales and molecular modalities are lacking. In this study, we generated multiplexed protein maps over a retinal organoid time course and primary adult human retinal tissue. We developed a toolkit to visualize progenitor and neuron location, the spatial arrangements of extracellular and subcellular components and global patterning in each organoid and primary tissue. In addition, we generated a single-cell transcriptome and chromatin accessibility timecourse dataset and inferred a gene regulatory network underlying organoid development. We integrated genomic data with spatially segmented nuclei into a multimodal atlas to explore organoid patterning and retinal ganglion cell (RGC) spatial neighborhoods, highlighting pathways involved in RGC cell death and showing that mosaic genetic perturbations in retinal organoids provide insight into cell fate regulation.

Maintenance of appropriate size scaling of the C. elegans pharynx by YAP-1

Even slight imbalance between the growth rate of different organs can accumulate to a large deviation from their appropriate size during development. Here, we use live imaging of the pharynx of C. elegans to ask if and how organ size scaling nevertheless remains uniform among individuals. Growth trajectories of hundreds of individuals reveal that pharynxes grow by a near constant volume per larval stage that is independent of their initial size, such that undersized pharynxes catch-up in size during development. Tissue-specific depletion of RAGA-1, an activator of mTOR and growth, shows that maintaining correct pharynx-to-body size proportions involves a bi-directional coupling between pharynx size and body growth. In simulations, this coupling cannot be explained by limitation of food uptake alone, and genetic experiments reveal an involvement of the mechanotransducing transcriptional co-regulator yap-1. Our data suggests that mechanotransduction coordinates pharynx growth with other tissues, ensuring body plan uniformity among individuals.

PhysiBoSS 2.0: a sustainable integration of stochastic Boolean and agent-based modelling frameworks

In systems biology, mathematical models and simulations play a crucial role in understanding complex biological systems. Different modelling frameworks are employed depending on the nature and scales of the system under study. For instance, signalling and regulatory networks can be simulated using Boolean modelling, whereas multicellular systems can be studied using agent-based modelling. Herein, we present PhysiBoSS 2.0, a hybrid agent-based modelling framework that allows simulating signalling and regulatory networks within individual cell agents. PhysiBoSS 2.0 is a redesign and reimplementation of PhysiBoSS 1.0 and was conceived as an add-on that expands the PhysiCell functionalities by enabling the simulation of intracellular cell signalling using MaBoSS while keeping a decoupled, maintainable and model-agnostic design. PhysiBoSS 2.0 also expands the set of functionalities offered to the users, including custom models and cell specifications, mechanistic submodels of substrate internalisation and detailed control over simulation parameters. Together with PhysiBoSS 2.0, we introduce PCTK, a Python package developed for handling and processing simulation outputs, and generating summary plots and 3D renders. PhysiBoSS 2.0 allows studying the interplay between the microenvironment, the signalling pathways that control cellular processes and population dynamics, suitable for modelling cancer. We show different approaches for integrating Boolean networks into multi-scale simulations using strategies to study the drug effects and synergies in models of cancer cell lines and validate them using experimental data. PhysiBoSS 2.0 is open-source and publicly available on GitHub with several repositories of accompanying interoperable tools.

General, open-source vertex modeling in biological applications using Tissue Forge

Vertex models are a widespread approach for describing the biophysics and behaviors of multicellular systems, especially of epithelial tissues. Vertex models describe a wide variety of developmental scenarios and behaviors like cell rearrangement and tissue folding. Often, these models are implemented as single-use or closed-source software, which inhibits reproducibility and decreases accessibility for researchers with limited proficiency in software development and numerical methods. We developed a physics-based vertex model methodology in Tissue Forge, an open-source, particle-based modeling and simulation environment. Our methodology describes the properties and processes of vertex model objects on the basis of vertices, which allows integration of vertex modeling with the particle-based formalism of Tissue Forge, enabling an environment for developing mixed-method models of multicellular systems. Our methodology in Tissue Forge inherits all features provided by Tissue Forge, delivering open-source, extensible vertex modeling with interactive simulation, real-time simulation visualization and model sharing in the C, C++ and Python programming languages and a Jupyter Notebook. Demonstrations show a vertex model of cell sorting and a mixed-method model of cell migration combining vertex- and particle-based models. Our methodology provides accessible vertex modeling for a broad range of scientific disciplines, and we welcome community-developed contributions to our open-source software implementation.

Game-theoretical description of the go-or-grow dichotomy in tumor development for various settings and parameter constellations

A medically important feature of several types of tumors is their ability to "decide" between staying at a primary site in the body or leaving it and forming metastases. The present theoretical study aims to provide a better understanding of the ultimate reasons for this so-called "go-or-grow" dichotomy. To that end, we use game theory, which has proven to be useful in analyzing the competition between tumors and healthy tissues or among different tumor cells. We begin by determining the game types in the Basanta-Hatzikirou-Deutsch model, depending on the parameter values. Thereafter, we suggest and analyze five modified variants of the model. For example, in the basic model, the deadlock game, Prisoner's Dilemma, and hawk-dove game can occur. The modified versions lead to several additional game types, such as battle of the sexes, route-choice, and stag-hunt games. For some game types, all cells are predicted to stay on their original site ("grow phenotype"), while for other types, only a certain fraction stay and the other cells migrate away ("go phenotype"). If the nutrient supply at a distant site is high, all the cells are predicted to go. We discuss our predictions in terms of the pros and cons of caloric restriction and limitations of the supply of vitamins or methionine. Our results may help devise treatments to prevent metastasis.

Topological data analysis of spatial patterning in heterogeneous cell populations: clustering and sorting with varying cell-cell adhesion

Different cell types aggregate and sort into hierarchical architectures during the formation of animal tissues. The resulting spatial organization depends (in part) on the strength of adhesion of one cell type to itself relative to other cell types. However, automated and unsupervised classification of these multicellular spatial patterns remains challenging, particularly given their structural diversity and biological variability. Recent developments based on topological data analysis are intriguing to reveal similarities in tissue architecture, but these methods remain computationally expensive. In this article, we show that multicellular patterns organized from two interacting cell types can be efficiently represented through persistence images. Our optimized combination of dimensionality reduction via autoencoders, combined with hierarchical clustering, achieved high classification accuracy for simulations with constant cell numbers. We further demonstrate that persistence images can be normalized to improve classification for simulations with varying cell numbers due to proliferation. Finally, we systematically consider the importance of incorporating different topological features as well as information about each cell type to improve classification accuracy. We envision that topological machine learning based on persistence images will enable versatile and robust classification of complex tissue architectures that occur in development and disease.

Inference of long-range cell-cell force transmission from ECM remodeling fluctuations

Cells sense, manipulate and respond to their mechanical microenvironment in a plethora of physiological processes, yet the understanding of how cells transmit, receive and interpret environmental cues to communicate with distant cells is severely limited due to lack of tools to quantitatively infer the complex tangle of dynamic cell-cell interactions in complicated environments. We present a computational method to systematically infer and quantify long-range cell-cell force transmission through the extracellular matrix (cell-ECM-cell communication) by correlating ECM remodeling fluctuations in between communicating cells and demonstrating that these fluctuations contain sufficient information to define unique signatures that robustly distinguish between different pairs of communicating cells. We demonstrate our method with finite element simulations and live 3D imaging of fibroblasts and cancer cells embedded in fibrin gels. While previous studies relied on the formation of a visible fibrous 'band' extending between cells to inform on mechanical communication, our method detected mechanical propagation even in cases where visible bands never formed. We revealed that while contractility is required, band formation is not necessary, for cell-ECM-cell communication, and that mechanical signals propagate from one cell to another even upon massive reduction in their contractility. Our method sets the stage to measure the fundamental aspects of intercellular long-range mechanical communication in physiological contexts and may provide a new functional readout for high content 3D image-based screening. The ability to infer cell-ECM-cell communication using standard confocal microscopy holds the promise for wide use and democratizing the method.

Organization of the human intestine at single-cell resolution

The intestine is a complex organ that promotes digestion, extracts nutrients, participates in immune surveillance, maintains critical symbiotic relationships with microbiota and affects overall health1. The intesting has a length of over nine metres, along which there are differences in structure and function2. The localization of individual cell types, cell type development trajectories and detailed cell transcriptional programs probably drive these differences in function. Here, to better understand these differences, we evaluated the organization of single cells using multiplexed imaging and single-nucleus RNA and open chromatin assays across eight different intestinal sites from nine donors. Through systematic analyses, we find cell compositions that differ substantially across regions of the intestine and demonstrate the complexity of epithelial subtypes, and find that the same cell types are organized into distinct neighbourhoods and communities, highlighting distinct immunological niches that are present in the intestine. We also map gene regulatory differences in these cells that are suggestive of a regulatory differentiation cascade, and associate intestinal disease heritability with specific cell types. These results describe the complexity of the cell composition, regulation and organization for this organ, and serve as an important reference map for understanding human biology and disease.

A pharmacoproteomic landscape of organotypic intervention responses in Gram-negative sepsis

Sepsis is the major cause of mortality across intensive care units globally, yet details of accompanying pathological molecular events remain unclear. This knowledge gap has resulted in ineffective biomarker development and suboptimal treatment regimens to prevent and manage organ dysfunction/damage. Here, we used pharmacoproteomics to score time-dependent treatment impact in a murine Escherichia coli sepsis model after administering beta-lactam antibiotic meropenem (Mem) and/or the immunomodulatory glucocorticoid methylprednisolone (Gcc). Three distinct proteome response patterns were identified, which depended on the underlying proteotype for each organ. Gcc enhanced some positive proteome responses of Mem, including superior reduction of the inflammatory response in kidneys and partial restoration of sepsis-induced metabolic dysfunction. Mem introduced sepsis-independent perturbations in the mitochondrial proteome that Gcc counteracted. We provide a strategy for the quantitative and organotypic assessment of treatment effects of candidate therapies in relationship to dosing, timing, and potential synergistic intervention combinations during sepsis.

Nucleocytoplasmic transport of active HER2 causes fractional escape from the DCIS-like state

Activation of HER2/ErbB2 coincides with escape from ductal carcinoma in situ (DCIS) premalignancy and disrupts 3D organization of cultured breast-epithelial spheroids. The 3D phenotype is infrequent, however, and mechanisms for its incomplete penetrance have been elusive. Using inducible HER2/ErbB2-EGFR/ErbB1 heterodimers, we match phenotype penetrance to the frequency of co-occurring transcriptomic changes and uncover a reconfiguration in the karyopherin network regulating ErbB nucleocytoplasmic transport. Induction of the exportin CSE1L inhibits nuclear accumulation of ErbBs, whereas nuclear ErbBs silence the importin KPNA1 by inducing miR-205. When these negative feedbacks are incorporated into a validated systems model of nucleocytoplasmic transport, steady-state localization of ErbB cargo becomes ultrasensitive to initial CSE1L abundance. Erbb2-driven carcinomas with Cse1l deficiency outgrow less irregularly from mammary ducts, and NLS-attenuating mutants or variants of HER2 favor escape in 3D culture. We conclude here that adaptive nucleocytoplasmic relocalization of HER2 creates a systems-level molecular switch at the premalignant-to-malignant transition.

The application of pancreatic cancer organoids for novel drug discovery

INTRODUCTION

Pancreatic ductal adenocarcinoma presents with a dismal prognosis. Personalized therapy is urgently warranted to overcome the treatment limitations of the "one-size-fits-all" scheme. Organoids have emerged as fundamental novel tools to study tumor biology and heterogeneity, hence overcoming limitations of other model systems by better-reflecting tissue heterogeneity and recapitulating in-vivo processes. Besides their crucial role in basic research, they have evolved as tools for translational drug discovery and patient stratification.

AREAS COVERED

This review highlights the achievements of an organoid-based drug investigation and discovery. The authors present an overview of studies using organoids for drug testing. Further, they pinpoint studies correlating the in vitro prediction of organoids to the actual patient`s response. Furthermore, the authors describe novel model systems and take a thorough overlook of microfluidic chips, synthetic matrices, multicellular systems, bioprinting, and stem cell-derived pancreatic organoid systems.

EXPERT OPINION

Organoid systems promise great potential for future clinical applications. Indeed, they may be implemented into informed decision-making for guiding therapies. However, validation by randomized trials is mandatory. Additionally, organoids in combination with other cellular compartments may be exploited for drug discovery by studying niche-tumor interaction. Yet, several precautions must be kept in mind, such as standardization and reproducibility.

3D multicellular systems in disease modelling: From organoids to organ-on-chip

Cell-cell interactions underlay organ formation and function during homeostasis. Changes in communication between cells and their surrounding microenvironment are a feature of numerous human diseases, including metabolic disease and neurological disorders. In the past decade, cross-disciplinary research has been conducted to engineer novel synthetic multicellular organ systems in 3D, including organoids, assembloids, and organ-on-chip models. These model systems, composed of distinct cell types, satisfy the need for a better understanding of complex biological interactions and mechanisms underpinning diseases. In this review, we discuss the emerging field of building 3D multicellular systems and their application for modelling the cellular interactions at play in diseases. We report recent experimental and computational approaches for capturing cell-cell interactions as well as progress in bioengineering approaches for recapitulating these complexities ex vivo. Finally, we explore the value of developing such multicellular systems for modelling metabolic, intestinal, and neurological disorders as major examples of multisystemic diseases, we discuss the advantages and disadvantages of the different approaches and provide some recommendations for further advancing the field.

Jag1-Notch cis-interaction determines cell fate segregation in pancreatic development

The Notch ligands Jag1 and Dll1 guide differentiation of multipotent pancreatic progenitor cells (MPCs) into unipotent pro-acinar cells (PACs) and bipotent duct/endocrine progenitors (BPs). Ligand-mediated trans-activation of Notch receptors induces oscillating expression of the transcription factor Hes1, while ligand-receptor cis-interaction indirectly represses Hes1 activation. Despite Dll1 and Jag1 both displaying cis- and trans-interactions, the two mutants have different phenotypes for reasons not fully understood. Here, we present a mathematical model that recapitulates the spatiotemporal differentiation of MPCs into PACs and BPs. The model correctly captures cell fate changes in Notch pathway knockout mice and small molecule inhibitor studies, and a requirement for oscillatory Hes1 expression to maintain the multipotent state. Crucially, the model entails cell-autonomous attenuation of Notch signaling by Jag1-mediated cis-inhibition in MPC differentiation. The model sheds light on the underlying mechanisms, suggesting that cis-interaction is crucial for exiting the multipotent state, while trans-interaction is required for adopting the bipotent fate.

Slowest possible replicative life at frigid temperatures for yeast

Determining whether life can progress arbitrarily slowly may reveal fundamental barriers to staying out of thermal equilibrium for living systems. By monitoring budding yeast's slowed-down life at frigid temperatures and with modeling, we establish that Reactive Oxygen Species (ROS) and a global gene-expression speed quantitatively determine yeast's pace of life and impose temperature-dependent speed limits - shortest and longest possible cell-doubling times. Increasing cells' ROS concentration increases their doubling time by elongating the cell-growth (G1-phase) duration that precedes the cell-replication (S-G2-M) phase. Gene-expression speed constrains cells' ROS-reducing rate and sets the shortest possible doubling-time. To replicate, cells require below-threshold concentrations of ROS. Thus, cells with sufficiently abundant ROS remain in G1, become unsustainably large and, consequently, burst. Therefore, at a given temperature, yeast's replicative life cannot progress arbitrarily slowly and cells with the lowest ROS-levels replicate most rapidly. Fundamental barriers may constrain the thermal slowing of other organisms' lives.

A lineage tree-based hidden Markov model quantifies cellular heterogeneity and plasticity

Individual cells can assume a variety of molecular and phenotypic states and recent studies indicate that cells can rapidly adapt in response to therapeutic stress. Such phenotypic plasticity may confer resistance, but also presents opportunities to identify molecular programs that could be targeted for therapeutic benefit. Approaches to quantify tumor-drug responses typically focus on snapshot, population-level measurements. While informative, these methods lack lineage and temporal information, which are particularly critical for understanding dynamic processes such as cell state switching. As new technologies have become available to measure lineage relationships, modeling approaches will be needed to identify the forms of cell-to-cell heterogeneity present in these data. Here we apply a lineage tree-based adaptation of a hidden Markov model that employs single cell lineages as input to learn the characteristic patterns of phenotypic heterogeneity and state transitions. In benchmarking studies, we demonstrated that the model successfully classifies cells within experimentally-tractable dataset sizes. As an application, we analyzed experimental measurements in cancer and non-cancer cell populations under various treatments. We find evidence of multiple phenotypically distinct states, with considerable heterogeneity and unique drug responses. In total, this framework allows for the flexible modeling of single cell heterogeneity across lineages to quantify, understand, and control cell state switching.

Mapping the epithelial-immune cell interactome upon infection in the gut and the upper airways

Increasing evidence points towards the key role of the epithelium in the systemic and over-activated immune response to viral infection, including SARS-CoV-2 infection. Yet, how viral infection alters epithelial-immune cell interactions regulating inflammatory responses, is not well known. Available experimental approaches are insufficient to properly analyse this complex system, and computational predictions and targeted data integration are needed as an alternative approach. In this work, we propose an integrated computational biology framework that models how infection alters intracellular signalling of epithelial cells and how this change impacts the systemic immune response through modified interactions between epithelial cells and local immune cell populations. As a proof-of-concept, we focused on the role of intestinal and upper-airway epithelial infection. To characterise the modified epithelial-immune interactome, we integrated intra- and intercellular networks with single-cell RNA-seq data from SARS-CoV-2 infected human ileal and colonic organoids as well as from infected airway ciliated epithelial cells. This integrated methodology has proven useful to point out specific epithelial-immune interactions driving inflammation during disease response, and propose relevant molecular targets to guide focused experimental analysis.

Single-cell analysis identifies a key role for Hhip in murine coronal suture development

Craniofacial development depends on formation and maintenance of sutures between bones of the skull. In sutures, growth occurs at osteogenic fronts along the edge of each bone, and suture mesenchyme separates adjacent bones. Here, we perform single-cell RNA-seq analysis of the embryonic, wild type murine coronal suture to define its population structure. Seven populations at E16.5 and nine at E18.5 comprise the suture mesenchyme, osteogenic cells, and associated populations. Expression of Hhip, an inhibitor of hedgehog signaling, marks a mesenchymal population distinct from those of other neurocranial sutures. Tracing of the neonatal Hhip-expressing population shows that descendant cells persist in the coronal suture and contribute to calvarial bone growth. In Hhip^-/- coronal sutures at E18.5, the osteogenic fronts are closely apposed and the suture mesenchyme is depleted with increased hedgehog signaling compared to those of the wild type. Collectively, these data demonstrate that Hhip is required for normal coronal suture development.

Evidence-Based Network Modelling to Simulate Nucleus Pulposus Multicellular Activity in Different Nutritional and Pro-Inflammatory Environments

Initiation of intervertebral disc degeneration is thought to be biologically driven. This reflects a process, where biochemical and mechanical stimuli affect cell activity (CA) that compromise the tissue strength over time. Experimental research enhanced our understanding about the effect of such stimuli on different CA, such as protein synthesis or mRNA expression. However, it is still unclear how cells respond to their native environment that consists of a “cocktail” of different stimuli that might locally vary. This work presents an interdisciplinary approach of experimental and in silico research to approximate Nucleus Pulposus CA within multifactorial biochemical environments. Thereby, the biochemical key stimuli glucose, pH, and the proinflammatory cytokines TNF-α and IL1β were considered that were experimentally shown to critically affect CA. To this end, a Nucleus Pulposus multicellular system was modelled. It integrated experimental findings from in vitro studies of human or bovine Nucleus Pulposus cells, to relate the individual effects of targeted stimuli to alterations in CA. Unknown stimulus-CA relationships were obtained through own experimental 3D cultures of bovine Nucleus Pulposus cells in alginate beads. Translation of experimental findings into suitable parameters for network modelling approaches was achieved thanks to a new numerical approach to estimate the individual sensitivity of a CA to each stimulus type. Hence, the effect of each stimulus type on a specific CA was assessed and integrated to approximate a multifactorial stimulus environment. Tackled CA were the mRNA expressions of Aggrecan, Collagen types I & II, MMP3, and ADAMTS4. CA was assessed for four different proinflammatory cell states; non-inflamed and inflamed for IL1β, TNF-α or both IL1β&TNF-α. Inflamed cell clusters were eventually predicted in a multicellular 3D agent-based model. Experimental results showed that glucose had no significant impact on proinflammatory cytokine or ADAMTS4 mRNA expression, whereas TNF-α caused a significant catabolic shift in most explored CA. In silico results showed that the presented methodology to estimate the sensitivity of a CA to a stimulus type importantly improved qualitative model predictions. However, more stimuli and/or further experimental knowledge need to be integrated, especially regarding predictions about the possible progression of inflammatory environments under adverse nutritional conditions. Tackling the multicellular level is a new and promising approach to estimate manifold responses of intervertebral disc cells. Such a top-down high-level network modelling approach allows to obtain information about relevant stimulus environments for a specific CA and could be shown to be suitable to tackle complex biological systems, including different proinflammatory cell states. The development of this methodology required a close interaction with experimental research. Thereby, specific experimental needs were derived from systematic in silico approaches and obtained results were directly used to enhance model predictions, which reflects a novelty in this research field. Eventually, the presented methodology provides modelling solutions suitable for multiscale approaches to contribute to a better understanding about dynamics over multiple spatial scales. Future work should focus on an amplification of the stimulus environment by integrating more key relevant stimuli, such as mechanical loading parameters, in order to better approximate native physiological environments.

Live cell tagging tracking and isolation for spatial transcriptomics using photoactivatable cell dyes

A cell's phenotype and function are influenced by dynamic interactions with its microenvironment. To examine cellular spatiotemporal activity, we developed SPACECAT-Spatially PhotoActivatable Color Encoded Cell Address Tags-to annotate, track, and isolate cells while preserving viability. In SPACECAT, samples are stained with photocaged fluorescent molecules, and cells are labeled by uncaging those molecules with user-patterned near-UV light. SPACECAT offers single-cell precision and temporal stability across diverse cell and tissue types. Illustratively, we target crypt-like regions in patient-derived intestinal organoids to enrich for stem-like and actively mitotic cells, matching literature expectations. Moreover, we apply SPACECAT to ex vivo tissue sections from four healthy organs and an autochthonous lung tumor model. Lastly, we provide a computational framework to identify spatially-biased transcriptome patterns and enriched phenotypes. This minimally perturbative and broadly applicable method links cellular spatiotemporal and/or behavioral phenotypes with diverse downstream assays, enabling insights into the connections between tissue microenvironments and (dys)function.

GMP compliant isolation of mucosal epithelial cells and fibroblasts from biopsy samples for clinical tissue engineering

Engineered epithelial cell sheets for clinical replacement of non-functional upper aerodigestive tract mucosa are regulated as medicinal products and should be manufactured to the standards of good manufacturing practice (GMP). The current gold standard for growth of epithelial cells for research utilises growth arrested murine 3T3 J2 feeder layers, which are not available for use as a GMP compliant raw material. Using porcine mucosal tissue, we demonstrate a new method for obtaining and growing non-keratinised squamous epithelial cells and fibroblast cells from a single biopsy, replacing the 3T3 J2 with a growth arrested primary fibroblast feeder layer and using pooled Human Platelet lysate (HPL) as the media serum supplement to replace foetal bovine serum (FBS). The initial isolation of the cells was semi-automated using an Octodissociator and the resultant cell suspension cryopreservation for future use. When compared to the gold standard of 3T3 J2 and FBS containing medium there was no reduction in growth, viability, stem cell population or ability to differentiate to mature epithelial cells. Furthermore, this method was replicated with Human buccal tissue, providing cells of sufficient quality and number to create a tissue engineered sheet.

Optogenetic current in myofibroblasts acutely alters electrophysiology and conduction of co-cultured cardiomyocytes

Interactions between cardiac myofibroblasts and myocytes may slow conduction and generate spontaneous beating in fibrosis, increasing the chance of life-threatening arrhythmia. While co-culture studies have shown that myofibroblasts can affect cardiomyocyte electrophysiology in vitro, the extent of myofibroblast-myocyte electrical conductance in a syncytium is unknown. In this neonatal rat study, cardiac myofibroblasts were transduced with Channelrhodopsin-2, which allowed acute and selective increase of myofibroblast current, and plated on top of cardiomyocytes. Optical mapping revealed significantly decreased conduction velocity (- 27 ± 6%, p < 10-3), upstroke rate (- 13 ± 4%, p = 0.002), and action potential duration (- 14 ± 7%, p = 0.004) in co-cultures when 0.017 mW/mm2 light was applied, as well as focal spontaneous beating in 6/7 samples and a decreased cycle length (- 36 ± 18%, p = 0.002) at 0.057 mW/mm2 light. In silico modeling of the experiments reproduced the experimental findings and suggested the light levels used in experiments produced excess current similar in magnitude to endogenous myofibroblast current. Fitting the model to experimental data predicted a tissue-level electrical conductance across the 3-D interface between myofibroblasts and cardiomyocytes of ~ 5 nS/cardiomyocyte, and showed how increased myofibroblast-myocyte conductance, increased myofibroblast/myocyte capacitance ratio, and increased myofibroblast current, which occur in fibrosis, can work in tandem to produce pro-arrhythmic increases in conduction and spontaneous beating.

Automated design of synthetic microbial communities

Microbial species rarely exist in isolation. In naturally occurring microbial systems there is strong evidence for a positive relationship between species diversity and productivity of communities. The pervasiveness of these communities in nature highlights possible advantages for genetically engineered strains to exist in cocultures as well. Building synthetic microbial communities allows us to create distributed systems that mitigate issues often found in engineering a monoculture, especially as functional complexity increases. Here, we demonstrate a methodology for designing robust synthetic communities that include competition for nutrients, and use quorum sensing to control amensal bacteriocin interactions in a chemostat environment. We computationally explore all two- and three- strain systems, using Bayesian methods to perform model selection, and identify the most robust candidates for producing stable steady state communities. Our findings highlight important interaction motifs that provide stability, and identify requirements for selecting genetic parts and further tuning the community composition.

A dynamic multi-tissue model to study human metabolism

Metabolic modeling enables the study of human metabolism in healthy and in diseased conditions, e.g., the prediction of new drug targets and biomarkers for metabolic diseases. To accurately describe blood and urine metabolite dynamics, the integration of multiple metabolically active tissues is necessary. We developed a dynamic multi-tissue model, which recapitulates key properties of human metabolism at the molecular and physiological level based on the integration of transcriptomics data. It enables the simulation of the dynamics of intra-cellular and extra-cellular metabolites at the genome scale. The predictive capacity of the model is shown through the accurate simulation of different healthy conditions (i.e., during fasting, while consuming meals or during exercise), and the prediction of biomarkers for a set of Inborn Errors of Metabolism with a precision of 83%. This novel approach is useful to prioritize new biomarkers for many metabolic diseases, as well as for the integration of various types of personal omics data, towards the personalized analysis of blood and urine metabolites.

Sound-induced morphogenesis of multicellular systems for rapid orchestration of vascular networks

Morphogenesis, a complex process, ubiquitous in developmental biology and many pathologies, is based on self-patterning of cells. Spatial patterns of cells, organoids, or inorganic particles can be forced on demand using acoustic surface standing waves, such as the Faraday waves. This technology allows tuning of parameters (sound frequency, amplitude, chamber shape) under contactless, fast and mild culture conditions, for morphologically relevant tissue generation. We call this method Sound Induced Morphogenesis (SIM). In this work, we use SIM to achieve tight control over patterning of endothelial cells and mesenchymal stem cells densities within a hydrogel, with the endpoint formation of vascular structures. Here, we first parameterize our system to produce enhanced cell density gradients. Second, we allow for vasculogenesis after SIM patterning control and compare our controlled technology against state-of-the-art microfluidic culture systems, the latter characteristic of pure self-organized patterning and uniform initial density. Our sound-induced cell density patterning and subsequent vasculogenesis requires less cells than the microfluidic chamber. We advocate for the use of SIM for rapid, mild, and reproducible morphogenesis induction and further explorations in the regenerative medicine and cell therapy fields.

Capturing Multicellular System Designs Using Synthetic Biology Open Language (SBOL)

Synthetic biology aims to develop novel biological systems and increase their reproducibility using engineering principles such as standardization and modularization. It is important that these systems can be represented and shared in a standard way to ensure they can be easily understood, reproduced, and utilized by other researchers. The Synthetic Biology Open Language (SBOL) is a data standard for sharing biological designs and information about their implementation and characterization. Previously, this standard has only been used to represent designs in systems where the same design is implemented in every cell; however, there is also much interest in multicellular systems, in which designs involve a mixture of different types of cells with differing genotype and phenotype. Here, we show how the SBOL standard can be used to represent multicellular systems, and, hence, how researchers can better share designs with the community and reliably document intended system functionality.

From Microbial Communities to Distributed Computing Systems

A distributed biological system can be defined as a system whose components are located in different subpopulations, which communicate and coordinate their actions through interpopulation messages and interactions. We see that distributed systems are pervasive in nature, performing computation across all scales, from microbial communities to a flock of birds. We often observe that information processing within communities exhibits a complexity far greater than any single organism. Synthetic biology is an area of research which aims to design and build synthetic biological machines from biological parts to perform a defined function, in a manner similar to the engineering disciplines. However, the field has reached a bottleneck in the complexity of the genetic networks that we can implement using monocultures, facing constraints from metabolic burden and genetic interference. This makes building distributed biological systems an attractive prospect for synthetic biology that would alleviate these constraints and allow us to expand the applications of our systems into areas including complex biosensing and diagnostic tools, bioprocess control and the monitoring of industrial processes. In this review we will discuss the fundamental limitations we face when engineering functionality with a monoculture, and the key areas where distributed systems can provide an advantage. We cite evidence from natural systems that support arguments in favor of distributed systems to overcome the limitations of monocultures. Following this we conduct a comprehensive overview of the synthetic communities that have been built to date, and the components that have been used. The potential computational capabilities of communities are discussed, along with some of the applications that these will be useful for. We discuss some of the challenges with building co-cultures, including the problem of competitive exclusion and maintenance of desired community composition. Finally, we assess computational frameworks currently available to aide in the design of microbial communities and identify areas where we lack the necessary tools.

Cellular Dialogues: Cell-Cell Communication through Diffusible Molecules Yields Dynamic Spatial Patterns

Cells form spatial patterns by coordinating their gene expressions. How a group of mesoscopic numbers (hundreds to thousands) of cells, without pre-existing morphogen gradients and spatial organization, self-organizes spatial patterns remains poorly understood. Of particular importance are dynamic spatial patterns such as spiral waves that perpetually move and transmit information. We developed an open-source software for simulating a field of cells that communicate by secreting any number of molecules. With this software and a theory, we identified all possible "cellular dialogues"-ways of communicating with two diffusing molecules-that yield diverse dynamic spatial patterns. These patterns emerge despite widely varying responses of cells to the molecules, gene-expression noise, spatial arrangements, and cell movements. A three-stage, "order-fluctuate-settle" process forms dynamic spatial patterns: cells form long-lived whirlpools of wavelets that, following erratic dynamics, settle into a dynamic spatial pattern. Our work helps in identifying gene-regulatory networks that underlie dynamic pattern formations.

A cellular automaton model for spheroid response to radiation and hyperthermia treatments

Thermo-radiosensitisation is a promising approach for treatment of radio-resistant tumours such as those containing hypoxic subregions. Response prediction and treatment planning should account for tumour response heterogeneity, e.g. due to microenvironmental factors, and quantification of the biological effects induced. 3D tumour spheroids provide a physiological in vitro model of tumour response and a systems oncology framework for simulating spheroid response to radiation and hyperthermia is presented. Using a cellular automaton model, 3D oxygen diffusion, delivery of radiation and/or hyperthermia were simulated for many ([Formula: see text]) individual cells forming a spheroid. The iterative oxygen diffusion model was compared to an analytical oxygenation model and simulations were calibrated and validated against experimental data for irradiated (0-10 Gy) and/or heated (0-240 CEM43) HCT116 spheroids. Despite comparable clonogenic survival, spheroid growth differed significantly following radiation or hyperthermia. This dynamic response was described well by the simulation ([Formula: see text] > 0.85). Heat-induced cell death was implemented as a fast, proliferation-independent process, allowing reoxygenation and repopulation, whereas radiation was modelled as proliferation-dependent mitotic catastrophe. This framework stands out both through its experimental validation and its novel ability to predict spheroid response to multimodality treatment. It provides a good description of response where biological dose-weighting based on clonogenic survival alone was insufficient.

Chemotactic behaviour of Escherichia coli at high cell density

At high cell density, swimming bacteria exhibit collective motility patterns, self-organized through physical interactions of a however still debated nature. Although high-density behaviours are frequent in natural situations, it remained unknown how collective motion affects chemotaxis, the main physiological function of motility, which enables bacteria to follow environmental gradients in their habitats. Here, we systematically investigate this question in the model organism Escherichia coli, varying cell density, cell length, and suspension confinement. The characteristics of the collective motion indicate that hydrodynamic interactions between swimmers made the primary contribution to its emergence. We observe that the chemotactic drift is moderately enhanced at intermediate cell densities, peaks, and is then strongly suppressed at higher densities. Numerical simulations reveal that this suppression occurs because the collective motion disturbs the choreography necessary for chemotactic sensing. We suggest that this physical hindrance imposes a fundamental constraint on high-density behaviours of motile bacteria, including swarming and the formation of multicellular aggregates and biofilms.

A multi-modal data resource for investigating topographic heterogeneity in patient-derived xenograft tumors

Patient-derived xenografts (PDXs) are an essential pre-clinical resource for investigating tumor biology. However, cellular heterogeneity within and across PDX tumors can strongly impact the interpretation of PDX studies. Here, we generated a multi-modal, large-scale dataset to investigate PDX heterogeneity in metastatic colorectal cancer (CRC) across tumor models, spatial scales and genomic, transcriptomic, proteomic and imaging assay modalities. To showcase this dataset, we present analysis to assess sources of PDX variation, including anatomical orientation within the implanted tumor, mouse contribution, and differences between replicate PDX tumors. A unique aspect of our dataset is deep characterization of intra-tumor heterogeneity via immunofluorescence imaging, which enables investigation of variation across multiple spatial scales, from subcellular to whole tumor levels. Our study provides a benchmark data resource to investigate PDX models of metastatic CRC and serves as a template for future, quantitative investigations of spatial heterogeneity within and across PDX tumor models.

In vivo measurements reveal a single 5'-intron is sufficient to increase protein expression level in Caenorhabditis elegans

Introns can increase gene expression levels using a variety of mechanisms collectively referred to as Intron Mediated Enhancement (IME). IME has been measured in cell culture and plant models by quantifying expression of intronless and intron-bearing reporter genes in vitro. We developed hardware and software to implement microfluidic chip-based gene expression quantification in vivo. We altered position, number and sequence of introns in reporter genes controlled by the hsp-90 promoter. Consistent with plant and mammalian studies, we determined a single, natural or synthetic, 5'-intron is sufficient for the full IME effect conferred by three synthetic introns, while a 3'-intron is not. We found coding sequence can affect IME; the same three synthetic introns that increase mcherry protein concentration by approximately 50%, increase mEGFP by 80%. We determined IME effect size is not greatly affected by the stronger vit-2 promoter. Our microfluidic imaging approach should facilitate screens for factors affecting IME and other intron-dependent processes.

Programming Morphogenesis through Systems and Synthetic Biology

Mammalian tissue development is an intricate, spatiotemporal process of self-organization that emerges from gene regulatory networks of differentiating stem cells. A major goal in stem cell biology is to gain a sufficient understanding of gene regulatory networks and cell-cell interactions to enable the reliable and robust engineering of morphogenesis. Here, we review advances in synthetic biology, single cell genomics, and multiscale modeling, which, when synthesized, provide a framework to achieve the ambitious goal of programming morphogenesis in complex tissues and organoids.

Cell-Cell Communication in Yeast Using Auxin Biosynthesis and Auxin Responsive CRISPR Transcription Factors

An engineering framework for synthetic multicellular systems requires a programmable means of cell-cell communication. Such a communication system would enable complex behaviors, such as pattern formation, division of labor in synthetic microbial communities, and improved modularity in synthetic circuits. However, it remains challenging to build synthetic cellular communication systems in eukaryotes due to a lack of molecular modules that are orthogonal to the host machinery, easy to reconfigure, and scalable. Here, we present a novel cell-to-cell communication system in Saccharomyces cerevisiae (yeast) based on CRISPR transcription factors and the plant hormone auxin that exhibits several of these features. Specifically, we engineered a sender strain of yeast that converts indole-3-acetamide (IAM) into auxin via the enzyme iaaH from Agrobacterium tumefaciens. To sense auxin and regulate transcription in a receiver strain, we engineered a reconfigurable library of auxin-degradable CRISPR transcription factors (ADCTFs). Auxin-induced degradation is achieved through fusion of an auxin-sensitive degron (from IAA corepressors) to the CRISPR TF and coexpression with an auxin F-box protein. Mirroring the tunability of auxin perception in plants, our family of ADCTFs exhibits a broad range of auxin sensitivities. We characterized the kinetics and steady-state behavior of the sender and receiver independently as well as in cocultures where both cell types were exposed to IAM. In the presence of IAM, auxin is produced by the sender cell and triggers deactivation of reporter expression in the receiver cell. The result is an orthogonal, rewireable, tunable, and, arguably, scalable cell-cell communication system for yeast and other eukaryotic cells.

A system for modelling cell-cell interactions during plant morphogenesis

BACKGROUND AND AIMS

During the development of multicellular organisms, cells are capable of interacting with each other through a range of biological and physical mechanisms. A description of these networks of cell-cell interactions is essential for an understanding of how cellular activity is co-ordinated in regionalized functional entities such as tissues or organs. The difficulty of experimenting on living tissues has been a major limitation to describing such systems, and computer modelling appears particularly helpful to characterize the behaviour of multicellular systems. The experimental difficulties inherent to the multitude of parallel interactions that underlie cellular morphogenesis have led to the need for computer models.

METHODS

A new generic model of plant cellular morphogenesis is described that expresses interactions amongst cellular entities explicitly: the plant is described as a multi-scale structure, and interactions between distinct entities is established through a topological neighbourhood. Tissues are represented as 2D biphasic systems where the cell wall responds to turgor pressure through a viscous yielding of the cell wall.

KEY RESULTS

This principle was used in the development of the CellModeller software, a generic tool dedicated to the analysis and modelling of plant morphogenesis. The system was applied to three contrasting study cases illustrating genetic, hormonal and mechanical factors involved in plant morphogenesis.

CONCLUSIONS

Plant morphogenesis is fundamentally a cellular process and the CellModeller software, through its underlying generic model, provides an advanced research tool to analyse coupled physical and biological morphogenetic mechanisms.