Publications

COVID-19 morbidity and mortality are increased via unknown mechanisms in patients with diabetes and kidney disease. SARS-CoV-2 uses angiotensin-converting enzyme 2 (ACE2) for entry into host cells. Because ACE2 is a susceptibility factor for infection, we investigated how diabetic kidney disease and medications alter ACE2 receptor expression in kidneys. Single cell RNA profiling of kidney biopsies from healthy living donors and patients with diabetic kidney disease revealed ACE2 expression primarily in proximal tubular epithelial cells. This cell-specific localization was confirmed by in situ hybridization. ACE2 expression levels were unaltered by exposures to renin-angiotensin-aldosterone system inhibitors in diabetic kidney disease. Bayesian integrative analysis of a large compendium of public -omics datasets identified molecular network modules induced in ACE2-expressing proximal tubular epithelial cells in diabetic kidney disease (searchable at hb.flatironinstitute.org/covid-kidney) that were linked to viral entry, immune activation, endomembrane reorganization, and RNA processing. The diabetic kidney disease ACE2-positive proximal tubular epithelial cell module overlapped with expression patterns seen in SARS-CoV-2–infected cells. Similar cellular programs were seen in ACE2-positive proximal tubular epithelial cells obtained from urine samples of 13 hospitalized patients with COVID-19, suggesting a consistent ACE2-coregulated proximal tubular epithelial cell expression program that may interact with the SARS-CoV-2 infection processes. Thus SARS-CoV-2 receptor networks can seed further research into risk stratification and therapeutic strategies for COVID-19–related kidney damage.

Although genetic alterations are initial drivers of disease, aberrantly activated transcriptional regulatory programs are often responsible for the maintenance and progression of cancer. CRLF2-overexpression in B-ALL patients leads to activation of JAK-STAT, PI3K and ERK/MAPK signaling pathways and is associated with poor outcome. Although inhibitors of these pathways are available, there remains the issue of treatment-associated toxicities, thus it is important to identify new therapeutic targets. Using a network inference approach, we reconstructed a B-ALL specific transcriptional regulatory network to evaluate the impact of CRLF2-overexpression on downstream regulatory interactions.

Comparing RNA-seq from CRLF2-High and other B-ALL patients (CRLF2-Low), we defined a CRLF2-High gene signature. Patient-specific chromatin accessibility was interrogated to identify altered putative regulatory elements that could be linked to transcriptional changes. To delineate these regulatory interactions, a B-ALL cancer-specific regulatory network was inferred using 868 B-ALL patient samples from the NCI TARGET database coupled with priors generated from ATAC-seq peak TF-motif analysis. CRISPRi, siRNA knockdown and ChIP-seq of nine TFs involved in the inferred network were analyzed to validate predicted TF-gene regulatory interactions.

In this study, a B-ALL specific regulatory network was constructed using ATAC-seq derived priors. Inferred interactions were used to identify differential patient-specific transcription factor activities predicted to control CRLF2-High deregulated genes, thereby enabling identification of new potential therapeutic targets.

The spindle shows remarkable diversity, and changes in an integrated fashion, as cells vary over evolution. Here, we provide a mechanistic explanation for variations in the first mitotic spindle in nematodes. We used a combination of quantitative genetics and biophysics to rule out broad classes of models of the regulation of spindle length and dynamics, and to establish the
importance of a balance of cortical pulling forces acting in different directions. These experiments led us to construct a model of cortical pulling forces in which the stoichiometric interactions of microtubules and force generators (each force generator can bind only one microtubule), is key to\ explaining the dynamics of spindle positioning and elongation, and spindle final length and scaling with cell size. This model accounts for variations in all the spindle traits we studied here, both within species and across nematode species spanning over 100 million years of evolution.

Precise discrimination of tumor from normal tissues remains a major roadblock for therapeutic efficacy of chimeric antigen receptor (CAR) T cells. Here, we perform a comprehensive in silico screen to identify multi-antigen signatures that improve tumor discrimination by CAR T cells engineered to integrate multiple antigen inputs via Boolean logic, e.g., AND and NOT. We screen >2.5 million dual antigens and ∼60 million triple antigens across 33 tumor types and 34 normal tissues. We find that dual antigens significantly outperform the best single clinically investigated CAR targets and confirm key predictions experimentally. Further, we identify antigen triplets that are predicted to show close to ideal tumor-versus-normal tissue discrimination for several tumor types. This work demonstrates the potential of 2- to 3-antigen Boolean logic gates for improving tumor discrimination by CAR T cell therapies. Our predictions are available on an interactive web server resource (antigen.princeton.edu).

