Statistical inference with a manifold-constrained RNA velocity model uncovers cell cycle speed modulations.

Across a range of biological processes, cells undergo coordinated changes in gene expression, resulting in transcriptome dynamics that unfold within a low-dimensional manifold. Single-cell RNA-sequencing (scRNA-seq) only measures temporal snapshots of gene expression. However, information on the underlying low-dimensional dynamics can be extracted using RNA velocity, which models unspliced and spliced RNA abundances to estimate the rate of change of gene expression. Available RNA velocity algorithms can be fragile and rely on heuristics that lack statistical control. Moreover, the estimated vector field is not dynamically consistent with the traversed gene expression manifold. Here, we develop a generative model of RNA velocity and a Bayesian inference approach that solves these problems. Our model couples velocity field and manifold estimation in a reformulated, unified framework, so as to coherently identify the parameters of an autonomous dynamical system. Focusing on the cell cycle, we implemented VeloCycle to study gene regulation dynamics on one-dimensional periodic manifolds and validated using live-imaging its ability to infer actual cell cycle periods. We benchmarked RNA velocity inference with sensitivity analyses and demonstrated one- and multiple-sample testing. We also conducted Markov chain Monte Carlo inference on the model, uncovering key relationships between gene-specific kinetics and our gene-independent velocity estimate. Finally, we applied VeloCycle to in vivo samples and in vitro genome-wide Perturb-seq, revealing regionally-defined proliferation modes in neural progenitors and the effect of gene knockdowns on cell cycle speed. Ultimately, VeloCycle expands the scRNA-seq analysis toolkit with a modular and statistically rigorous RNA velocity inference framework.

Self-organizing in vitro mouse neural tube organoids mimic embryonic development.

The embryonic neural tube is the origin of the entire adult nervous system, and disturbances in its development cause life-threatening birth defects. However, the study of mammalian neural tube development is limited by the lack of physiologically realistic three-dimensional (3D) in vitro models. Here, we report a self-organizing 3D neural tube organoid model derived from single mouse embryonic stem cells that exhibits an in vivo-like tissue architecture, cell type composition and anterior-posterior (AP) patterning. Moreover, maturation of the neural tube organoids showed the emergence of multipotent neural crest cells and mature neurons. Single-cell transcriptome analyses revealed the sequence of transcriptional events in the emergence of neural crest cells and neural differentiation. Thanks to the accessibility of this model, phagocytosis of migrating neural crest cells could be observed in real time for the first time in a mammalian model. We thus introduce a tractable in vitro model to study some of the key morphogenetic and cell type derivation events during early neural development.

Molecular profiling of stem cell-derived retinal pigment epithelial cell differentiation established for clinical translation.

Human embryonic stem cell-derived retinal pigment epithelial cells (hESC-RPE) are a promising cell source to treat age-related macular degeneration (AMD). Despite several ongoing clinical studies, a detailed mapping of transient cellular states during in vitro differentiation has not been performed. Here, we conduct single-cell transcriptomic profiling of an hESC-RPE differentiation protocol that has been developed for clinical use. Differentiation progressed through a culture diversification recapitulating early embryonic development, whereby cells rapidly acquired a rostral embryo patterning signature before converging toward the RPE lineage. At intermediate steps, we identified and examined the potency of an NCAM1+ retinal progenitor population and showed the ability of the protocol to suppress non-RPE fates. We demonstrated that the method produces a pure RPE pool capable of maturing further after subretinal transplantation in a large-eyed animal model. Our evaluation of hESC-RPE differentiation supports the development of safe and efficient pluripotent stem cell-based therapies for AMD.

Sphingolipids control dermal fibroblast heterogeneity.

