Morpheus: a user-friendly modeling environment for multiscale and multicellular systems biology

Understanding Insulin in the Age of Precision Medicine and Big Data: Under-Explored Nature of Genomics

Insulin is amongst the human genome's most well-studied genes/proteins due to its connection to metabolic health. Within this article, we review literature and data to build a knowledge base of Insulin (INS) genetics that influence transcription, transcript processing, translation, hormone maturation, secretion, receptor binding, and metabolism while highlighting the future needs of insulin research. The INS gene region has 2076 unique variants from population genetics. Several variants are found near the transcriptional start site, enhancers, and following the INS transcripts that might influence the readthrough fusion transcript INS-IGF2. This INS-IGF2 transcript splice site was confirmed within hundreds of pancreatic RNAseq samples, lacks drift based on human genome sequencing, and has possible elevated expression due to viral regulation within the liver. Moreover, a rare, poorly characterized African population-enriched variant of INS-IGF2 results in a loss of the stop codon. INS transcript UTR variants rs689 and rs3842753, associated with type 1 diabetes, are found in many pancreatic RNAseq datasets with an elevation of the 3'UTR alternatively spliced INS transcript. Finally, by combining literature, evolutionary profiling, and structural biology, we map rare missense variants that influence preproinsulin translation, proinsulin processing, dimer/hexamer secretory storage, receptor activation, and C-peptide detection for quasi-insulin blood measurements.

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.

Single-cell and bulk RNA sequencing reveal ligands and receptors associated with worse overall survival in serous ovarian cancer

Background

Serous ovarian carcinoma is the most frequent histological subgroup of ovarian cancer and the leading cause of death among gynecologic tumors. The tumor microenvironment and cancer-associated fibroblasts (CAFs) have a critical role in the origin and progression of cancer. We comprehensively characterized the crosstalk between CAFs and ovarian cancer cells from malignant fluids to identify specific ligands and receptors mediating intercellular communications and disrupted pathways related to prognosis and therapy response.

Methods

Malignant fluids of serous ovarian cancer, including tumor-derived organoids, CAFs-enriched (eCAFs), and malignant effusion cells (no cultured) paired with normal ovarian tissues, were explored by RNA-sequencing. These data were integrated with single-cell RNA-sequencing data of ascites from ovarian cancer patients. The most relevant ligand and receptor interactions were used to identify differentially expressed genes with prognostic values in ovarian cancer.

Results

CAF ligands and epithelial cancer cell receptors were enriched for PI3K-AKT, focal adhesion, and epithelial-mesenchymal transition signaling pathways. Collagens, MIF, MDK, APP, and laminin were detected as the most significant signaling, and the top ligand-receptor interactions THBS2/THBS3 (CAFs)—CD47 (cancer cells), MDK (CAFs)—NCL/SDC2/SDC4 (cancer cells) as potential therapeutic targets. Interestingly, 34 genes encoding receptors and ligands of the PI3K pathway were associated with the outcome, response to treatment, and overall survival in ovarian cancer. Up-regulated genes from this list consistently predicted a worse overall survival (hazard ratio > 1.0 and log-rank P < 0.05) in two independent validation cohorts.

Conclusions

This study describes critical signaling pathways, ligands, and receptors involved in the communication between CAFs and cancer cells that have prognostic and therapeutic significance in ovarian cancer.

AIM2 as a putative target in acute kidney graft rejection

Acute rejection (AR) is a process triggered via the recognition of grafted organ-derived antigens by the immune system, which could present as a life-threatening condition. In the context of a kidney transplant, despite improvement with immunosuppressive therapies, AR maintains a significant incidence of 10%, and currently available drugs generally act in similar and canonical pathways of lymphocyte activation. This prompted the research for different approaches to identify potential novel targets that could improve therapeutic interventions. Here, we conducted a transcriptome analysis comparing groups of acute rejection (including T cell-mediated rejection and antibody-mediated rejection) to stable grafts that included differentially expressed genes, transcription factor and kinase enrichment, and Gene Set Enrichment Analysis. These analyses revealed inflammasome enhancement in rejected grafts and AIM2 as a potential component linked to acute rejection, presenting a positive correlation to T-cell activation and a negative correlation to oxidative phosphorylation metabolism. Also, the AIM2 expression showed a global accuracy in discerning acute rejection grafts (area under the curve (AUC) = 0.755 and 0.894, p &lt; 0.0001), and meta-analysis comprising different studies indicated a considerable enhancement of AIM2 in rejection (standardized mean difference (SMD) = 1.45, [CI 95%, 1.18 to 1.71]), especially for T cell-mediated rejection (TCMR) (SMD = 2.01, [CI 95%, 1.58 to 2.45]). These findings could guide future studies of AIM2 as either an adjuvant target for immunosuppression or a potential biomarker for acute rejection and graft survival.

In Silico Modeling Study of Curcumin Diffusion and Cellular Growth

Curcumin can enhance cutaneous wound healing by improving fibroblast proliferation. However, its therapeutic properties are dose-dependent: high concentrations produce cytotoxic effects, whereas low concentrations benefit cell proliferation. Similarly, the type of administration and its moderation are key aspects, as an erroneous distribution may result in null or noxious activity to the organism. In silico models for curcumin diffusion work as predictive tools for evaluating curcumin’s cytotoxic effects and establishing therapeutic windows. A 2D fibroblast culture growth model was created based on a model developed by Gérard and Goldbeter. Similarly, a curcumin diffusion model was developed by adjusting experimental release values obtained from Aguilar-Rabiela et al. and fitted to Korsmeyer–Peppas and Peleg’s hyperbolic models. The release of six key curcumin concentrations was achieved. Both models were integrated using Morpheus software, and a scratch-wound assay simulated curcumin’s dose-dependent effects on wound healing. The most beneficial effect was achieved at 0.25 μM, which exhibited the lowest cell-division period, the highest confluence (~60% for both release models, 447 initial cells), and the highest final cell population. The least beneficial effect was found at 20 μM, which inhibited cell division and achieved the lowest confluence (~34.30% for both release models, 447 initial cells). Confluence was shown to decrease as curcumin concentration increased, since higher concentrations of curcumin have inhibitory and cytotoxic effects.

Evolutionary origin of vertebrate OCT4/POU5 functions in supporting pluripotency

The support of pluripotent cells over time is an essential feature of development. In eutherian embryos, pluripotency is maintained from naïve states in peri-implantation to primed pluripotency at gastrulation. To understand how these states emerged, we reconstruct the evolutionary trajectory of the Pou5 gene family, which contains the central pluripotency factor OCT4. By coupling evolutionary sequence analysis with functional studies in mouse embryonic stem cells, we find that the ability of POU5 proteins to support pluripotency originated in the gnathostome lineage, prior to the generation of two paralogues, Pou5f1 and Pou5f3 via gene duplication. In osteichthyans, retaining both genes, the paralogues differ in their support of naïve and primed pluripotency. The specialization of these duplicates enables the diversification of function in self-renewal and differentiation. By integrating sequence evolution, cell phenotypes, developmental contexts and structural modelling, we pinpoint OCT4 regions sufficient for naïve pluripotency and describe their adaptation over evolutionary time.

Melatonin Regulates the Daily Levels of Plasma Amino Acids, Acylcarnitines, Biogenic Amines, Sphingomyelins, and Hexoses in a Xenograft Model of Triple Negative Breast Cancer

