Data-driven model discovery and model selection for noisy biological systems.

Biological systems exhibit complex dynamics that differential equations can often adeptly represent. Ordinary differential equation models are widespread; until recently their construction has required extensive prior knowledge of the system. Machine learning methods offer alternative means of model construction: differential equation models can be learnt from data via model discovery using sparse identification of nonlinear dynamics (SINDy).

Transition paths across the EMT landscape are dictated by network logic.

During development and cancer metastasis, cells transition reversibly from epithelial to mesenchymal via intermediate cell states during epithelial-mesenchymal transition (EMT). EMT is controlled by gene regulatory networks (GRNs) and can be described by a three-node GRN that permits tristable EMT landscapes. In this GRN, multiple inputs regulate the transcription factor ZEB that induces EMT.

Developmental hematopoietic stem cell variation explains clonal hematopoiesis later in life.

Clonal hematopoiesis becomes increasingly common with age, but its cause is enigmatic because driver mutations are often absent. Serial observations infer weak selection indicating variants are acquired much earlier in life with unexplained initial growth spurts. Here we use fluctuating CpG methylation as a lineage marker to track stem cell clonal dynamics of hematopoiesis.

The RNA exosome maintains cellular RNA homeostasis by controlling transcript abundance in the brain.

Intracellular ribonucleases (RNases) are essential in all aspects of RNA metabolism, including maintaining accurate RNA levels. Inherited mutations in genes encoding ubiquitous RNases are associated with human diseases, primarily affecting the nervous system. Recessive mutations in genes encoding an evolutionarily conserved RNase complex, the RNA exosome, lead to syndromic neurodevelopmental disorders characterized by progressive neurodegeneration, such as Pontocerebellar Hypoplasia Type 1b (PCH1b).

Intratumoral CXCL12 Gradients Contextualize Tumor Cell Invasion, Migration and Immune Suppression in Breast Cancer.

Although the CXCL12/CXCR4 pathway has been prior investigated for its prometastatic and immuno- suppressive roles in the tumor microenvironment, evidence on the spatiotemporal regulation of these hallmarks has been lacking. Here, we demonstrate that CXCL12 forms a gradient specifically around cancer cell intravasation doorways, also known as Tumor Microenvironment of Metastasis (TMEM) doorways, thus facilitating the chemotactic translocation of prometastatic tumor cells expressing CXCR4 toward the perivascular TMEM doorways for subsequent entry into peripheral circulation. Fur- thermore, we demonstrate that the CXCL12-rich micro-environment around TMEM doorways may cre- ate immunosuppressive niches, whereby CD8 T cells, despite being attracted to these regions, often exhibit reduced effector functions, limiting their efficacy.

Gastrulation-stage gene expression in mouse embryos foreshadows the development of syndromic birth defects.

In animal models, deficiency phenocopies gene expression changes and birth defects seen in Cornelia de Lange syndrome, the most common cause of which is haploinsufficiency. Previous studies in mice suggested that heart development is abnormal as soon as cardiogenic tissue is formed. To investigate this, we performed single-cell RNA sequencing on wild-type and mouse embryos at gastrulation and early cardiac crescent stages.

Developmental hematopoietic stem cell variation explains clonal hematopoiesis later in life.

Clonal hematopoiesis becomes increasingly common with age, but its cause is enigmatic because driver mutations are often absent. Serial observations infer weak selection indicating variants are acquired much earlier in life with unexplained initial growth spurts. Here we use fluctuating CpG methylation as a lineage marker to track stem cell clonal dynamics of hematopoiesis.

Gastrulation-stage gene expression in mouse embryos foreshadows the development of syndromic birth defects.

Unlabelled: In animal models, -deficiency phenocopies gene expression changes and birth defects seen in Cornelia de Lange Syndrome (CdLS), the most common cause of which is -haploinsufficiency. Previous studies in mice suggested that heart development is abnormal as soon as cardiogenic tissue is formed. To investigate this, we performed single-cell RNA-sequencing on wildtype (WT) and mouse embryos at gastrulation and early cardiac crescent stages.

Direct androgen receptor control of sexually dimorphic gene expression in the mammalian kidney.

