Publications

During epithelial morphogenesis, cell polarity aligns individual cell behaviors into collective motions that shape developing tissues. Here, we combine experiments with computational modeling to investigate how cell-scale forces oriented by Planar Cell Polarity (PCP) direct the collective, counter-rotational cell flows that occur during hair placode morphogenesis. We rule out that PCP directs apical neighbor exchanges, as junctional myosin and PCP protein localization are not co-correlated with junction shrinkage. Instead, we find that PCP directs anterior-directed crawling of placode cells along the basal surface of the tissue through a mechanism that requires cell crawling regulator Rac1. Modeling the placode as a continuum viscoelastic fluid, we find that active forces from cell crawling at the basal surface is sufficient to generate the experimentally observed counter-rotational cell motion at the apical surface. Our results show an unexpected role for PCP in epithelial morphogenesis, centering the basal surface as the site of force generation.

Elliptic partial differential equations (PDEs) can model many physical phenomena, such as electrostatics, acoustics, wave propagation, and diffusion. In scientific machine learning settings, a high-throughput PDE solver may be required to generate a training dataset, run in the inner loop of an iterative algorithm, or interface directly with a deep neural network. To provide value to machine learning users, such a PDE solver must be compatible with standard automatic differentiation frameworks, scale efficiently when run on graphics processing units (GPUs), and maintain high accuracy for a large range of input parameters. We have designed the jaxhps package with these use-cases in mind by implementing a highly efficient and accurate solver for elliptic problems with native hardware acceleration and automatic differentiation support. This is achieved by expressing a highly-efficient solution method for elliptic PDEs in JAX (Bradbury et al., 2018). This software implements algorithms specifically designed for fast GPU execution of a family of elliptic PDE solvers, which are described in full in Melia et al. (2025).

3D chromatin structure is critical for the regulation of gene expression during development. Here we used Micro-C assays at 100-bp resolution to map genome organization in Drosophila melanogaster throughout the first half of embryogenesis. These high-resolution contact maps reveal fine-scale features such as loops and boundaries delineating topologically associating domains. Notably, we observe that 3D chromatin structures form prior to zygotic genome activation and persist during successive mitotic cycles. Integrative analysis with 149 public chromatin immunoprecipitation sequencing (ChIP-seq) datasets identifies four classes of chromatin structuring elements, including a distinct group enriched for GAGA-associated factor (GAF) and Zelda binding, associated with developmental-gene regulation. These elements are mitotically retained and exhibit sequence and structure similarity between D. melanogaster and D. virilis. We propose that 3D chromatin organization in the pre-cellular embryo facilitates deployment of developmentally regulated genes during Drosophila embryogenesis.

Despite the rise of single particle cryo-electron microscopy (cryo-EM) as a premier method for resolving macromolecular structures at atomic resolution, methods to address molecular heterogeneity in vitrified samples have yet to reach maturity. With an increasing number of new methods to analyze the multitude of heterogeneous states captured in single particle images, a systematic approach to validation in this field is needed. With this motivation, we issued a challenge to the community to analyze two cryo-EM particle image sets of thyroglobulin that exhibit continuous conformational heterogeneity. The first dataset was experimental and the second was generated with a simulator, allowing control over the distribution of molecular structures and enabled direct comparison between participants’ submissions and the ground truth molecular structures and distributions. Participants were asked to submit 80 volumes representing the heterogeneous ensemble and estimate their respective populations in the image sets provided. Participation of the research community in the challenge was strong, with submissions from nearly all developers of heterogeneity methods, resulting in 41 submissions across both datasets. Submissions qualitatively exceeded expectations, with the molecular motions identified by methods resembling both each other and the ground truth motion. However, quantitatively assessing these similarities was a challenge in and of itself. In the process of assessing the submissions, we developed several validation metrics, most of which require reference to the underlying ground truth volumes. However, we have also explored the use of metrics that do not necessarily reference ground truth. This is particularly apt for experimental datasets where ground truth is inaccessible. These approaches allowed us to assess the similarity and accuracy in volume quality, molecular motions, and conformational distribution of di!erent submissions. These metrics and the e!orts of all participants help chart a path forward for the improvements of heterogeneity methods for cryo-EM and for future challenges to validate these new methods as they continue to be developed by the community.

Septins can assemble into scaffolds at the plasma membrane to regulate cell morphology. While septins preferentially bind convex membranes via amphipathic helices, their assembly on varied geometries in cells suggests additional localization cues. We tested the hypothesis that lipid composition directs septin assembly through the property of lipid packing. We used pharmacological perturbations that alter fatty acid chain saturation to manipulate lipid packing and found septin structures were selectively disrupted at flat regions of the plasma membrane. To determine whether lipid packing is sufficient to impact septin assembly, molecular dynamics simulations were used to design lipid mixtures with varied packing to monitor septin adsorption in vitro. Septins strongly favored loosely packed lipid bilayers, but additional geometrical cues act in conjunction with this membrane property. This work demonstrates that packing defects and geometry jointly regulate septin localization, highlighting how distinct membrane properties are integrated to organize the septin cytoskeleton.