Metabolic dysregulation as a reflection of specific metabolite production and its utilization is a common feature of many human neoplasms. Melatonin, an indoleamine that is highly available during darkness, has a variety of metabolic functions in solid tumors. Because plasma metabolites undergo circadian changes, we investigated the role of melatonin on the profile of amino acids (AAs), biogenic amines, carnitines, sphingolipids, and hexoses present in the plasma of mice bearing xenograft triple negative breast cancer (MDA-MB-231 cells) over 24 h. Plasma concentrations of nine AAs were reduced by melatonin, especially during the light phase, with a profile closer to that of non-breast cancer (BC) animals. With respect to acylcarnitine levels, melatonin reduced 12 out of 24 molecules in BC-bearing animals compared to their controls, especially at 06:00 h and 15:00 h. Importantly, melatonin reduced the concentrations of asymmetric dimethylarginine, carnosine, histamine, kynurenine, methionine sulfoxide, putrescine, spermidine, spermine, and symmetric dimethylarginine, which are associated with the BC metabolite sets. Melatonin also led to reduced levels of sphingomyelins and hexoses, which showed distinct daily variations over 24 h. These results highlight the role of melatonin in controlling the levels of plasma metabolites in human BC xenografts, which may impact cancer bioenergetics, in addition to emphasizing the need for a more accurate examination of its metabolomic changes at different time points.

Unraveling the effect of intra- and intercellular processes on acetaminophen-induced liver injury

In high dosages, acetaminophen (APAP) can cause severe liver damage, but susceptibility to liver failure varies across individuals and is influenced by factors such as health status. Because APAP-induced liver injury and recovery is regulated by an intricate system of intra- and extracellular molecular signaling, we here aim to quantify the importance of specific modules in determining the outcome after an APAP insult and of potential targets for therapies that mitigate adversity. For this purpose, we integrated hepatocellular acetaminophen metabolism, DNA damage response induction and cell fate into a multiscale mechanistic liver lobule model which involves various cell types, such as hepatocytes, residential Kupffer cells and macrophages. Our model simulations show that zonal differences in metabolism and detoxification efficiency are essential determinants of necrotic damage. Moreover, the extent of senescence, which is regulated by intracellular processes and triggered by extracellular signaling, influences the potential to recover. In silico therapies at early and late time points after APAP insult indicated that prevention of necrotic damage is most beneficial for recovery, whereas interference with regulation of senescence promotes regeneration in a less pronounced way.

ChemChaste: Simulating spatially inhomogeneous biochemical reaction-diffusion systems for modeling cell-environment feedbacks

BACKGROUND

Spatial organization plays an important role in the function of many biological systems, from cell fate specification in animal development to multistep metabolic conversions in microbial communities. The study of such systems benefits from the use of spatially explicit computational models that combine a discrete description of cells with a continuum description of one or more chemicals diffusing within a surrounding bulk medium. These models allow the in silico testing and refinement of mechanistic hypotheses. However, most existing models of this type do not account for concurrent bulk and intracellular biochemical reactions and their possible coupling.

CONCLUSIONS

Here, we describe ChemChaste, an extension for the open-source C++ computational biology library Chaste. ChemChaste enables the spatial simulation of both multicellular and bulk biochemistry by expanding on Chaste's existing capabilities. In particular, ChemChaste enables (i) simulation of an arbitrary number of spatially diffusing chemicals, (ii) spatially heterogeneous chemical diffusion coefficients, and (iii) inclusion of both bulk and intracellular biochemical reactions and their coupling. ChemChaste also introduces a file-based interface that allows users to define the parameters relating to these functional features without the need to interact directly with Chaste's core C++ code. We describe ChemChaste and demonstrate its functionality using a selection of chemical and biochemical exemplars, with a focus on demonstrating increased ability in modeling bulk chemical reactions and their coupling with intracellular reactions.

AVAILABILITY AND IMPLEMENTATION

ChemChaste version 1.0 is a free, open-source C++ library, available via GitHub at https://github.com/OSS-Lab/ChemChaste under the BSD license, on the Zenodo archive at zendodo doi, as well as on BioTools (biotools:chemchaste) and SciCrunch (RRID:SCR022208) databases.

Guidance by followers ensures long-range coordination of cell migration through α-catenin mechanoperception

Morphogenesis, wound healing, and some cancer metastases depend upon the migration of cell collectives that need to be guided to their destination as well as coordinated with other cell movements. During zebrafish gastrulation, the extension of the embryonic axis is led by the mesendodermal polster that migrates toward the animal pole, followed by the axial mesoderm that undergoes convergence and extension. Here, we investigate how polster cells are guided toward the animal pole. Using a combination of precise laser ablations, advanced transplants, and functional as well as in silico approaches, we establish that each polster cell is oriented by its immediate follower cells. Each cell perceives the migration of followers, through E-cadherin/α-catenin mechanotransduction, and aligns with them. Therefore, directional information propagates from cell to cell over the whole tissue. Such guidance of migrating cells by followers ensures long-range coordination of movements and developmental robustness.

Density-Dependent Migration Characteristics of Cancer Cells Driven by Pseudopod Interaction

The ability of cancer cells to invade neighboring tissue from primary tumors is an important determinant of metastatic behavior. Quantification of cell migration characteristics such as migration speed and persistence helps to understand the requirements for such invasiveness. One factor that may influence invasion is how local tumor cell density shapes cell migration characteristics, which we here investigate with a combined experimental and computational modeling approach. First, we generated and analyzed time-lapse imaging data on two aggressive Triple-Negative Breast Cancer (TNBC) cell lines, HCC38 and Hs578T, during 2D migration assays at various cell densities. HCC38 cells exhibited a counter-intuitive increase in speed and persistence with increasing density, whereas Hs578T did not exhibit such an increase. Moreover, HCC38 cells exhibited strong cluster formation with active pseudopod-driven migration, especially at low densities, whereas Hs578T cells maintained a dispersed positioning. In order to obtain a mechanistic understanding of the density-dependent cell migration characteristics and cluster formation, we developed realistic spatial simulations using a Cellular Potts Model (CPM) with an explicit description of pseudopod dynamics. Model analysis demonstrated that pseudopods exerting a pulling force on the cell and interacting via increased adhesion at pseudopod tips could explain the experimentally observed increase in speed and persistence with increasing density in HCC38 cells. Thus, the density-dependent migratory behavior could be an emergent property of single-cell characteristics without the need for additional mechanisms. This implies that pseudopod dynamics and interaction may play a role in the aggressive nature of cancers through mediating dispersal.

BioSANS: A software package for symbolic and numeric biological simulation

Modeling biochemical systems can provide insights into behaviors that are difficult to observe or understand. It requires software, programming, and understanding of the system to build a model and study it. Softwares exist for systems biology modeling, but most support only certain types of modeling tasks. Desirable features including ease in preparing input, symbolic or analytical computation, parameter estimation, graphical user interface, and systems biology markup language (SBML) support are not seen concurrently in one software package. In this study, we developed a python-based software that supports these features, with both deterministic and stochastic propagations. The software can be used by graphical user interface, command line, or as a python import. We also developed a semi-programmable and intuitively easy topology input method for the biochemical reactions. We tested the software with semantic and stochastic SBML test cases. Tests on symbolic solution and parameter estimation were also included. The software we developed is reliable, well performing, convenient to use, and compliant with most of the SBML tests. So far it is the only systems biology software that supports symbolic, deterministic, and stochastic modeling in one package that also features parameter estimation and SBML support. This work offers a comprehensive set of tools and allows for better availability and accessibility for studying kinetics and dynamics in biochemical systems.

The Role of in silico Research in Developing Nanoparticle-Based Therapeutics

