Engineering programmable material-to-cell pathways via synthetic notch receptors to spatially control differentiation in multicellular constructs.

Synthetic Notch (synNotch) receptors are genetically encoded, modular synthetic receptors that enable mammalian cells to detect environmental signals and respond by activating user-prescribed transcriptional programs. Although some materials have been modified to present synNotch ligands with coarse spatial control, applications in tissue engineering generally require extracellular matrix (ECM)-derived scaffolds and/or finer spatial positioning of multiple ligands. Thus, we develop here a suite of materials that activate synNotch receptors for generalizable engineering of material-to-cell signaling. We genetically and chemically fuse functional synNotch ligands to ECM proteins and ECM-derived materials. We also generate tissues with microscale precision over four distinct reporter phenotypes by culturing cells with two orthogonal synNotch programs on surfaces microcontact-printed with two synNotch ligands. Finally, we showcase applications in tissue engineering by co-transdifferentiating fibroblasts into skeletal muscle or endothelial cell precursors in user-defined micropatterns. These technologies provide avenues for spatially controlling cellular phenotypes in mammalian tissues.

Patterning and folding of intestinal villi by active mesenchymal dewetting.

Tissue folds are structural motifs critical to organ function. In the intestine, bending of a flat epithelium into a periodic pattern of folds gives rise to villi, finger-like protrusions that enable nutrient absorption. However, the molecular and mechanical processes driving villus morphogenesis remain unclear. Here, we identify an active mechanical mechanism that simultaneously patterns and folds the intestinal epithelium to initiate villus formation. At the cellular level, we find that PDGFRA+ subepithelial mesenchymal cells generate myosin II-dependent forces sufficient to produce patterned curvature in neighboring tissue interfaces. This symmetry-breaking process requires altered cell and extracellular matrix interactions that are enabled by matrix metalloproteinase-mediated tissue fluidization. Computational models, together with in vitro and in vivo experiments, revealed that these cellular features manifest at the tissue level as differences in interfacial tensions that promote mesenchymal aggregation and interface bending through a process analogous to the active dewetting of a thin liquid film.

Commensal bacteria promote type I interferon signaling to maintain immune tolerance in mice.

Type I interferons (IFNs) exert a broad range of biological effects important in coordinating immune responses, which have classically been studied in the context of pathogen clearance. Yet, whether immunomodulatory bacteria operate through IFN pathways to support intestinal immune tolerance remains elusive. Here, we reveal that the commensal bacterium, Bacteroides fragilis, utilizes canonical antiviral pathways to modulate intestinal dendritic cells (DCs) and regulatory T cell (Treg) responses. Specifically, IFN signaling is required for commensal-induced tolerance as IFNAR1-deficient DCs display blunted IL-10 and IL-27 production in response to B. fragilis. We further establish that IFN-driven IL-27 in DCs is critical in shaping the ensuing Foxp3+ Treg via IL-27Rα signaling. Consistent with these findings, single-cell RNA sequencing of gut Tregs demonstrated that colonization with B. fragilis promotes a distinct IFN gene signature in Foxp3+ Tregs during intestinal inflammation. Altogether, our findings demonstrate a critical role of commensal-mediated immune tolerance via tonic type I IFN signaling.

Spatial transcriptomics defines injury-specific microenvironments in the adult mouse kidney and novel cellular interactions in regeneration and disease.

Kidney injury disrupts the intricate renal architecture and triggers limited regeneration, and injury-invoked inflammation and fibrosis. Deciphering molecular pathways and cellular interactions driving these processes is challenging due to the complex renal architecture. Here, we applied single cell spatial transcriptomics to examine ischemia-reperfusion injury in the mouse kidney. Spatial transcriptomics revealed injury-specific and spatially-dependent gene expression patterns in distinct cellular microenvironments within the kidney and predicted Clcf1-Crfl1 in a molecular interplay between persistently injured proximal tubule cells and neighboring fibroblasts. Immune cell types play a critical role in organ repair. Spatial analysis revealed cellular microenvironments resembling early tertiary lymphoid structures and identified associated molecular pathways. Collectively, this study supports a focus on molecular interactions in cellular microenvironments to enhance understanding of injury, repair and disease.

TRILL: ORCHESTRATING MODULAR DEEP-LEARNING WORKFLOWS FOR DEMOCRATIZED, SCALABLE PROTEIN ANALYSIS AND ENGINEERING.