Human cells produce thousands of lipids that change during cell differentiation and can vary across individual cells of the same type. However, we are only starting to characterize the function of these cell-to-cell differences in lipid composition. Here, we measured the lipidomes and transcriptomes of individual human dermal fibroblasts by coupling high-resolution mass spectrometry imaging with single-cell transcriptomics. We found that the cell-to-cell variations of specific lipid metabolic pathways contribute to the establishment of cell states involved in the organization of skin architecture. Sphingolipid composition is shown to define fibroblast subpopulations, with sphingolipid metabolic rewiring driving cell-state transitions. Therefore, cell-to-cell lipid heterogeneity affects the determination of cell states, adding a new regulatory component to the self-organization of multicellular systems.

Molecular architecture of the developing mouse brain.

The mammalian brain develops through a complex interplay of spatial cues generated by diffusible morphogens, cell-cell interactions and intrinsic genetic programs that result in probably more than a thousand distinct cell types. A complete understanding of this process requires a systematic characterization of cell states over the entire spatiotemporal range of brain development. The ability of single-cell RNA sequencing and spatial transcriptomics to reveal the molecular heterogeneity of complex tissues has therefore been particularly powerful in the nervous system. Previous studies have explored development in specific brain regions1-8, the whole adult brain9 and even entire embryos10. Here we report a comprehensive single-cell transcriptomic atlas of the embryonic mouse brain between gastrulation and birth. We identified almost eight hundred cellular states that describe a developmental program for the functional elements of the brain and its enclosing membranes, including the early neuroepithelium, region-specific secondary organizers, and both neurogenic and gliogenic progenitors. We also used in situ mRNA sequencing to map the spatial expression patterns of key developmental genes. Integrating the in situ data with our single-cell clusters revealed the precise spatial organization of neural progenitors during the patterning of the nervous system.

Structure of ribosome-bound azole-modified peptide phazolicin rationalizes its species-specific mode of bacterial translation inhibition.

Ribosome-synthesized post-translationally modified peptides (RiPPs) represent a rapidly expanding class of natural products with various biological activities. Linear azol(in)e-containing peptides (LAPs) comprise a subclass of RiPPs that display outstanding diversity of mechanisms of action while sharing common structural features. Here, we report the discovery of a new LAP biosynthetic gene cluster in the genome of Rhizobium Pop5, which encodes the precursor peptide and modification machinery of phazolicin (PHZ) - an extensively modified peptide exhibiting narrow-spectrum antibacterial activity against some symbiotic bacteria of leguminous plants. The cryo-EM structure of the Escherichia coli 70S-PHZ complex reveals that the drug interacts with the 23S rRNA and uL4/uL22 proteins and obstructs ribosomal exit tunnel in a way that is distinct from other compounds. We show that the uL4 loop sequence determines the species-specificity of antibiotic action. PHZ expands the known diversity of LAPs and may be used in the future as biocontrol agent for agricultural needs.

Sorting Out Antibiotics' Mechanisms of Action: a Double Fluorescent Protein Reporter for High-Throughput Screening of Ribosome and DNA Biosynthesis Inhibitors.

In order to accelerate drug discovery, a simple, reliable, and cost-effective system for high-throughput identification of a potential antibiotic mechanism of action is required. To facilitate such screening of new antibiotics, we created a double-reporter system for not only antimicrobial activity detection but also simultaneous sorting of potential antimicrobials into those that cause ribosome stalling and those that induce the SOS response due to DNA damage. In this reporter system, the red fluorescent protein gene rfp was placed under the control of the SOS-inducible sulA promoter. The gene of the far-red fluorescent protein, katushka2S, was inserted downstream of the tryptophan attenuator in which two tryptophan codons were replaced by alanine codons, with simultaneous replacement of the complementary part of the attenuator to preserve the ability to form secondary structures that influence transcription termination. This genetically modified attenuator makes possible Katushka2S expression only upon exposure to ribosome-stalling compounds. The application of red and far-red fluorescent proteins provides a high signal-to-background ratio without any need of enzymatic substrates for detection of the reporter activity. This reporter was shown to be efficient in high-throughput screening of both synthetic and natural chemicals.