Nanoparticles (NPs) hold great potential as therapeutics, particularly in the realm of drug delivery. They are effective at functional cargo delivery and offer a great degree of amenability that can be used to offset toxic side effects or to target drugs to specific regions in the body. However, there are many challenges associated with the development of NP-based drug formulations that hamper their successful clinical translation. Arguably, the most significant barrier in the way of efficacious NP-based drug delivery systems is the tedious and time-consuming nature of NP formulation—a process that needs to account for downstream effects, such as the onset of potential toxicity or immunogenicity, in vivo biodistribution and overall pharmacokinetic profiles, all while maintaining desirable therapeutic outcomes. Computational and AI-based approaches have shown promise in alleviating some of these restrictions. Via predictive modeling and deep learning, in silico approaches have shown the ability to accurately model NP-membrane interactions and cellular uptake based on minimal data, such as the physicochemical characteristics of a given NP. More importantly, machine learning allows computational models to predict how specific changes could be made to the physicochemical characteristics of a NP to improve functional aspects, such as drug retention or endocytosis. On a larger scale, they are also able to predict the in vivo pharmacokinetics of NP-encapsulated drugs, predicting aspects such as circulatory half-life, toxicity, and biodistribution. However, the convergence of nanomedicine and computational approaches is still in its infancy and limited in its applicability. The interactions between NPs, the encapsulated drug and the body form an intricate network of interactions that cannot be modeled with absolute certainty. Despite this, rapid advancements in the area promise to deliver increasingly powerful tools capable of accelerating the development of advanced nanoscale therapeutics. Here, we describe computational approaches that have been utilized in the field of nanomedicine, focusing on approaches for NP design and engineering.

Role of senescent cells in the motile behavior of active, non-senescent cells in confluent populations

Characteristics of cell migration in a confluent population depend on the nature of cell-to-cell interactions as well as cell-intrinsic properties such as the directional persistence in crawling. In addition, biological tissues (or cell cultures) almost always carry anisotropies and they too can significantly affect cell motility. In the light of this viewpoint, the emergence of cellular senescences in a confluent population of active cells raises an interesting question. Cellular senescence is a process through which a cell enters a permanent growth-arrest state and generally exhibits a dramatic body expansion. Therefore, randomly emerging senescent cells transform an initially homogeneous cell population to a “binary mixture” of two distinct cell types. Here, using in vitro cultures of MDA-MB-231 cells we investigate how spatially localized cellular senescence affect the motility of active cells within a confluent population. Importantly, we estimate the intercellular surface energy of the interface between non-senescent and senescent MDA-MB-231 cells by combining the analysis on the motile behaviors of non-senescent cells encircling senescent cells and the result of extensive numerical simulations of a cellular Potts model. We find that the adhesion of normal cells to senescent cells is much weaker than that among normal cells and that the ‘arclength’ traveled by a normal cell along the boundary of a senescent cell, on average, is several times greater than the persistence length of normal cell in a densely packed homogeneous population. The directional persistent time of normal cell during its contact with a senescent cell also increases significantly. We speculate that the phenomenon could be a general feature associated with senescent cells as the enormous expansion of senescent cell’s membrane would inevitably decrease the density of cell adhesion molecules.

Severe COVID-19 Shares a Common Neutrophil Activation Signature with Other Acute Inflammatory States

Severe COVID-19 patients present a clinical and laboratory overlap with other hyperinflammatory conditions such as hemophagocytic lymphohistiocytosis (HLH). However, the underlying mechanisms of these conditions remain to be explored. Here, we investigated the transcriptome of 1596 individuals, including patients with COVID-19 in comparison to healthy controls, other acute inflammatory states (HLH, multisystem inflammatory syndrome in children [MIS-C], Kawasaki disease [KD]), and different respiratory infections (seasonal coronavirus, influenza, bacterial pneumonia). We observed that COVID-19 and HLH share immunological pathways (cytokine/chemokine signaling and neutrophil-mediated immune responses), including gene signatures that stratify COVID-19 patients admitted to the intensive care unit (ICU) and COVID-19_nonICU patients. Of note, among the common differentially expressed genes (DEG), there is a cluster of neutrophil-associated genes that reflects a generalized hyperinflammatory state since it is also dysregulated in patients with KD and bacterial pneumonia. These genes are dysregulated at the protein level across several COVID-19 studies and form an interconnected network with differentially expressed plasma proteins that point to neutrophil hyperactivation in COVID-19 patients admitted to the intensive care unit. scRNAseq analysis indicated that these genes are specifically upregulated across different leukocyte populations, including lymphocyte subsets and immature neutrophils. Artificial intelligence modeling confirmed the strong association of these genes with COVID-19 severity. Thus, our work indicates putative therapeutic pathways for intervention.

Computational modelling of cell motility modes emerging from cell-matrix adhesion dynamics

Lymphocytes have been described to perform different motility patterns such as Brownian random walks, persistent random walks, and Lévy walks. Depending on the conditions, such as confinement or the distribution of target cells, either Brownian or Lévy walks lead to more efficient interaction with the targets. The diversity of these motility patterns may be explained by an adaptive response to the surrounding extracellular matrix (ECM). Indeed, depending on the ECM composition, lymphocytes either display a floating motility without attaching to the ECM, or sliding and stepping motility with respectively continuous or discontinuous attachment to the ECM, or pivoting behaviour with sustained attachment to the ECM. Moreover, on the long term, lymphocytes either perform a persistent random walk or a Brownian-like movement depending on the ECM composition. How the ECM affects cell motility is still incompletely understood. Here, we integrate essential mechanistic details of the lymphocyte-matrix adhesions and lymphocyte intrinsic cytoskeletal induced cell propulsion into a Cellular Potts model (CPM). We show that the combination of de novo cell-matrix adhesion formation, adhesion growth and shrinkage, adhesion rupture, and feedback of adhesions onto cell propulsion recapitulates multiple lymphocyte behaviours, for different lymphocyte subsets and various substrates. With an increasing attachment area and increased adhesion strength, the cells’ speed and persistence decreases. Additionally, the model predicts random walks with short-term persistent but long-term subdiffusive properties resulting in a pivoting type of motility. For small adhesion areas, the spatial distribution of adhesions emerges as a key factor influencing cell motility. Small adhesions at the front allow for more persistent motility than larger clusters at the back, despite a similar total adhesion area. In conclusion, we present an integrated framework to simulate the effects of ECM proteins on cell-matrix adhesion dynamics. The model reveals a sufficient set of principles explaining the plasticity of lymphocyte motility.

During immunosurveillance, lymphocytes patrol through tissues to interact with cancer cells, other immune cells, and pathogens. The efficiency of this process depends on the kinds of trajectories taken, ranging from simple Brownian walks to Lévy walks. The composition of the extracellular matrix (ECM), a network of macromolecules, affects the formation of cell-matrix adhesions, thus strongly influencing the way lymphocytes move. Here, we present a model of lymphocyte motility driven by adhesions that grow, shrink and rupture in response to the ECM and cellular forces. Compared to other models, our model is computationally light making it suitable for generating long term cell track data, while still capturing actin dynamics and adhesion turnover. Our model suggests that cell motility is affected by the force required to break adhesions and the rate at which new adhesions form. Adhesions can promote cell protrusion by inhibiting retrograde actin flow. After introducing this effect into the model, we found that it reduces the cellular diffusivity and that it promotes stick-slip behaviour. Furthermore, location and size of adhesion clusters determined cell persistence. Overall, our model explains the plasticity of lymphocyte behaviour in response to the ECM.

Transcriptomic Response to Acidosis Reveals Its Contribution to Bone Metastasis in Breast Cancer Cells

