Toward a Common Coordinate Framework for the Human Body.

Understanding the genetic and molecular drivers of phenotypic heterogeneity across individuals is central to biology. As new technologies enable fine-grained and spatially resolved molecular profiling, we need new computational approaches to integrate data from the same organ across different individuals into a consistent reference and to construct maps of molecular and cellular organization at histological and anatomical scales. Here, we review previous efforts and discuss challenges involved in establishing such a common coordinate framework, the underlying map of tissues and organs. We focus on strategies to handle anatomical variation across individuals and highlight the need for new technologies and analytical methods spanning multiple hierarchical scales of spatial resolution.

Molecular logic of cellular diversification in the mouse cerebral cortex.

The mammalian cerebral cortex has an unparalleled diversity of cell types, which are generated during development through a series of temporally orchestrated events that are under tight evolutionary constraint and are critical for proper cortical assembly and function1,2. However, the molecular logic that governs the establishment and organization of cortical cell types remains unknown, largely due to the large number of cell classes that undergo dynamic cell-state transitions over extended developmental timelines. Here we generate a comprehensive atlas of the developing mouse neocortex, using single-cell RNA sequencing and single-cell assay for transposase-accessible chromatin using sequencing. We sampled the neocortex every day throughout embryonic corticogenesis and at early postnatal ages, and complemented the sequencing data with a spatial transcriptomics time course. We computationally reconstruct developmental trajectories across the diversity of cortical cell classes, and infer their spatial organization and the gene regulatory programs that accompany their lineage bifurcation decisions and differentiation trajectories. Finally, we demonstrate how this developmental map pinpoints the origin of lineage-specific developmental abnormalities that are linked to aberrant corticogenesis in mutant mice. The data provide a global picture of the regulatory mechanisms that govern cellular diversification in the neocortex.

A transcriptomic and epigenomic cell atlas of the mouse primary motor cortex.

Single-cell transcriptomics can provide quantitative molecular signatures for large, unbiased samples of the diverse cell types in the brain1-3. With the proliferation of multi-omics datasets, a major challenge is to validate and integrate results into a biological understanding of cell-type organization. Here we generated transcriptomes and epigenomes from more than 500,000 individual cells in the mouse primary motor cortex, a structure that has an evolutionarily conserved role in locomotion. We developed computational and statistical methods to integrate multimodal data and quantitatively validate cell-type reproducibility. The resulting reference atlas-containing over 56 neuronal cell types that are highly replicable across analysis methods, sequencing technologies and modalities-is a comprehensive molecular and genomic account of the diverse neuronal and non-neuronal cell types in the mouse primary motor cortex. The atlas includes a population of excitatory neurons that resemble pyramidal cells in layer 4 in other cortical regions4. We further discovered thousands of concordant marker genes and gene regulatory elements for these cell types. Our results highlight the complex molecular regulation of cell types in the brain and will directly enable the design of reagents to target specific cell types in the mouse primary motor cortex for functional analysis.

Deep learning and alignment of spatially resolved single-cell transcriptomes with Tangram.

Charting an organs' biological atlas requires us to spatially resolve the entire single-cell transcriptome, and to relate such cellular features to the anatomical scale. Single-cell and single-nucleus RNA-seq (sc/snRNA-seq) can profile cells comprehensively, but lose spatial information. Spatial transcriptomics allows for spatial measurements, but at lower resolution and with limited sensitivity. Targeted in situ technologies solve both issues, but are limited in gene throughput. To overcome these limitations we present Tangram, a method that aligns sc/snRNA-seq data to various forms of spatial data collected from the same region, including MERFISH, STARmap, smFISH, Spatial Transcriptomics (Visium) and histological images. Tangram can map any type of sc/snRNA-seq data, including multimodal data such as those from SHARE-seq, which we used to reveal spatial patterns of chromatin accessibility. We demonstrate Tangram on healthy mouse brain tissue, by reconstructing a genome-wide anatomically integrated spatial map at single-cell resolution of the visual and somatomotor areas.

Disentangling bacterial invasiveness from lethality in an experimental host-pathogen system.

