Controlling organoid symmetry breaking uncovers an excitable system underlying human axial elongation.

The human embryo breaks symmetry to form the anterior-posterior axis of the body. As the embryo elongates along this axis, progenitors in the tail bud give rise to tissues that generate spinal cord, skeleton, and musculature. This raises the question of how the embryo achieves axial elongation and patterning. While ethics necessitate in vitro studies, the variability of organoid systems has hindered mechanistic insights. Here, we developed a bioengineering and machine learning framework that optimizes organoid symmetry breaking by tuning their spatial coupling. This framework enabled reproducible generation of axially elongating organoids, each possessing a tail bud and neural tube. We discovered that an excitable system composed of WNT/FGF signaling drives elongation by inducing a neuromesodermal progenitor-like signaling center. We discovered that instabilities in the excitable system are suppressed by secreted WNT inhibitors. Absence of these inhibitors led to ectopic tail buds and branches. Our results identify mechanisms governing stable human axial elongation.

Multimodal charting of molecular and functional cell states via in situ electro-sequencing.

Paired mapping of single-cell gene expression and electrophysiology is essential to understand gene-to-function relationships in electrogenic tissues. Here, we developed in situ electro-sequencing (electro-seq) that combines flexible bioelectronics with in situ RNA sequencing to stably map millisecond-timescale electrical activity and profile single-cell gene expression from the same cells across intact biological networks, including cardiac and neural patches. When applied to human-induced pluripotent stem-cell-derived cardiomyocyte patches, in situ electro-seq enabled multimodal in situ analysis of cardiomyocyte electrophysiology and gene expression at the cellular level, jointly defining cell states and developmental trajectories. Using machine-learning-based cross-modal analysis, in situ electro-seq identified gene-to-electrophysiology relationships throughout cardiomyocyte development and accurately reconstructed the evolution of gene expression profiles based on long-term stable electrical measurements. In situ electro-seq could be applicable to create spatiotemporal multimodal maps in electrogenic tissues, potentiating the discovery of cell types and gene programs responsible for electrophysiological function and dysfunction.

Spatial atlas of the mouse central nervous system at molecular resolution.

Spatially charting molecular cell types at single-cell resolution across the 3D volume is critical for illustrating the molecular basis of brain anatomy and functions. Single-cell RNA sequencing has profiled molecular cell types in the mouse brain1,2, but cannot capture their spatial organization. Here we used an in situ sequencing method, STARmap PLUS3,4, to profile 1,022 genes in 3D at a voxel size of 194 × 194 × 345 nm3, mapping 1.09 million high-quality cells across the adult mouse brain and spinal cord. We developed computational pipelines to segment, cluster and annotate 230 molecular cell types by single-cell gene expression and 106 molecular tissue regions by spatial niche gene expression. Joint analysis of molecular cell types and molecular tissue regions enabled a systematic molecular spatial cell-type nomenclature and identification of tissue architectures that were undefined in established brain anatomy. To create a transcriptome-wide spatial atlas, we integrated STARmap PLUS measurements with a published single-cell RNA-sequencing atlas1, imputing single-cell expression profiles of 11,844 genes. Finally, we delineated viral tropisms of a brain-wide transgene delivery tool, AAV-PHP.eB5,6. Together, this annotated dataset provides a single-cell resource that integrates the molecular spatial atlas, brain anatomy and the accessibility to genetic manipulation of the mammalian central nervous system.

Spatiotemporally resolved transcriptomics reveals the subcellular RNA kinetic landscape.

Spatiotemporal regulation of the cellular transcriptome is crucial for proper protein expression and cellular function. However, the intricate subcellular dynamics of RNA remain obscured due to the limitations of existing transcriptomics methods. Here, we report TEMPOmap-a method that uncovers subcellular RNA profiles across time and space at the single-cell level. TEMPOmap integrates pulse-chase metabolic labeling with highly multiplexed three-dimensional in situ sequencing to simultaneously profile the age and location of individual RNA molecules. Using TEMPOmap, we constructed the subcellular RNA kinetic landscape in various human cells from transcription and translocation to degradation. Clustering analysis of RNA kinetic parameters across single cells revealed 'kinetic gene clusters' whose expression patterns were shaped by multistep kinetic sculpting. Importantly, these kinetic gene clusters are functionally segregated, suggesting that subcellular RNA kinetics are differentially regulated in a cell-state- and cell-type-dependent manner. Spatiotemporally resolved transcriptomics provides a gateway to uncovering new spatiotemporal gene regulation principles.

Cyborg Organoids: Implantation of Nanoelectronics via Organogenesis for Tissue-Wide Electrophysiology.