The formation and maintenance of the finite mammalian ovarian reserve are critical for fertility and species survival. Genetic and developmental studies have uncovered various mechanisms underlying oocyte development and maturation, revealing two curious features of the ovarian germline: (a) The establishment of the follicle reserve involves an initial massive overproduction of oocyte precursors, and (b) the total number of ovulated oocytes across an animal's fertile lifetime is a very small proportion of the initial ovarian reserve. Many have proposed that this indicates the existence of selective quality control to ensure gamete fitness. Here, we review the findings underlying the hypotheses for germline quality control during prepubertal development, homeostatic fertility, and reproductive aging. We evaluate whether the existing evidence base distinguishes the active selection of specific germ cell subsets from neutral dynamics. Throughout, we discuss strategies for applying statistical frameworks to evaluate selection in oogenesis and the implications of neutrality versus selection at various points in oocyte development.

Two-dimensional moiré materials are formed by artificially stacking atomically thin monolayers. Correlated and topological quantum phases can be engineered by precise choice of stacking geometry1--3. These designer electronic properties depend crucially on interlayer coupling and atomic registry4,5. An open question is how the atomic registry responds on ultrafast timescales to optical excitation and whether the moiré geometry can be dynamically reconfigured to tune emergent phenomena in real time. Here we show that femtosecond photoexcitation drives a coherent twist--untwist motion of the moiré superlattice in 2° and 57° twisted WSe2/MoSe2 heterobilayers, resolved directly by ultrafast electron diffraction. On above-band-gap photoexcitation, the moiré superlattice diffraction features are enhanced within 1ps and subsequently suppressed several picoseconds after, deviating markedly from typical photoinduced lattice heating. Kinetic diffraction analysis, supported by simulations of the sample dynamics, indicates a peak-to-trough local twist angle modulation of 0.6°, correlated with a sub-THz frequency moiré phonon. This motion is driven by ultrafast charge transfer that transiently increases interlayer attraction. Our results could lead to ultrafast control of moiré periodic lattice distortions and, by extension, the local moiré potential that shapes excitons, polarons and correlation-driven behaviours.

We model sensory streams as observations from high-dimensional stochastic dynamical systems and conceptualize sensory neurons as self-supervised learners of compact representations of such dynamics. From prior experience, neurons learn coherent sets-regions of stimulus state space whose trajectories evolve cohesively over finite times-and assign membership indices to new stimuli. Coherent sets are identified via spectral clustering of the stochastic Koopman operator (SKO), where the sign pattern of a subdominant singular function partitions the state space into minimally coupled regions. For multivariate Ornstein-Uhlenbeck processes, this singular function reduces to a linear projection onto the dominant singular vector of the whitened state-transition matrix. Encoding this singular vector as a receptive field enables neurons to compute membership indices via the projection sign in a biologically plausible manner. Each neuron detects either a predictive coherent set (stimuli with common futures) or a retrospective coherent set (stimuli with common pasts), suggesting a functional dichotomy among neurons. Since neurons lack access to explicit dynamical equations, the requisite singular vectors must be estimated directly from data, for example, via past-future canonical correlation analysis on lag-vector representations-an approach that naturally extends to nonlinear dynamics. This framework provides a novel account of neuronal temporal filtering, the ubiquity of rectification in neural responses, and known functional dichotomies. Coherent-set clustering thus emerges as a fundamental computation underlying sensory processing and transferable to bio-inspired artificial systems.

The first-order continuum partial differential equation (PDE) model proposed by Bistritzer and MacDonald [Proc. Natl. Acad. Sci. U. S. A. 108, 12233–12237 (2011)] accurately describes the single-particle electronic properties of twisted bilayer graphene at small twist angles. In this paper, we obtain higher-order corrections to the Bistritzer–MacDonald (BM) model via a systematic multiple-scales expansion. We prove that the solution of the resulting higher-order PDE model accurately approximates the corresponding tight-binding wave function under a natural choice of parameters and given initial conditions that are spectrally localized to the monolayer Dirac points. Numerical simulations of tight-binding and continuum dynamics demonstrate the validity of the higher-order continuum model. Symmetries of the higher-order models are also discussed. This work extends the analysis from Watson et al., J. Math. Phys. 64, 031502 (2023), which rigorously established the validity of the (first-order) BM model.

Viral infections can induce prolonged changes in innate immunity. Here, we use blood samples from a human influenza H3N2 challenge study (NCT03883113) to perform comprehensive multi-omics analyses. We detect remodeling of immune programs in circulating innate immune cells that persist after resolution of the infection. We find changes associated with suppressed inflammation, including decreased cytokine and AP-1 gene expression as well as decreased accessibility at AP-1 targets and interleukin-related gene promoter regions. We also find decreased histone deacetylase gene expression, increased MAP kinase gene expression, and increased accessibility at interferon-related gene promoter regions. Genes involved in inflammation and methylation remodeling show modulation of gene-chromatin site regulatory circuit activity. These results reveal a coordinated rewiring of the molecular landscape in innate immune cells induced by mild influenza virus infection.