2026-11-02 | | Total: 27
Phylogenetic networks provide a natural framework for representing reticulate evolutionary events such as horizontal gene transfer and hybridization. In the present manuscript, we investigate the identifiability of phylogenetic networks under (hardwired) parsimony, a natural extension of Maximum Parsimony from trees to networks. More precisely, we ask which structural features of a network can be uniquely inferred from so-called $A_k$-alignments under this criterion. While identifiability has been studied extensively for phylogenetic trees under parsimony, little is known about the identifiability of networks under hardwired parsimony. We take a first step towards addressing this gap by developing a general framework for studying identifiability under hardwired parsimony. We then use this framework to show that the tree of blobs of a phylogenetic network can be reconstructed from its $A_1$-alignment. This result can be viewed as a natural extension of Buneman's well-known splits equivalence theorem. Finally, we show that, under hardwired parsimony, cycles of length at least $5$ cannot in general be uniquely identified from character data, establishing a fundamental limitation of hardwired parsimony.
Many biological interactions involve groups rather than pairs, as when several proteins assemble into a multi-protein complex. Simplicial complexes naturally capture these collective interactions. Existing hyperbolic embeddings represent network hierarchy and similarity through a few interpretable coordinates, but rely only on pairwise links: they cannot distinguish three units engaged in a common interaction from three units connected only in pairs. We introduce here the hyperbolic higher-order (HHO) embedding, which embeds a simplicial complex into the Poincaré disk, while incorporating its group interactions. The hyperbolic area of the triangle spanned by three interacting nodes then measures how compactly their interaction is represented. Testson synthetic and real-world networks, including a co-authorship network and the protein interaction network of \textit{Escherichia coli}, show that HHO reduces the hyperbolic area of triangles representing simplices relative to those obtained from pairwise network embeddings. HHO also recovers node groups more accurately than pairwise embeddings and a weighted-clique alternative, with performance competitive with Euclidean embeddings, while providing a low-dimensional representation that can be analyzed with the tools of network geometry.
Biochemical networks propagate the signals that regulate cellular behaviour, and which routes carry them depends on cellular context: a shift in context rewires the network's effective dynamics. Mechanistic simulators are misspecified here because the effective dynamics emerge from an interaction structure that is only partially known. Model discovery instead learns the dynamics from data, but sparse, correlated measurements poorly identify a full network model. What matters is narrower: the portion the context rewires. We introduce conditionally coarse-grained ODEs (Co$^2$ODE), which model that portion as a sparse neural or symbolic ODE: a conditional mask selects the markers carrying the rewiring, and a variational autoencoder compresses the remaining markers into a learnt closure term. On single-cell phosphoproteomic data spanning 30 cancer-relevant overexpressions, Co$^2$ODE compresses a 30-marker panel into a 10-dimensional state while matching the predictive accuracy of a neural ODE fit to the full panel. The closure absorbs what the mask drops, hence the accuracy is insensitive to which markers are kept. Adding a prior score to the loss shifts the selected markers at no cost in fit, indicating the data are not interventional enough to identify the underlying biology.
How stochastic rhythms in upstream regulatory signals are transmitted to downstream gene expression is a fundamental question in systems biology. Although experiments have shown that molecular degradation strongly influences downstream rhythmic responses, a general quantitative theory explaining how degradation, stochastic bursting, and upstream timing variability jointly govern rhythm transmission remains lacking. Here we develop a general analytical framework based on a martingale-based orthogonal decomposition, applicable to arbitrarily complex upstream dynamics, including memory effects. We derive exact autocorrelation and power-spectrum transmission laws, establish an oscillation-efficiency contraction principle, and identify a critical degradation rate that provides a necessary and sufficient condition for stochastic rhythm transmission. An explicit rhythm-transmission relation reveals that increasing burst frequency and modulation depth promotes rhythm transmission, whereas greater variability in upstream cycle duration suppresses it. In particular, increasing burst size leaves oscillation efficiency unchanged. The predicted degradation threshold correctly separates all ten pulsing and nonpulsing p53 target proteins reported in previous experiments. These results provide a theoretical foundation for understanding how cells selectively preserve or integrate rhythmic regulatory signals.
Direct-to-biology (D2B) screening enables rapid testing of newly synthesized molecules without purification. However, uncertain concentrations in crude mixtures and limited purified-compound measurements complicate the use of these data in machine learning. Here we show that models tailored to the structure of D2B data can improve prediction and guide the discovery of more potent and selective compounds. Our assay-aware models predict activity from single-point screens, relative potency across targets from crude dose-response assays, and absolute potency from purified-compound measurements. For crude assays, potency ratios between targets tested from the same well can retain selectivity information despite uncertain concentrations. A joint model that shares information across stages and balances their contributions to training improved predictions over separate neural networks. The largest gain was in crude dose-response selectivity prediction, where rank correlation increased from 0.27 to 0.40. In a separate experiment, independently trained models guided the selection of 84 new compounds to improve ROS1 potency while limiting activity against the anti-target TRKB. Synthesis and testing identified compounds that extended the measured Pareto front, achieving better combinations of potency and selectivity than those in the training data. The most selective new compound achieved a ROS1 IC50 of 14.7 nM and 233-fold selectivity over TRKB. By comparison, the most selective training compound had a ROS1 IC50 of 20.3 nM and 79-fold selectivity. ATLAS provides the experimental setting for this approach, supporting synthesis and screening of 1,000-10,000 model-selected compounds in two- to three-week cycles. Together, these results show how assay-aware learning can use the full D2B assay sequence to guide multi-objective molecular optimization.
Molecular mechanisms governing protein-protein interactions (PPI) contributing to the biogenesis of membraneless condensates remain poorly understood. Observations reported in this contribution demonstrate that liquid-liquid phase-separated (LLPS) condensate-forming proteins and their variants harboring condensate-disrupting human pathogenic mutations are significantly enriched among PPI cascade constituents of the Hybrid Assembly Pathway (HAP) of endogenous multiprotein complexes (EMC). In contrast, no structured or disordered proteins defined as condensate-negative molecules were detected among PPI cascade constituents. Cumulative evidence derived from different lines of analytical inquiries support the concept of the significant contributions of PPI cascade constituents of HAP to the biogenesis of intracellular pools of cell type-specific protein heterodimers assembled with the engagement of at least one protein monomer selected from condensate-forming protein molecules. Analytical evidence documented in this study support the concept that genome-wide chromosome-naïve networks of transposon-seeded DNA-templates nucleate the initiation of PPI cascade chain reactions with the capacity to facilitate the concurrent iterative assembly in human cells of both the membraneless biomolecular condensates and the structured multiprotein complexes. AlphaFold 3D modelling of mutation patterns and analysis of partner-bound experimental structures identified variants located at a measured protein dimerization interface. This contribution nominates numerous protein heterodimers of potential relevance to human pathophysiology for follow-up mechanistic studies of potential fundamental and translational impacts.
Rhythmic behaviours such as walking and running are commonly modelled as dynamical systems around a stable limit cycle. To keep such models tractable, it is usually assumed that each step depends only on the previous one, that the limit cycle is fixed, and that the sagittal and frontal planes can be studied separately. Here, we test these assumptions in seven public datasets of human gait (71 participants), including 16 people who both walked and ran. We fit linear step-to-step maps and validate them by cross-validation, holding out more data where more data exist. At first sight, gait has memory: the state two steps back lowers the prediction error by up to 29\%. Yet most of this temporal coupling vanishes once the operating point may drift slowly and measurement noise is modelled; gait remains limit-cycle-like, around an operating point that moves. Running keeps more memory than walking in the same people. Spatial coupling between the planes is present in most participants, mainly from sagittal to frontal, but weak. This supports planar models of healthy gait, while a repeatable, task-dependent coupling index may serve as a measure of gait quality. The simple models are wrong, but useful enough.
Summary statistics that measure the shapes of phylogenetic trees contributed to quantitatively understanding evolutionary processes. Among phylogenetic summary statistics, the tree balance indices have been designed to capture intuitive notions of balance or imbalance in trees. Here, we characterize the tree balance indices in terms of basis statistics, mapping the exponentially large space of rooted, binary, unlabeled trees to lower-dimensional spaces. We observe that many tree balance indices can be represented in terms of linear functions of these low-dimensional spaces, enabling us to compute exact formulas for the means, variances, and covariances of many balance indices under the Yule--Harding model. The new framework recapitulates existing formulas for the expectations of 8 indices: Sackin index, average leaf depth, total cophenetic index, $\hat{s}$-shape statistic, Colless index, corrected Colless index, quadratic Colless index, and rooted quartet index for binary trees. It produces new formulas for the expectations of 9 additional indices: equal weights Colless index, total $I$ index, total $I'$ index, Rogers $J$ index, total internal path length, total path length, average vertex depth, stairs1 index, and stairs2 index. Moreover, the framework also allows us to calculate the exact variance of all 17 indices that are represented by the bases as well as the covariance and correlation for all $\binom{17}{2} = 136$ pairs among the 17 indices.
Transcription factors (TFs) play an important role in regulating gene expression, making their identification a fundamental task in computational biology. In this work, we propose GraphTFactor, a novel end-to-end framework that models protein sequences as graphs and directly classifies transcription factors. By encoding sequences into graph structures, the proposed approach naturally handles variable-length inputs while preserving relevant relational information, without imposing length constraints. Experimental results show that GraphTFactor outperforms state-of-the-art methods on benchmark datasets, demonstrating the effectiveness of graph-based representations for TF prediction while providing interpretable insights into the sequence regions contributing to the model's predictions. To support reproducibility and practical adoption, we release an open-source software package and an interactive web application that let users easily perform TF prediction on custom sequences.
Functional connectivity studies have provided valuable insights into altered brain organization in autism spectrum disorder (ASD). However, most functional connectivity studies focus exclusively on either gray matter (GM) regions or white matter (WM) tracts and examine specific age groups, rather than investigating their functional interactions across brain development. This limits our understanding of how GM-WM functional interactions evolve developmentally and may obscure important aspects of brain network organization in ASD. We combined integrated WM+GM functional network connectivity (FNC) with brain-wide risk scores (BRSs) to quantify ASD-associated pattern expression in ABIDE I/II and an independent infant cohort. We employed a NeuroMark framework with predefined WM and GM intrinsic connectivity networks templates, to estimate subject-level FNC. First, we identified two distinct patterns of altered WM+GM FNC in ASD: GM-related hypoconnectivity centered on the paralimbic network and WM-GM hyperconnectivity centered on the frontal network. These findings highlight the involvement of higher-order cognitive networks and the importance of WM-GM functional interactions in characterizing ASD-related connectivity alterations. Second, infants with a higher familial likelihood of autism exhibited greater mean expression of ASD-associated FNC patterns than those with a lower familial likelihood, as quantified by BRS. These group-level findings do not establish developmental persistence or predict diagnostic outcomes at the individual level.
As climate change reshapes grape-growing conditions, predicting grapevine phenotypes from genomic variation is increasingly important for identifying and breeding cultivars suited to emerging environmental pressures. Existing research on the grapevine genotype-phenotype relationship using machine learning has often focused on individual traits, isolated regions, or a limited set of predictive models. We introduce VitisGP, a multi-phenotype curated dataset, integrating genome-wide SNP data with 26 heterogeneous phenotype prediction tasks spanning sommelier-curated aroma traits (N=11), botanical traits (N=14) and wine color---across 1,197 Vitis vinifera cultivars. We employ VitisGP to benchmark phenotype prediction, comparing machine-learning models with differing inductive biases---L2-regularized logistic regression, Random Forest, XGBoost, CatBoost, GBLUP and Transformer-based models, providing the first systematic comparison of this breadth on harmonized grapevine data. Across the benchmark, CatBoost achieves the strongest overall performance (macro-average AUROC=0.696 across 26 tasks), although model performance varies substantially across phenotype classes, with sensory and color traits showing stronger predictability than agronomic characteristics. VitisGP provides the community with a resource for investigating the genomic basis of phenotypic traits and screening grapevine diversity for climate-relevant characteristics.
Median genomes serve as candidates for ancestral reconstruction and help study the diversity of parsimonious rearrangement solutions. However, the breakpoint median problem is NP-hard for unichromosomal genomes, and medians are non-unique, making their diversity important to study. In this paper, we develop tree-based algorithms for studying breakpoint medians of linear unichromosomal genomes represented as unsigned permutations. The algorithms construct rooted labeled trees whose leaves encode candidate median genomes satisfying local adjacency constraints derived from the inputs. We show that the resulting reduced search spaces contain all breakpoint medians while being much smaller than the full permutation space. We investigate these reduced spaces experimentally for three-genome instances. Our method finds exact medians and many near-optimal approximations, while early sampling of leaves often produces near-medians quickly. We also introduce constrained sampling procedures for targeting balanced near-medians, whose distances to the input genomes are more evenly distributed and hence are particularly relevant for ancestral reconstruction. The resulting framework provides a unified method for exact enumeration, efficient sampling, and targeted exploration of median and near-median genomes. In contrast to Traveling Salesman Problem-based approaches, which typically return a single solution and do not reveal the diversity or balance of the median set, our method can enumerate all medians, sample near-medians early, and target structurally distinct candidates under the breakpoint distance.
Quantitative identification of molecular products and reaction pathways from scanning tunneling microscopy images is difficult when labeled images are unavailable for a new chemical system. We develop a simulation-assisted framework for zero-shot analysis that requires only precursor information for a new target system, with no experimental images or annotations from that system used during model training. The framework couples two functional components: NanoSim generates physics-grounded synthetic images from predicted reaction products and sampled surface scenes, while NanoAnalysis translates molecular detections into product populations, spatial arrangements, and mechanistic interpretations. Applied to a representative on-surface reaction, the detector achieves a micro-averaged F1 score of 0.92 on experimental images. Analysis of the detected molecules reveals annealing-dependent population changes and reaction-site preferences consistent with experimental observations. Beyond this capability, its modular organization allows heterogeneous data from other systems to be used to improve each component independently, providing a scalable route for continued generalization.
AlphaFold 3-like models (e.g., Boltz, Protenix, etc.) have pushed the frontier of protein binder design. This architecture is used both for the generation of binder candidates (e.g., by hallucination or converting these models into backbone or full atom generative models) as well as for the screening of these candidates for valid binders. In this paper, we start from the observation that AlphaFold 3's confidence metrics (ipSAE, ipTM) as well as structural consistency metrics can be predicted directly from the trunk even at early recycling. Because our proposed predictor is cheap, differentiable, and available at every cycle, it can be used in two distinct ways: as a decision signal to wisely spend the recycling budget for screening and as an objective for hallucination-based design. We show that both applications lead to increased performance at significantly reduced computational cost.
Chronic diseases evolve through long, heterogeneous trajectories shaped by interacting mediators, confounders, and treatment histories. Estimating long-term patient trajectories under alternative interventions is essential for in silico clinical trial design. However, most existing models overlook irregular visit intervals and cross-mediator dependencies, both of which are crucial for reliable long-term longitudinal prediction. Moreover, these models typically focus on outcome prediction rather than generating complete patient-state trajectories. To address these challenges, we propose CHDT, an intervention-conditioned Human Digital Twin generator under standard identification assumptions for multimodal longitudinal trajectory estimation. Its main empirical advantage arises from explicit modeling of irregular visit intervals, while cross-mediator refinement provides additional gains in selected intervention-conditioned settings. Across real-world ADNI and simulator-generated ALPACA cohorts, CHDT shows improved long-horizon patient-state fidelity relative to competing baselines. Ablation results indicate that irregular-time conditioning provides the primary empirical gain, and cross-mediator refinement provides additional improvements.
In traditional Chinese medicine (TCM), treatment follows syndrome differentiation (bianzheng), the step that maps a patient's findings to a TCM syndrome. A misdiagnosed syndrome therefore leads directly to inappropriate treatment, and the error usually surfaces only after that treatment fails. The physician who sees the patient next must diagnose a condition already altered by the wrong therapy, identify the initial physician's error, and avoid continuing the failed therapeutic direction. We introduce TCM-MisDx, an interactive benchmark of 820 consultations after treatment failure built from three published monographs of TCM misdiagnosis cases. An agent gathers symptoms, tongue and pulse findings through the four diagnostic methods before reporting the current syndrome, the initial error, and a corrective treatment plan. We further present TCM-SEA, a self-evolving consultation agent that rewrites its own clinical skill document: failed consultations are compared with the expert physician's documented analysis, the agent proposes a minimal edit, and the edit is kept only if it passes an acceptance gate on held-in and held-out cases. The underlying language model is never updated. On the sealed test set, TCM-SEA raises the composite score from 61.8 to 65.5, syndrome accuracy from 53.9 to 57.2 and treatment-principle accuracy from 44.7 to 48.9, and lowers the therapeutic momentum rate, the fraction of corrective plans that keep the initial physician's failed direction, from 44.1 to 38.2. Learning from expert cases outperforms learning from outcome feedback alone, with the largest gains on the monograph of inappropriate-treatment cases. The two accepted edits are short, human-readable clinical checklists that physicians can audit directly.
HER2 status decides who receives HER2-targeted therapy. Pathologists score HER2 protein by immunohistochemistry (IHC); equivocal IHC 2+ tumours, their hardest calls, wait for in situ hybridisation (ISH). We built a HER2 reader for routine H&E slides, owned end to end, as a baseline for pathologists. A 22M-parameter ViT-S/16 tile encoder was trained on public TCGA slides by self-supervision, distillation from H-optimus-0 and two rounds of JEPA refinement; five cross-validated attention-based slide classifiers average their probabilities; and an operating point is fixed on training data. Compared under one paired ten-seed protocol, readers on our encoder rank first of 7 on external validation from TCGA-BRCA to Yale-HER2 (AUC 0.812), level with readers on encoders 4 and 50 times larger. We also rebuilt PathBench's encoder benchmark on public cohorts; our replica reproduces its published order of reference encoders, and under it our encoder ranks first of 7 on TCGA-BRCA (0.689) and sits level with H-optimus-0 and Phikon on Yale-HER2 (0.897). On three external hospitals our deployed reader reaches AUC 0.82 on Yale-HER2 (about 0.81 with training slides balanced for tissue size), 0.75 on HEROHE and 0.63 on BCNB biopsies, and among HEROHE's 85 IHC 2+ cases it ranks ISH-amplified above non-amplified tumours at 0.76, best among our TCGA-trained readers. Per-laboratory score normalisation and a laboratory's own slides improve it further on a new site. Without a GPU our reader needs 5-6 s per slide. It serves as a research-use pre-read in a pathology platform, where it proposes a probability and a call beside IHC and helps prioritise ISH for equivocal cases. Weights, per-slide predictions and the data behind every figure are public at https://huggingface.co/cloudkites/her2-reader.
Objective: To distinguish the transferable predictive value of physiological trajectories from that of observation patterns in ICU mortality prediction, and to assess whether measurement timing affects performance independently of observation quantity. Materials and Methods: In-hospital mortality was predicted at 48 hours using MIMIC-III, MIMIC-IV, and eICU. Nested representations progressively integrated latest-value snapshots (R1), longitudinal physiological summaries (R2), and measurement frequency and timing features (R3). Models fitted on 7,762 MetaVision patients were evaluated on 21,682 MIMIC-IV patients and 43,878 eICU patients across 204 hospitals. The effect of measurement timing was assessed while holding the number of retained observations constant. Results: Physiological trajectories improved external prediction in eICU (R2-R1 log-loss difference -0.00922, 95% CI -0.01148 to -0.00704). Observation-pattern features improved locally developed eICU models (R3-R2 -0.00257, 95% CI -0.00370 to -0.00146), but not clearly after transfer from MetaVision (-0.00012, 95% CI -0.00206 to +0.00208). R3 was favoured in 10 of 30 source bootstrap refits. At matched observation counts, retaining measurements spread throughout the window reduced log loss relative to random retention by 0.03280 for R3 versus 0.00703 for R2. Removing the final 12 hours increased log loss relative to count-matched random removal by 0.01947 for R3 versus 0.00822 for R2. Under reduced availability, the local R3 increment changed from beneficial to harmful. Discussion: Physiological trajectories transferred more consistently, while observation-pattern features were particularly sensitive to the timing and availability of measurements. Conclusion: Clinical prediction models should be evaluated not only in new patient populations but also under the measurement conditions that supply their inputs.
Acoustic indices are widely used as rapid proxies for biodiversity, but they summarise the distribution of sound energy rather than the species producing it and often track richness poorly. We revise and test an Acoustic Biodiversity Index (ABI) built from species-level detections. Its primary form, ABI_E, is Shannon diversity of detections normalised by the regional species pool expected at the site and season (H / ln S_ref); temporal stability T is reported as a separate component. We first benchmarked four classifiers on a strongly labelled tropical dataset (1,216 calls, 24 classes) with recording-grouped cross-validation. BirdNET embeddings + logistic regression performed best (macro-F1 0.750 +/- 0.029), and random frame-level splits inflated accuracy from 0.41 to 0.70. We then applied BirdNET v2.4 with a location and season filter to recordings from a riverine savanna in the Zambezi Valley, Zimbabwe, and from three habitats at Emirates Park Zoo and Resort, Abu Dhabi, United Arab Emirates. A stricter post-hoc location filter removed 517 high-confidence detections of 9 implausible species, all at the zoo. After diel-matched rarefaction to equal effort (6 h), the savanna site held 21 species (Hill number 1D = 9.83; ABI_E = 0.38) against 1.3-13 species (ABI_E = 0.03-0.21) at the zoo sites. Across 149 recorded hours, ABI_E correlated strongly with richness (Spearman rho = 0.93), whereas six standard indices did not (best: ACI, rho = 0.48). The stability term was lower at the savanna site (0.39) than at the zoo (0.65-0.67), showing that it penalises the natural diel cycle of bird activity. Species-level indices with explicit false-positive control, a declared regional reference pool and effort standardisation provide a more interpretable biodiversity indicator than energy-based acoustic indices.
At low Reynolds number, microswimmers are coupled by the flows they create, even when they are apart. The best way for each swimmer to push itself forward then depends on its neighbors, turning swimming together into a game. To study this, we optimize the shape changes of soft swimmers, alone and in pairs, with a differentiable Stokes solver. From random initial membrane motions, a lone swimmer converges to a traveling membrane wave. In a pair, the stroke of a swimmer affects its neighbor far more than the neighbor affects its own best stroke. For swimmers in line, this influence decays as $1/\ell$ with their separation $\ell$, whereas the incentive to respond decays as $1/\ell^2$. Selfish swimmers, which maximize their own payoff, respond only to the weak incentive. The cooperative gain they miss decays only as $1/\ell^2$, and their selfish adaptations leave the pair worse off than if both had kept their solitary strokes. Related swimmers, which also value the payoff of their neighbor, recover this gain, whereas competitors pay to hinder each other, both following a law that we derive analytically. Hydrodynamics alone thus makes swimming in company a social dilemma, whose outcome depends on how much the swimmers value each other.
Analysis of mathematical models in ecology and epidemiology often emphasizes asymptotic dynamics, such as stable equilibria, periodic orbits, and attractors. However, many biologically relevant trajectories exhibit long transient dynamics in which a relevant observable evolves slowly for an extended period before a comparatively rapid transition occurs. On finite time scales relevant to data collection, management, and intervention, such transient behavior can be more informative than the eventual asymptotic outcome. In this work, we study long-lasting and slowly varying transient dynamics in continuous-time systems through the framework of nonequilibrium transient centers. These centers are points in the state space near which the observable can remain nearly constant for a long time, even when the full state is not near an equilibrium. We develop finite-time escape criteria for identifying such centers and their reachable counterparts. These criteria provide practical conditions under which nearby trajectories can remain in a slow regime for a long period of time before eventually leaving it. We apply the theory to an epidemic model with delay of stability loss. In this model, trajectories may remain close to the disease-free manifold for a long time before the infectious population resurges. The framework developed here describes this behavior as slow evolution near a nonequilibrium transient center followed by finite-time escape from the slow regime.
Morphogenesis of diatom silica shells involves a sequence of topological transitions: initially separate ribs grow and branch, then merge through bridges, leaving a regular network of pores. Existing models based on diffusion-limited aggregation or reaction-diffusion reproduce either branching or pore formation, but not both within a single framework. Here we propose a two-stage model of diffusion-driven morphogenesis with evolving connectivity. The model is formulated on a fixed finite element mesh with a density field on an active layer. Stage 1 reproduces dendritic growth, where a stimulus based on the local diffusive flux exceeds a threshold. Stage 2 uses an inverted stimulus, where high local flux suppresses precipitation, stabilizing pores between the dendrites. The update rules are analogous to sensitivity in density-based topology optimization and to strain energy density in bone remodeling. The model predicts distinct morphological regimes on each stage, controlled by the respective thresholds. For a reference calculation, the predicted pore network is regular, with approximately uniform spacing, and pores are absent near the domain edges, consistent with observations on diatom valves. The model provides a unified framework for describing both stages of silica morphogenesis and suggests a route for comparing computational predictions with cryo-electron tomography data.
Designing regulatory DNA sequences for targeted gene-expression control is an important challenge. Recent approaches frame this task as property-conditioned generation, where diffusion models trained on natural sequences are guided toward high predicted activity through reward-based fine-tuning or sampling-time guidance. We show that single-reward evaluation can obscure reward-hacking behavior. By decomposing evaluation into property, novelty, and plausibility, we identify an empirically underexplored region of the evaluation space in which all three objectives are satisfied, but no existing baseline succeeds. We postulate that this gap arises because prior methods treat guidance primarily as a matter of magnitude. In regulatory DNA, however, sequence grammar is spatially nonuniform: motif positions are sensitive to perturbation, while background regions allow more exploration. We therefore propose LOCOGEN by reformulating guidance as a problem of locality. LOCOGEN learns an invariant motif mask through sequence-level invariance learning to protect motif positions, while applying OOD damping only to background positions. It operates at sampling time on a frozen diffusion backbone, without retraining the backbone. Across HepG2, K562, and SK-N-SH, LOCOGEN reaches this previously underexplored region of the property-novelty-plausibility space. Code is available at https://github.com/Lotfollahi-lab/LocoGen.
We study, after Coron Le Jan (2024), the genetic contribution of a single strongly advantageous mutant in a biparental Moran model under infinite selection. We prove a functional limit theorem for the entire beneficial mutation spreading period: after normalization by \(\sqrt N\), the contributions carried by the advantaged and disadvantaged subpopulations converge jointly in law to a non-degenerate random multiple of an explicit deterministic profile. Randomness is generated during the initial microscopic growth of the mutant lineage and subsequently propagated through the population along a deterministic macroscopic trajectory. At the heart of this phenomenon is an eigenmode martingale that measures the imbalance between the mean genetic weight of the founding ancestor within the two classes; its limiting value provides precisely the random amplitude of the macroscopic profile and is characterized by a smoothing equation.
Understanding which regulatory sequence motifs a genomic deep learning model has learned, with statistical error guarantees, is essential for trusting its use in downstream tasks such as variant effect prediction. Existing approaches to identifying which units in a trained network encode human- understandable concepts fall into two groups. Exploratory methods rank and visualize units but offer no statistical guarantee that the selected units are associated with the prediction target across the population of inputs. Testing-based methods test only single units or one pre-defined set of them, and do not control the error across the many sets of units. We introduce CUTTIE (Cluster-tree Unit Top-down Testing of Independence with Error control), which clusters the units of each layer into a tree and tests the clusters top-down with conditional independence tests, localizing information about the target while controlling the family-wise error rate. Thereby, we identify clusters that carry information about the target beyond the rest of the layer. For the units in the identified clusters, conditional attribution maps then reveal what motifs they encode. In a controlled simulation, the procedure identifies the clusters whose concepts affect the target. We apply the procedure to a convolutional model predicting cell-type-specific chromatin accessibility from DNA sequence, where the concepts encoded in the identified clusters map onto known regulatory motifs.