
Conference Dates Location Organizing Committee
August 11-14, 2025
Share:
In response to the recent financial uncertainty associated with federal awards,
we have significantly reduced the conference registration fee.
Overview
The NITMB MathBio Convergence Conference will be an international forum for the presentation and discussion of research at the frontier between biology and mathematics. The Conference is sponsored and hosted by a new national institute charged with enhancing greater integration between these disciplines, the NSF-Simons National Institute for Theory and Mathematics in Biology (NITMB). Located in downtown Chicago, the NITMB was founded in 2023 by the National Science Foundation and the Simons Foundation and is a working partnership between Northwestern University and the University of Chicago.
The conference will be a forum to discuss challenging problems in biology that are inspiring new mathematics, statistics, theory and computation and to explore mathematics that is generating important new insights into biology. A half-dozen plenary talks from distinguished international speakers are scheduled at which theory and experiment will be discussed.
There will be 20 speaking slots and two poster sessions selected from submitted abstract applications. The schedule provides ample time for discussion, professional networking and career development activities.
Those who work at or who are interested in the math-bio interface are strongly
encouraged to register for a modest fee to hold their place.
Location
The conference will be held near to the headquarters of the National Institute for Theory and Mathematics in Biology at the Drake Hotel, 140 E. Walton Place, Chicago, Illinois.
Hotel Room Block
A block of hotel rooms has been reserved at The Drake at $199.99 a night. A link to book your room will be sent after registration starting May 22, 2025. Rooms are automatically reserved for invited and selected speakers.
Attendees applying for and receiving travel awards will have their lodging at the conference hotel, meals and travel covered by the NITMB within reimbursement policy. See Travel Awards Tab for more details.
Conference Amenities
Meals
Breakfast and Lunch will be provided by the conference.
Workspace
All registrants will be welcome to work out of NITMB the days of the workshop.
Questions?
Email nitmb@nitmb.org
Confirmed Plenary Speakers
Alexander Aulehla
Head of Developmental Biology Unit
The European Molecular Biology Laboratory (EMBL)
Dr. Aulehla uses dynamical systems to study how a vertebrate embryo forms the unique number of body segments that make up the vertebrae.
Carina Curto
Pablo J. Salame Goldman Sachs Professor of Computational Neuroscience
Brown University
Dr. Curto combines mathematical tools from algebra, geometry, topology and dynamical systems to elucidate complex patterns of activity in networks of neurons.
Lorin Crawford
Principal Researcher & Distinguished Senior Fellow in Biostatistics
Microsoft Research & Brown University
Dr. Crawford uses machine learning and statistical tools to study the genotype-to-phenotype map for complex traits and diseases.
Michael Hawrylycz
Senior Investigator
The Allen Institute
Dr. Hawrylycz studies the cellular constitution and organization of the brain using single-cell transcriptomics and cutting-edge computational and statistical tools.
Anna-Liisa Laine
Professor, Organismal and Evolutionary Biology Research Programme
University of Helsinki
Dr. Laine studies complex problems in ecology across a spectrum of length and time scales.
Reidun Twarock
Professor, Department of Mathematics
University of York
Dr Twarock uses mathematics to study the symmetries of animal viruses and how they self-assemble inside host cells.
Deadlines
Abstract Submission
Travel Fund Application
Early Registration (to be considered for a presentation) ($40)
(Early Registration for Undergraduates: $20)
Notification of Selection for a Talk or Poster
Regular Registration ($40)
(Regular Registration for Undergraduates: $20)
Late Registration ($80)
Hotel Registration*
April 15, 2025 - now closed
April 15, 2025 - now closed
February 15 - April 15, 2025
Late May
February 15 - July 14, 2025
July 15 - August 4, 2025
May 22 - July 11, 2025
*see FAQ for hotel details
Abstract Submission
NITMB MathBio Convergence Conference Abstract Submission Instructions
Abstract Submissions were due by April 15, 2025.
After confirmation of payment, you will be sent a link to submit an abstract via the University of Chicago InfoReady site. You will need to make an account if you don’t already have one. Submissions will be reviewed and presenters notified in May.
The conference will have a single session. There will be approximately 20 speaking slots selected from submissions, and two poster sessions. Priority for the speaking slots will be given to early career researchers.
Abstract submissions should be made in PDF format, and one-page. (If you feel you need more space, a second page is allowed. But please, no more than 2 pages.) On the submission page indicate your choice: Oral talk preferred; Poster presentation preferred; No preference.
Abstract templates will be provided on this submission page in LaTeX (Euclid.tex and ltexpprt.sty) and Word (Darwin.docx) are provided. The LaTeX template does allow a two-column format if desired. Please see the comments in the LaTeX source.
Travel Awards
Travel Awards Submissions were due by April 15, 2025.
After confirmation of registration payment, you will be sent a link to submit an abstract via the University of Chicago InfoReady site. You will need to make an account if you don’t already have one.
NITMB Student Travel Award
Any full-time student in good standing is eligible to receive an award to assist their participation at the NITMB MathBio Convergence Conference.
NITMB Student Travel Awards fund up to $650 for domestic travel and $800 for intercontinental travel, plus reimbursement for local hotel costs at the conference hotel and meals not provided by the conference (gsa.gov rates apply).
Please note, award is subject to all NITMB reimbursement policies found at www.nitmb.org/reimbursement. Registrants receiving a travel award must use The Drake Hotel (140 East Walton Place Chicago, IL 60611) to be eligible for reimbursement.
Required Documents for Application:
● Applicant Statement (1 page max) - Short summary of why you are requesting conference support, the rationale for your request (financial need, scholarship, research, etc.), and your interest in the conference.
● CV - An updated copy of your curriculum vitae.
● Letter of Recommendation - A letter from your advisor or senior colleague noting your achievements, potential, and any special circumstances that may justify receiving an award.
● Abstract - If you have already submitted an abstract or intend to do so, please indicate this on the application portal. A code to use for abstract submission nominally will be sent within one business day of receipt of your application.
NITMB Early Career Travel Award
To qualify as "early career" an applicant must be a postdoctoral or non-tenured faculty member who obtained their Ph.D. in the ten years preceding the conference. NITMB Early Career Travel Awards fund up to $650 for domestic travel and $800 for intercontinental travel, plus reimbursement for local hotel costs and meals not provided by the conference (gsa.gov rates apply)
and
NITMB Teaching Intensive Institutions (TII) and Primarily Undergraduate Institutions (PUI) Faculty and Lecturer Travel Award
To qualify, an applicant must be faculty, a lecturer, or adjunct in mathematics, biology, engineering, physics, or related discipline at an institution within the U.S. that does not offer Ph.D. degrees (specifically, an institution not classified as “Doctoral” by the Carnegie Classification of Institutions of Higher Education, https://carnegieclassifications.acenet.edu/ ).
These awards fund up to $650 for domestic travel and $800 for intercontinental travel, plus reimbursement for local hotel costs and meals not provided by the conference (gsa.gov rates apply).
Required Documents for Application for both Early Career, TII/PUI Faculty, and Lecturer categories:
● Applicant Statement (1 page max) - Short summary of why you are requesting conference support, the rationale for your request (financial need, scholarship, research, etc.), and your interest in the conference.
● CV - An updated copy of your curriculum vitae.
● Abstract - If you have already submitted an abstract or intend to do so, please indicate this on the application portal. A code to use for abstract submission nominally will be sent within one business day of receipt of your application.
Liam O’Brien
Recorded on 8/11/2025
Title: Structural causes of pattern formation and loss through model-independent bifurcation theory
Abstract: During development, precise cellular patterning is essential for the formation of functional tissues and organs. These patterns arise from conserved signaling networks that regulate communication both within and between cells. In this talk, we will develop and present a model-independent ordinary differential equation (ODE) framework for analyzing pattern formation in a homogeneous cell array. In contrast to traditional approaches that focus on specific equations, our method relies solely on general assumptions about global intercellular communication (between cells) and qualitative properties of local intracellular biochemical signaling (within cells). Prior work has shown that global intercellular communication networks alone determine the possible emergent patterns in a generic system. We build on these results by demonstrating that additional constraints on the local intracellular signaling network lead to a single stable pattern which depends on the qualitative features of the network. Our framework enables the prediction of cell fate patterns with minimal modeling assumptions, and provides a powerful tool for inferring unknown interactions within signaling networks by analyzing tissue-level patterns.
Amlan Nayak
Recorded on 8/11/2025
Title: Unraveling prey evasion mechanisms in zebrafish
Abstract: We developed an interactive platform using computer vision and feedback control to study zebrafish (Danio rerio) escape behavior in response to virtual predators and support the development of predictive models grounded in experimental data. Our results show that zebrafish adapt escape strategies based on neural impairments and social context, highlighting the role of internal models and predictive control. Embedding theory within our closed-loop experiments provides a testable framework for understanding the neural and computational basis of survival behavior.
Anna-Liisa Laine
Recorded on 8/11/2025
Title: Uncovering the determinants of disease risk and within host-pathogen diversity in natural systems
Abstract: Disease risk in natural plant populations emerges from complex interactions between host diversity, community structure, and environmental variation. We integrate field experiments, long-term datasets, and molecular tools to explore how both ecological and evolutionary dynamics shape pathogen prevalence and community composition. We show that genetic and species-level host diversity can modulate pathogen spread — sometimes buffering, other times amplifying infection risk-depending on landscape structure and host traits. Our findings emphasize the need to bridge scales of biological organization-from genes to ecosystemsto predict when and where outbreaks will occur. Understanding these dynamics is increasingly important as ecosystems are reshaped by land-use change and climate variability.
Haina Wang
Recorded on 8/11/2025
Title: Rigidity homeostasis of the actin cortex via tensionsensitive filament and crosslinker dynamics
Abstract: The actin cortex is a biopolymer network that maintains mechanical rigidity despite constant structural changes through assembly and disassembly of filaments and crosslinkers. The biological advantage of this energy-intensive remodeling, as well as the microscopic mechanisms underlying rigidity homeostasis, remains unclear. To address these questions, we have developed two independent elastic network models that achieve rigidity homeostasis while undergoing complete architectural turnover, capturing filament and crosslinkers dynamics, respectively. We find that both models require the following minimal ingredients: (1) preferential disassembly of edges or nodes that are under small tension or force, (2) a small fraction of random disassembly, and (3) energy injection upon assembly. The trajectories toward steady states are analyzed with graph-theoretical metrics to reveal the relationship between rigidity and network topology. Remarkably, we find that the edges in our networks undergo diffusion, whereas the elastic moduli and structural correlations remain statistically invariant, analogous to the representational drift found in learning systems in which the architecture self-adjust via local rules. We hypothesize that the cortex functions as a tunable matter, where remodeling dynamics enables the cell to balance robustness and flexibility in fluctuating mechanical environments, creating survival advantages that justify the energy consumption. Our models provide platforms to design and study bio-inspired tunable metamaterials.
Ryan Robinett
Recorded on 8/11/2025
Title: Novel quantification of morphological variation in D. Melanogaster wings using the Wasserstein-FisherRao metric
Abstract: Living systems exhibit remarkable fidelity in form (morphology) when comparing individual organisms from the same species or closely related species. To understand the constraints underlying morphological variation requires geometric analysis of form. However, historical study of morphological geometry has relied on a sparse number of handpicked anatomical landmarks within the given structures. In this talk, we present methods for quantifying morphological differences between anatomical structures of two individuals, as well as for quantifying modes of variation in morphology for a large number of specimens. These methods leverage the Riemannian manifold structure of the WassersteinFisher-Rao metric introduced by Chizat, Peyr´e, Schmitzer, and Vialard1 . We demonstrate the power of these methods by applying them to wing morphology in Drosophila melanogaster, capturing modes of intra- and inter-strain morphological variation
Palash Sashittal
Recorded on 8/11/2025
Title: Inferring cell lineage trees and differentiation maps from lineage tracing data
Abstract: Organismal development is a complex process that involves repeated cell divisions and gradual differentiation of cells through a hierarchy of progenitor cell types, each with progressively restricted differentiation potential. A central goal in developmental biology is to reconstruct the cell lineage tree – which describes the complete history of cell divisions – and the cell differentiation map – which describes the hierarchy of progenitor cell types and cell type transitions that occur during development. Even for simple organisms such as Caenorhabditis elegans, uncovering the cell lineage tree and the differentiation map of has led to key insights, including the critical role of programmed cell death during development [23]. Creating comparable cell lineage trees and differentiation maps for mammalian development will have far-reaching implications in regenerative medicine and advancement of therapies for developmental disorders. Recent advances in genome editing and single-cell sequencing have enabled high-throughput lineage tracing of cells in complex developmental systems [1, 16, 17, 21]. In these technologies, heritable barcodes are induced in cells, either at specific stages of development [4, 7, 13, 14] or dynamically as cells divide and differentiate [8, 11, 18–20, 22], using genome editing tools such as CRISPR-Cas9. Single-cell RNA sequencing is subsequently used to simultaneously measure barcodes (revealing the lineage of cells) and gene expression (revealing cell types) for thousands of individual cells [2, 5, 6, 9, 10, 12, 15]. While these technologies offer the scalability to investigate development in complex organisms, they do not observe every cell division and the differentiation decisions made by each dividing cell during development. As such, we require computational methods to infer cell lineage trees and differentiation maps from single-cell lineage tracing data. In this talk, we will present two algorithms to infer cell lineage trees and differentiation maps from lineage tracing data. First, we will present Startle [6], an algorithm to infer cell lineage trees from CRISPR-Cas9-based lineage tracing. Startle introduces the star homoplasy evolutionary model that constrains a phylogenetic character to mutate at most once along a lineage, capturing the non-modifiability property of CRISPR-Cas9 mutations. We derive a combinatorial characterization of star homoplasy phylogenies and use this characterization in Startle to compute a maximum parsimony star homoplasy phylogeny from lineage tracing data. The second algorithm, CARTA [3], infers cell differentiation maps from lineage tracing data in which observed and unobserved progenitor cell types are represented by their potency. CARTA introduces a novel mathematical framework that balances the trade-off between the complexity of the differentiation map and the fit of the map with the cell lineage trees derived from lineage tracing data. We also introduce a formal quantitative framework to systematically assess different models of cell differentiation maps with varying numbers of progenitors. Throughout the talk, we will demonstrate that Startle and CARTA outperform existing methods on both simulated and real data. Overall, my talk will underscore the crucial role of specialized models and algorithms to effectively analyze lineage tracing data and uncover the dynamics of developmental biological systems.
Eleanor Degen
Recorded on 8/11/2025
Title: Bicoid-nucleosome competition sets a concentration threshold for transcription constrained by genome replication
Abstract: Concentration gradients of transcription factors establish patterns of gene expression during a time in early embryonic development characterized by frequent mitotic divisions, rounds of genome replication, and chromatin reorganization. In the Drosophila embryo, an exponential gradient of the transcription factor Bicoid (Bcd) activates the earliest expressed patterning genes across the anterior-posterior axis. Bcd successfully navigates the chromatin of the replicating genome to bind its sites and facilitate transcription. However, we do not fully understand how genomic context leads to differential transcriptional outputs across Bcd concentrations. We aim to model how enhancer sequence, chromatin, and DNA replication together determine the Bcd concentrationsensitivity of the transcriptional process. By live-imaging an MS2-MCP transcriptional reporter for the hunchback P2 (hbP2) enhancer, we have found that the length of the delay between mitosis and the initiation of transcription uniquely reflects Bcd concentration-sensitive regulation. We have defined a stochastic model of transcriptional regulation that accurately predicts transcriptional onset times using mathematical descriptions of Bcd-nucleosome competition for occupancy and the probability of DNA replication at hbP2. Our work suggests that Bcd’s ability to outcompete nucleosomes dictates a Bcd concentration threshold for expression, while DNA replication limits the rate of transcriptional activation at high concentrations where Bcd readily outcompetes nucleosomes. Disrupting nucleosome stability by promoting pioneer factor binding to hbP2 both expands the MS2 expression domain and preserves onset time variance at high concentrations, supporting the model’s predictions. We provide a theoretical framework for understanding how the critical parameters of nucleosome positioning, stability, and replication dynamics can govern the sensitivity of gene expression to transcription factor concentration.
Boya Hou
Recorded on 8/11/2025
Title: Multi-network RNA velocity
Abstract: RNA velocity is a powerful model for cellular RNA splicing events which can be used to infer dynamical properties of various regulatory functions. Despite its wide applicability, the model has not been adequately adjusted to account for multiway gene regulatory or intercellular interactions. Here, we propose to revise the RNA velocity approach by incorporating two new network components: an intracellular gene regulatory network (GRN) that regulates gene expression, and an intercellular interaction network that captures interactions between (neighboring) cells. By combining GRN and intercellular communication models, we aim to provide a more holistic understanding of single-cell gene expression dynamics through the lens of systems and control theory. In particular, we investigate network steady states, stability and targeted drug intervention.
Michael Hawrylycz
Recorded on 8/11/2025
Title: Discrete Morse graph construction for highdimensional transcriptomic data
Abstract: Discrete Morse theory, a combinatorial adaptation of classical Morse theory, is instrumental in analyzing topological features of discrete spaces like cell complexes. When combined with persistent homology in topological data analysis, it is particularly effective in simplifying complex data structures and extracting meaningful features from high-dimensional or noisy datasets. Single-cell transcriptomic data such as RNA-seq provides an interesting application of discrete Morse theory, where measurements of thousands of genes reflecting differential expression patterns are used to determine cell type and cell state. We introduce a novel algorithm, scDMGC, for discrete Morse graph reconstruction and demonstrate its efficacy in validating brain cell type taxonomies and gradients. Additionally, discrete Morse graph construction is employed to distinguish cell loss and altered transcription in early to late Alzheimer’s disease transcriptomic data.
Juliana Londono Alvarez
Recorded on 8/12/2025
Title: Attractor-based models for sequences and pattern generation in neural circuits
Abstract: Brain rhythms that generate walking, breathing, or swimming are traditionally modeled using coupled oscillators—intrinsically oscillating (“pacemaker”) neurons. Alternatively, attractor networks—a framework in which cognitive processes are modeled as attractors of a dynamical system, typically implemented by a recurrent network—have been a central tool in theoretical neuroscience for modeling memory and pattern completion [1]. In this work, we show that attractor networks can also support rhythmic pattern generation. We present a single attractor network that can robustly encode five distinct quadruped gaits (bound, pace, trot, walk, and pronk) as coexisting attractors, without requiring parameter changes. Transitions between gaits can be triggered by simple external pulses. This contrasts with most existing locomotion models, which rely on finely tuned coupled oscillators whose parameters must be adjusted to switch between gaits. In addition, we introduce a novel method for encoding ordered sequences of gaits using fusion attractors, allowing the network to flexibly reuse existing patterns in different combinations (as in a sequence of dance moves).
Reidun Twarock
Recorded on 8/12/2025
Title: Viruses under the mathematical microscope: viral geometry as a key to understanding viral infections
Abstract: Most viruses have protein shells, called viral capsids, that surround, and thus protect, their genetic material. Due to the highly symmetric nature of these capsids, mathematical techniques from group, graph and tiling theory can be used to model and classify virus architecture. By combining these geometric, and related topological, descriptors of virus architecture with stochastic simulations, I will demonstrate how viral geometry provides insights into viral life cycles that pave the way to innovation in antiviral therapy and virus nanotechnology.
Ning Wei
Recorded on 8/12/2025
Title: The impact of ephaptic coupling and ionic electrodiffusion on arrhythmogenesis in the heart
Abstract: The synchronized contraction of the heart relies on the ordered propagation of the cardiac action potential, which occurs through gap junctions (GJs)-rich intercalated disc (ID). GJs are low-resistance pathways that directly connect the cytoplasm of adjacent cardiac cells for the exchange of ions and small molecules. However, ongoing theoretical and experimental studies have suggested the potential incompatibilities in applying GJ-mediated microscale coupling to explain macroscopic propagation. Notably, GJ knockout mice are still able to maintain heart structure and function, albeit GJs are undetectable at IDs. Ephaptic coupling (EpC) is an electric field effect developed at the ID, which is a narrow and torturous cleft connecting the individual cardiac cells to form a functional syncytium. Recent research has substantiated EpC as an alternative mechanism for cellular communication when GJs are impaired. Given the current absence of direct experimental evidence for the existence of EpC, modeling studies play a vital role in elucidating its physiological and pathological functions in the heart. We developed a two-dimensional (2D) model of EpC that incorporates multidomain electrodiffusion of multiple ions, providing the first physiologically detailed framework for studying ephaptic conduction in both healthy and ischemic heart. Our results showed that strong EpC, combined with ionic electrodiffusion, significantly decreases the occurrence of conduction failure into ischemic region, indicating the potential benefits of EpC. In addition, strong EpC can markedly influence ionic concentrations within the cleft, notably elevating K+ levels and nearly depleting Ca2+ , while changes in Na+ remain moderate. These findings shed light on the underlying mechanism of EpC. Moreover, our findings indicate that sufficiently strong EpC tends to suppress the initiation of reentry, leading to absent or nonsustained reentrant activity, while could also introduce instability and heterogeneity into the cardiac dynamics. In contrast, relatively weak EpC supports sustained reentry with a stable rotor. Furthermore, we demonstrated that strong EpC can terminate reentrant arrhythmias in ischemic hearts with complex anatomical structures. Our findings revealed two distinct termination mechanisms: (1) sufficiently strong EpC leads to rapid self-attenuation of reentry, and (2) moderate EpC results in termination via selfcollision over a longer period, driven by increased conduction velocity and anisotropy. In summary, our research is the first to demonstrate that EpC exerts anti-arrhythmic effects by both preventing the initiation of reentry and terminating established reentrant circuits, fundamentally shaping our understanding of EpC.
Nonthakorn Olaranont
Recorded on 8/12/2025
Title: A cell-based mechanical model captures stress relaxation and flow in proliferating tissues with subcellular elasticity
Abstract: Many developmental shape changes, such as brain convolution and gut folding, are governed by nonlinear mechanics under tissue growth. Yet living tissues also exhibit fluid-like behavior, dynamically rearranging to relieve stress and form complex shapes. It remains unclear how proliferating tissues simultaneously generate elastic forces, relax them, and exhibit fluidity at the intercellular level. We present a vertex-based computational model that captures these multiscale behaviors by embedding subcellular shear elasticity via triangulated reference states and evolving deformation gradients. Cell pressure is tracked through deviations from reference cell size, with proliferation modeled as both a change in cell number and reference configuration. This framework enables simulation of tissue-scale flows with built-in growth and mechanical heterogeneity. Our model reproduces buckling behavior consistent with continuum nonlinear growth-elasticity theory and further reveals how local intra- and intercellular rearrangements mediate stress relaxation and shape change. We also demonstrate how microscopic cellular architecture impacts macroscopic growth-driven flow and morphogenesis. Mathematically, the model leverages smooth optimization techniques, incorporating gradients and Hessians to improve computational performance and allow local linear stability analysis of solutions. We apply the framework to simulate large-scale tissue flows during Drosophila wing disc and dorsal thorax morphogenesis, as well as in vitro circular wound closure in partial EMT monolayers. This predictive tool bridges cellular mechanics with emergent tissue-scale dynamics, with applications to development, regeneration, and cancer.
Orlando Arguello Miranda
Recorded on 8/12/2025
Title: Biomolecular network discovery through GenAIdriven time series analysis in living cells
Abstract: Biological processes are driven by networks of proteins coordinated in time and space that define cell behavior in health and disease. Although protein networks can be inferred from sequencing-based techniques and fluorescence microscopy, these methods often destroy cells to recover information or use biochemical sensors that alter cell physiology. These limitations prevent the direct visualization of complex protein networks in single living cells. For instance, under optimal conditions, up to six fluorescent markers can be reliably tracked in living cells; however, even small molecular networks, such as signaling cascades, can surpass dozens of proteins. Therefore, new methods are required to visualize large protein networks in single living cells and produce more accurate models of cellular behavior. Most methods to increase the number of proteins measured in cells focus on improving microscopy hardware or optimizing biochemical sensors. Here, we propose an alternative approach based on generative AI (GenAI) to produce bio-realistic visualizations using customized algorithms for image-to-image translation. We created generative models that take simple microscopy images, such as phase contrast micrographs, as input and produce a multidimensional array of synthetic images depicting multiple proteins in living cells over time (Fig. 1). To achieve this, we produced a unique ground truth dataset of multidimensional fluorescent images with dozens of proteins encompassing signaling networks in the model organism Saccharomyces cerevisiae. To test whether the time series derived from synthetic images could be used for biological inference, we compared ordinary differential equation models for networks based on time series measurements from real or synthetic images. Our results showed that time series derived from synthetic images resembled real time series, enabling the discovery of protein correlations, cross-correlations, hypothetical network motifs, and protein-protein interactions in single living cells. We envision that biomolecular network discovery through GenAI-driven time series analysis will accelerate research by expanding the information obtained from single images and allowing faster parameter optimizations for network models of complex biological systems.
Jason Kim
Recorded on 8/12/2025
Title: Geometric model manifold of space, time, and belief in hippocampal cognitive maps
Abstract: Cognitive maps are mental representations of spatial and conceptual relationships in an environment. The hippocampus, a key brain region for these maps, exhibits neural population activity that often evolves along nonlinear, low-dimensional manifolds. While cognitive maps have been characterized at the level of cell tunings to spatial location and decorrelations over time, a precise, unsupervised, and quantitative model of how this map forms and morphs at the population level remains elusive due to the large dimensionality and nonlinearity of neural activity. Current dimensionality reduction techniques have difficulty capturing this manifold in an interpretable way due to two shortcomings: they do not model regions of neural activity where there are no data, and the map from the embedding back to neural activity is too nonlinear. We solve these problems in a novel unsupervised technique using differential geometry to model the low-dimensional space in which the data reside as a gently curved manifold. We apply this method to activity from thousands of neurons across several days in the CA1 region of the hippocampus of mice learning to collect rewards under two task conditions on a linear maze. We show that the 3 relevant variables of position (space), trial (time), and task belief emerge as local directions on the population-level manifold that captures track topology, a consistent direction of representational drift within and between days, and task learning that evolves orthogonally to baseline drift (resulting in tuning decorrelation). Further, the manifold quantifies neuron-level features—such as changes in positional tuning across trials—at the population level, and we discover neurons tuned for generalized states including space, time and belief. We provide a novel technique and a unifying theory for how populationlevel neural manifolds encode generalized state-time information along gently curved manifolds, which we hypothesize is for efficient recall and co-localization in downstream regions.
Lorin Crawford
Recorded on 8/12/2025
Title: Statistical opportunities in defining, modeling, and targeting cell state in cancer
Abstract: Project Ex Vivo is a joint cancer research collaboration between Microsoft and the Broad Institute of MIT and Harvard. Our group views cancers as complex (eco)systems, beyond just mutational variation, that necessitate systems-level understanding and intervention. In this talk, I will discuss a series of multimodal statistical and deep learning approaches to understand accurate representations of tumors by integrating genetic markers, expression state, and microenvironmental interactions. These representations help us precisely define and quantify the trajectory of each tumor in each patient. Our ultimate objective is to more effectively model cancer ex vivo – outside the body – in a patient-specific manner. In doing so, we aim to unlock the ability to better stratify patient populations and identify therapies that target diverse aspects of human cancers.
Selimzhan Chalyshkan
Recorded on 8/13/2025
Title: An ensemble modeling framework to predict synapses from optokinetic stimuli in larval zebrafish
Abstract: Studying the link between neural network structure and function is crucial in understanding the principles underlying brain activity. With the growing use of large-scale activity recordings and advances in anatomical circuitlevel reconstructions, researchers can now posit numerous neural network models that tie neural recordings with synaptic connectivity. However, understanding how well various neural network model classes relate to and predict observed neural activity with structural connections remains a significant challenge in systems and computational neuroscience. Here, rather than fitting a single model to find structural parameters matching a network's observed functional profile, we examined an ensemble modeling approach that finds consistent structure across models (Biswas et al., 2024). We applied ensemble modeling in a threshold-linear network to predict connectivity in larval zebrafish pretectum from the most common neuronal response types to optokinetic stimuli. In line with the previously hypothesized model, the strongest predictions included excitatory connections from the retina to monocular response types and connections from the monocular response types to binocular response types. We also tested the ensemble modeling framework by predicting neuronal activity for a given stimulus condition from the ensemble of models that account for the other stimulus conditions. We found that ensemble modeling produced up to five correct predictions for each response type before making an error. In contrast, the fit model that finds the solution that minimizes the weight norm predicted only one or two correct activity responses. Together, we have shown that ensemble modeling is a robust approach that makes accurate activity predictions and compelling structural predictions from thresholdlinear recurrent neural network models. Forthcoming comparisons with the anatomically reconstructed circuit will further evaluate accuracy of the theoretical framework predictions.
Harrison Oatman
Recorded on 8/13/2025
Title: Modeling mitotic wave origins in Drosophila Melanogaster
Abstract: In the first few hours of embryonic development in Drosophila melanogaster, 13 rounds of rapid nuclear divisions bring the count of nuclei from one to nearly six thousand. Such quick advancement offers a selective advantage, but synchronizing cell cycle timing across the ~500 µm embryo with a precision better than five minutes presents a significant organizational challenge. The embryo's solution—studied for over four decades—is a series of mitotic waves originating from the poles of the egg and cascading towards the equator. Biologists and physicists, inspired by this feat of large-scale organization, have modeled these waves as a reaction-diffusion system based on careful observation of cell cycle kinase activity. Still, the precise reasoning for mitotic wave slowdown and the factors influencing mitotic wave origin are unclear. Using high-resolution in toto lightsheet live imaging combined with newly developed machine learning methods for nuclear segmentation and lineage tracking, I provide a single-cell perspective on this phenomenon. By tracking cell lineages and monitoring cell cycle progression through the classification of nuclear morphology, I have determined that mitotic wave slowdown can be explained through small but persistent delays in interphase. Comparing healthy embryos with embryos lacking terminal signalling pattern, I have shown that Erk patterning speeds interphase progression and pushes the wave origins towards the anterior and posterior poles. This work highlights how advances in high-throughput imaging and modern machine learning can yield new insights into developmental dynamics and enable data-driven modeling of complex biological systems.
Alexander Aulehla
Recorded on 8/13/2025
Title: Cycles upon cycles - collective rhythms during embryonic development
Abstract: We study the origin and function of collective signaling oscillations in embryonic development. Oscillatory signaling is linked to the sequential segmentation of the vertebrate embryo body axis and the formation of pre-vertebrae, somites. Most strikingly, signaling oscillations are coordinated between neighboring cells and result in spatio-temporal wave patterns that traverse the embryo periodically. I will discuss how we employ general synchronisation principles and experimental entrainment to reveal the hidden dynamical properties of this embryonic coupled oscillator network.
Connor Shrader
Recorded on 8/13/2025
Title: Modeling genetic drift and selection in spermatogonial stem cell dynamics
Abstract: Stem cells maintain and repair our tissues, but not all stem cells are identical. As organisms age, distinct stem cell "clones" can begin to dominate the cell population. While this behavior has been observed across multiple species and organs, the mechanisms and consequences of stem cell clonality are still poorly understood. We have developed a novel experimental approach using a CRISPR-Cas9 system to permanently “barcode” the spermatogonial stem cell clones in the zebrafish testis. Once these fish reach sexual maturity, we sample sperm each month to determine the contribution of each uniquely labeled stem cell clone to the sperm pool. To better understand the factors that drive clonal dynamics, we have also designed stochastic models of stem cell population dynamics. These models are formulated as hidden Markov models that describe rules for the division and differentiation of stem cells within the testis. We then use these models to quantify evidence of genetic drift and selection in our experimental data. Our models provide insight into how individual stem cell behavior can lead to population-level mosaicism and inform experimental efforts to verify our hypotheses.
Daniel Cruz
Recorded on 8/13/2025
Title: Personalizing agent-based models to construct medical digital twins
Abstract: Digital twin technology, originally developed for engineering, is being adapted to biomedicine and healthcare. A key challenge in this process is dynamically calibrating computational models to individual patients using data collected over time. This calibration is vital for improving model-based predictions and enabling personalized medicine. Biomedical models are often complex, incorporating multiple scales of biology and both stochastic and spatially heterogeneous elements. Agent-based models (ABMs), which simulate autonomous agents such as cells, are commonly used to capture how local interactions affect system-level behavior. However, no standard personalization methods exist for these models. The main challenge is bridging the gap between clinically measurable macrostates (e.g., blood pressure, heart rate) and the detailed microstate data (e.g., cellular processes) needed to run the model. In this work, we propose an algorithm that applies the ensemble Kalman filter (EnKF), a classic data assimilation technique, at the macrostate level. We then link the Kalman update at the macrostate to corresponding updates at the microstate level, ensuring that the resulting microstates are compatible with the desired macrostates and consistent with the model’s dynamics. This approach improves the personalization of complex biomedical models and enhances model-based forecasts for individual patients.
Sidney Holden
Recorded on 8/13/2025
Title: A continuum limit for dense spatial networks
Abstract: Many physical systems–such as dense neuronal or vascular networks and optical waveguide lattices– can be modeled by spatial networks, where slender “wires” (edges) support wave or diffusion equations subject to conservation conditions at nodes. We propose a continuum-limit framework which replaces edgewise differential equations with a coarse-grained partial differential equation (PDE) defined on the continuous space occupied by the network. The derivation naturally introduces an edge-conductivity tensor, an edge-capacity function, and a vertex number density to encode how each microscopic patch of the graph contributes to the macroscopic phenomena. We calculate all macroscopic parameters from first principles via a systematic discrete-to-continuous local homogenization, finding an anomalous effective embedding dimension resulting from a homogenized diffusivity. These high-density networks encode emergent material and functional properties. They reflect the ability of many real-world, space-filling networks to operate simultaneously at multiple scales–both at the system-wide and local levels–using the continuum as a feature. Numerical examples—including periodic lattices and random graphs (figure left, center)—demonstrate that each finite model converges to its corresponding PDE (posed on different manifolds like tori, disks, and spheres) in the limit of increasing vertex density. We expect our results to be useful in modelling biological network growth and function (e.g. quail embryo vasculature–figure right).
Carina Curto
Recorded on 8/13/2025
Title: Graphical domination and inhibitory control in recurrent networks
Abstract: Recurrent neural networks can be modeled as dynamical systems on directed graphs. What graph features are important for shaping the emergent dynamics? In this talk we will introduce the concept graphical domination and present key theorems about domination that help us understand the associated nonlinear dynamics. In particular, domination can be used to reduce graphs to smaller equivalent networks. We also show how reducible graph modules can be chained together to produce larger networks with predictable dynamics. These networks are amenable to control via inhibitory pulses. While these results were originally developed for a special class of effectively inhibitory threshold-linear networks, we will show how they apply equally well to E-I networks with global inhibition.