Deep-learning models have been rapidly adopted by many fields, partly due to the deluge of data humanity has amassed. In particular, the petabases of biological sequencing data enable the unsupervised training of protein language models that learn the "language of life." However, due to their prohibitive size and complexity, contemporary deep-learning models are often unwieldy, especially for scientists with limited machine learning backgrounds. TRILL (TRaining and Inference using the Language of Life) is a platform for creative protein design and discovery. Leveraging several state-of-the-art models such as ESM-2, DiffDock, and RFDiffusion, TRILL allows researchers to generate novel proteins, predict 3-D structures, extract high-dimensional representations of proteins, functionally classify proteins and more. What sets TRILL apart is its ability to enable complex pipelines by chaining together models and effectively merging the capabilities of different models to achieve a sum greater than its individual parts. Whether using Google Colab with one GPU or a supercomputer with hundreds, TRILL allows scientists to effectively utilize models with millions to billions of parameters by using optimized training strategies such as ZeRO-Offload and distributed data parallel. Therefore, TRILL not only bridges the gap between complex deep-learning models and their practical application in the field of biology, but also simplifies the orchestration of these models into comprehensive workflows, democratizing access to powerful methods. Documentation: https://trill.readthedocs.io/en/latest/home.html.

Recovery of missing single-cell RNA-sequencing data with optimized transcriptomic references.

Single-cell RNA-sequencing (scRNA-seq) is an indispensable tool for characterizing cellular diversity and generating hypotheses throughout biology. Droplet-based scRNA-seq datasets often lack expression data for genes that can be detected with other methods. Here we show that the observed sensitivity deficits stem from three sources: (1) poor annotation of 3' gene ends; (2) issues with intronic read incorporation; and (3) gene overlap-derived read loss. We show that missing gene expression data can be recovered by optimizing the reference transcriptome for scRNA-seq through recovering false intergenic reads, implementing a hybrid pre-mRNA mapping strategy and resolving gene overlaps. We demonstrate, with a diverse collection of mouse and human tissue data, that reference optimization can substantially improve cellular profiling resolution and reveal missing cell types and marker genes. Our findings argue that transcriptomic references need to be optimized for scRNA-seq analysis and warrant a reanalysis of previously published datasets and cell atlases.

Patterning and folding of intestinal villi by active mesenchymal dewetting.

Tissue folding generates structural motifs critical to organ function. In the intestine, bending of a flat epithelium into a periodic pattern of folds gives rise to villi, the numerous finger-like protrusions that are essential for nutrient absorption. However, the molecular and mechanical mechanisms driving the initiation and morphogenesis of villi remain a matter of debate. Here, we identify an active mechanical mechanism that simultaneously patterns and folds intestinal villi. We find that PDGFRA+ subepithelial mesenchymal cells generate myosin II-dependent forces sufficient to produce patterned curvature in neighboring tissue interfaces. At the cell-level, this occurs through a process dependent upon matrix metalloproteinase-mediated tissue fluidization and altered cell-ECM adhesion. By combining computational models with in vivo experiments, we reveal these cellular features manifest at the tissue-level as differences in interfacial tensions that promote mesenchymal aggregation and interface bending through a process analogous to the active de-wetting of a thin liquid film.

Configurational entropy is an intrinsic driver of tissue structural heterogeneity.

Tissues comprise ordered arrangements of cells that can be surprisingly disordered in their details. How the properties of single cells and their microenvironment contribute to the balance between order and disorder at the tissue-scale remains poorly understood. Here, we address this question using the self-organization of human mammary organoids as a model. We find that organoids behave like a dynamic structural ensemble at the steady state. We apply a maximum entropy formalism to derive the ensemble distribution from three measurable parameters - the degeneracy of structural states, interfacial energy, and tissue activity (the energy associated with positional fluctuations). We link these parameters with the molecular and microenvironmental factors that control them to precisely engineer the ensemble across multiple conditions. Our analysis reveals that the entropy associated with structural degeneracy sets a theoretical limit to tissue order and provides new insight for tissue engineering, development, and our understanding of disease progression.

Daily SARS-CoV-2 Nasal Antigen Tests Miss Infected and Presumably Infectious People Due to Viral Load Differences among Specimen Types.

