A new Neu-a syngeneic model of spontaneously metastatic HER2-positive breast cancer.

Metastatic disease results from the dissemination of tumor cells beyond their organ of origin to grow in distant organs and is the primary cause of death in patients with advanced breast cancer. Preclinical murine models in which primary tumors spontaneously metastasize are valuable tools for studying metastatic progression and novel cancer treatment combinations. Here, we characterize a novel syngeneic murine breast tumor cell line that provides a model of spontaneously metastatic neu-expressing breast cancer with quicker onset of widespread metastases after orthotopic mammary implantation in immune-competent NeuN mice. The NT2.5-lung metastasis (-LM) cell line was derived from serial passaging of tumor cells that were macro-dissected from spontaneous lung metastases after orthotopic mammary implantation of parental NT2.5 cells. Within one week of NT2.5-LM implantation, metastases are observed in the lungs. Within four weeks, metastases are also observed in the bones, spleen, colon, and liver. We demonstrate that NT2.5-LM metastases are positive for NeuN-the murine equivalent of human epidermal growth factor 2 (HER2). We further demonstrate altered expression of markers of epithelial-to-mesenchymal transition (EMT), suggestive of their enhanced metastatic potential. Genomic analyses support these findings and reveal enrichment in EMT-regulating pathways. In addition, the metastases are rapidly growing, proliferative, and responsive to HER2-directed therapy. The new NT2.5-LM model provides certain advantages over the parental NT2/NT2.5 model, given its more rapid and spontaneous development of metastases. Besides investigating mechanisms of metastatic progression, this new model may be used for the rationalized development of novel therapeutic interventions and assessment of therapeutic responses.

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

In animal models, Nipbl deficiency phenocopies gene expression changes and birth defects seen in Cornelia de Lange syndrome, the most common cause of which is Nipbl haploinsufficiency. Previous studies in Nipbl+/- 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 Nipbl+/- mouse embryos at gastrulation and early cardiac crescent stages. Nipbl+/- embryos had fewer mesoderm cells than wild-type and altered proportions of mesodermal cell subpopulations. These findings were associated with underexpression of genes implicated in driving specific mesodermal lineages. In addition, Nanog was found to be overexpressed in all germ layers, and many gene expression changes observed in Nipbl+/- embryos could be attributed to Nanog overexpression. These findings establish a link between Nipbl deficiency, Nanog overexpression, and gene expression dysregulation/lineage misallocation, which ultimately manifest as birth defects in Nipbl+/- animals and Cornelia de Lange syndrome.

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. We show, via the shared prenatal circulation of monozygotic twins, that weak selection conferred by stem cell variation created before birth can reliably yield clonal hematopoiesis later in life. Theory indicates weak selection will lead to dominance given enough time and large enough population sizes. Human hematopoiesis satisfies both these conditions. Stochastic loss of weakly selected variants is naturally prevented by the expansion of stem cell lineages during development. The dominance of stem cell clones created before birth is supported by blood fluctuating CpG methylation patterns that exhibit low correlation between unrelated individuals but are highly correlated between many elderly monozygotic twins. Therefore, clonal hematopoiesis driven by weak selection in later life appears to reflect variation created before birth.

Mimicking the breast metastatic microenvironment: characterization of a novel syngeneic model of HER2+ breast cancer.

Preclinical murine models in which primary tumors spontaneously metastasize to distant organs are valuable tools to study metastatic progression and novel cancer treatment combinations. Here, we characterize a novel syngeneic murine breast tumor cell line, NT2.5-lung metastasis (-LM), that provides a model of spontaneously metastatic neu-expressing breast cancer with quicker onset of widespread metastases after orthotopic mammary implantation in immune-competent NeuN mice. Within one week of orthotopic implantation of NT2.5-LM in NeuN mice, distant metastases can be observed in the lungs. Within four weeks, metastases are also observed in the bones, spleen, colon, and liver. Metastases are rapidly growing, proliferative, and responsive to HER2-directed therapy. We demonstrate altered expression of markers of epithelial-to-mesenchymal transition (EMT) and enrichment in EMT-regulating pathways, suggestive of their enhanced metastatic potential. The new NT2.5-LM model provides more rapid and spontaneous development of widespread metastases. Besides investigating mechanisms of metastatic progression, this new model may be used for the rationalized development of novel therapeutic interventions and assessment of therapeutic responses targeting distant visceral metastases.

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

