Explore the vast range of titles available.

Spatial transcriptomic technologies and spatially annotated single-cell RNA sequencing datasets provide unprecedented opportunities to dissect cell–cell communication (CCC). However, incorporation of the spatial information and complex biochemical processes required in the reconstruction of CCC remains a major challenge. Here, we present COMMOT (COMMunication analysis by Optimal Transport) to infer CCC in spatial transcriptomics, which accounts for the competition between different ligand and receptor species as well as spatial distances between cells. A collective optimal transport method is developed to handle complex molecular interactions and spatial constraints. Furthermore, we introduce downstream analysis tools to infer spatial signaling directionality and genes regulated by signaling using machine learning models. We apply COMMOT to simulation data and eight spatial datasets acquired with five different technologies to show its effectiveness and robustness in identifying spatial CCC in data with varying spatial resolutions and gene coverages. Finally, COMMOT identifies new CCCs during skin morphogenesis in a case study of human epidermal development. This work presents a computational framework, COMMOT, to spatially infer cell–cell communication from transcriptomics data based on a variant of optimal transport (OT).

Biomaterial wound dressings, such as hydrogels, interact with host cells to regulate tissue repair. This study investigates how crosslinking of gelatin-based hydrogels influences immune and stromal cell behavior and wound healing in female mice. We observe that softer, lightly crosslinked hydrogels promote greater cellular infiltration and result in smaller scars compared to stiffer, heavily crosslinked hydrogels. Using single-cell RNA sequencing, we further show that heavily crosslinked hydrogels increase inflammation and lead to the formation of a distinct macrophage subpopulation exhibiting signs of oxidative activity and cell fusion. Conversely, lightly crosslinked hydrogels are more readily taken up by macrophages and integrated within the tissue. The physical properties differentially affect macrophage and fibroblast interactions, with heavily crosslinked hydrogels promoting pro-fibrotic fibroblast activity that drives macrophage fusion through RANKL signaling. These findings suggest that tuning the physical properties of hydrogels can guide cellular responses and improve healing, offering insights for designing better biomaterials for wound treatment. Impaired wound healing that leads to scar remains a clinical challenge. Here, the authors study the effects of hydrogel crosslinking on cellular behavior in skin wounds and its effect on immune and stromal cell activity.

Cell-cell communication (CCC) occurs across different biological scales, ranging from interactions between large groups of cells to interactions between individual cells, forming a hierarchical structure. Globally, CCC may exist between clusters or only subgroups of a cluster with varying size, while locally, a group of cells as sender or receiver may exhibit distinct signaling properties. Current existing methods infer CCC from single-cell RNA-seq or Spatial Transcriptomics only between predefined cell groups, neglecting the existing hierarchical structure within CCC that are determined by signaling molecules, in particular, ligands and receptors. Here, we develop CrossChat, a novel computational framework designed to infer and analyze the hierarchical cell-cell communication structures using two complementary approaches: a global hierarchical structure using a multi-resolution clustering method, and multiple local hierarchical structures using a tree detection method. This framework provides a comprehensive approach to understand the hierarchical relationships within CCC that govern complex tissue functions. By applying our method to two nonspatial scRNA-seq datasets sampled from COVID-19 patients and mouse embryonic skin, and two spatial transcriptomics datasets generated from Stereo-seq of mouse embryo and 10x Visium of mouse wounded skin, we showcase CrossChat’s functionalities for analyzing both global and local hierarchical structures within cell-cell communication. Cell-cell communication (CCC) plays a critical role in biological processes across different scales. Here, authors develop CrossChat, a computational framework that identifies and analyses hierarchical CCC structures using single-cell RNA-seq and spatial transcriptomics data, offering insights into complex tissue functions.

The crypts of the intestinal epithelium house the stem cells that ensure the continual renewal of the epithelial cells that line the intestinal tract. Crypt number increases by a process called crypt fission, the division of a single crypt into two daughter crypts. Fission drives normal tissue growth and maintenance. Correspondingly, it becomes less frequent in adulthood. Importantly, fission is reactivated to drive adenoma growth. The mechanisms governing fission are poorly understood. However, only by knowing how normal fission operates can cancer-associated changes be elucidated. We studied normal fission in tissue in three dimensions using high-resolution imaging and used intestinal organoids to identify underlying mechanisms. We discovered that both the number and relative position of Paneth cells and Lgr5+ cells are important for fission. Furthermore, the higher stiffness and increased adhesion of Paneth cells are involved in determining the site of fission. Formation of a cluster of Lgr5+ cells between at least two Paneth-cell-rich domains establishes the site for the upward invagination that initiates fission.

Background Many cancers arise from mutations in cells within epithelial tissues. Mutations manifesting at the subcellular level influence the structure and function of the tissue resulting in cancer. Previous work has proposed how cell level properties can lead to mutant cell invasion, but has not incorporated detailed subcellular modelling Results We present a framework that allows the straightforward integration and simulation of SBML representations of subcellular dynamics within multiscale models of epithelial tissues. This allows us to investigate the effect of mutations in subcellular pathways on the migration of cells within the colorectal crypt. Using multiple models we find that mutations in APC, a key component in the Wnt signalling pathway, can bias neutral drift and can also cause downward invasion of mutant cells in the crypt. Conclusions Our framework allows us to investigate how subcellular mutations, i.e. knockouts and knockdowns, affect cell-level properties and the resultant migration of cells within epithelial tissues. In the context of the colorectal crypt, we see that mutations in APC can lead directly to mutant cell invasion.