Quantifying virulence remains a central problem in human health, pest control, disease ecology, and evolutionary biology. Bacterial virulence is typically quantified by the LT50 (i.e., the time taken to kill 50% of infected hosts); however, such an indicator cannot account for the full complexity of the infection process, such as distinguishing between the pathogen's ability to colonize versus kill the hosts. Indeed, the pathogen needs to breach the primary defenses in order to colonize, find a suitable environment to replicate, and finally express the virulence factors that cause disease. Here, we show that two virulence attributes, namely pathogen lethality and invasiveness, can be disentangled from the survival curves of a laboratory population of Caenorhabditis elegans nematodes exposed to three bacterial pathogens: Pseudomonas aeruginosa, Serratia marcescens, and Salmonella enterica We first show that the host population eventually experiences a constant mortality rate, which quantifies the lethality of the pathogen. We then show that the time necessary to reach this constant mortality rate regime depends on the pathogen growth rate and colonization rate, and thus determines the pathogen invasiveness. Our framework reveals that Serratia marcescens is particularly good at the initial colonization of the host, whereas Salmonella enterica is a poor colonizer yet just as lethal once established. Pseudomonas aeruginosa, on the other hand, is both a good colonizer and highly lethal after becoming established. The ability to quantitatively characterize the ability of different pathogens to perform each of these steps has implications for treatment and prevention of disease and for the evolution and ecology of pathogens.

Correction: Resource Availability Modulates the Cooperative and Competitive Nature of a Microbial Cross-Feeding Mutualism.

[This corrects the article DOI: 10.1371/journal.pbio.1002540.].

Resource Availability Modulates the Cooperative and Competitive Nature of a Microbial Cross-Feeding Mutualism.

Mutualisms between species play an important role in ecosystem function and stability. However, in some environments, the competitive aspects of an interaction may dominate the mutualistic aspects. Although these transitions could have far-reaching implications, it has been difficult to study the causes and consequences of this mutualistic-competitive transition in experimentally tractable systems. Here, we study a microbial cross-feeding mutualism in which each yeast strain supplies an essential amino acid for its partner strain. We find that, depending upon the amount of freely available amino acid in the environment, this pair of strains can exhibit an obligatory mutualism, facultative mutualism, competition, parasitism, competitive exclusion, or failed mutualism leading to extinction of the population. A simple model capturing the essential features of this interaction explains how resource availability modulates the interaction and predicts that changes in the dynamics of the mutualism in deteriorating environments can provide advance warning that collapse of the mutualism is imminent. We confirm this prediction experimentally by showing that, in the high nutrient competitive regime, the strains rapidly reach a common carrying capacity before slowly reaching the equilibrium ratio between the strains. However, in the low nutrient regime, before collapse of the obligate mutualism, we find that the ratio rapidly reaches its equilibrium and it is the total abundance that is slow to reach equilibrium. Our results provide a general framework for how mutualisms may transition between qualitatively different regimes of interaction in response to changes in nutrient availability in the environment.

Giant Amplification of Noise in Fluctuation-Induced Pattern Formation.

The amplitude of fluctuation-induced patterns might be expected to be proportional to the strength of the driving noise, suggesting that such patterns would be difficult to observe in nature. Here, we show that a large class of spatially extended dynamical systems driven by intrinsic noise can exhibit giant amplification, yielding patterns whose amplitude is comparable to that of deterministic Turing instabilities. The giant amplification results from the interplay between noise and nonorthogonal eigenvectors of the linear stability matrix, yielding transients that grow with time, and which, when driven by the ever-present intrinsic noise, lead to persistent large amplitude patterns. This mechanism shows that fluctuation-induced Turing patterns are observable, and are not strongly limited by the amplitude of demographic stochasticity nor by the value of the diffusion coefficients.

Universal biology and the statistical mechanics of early life.

All known life on the Earth exhibits at least two non-trivial common features: the canonical genetic code and biological homochirality, both of which emerged prior to the Last Universal Common Ancestor state. This article describes recent efforts to provide a narrative of this epoch using tools from statistical mechanics. During the emergence of self-replicating life far from equilibrium in a period of chemical evolution, minimal models of autocatalysis show that homochirality would have necessarily co-evolved along with the efficiency of early-life self-replicators. Dynamical system models of the evolution of the genetic code must explain its universality and its highly refined error-minimization properties. These have both been accounted for in a scenario where life arose from a collective, networked phase where there was no notion of species and perhaps even individuality itself. We show how this phase ultimately terminated during an event sometimes known as the Darwinian transition, leading to the present epoch of tree-like vertical descent of organismal lineages. These examples illustrate concrete examples of universal biology: the quest for a fundamental understanding of the basic properties of living systems, independent of precise instantiation in chemistry or other media.This article is part of the themed issue 'Reconceptualizing the origins of life'.

Genetic Toggle Switch in the Absence of Cooperative Binding: Exact Results.

We present an analytical treatment of a genetic switch model consisting of two mutually inhibiting genes operating without cooperative binding of the corresponding transcription factors. Previous studies have numerically shown that these systems can exhibit bimodal dynamics without possessing two stable fixed points at the deterministic level. We analytically show that bimodality is induced by the noise and find the critical repression strength that controls a transition between the bimodal and nonbimodal regimes. We also identify characteristic polynomial scaling laws of the mean switching time between bimodal states. These results, independent of the model under study, reveal essential differences between these systems and systems with cooperative binding, where there is no critical threshold for bimodality and the mean switching time scales exponentially with the system size.