In animal models, Nipbl -deficiency phenocopies gene expression changes and birth defects seen in Cornelia de Lange Syndrome (CdLS), the most common cause of which is Nipbl -haploinsufficiency. Previous studies in Nipbl +/- 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 Nipbl +/- mouse embryos at gastrulation and early cardiac crescent stages. Nipbl +/- embryos had fewer mesoderm cells than WT and altered proportions of mesodermal cell subpopulations. These findings were associated with underexpression of genes implicated in driving specific mesodermal lineages. In addition, Nanog was found to be overexpressed in all germ layers, and many gene expression changes observed in Nipbl +/- embryos could be attributed to Nanog overexpression. These findings establish a link between Nipbl -deficiency, Nanog overexpression, and gene expression dysregulation/lineage misallocation, which ultimately manifest as birth defects in Nipbl +/- animals and CdLS. Teaser: Gene expression changes during gastrulation of Nipbl -deficient mice shed light on early origins of structural birth defects.

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. Hormone injection studies and genetic removal of androgen and estrogen receptors demonstrated androgen receptor (AR)-mediated regulation of gene activity in PT cells as the regulatory mechanism. Interestingly, caloric restriction feminizes the male kidney. Single-nuclear multiomic analysis identified putative cis-regulatory regions and cooperating factors mediating PT responses to AR activity in the mouse kidney. In the human kidney, a limited set of genes showed conserved sex-linked regulation, whereas analysis of the mouse liver underscored organ-specific differences in the regulation of sexually dimorphic gene expression. These findings raise interesting questions on the evolution, physiological significance, disease, and metabolic linkage of sexually dimorphic gene activity.

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. Rather, the predominant macrophage signature in dystrophic muscle was characterized by high expression of fibrotic factors, galectin-3 (gal-3) and osteopontin (Spp1). Spatial transcriptomics, computational inferences of intercellular communication, and in vitro assays indicated that macrophage-derived Spp1 regulates stromal progenitor differentiation. Gal-3+ macrophages were chronically activated in dystrophic muscle, and adoptive transfer assays showed that the gal-3+ phenotype was the dominant molecular program induced within the dystrophic milieu. Gal-3+ macrophages were also elevated in multiple human myopathies. These studies advance our understanding of macrophages in muscular dystrophy by defining their transcriptional programs and reveal Spp1 as a major regulator of macrophage and stromal progenitor interactions.

Single-cell Ca2+ parameter inference reveals how transcriptional states inform dynamic cell responses.