Bone is the most common site of metastasis in breast cancer. Metastasis is promoted by acidosis, which is associated with osteoporosis. To investigate how acidosis could promote bone metastasis, we compared differentially expressed genes (DEGs) in MDA-MB-231 cancer cells in acidosis, bone metastasis, and bone metastatic tumors. The DEGs were identified using Biojupies and GEO2R. The expression profiles were assessed with Morpheus. The overlapping DEGs between acidosis and bone metastasis were compared to the bulk of the DEGs in terms of the most important genes and enriched terms using CytoHubba and STRING. The expression of the genes in this overlap filtered by secreted proteins was assessed in the osteoporosis secretome. The analysis revealed that acidosis-associated transcriptomic changes were more similar to bone metastasis than bone metastatic tumors. Extracellular matrix (ECM) organization would be the main biological process shared between acidosis and bone metastasis. The secretome genes upregulated in acidosis, bone metastasis, and osteoporosis-associated mesenchymal stem cells are enriched for ECM organization and angiogenesis. Therefore, acidosis may be more important in the metastatic niche than in the primary tumor. Acidosis may contribute to bone metastasis by promoting ECM organization. Untreated osteoporosis could favor bone metastasis through the increased secretion of ECM organization proteins.

The proteomic landscape of ovarian cancer cells in response to melatonin

Ovarian cancer (OC) is the most lethal gynecological malignancy with a highly negative prognosis. Melatonin is an indoleamine secreted by the pineal gland during darkness and has shown antitumor activity in both in vitro and in vivo experiments. Herein, we investigated the influence of melatonin on the proteome of human ovarian carcinoma cells (SKOV-3 cell line) using the Ultimate 3000 LC Liquid NanoChromatography equipment coupled to a Q-Exactive mass spectrometry. After 48 h of treatment, melatonin induced a significant cytotoxicity especially with the highest melatonin concentration. The proteomic profile revealed 639 proteins in the control group, and 98, 110, and 128 proteins were altered by melatonin at the doses of 0.8, 1.6, and 2.4 mM, respectively. Proteins associated with the immune system and tricarboxylic acid cycle were increased in the three melatonin-exposed groups of cells. Specifically, the dose of 2.4 mM led to a reduction in molecules associated with protein synthesis, especially those of the ribosomal protein family. We also identified 28 potential genes shared between normal ovarian tissue and OC in all experimental groups, and melatonin was predicted to alter genes encoding ribosomal proteins. Notably, the set of proteins changed by melatonin was linked to a better prognosis for OC patients. We conclude that melatonin significantly alters the proteome of SKOV-3 cells by changing proteins involved with the immune response and mitochondrial metabolism. The concentration of 2.4 mM of melatonin promoted the largest number of protein changes. The evidence suggests that melatonin may be an effective therapeutic strategy against OC.

Amino Acids and IGF1 Regulation of Fish Muscle Growth Revealed by Transcriptome and microRNAome Integrative Analyses of Pacu (Piaractus mesopotamicus) Myotubes

Amino acids (AA) and IGF1 have been demonstrated to play essential roles in protein synthesis and fish muscle growth. The myoblast cell culture is useful for studying muscle regulation, and omics data have contributed enormously to understanding its molecular biology. However, to our knowledge, no study has performed the large-scale sequencing of fish-cultured muscle cells stimulated with pro-growth signals. In this work, we obtained the transcriptome and microRNAome of pacu (Piaractus mesopotamicus)-cultured myotubes treated with AA or IGF1. We identified 1228 and 534 genes differentially expressed by AA and IGF1. An enrichment analysis showed that AA treatment induced chromosomal changes, mitosis, and muscle differentiation, while IGF1 modulated IGF/PI3K signaling, metabolic alteration, and matrix structure. In addition, potential molecular markers were similarly modulated by both treatments. Muscle-miRNAs (miR-1, -133, -206 and -499) were up-regulated, especially in AA samples, and we identified molecular networks with omics integration. Two pairs of genes and miRNAs demonstrated a high-level relationship, and involvement in myogenesis and muscle growth: marcksb and miR-29b in AA, and mmp14b and miR-338-5p in IGF1. Our work helps to elucidate fish muscle physiology and metabolism, highlights potential molecular markers, and creates a perspective for improvements in aquaculture and in in vitro meat production.

A scalable solver for a stochastic, hybrid cellular automaton model of personalized breast cancer therapy