In a recent household transmission study of SARS-CoV-2, we found extreme differences in SARS-CoV-2 viral loads among paired saliva, anterior nares swab (ANS), and oropharyngeal swab specimens collected from the same time point. We hypothesized these differences may hinder low-analytical-sensitivity assays (including antigen rapid diagnostic tests [Ag-RDTs]) by using a single specimen type (e.g., ANS) from reliably detecting infected and infectious individuals. We evaluated daily at-home ANS Ag-RDTs (Quidel QuickVue) in a cross-sectional analysis of 228 individuals and a longitudinal analysis (throughout infection) of 17 individuals enrolled early in the course of infection. Ag-RDT results were compared to reverse transcription-quantitative PCR (RT-qPCR) results and high, presumably infectious viral loads (in each, or any, specimen type). The ANS Ag-RDT correctly detected only 44% of time points from infected individuals on cross-sectional analysis, and this population had an inferred limit of detection of 7.6 × 106 copies/mL. From the longitudinal cohort, daily Ag-RDT clinical sensitivity was very low (<3%) during the early, preinfectious period of the infection. Further, the Ag-RDT detected ≤63% of presumably infectious time points. The poor observed clinical sensitivity of the Ag-RDT was similar to what was predicted based on quantitative ANS viral loads and the inferred limit of detection of the ANS Ag-RDT being evaluated, indicating high-quality self-sampling. Nasal Ag-RDTs, even when used daily, can miss individuals infected with the Omicron variant and even those presumably infectious. Evaluations of Ag-RDTs for detection of infected or infectious individuals should be compared with a composite (multispecimen) infection status to correctly assess performance. IMPORTANCE We reveal three findings from a longitudinal study of daily nasal antigen rapid diagnostic test (Ag-RDT) evaluated against SARS-CoV-2 viral load quantification in three specimen types (saliva, nasal swab, and throat swab) in participants enrolled at the incidence of infection. First, the evaluated Ag-RDT showed low (44%) clinical sensitivity for detecting infected persons at all infection stages. Second, the Ag-RDT poorly detected (≤63%) time points that participants had high and presumably infectious viral loads in at least one specimen type. This poor clinical sensitivity to detect infectious individuals is inconsistent with the commonly held view that daily Ag-RDTs have near-perfect detection of infectious individuals. Third, use of a combination nasal-throat specimen type was inferred by viral loads to significantly improve Ag-RDT performance to detect infectious individuals.

Spatio-temporal patterning of extensile active stresses in microtubule-based active fluids.

Microtubule-based active fluids exhibit turbulent-like autonomous flows, which are driven by the molecular motor powered motion of filamentous constituents. Controlling active stresses in space and time is an essential prerequisite for controlling the intrinsically chaotic dynamics of extensile active fluids. We design single-headed kinesin molecular motors that exhibit optically enhanced clustering and thus enable precise and repeatable spatial and temporal control of extensile active stresses. Such motors enable rapid, reversible switching between flowing and quiescent states. In turn, spatio-temporal patterning of the active stress controls the evolution of the ubiquitous bend instability of extensile active fluids and determines its critical length dependence. Combining optically controlled clusters with conventional kinesin motors enables one-time switching from contractile to extensile active stresses. These results open a path towards real-time control of the autonomous flows generated by active fluids.

D-SPIN constructs gene regulatory network models from multiplexed scRNA-seq data revealing organizing principles of cellular perturbation response.

Gene regulatory networks within cells modulate the expression of the genome in response to signals and changing environmental conditions. Reconstructions of gene regulatory networks can reveal the information processing and control principles used by cells to maintain homeostasis and execute cell-state transitions. Here, we introduce a computational framework, D-SPIN, that generates quantitative models of gene-regulatory networks from single-cell mRNA-seq data sets collected across thousands of distinct perturbation conditions. D-SPIN models the cell as a collection of interacting gene-expression programs, and constructs a probabilistic model to infer regulatory interactions between gene-expression programs and external perturbations. Using large Perturb-seq and drug-response datasets, we demonstrate that D-SPIN models reveal the organization of cellular pathways, sub-functions of macromolecular complexes, and the logic of cellular regulation of transcription, translation, metabolism, and protein degradation in response to gene knockdown perturbations. D-SPIN can also be applied to dissect drug response mechanisms in heterogeneous cell populations, elucidating how combinations of immunomodulatory drugs can induce novel cell states through additive recruitment of gene expression programs. D-SPIN provides a computational framework for constructing interpretable models of gene-regulatory networks to reveal principles of cellular information processing and physiological control.