Tissue-wide electrophysiology with single-cell and millisecond spatiotemporal resolution is critical for heart and brain studies. Issues arise, however, from the invasive, localized implantation of electronics that destroys well-connected cellular networks within matured organs. Here, we report the creation of cyborg organoids: the three-dimensional (3D) assembly of soft, stretchable mesh nanoelectronics across the entire organoid by the cell-cell attraction forces from 2D-to-3D tissue reconfiguration during organogenesis. We demonstrate that stretchable mesh nanoelectronics can migrate with and grow into the initial 2D cell layers to form the 3D organoid structure with minimal impact on tissue growth and differentiation. The intimate contact between the dispersed nanoelectronics and cells enables us to chronically and systematically observe the evolution, propagation, and synchronization of the bursting dynamics in human cardiac organoids through their entire organogenesis.

Scalable spatial single-cell transcriptomics and translatomics in 3D thick tissue blocks.

Characterizing the transcriptional and translational gene expression patterns at the single-cell level within their three-dimensional (3D) tissue context is essential for revealing how genes shape tissue structure and function in health and disease. However, most existing spatial profiling techniques are limited to 5-20 µm thin tissue sections. Here, we developed Deep-STARmap and Deep-RIBOmap, which enable 3D in situ quantification of thousands of gene transcripts and their corresponding translation activities, respectively, within 200-µm thick tissue blocks. This is achieved through scalable probe synthesis, hydrogel embedding with efficient probe anchoring, and robust cDNA crosslinking. We first utilized Deep-STARmap in combination with multicolor fluorescent protein imaging for simultaneous molecular cell typing and 3D neuron morphology tracing in the mouse brain. We also demonstrate that 3D spatial profiling facilitates comprehensive and quantitative analysis of tumor-immune interactions in human skin cancer.

Cyborg islets: implanted flexible electronics reveal principles of human islet electrical maturation.

Flexible electronics implanted during tissue formation enable chronic studies of tissue-wide electrophysiology. Here, we integrate tissue-like stretchable electronics during organogenesis of human stem cell-derived pancreatic islets, stably tracing single-cell extracellular spike bursting dynamics over months of functional maturation. Adapting spike sorting methods from neural studies reveals maturation-dependent electrical patterns of α and ß-like (SC-α and ß) cells, and their stimulus-coupled dynamics. We identified two major electrical states for both SC-α and ß cells, distinguished by their glucose threshold for action potential firing. We find that improved hormone stimulation capacity during extended culture reflects increasing numbers of SC-α/ß cells in low basal firing states, linked to energy and hormone metabolism gene upregulation. Continuous recording during further maturation by entrainment to daily feeding cycles reveals that circadian islet-level hormone secretion rhythms reflect sustained and coordinate oscillation of cell-level SC-α and ß electrical activities. We find that this correlates with cell-cell communication and exocytic network induction, indicating a role for circadian rhythms in coordinating system-level stimulus-coupled responses. Cyborg islets thus reveal principles of electrical maturation that will be useful to build fully functional in vitro islets for research and therapeutic applications.

Brain implantation of tissue-level-soft bioelectronics via embryonic development.

The design of bioelectronics capable of stably tracking brain-wide, single-cell, and millisecond-resolved neural activities in the developing brain is critical to the study of neuroscience and neurodevelopmental disorders. During development, the three-dimensional (3D) structure of the vertebrate brain arises from a 2D neural plate 1,2 . These large morphological changes previously posed a challenge for implantable bioelectronics to track neural activity throughout brain development 3-9 . Here, we present a tissue-level-soft, sub-micrometer-thick, stretchable mesh microelectrode array capable of integrating into the embryonic neural plate of vertebrates by leveraging the 2D-to-3D reconfiguration process of the tissue itself. Driven by the expansion and folding processes of organogenesis, the stretchable mesh electrode array deforms, stretches, and distributes throughout the entire brain, fully integrating into the 3D tissue structure. Immunostaining, gene expression analysis, and behavioral testing show no discernible impact on brain development or function. The embedded electrode array enables long-term, stable, brain-wide, single-unit-single-spike-resolved electrical mapping throughout brain development, illustrating how neural electrical activities and population dynamics emerge and evolve during brain development.

Tracking neural activity from the same cells during the entire adult life of mice.

Stably recording the electrical activity of the same neurons over the adult life of an animal is important to neuroscience research and biomedical applications. Current implantable devices cannot provide stable recording on this timescale. Here, we introduce a method to precisely implant electronics with an open, unfolded mesh structure across multiple brain regions in the mouse. The open mesh structure forms a stable interwoven structure with the neural network, preventing probe drifting and showing no immune response and neuron loss during the year-long implantation. Rigorous statistical analysis, visual stimulus-dependent measurement and unbiased, machine-learning-based analysis demonstrated that single-unit action potentials have been recorded from the same neurons of behaving mice in a very long-term stable manner. Leveraging this stable structure, we demonstrated that the same neurons can be recorded over the entire adult life of the mouse, revealing the aging-associated evolution of single-neuron activities.