Mathematical modeling and simulation is a promising approach to personalized cancer medicine. Yet, the complexity, heterogeneity and multi-scale nature of cancer pose significant computational challenges. Coupling discrete cell-based models with continuous models using hybrid cellular automata (CA) is a powerful approach for mimicking biological complexity and describing the dynamical exchange of information across different scales. However, when clinically relevant cancer portions are taken into account, such models become computationally very expensive. While efficient parallelization techniques for continuous models exist, their coupling with discrete models, particularly CA, necessitates more elaborate solutions. Building upon FEniCS, a popular and powerful scientific computing platform for solving partial differential equations, we developed parallel algorithms to link stochastic CA with differential equations (https://bitbucket.org/HTasken/cansim). The algorithms minimize the communication between processes that share CA neighborhood values while also allowing for reproducibility during stochastic updates. We demonstrated the potential of our solution on a complex hybrid cellular automaton model of breast cancer treated with combination chemotherapy. On a single-core processor, we obtained nearly linear scaling with an increasing problem size, whereas weak parallel scaling showed moderate growth in solving time relative to increase in problem size. Finally, we applied the algorithm to a problem that is 500 times larger than previous work, allowing us to run personalized therapy simulations based on heterogeneous cell density and tumor perfusion conditions estimated from magnetic resonance imaging data on an unprecedented scale.

Size matters: tissue size as a marker for a transition between reaction-diffusion regimes in spatio-temporal distribution of morphogens

The reaction-diffusion model constitutes one of the most influential mathematical models to study distribution of morphogens in tissues. Despite its widespread use, the effect of finite tissue size on model-predicted spatio-temporal morphogen distributions has not been completely elucidated. In this study, we analytically investigated the spatio-temporal distributions of morphogens predicted by a reaction-diffusion model in a finite one-dimensional domain, as a proxy for a biological tissue, and compared it with the solution of the infinite-domain model. We explored the reduced parameter, the tissue length in units of a characteristic reaction-diffusion length, and identified two reaction-diffusion regimes separated by a crossover tissue size estimated in approximately three characteristic reaction-diffusion lengths. While above this crossover the infinite-domain model constitutes a good approximation, it breaks below this crossover, whereas the finite-domain model faithfully describes the entire parameter space. We evaluated whether the infinite-domain model renders accurate estimations of diffusion coefficients when fitted to finite spatial profiles, a procedure typically followed in fluorescence recovery after photobleaching (FRAP) experiments. We found that the infinite-domain model overestimates diffusion coefficients when the domain is smaller than the crossover tissue size. Thus, the crossover tissue size may be instrumental in selecting the suitable reaction-diffusion model to study tissue morphogenesis.

Seven challenges in the multiscale modeling of multicellular tissues

The growth and dynamics of multicellular tissues involve tightly regulated and coordinated morphogenetic cell behaviors, such as shape changes, movement, and division, which are governed by subcellular machinery and involve coupling through short- and long-range signals. A key challenge in the fields of developmental biology, tissue engineering and regenerative medicine is to understand how relationships between scales produce emergent tissue-scale behaviors. Recent advances in molecular biology, live-imaging and ex vivo techniques have revolutionized our ability to study these processes experimentally. To fully leverage these techniques and obtain a more comprehensive understanding of the causal relationships underlying tissue dynamics, computational modeling approaches are increasingly spanning multiple spatial and temporal scales, and are coupling cell shape, growth, mechanics, and signaling. Yet such models remain challenging: modeling at each scale requires different areas of technical skills, while integration across scales necessitates the solution to novel mathematical and computational problems. This review aims to summarize recent progress in multiscale modeling of multicellular tissues and to highlight ongoing challenges associated with the construction, implementation, interrogation, and validation of such models. This article is categorized under: Reproductive System Diseases > Computational Models Metabolic Diseases > Computational Models Cancer > Computational Models.

Decreased miR-497-5p Suppresses IL-6 Induced Atrophy in Muscle Cells

Interleukin-6 (IL-6) is a pro-inflammatory cytokine associated with skeletal muscle wasting in cancer cachexia. The control of gene expression by microRNAs (miRNAs) in muscle wasting involves the regulation of thousands of target transcripts. However, the miRNA-target networks associated with IL6-induced muscle atrophy remain to be characterized. Here, we show that IL-6 promotes the atrophy of C2C12 myotubes and changes the expression of 20 miRNAs (5 up-regulated and 15 down-regulated). Gene Ontology analysis of predicted miRNAs targets revealed post-transcriptional regulation of genes involved in cell differentiation, apoptosis, migration, and catabolic processes. Next, we performed a meta-analysis of miRNA-published data that identified miR-497-5p, a down-regulated miRNAs induced by IL-6, also down-regulated in other muscle-wasting conditions. We used miR-497-5p mimics and inhibitors to explore the function of miR-497-5p in C2C12 myoblasts and myotubes. We found that miR-497-5p can regulate the expression of the cell cycle genes CcnD2 and CcnE1 without affecting the rate of myoblast cellular proliferation. Notably, miR-497-5p mimics induced myotube atrophy and reduced Insr expression. Treatment with miR-497-5p inhibitors did not change the diameter of the myotubes but increased the expression of its target genes Insr and Igf1r. These genes are known to regulate skeletal muscle regeneration and hypertrophy via insulin-like growth factor pathway and were up-regulated in cachectic muscle samples. Our miRNA-regulated network analysis revealed a potential role for miR-497-5p during IL6-induced muscle cell atrophy and suggests that miR-497-5p is likely involved in a compensatory mechanism of muscle atrophy in response to IL-6.

Derivation of continuum models from discrete models of mechanical forces in cell populations

In certain discrete models of populations of biological cells, the mechanical forces between the cells are center based or vertex based on the microscopic level where each cell is individually represented. The cells are circular or spherical in a center based model and polygonal or polyhedral in a vertex based model. On a higher, macroscopic level, the time evolution of the density of the cells is described by partial differential equations (PDEs). We derive relations between the modelling on the micro and macro levels in one, two, and three dimensions by regarding the micro model as a discretization of a PDE for conservation of mass on the macro level. The forces in the micro model correspond on the macro level to a gradient of the pressure scaled by quantities depending on the cell geometry. The two levels of modelling are compared in numerical experiments in one and two dimensions.

Surface Roughness of Titanium Orthopedic Implants Alters the Biological Phenotype of Human Mesenchymal Stromal Cells

Metal orthopedic implants are largely biocompatible and generally achieve long-term structural fixation. However, some orthopedic implants may loosen over time even in the absence of infection. In vivo fixation failure is multifactorial, but the fundamental biological defect is cellular dysfunction at the host-implant interface. Strategies to reduce the risk of short- and long-term loosening include surface modifications, implant metal alloy type, and adjuvant substances such as polymethylmethacrylate cement. Surface modifications (e.g., increased surface rugosity) can increase osseointegration and biological ingrowth of orthopedic implants. However, the localized responses of cells to implant surface modifications need to be better characterized. As an in vitro model for investigating cellular responses to metallic orthopedic implants, we cultured mesenchymal stromal/stem cells on clinical-grade titanium disks (Ti6Al4V) that differed in surface roughness as high (porous structured), medium (grit blasted), and low (bead blasted). Topological characterization of clinically relevant titanium (Ti) materials combined with differential mRNA expression analyses (RNA-seq and real-time quantitative polymerase chain reaction) revealed alterations to the biological phenotype of cells cultured on titanium structures that favor early extracellular matrix production and observable responses to oxidative stress and heavy metal stress. These results provide a descriptive model for the interpretation of cellular responses at the interface between native host tissues and three-dimensionally printed modular orthopedic implants, and will guide future studies aimed at increasing the long-term retention of such materials after total joint arthroplasty. Impact statement Using an in vitro model of implant-to-cell interactions by culturing mesenchymal stromal cells (MSCs) on clinically relevant titanium materials of varying topological roughness, we identified mRNA expression patterns consistent with early extracellular matrix (ECM) production and responses to oxidative/heavy metal stress. Implants with high surface roughness may delay the differentiation and ECM formation of MSCs and alter the expression of genes sensitive to reactive oxygen species and protein kinases. In combination with ongoing animal studies, these results will guide future studies aimed at increasing the long-term retention of widely used titanium materials after total joint arthroplasty.

Navigating Multi-Scale Cancer Systems Biology Towards Model-Driven Clinical Oncology and Its Applications in Personalized Therapeutics

Rapid advancements in high-throughput omics technologies and experimental protocols have led to the generation of vast amounts of scale-specific biomolecular data on cancer that now populates several online databases and resources. Cancer systems biology models built using this data have the potential to provide specific insights into complex multifactorial aberrations underpinning tumor initiation, development, and metastasis. Furthermore, the annotation of these single- and multi-scale models with patient data can additionally assist in designing personalized therapeutic interventions as well as aid in clinical decision-making. Here, we have systematically reviewed the emergence and evolution of (i) repositories with scale-specific and multi-scale biomolecular cancer data, (ii) systems biology models developed using this data, (iii) associated simulation software for the development of personalized cancer therapeutics, and (iv) translational attempts to pipeline multi-scale panomics data for data-driven in silico clinical oncology. The review concludes that the absence of a generic, zero-code, panomics-based multi-scale modeling pipeline and associated software framework, impedes the development and seamless deployment of personalized in silico multi-scale models in clinical settings.

Iterative community-driven development of a SARS-CoV-2 tissue simulator

The 2019 novel coronavirus, SARS-CoV-2, is a pathogen of critical significance to international public health. Knowledge of the interplay between molecular-scale virus-receptor interactions, single-cell viral replication, intracellular-scale viral transport, and emergent tissue-scale viral propagation is limited. Moreover, little is known about immune system-virus-tissue interactions and how these can result in low-level (asymptomatic) infections in some cases and acute respiratory distress syndrome (ARDS) in others, particularly with respect to presentation in different age groups or pre-existing inflammatory risk factors. Given the nonlinear interactions within and among each of these processes, multiscale simulation models can shed light on the emergent dynamics that lead to divergent outcomes, identify actionable "choke points" for pharmacologic interventions, screen potential therapies, and identify potential biomarkers that differentiate patient outcomes. Given the complexity of the problem and the acute need for an actionable model to guide therapy discovery and optimization, we introduce and iteratively refine a prototype of a multiscale model of SARS-CoV-2 dynamics in lung tissue. The first prototype model was built and shared internationally as open source code and an online interactive model in under 12 hours, and community domain expertise is driving regular refinements. In a sustained community effort, this consortium is integrating data and expertise across virology, immunology, mathematical biology, quantitative systems physiology, cloud and high performance computing, and other domains to accelerate our response to this critical threat to international health. More broadly, this effort is creating a reusable, modular framework for studying viral replication and immune response in tissues, which can also potentially be adapted to related problems in immunology and immunotherapy.

Multidimensional Single-Nuclei RNA-Seq Reconstruction of Adipose Tissue Reveals Adipocyte Plasticity Underlying Thermogenic Response

Adipose tissue has been classified based on its morphology and function as white, brown, or beige/brite. It plays an essential role as a regulator of systemic metabolism through paracrine and endocrine signals. Recently, multiple adipocyte subtypes have been revealed using RNA sequencing technology, going beyond simply defined morphology but also by their cellular origin, adaptation to metabolic stress, and plasticity. Here, we performed an in-depth analysis of publicly available single-nuclei RNAseq from adipose tissue and utilized a workflow template to characterize adipocyte plasticity, heterogeneity, and secretome profiles. The reanalyzed dataset led to the identification of different subtypes of adipocytes including three subpopulations of thermogenic adipocytes, and provided a characterization of distinct transcriptional profiles along the adipocyte trajectory under thermogenic challenges. This study provides a useful resource for further investigations regarding mechanisms related to adipocyte plasticity and trans-differentiation.

Aging whole blood transcriptome reveals candidate genes for SARS-CoV-2-related vascular and immune alterations

Abstract

The risk of severe COVID-19 increases with age as older patients are at highest risk. Thus, there is an urgent need to identify how severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) interacts with blood components during aging. We investigated the whole blood transcriptome from the Genotype-Tissue Expression (GTEx) database to explore differentially expressed genes (DEGs) translated into proteins interacting with viral proteins during aging. From 22 DEGs in aged blood, FASLG, CTSW, CTSE, VCAM1, and BAG3 were associated with immune response, inflammation, cell component and adhesion, and platelet activation/aggregation. Males and females older than 50 years old overexpress FASLG, possibly inducing a hyperinflammatory cascade. The expression of cathepsins (CTSW and CTSE) and the anti-apoptotic co-chaperone molecule BAG3 also increased throughout aging in both genders. By exploring single-cell RNA-sequencing data from peripheral blood of SARS-CoV-2-infected patients, we found FASLG and CTSW expressed in natural killer cells and CD8 + T lymphocytes, whereas BAG3 was expressed mainly in CD4 + T cells, naive T cells, and CD14 + monocytes. In addition, T cell exhaustion was associated with increased expression of CCL4L2 and DUSP4 over blood aging. LAG3, PDCD1, TIGIT, VCAM1, HLA-DRA, and TOX also increased in individuals aged 60–69 years old; conversely, the RGS2 gene decreased with aging. We further identified a distinct gene expression profile associated with type I interferon signaling following blood aging. These results revealed changes in blood molecules potentially related to SARS-CoV-2 infection throughout aging, emphasizing them as therapeutic candidates for aggressive clinical manifestation of COVID-19.