Noise-Induced Mechanism for Biological Homochirality of Early Life Self-Replicators.

The observed single-handedness of biological amino acids and sugars has long been attributed to autocatalysis. However, the stability of homochiral states in deterministic autocatalytic systems relies on cross inhibition of the two chiral states, an unlikely scenario for early life self-replicators. Here, we present a theory for a stochastic individual-level model of autocatalysis due to early life self-replicators. Without chiral inhibition, the racemic state is the global attractor of the deterministic dynamics, but intrinsic multiplicative noise stabilizes the homochiral states, in both well-mixed and spatially extended systems. We conclude that autocatalysis is a viable mechanism for homochirality, without imposing additional nonlinearities such as chiral inhibition.

Noise-induced bistable states and their mean switching time in foraging colonies.

We investigate a type of bistability occurring in population systems where noise not only causes transitions between stable states, but also constructs the states themselves. We focus on the experimentally well-studied system of ants choosing between two food sources to illustrate the essential points, but the ideas are more general. The mean time for switching between the two bistable states of the system is calculated. This suggests a procedure for estimating, in a real system, the critical population size above which bistability ceases to occur.

Stochastic pattern formation and spontaneous polarisation: the linear noise approximation and beyond.

We review the mathematical formalism underlying the modelling of stochasticity in biological systems. Beginning with a description of the system in terms of its basic constituents, we derive the mesoscopic equations governing the dynamics which generalise the more familiar macroscopic equations. We apply this formalism to the analysis of two specific noise-induced phenomena observed in biologically inspired models. In the first example, we show how the stochastic amplification of a Turing instability gives rise to spatial and temporal patterns which may be understood within the linear noise approximation. The second example concerns the spontaneous emergence of cell polarity, where we make analytic progress by exploiting a separation of time-scales.

Gut microbial species and endotypes associate with remission in ulcerative colitis patients treated with anti-TNF or anti-integrin therapy.

BACKGROUND AND AIMS: The gut microbiota contributes to aberrant inflammation in inflammatory bowel disease, but the bacterial factors causing or exacerbating inflammation are not fully understood. Further, the predictive or prognostic value of gut microbial biomarkers for remission in response to biologic therapy is unclear. METHODS: We perform whole metagenomic sequencing of 550 stool samples from 287 ulcerative colitis patients from a large phase 3 head-to-head study of infliximab and etrolizumab. RESULTS: We identify several bacterial species in baseline and/or post-treatment samples that associate with clinical remission. These include previously described associations (Faecalibacterium prausnitzii_F) as well as new associations with remission to biologic therapy (Flavonifractor plautii). We build multivariate models and find that gut microbial species are better predictors for remission than clinical variables alone. Finally, we describe patient groups that differ in microbiome composition and remission rate after induction therapy, suggesting the potential utility of microbiome-based endotyping. CONCLUSIONS: In this large study of ulcerative colitis patients, we show that few individual species associate strongly with clinical remission, but multivariate models including microbiome can predict clinical remission and have better predictive power compared to clinical data alone.

Prediction of single-cell RNA expression profiles in live cells by Raman microscopy with Raman2RNA.

Single-cell RNA sequencing and other profiling assays have helped interrogate cells at unprecedented resolution and scale, but are inherently destructive. Raman microscopy reports on the vibrational energy levels of proteins and metabolites in a label-free and nondestructive manner at subcellular spatial resolution, but it lacks genetic and molecular interpretability. Here we present Raman2RNA (R2R), a method to infer single-cell expression profiles in live cells through label-free hyperspectral Raman microscopy images and domain translation. We predict single-cell RNA sequencing profiles nondestructively from Raman images using either anchor-based integration with single molecule fluorescence in situ hybridization, or anchor-free generation with adversarial autoencoders. R2R outperformed inference from brightfield images (cosine similarities: R2R >0.85 and brightfield <0.15). In reprogramming of mouse fibroblasts into induced pluripotent stem cells, R2R inferred the expression profiles of various cell states. With live-cell tracking of mouse embryonic stem cell differentiation, R2R traced the early emergence of lineage divergence and differentiation trajectories, overcoming discontinuities in expression space. R2R lays a foundation for future exploration of live genomic dynamics.

Reference-based cell type matching of in situ image-based spatial transcriptomics data on primary visual cortex of mouse brain.