Tissue-embedded stretchable nanoelectronics reveal endothelial cell-mediated electrical maturation of human 3D cardiac microtissues.

Clinical translation of stem cell therapies for heart disease requires electrical integration of transplanted cardiomyocytes. Generation of electrically matured human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) is critical for electrical integration. Here, we found that hiPSC-derived endothelial cells (hiPSC-ECs) promoted the expression of selected maturation markers in hiPSC-CMs. Using tissue-embedded stretchable mesh nanoelectronics, we achieved a long-term stable map of human three-dimensional (3D) cardiac microtissue electrical activity. The results revealed that hiPSC-ECs accelerated the electrical maturation of hiPSC-CMs in 3D cardiac microtissues. Machine learning-based pseudotime trajectory inference of cardiomyocyte electrical signals further revealed the electrical phenotypic transition path during development. Guided by the electrical recording data, single-cell RNA sequencing identified that hiPSC-ECs promoted cardiomyocyte subpopulations with a more mature phenotype, and multiple ligand-receptor interactions were up-regulated between hiPSC-ECs and hiPSC-CMs, revealing a coordinated multifactorial mechanism of hiPSC-CM electrical maturation. Collectively, these findings show that hiPSC-ECs drive hiPSC-CM electrical maturation via multiple intercellular pathways.

Explainable multi-task learning for multi-modality biological data analysis.

Current biotechnologies can simultaneously measure multiple high-dimensional modalities (e.g., RNA, DNA accessibility, and protein) from the same cells. A combination of different analytical tasks (e.g., multi-modal integration and cross-modal analysis) is required to comprehensively understand such data, inferring how gene regulation drives biological diversity and functions. However, current analytical methods are designed to perform a single task, only providing a partial picture of the multi-modal data. Here, we present UnitedNet, an explainable multi-task deep neural network capable of integrating different tasks to analyze single-cell multi-modality data. Applied to various multi-modality datasets (e.g., Patch-seq, multiome ATAC + gene expression, and spatial transcriptomics), UnitedNet demonstrates similar or better accuracy in multi-modal integration and cross-modal prediction compared with state-of-the-art methods. Moreover, by dissecting the trained UnitedNet with the explainable machine learning algorithm, we can directly quantify the relationship between gene expression and other modalities with cell-type specificity. UnitedNet is a comprehensive end-to-end framework that could be broadly applicable to single-cell multi-modality biology. This framework has the potential to facilitate the discovery of cell-type-specific regulation kinetics across transcriptomics and other modalities.

Stretchable Mesh Nanoelectronics for 3D Single-Cell Chronic Electrophysiology from Developing Brain Organoids.

Human induced pluripotent stem cell derived brain organoids have shown great potential for studies of human brain development and neurological disorders. However, quantifying the evolution of the electrical properties of brain organoids during development is currently limited by the measurement techniques, which cannot provide long-term stable 3D bioelectrical interfaces with developing brain organoids. Here, a cyborg brain organoid platform is reported, in which "tissue-like" stretchable mesh nanoelectronics are designed to match the mechanical properties of brain organoids and to be folded by the organogenetic process of progenitor or stem cells, distributing stretchable electrode arrays across the 3D organoids. The tissue-wide integrated stretchable electrode arrays show no interruption to brain organoid development, adapt to the volume and morphological changes during brain organoid organogenesis, and provide long-term stable electrical contacts with neurons within brain organoids during development. The seamless and noninvasive coupling of electrodes to neurons enables long-term stable, continuous recording and captures the emergence of single-cell action potentials from early-stage brain organoid development.

ClusterMap for multi-scale clustering analysis of spatial gene expression.

Quantifying RNAs in their spatial context is crucial to understanding gene expression and regulation in complex tissues. In situ transcriptomic methods generate spatially resolved RNA profiles in intact tissues. However, there is a lack of a unified computational framework for integrative analysis of in situ transcriptomic data. Here, we introduce an unsupervised and annotation-free framework, termed ClusterMap, which incorporates the physical location and gene identity of RNAs, formulates the task as a point pattern analysis problem, and identifies biologically meaningful structures by density peak clustering (DPC). Specifically, ClusterMap precisely clusters RNAs into subcellular structures, cell bodies, and tissue regions in both two- and three-dimensional space, and performs consistently on diverse tissue types, including mouse brain, placenta, gut, and human cardiac organoids. We demonstrate ClusterMap to be broadly applicable to various in situ transcriptomic measurements to uncover gene expression patterns, cell niche, and tissue organization principles from images with high-dimensional transcriptomic profiles.