Key messages

• Prediction of host-viral interactions in the whole blood transcriptome during aging.

• Expression levels of FASLG, CTSW, CTSE, VCAM1, and BAG3 increase in aged blood.

• Blood interactome reveals targets involved with immune response, inflammation, and blood clots.

• SARS-CoV-2-infected patients with high viral load showed FASLG overexpression.

• Gene expression profile associated with T cell exhaustion and type I interferon signaling were affected with blood aging.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00109-021-02161-4.

Drug Repositioning Based on the Reversal of Gene Expression Signatures Identifies TOP2A as a Therapeutic Target for Rectal Cancer

Rectal cancer is a common disease with high mortality rates and limited therapeutic options. Here we combined the gene expression signatures of rectal cancer patients with the reverse drug-induced gene-expression profiles to identify drug repositioning candidates for cancer therapy. Among the predicted repurposable drugs, topoisomerase II inhibitors (doxorubicin, teniposide, idarubicin, mitoxantrone, and epirubicin) presented a high potential to reverse rectal cancer gene expression signatures. We showed that these drugs effectively reduced the growth of colorectal cancer cell lines closely representing rectal cancer signatures. We also found a clear correlation between topoisomerase 2A (TOP2A) gene copy number or expression levels with the sensitivity to topoisomerase II inhibitors. Furthermore, CRISPR-Cas9 and shRNA screenings confirmed that loss-of-function of the TOP2A has the highest efficacy in reducing cellular proliferation. Finally, we observed significant TOP2A copy number gains and increased expression in independent cohorts of rectal cancer patients. These findings can be translated into clinical practice to evaluate TOP2A status for targeted and personalized therapies based on topoisomerase II inhibitors in rectal cancer patients.

Cellular Tango: how extracellular matrix adhesion choreographs Rac-Rho signaling and cell movement

The small GTPases Rac and Rho are known to regulate eukaryotic cell shape, promoting front protrusion (Rac) or rear retraction (Rho) of the cell edge. Such cell deformation changes the contact and adhesion of cell to the extracellular matrix (ECM), while ECM signaling through integrin receptors also affects GTPase activity. We develop and investigate a model for this three-way feedback loop in 1D and 2D spatial domains, as well as in a fully deforming 2D cell shapes with detailed adhesion-bond biophysics. The model consists of reaction-diffusion equations solved numerically with open-source software, Morpheus, and with custom-built cellular Potts model simulations. We find a variety of patterns and cell behaviors, including persistent polarity, flipped front-back cell polarity oscillations, spiral waves, and random protrusion-retraction. We show that the observed spatial patterns depend on the cell shape, and vice versa.

The network interplay of interferon and Toll-like receptor signaling pathways in the anti-Candida immune response

Fungal infections represent a major global health problem affecting over a billion people that kills more than 1.5 million annually. In this study, we employed an integrative approach to reveal the landscape of the human immune responses to Candida spp. through meta-analysis of microarray, bulk, and single-cell RNA sequencing (scRNA-seq) data for the blood transcriptome. We identified across these different studies a consistent interconnected network interplay of signaling molecules involved in both Toll-like receptor (TLR) and interferon (IFN) signaling cascades that is activated in response to different Candida species (C. albicans, C. auris, C. glabrata, C. parapsilosis, and C. tropicalis). Among these molecules are several types I IFN, indicating an overlap with antiviral immune responses. scRNA-seq data confirmed that genes commonly identified by the three transcriptomic methods show cell type-specific expression patterns in various innate and adaptive immune cells. These findings shed new light on the anti-Candida immune response, providing putative molecular pathways for therapeutic intervention.

Neighbor-enhanced diffusivity in dense, cohesive cell populations

The dispersal or mixing of cells within cellular tissue is a crucial property for diverse biological processes, ranging from morphogenesis, immune action, to tumor metastasis. With the phenomenon of ‘contact inhibition of locomotion,’ it is puzzling how cells achieve such processes within a densely packed cohesive population. Here we demonstrate that a proper degree of cell-cell adhesiveness can, intriguingly, enhance the super-diffusive nature of individual cells. We systematically characterize the migration trajectories of crawling MDA-MB-231 cell lines, while they are in several different clustering modes, including freely crawling singles, cohesive doublets of two cells, quadruplets, and confluent population on two-dimensional substrate. Following data analysis and computer simulation of a simple cellular Potts model, which faithfully recapitulated all key experimental observations such as enhanced diffusivity as well as periodic rotation of cell-doublets and cell-quadruplets with mixing events, we found that proper combination of active self-propelling force and cell-cell adhesion is sufficient for generating the observed phenomena. Additionally, we found that tuning parameters for these two factors covers a variety of different collective dynamic states.

Dispersal or movement of cells within dense biological tissue is essential for diverse biological processes, ranging from pattern formation, immune action, to tumor metastasis. However, it is quite puzzling how cells acquire such ability when they are supposedly “caged” by neighboring cells. Here, we report an unusual property of (MDA-MB-231) breast cancer cells that diffuse more persistently within a densely packed population than when they are free to crawl around with little interference. This property is rather surprising since they prefer to stick together, forming clusters. Interestingly, however, we find that having sticky neighbors not only makes two active cells in contact periodically rotate, reminiscent of a ballroom dance, but also enhances the persistence of the cells within a dense population. These intriguing phenomena appear to be universal as they can be generated by a simple cellular Potts model with appropriate combination of active self-propulsion and cell-cell adhesion force.

Generation of multicellular spatiotemporal models of population dynamics from ordinary differential equations, with applications in viral infection

BACKGROUND

