Publications

Comprehensive and spatially mapped molecular atlases of organs at a cellular level are a critical resource to gain insights into pathogenic mechanisms and personalized therapies for diseases. The Kidney Precision Medicine Project (KPMP) is an endeavor to generate 3-dimensional (3D) molecular atlases of healthy and diseased kidney biopsies using multiple state-of-the-art OMICS and imaging technologies across several institutions. Obtaining rigorous and reproducible results from disparate methods and at different sites to interrogate biomolecules at a single cell level or in 3D space is a significant challenge that can be a futile exercise if not well controlled. We describe a "follow the tissue" pipeline for generating a reliable and authentic single cell/region 3D molecular atlas of human adult kidney. Our approach emphasizes quality assurance, quality control, validation and harmonization across different OMICS and imaging technologies from sample procurement, processing, storage, shipping to data generation, analysis and sharing. We established benchmarks for quality control, rigor, reproducibility and feasibility across multiple technologies through a pilot experiment using common source tissue that was processed and analyzed at different institutions and different technologies. A peer review system was established to critically review quality control measures and the reproducibility of data generated by each technology before being approved to interrogate clinical biopsy specimens. The process established economizes the use of valuable biopsy tissue for multi-OMICS and imaging analysis with stringent quality control to ensure rigor and reproducibility of results and serves as a model for precision medicine projects across laboratories, institutions and consortia.

The early Drosophila embryo provides unique experimental advantages for addressing fundamental questions of gene regulation at multiple levels of organization, from individual gene loci to the whole genome. Using Drosophila embryos undergoing the first wave of genome activation, we detected discrete “speckles” of RNA Polymerase II (Pol II), and showed that they overlap with transcribing loci. We characterized the spatial distribution of Pol II speckles and quantified how this distribution changes in the absence of the primary driver of Drosophila genome activation, the pioneer factor Zelda. Although the number and size of Pol II speckles were reduced, indicating that Zelda promotes Pol II speckle formation, we observed a uniform distribution of distances between active genes in the nuclei of both wildtype and zelda mutant embryos. This suggests that the topologically associated domains identified by Hi-C studies do little to spatially constrain groups of transcribed genes at this time. We provide evidence that linear genomic distance between transcribed genes is the primary determinant of measured physical distance between the active loci. Furthermore, we show active genes can have distinct Pol II pools even if the active loci are in close proximity. In contrast to the emerging model whereby active genes are clustered to facilitate co-regulation and sharing of transcriptional resources, our data support an “individualist” model of gene control at early genome activation in Drosophila. This model is in contrast to a “collectivist” model where active genes are spatially clustered and share transcriptional resources, motivating rigorous tests of both models in other experimental systems.

The three-dimensional structure of human chromosomes is tied to gene regulation and replication timing, but there is still a lack of consensus on the computational and biological definitions for chromosomal substructures such as topologically associating domains (TADs). TADs are described and identified by various computational properties leading to different TAD sets with varying compatibility with biological properties such as boundary occupancy of structural proteins. We unify many of these computational and biological targets into one algorithmic framework that jointly maximizes several computational TAD definitions and optimizes TAD selection for a quantifiable biological property. Using this framework, we explore the variability of TAD sets optimized for six different desirable properties of TAD sets: high occupancy of CTCF, RAD21, and H3K36me3 at boundaries, reproducibility between replicates, high intra- vs inter-TAD difference in contact frequencies, and many CTCF binding sites at boundaries. The compatibility of these biological targets varies by cell type, and our results suggest that these properties are better reflected as subpopulations or families of TADs rather than a singular TAD set fitting all TAD definitions and properties. We explore the properties that produce similar TAD sets (reproducibility and inter- vs intra-TAD difference, for example) and those that lead to very different TADs (such as CTCF binding sites and inter- vs intra-TAD contact frequency difference).

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 toward a theory of how molecular perturbations impact cytoskeletal organization, and provides a framework for designing cytoskeletal networks with desirable properties in the lab.

Commonly prescribed antidepressants, such as selective serotonin reuptake inhibitors (SSRIs) take weeks to achieve therapeutic benefits1, 2. The underlying mechanisms of why antidepressants take weeks or months to reverse depressed mood are not understood. Using a single cell sequencing approach, we analyzed gene expression changes in mice subjected to stress-induced depression and determined their temporal response to antidepressant treatment in the cerebral cortex. We discovered that both glial and neuronal cell populations elicit gene expression changes in response to stress, and that these changes are reversed upon treatment with fluoxetine (Prozac), a widely prescribed selective serotonin reuptake inhibitor (SSRI). Upon reproducing the molecular signaling events regulated by fluoxetine3 in a cortical culture system, we found that these transcriptional changes are serotonin-dependent, require reciprocal neuron-glia communication, and involve temporally-specified sequences of autoregulation and cross-regulation between FGF2 and BDNF signaling pathways. Briefly, stimulation of Fgf2 synthesis and signaling directly regulates Bdnf synthesis and secretion cell-non-autonomously requiring neuron-glia interactions, which then activates neuronal BDNF-TrkB signaling to drive longer-term neuronal adaptations4–6 leading to improved mood. Our studies highlight temporal and cell type specific mechanisms promoting the onset of the antidepressant response, that we propose could offer novel avenues for mitigating delayed onset of antidepressant therapies.