Single-cell genomic technologies offer vast new resources with which to study cells, but their potential to inform parameter inference of cell dynamics has yet to be fully realized. Here we develop methods for Bayesian parameter inference with data that jointly measure gene expression and Ca2+ dynamics in single cells. We propose to share information between cells via transfer learning: for a sequence of cells, the posterior distribution of one cell is used to inform the prior distribution of the next. In application to intracellular Ca2+ signalling dynamics, we fit the parameters of a dynamical model for thousands of cells with variable single-cell responses. We show that transfer learning accelerates inference with sequences of cells regardless of how the cells are ordered. However, only by ordering cells based on their transcriptional similarity can we distinguish Ca2+ dynamic profiles and associated marker genes from the posterior distributions. Inference results reveal complex and competing sources of cell heterogeneity: parameter covariation can diverge between the intracellular and intercellular contexts. Overall, we discuss the extent to which single-cell parameter inference informed by transcriptional similarity can quantify relationships between gene expression states and signalling dynamics in single cells.

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. We developed popInfer to infer networks that characterize lineage-specific dynamic cell state transitions from joint gene expression and chromatin accessibility data. Benchmarking against alternative methods for GRN inference, we showed that popInfer achieves higher accuracy in the GRNs inferred. popInfer was applied to study single-cell multiomics data characterizing hematopoietic stem cells (HSCs) and the transition from HSC to a multipotent progenitor cell state during murine hematopoiesis across age and dietary conditions. From networks predicted by popInfer, we discovered gene interactions controlling entry to/exit from HSC quiescence that are perturbed in response to diet or aging.

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. We identified six novel clusters. Unexpectedly, none corresponded to traditional definitions of M1 or M2 macrophage activation. Rather, the predominant macrophage signature in dystrophic muscle was characterized by high expression of fibrotic factors, galectin-3 and spp1. Spatial transcriptomics and computational inferences of intercellular communication indicated that spp1 regulates stromal progenitor and macrophage interactions during muscular dystrophy. Galectin-3 + macrophages were chronically activated in dystrophic muscle and adoptive transfer assays showed that the galectin-3 + phenotype was the dominant molecular program induced within the dystrophic milieu. Histological examination of human muscle biopsies revealed that galectin-3 + macrophages were also elevated in multiple myopathies. These studies advance our understanding of macrophages in muscular dystrophy by defining the transcriptional programs induced in muscle macrophages, and reveal spp1 as a major regulator of macrophage and stromal progenitor interactions.

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. Hormone injection studies and genetic removal of androgen and estrogen receptors demonstrated androgen receptor (AR) mediated regulation of gene activity in PT cells as the regulatory mechanism. Interestingly, caloric restriction feminizes the male kidney. Single-nuclear multiomic analysis identified putative cis-regulatory regions and cooperating factors mediating PT responses to AR activity in the mouse kidney. In the human kidney, a limited set of genes showed conserved sex-linked regulation while analysis of the mouse liver underscored organ-specific differences in the regulation of sexually dimorphic gene expression. These findings raise interesting questions on the evolution, physiological significance, and disease and metabolic linkage, of sexually dimorphic gene activity.

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. We highlight how Julia's design is already enabling new ways of analyzing biological data and systems, and we provide a list of resources that can facilitate the transition into Julian computing.

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. We modeled the tumor-immune dynamics with stochastic delay differential equations and studied the impact of delays in MDSC activation/recruitment on tumor growth outcomes. In the lung environment, when the circulating level of MDSCs was low, the MDSC delay had a pronounced impact on the probability of new metastatic establishment: blocking MDSC recruitment could reduce the probability of metastasis by as much as 50%. To predict patient-specific MDSC responses, we fit to the model individual tumors treated with immune checkpoint inhibitors via Bayesian parameter inference. We reveal that control of the inhibition rate of natural killer (NK) cells by MDSCs had a larger influence on tumor outcomes than controlling the tumor growth rate directly. Posterior classification of tumor outcomes demonstrates that incorporating knowledge of the MDSC responses improved predictive accuracy from 63% to 82%. Investigation of the MDSC dynamics in an environment low in NK cells and abundant in cytotoxic T cells revealed, in contrast, that small MDSC delays no longer impacted metastatic growth dynamics. Our results illustrate the importance of MDSC dynamics in the tumor microenvironment overall and predict interventions promoting shifts toward less immune-suppressed states. We propose that there is a pressing need to consider MDSCs more often in analyses of 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. The role of Igf2bp2 in HSC aging is unknown. In this study, an analysis of wild-type and Igf2bp2 knockout mice showed that Igf2bp2 regulates oxidative metabolism in HSPCs and the expression of metabolism, protein synthesis, and stemness-related genes in HSCs of young mice. Interestingly, Igf2bp2 expression and function strongly declined in aging HSCs. In young mice, Igf2bp2 deletion mimicked aging-related changes in HSCs, including changes in Igf2bp2 target gene expression and impairment of colony formation and repopulation capacity. In aged mice, Igf2bp2 gene status had no effect on these parameters in HSCs. Unexpectedly, Igf2bp2-deficient mice exhibited an amelioration of the aging-associated increase in HSCs and myeloid-skewed differentiation. The results suggest that Igf2bp2 controls mitochondrial metabolism, protein synthesis, growth, and stemness of young HSCs, which is necessary for full HSC function during young adult age. However, Igf2bp2 gene function is lost during aging, and it appears to contribute to HSC aging in 2 ways: the aging-related loss of Igf2bp2 gene function impairs the growth and repopulation capacity of aging HSCs, and the activity of Igf2bp2 at a young age contributes to aging-associated HSC expansion and myeloid skewing.

A single-cell resolved cell-cell communication model explains lineage commitment in hematopoiesis.

Cells do not make fate decisions independently. Arguably, every cell-fate decision occurs in response to environmental signals. In many cases, cell-cell communication alters the dynamics of the internal gene regulatory network of a cell to initiate cell-fate transitions, yet models rarely take this into account. Here, we have developed a multiscale perspective to study the granulocyte-monocyte versus megakaryocyte-erythrocyte fate decisions. This transition is dictated by the GATA1-PU.1 network: a classical example of a bistable cell-fate system. We show that, for a wide range of cell communication topologies, even subtle changes in signaling can have pronounced effects on cell-fate decisions. We go on to show how cell-cell coupling through signaling can spontaneously break the symmetry of a homogenous cell population. Noise, both intrinsic and extrinsic, shapes the decision landscape profoundly, and affects the transcriptional dynamics underlying this important hematopoietic cell-fate decision-making system. This article has an associated 'The people behind the papers' interview.

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. In addition, by computationally integrating these data with RNA-seq data from a transgenic Wnt reporter line, we identify gene cohorts with similar temporal responses to Wnts during migration and show that one, Atp6ap2, is required for melanocyte differentiation. Together, our results show that cranial NC cell lineages arise progressively and uncover a series of spatially restricted cell interactions likely to regulate such cell-fate decisions.