The evolution of landscapes, landforms, and other natural structures involves highly interactive physical and chemical processes that often lead to intriguing shapes and recurring motifs. Particularly intricate and fine-scale features characterize the so-called karst morphologies formed by mineral dissolution into water. An archetypal form is the tall, slender, and sharply tipped karst pinnacle or rock spire that appears in multitudes in striking landforms called stone forests, but whose formative mechanisms remain unclear due to complex, fluctuating, and incompletely understood developmental conditions. Here, we demonstrate that exceedingly sharp spires also form under the far-simpler conditions of a solid dissolving into a surrounding liquid. Laboratory experiments on solidified sugars in water show that needlelike pinnacles, as well as bed-of-nails-like arrays of pinnacles, emerge robustly from the dissolution of solids with smooth initial shapes. Although the liquid is initially quiescent and no external flow is imposed, persistent flows are generated along the solid boundary as dense, solute-laden fluid descends under gravity. We use these observations to motivate a mathematical model that links such boundary-layer flows to the shape evolution of the solid. Dissolution induces these natural convective flows that, in turn, enhance dissolution rates, and simulations show that this feedback drives the shape toward a finite-time singularity or blow-up of apex curvature that is cut off once the pinnacle tip reaches microscales. This autogenic mechanism produces ultra-fine structures as an attracting state or natural consequence of the coupled processes at work in the closed solid-fluid system.

Registration of histology images from multiple sources is a pressing problem in large-scale studies of spatial -omics data. Researchers often perform “common coordinate registration,” akin to segmentation, in which samples are partitioned based on tissue type to allow for quantitative comparison of similar regions across samples. Accuracy in such registration requires both high image resolution and global awareness, which mark a difficult balancing act for contemporary deep learning architectures. We present a novel convolutional neural network (CNN) architecture that combines (1) a local classification CNN that extracts features from image patches sampled sparsely across the tissue surface, and (2) a global segmentation CNN that operates on these extracted features. This hybrid network can be trained in an end-to-end manner, and we demonstrate its relative merits over competing approaches on a reference histology dataset as well as two published spatial transcriptomics datasets. We believe that this paradigm will greatly enhance our ability to process spatial -omics data, and has general purpose applications for the processing of high-resolution histology images on commercially available GPUs.

Living matter moves, deforms, and organizes itself. In cells this is made possible by networks of polymer filaments and crosslinking molecules that connect filaments to each other and that act as motors to do mechanical work on the network. For the case of highly cross-linked filament networks, we discuss how the material properties of assemblies emerge from the forces exerted by microscopic agents. First, we introduce a phenomenological model that characterizes the forces that crosslink populations exert between filaments. Second, we derive a theory that predicts the material properties of highly crosslinked filament networks, given the crosslinks present. Third, we discuss which properties of crosslinks set the material properties and behavior of highly crosslinked cytoskeletal networks. The work presented here, will enable the better understanding of cytoskeletal mechanics and its molecular underpinnings. This theory is also a first step towards a theory of how molecular perturbations impact cytoskeletal organization, and provides a framework for designing cytoskeletal networks with desirable properties in the lab.

A major obstacle to treating Alzheimer’s disease (AD) is our lack of understanding of the molecular mechanisms underlying selective neuronal vulnerability, a key characteristic of the disease. Here, we present a framework integrating high-quality neuron-type-specific molecular profiles across the lifetime of the healthy mouse, which we generated using bacTRAP, with postmortem human functional genomics and quantitative genetics data. We demonstrate human-mouse conservation of cellular taxonomy at the molecular level for neurons vulnerable and resistant in AD, identify specific genes and pathways associated with AD neuropathology, and pinpoint a specific functional gene module underlying selective vulnerability, enriched in processes associated with axonal remodeling, and affected by amyloid accumulation and aging. We have made all cell-type-specific profiles and functional networks available at http://alz.princeton.edu. Overall, our study provides a molecular framework for understanding the complex interplay between Aβ, aging, and neurodegeneration within the most vulnerable neurons in AD.

A generic flow distribution network typically does not deliver its load at a uniform rate across a service area, instead oversupplying regions near the nutrient source while leaving downstream regions undersupplied. In this work we demonstrate how a local adaptive rule coupling tissue growth with nutrient density results in a flow network that self-organizes to deliver nutrients uniformly. This geometric adaptive rule can be generalized and imported to mechanics-based adaptive models to address the effects spatial gradients in nutrients or growth factors in tissues.

Loss of normal tissue architecture is a hallmark of oncogenic transformation1. In developing organisms, tissues architectures are sculpted by mechanical forces during morphogenesis2. However, the origins and consequences of tissue architecture during tumorigenesis remain elusive. In skin, premalignant basal cell carcinomas form ‘buds’, while invasive squamous cell carcinomas initiate as ‘folds’. Here, using computational modelling, genetic manipulations and biophysical measurements, we identify the biophysical underpinnings and biological consequences of these tumour architectures. Cell proliferation and actomyosin contractility dominate tissue architectures in monolayer, but not multilayer, epithelia. In stratified epidermis, meanwhile, softening and enhanced remodelling of the basement membrane promote tumour budding, while stiffening of the basement membrane promotes folding. Additional key forces stem from the stratification and differentiation of progenitor cells. Tumour-specific suprabasal stiffness gradients are generated as oncogenic lesions progress towards malignancy, which we computationally predict will alter extensile tensions on the tumour basement membrane. The pathophysiologic ramifications of this prediction are profound. Genetically decreasing the stiffness of basement membranes increases membrane tensions in silico and potentiates the progression of invasive squamous cell carcinomas in vivo. Our findings suggest that mechanical forces—exerted from above and below progenitors of multilayered epithelia—function to shape premalignant tumour architectures and influence tumour progression.