We develop a continuum mechanics model of blastocyst hatching. The blastocyst and the zona pellucida are modeled as concentric thick-walled initially spherical shells embedded in a viscous medium. Each shell is characterized by a nonlinear elastic–viscous–constitutive relation. The stiffer outer shell (the zona pellucida) contains an opening. The softer inner shell (the blastocyst) is subject to a continually increasing pressure, which can eventually drive the escape of the inner shell from the outer shell (“hatching”). The focus is on the continuum mechanics modeling framework and illustrating the sort of quantitative predictions that can be made. Numerical examples are presented for the predicted dependence of the evolution of the escape process on values of parameters characterizing the constitutive response of the shells, on the viscosity of the external medium and on the size of the opening in the zona pellucida.

Understanding the changes in diverse molecular pathways underlying the development of breast tumors is critical for improving diagnosis, treatment, and drug development. Here, we used RNA-profiling of canine mammary tumors (CMTs) coupled with a robust analysis framework to model molecular changes in human breast cancer. Our study leveraged a key advantage of the canine model, the frequent presence of multiple naturally occurring tumors at diagnosis, thus providing samples spanning normal tissue, benign and malignant tumors from each patient. We demonstrated human breast cancer signals, at both expression and mutation level, are evident in CMTs. Profiling multiple tumors per patient enabled by the CMT model allowed us to resolve statistically robust transcription patterns and biological pathways specific to malignant tumors versus those arising in benign tumors or shared with normal tissues. We demonstrated that multiple-histological-samples per patient is necessary to effectively capture these progression-related signatures, and that carcinoma-specific signatures are predictive of survival for human breast cancer patients. To catalyze and support similar analyses and use of the CMT model by other biomedical researchers, we provide FREYA, a robust data processing pipeline and statistical analyses framework.

Understanding the changes in diverse molecular pathways underlying the development of breast tumors is critical for improving diagnosis, treatment, and drug development. Here, we used RNA-profiling of canine mammary tumors (CMTs) coupled with a robust analysis framework to model molecular changes in human breast cancer. Our study leveraged a key advantage of the canine model, the frequent presence of multiple naturally occurring tumors at diagnosis, thus providing samples spanning normal tissue and benign and malignant tumors from each patient. We showed human breast cancer signals, at both expression and mutation level, are evident in CMTs. Profiling multiple tumors per patient enabled by the CMT model allowed us to resolve statistically robust transcription patterns and biological pathways specific to malignant tumors versus those arising in benign tumors or shared with normal tissues. We showed that multiple histological samples per patient is necessary to effectively capture these progression-related signatures, and that carcinoma-specific signatures are predictive of survival for human breast cancer patients. To catalyze and support similar analyses and use of the CMT model by other biomedical researchers, we provide FREYA, a robust data processing pipeline and statistical analyses framework.

A classical problem in fluid mechanics is the motion of an axisymmetric vor-tex sheet evolving under the action of surface tension, surrounded by an invis-cid fluid. Lagrangian descriptions of these dynamics are well-known, involv-ing complex nonlocal expressions for the radial and longitudinal velocities interms of elliptic integrals. Here we use these prior results to arrive at a remark-ably compact and exact Eulerian evolution equation for the sheet radius r.´; t/in an explicit flux form associated with the conservation of enclosed volume.The flux appears as an integral involving the pairwise mutual induction formulafor vortex loop pairs first derived by Helmholtz and Maxwell. We show howthe well-known linear stability results for cylindrical vortex sheets in the pres-ence of surface tension and streaming flows [A. M. Sterling and C. A. Sleicher,J. Fluid Mech. 68, 477 (1975)] can be obtained directly from this formulation.Furthermore, the inviscid limit of the empirical model of Eggers and Dupont[J. Fluid Mech. 262 205 (1994); SIAM J. Appl. Math. 60, 1997 (2000)], whichhas served as the basis for understanding singularity formation in droplet pin-choff, is derived within the present formalism as the leading-order term in anasymptotic analysis for long slender axisymmetric vortex sheets and should pro-vide the starting point for a rigorous analysis of singularity formation.

With the rise of commercial spaceflight and prospective human missions to Mars, a wider health range of humans will enter space for longer spans and at higher exposure to environmental stressors than ever before. Numerous adverse health effects have been observed in space, including bone demineralization and skeletal muscle atrophy, among others. Scientists across the world are conducting space omics studies to develop countermeasures for safe and effective crewed space missions. However, optimal extraction of scientific insight from such data is contingent on improved standardization. In response, we founded ISSOP (International Standards for Space Omics Processing), an international consortium of scientists who aim to enhance guidelines between space biologists globally. This paper informs scientists and data scientists from many fields about the challenges and future avenues of space omics and can serve as an introductory reference for new members in the space biology discipline.