The biophysics of an organism span multiple scales from subcellular to organismal and include processes characterized by spatial properties, such as the diffusion of molecules, cell migration, and flow of intravenous fluids. Mathematical biology seeks to explain biophysical processes in mathematical terms at, and across, all relevant spatial and temporal scales, through the generation of representative models. While non-spatial, ordinary differential equation (ODE) models are often used and readily calibrated to experimental data, they do not explicitly represent the spatial and stochastic features of a biological system, limiting their insights and applications. However, spatial models describing biological systems with spatial information are mathematically complex and computationally expensive, which limits the ability to calibrate and deploy them and highlights the need for simpler methods able to model the spatial features of biological systems.

RESULTS

In this work, we develop a formal method for deriving cell-based, spatial, multicellular models from ODE models of population dynamics in biological systems, and vice versa. We provide examples of generating spatiotemporal, multicellular models from ODE models of viral infection and immune response. In these models, the determinants of agreement of spatial and non-spatial models are the degree of spatial heterogeneity in viral production and rates of extracellular viral diffusion and decay. We show how ODE model parameters can implicitly represent spatial parameters, and cell-based spatial models can generate uncertain predictions through sensitivity to stochastic cellular events, which is not a feature of ODE models. Using our method, we can test ODE models in a multicellular, spatial context and translate information to and from non-spatial and spatial models, which help to employ spatiotemporal multicellular models using calibrated ODE model parameters. We additionally investigate objects and processes implicitly represented by ODE model terms and parameters and improve the reproducibility of spatial, stochastic models.

CONCLUSION

We developed and demonstrate a method for generating spatiotemporal, multicellular models from non-spatial population dynamics models of multicellular systems. We envision employing our method to generate new ODE model terms from spatiotemporal and multicellular models, recast popular ODE models on a cellular basis, and generate better models for critical applications where spatial and stochastic features affect outcomes.

Disruption of RING and PHD Domains of TRIM28 Evokes Differentiation in Human iPSCs

TRIM28, a multi-domain protein, is crucial in the development of mouse embryos and the maintenance of embryonic stem cells’ (ESC) self-renewal potential. As the epigenetic factor modulating chromatin structure, TRIM28 regulates the expression of numerous genes and is associated with progression and poor prognosis in many types of cancer. Because of many similarities between highly dedifferentiated cancer cells and normal pluripotent stem cells, we applied human induced pluripotent stem cells (hiPSC) as a model for stemness studies. For the first time in hiPSC, we analyzed the function of individual TRIM28 domains. Here we demonstrate the essential role of a really interesting new gene (RING) domain and plant homeodomain (PHD) in regulating pluripotency maintenance and self-renewal capacity of hiPSC. Our data indicate that mutation within the RING or PHD domain leads to the loss of stem cell phenotypes and downregulation of the FGF signaling. Moreover, impairment of RING or PHD domain results in decreased proliferation and impedes embryoid body formation. In opposition to previous data indicating the impact of phosphorylation on TRIM28 function, our data suggest that TRIM28 phosphorylation does not significantly affect the pluripotency and self-renewal maintenance of hiPSC. Of note, iPSC with disrupted RING and PHD functions display downregulation of genes associated with tumor metastasis, which are considered important targets in cancer treatment. Our data suggest the potential use of RING and PHD domains of TRIM28 as targets in cancer therapy.

An insight on the impact of teleost whole genome duplication on the regulation of the molecular networks controlling skeletal muscle growth

Fish muscle growth is a complex process regulated by multiple pathways, resulting on the net accumulation of proteins and the activation of myogenic progenitor cells. Around 350–320 million years ago, teleost fish went through a specific whole genome duplication (WGD) that expanded the existent gene repertoire. Duplicated genes can be retained by different molecular mechanisms such as subfunctionalization, neofunctionalization or redundancy, each one with different functional implications. While the great majority of ohnolog genes have been identified in the teleost genomes, the effect of gene duplication in the fish physiology is still not well characterized. In the present study we studied the effect of WGD on the transcription of the duplicated components controlling muscle growth. We compared the expression of lineage-specific ohnologs related to myogenesis and protein balance in the fast-skeletal muscle of pacus (Piaractus mesopotamicus—Ostariophysi) and Nile tilapias (Oreochromis niloticus—Acanthopterygii) fasted for 4 days and refed for 3 days. We studied the expression of 20 ohnologs and found that in the great majority of cases, duplicated genes had similar expression profiles in response to fasting and refeeding, indicating that their functions during growth have been conserved during the period after the WGD. Our results suggest that redundancy might play a more important role in the retention of ohnologs of regulatory pathways than initially thought. Also, comparison to non-duplicated orthologs showed that it might not be uncommon for the duplicated genes to gain or loss new regulatory elements simultaneously. Overall, several of duplicated ohnologs have similar transcription profiles in response to pro-growth signals suggesting that evolution tends to conserve ohnolog regulation during muscle development and that in the majority of ohnologs related to muscle growth their functions might be very similar.

HCV Spread Kinetics Reveal Varying Contributions of Transmission Modes to Infection Dynamics

The hepatitis C virus (HCV) is capable of spreading within a host by two different transmission modes: cell-free and cell-to-cell. However, the contribution of each of these transmission mechanisms to HCV spread is unknown. To dissect the contribution of these different transmission modes to HCV spread, we measured HCV lifecycle kinetics and used an in vitro spread assay to monitor HCV spread kinetics after a low multiplicity of infection in the absence and presence of a neutralizing antibody that blocks cell-free spread. By analyzing these data with a spatially explicit mathematical model that describes viral spread on a single-cell level, we quantified the contribution of cell-free, and cell-to-cell spread to the overall infection dynamics and show that both transmission modes act synergistically to enhance the spread of infection. Thus, the simultaneous occurrence of both transmission modes represents an advantage for HCV that may contribute to viral persistence. Notably, the relative contribution of each viral transmission mode appeared to vary dependent on different experimental conditions and suggests that viral spread is optimized according to the environment. Together, our analyses provide insight into the spread dynamics of HCV and reveal how different transmission modes impact each other.

A Keller-Segel model for C elegans L1 aggregation

We describe a mathematical model for the aggregation of starved first-stage C elegans larvae (L1s). We propose that starved L1s produce and respond chemotactically to two labile diffusible chemical signals, a short-range attractant and a longer range repellent. This model takes the mathematical form of three coupled partial differential equations, one that describes the movement of the worms and one for each of the chemical signals. Numerical solution of these equations produced a pattern of aggregates that resembled that of worm aggregates observed in experiments. We also describe the identification of a sensory receptor gene, srh-2, whose expression is induced under conditions that promote L1 aggregation. Worms whose srh-2 gene has been knocked out form irregularly shaped aggregates. Our model suggests this phenotype may be explained by the mutant worms slowing their movement more quickly than the wild type.

The relationship between cytokine and neutrophil gene network distinguishes SARS-CoV-2-infected patients by sex and age

The fact that the COVID-19 fatality rate varies by sex and age is poorly understood. Notably, the outcome of SARS-CoV-2 infections mostly depends on the control of cytokine storm and the increasingly recognized pathological role of uncontrolled neutrophil activation. Here, we used an integrative approach with publicly available RNA-Seq data sets of nasopharyngeal swabs and peripheral blood leukocytes from patients with SARS-CoV-2, according to sex and age. Female and young patients infected by SARS-CoV-2 exhibited a larger number of differentially expressed genes (DEGs) compared with male and elderly patients, indicating a stronger immune modulation. Among them, we found an association between upregulated cytokine/chemokine- and downregulated neutrophil-related DEGs. This was correlated with a closer relationship between female and young subjects, while the relationship between male and elderly patients was closer still. The association between these cytokine/chemokines and neutrophil DEGs is marked by a strongly correlated interferome network. Here, female patients exhibited reduced transcriptional levels of key proinflammatory/neutrophil-related genes, such as CXCL8 receptors (CXCR1 and CXCR2), IL-1β, S100A9, ITGAM, and DBNL, compared with male patients. These genes are well known to be protective against inflammatory damage. Therefore, our work suggests specific immune-regulatory pathways associated with sex and age of patients infected with SARS-CoV-2 and provides a possible association between inverse modulation of cytokine/chemokine and neutrophil transcriptional signatures.