Mammalian organs exhibit distinct physiology, disease susceptibility, and injury responses between the sexes. In the mouse kidney, sexually dimorphic gene activity maps predominantly to proximal tubule (PT) segments. Bulk RNA sequencing (RNA-seq) data demonstrated that sex differences were established from 4 and 8 weeks after birth under gonadal control.

Single-cell and spatial transcriptomics identify a macrophage population associated with skeletal muscle fibrosis.

Macrophages are essential for skeletal muscle homeostasis, but how their dysregulation contributes to the development of fibrosis in muscle disease remains unclear. Here, we used single-cell transcriptomics to determine the molecular attributes of dystrophic and healthy muscle macrophages. We identified six clusters and unexpectedly found that none corresponded to traditional definitions of M1 or M2 macrophages.

Direct androgen receptor regulation of sexually dimorphic gene expression in the mammalian kidney.

Mammalian organs exhibit distinct physiology, disease susceptibility and injury responses between the sexes. In the mouse kidney, sexually dimorphic gene activity maps predominantly to proximal tubule (PT) segments. Bulk RNA-seq data demonstrated sex differences were established from 4 and 8 weeks after birth under gonadal control.

Single-cell and spatial transcriptomics identify a macrophage population associated with skeletal muscle fibrosis.

The monocytic/macrophage system is essential for skeletal muscle homeostasis, but its dysregulation contributes to the pathogenesis of muscle degenerative disorders. Despite our increasing knowledge of the role of macrophages in degenerative disease, it still remains unclear how macrophages contribute to muscle fibrosis. Here, we used single-cell transcriptomics to determine the molecular attributes of dystrophic and healthy muscle macrophages.

Gene regulatory network inference with popInfer reveals dynamic regulation of hematopoietic stem cell quiescence upon diet restriction and aging.

Inference of gene regulatory networks (GRNs) can reveal cell state transitions from single-cell genomics data. However, obstacles to temporal inference from snapshot data are difficult to overcome. Single-nuclei multiomics data offer means to bridge this gap and derive temporal information from snapshot data using joint measurements of gene expression and chromatin accessibility in the same single cells.

Julia for biologists.

Major computational challenges exist in relation to the collection, curation, processing and analysis of large genomic and imaging datasets, as well as the simulation of larger and more realistic models in systems biology. Here we discuss how a relative newcomer among programming languages-Julia-is poised to meet the current and emerging demands in the computational biosciences and beyond. Speed, flexibility, a thriving package ecosystem and readability are major factors that make high-performance computing and data analysis available to an unprecedented degree.

Myeloid-Derived Suppressor-Cell Dynamics Control Outcomes in the Metastatic Niche.

Myeloid-derived suppressor cells (MDSC) play a prominent role in the tumor microenvironment. A quantitative understanding of the tumor-MDSC interactions that influence disease progression is critical, and currently lacking. We developed a mathematical model of metastatic growth and progression in immune-rich tumor microenvironments.

Profiling intermediate cell states in high resolution.

Kong et al. present Capybara, a computational method to identify cell states from single-cell gene expression data. Notably, Capybara can identify intermediate cell states and cell state transitions, offering biologists new means with which to interrogate the states and fates of cells.

Age-dependent effects of Igf2bp2 on gene regulation, function, and aging of hematopoietic stem cells in mice.

Increasing evidence links metabolism, protein synthesis, and growth signaling to impairments in the function of hematopoietic stem and progenitor cells (HSPCs) during aging. The Lin28b/Hmga2 pathway controls tissue development, and the postnatal downregulation of this pathway limits the self-renewal of adult vs fetal hematopoietic stem cells (HSCs). Igf2bp2 is an RNA binding protein downstream of Lin28b/Hmga2, which regulates messenger RNA stability and translation.

Single-cell transcriptomic analysis of zebrafish cranial neural crest reveals spatiotemporal regulation of lineage decisions during development.

Neural crest (NC) cells migrate throughout vertebrate embryos to give rise to a huge variety of cell types, but when and where lineages emerge and their regulation remain unclear. We have performed single-cell RNA sequencing (RNA-seq) of cranial NC cells from the first pharyngeal arch in zebrafish over several stages during migration. Computational analysis combining pseudotime and real-time data reveals that these NC cells first adopt a transitional state, becoming specified mid-migration, with the first lineage decisions being skeletal and pigment, followed by neural and glial progenitors.