With the advent of multiplex fluorescence in situ hybridization (FISH) and in situ RNA sequencing technologies, spatial transcriptomics analysis is advancing rapidly, providing spatial location and gene expression information about cells in tissue sections at single cell resolution. Cell type classification of these spatially-resolved cells can be inferred by matching the spatial transcriptomics data to reference atlases derived from single cell RNA-sequencing (scRNA-seq) in which cell types are defined by differences in their gene expression profiles. However, robust cell type matching of the spatially-resolved cells to reference scRNA-seq atlases is challenging due to the intrinsic differences in resolution between the spatial and scRNA-seq data. In this study, we systematically evaluated six computational algorithms for cell type matching across four image-based spatial transcriptomics experimental protocols (MERFISH, smFISH, BaristaSeq, and ExSeq) conducted on the same mouse primary visual cortex (VISp) brain region. We find that many cells are assigned as the same type by multiple cell type matching algorithms and are present in spatial patterns previously reported from scRNA-seq studies in VISp. Furthermore, by combining the results of individual matching strategies into consensus cell type assignments, we see even greater alignment with biological expectations. We present two ensemble meta-analysis strategies used in this study and share the consensus cell type matching results in the Cytosplore Viewer ( https://viewer.cytosplore.org ) for interactive visualization and data exploration. The consensus matching can also guide spatial data analysis using SSAM, allowing segmentation-free cell type assignment.

Genetic vulnerability to Crohn's disease reveals a spatially resolved epithelial restitution program.

Effective tissue repair requires coordinated intercellular communication to sense damage, remodel the tissue, and restore function. Here, we dissected the healing response in the intestinal mucosa by mapping intercellular communication at single-cell resolution and integrating with spatial transcriptomics. We demonstrated that a risk variant for Crohn's disease, hepatocyte growth factor activator (HGFAC) Arg509His (R509H), disrupted a damage-sensing pathway connecting the coagulation cascade to growth factors that drive the differentiation of wound-associated epithelial (WAE) cells and production of a localized retinoic acid (RA) gradient to promote fibroblast-mediated tissue remodeling. Specifically, we showed that HGFAC R509H was activated by thrombin protease activity but exhibited impaired proteolytic activation of the growth factor macrophage-stimulating protein (MSP). In Hgfac R509H mice, reduced MSP activation in response to wounding of the colon resulted in impaired WAE cell induction and delayed healing. Through integration of single-cell transcriptomics and spatial transcriptomics, we demonstrated that WAE cells generated RA in a spatially restricted region of the wound site and that mucosal fibroblasts responded to this signal by producing extracellular matrix and growth factors. We further dissected this WAE cell-fibroblast signaling circuit in vitro using a genetically tractable organoid coculture model. Collectively, these studies exploited a genetic perturbation associated with human disease to disrupt a fundamental biological process and then reconstructed a spatially resolved mechanistic model of tissue healing.

Stepwise-edited, human melanoma models reveal mutations' effect on tumor and microenvironment.

Establishing causal relationships between genetic alterations of human cancers and specific phenotypes of malignancy remains a challenge. We sequentially introduced mutations into healthy human melanocytes in up to five genes spanning six commonly disrupted melanoma pathways, forming nine genetically distinct cellular models of melanoma. We connected mutant melanocyte genotypes to malignant cell expression programs in vitro and in vivo, replicative immortality, malignancy, rapid tumor growth, pigmentation, metastasis, and histopathology. Mutations in malignant cells also affected tumor microenvironment composition and cell states. Our melanoma models shared genotype-associated expression programs with patient melanomas, and a deep learning model showed that these models partially recapitulated genotype-associated histopathological features as well. Thus, a progressive series of genome-edited human cancer models can causally connect genotypes carrying multiple mutations to phenotype.

Epigenetic encoding, heritability and plasticity of glioma transcriptional cell states.

Single-cell RNA sequencing has revealed extensive transcriptional cell state diversity in cancer, often observed independently of genetic heterogeneity, raising the central question of how malignant cell states are encoded epigenetically. To address this, here we performed multiomics single-cell profiling-integrating DNA methylation, transcriptome and genotype within the same cells-of diffuse gliomas, tumors characterized by defined transcriptional cell state diversity. Direct comparison of the epigenetic profiles of distinct cell states revealed key switches for state transitions recapitulating neurodevelopmental trajectories and highlighted dysregulated epigenetic mechanisms underlying gliomagenesis. We further developed a quantitative framework to directly measure cell state heritability and transition dynamics based on high-resolution lineage trees in human samples. We demonstrated heritability of malignant cell states, with key differences in hierarchal and plastic cell state architectures in IDH-mutant glioma versus IDH-wild-type glioblastoma, respectively. This work provides a framework anchoring transcriptional cancer cell states in their epigenetic encoding, inheritance and transition dynamics.