Extreme differences in SARS-CoV-2 viral loads among respiratory specimen types during presumed pre-infectious and infectious periods.

SARS-CoV-2 viral-load measurements from a single-specimen type are used to establish diagnostic strategies, interpret clinical-trial results for vaccines and therapeutics, model viral transmission, and understand virus-host interactions. However, measurements from a single-specimen type are implicitly assumed to be representative of other specimen types. We quantified viral-load timecourses from individuals who began daily self-sampling of saliva, anterior-nares (nasal), and oropharyngeal (throat) swabs before or at the incidence of infection with the Omicron variant. Viral loads in different specimen types from the same person at the same timepoint exhibited extreme differences, up to 109 copies/mL. These differences were not due to variation in sample self-collection, which was consistent. For most individuals, longitudinal viral-load timecourses in different specimen types did not correlate. Throat-swab and saliva viral loads began to rise as many as 7 days earlier than nasal-swab viral loads in most individuals, leading to very low clinical sensitivity of nasal swabs during the first days of infection. Individuals frequently exhibited presumably infectious viral loads in one specimen type while viral loads were low or undetectable in other specimen types. Therefore, defining an individual as infectious based on assessment of a single-specimen type underestimates the infectious period, and overestimates the ability of that specimen type to detect infectious individuals. For diagnostic COVID-19 testing, these three single-specimen types have low clinical sensitivity, whereas a combined throat-nasal swab, and assays with high analytical sensitivity, was inferred to have significantly better clinical sensitivity to detect presumed pre-infectious and infectious individuals.

Spin glasses, error correcting codes, and synchronization of human stem cell organoids.

Magnetic spins, pendulum clocks, and fireflies all self-organize into coherent collectives when arranged into groups of spatially coupled and interacting individuals. Ramanathan and colleagues demonstrate that spatial coupling of human stem cell organoids induces coherent progression through developmental transitions, allowing the dissection of molecular circuits underlying human development.

Motor processivity and speed determine structure and dynamics of microtubule-motor assemblies.

Active matter systems can generate highly ordered structures, avoiding equilibrium through the consumption of energy by individual constituents. How the microscopic parameters that characterize the active agents are translated to the observed mesoscopic properties of the assembly has remained an open question. These active systems are prevalent in living matter; for example, in cells, the cytoskeleton is organized into structures such as the mitotic spindle through the coordinated activity of many motor proteins walking along microtubules. Here, we investigate how the microscopic motor-microtubule interactions affect the coherent structures formed in a reconstituted motor-microtubule system. This question is of deeper evolutionary significance as we suspect motor and microtubule type contribute to the shape and size of resulting structures. We explore key parameters experimentally and theoretically, using a variety of motors with different speeds, processivities, and directionalities. We demonstrate that aster size depends on the motor used to create the aster, and develop a model for the distribution of motors and microtubules in steady-state asters that depends on parameters related to motor speed and processivity. Further, we show that network contraction rates scale linearly with the single-motor speed in quasi-one-dimensional contraction experiments. In all, this theoretical and experimental work helps elucidate how microscopic motor properties are translated to the much larger scale of collective motor-microtubule assemblies.

Unexplored regions of the protein sequence-structure map revealed at scale by a library of foldtuned language models.

Nature has likely sampled only a fraction of all protein sequences and structures allowed by the laws of biophysics. However, the combinatorial scale of amino-acid sequence-space has traditionally precluded substantive study of the full protein sequence-structure map. In particular, it remains unknown how much of the vast uncharted landscape of far-from-natural sequences consists of alternate ways to encode the familiar ensemble of natural folds; proteins in this category also represent an opportunity to diversify candidates for downstream applications. Here, we characterize sequence-structure mapping in far-from-natural regions of sequence-space guided by the capacity of protein language models (pLMs) to explore sequences outside their natural training data through generation. We demonstrate that pretrained generative pLMs sample a limited structural snapshot of the natural protein universe, including >350 common (sub)domain elements. Incorporating pLM, structure prediction, and structure-based search techniques, we surpass this limitation by developing a novel "foldtuning" strategy that pushes a pretrained pLM into a generative regime that maintains structural similarity to a target protein fold (e.g. TIM barrel, thioredoxin, etc) while maximizing dissimilarity to natural amino-acid sequences. We apply "foldtuning" to build a library of pLMs for >700 naturally-abundant folds in the SCOP database, accessing swaths of proteins that take familiar structures yet lie far from known sequences, spanning targets that include enzymes, immune ligands, and signaling proteins. By revealing protein sequence-structure information at scale outside of the context of evolution, we anticipate that this work will enable future systematic searches for wholly novel folds and facilitate more immediate protein design goals in catalysis and medicine.