From single-cell RNA-sequencing (scRNA-seq) and spatial transcriptomics (ST), one can extract high-dimensional gene expression patterns that can be described by intercellular communication networks or decoupled gene modules. These two descriptions of information flow are often assumed to occur independently. However, intercellular communication drives directed flows of information that are mediated by intracellular gene modules, in turn triggering outflows of other signals. Methodologies to describe such intercellular flows are lacking. We present FlowSig, a method that infers communication-driven intercellular flows from scRNA-seq or ST data using graphical causal modeling and conditional independence. We benchmark FlowSig using newly generated experimental cortical organoid data and synthetic data generated from mathematical modeling. We demonstrate FlowSig’s utility by applying it to various studies, showing that FlowSig can capture stimulation-induced changes to paracrine signaling in pancreatic islets, demonstrate shifts in intercellular flows due to increasing COVID-19 severity and reconstruct morphogen-driven activator–inhibitor patterns in mouse embryogenesis. Using single-cell and spatial transcriptomics data, FlowSig provides a unified signaling modeling framework by connecting intercellular communication mediated by ligand–receptor interactions and intracellular gene expression modules.

From mathematical models of growth to computer simulations of pigmentation, the study of shell formation has given rise to an abundant number of models, working at various scales. Yet, attempts to combine those models have remained sparse, due to the challenge of combining categorically different approaches. In this paper, we propose a framework to streamline the process of combining the molecular and tissue scales of shell formation. We choose these levels as a proxy to link the genotype level, which is better described by molecular models, and the phenotype level, which is better described by tissue-level mechanics. We also show how to connect observations on shell populations to the approach, resulting in collections of molecular parameters that may be associated with different populations of real shell specimens. The approach is as follows: we use a Quality-Diversity algorithm, a type of black-box optimization algorithm, to explore the range of concentration profiles emerging as solutions of a molecular model, and that define growth patterns for the mechanical model. At the same time, the mechanical model is simulated over a wide range of growth patterns, resulting in a variety of spine shapes. While time-consuming, these steps only need to be performed once and then function as look-up tables. Actual pictures of shell spines can then be matched against the list of existing spine shapes, yielding a potential growth pattern which, in turn, gives us matching molecular parameters. The framework is modular, such that models can be easily swapped without changing the overall working of the method. As a demonstration of the approach, we solve specific molecular and mechanical models, adapted from available theoretical studies on molluscan shells, and apply the multiscale framework to evaluate the characteristics of spines from three distinct populations of Turbo sazae .

Phenotypic variation among species is a product of evolutionary changes to developmental programs 1 , 2 . However, how these changes generate novel morphological traits remains largely unclear. Here we studied the genomic and developmental basis of the mammalian gliding membrane, or patagium—an adaptative trait that has repeatedly evolved in different lineages, including in closely related marsupial species. Through comparative genomic analysis of 15 marsupial genomes, both from gliding and non-gliding species, we find that the Emx2 locus experienced lineage-specific patterns of accelerated cis -regulatory evolution in gliding species. By combining epigenomics, transcriptomics and in-pouch marsupial transgenics, we show that Emx2 is a critical upstream regulator of patagium development. Moreover, we identify different cis -regulatory elements that may be responsible for driving increased Emx2 expression levels in gliding species. Lastly, using mouse functional experiments, we find evidence that Emx2 expression patterns in gliders may have been modified from a pre-existing program found in all mammals. Together, our results suggest that patagia repeatedly originated through a process of convergent genomic evolution, whereby regulation of Emx2 was altered by distinct cis -regulatory elements in independently evolved species. Thus, different regulatory elements targeting the same key developmental gene may constitute an effective strategy by which natural selection has harnessed regulatory evolution in marsupial genomes to generate phenotypic novelty. Patagia—the mammalian gliding membrane—repeatedly originated through a process of convergent genomic evolution, whereby the regulation of Emx2 was altered by distinct cis -regulatory elements in independently evolved species.

We consider mechanically-induced pattern formation within the framework of a growing, planar, elastic rod attached to an elastic foundation. Through a combination of weakly nonlinear analysis and numerical methods, we identify how the shape and type of buckling (super- or subcritical) depend on material parameters, and a complex phase-space of transition from super- to subcritical is uncovered. We then examine the effect of heterogeneity on buckling and post-buckling behaviour, in the context of a heterogeneous substrate adhesion, elastic stiffness, or growth. We show how the same functional form of heterogeneity in different properties is manifest in a vastly differing post-buckled shape. Finally, a fourth form of heterogeneity, an imperfect foundation, is incorporated and shown to have a more dramatic impact on the buckling instability, a difference that can be qualitatively understood via the weakly nonlinear analysis.

Crypt fission is an in vivo tissue deformation process that is involved in both intestinal homeostasis and colorectal tumourigenesis. Despite its importance, the mechanics underlying crypt fission are currently poorly understood. Recent experimental development of organoids, organ-like buds cultured from crypt stem cells in vitro, has shown promise in shedding light on crypt fission. Drawing inspiration from observations of organoid growth and fission in vivo, we develop a computational model of a deformable epithelial tissue layer. Results from in silico experiments show the stiffness of cells and the proportions of cell subpopulations affect the nature of deformation in the epithelial layer. In particular, we find that increasing the proportion of stiffer cells in the layer increases the likelihood of crypt fission occurring. This is in agreement with and helps explain recent experimental work.