Artistoo, a library to build, share, and explore simulations of cells and tissues in the web browser

The cellular Potts model (CPM) is a powerful in silico method for simulating biological processes at tissue scale. Their inherently graphical nature makes CPMs very accessible in theory, but in practice, they are mostly implemented in specialised frameworks users need to master before they can run simulations. We here present Artistoo (Artificial Tissue Toolbox), a JavaScript library for building ‘explorable’ CPM simulations where viewers can change parameters interactively, exploring their effects in real time. Simulations run directly in the web browser and do not require third-party software, plugins, or back-end servers. The JavaScript implementation imposes no major performance loss compared to frameworks written in C++; Artistoo remains sufficiently fast for interactive, real-time simulations. Artistoo provides an opportunity to unlock CPM models for a broader audience: interactive simulations can be shared via a URL in a zero-install setting. We discuss applications in CPM research, science dissemination, open science, and education.

Modeling the Dynamics of T-Cell Development in the Thymus

The thymus hosts the development of a specific type of adaptive immune cells called T cells. T cells orchestrate the adaptive immune response through recognition of antigen by the highly variable T-cell receptor (TCR). T-cell development is a tightly coordinated process comprising lineage commitment, somatic recombination of Tcr gene loci and selection for functional, but non-self-reactive TCRs, all interspersed with massive proliferation and cell death. Thus, the thymus produces a pool of T cells throughout life capable of responding to virtually any exogenous attack while preserving the body through self-tolerance. The thymus has been of considerable interest to both immunologists and theoretical biologists due to its multi-scale quantitative properties, bridging molecular binding, population dynamics and polyclonal repertoire specificity. Here, we review experimental strategies aimed at revealing quantitative and dynamic properties of T-cell development and how they have been implemented in mathematical modeling strategies that were reported to help understand the flexible dynamics of the highly dividing and dying thymic cell populations. Furthermore, we summarize the current challenges to estimating in vivo cellular dynamics and to reaching a next-generation multi-scale picture of T-cell development.

Prediction of SARS-CoV Interaction with Host Proteins during Lung Aging Reveals a Potential Role for TRIB3 in COVID-19

COVID-19 is prevalent in the elderly. Old individuals are more likely to develop pneumonia and respiratory failure due to alveolar damage, suggesting that lung senescence may increase the susceptibility to SARS-CoV-2 infection and replication. Considering that human coronavirus (HCoVs; SARS-CoV-2 and SARS-CoV) require host cellular factors for infection and replication, we analyzed Genotype-Tissue Expression (GTEx) data to test whether lung aging is associated with transcriptional changes in human protein-coding genes that potentially interact with these viruses. We found decreased expression of the gene tribbles homolog 3 (TRIB3) during aging in male individuals, and its protein was predicted to interact with HCoVs nucleocapsid protein and RNA-dependent RNA polymerase. Using publicly available lung single-cell data, we found TRIB3 expressed mainly in alveolar epithelial cells that express SARS-CoV-2 receptor ACE2. Functional enrichment analysis of age-related genes, in common with SARS-CoV-induced perturbations, revealed genes associated with the mitotic cell cycle and surfactant metabolism. Given that TRIB3 was previously reported to decrease virus infection and replication, the decreased expression of TRIB3 in aged lungs may help explain why older male patients are related to more severe cases of the COVID-19. Thus, drugs that stimulate TRIB3 expression should be evaluated as a potential therapy for the disease.

Tools for time-course simulation in systems biology: a brief overview

Dynamic modeling of biological systems is essential for understanding all properties of a given organism as it allows us to look not only at the static picture of an organism but also at its behavior under various conditions. With the increasing amount of experimental data, the number of tools that enable dynamic analysis also grows. However, various tools are based on different approaches, use different types of data and offer different functions for analyses; so it can be difficult to choose the most suitable tool for a selected type of model. Here, we bring a brief overview containing descriptions of 50 tools for the reconstruction of biological models, their time-course simulation and dynamic analysis. We examined each tool using test data and divided them based on the qualitative and quantitative nature of the mathematical apparatus they use.

Towards understanding the messengers of extracellular space: Computational models of outside-in integrin reaction networks

The interactions between cells and their extracellular matrix (ECM) are critically important for homeostatic control of cell growth, proliferation, differentiation and apoptosis. Transmembrane integrin molecules facilitate the communication between ECM and the cell. Since the characterization of integrins in the late 1980s, there has been great advancement in understanding the function of integrins at different subcellular levels. However, the versatility in molecular pathways integrins are involved in, the high diversity in their interaction partners both outside and inside the cell as well as on the cell membrane and the short lifetime of events happening at the cell-ECM interface make it difficult to elucidate all the details regarding integrin function experimentally. To overcome the experimental challenges and advance the understanding of integrin biology, computational modeling tools have been used extensively. In this review, we summarize the computational models of integrin signaling while we explain the function of integrins at three main subcellular levels (outside the cell, cell membrane, cytosol). We also discuss how these computational modeling efforts can be helpful in other disciplines such as biomaterial design. As such, this review is a didactic modeling summary for biomaterial researchers interested in complementing their experimental work with computational tools or for seasoned computational scientists that would like to advance current in silico integrin models.

Reoccurring neural stem cell divisions in the adult zebrafish telencephalon are sufficient for the emergence of aggregated spatiotemporal patterns

Regulation of quiescence and cell cycle entry is pivotal for the maintenance of stem cell populations. Regulatory mechanisms, however, are poorly understood. In particular, it is unclear how the activity of single stem cells is coordinated within the population or if cells divide in a purely random fashion. We addressed this issue by analyzing division events in an adult neural stem cell (NSC) population of the zebrafish telencephalon. Spatial statistics and mathematical modeling of over 80,000 NSCs in 36 brain hemispheres revealed weakly aggregated, nonrandom division patterns in space and time. Analyzing divisions at 2 time points allowed us to infer cell cycle and S-phase lengths computationally. Interestingly, we observed rapid cell cycle reentries in roughly 15% of newly born NSCs. In agent-based simulations of NSC populations, this redividing activity sufficed to induce aggregated spatiotemporal division patterns that matched the ones observed experimentally. In contrast, omitting redivisions leads to a random spatiotemporal distribution of dividing cells. Spatiotemporal aggregation of dividing stem cells can thus emerge solely from the cells’ history.

An interdisciplinary study of the rules governing cell divisions in a population of neural stem cells in the zebrafish brain reveals the existence of aggregated spatio-temporal division patterns of rapid cell cycles in stem cells, and shows that these patterns can be explained by a simple agent-based model relying solely on the cells‘ division history.

Bridging from single to collective cell migration: A review of models and links to experiments

Mathematical and computational models can assist in gaining an understanding of cell behavior at many levels of organization. Here, we review models in the literature that focus on eukaryotic cell motility at 3 size scales: intracellular signaling that regulates cell shape and movement, single cell motility, and collective cell behavior from a few cells to tissues. We survey recent literature to summarize distinct computational methods (phase-field, polygonal, Cellular Potts, and spherical cells). We discuss models that bridge between levels of organization, and describe levels of detail, both biochemical and geometric, included in the models. We also highlight links between models and experiments. We find that models that span the 3 levels are still in the minority.