Stem cell-derived synthetic embryos self-assemble by exploiting cadherin codes and cortical tension.

Mammalian embryos sequentially differentiate into trophectoderm and an inner cell mass, the latter of which differentiates into primitive endoderm and epiblast. Trophoblast stem (TS), extraembryonic endoderm (XEN) and embryonic stem (ES) cells derived from these three lineages can self-assemble into synthetic embryos, but the mechanisms remain unknown. Here, we show that a stem cell-specific cadherin code drives synthetic embryogenesis. The XEN cell cadherin code enables XEN cell sorting into a layer below ES cells, recapitulating the sorting of epiblast and primitive endoderm before implantation. The TS cell cadherin code enables TS cell sorting above ES cells, resembling extraembryonic ectoderm clustering above epiblast following implantation. Whereas differential cadherin expression drives initial cell sorting, cortical tension consolidates tissue organization. By optimizing cadherin code expression in different stem cell lines, we tripled the frequency of correctly formed synthetic embryos. Thus, by exploiting cadherin codes from different stages of development, lineage-specific stem cells bypass the preimplantation structure to directly assemble a postimplantation embryo.

Localization of signaling receptors maximizes cellular information acquisition in spatially structured natural environments.

Cells in natural environments, such as tissue or soil, sense and respond to extracellular ligands with intricately structured and non-monotonic spatial distributions, sculpted by processes such as fluid flow and substrate adhesion. In this work, we show that spatial sensing and navigation can be optimized by adapting the spatial organization of signaling pathways to the spatial structure of the environment. We develop an information-theoretic framework for computing the optimal spatial organization of a sensing system for a given signaling environment. We find that receptor localization previously observed in cells maximizes information acquisition in simulated natural contexts, including tissue and soil. Specifically, information acquisition is maximized when receptors form localized patches at regions of maximal ligand concentration. Receptor localization extends naturally to produce a dynamic protocol for continuously redistributing signaling receptors, which when implemented using simple feedback, boosts cell navigation efficiency by 30-fold.

Reinforcement learning reveals fundamental limits on the mixing of active particles.

The control of far-from-equilibrium physical systems, including active materials, requires advanced control strategies due to the non-linear dynamics and long-range interactions between particles, preventing explicit solutions to optimal control problems. In such situations, Reinforcement Learning (RL) has emerged as an approach to derive suitable control strategies. However, for active matter systems, it is an important open question how the mathematical structure and the physical properties determine the tractability of RL. In this paper, we demonstrate that RL can only find good mixing strategies for active matter systems that combine attractive and repulsive interactions. Using analytic results from dynamical systems theory, we show that combining both interaction types is indeed necessary for the existence of mixing-inducing hyperbolic dynamics and therefore the ability of RL to find homogeneous mixing strategies. In particular, we show that for drag-dominated translational-invariant particle systems, mixing relies on combined attractive and repulsive interactions. Therefore, our work demonstrates which experimental developments need to be made to make protein-based active matter applicable, and it provides some classification of microscopic interactions based on macroscopic behavior.

Minimal gene set discovery in single-cell mRNA-seq datasets with ActiveSVM.

Sequencing costs currently prohibit the application of single-cell mRNA-seq to many biological and clinical analyses. Targeted single-cell mRNA-sequencing reduces sequencing costs by profiling reduced gene sets that capture biological information with a minimal number of genes. Here we introduce an active learning method that identifies minimal but highly informative gene sets that enable the identification of cell types, physiological states and genetic perturbations in single-cell data using a small number of genes. Our active feature selection procedure generates minimal gene sets from single-cell data by employing an active support vector machine (ActiveSVM) classifier. We demonstrate that ActiveSVM feature selection identifies gene sets that enable ~90% cell-type classification accuracy across, for example, cell atlas and disease-characterization datasets. The discovery of small but highly informative gene sets should enable reductions in the number of measurements necessary for application of single-cell mRNA-seq to clinical tests, therapeutic discovery and genetic screens.