Publications

2024

1. Spotiflow: accurate and efficient spot detection for fluorescence microscopy with deep stereographic flow regression

   Albert Dominguez Mantes ,  Antonio Herrera ,  Irina Khven , and 14 more authors

   bioRxiv, 2024

   Abs PDF

   Identifying spot-like structures in large and noisy microscopy images is a crucial step to produce high quality results in various life-science applications. Imaging-based spatial transcriptomics (iST) methods, in particular, critically depend on the precise detection of millions of transcripts in images with low signal-to-noise ratio. Despite advances in computer vision that have revolutionized many biological imaging tasks, currently adopted spot detection techniques are mostly still based on classical signal processing methods that are fragile and require tedious manual tuning per dataset. In this work, we introduce Spotiflow, a deep learning method that achieves subpixel-accurate localizations by formulating the spot detection task as a multiscale heatmap and stereographic flow regression problem. Spotiflow can be used for 2D images and 3D volumetric stacks and can be trained to generalize across different imaging conditions, tissue types and chemical preparations, while being substantially more time- and memory-e cient than existing methods. We show the eficacy of Spotiflow via extensive quantitative experiments on a variety of diverse datasets and demonstrate that the enhanced accuracy of Spotiflow leads to meaningful improvements in the biological insights obtained from iST and live imaging experiments. Spotiflow is available as an easy-to-use Python library as well as a napari plugin at https://github.com/weigertlab/spotiflow.

2. Unified Mass Imaging Maps the Lipidome of Vertebrate Development

   Halima Hannah Schede ,  Leila Haj Abdullah Alieh ,  Laurel Ann Rohde , and 12 more authors

   bioRxiv, 2024

   Abs

   Embryo development entails the formation of anatomical structures with distinct biochemical compositions. Compared with the wealth of knowledge on gene regulation, our understanding of metabolic programs operating during embryogenesis is limited. Mass spectrometry imaging (MSI) has the potential to map the distribution of metabolites across embryo development. Here, we established an analytical framework for the joint analysis of large MSI datasets that allows for the construction of multi-dimensional metabolomic atlases. Employing this framework, we mapped the 4D distribution of over a hundred lipids at quasi-single-cell resolution in Danio rerio embryos. We discovered metabolic trajectories that unfold in concert with morphogenesis and revealed spatially organized biochemical coordination overlooked by bulk measurements. Interestingly, lipid mapping revealed unexpected distributions of sphingolipid and triglyceride species, suggesting their involvement in pattern establishment and organ development. Our approach empowers a new generation of whole-organism metabolomic atlases and enables the discovery of spatially organized metabolic circuits.Competing Interest StatementThe authors have declared no competing interest.

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

   Alex R. Lederer ,  Maxine Leonardi ,  Lorenzo Talamanca , and 11 more authors

   Nature Methods, 2024

   Abs

   Across biological systems, cells undergo coordinated changes in gene expression, resulting in transcriptome dynamics that unfold within a low-dimensional manifold. While low-dimensional dynamics can be extracted using RNA velocity, these algorithms can be fragile and rely on heuristics lacking statistical control. Moreover, the estimated vector field is not dynamically consistent with the traversed gene expression manifold. To address these challenges, we introduce a Bayesian model of RNA velocity that couples velocity field and manifold estimation in a reformulated, unified framework, identifying the parameters of an explicit dynamical system. Focusing on the cell cycle, we implement VeloCycle to study gene regulation dynamics on one-dimensional periodic manifolds and validate its ability to infer cell cycle periods using live imaging. We also apply VeloCycle to reveal speed differences in regionally defined progenitors and Perturb-seq gene knockdowns. Overall, VeloCycle expands the single-cell RNA sequencing analysis toolkit with a modular and statistically consistent RNA velocity inference framework.

4. Specialized signaling centers direct cell fate and spatial organization in a limb organoid model

   Evangelia Skoufa ,  Jixing Zhong ,  Oliver Kahre , and 12 more authors

   bioRxiv, 2024

   Abs

   Specialized signaling centers orchestrate robust development and regeneration. Limb morphogenesis, for instance, requires interactions between the mesoderm and the signaling center apical-ectodermal ridge (AER), whose properties and role in cell fate decisions have remained challenging to dissect. To tackle this, we developed mouse embryonic stem cells (mESCs)-based heterogeneous cultures and a limb organoid model, termed budoids, comprising cells with AER, surface ectoderm, and mesoderm properties. mESCs were first induced into heterogeneous cultures that self-organized into domes in 2D. Aggregating these cultures resulted in formation of limb bud-like structures in 3D, exhibiting chondrogenesis-based symmetry breaking and elongation. Using our organoids and quantitative in situ expression profiling, we uncovered that AER-like cells support nearby limb mesoderm and fibroblast identities while enhancing tissue polarization that permits distant cartilage formation. Together, our findings provide a powerful model to study aspects of limb morphogenesis, and reveal the ability of signaling center AER cells to concurrently modulate cell fate and spatial organization.

2023

1. Neural ADMIXTURE for rapid genomic clustering

   Albert Dominguez Mantes ,  Daniel Mas Montserrat ,  Carlos D. Bustamante , and 2 more authors

   Nature Computational Science, 2023

   Abs PDF

   Characterizing the genetic structure of large cohorts has become increasingly important as genetic studies extend to massive, increasingly diverse biobanks. Popular methods decompose individual genomes into fractional cluster assignments with each cluster representing a vector of DNA variant frequencies. However, with rapidly increasing biobank sizes, these methods have become computationally intractable. Here we present Neural ADMIXTURE, a neural network autoencoder that follows the same modeling assumptions as the current standard algorithm, ADMIXTURE, while reducing the compute time by orders of magnitude surpassing even the fastest alternatives. One month of continuous compute using ADMIXTURE can be reduced to just hours with Neural ADMIXTURE. A multi-head approach allows Neural ADMIXTURE to offer even further acceleration by computing multiple cluster numbers in a single run. Furthermore, the models can be stored, allowing cluster assignment to be performed on new data in linear time without needing to share the training samples.

2022

1. Archetypal Analysis for population genetics

   Julia Gimbernat-Mayol ,  Albert Dominguez Mantes ,  Carlos D. Bustamante , and 2 more authors

   PLOS Computational Biology, Aug 2022

   Abs

   The estimation of genetic clusters using genomic data has application from genome-wide association studies (GWAS) to demographic history to polygenic risk scores (PRS) and is expected to play an important role in the analyses of increasingly diverse, large-scale cohorts. However, existing methods are computationally-intensive, prohibitively so in the case of nationwide biobanks. Here we explore Archetypal Analysis as an efficient, unsupervised approach for identifying genetic clusters and for associating individuals with them. Such unsupervised approaches help avoid conflating socially constructed ethnic labels with genetic clusters by eliminating the need for exogenous training labels. We show that Archetypal Analysis yields similar cluster structure to existing unsupervised methods such as ADMIXTURE and provides interpretative advantages. More importantly, we show that since Archetypal Analysis can be used with lower-dimensional representations of genetic data, significant reductions in computational time and memory requirements are possible. When Archetypal Analysis is run in such a fashion, it takes several orders of magnitude less compute time than the current standard, ADMIXTURE. Finally, we demonstrate uses ranging across datasets from humans to canids.