q-bio updates on arXiv.org
q-bio updates on the arXiv.org e-print archive.
Vision Foundry: A System for Training Foundational Vision AI Models
oai:arXiv.org:2512.11837v1
arXiv:2512.11837v1 Announce Type: new
Abstract: Self-supervised learning (SSL) leverages vast unannotated medical datasets, yet steep technical barriers limit adoption by clinical researchers. We introduce Vision Foundry, a code-free, HIPAA-compliant platform that democratizes pre-training, adaptation, and deployment of foundational vision models. The system integrates the DINO-MX framework, abstracting distributed infrastructure complexities while implementing specialized strategies like Magnification-Aware Distillation (MAD) and Parameter-Efficient Fine-Tuning (PEFT). We validate the platform across domains, including neuropathology segmentation, lung cellularity estimation, and coronary calcium scoring. Our experiments demonstrate that models trained via Vision Foundry significantly outperform generic baselines in segmentation fidelity and regression accuracy, while exhibiting robust zero-shot generalization across imaging protocols. By bridging the gap between advanced representation learning and practical application, Vision Foundry enables domain experts to develop state-of-the-art clinical AI tools with minimal annotation overhead, shifting focus from engineering optimization to clinical discovery.
Gene regulatory network inference algorithm based on spectral signed directed graph convolution
oai:arXiv.org:2512.11927v1
arXiv:2512.11927v1 Announce Type: new
Abstract: Accurately reconstructing Gene Regulatory Networks (GRNs) is crucial for understanding gene functions and disease mechanisms. Single-cell RNA sequencing (scRNA-seq) technology provides vast data for computational GRN reconstruction. Since GRNs are ideally modeled as signed directed graphs to capture activation/inhibition relationships, the most intuitive and reasonable approach is to design feature extractors based on the topological structure of GRNs to extract structural features, then combine them with biological characteristics for research. However, traditional spectral graph convolution struggles with this representation. Thus, we propose MSGRNLink, a novel framework that explicitly models GRNs as signed directed graphs and employs magnetic signed Laplacian convolution. Experiments across simulated and real datasets demonstrate that MSGRNLink outperforms all baseline models in AUROC. Parameter sensitivity analysis and ablation studies confirmed its robustness and the importance of each module. In a bladder cancer case study, MSGRNLink predicted more known edges and edge signs than benchmark models, further validating its biological relevance.
GNN-Based Deep Surrogate Modeling of Knee Contact Mechanics: Generalizing Neuromuscular Control Patterns Across Subjects
oai:arXiv.org:2512.11936v1
arXiv:2512.11936v1 Announce Type: new
Abstract: Background: Accumulation of abnormal contact stress is a primary biomechanical driver of acute meniscal tears and chronic osteoarthritis. While Finite Element Analysis (FEA) provides the necessary fidelity to quantify these injury-inducing loads, its high computational cost precludes clinical utility. Emerging deep surrogate models promise real-time assessment but suffer a critical blind spot: they predominantly focus on learning anatomical variations, largely overlooking the neuromuscular control patterns. These dynamic, subject-specific motor strategies fundamentally dictate potentially injurious stress distributions inside the knee. Methods: This study investigates the generalization capability of the topology-aware MeshGraphNet regarding cross-subject neuromuscular control patterns under fixed anatomical conditions. We constructed a dataset using gait data from nine subjects via an OpenSim-FEBio co-simulation platform. The MGN was compared against a structure-agnostic Node-wise MLP using a rigorous grouped 3-fold cross-validation on unseen subjects. Results: The MGN demonstrated superior fidelity, achieving a correlation of 0.94 with ground truth (vs. 0.88 for MLP). In contrast to the MLP, which exhibited the "peak shaving" defect common in deep learning, MGN significantly reduced peak-stress prediction errors and achieved higher spatial overlap in high-risk regions. This indicates that MGN effectively captured the non-local force-transmission pathways unique to each subject's movement strategy. Conclusion: By mimicking the propagation of physical stress through message passing, MGN successfully decodes the heterogeneity of human neuromuscular control, even under fixed anatomy. This establishes GNNs as robust clinical tools capable of identifying functional injury risks that are invisible to purely geometry-based surrogate models.
An algorithm to align a chain of sequences to paths in a pangenome graph
oai:arXiv.org:2512.12052v1
arXiv:2512.12052v1 Announce Type: new
Abstract: Affordable, high-quality whole-genome assemblies have made it possible to construct rich pangenomes that capture haplotype diversity across many species. As these datasets grow, they motivate the development of specialized techniques capable of handling the dense sequence variation found in large groups of related genomes. A common strategy is to encode pangenomic information in graph form, which provides a flexible substrate for improving algorithms in areas such as alignment, visualization, and functional analysis. Methods built on these graph models have already shown clear advantages in core bioinformatics workflows, including read mapping, variant discovery, and genotyping. By integrating multiple sequence and coordinate representations into a single structure, pangenome graphs offer a unified and expressive framework for comparative genomics. Although it remains unclear whether graph-based references will ultimately supplant traditional linear genomes, their versatility ensures that they will play a central role in emerging pangenomic approaches. This paper introduces an algorithm to mine a chain of sequences in pangenome graphs that might be useful in the functional analysis of pangenome graphs. Specifically, the algorithm calculates all maximal paths in a pangenome graph aligning with a given chain of sequences in the segments of the path vertices, possibly with some maximal gap as specified by the user.
Prediction of PLX-4720 Sensitivity in Cancer Cell Lines through Multi-Omics Integration and Attention-Based Fusion Modeling
oai:arXiv.org:2512.12113v1
arXiv:2512.12113v1 Announce Type: new
Abstract: Predicting the sensitivity of cancer cell lines to PLX-4720, a preclinical BRAF inhibitor, requires models capable of capturing the multilayered regulation of oncogenic signaling. Single-omics predictors are often insufficient because drug response is shaped by interactions among genomic alterations, epigenetic regulation, transcriptional activity, protein signaling, metabolic state, and network-level context. In this study we develop an attention-based multi-omics integration framework using genomic, epigenomic, transcriptomic, proteomic, metabolomic, and protein interaction data from the GDSC1 panel. Each modality is encoded into a latent representation using feed-forward neural networks or graph convolutional networks, and fused through an attention mechanism that assigns modality-specific importance weights. A regression model is then used to predict PLX-4720 response. Across single- and multi-omics configurations, the best performance is achieved by integrating genomics and transcriptomics, which yields validation R2 values above 0.92. This reflects the complementary roles of mutational status and downstream transcriptional activation in shaping sensitivity to BRAF inhibition. Epigenomics is the strongest single-omics predictor, while metabolomics and PPI data contribute additional context when combined with other modalities. Integration of three to five omics layers improves stability but does not surpass the accuracy of the best two-modality combinations, likely due to information redundancy and sample-size imbalance. These findings highlight the importance of modality selection rather than maximal data depth. The proposed framework provides an efficient and biologically grounded strategy for drug response prediction and supports the development of precision pharmacogenomics.
Modeling Dabrafenib Response Using Multi-Omics Modality Fusion and Protein Network Embeddings Based on Graph Convolutional Networks
oai:arXiv.org:2512.12134v1
arXiv:2512.12134v1 Announce Type: new
Abstract: Cancer cell response to targeted therapy arises from complex molecular interactions, making single omics insufficient for accurate prediction. This study develops a model to predict Dabrafenib sensitivity by integrating multiple omics layers (genomics, transcriptomics, proteomics, epigenomics, and metabolomics) with protein network embeddings generated using Graph Convolutional Networks (GCN). Each modality is encoded into low dimensional representations through neural network preprocessing. Protein interaction information from STRING is incorporated using GCN to capture biological topology. An attention based fusion mechanism assigns adaptive weights to each modality according to its relevance. Using GDSC cancer cell line data, the model shows that selective integration of two modalities, especially proteomics and transcriptomics, achieves the best test performance (R2 around 0.96), outperforming all single omics and full multimodal settings. Genomic and epigenomic data were less informative, while proteomic and transcriptomic layers provided stronger phenotypic signals related to MAPK inhibitor activity. These results show that attention guided multi omics fusion combined with GCN improves drug response prediction and reveals complementary molecular determinants of Dabrafenib sensitivity. The approach offers a promising computational framework for precision oncology and predictive modeling of targeted therapies.
Prevalence of Upper Extremity Distal Predominant Weakness Pattern in Chronic Stroke
oai:arXiv.org:2512.12147v1
arXiv:2512.12147v1 Announce Type: new
Abstract: Background: Hemiparesis after subcortical stroke is classically described as distal upper-extremity (UE) predominant, but prevalence data in chronic stroke is limited. Objective: Determine the prevalence of distal predominant UE weakness in exclusively subcortical chronic stroke versus other stroke distributions, characterize cohort differences, and describe UE weakness patterns in chronic stroke overall. Methods: Outpatient records were retrospectively reviewed to identify chronic stroke subjects. Lesion locations were classified from radiographic reports as exclusively subcortical or not (using a whole brain and supratentorial definition). UE weakness was categorized as distal predominant or not. Prevalence was compared with $\chi$-squared testing and odds ratios (OR). Results: 250 subjects were included (mean 861 days post-stroke). Using the whole-brain definition, distal predominant weakness occurred in 30.6% of exclusively subcortical versus 17.4% of non-exclusively subcortical strokes (OR 2.09, 95% CI 1.15-3.81; p=0.014). Using the supratentorial definition, distal predominant weakness occurred in 27.9% versus 17.9%, respectively (OR 2.16, 95% CI 1.17-3.96; p=0.012). Across all chronic strokes, 60% had no UE weakness; distal predominant weakness was the most common weakness pattern (23%), followed by uniform UE weakness (12%); proximal predominant weakness was rare (3%). Conclusions: Distal predominant UE weakness is more prevalent in chronic exclusively subcortical stroke than in non-subcortical stroke. These prevalence estimates may help predict long-term outcomes based on lesion location, support rehabilitation planning, and aid clinical lesion localization and research prioritization.
DCAF-Net: Dual-Channel Attentive Fusion Network for Lower Limb Motion Intention Prediction in Stroke Rehabilitation Exoskeletons
oai:arXiv.org:2512.12184v1
arXiv:2512.12184v1 Announce Type: new
Abstract: Rehabilitation exoskeletons have shown promising results in promoting recovery for stroke patients. Accurately and timely identifying the motion intentions of patients is a critical challenge in enhancing active participation during lower limb exoskeleton-assisted rehabilitation training. This paper proposes a Dual-Channel Attentive Fusion Network (DCAF-Net) that synergistically integrates pre-movement surface electromyography (sEMG) and inertial measurement unit (IMU) data for lower limb intention prediction in stroke patients. First, a dual-channel adaptive channel attention module is designed to extract discriminative features from 48 time-domain and frequency-domain features derived from bilateral gastrocnemius sEMG signals. Second, an IMU encoder combining convolutional neural network (CNN) and attention-based long short-term memory (attention-LSTM) layers is designed to decode temporal-spatial movement patterns. Third, the sEMG and IMU features are fused through concatenation to enable accurate recognition of motion intention. Extensive experiment on 11 participants (8 stroke subjects and 3 healthy subjects) demonstrate the effectiveness of DCAF-Net. It achieved a prediction accuracies of 97.19% for patients and 93.56% for healthy subjects. This study provides a viable solution for implementing intention-driven human-in-the-loop assistance control in clinical rehabilitation robotics.
Scalable branch-and-bound model selection with non-monotonic criteria including AIC, BIC and Mallows's $\mathit{C_p}$
oai:arXiv.org:2512.12221v1
arXiv:2512.12221v1 Announce Type: new
Abstract: Model selection is a pivotal process in the quantitative sciences, where researchers must navigate between numerous candidate models of varying complexity. Traditional information criteria, such as the corrected Akaike Information Criterion (AICc), Bayesian Information Criterion (BIC), and Mallows's $\mathit{C_p}$, are valuable tools for identifying optimal models. However, the exponential increase in candidate models with each additional model parameter renders the evaluation of these criteria for all models -- a strategy known as exhaustive, or brute-force, searches -- computationally prohibitive. Consequently, heuristic approaches like stepwise regression are commonly employed, albeit without guarantees of finding the globally-optimal model.
In this study, we challenge the prevailing notion that non-monotonicity in information criteria precludes bounds on the search space. We introduce a simple but novel bound that enables the development of branch-and-bound algorithms tailored for these non-monotonic functions. We demonstrate that our approach guarantees identification of the optimal model(s) across diverse model classes, sizes, and applications, often with orders of magnitude computational speedups. For instance, in one previously-published model selection task involving $2^{32}$ (approximately 4 billion) candidate models, our method achieves a computational speedup exceeding 6,000. These findings have broad implications for the scalability and effectiveness of model selection in complex scientific domains.
Advancements in Hematology Analyzers: Next-Generation Technologies for Precision Diagnostics and Personalized Medicine
oai:arXiv.org:2512.12248v1
arXiv:2512.12248v1 Announce Type: new
Abstract: Hematology analyzers are essential diagnostic and monitoring tools for detecting blood diseases. Although contemporary analyzers produce only basic insights, they are often not as detailed as required under the personalized medicine paradigm. Next-Generation Hematology Analyzers (NGHAs) are revolutionary newcomers in the field, with significant advantages over regular hematology analyzers. They provide deeper insights into cellular morphology, function, and genetic profiles. This detailed information opens up possibilities for tailor-made diagnostic and therapeutic approaches in precision medicine. This review presents some revolutionary technologies that have changed hematology analyzers and provides an overview of their limitations, basic functions, and influence on clinical practice. It focuses on the integration of state-of-the-art technologies, such as microfluidics, advanced optics, artificial intelligence, flow cytometry, and digital imaging, empowering NGHAs to improve diagnostic accuracy, rapidly detect diseases, and support flexible, targeted therapy. Hints regarding point-of-care hematology testing are also provided to discuss its implications for transforming healthcare patterns. This review highlights the data management, standardization, regulatory, and ethical challenges associated with these technologies. A review tracking the current state-of-the-art and trends for the future is provided to show how these advancements may reconfigure hematology analyzer design and act as a stepping stone for future therapeutic reforms.
Accurate de novo sequencing of the modified proteome with OmniNovo
oai:arXiv.org:2512.12272v1
arXiv:2512.12272v1 Announce Type: new
Abstract: Post-translational modifications (PTMs) serve as a dynamic chemical language regulating protein function, yet current proteomic methods remain blind to a vast portion of the modified proteome. Standard database search algorithms suffer from a combinatorial explosion of search spaces, limiting the identification of uncharacterized or complex modifications. Here we introduce OmniNovo, a unified deep learning framework for reference-free sequencing of unmodified and modified peptides directly from tandem mass spectra. Unlike existing tools restricted to specific modification types, OmniNovo learns universal fragmentation rules to decipher diverse PTMs within a single coherent model. By integrating a mass-constrained decoding algorithm with rigorous false discovery rate estimation, OmniNovo achieves state-of-the-art accuracy, identifying 51\% more peptides than standard approaches at a 1\% false discovery rate. Crucially, the model generalizes to biological sites unseen during training, illuminating the dark matter of the proteome and enabling unbiased comprehensive analysis of cellular regulation.
Modeling the Prey-Predator Dynamics of Habu Snakes and Mongooses Leading to Ecological Disaster on Amami Oshima Island in Japan
oai:arXiv.org:2512.12388v1
arXiv:2512.12388v1 Announce Type: new
Abstract: The introduction of mongooses from Indian subcontinent to Amami Oshima Island, Japan, aimed at controlling the population of venomous Habu snakes, has led to significant ecological disruptions, raising concerns about the long-term sustainability of the islands biodiversity. To highlight the unintended consequences of such interventions and the necessity of understanding predator-prey dynamics in preserving ecological balance, a mathematical model incorporating snake, mongooses, mouse and natural resources has been proposed to explore their role in the ongoing ecological disaster and analysis the other scenarios if the authorities applied different approaches in place of already implemented strategy. Determining the model's existence and uniqueness, stability at equilibrium points, and state variable characteristics are some of the parts of the analytical analysis of the model. Additionally, sensitivity analysis is conducted to identify sensitive factors. In addition, the Runge-Kutta 4th order has been used to execute the numerical simulations. Our research reveals that although the government began killing and trapping mongooses almost 20 years after their introduction, but if trapping had started just 10 years after their introduction, the outcome could have been drastically different. This time, mongooses would not have been extinct, and their coexistence with other native species would have helped to preserve the ecological balance and prevent the severe ecological damage that is presently being seen. Thus, It is recommended to use mathematical modeling to explore alternatives before decision-making, ensuring sustainable ecosystem management while preventing irreversible impacts of invasive species.
Reduced rank regression for neural communication: a tutorial for neuroscientists
oai:arXiv.org:2512.12467v1
arXiv:2512.12467v1 Announce Type: new
Abstract: Reduced rank regression (RRR) is a statistical method for finding a low-dimensional linear mapping between a set of high-dimensional inputs and outputs. In recent years, RRR has found numerous applications in neuroscience, in particular for identifying "communication subspaces" governing the interactions between brain regions. This tutorial article seeks to provide an introduction to RRR and its mathematical foundations, with a particular emphasis on neural communication. We discuss RRR's relationship to alternate dimensionality reduction techniques such as singular value decomposition (SVD), principal components analysis (PCA), principal components regression (PCR), and canonical correlation analysis (CCA). We also derive important extensions to RRR, including ridge regularization and non-spherical noise. Finally, we introduce new metrics for quantifying communication strength as well as the alignment between communication axes and the principal modes of neural activity. By the end of this article, readers should have a clear understanding of RRR and the practical considerations involved in applying it to their own data.
Dual-Model Framework for CHIKV Transmission Modeling: ODE and Petri Net Analysis of the 2025 Foshan Outbreak
oai:arXiv.org:2512.12577v1
arXiv:2512.12577v1 Announce Type: new
Abstract: This study constructs a dual-model framework integrating Ordinary Differential Equations (ODE) and Petri Nets (PN) to analyze the 2025 Chikungunya outbreak in Foshan City, China. We employ SEICR compartmental modeling to compare two distinct approaches under identical epidemiological scenarios and evaluate intervention effectiveness through three-phase fitting protocols. Both models demonstrate excellent accuracy with MAE of 18.77-18.91 cases and RMSE of 36.52-36.54 cases. Models predicted epidemic peaks at day 32 (406 cases), 3 days earlier than observed (day 35, 432 cases), with 6.0% peak value error. Reproduction number analysis revealed initial R0 of 14.67 (ODE)/13.90 (PN), with effective reproduction numbers decreasing through intervention phases: 7.85/7.86 after Phase 1, 7.59/7.56 after Phase 2, and 0.059 in Phase 3, achieving transmission blockade. Sensitivity analysis showed recovery rate) as the most sensitive parameter (Sobol index 0.9672), explaining 96.72% of R0 variation. This study presents the first systematic ODE-Petri Net comparison, providing a novel dual-model framework for vector-borne disease modeling with significant theoretical and practical value for epidemic control strategy formulation.
Investigating High-Order Behaviors in Multivariate Cardiovascular Interactions via Nonlinear Prediction and Information-Theoretic Tools
oai:arXiv.org:2512.12709v1
arXiv:2512.12709v1 Announce Type: new
Abstract: Assessing the synergistic high-order behaviors (HOBs) that emerge from underlying structural mechanisms is crucial to characterize complex systems. This work leverages the combined use of predictability and information measures to detect and quantify HOBs in synthetic and physiological network systems. After providing formal definitions of mechanisms and behaviors in a complex system, measures of statistical synergy are defined as the whole-minus-sum excess of mutual predictability ($\Delta_\textrm{MP}$) or mutual information ($\Delta_\textrm{MI}$) obtained when considering the system as a whole rather than as a combination of its units. The two measures are computed using model-free methods based on nonlinear prediction and entropy estimation. The application to simulated linear Gaussian systems and nonlinear deterministic and stochastic dynamic systems shows that $\Delta_\textrm{MP}$ tends to vanish for target variables influenced by additive effects of single independent source variables and is positive in the presence of group interactions between sources, while $\Delta_\textrm{MI}$ exhibits a higher propensity to display positive values. The analysis of physiological variables shows significant values of $\Delta_\textrm{MI}$ when investigating the additive effect of systolic and diastolic arterial pressure on mean arterial pressure, and of both $\Delta_\textrm{MP}$ and $\Delta_\textrm{MI}$ when assessing how diastolic pressure is modulated by pre-ejection and left-ventricular ejection times. HOBs can be more clearly identified by information-theoretic measures, while prediction measures are more sensitive to synergy arising from the governing rules of the system analyzed rather than from pure statistical dependencies. Quantifying HOBs through measures sensitive to structural mechanisms can provide biomarkers to assess physio-pathological alterations of cardiovascular networks.
Random matrix theory of sparse neuronal networks with heterogeneous timescales
oai:arXiv.org:2512.12767v1
arXiv:2512.12767v1 Announce Type: new
Abstract: Training recurrent neuronal networks consisting of excitatory (E) and inhibitory (I) units with additive noise for working memory computation slows and diversifies inhibitory timescales, leading to improved task performance that is attributed to emergent marginally stable equilibria [PNAS 122 (2025) e2316745122]. Yet the link between trained network characteristics and their roles in shaping desirable dynamical landscapes remains unexplored. Here, we investigate the Jacobian matrices describing the dynamics near these equilibria and show that they are sparse, non-Hermitian rectangular-block matrices modified by heterogeneous synaptic decay timescales and activation-function gains. We specify a random matrix ensemble that faithfully captures the spectra of trained Jacobian matrices, arising from the inhibitory core - excitatory periphery network motif (pruned E weights, broadly distributed I weights) observed post-training. An analytic theory of this ensemble is developed using statistical field theory methods: a Hermitized resolvent representation of the spectral density processed with a supersymmetry-based treatment in the style of Fyodorov and Mirlin. In this manner, an analytic description of the spectral edge is obtained, relating statistical parameters of the Jacobians (sparsity, weight variances, E/I ratio, and the distributions of timescales and gains) to near-critical features of the equilibria essential for robust working memory computation.
A Disproof of Large Language Model Consciousness: The Necessity of Continual Learning for Consciousness
oai:arXiv.org:2512.12802v1
arXiv:2512.12802v1 Announce Type: new
Abstract: The requirements for a falsifiable and non-trivial theory of consciousness significantly constrain such theories. Specifically, recent research on the Unfolding Argument and the Substitution Argument has given us formal tools to analyze requirements for a theory of consciousness. I show via a new Proximity Argument that these requirements especially constrain the potential consciousness of contemporary Large Language Models (LLMs) because of their proximity to systems that are equivalent to LLMs in terms of input/output function; yet, for these functionally equivalent systems, there cannot be any non-trivial theory of consciousness that judges them conscious. This forms the basis of a disproof of contemporary LLM consciousness. I then show a positive result, which is that theories of consciousness based on (or requiring) continual learning do satisfy the stringent formal constraints for a theory of consciousness in humans. Intriguingly, this work supports a hypothesis: If continual learning is linked to consciousness in humans, the current limitations of LLMs (which do not continually learn) are intimately tied to their lack of consciousness.
Cycles Communities from the Perspective of Dendrograms and Gradient Sampling
oai:arXiv.org:2512.12974v1
arXiv:2512.12974v1 Announce Type: new
Abstract: Identifying and comparing topological features, particularly cycles, across different topological objects remains a fundamental challenge in persistent homology and topological data analysis. This work introduces a novel framework for constructing cycle communities through two complementary approaches. First, a dendrogram-based methodology leverages merge-tree algorithms to construct hierarchical representations of homology classes from persistence intervals. The Wasserstein distance on merge trees is introduced as a metric for comparing dendrograms, establishing connections to hierarchical clustering frameworks. Through simulation studies, the discriminative power of dendrogram representations for identifying cycle communities is demonstrated. Second, an extension of Stratified Gradient Sampling simultaneously learns multiple filter functions that yield cycle barycenter functions capable of faithfully reconstructing distinct sets of cycles. The set of cycles each filter function can reconstruct constitutes cycle communities that are non-overlapping and partition the space of all cycles. Together, these approaches transform the problem of cycle matching into both a hierarchical clustering and topological optimization framework, providing principled methods to identify similar topological structures both within and across groups of topological objects.
Macular: a multi-scale simulation platform for the retina and the primary visual system
oai:arXiv.org:2512.13052v1
arXiv:2512.13052v1 Announce Type: new
Abstract: We developed Macular, a simulation platform with a graphical interface, designed to produce in silico experiment scenarios for the retina and the primary visual system. A scenario consists of generating a three-dimensional structure with interconnected layers, each layer corresponding to a type of 'cell' in the retina or visual cortex. The cells can correspond to neurons or more complex structures (such as cortical columns). The inputs are arbitrary videos. The user can use the cells and synapses provided with the software, or create their own using a graphical interface where they enter the constituent equations in text format (e.g., LaTeX). They also create the three-dimensional structure via the graphical interface. Macular then automatically generates and compiles the C++ code and generates the simulation interface. This allows the user to view the input video and the three-dimensional structure in layers. It also allows the user to select cells and synapses in each layer and view the activity of their state variables. Finally, the user can adjust the phenomenological parameters of the cells or synapses via the interface. We provide several example scenarios, corresponding to published articles, including an example of a retino-cortical model. Macular was designed for neurobiologists and modelers, specialists in the primary visual system, who want to test hypotheses in silico without the need for programming. By design, this tool allows natural or altered conditions (pharmacology, pathology, development) to be simulated.
Ecological interactions and spatial dynamics in microbial aggregates: A novel modelling framework
oai:arXiv.org:2512.13156v1
arXiv:2512.13156v1 Announce Type: new
Abstract: We present a mathematical model based on a system of partial differential equations (PDEs) with cross-diffusion and reaction terms to describe ecological interactions between multiple bacterial species and substrates within microaggregates, where bacteria proliferate in response to substrate availability and undergo passive dispersal driven by population pressure gradients. The ecological interactions include interspecific competition for shared substrates, and commensalism, whereby one species benefits from the metabolic by-products of another. The main motivation comes from individual-based models (IBMs) of microbial aggregates, where simulations reveal that substrate-limited conditions can give rise to rich spatial patterns. Our numerical experiments demonstrate that our PDE-based model captures the key qualitative features of three verification scenarios that have previously been investigated with IBMs. Moreover, we formally derive a competition system from an on-lattice biased random walk, and establish local well-posedness for a parameter-symmetric subcase of it. We then formally analyse the travelling wave behaviour of this case in one spatial dimension and compare the minimal travelling wave speed with the wave speed measured in the simulations.
Stable equilibria in the Lotka-Volterra equations
oai:arXiv.org:2512.13347v1
arXiv:2512.13347v1 Announce Type: new
Abstract: We consider the Lotka-Volterra system and provide necessary conditions for an equilibrium to be stable. Our results naturally complement earlier fundamental results by N. Adachi, Y. Takeuchi, and H. Tokumaru, who, in a series of papers, give sufficient (and for some cases necessary) conditions for the existence of a stable equilibrium point.
Nondimensionalization is more science than art
oai:arXiv.org:2512.13455v1
arXiv:2512.13455v1 Announce Type: new
Abstract: When faced with a mathematical model, often the first step is to reduce the complexity of the model by turning variables and parameters into dimensionless quantities. This process is often performed by hand, relying on a skill practiced over many years, and attempted for small models. Nondimensionalization is often considered an art, as there is no formal method accessible to applied scientists. Here we show how to systematically perform nondimensionalization for arbitrarily sized models described by rational first order ordinary differential equations. We translate and extend an existing approach for computing rational invariants of the maximal scaling symmetry, which combines ideas from differential algebra, invariant theory and linear algebra, to the setting arising in biological models. The modeler inputs the system of equations and our implemented algorithm outputs the nondimensional quantities for the corresponding nondimensionalized model. We extend the algorithm to include initial conditions, and the modeler's choice of invariants, thereby including a larger class of nondimensionalizations. We further prove that any dimensionally consistent change of variables preserves the dimension of the maximal scaling symmetry. We showcase the framework on various models, including the classical Michaelis-Menten equations, which serves as a benchmark for asking and answering specific modeling questions.
A Deep Learning Model of Mental Rotation Informed by Interactive VR Experiments
oai:arXiv.org:2512.13517v1
arXiv:2512.13517v1 Announce Type: new
Abstract: Mental rotation -- the ability to compare objects seen from different viewpoints -- is a fundamental example of mental simulation and spatial world modelling in humans. Here we propose a mechanistic model of human mental rotation, leveraging advances in deep, equivariant, and neuro-symbolic learning. Our model consists of three stacked components: (1) an equivariant neural encoder, taking images as input and producing 3D spatial representations of objects, (2) a neuro-symbolic object encoder, deriving symbolic descriptions of objects from these spatial representations, and (3) a neural decision agent, comparing these symbolic descriptions to prescribe rotation simulations in 3D latent space via a recurrent pathway. Our model design is guided by the abundant experimental literature on mental rotation, which we complemented with experiments in VR where participants could at times manipulate the objects to compare, providing us with additional insights into the cognitive process of mental rotation. Our model captures well the performance, response times and behavior of participants in our and others' experiments. The necessity of each model component is shown through systematic ablations. Our work adds to a recent collection of deep neural models of human spatial reasoning, further demonstrating the potency of integrating deep, equivariant, and symbolic representations to model the human mind.
Altered oscillatory brain networks during emotional face processing in ADHD: an eLORETA and functional ICA study
oai:arXiv.org:2512.13539v1
arXiv:2512.13539v1 Announce Type: new
Abstract: Attention-deficit/hyperactivity disorder (ADHD) is characterized by executive dysfunction and difficulties in processing emotional facial expressions, yet the large-scale neural dynamics underlying these impairments remain insufficiently understood. This study applied network-based EEG source analysis to examine oscillatory cortical activity during cognitive and emotional Go/NoGo tasks in individuals with ADHD. EEG data from 272 participants (ADHD n equals 102, controls n equals 170, age range 6 to 60 years) were analyzed using exact low-resolution brain electromagnetic tomography combined with functional independent component analysis, yielding ten frequency-resolved cortical networks. Mixed-effects ANCOVAs were conducted on independent component loadings with Group, Task, and Condition as factors and age and sex as covariates. ADHD participants showed statistically significant but small increases in activation across several networks, including a gamma-dominant inferior temporal component showing a Group effect and a Group by Condition interaction with stronger NoGo-related activation in ADHD. Two additional components showed similar but weaker NoGo-selective patterns. A main effect of Task emerged only for one temporal delta component, with higher activation during the VCPT than the ECPT. No Group by Task interactions were observed. Behavioral results replicated the established ADHD performance profile, with slower responses, greater variability, and higher error rates, particularly during the emotional ECPT. Overall, the findings reveal subtle alterations in oscillatory brain networks during inhibitory processing in ADHD, with modest effect sizes embedded within substantial within-group variability. These results support a dimensional view of ADHD neurobiology and highlight the limited discriminative power of network-level EEG markers.
BaRISTA: Brain Scale Informed Spatiotemporal Representation of Human Intracranial Neural Activity
oai:arXiv.org:2512.12135v1
arXiv:2512.12135v1 Announce Type: cross
Abstract: Intracranial recordings have opened a unique opportunity to simultaneously measure activity across multiregional networks in the human brain. Recent works have focused on developing transformer-based neurofoundation models of such recordings that can generalize across subjects and datasets. However, these recordings exhibit highly complex spatiotemporal interactions across diverse spatial scales, from the single-channel scale to the scale of brain regions. As such, there remain critical open questions regarding how best to encode spatial information and how to design self-supervision tasks that enable the learning of brain network patterns and enhance downstream decoding performance using such high-dimensional, multiregional recordings. To allow for exploring these questions, we propose a new spatiotemporal transformer model of multiregional neural activity and a corresponding self-supervised masked latent reconstruction task, designed to enable flexibility in the spatial scale used for token encoding and masking. Applying this model on publicly available multiregional intracranial electrophysiology (iEEG) data, we demonstrate that adjusting the spatial scale for both token encoding and masked reconstruction significantly impacts downstream decoding. Further, we find that spatial encoding at larger scales than channel-level encoding, which is commonly used in existing iEEG transformer models, improves downstream decoding performance. Finally, we demonstrate that our method allows for region-level token encoding while also maintaining accurate channel-level neural reconstruction. Taken together, our modeling framework enables exploration of the spatial scales used for token encoding and masking, reveals their importance towards self-supervised pretraining of neurofoundation models of multiregional human brain activity, and enhances downstream decoding performance.
MolGuidance: Advanced Guidance Strategies for Conditional Molecular Generation with Flow Matching
oai:arXiv.org:2512.12198v1
arXiv:2512.12198v1 Announce Type: cross
Abstract: Key objectives in conditional molecular generation include ensuring chemical validity, aligning generated molecules with target properties, promoting structural diversity, and enabling efficient sampling for discovery. Recent advances in computer vision introduced a range of new guidance strategies for generative models, many of which can be adapted to support these goals. In this work, we integrate state-of-the-art guidance methods -- including classifier-free guidance, autoguidance, and model guidance -- in a leading molecule generation framework built on an SE(3)-equivariant flow matching process. We propose a hybrid guidance strategy that separately guides continuous and discrete molecular modalities -- operating on velocity fields and predicted logits, respectively -- while jointly optimizing their guidance scales via Bayesian optimization. Our implementation, benchmarked on the QM9 and QMe14S datasets, achieves new state-of-the-art performance in property alignment for de novo molecular generation. The generated molecules also exhibit high structural validity. Furthermore, we systematically compare the strengths and limitations of various guidance methods, offering insights into their broader applicability.
Morphogenesis of bacterial colonies in liquid crystalline environments
oai:arXiv.org:2512.12406v1
arXiv:2512.12406v1 Announce Type: cross
Abstract: Natural bacterial habitats are often complex fluids with viscoelastic and anisotropic responses to stress; for example, they can take the form of liquid crystals (LCs), with elongated microscopic constituents that collectively align while still retaining the ability to flow. However, laboratory studies typically focus on cells in simple liquids or complex fluids with randomly-oriented constituents. Here, we show how interactions with LCs shape bacterial proliferation in multicellular colonies. Using experiments, we find that in a nematic LC, cells generically form aligned single-cell-wide "chains" as they reproduce. As these chains lengthen, they eventually buckle in a highly localized manner. By combining our measurements with a continuum mechanical theory, we demonstrate that this distinctive morphogenetic program emerges because cells are kept in alignment due to the LC's elasticity; as each chain lengthens, growth-induced viscous stresses along its contour eventually overcome the elasticity of the surrounding nematic, leading to buckling. Our work thus reveals and provides mechanistic insight into the previously-overlooked role of LCs in sculpting bacterial life in complex environments.
Cross-Modal Representational Knowledge Distillation for Enhanced Spike-Informed LFP Modeling
oai:arXiv.org:2512.12461v1
arXiv:2512.12461v1 Announce Type: cross
Abstract: Local field potentials (LFPs) can be routinely recorded alongside spiking activity in intracortical neural experiments, measure a larger complementary spatiotemporal scale of brain activity for scientific inquiry, and can offer practical advantages over spikes, including greater long-term stability, robustness to electrode degradation, and lower power requirements. Despite these advantages, recent neural modeling frameworks have largely focused on spiking activity since LFP signals pose inherent modeling challenges due to their aggregate, population-level nature, often leading to lower predictive power for downstream task variables such as motor behavior. To address this challenge, we introduce a cross-modal knowledge distillation framework that transfers high-fidelity representational knowledge from pretrained multi-session spike transformer models to LFP transformer models. Specifically, we first train a teacher spike model across multiple recording sessions using a masked autoencoding objective with a session-specific neural tokenization strategy. We then align the latent representations of the student LFP model to those of the teacher spike model. Our results show that the Distilled LFP models consistently outperform single- and multi-session LFP baselines in both fully unsupervised and supervised settings, and can generalize to other sessions without additional distillation while maintaining superior performance. These findings demonstrate that cross-modal knowledge distillation is a powerful and scalable approach for leveraging high-performing spike models to develop more accurate LFP models.
Dynamical modeling of nonlinear latent factors in multiscale neural activity with real-time inference
oai:arXiv.org:2512.12462v1
arXiv:2512.12462v1 Announce Type: cross
Abstract: Real-time decoding of target variables from multiple simultaneously recorded neural time-series modalities, such as discrete spiking activity and continuous field potentials, is important across various neuroscience applications. However, a major challenge for doing so is that different neural modalities can have different timescales (i.e., sampling rates) and different probabilistic distributions, or can even be missing at some time-steps. Existing nonlinear models of multimodal neural activity do not address different timescales or missing samples across modalities. Further, some of these models do not allow for real-time decoding. Here, we develop a learning framework that can enable real-time recursive decoding while nonlinearly aggregating information across multiple modalities with different timescales and distributions and with missing samples. This framework consists of 1) a multiscale encoder that nonlinearly aggregates information after learning within-modality dynamics to handle different timescales and missing samples in real time, 2) a multiscale dynamical backbone that extracts multimodal temporal dynamics and enables real-time recursive decoding, and 3) modality-specific decoders to account for different probabilistic distributions across modalities. In both simulations and three distinct multiscale brain datasets, we show that our model can aggregate information across modalities with different timescales and distributions and missing samples to improve real-time target decoding. Further, our method outperforms various linear and nonlinear multimodal benchmarks in doing so.
Pattern Formation Beyond Turing: Physical Principles of Mass-Conserving Reaction--Diffusion Systems
oai:arXiv.org:2512.12558v1
arXiv:2512.12558v1 Announce Type: cross
Abstract: Intracellular protein patterns govern essential cellular functions by dynamically redistributing proteins between membrane-bound and cytosolic states, conserving their total numbers. This review presents a theoretical framework for understanding such patterns based on mass-conserving reaction--diffusion systems. The emergence, selection, and evolution of patterns are analyzed in terms of mass redistribution and interface motion, resulting in mesoscale laws of coarsening and wavelength selection. A geometric phase-space perspective provides a conceptual tool to link local reactive equilibria with global pattern dynamics through conserved mass fluxes. The Min protein system of \emph{Escherichia coli} provides a paradigmatic example, enabling direct comparison between theory and experiment. Successive model refinements capture both the robustness of pattern formation and the diversity of dynamic regimes observed \emph{in vivo} and \emph{in vitro}. The Min system thus illustrates how to extract predictive, multiscale theory from biochemical detail, providing a foundation for understanding pattern formation in more complex and synthetic systems.
Unsupervised learning of multiscale switching dynamical system models from multimodal neural data
oai:arXiv.org:2512.12881v1
arXiv:2512.12881v1 Announce Type: cross
Abstract: Neural population activity often exhibits regime-dependent non-stationarity in the form of switching dynamics. Learning accurate switching dynamical system models can reveal how behavior is encoded in neural activity. Existing switching approaches have primarily focused on learning models from a single neural modality, either continuous Gaussian signals or discrete Poisson signals. However, multiple neural modalities are often recorded simultaneously to measure different spatiotemporal scales of brain activity, and all these modalities can encode behavior. Moreover, regime labels are typically unavailable in training data, posing a significant challenge for learning models of regime-dependent switching dynamics. To address these challenges, we develop a novel unsupervised learning algorithm that learns the parameters of switching multiscale dynamical system models using only multiscale neural observations. We demonstrate our method using both simulations and two distinct experimental datasets with multimodal spike-LFP observations during different motor tasks. We find that our switching multiscale dynamical system models more accurately decode behavior than switching single-scale dynamical models, showing the success of multiscale neural fusion. Further, our models outperform stationary multiscale models, illustrating the importance of tracking regime-dependent non-stationarity in multimodal neural data. The developed unsupervised learning framework enables more accurate modeling of complex multiscale neural dynamics by leveraging information in multimodal recordings while incorporating regime switches. This approach holds promise for improving the performance and robustness of brain-computer interfaces over time and for advancing our understanding of the neural basis of behavior.
Binary normal networks without near reticulations can be reconstructed from their rooted triples
oai:arXiv.org:2512.12969v1
arXiv:2512.12969v1 Announce Type: cross
Abstract: Normal networks are an important class of phylogenetic networks that have compelling mathematical properties which align with intuition about inference from genetic data. While tools enabling widespread use of phylogenetic networks in the biological literature are still under mathematical, statistical, and computational development, many such results are being assembled, and in particular for normal phylogenetic networks. For instance, it has been shown that binary normal networks can be reconstructed from the sets of three- and four-leaf rooted phylogenetic trees that they display. It is also known that one can reconstruct particular subclasses of normal networks from just the displayed rooted triples. This applies, for instance, to rooted binary phylogenetic trees and to binary level-$1$ normal networks. In this paper we address the question of how much of the class of binary normal networks can be reconstructed from just the rooted triples that they display. We find that all except those with substructures that we call ``near-sibling reticulations'' and ``near-stack reticulations'' can be reconstructed just from their rooted triples. This goes some way to answering the natural question of how much information can be extracted from a set of displayed rooted triples, which are arguably the simplest substructure that one may hope for in a phylogenetic object.
Deep Learning-Driven Inversion Framework for Shear Modulus Estimation in Magnetic Resonance Elastography (DIME)
oai:arXiv.org:2512.13010v1
arXiv:2512.13010v1 Announce Type: cross
Abstract: The Multimodal Direct Inversion (MMDI) algorithm is widely used in Magnetic Resonance Elastography (MRE) to estimate tissue shear stiffness. However, MMDI relies on the Helmholtz equation, which assumes wave propagation in a uniform, homogeneous, and infinite medium. Furthermore, the use of the Laplacian operator makes MMDI highly sensitive to noise, which compromises the accuracy and reliability of stiffness estimates. In this study, we propose the Deep-Learning driven Inversion Framework for Shear Modulus Estimation in MRE (DIME), aimed at enhancing the robustness of inversion. DIME is trained on the displacement fields-stiffness maps pair generated through Finite Element Modelling (FEM) simulations. To capture local wave behavior and improve robustness to global image variations, DIME is trained on small image patches. We first validated DIME using homogeneous and heterogeneous datasets simulated with FEM, where DIME produced stiffness maps with low inter-pixel variability, accurate boundary delineation, and higher correlation with ground truth (GT) compared to MMDI. Next, DIME was evaluated in a realistic anatomy-informed simulated liver dataset with known GT and compared directly to MMDI. DIME reproduced ground-truth stiffness patterns with high fidelity (r = 0.99, R^2 = 0.98), while MMDI showed greater underestimation. After validating DIME on synthetic data, we tested the model in in vivo liver MRE data from eight healthy and seven fibrotic subjects. DIME preserved physiologically consistent stiffness patterns and closely matched MMDI, which showed directional bias. Overall, DIME showed higher correlation with ground truth and visually similar stiffness patterns, whereas MMDI displayed a larger bias that can potentially be attributed to directional filtering. These preliminary results highlight the feasibility of DIME for clinical applications in MRE.
FlowClass.jl: Classifying Dynamical Systems by Structural Properties in Julia
oai:arXiv.org:2512.13084v1
arXiv:2512.13084v1 Announce Type: cross
Abstract: FlowClass.jl is a Julia package for classifying continuous-time dynamical systems into a hierarchy of structural classes: Gradient, Gradient-like, Morse-Smale, Structurally Stable, and General. Given a vector field \(\mathbf{F}(\mathbf{x})\) defining the system \(\mathrm{d}\mathbf{x}/\mathrm{d}t = \mathbf{F}(\mathbf{x})\), the package performs a battery of computational tests -- Jacobian symmetry analysis, curl magnitude estimation, fixed point detection and stability classification, periodic orbit detection, and stable/unstable manifold computation -- to determine where the system sits within the classification hierarchy. This classification has direct implications for qualitative behaviour: gradient systems cannot oscillate, Morse-Smale systems are structurally stable in less than 3 dimensions, and general systems may exhibit chaos. Much of classical developmental theory going back to Waddington's epigenetic landscape rests on an implicit assumption of gradient dynamics.
The package is designed with applications in systems and developmental biology in mind, particularly the analysis of gene regulatory networks and cell fate decision models in the context of Waddington's epigenetic landscape. It provides tools to assess whether a landscape metaphor is appropriate for a given dynamical model, and to quantify the magnitude of non-gradient (curl) dynamics.
Vertex Model Mechanics Explain the Emergence of Centroidal Voronoi Tiling in Epithelia
oai:arXiv.org:2512.13116v1
arXiv:2512.13116v1 Announce Type: cross
Abstract: Epithelia are confluent cell layers that self-organize into polygonal networks whose geometry encodes their mechanical state. A principal driver is the tunable contractility of the actomyosin cortex, which links cell-junction tension to tissue architecture. Notably, epithelial tilings frequently resemble centroidal Voronoi tessellations (CVTs), yet the physical origin of this resemblance has remained unclear. Here, using a minimal vertex model that relates cell shape to a mechanical energy, we show that CVT-like patterns arise naturally in the solid (rigid) regime of tissues. Analytical theory reveals that isotropic strain minimization drives cell centroids toward Voronoi configurations, a result we corroborate with a analytical mean-field formulation of the vertex model. We further demonstrate that physiologically relevant perturbations -- such as cyclic stretch -- shift tissues into distinct, geometrically disordered CVT states, and that these shifts provide quantitative, image-based readouts of mechanical state. Together, our results identify a mechanical origin for CVT-like organization in epithelia and establish a geometric framework that infers tissue stresses directly from morphology, offering broadly applicable metrics for assessing rigidity and remodeling in living tissues.
Large language models are not about language
oai:arXiv.org:2512.13441v1
arXiv:2512.13441v1 Announce Type: cross
Abstract: Large Language Models are useless for linguistics, as they are probabilistic models that require a vast amount of data to analyse externalized strings of words. In contrast, human language is underpinned by a mind-internal computational system that recursively generates hierarchical thought structures. The language system grows with minimal external input and can readily distinguish between real language and impossible languages.
Understanding Cellular Noise with Optical Perturbation and Deep Learning
oai:arXiv.org:2401.12498v2
arXiv:2401.12498v2 Announce Type: replace
Abstract: Noise plays a crucial role in the regulation of cellular and organismal function and behavior.
Exploring noise's impact is key to understanding fundamental biological processes, such as gene expression, signal transduction, and the mechanisms of development and evolution.
Currently, a comprehensive method to quantify dynamical behavior of cellular noise within these biochemical systems is lacking.
In this study, we introduce an optically-controlled perturbation system utilizing the light-sensitive Phytochrome B (PhyB) from \textit{Arabidopsis thaliana}, which enables precise noise modulation with high spatial-temporal resolution.
Our system exhibits exceptional sensitivity to light, reacting consistently to pulsed light signals, distinguishing it from other photoreceptor-based promoter systems that respond to a single light wavelength.
To characterize our system, we developed a stochastic model for phytochromes that accounts for photoactivation/deactivation, thermal reversion, and the dynamics of the light-activated gene promoter system.
To precisely control our system, we determined the rate constants for this model using an omniscient deep neural network that can directly map rate constant combinations to time-dependent state joint distributions.
By adjusting the activation rates through light intensity and degradation rates via N-terminal mutagenesis, we illustrate that out optical-controlled perturbation can effectively modulate molecular expression level as well as noise.
Our results highlight the potential of employing an optically-controlled gene perturbation system as a noise-controlled stimulus source.
This approach, when combined with the analytical capabilities of a sophisticated deep neural network, enables the accurate estimation of rate constants from observational data in a broad range of biochemical reaction networks.
Dy-mer: An Explainable DNA Sequence Representation Scheme using Dictionary Learning
oai:arXiv.org:2407.12051v2
arXiv:2407.12051v2 Announce Type: replace
Abstract: DNA sequences encode critical genetic information, yet their variable length and discrete nature impede direct utilization in deep learning models. Existing DNA representation schemes convert sequences into numerical vectors but fail to capture structural features of local subsequences and often suffer from limited interpretability and poor generalization on small datasets. To address these limitations, we propose Dy-mer, an interpretable and robust DNA representation scheme based on dictionary learning. Dy-mer formulates an optimization problem in tensor format, which ensures computational efficiency in batch processing. Our scheme reconstructs DNA sequences as concatenations of dynamic-length subsequences (dymers) through a convolution operation and simultaneously optimize a learnable dymer dictionary and sparse representations. Our method achieves state-of-the-art performance in downstream tasks such as DNA promoter classification and motif detection. Experiments further show that the learned dymers match known DNA motifs and clustering using Dy-mer yields semantically meaningful phylogenetic trees. These results demonstrate that the proposed approach achieves both strong predictive performance and high interpretability, making it well suited for biological research applications.
Learning in Focus: Detecting Behavioral and Collaborative Engagement Using Vision Transformers
oai:arXiv.org:2508.15782v2
arXiv:2508.15782v2 Announce Type: replace
Abstract: In early childhood education, accurately detecting collaborative and behavioral engagement is essential to foster meaningful learning experiences. This paper presents an AI driven approach that leverages Vision Transformers (ViTs) to automatically classify children s engagement using visual cues such as gaze direction, interaction, and peer collaboration. Utilizing the ChildPlay gaze dataset, our method is trained on annotated video segments to classify behavioral and collaborative engagement states (e.g., engaged, not engaged, collaborative, not collaborative). We evaluated six state of the art transformer models: Vision Transformer (ViT), Data efficient Image Transformer (DeiT), Swin Transformer, VitGaze, APVit and GazeTR. Among these, the Swin Transformer achieved the highest classification performance with an accuracy of 97.58 percent, demonstrating its effectiveness in modeling local and global attention. Our results highlight the potential of transformer based architectures for scalable, automated engagement analysis in real world educational settings.
A self-organized compression network arrests epithelial proliferation
oai:arXiv.org:2509.16661v3
arXiv:2509.16661v3 Announce Type: replace
Abstract: As epithelial development or wound closure approaches completion, cell proliferation progressively slows via contact inhibition of proliferation (CIP) - a mechanism understood as being strictly local. Here we report the discovery of inhibition of proliferation through an unanticipated mechanism that is non-local. As a confluent epithelial layer becomes progressively more jammed, two interpenetrating networks emerge: islands of mechanically compressed non-cycling cells percolating within an ocean of mechanically tensed cycling cells. The evolution of the compression network was found to be susceptible to both specific molecular stimulus and to injury-induced unjamming. Yet, in all circumstances, the size of compressed islands followed a power-law distribution that was well-captured by preferential network theory. Together, these findings demonstrate the existence of a network-based inhibition of proliferation (NIP) that is self-organizing and poised in proximity to criticality.
Human-computer interactions predict mental health
oai:arXiv.org:2511.20179v2
arXiv:2511.20179v2 Announce Type: replace
Abstract: Scalable assessments of mental illness, the leading driver of disability worldwide, remain a critical roadblock toward accessible and equitable care. Here, we show that human-computer interactions encode mental health with state-of-the-art biomarker precision. We introduce MAILA, a MAchine-learning framework for Inferring Latent mental states from digital Activity. We trained MAILA to predict 1.3 million mental-health self-reports from 20,000 cursor and touchscreen recordings recorded in 9,000 online participants. The dataset includes 2,000 individuals assessed longitudinally, 1,500 diagnosed with depression, and 500 with obsessive-compulsive disorder. MAILA tracks dynamic mental states along three orthogonal dimensions, identifies individuals living with mental illness, and achieves near-ceiling accuracy when predicting group-level mental health. By extracting non-verbal signatures of psychological function that have so far remained untapped, MAILA represents a key step toward foundation models for mental health. The ability to decode mental states at zero marginal cost creates new opportunities in neuroscience, medicine, and public health, while raising urgent questions about privacy, agency, and autonomy online.
Automated Histopathologic Assessment of Hirschsprung Disease Using a Multi-Stage Vision Transformer Framework
oai:arXiv.org:2511.20734v2
arXiv:2511.20734v2 Announce Type: replace
Abstract: Hirschsprung Disease is characterized by the absence of ganglion cells in the myenteric plexus. Therefore, the correct identification of ganglion cells is crucial for diagnosing Hirschsprung disease. We introduce a three-stage analysis framework that mimics the pathologist's diagnostic approach. The framework, based on a Vision Transformer model (ViT-B/16), sequentially segments the muscularis propria, segments the myenteric plexus, and detects ganglion cells within anatomically valid regions. 30 whole-slide images of colon tissue were used, each containing manual annotations of muscularis, plexus, and ganglion cells. A 5-fold cross-validation scheme was applied to each stage, along with resolution-specific tiling strategies and tailored postprocessing to ensure anatomical consistency. The proposed method achieved a Dice coefficient of 89.9% and a Plexus Inclusion Rate of 100% for muscularis segmentation. Plexus segmentation reached a recall of 94.8%, a precision of 84.2% and a Ganglia Inclusion Rate of 99.7%. For ganglion cells annotated with high certainty, the model achieved 62.1\% precision and 89.1% recall. When considering all annotated ganglion cells, regardless of certainty level, the overall precision was 67.0%. These results indicate that ViT-based models are effective at leveraging global tissue context and capturing cellular morphology at small scales, even within complex histological tissue structures. This multi-stage methodology has great potential to support digital pathology workflows by reducing inter-observer variability and assisting in the evaluation of Hirschsprung disease. The clinical impact will be evaluated in future work with larger multi-center datasets and additional expert annotations.
Stabilizing Fractional Dynamical Networks Suppresses Epileptic Seizures
oai:arXiv.org:2511.20950v2
arXiv:2511.20950v2 Announce Type: replace
Abstract: Medically uncontrolled epileptic seizures affect nearly 15 million people worldwide, resulting in enormous economic and psychological burdens. Treatment of medically refractory epilepsy is essential for patients to achieve remission, improve psychological functioning, and enhance social and vocational outcomes. Here, we show a state-of-the-art method that stabilizes fractional dynamical networks modeled from intracranial EEG data, effectively suppressing seizure activity in 34 out of 35 total spontaneous episodes from patients at the University of Pennsylvania and the Mayo Clinic. We perform a multi-scale analysis and show that the fractal behavior and stability properties of these data distinguish between four epileptic states: interictal, pre-ictal, ictal, and post-ictal. Furthermore, the simulated controlled signals exhibit substantial amplitude reduction ($49\%$ average). These findings highlight the potential of fractional dynamics to characterize seizure-related brain states and demonstrate its capability to suppress epileptic activity.
Hierarchical Molecular Language Models (HMLMs)
oai:arXiv.org:2512.00696v3
arXiv:2512.00696v3 Announce Type: replace
Abstract: Artificial intelligence (AI) is reshaping computational and network biology by enabling new approaches to decode cellular communication networks. We introduce Hierarchical Molecular Language Models (HMLMs), a novel framework that models cellular signaling as a specialized molecular language, where signaling molecules function as tokens, protein interactions define syntax, and functional consequences constitute semantics. HMLMs employ a transformer-based architecture adapted to accommodate graph-structured signaling networks through information transducers, mathematical entities that capture how molecules receive, process, and transmit signals. The architecture integrates multi-modal data sources across molecular, pathway, and cellular scales through hierarchical attention mechanisms and scale-bridging operators that enable information flow across biological hierarchies. Applied to a complex network of cardiac fibroblast signaling, HMLMs outperformed traditional approaches in temporal dynamics prediction, particularly under sparse sampling conditions. Attention-based analysis revealed biologically meaningful crosstalk patterns, including previously uncharacterized interactions between signaling pathways. By bridging molecular mechanisms with cellular phenotypes through AI-driven molecular language representation, HMLMs establish a foundation for biology-oriented large language models (LLMs) that could be pre-trained on comprehensive pathway datasets and applied across diverse signaling systems and tissues, advancing precision medicine and therapeutic discovery.
Tracking large chemical reaction networks and rare events by neural networks
oai:arXiv.org:2512.10309v2
arXiv:2512.10309v2 Announce Type: replace
Abstract: Chemical reaction networks are widely used to model stochastic dynamics in chemical kinetics, systems biology and epidemiology. Solving the chemical master equation that governs these systems poses a significant challenge due to the large state space exponentially growing with system sizes. The development of autoregressive neural networks offers a flexible framework for this problem; however, its efficiency is limited especially for high-dimensional systems and in scenarios with rare events. Here, we push the frontier of neural-network approach by exploiting faster optimizations such as natural gradient descent and time-dependent variational principle, achieving a 5- to 22-fold speedup, and by leveraging enhanced-sampling strategies to capture rare events. We demonstrate reduced computational cost and higher accuracy over the previous neural-network method in challenging reaction networks, including the mitogen-activated protein kinase (MAPK) cascade network, the hitherto largest biological network handled by the previous approaches of solving the chemical master equation. We further apply the approach to spatially extended reaction-diffusion systems, the Schl\"ogl model with rare events, on two-dimensional lattices, beyond the recent tensor-network approach that handles one-dimensional lattices. The present approach thus enables efficient modeling of chemical reaction networks in general.
Foveated Retinotopy Improves Classification and Localization in Convolutional Neural Networks
oai:arXiv.org:2402.15480v5
arXiv:2402.15480v5 Announce Type: replace-cross
Abstract: From falcons spotting preys to humans recognizing faces, rapid visual abilities depend on a foveated retinal organization which delivers high-acuity central vision while preserving low-resolution periphery. This organization is conserved along early visual pathways but remains underexplored in machine learning. Here we examine how embedding a foveated retinotopic transformation as a preprocessing layer impacts convolutional neural networks (CNNs) for image classification. By applying a log-polar mapping to off-the-shelf models and retraining them, we retain comparable accuracy while improving robustness to scale and rotation. We show that this architecture becomes highly sensitive to fixation-point shifts, and that this sensitivity yields a proxy for defining saliency maps that effectively facilitates object localization. Our results show that foveated retinotopy encodes prior geometric knowledge, offering a solution to visual-search and enhancing both classification and localization. These findings connect biological vision principles with artificial networks, pointing to new, robust and efficient directions for computer-vision systems.
Wavespeed selection and interstitial gap formation in an acid-mediated cancer invasion model
oai:arXiv.org:2411.12232v2
arXiv:2411.12232v2 Announce Type: replace-cross
Abstract: We consider a two-component reaction-diffusion system that has previously been developed to model invasion of cells into a resident cell population. The system is an idealised version of models of tumour growth in which tumour cells degrade the surrounding tissue by increasing the acidity of the local environment. By numerically computing families of travelling wave solutions to this problem, we observe that a general initial condition with either compact support, or sufficiently large exponential decay in the far field, tends to the travelling wave solution that has the largest possible decay at its front. Initial conditions with sufficiently slow exponential decay tend to those travelling wave solutions that have the same exponential decay as their initial conditions. We also show that in the limit that the (nondimensional) degradation rate of resident cells is large, the system has similar asymptotic structure as previously observed in perturbed Fisher--KPP models. The asymptotic analysis in this limit explains the formation of an interstitial gap (a region between the invading and receding fronts, in which both cell populations are small), the width of which is logarithmically large in the limit of large degradation rate. These results show that the general mechanism behind the formation of the interstitial gap in reaction-diffusion tumour models is connected to perturbations of the Fisher-KPP system. Biologically, this implies that order of magnitude difference in degradation rate is required to produce appreciably different gap sizes, while the velocity of the invading front is largely determined by the Fisher-KPP velocity, and only very weakly affected by the presence of the interstitial gap.
PathRWKV: Enabling Whole Slide Prediction with Recurrent-Transformer
oai:arXiv.org:2503.03199v2
arXiv:2503.03199v2 Announce Type: replace-cross
Abstract: Pathological diagnosis is essential for cancer diagnosis, with whole slide image (WSI) providing histopathological and cellular information. Recent deep learning advancements have improved WSI analysis through a two-stage paradigm: tile-level feature extraction followed by slide-level modeling. In this paradigm, Transformer-based models surpass traditional multiple instance learning approaches in accuracy, yet still face four core limitations: (1) inadequate handling of variable tissue sizes across slides, (2) inability to effectively infer from all tiles for slide-level conclusions, (3) challenges in balancing model complexity with limited training data, and (4) difficulty balancing training efficiency and inference performance. Consequently, these issues limit whole-slide perception for diagnosis with restricted WSI training scales. To address them, we introduce PathRWKV, a novel state space model for slide-level feature modeling. To handle variable tissue sizes, PathRWKV employs two modules: Time Mix and Channel Mix, enabling dynamic perception of tiles for improved slide-level modeling. To draw effective conclusions, we propose an asymmetric design that samples tiles during training and iterates over all tiles at inference, scaling up to cover the entire slide. To balance model complexity and data size, we adopt linear attention and state space architecture with a Recurrent module. To balance training efficiency and inference, we design a tailored multi-task learning module handling versatile tasks simultaneously, enhancing model ability via multiple clinical indicators in slide reading. Experimental results show PathRWKV outperforms nine recent state-of-the-art methods across 10 downstream tasks on 11 datasets with 17,292 WSIs, paving its way for efficient slide-level pathological inference. The project is open-sourced.
Decoding and Engineering the Phytobiome Communication for Smart Agriculture
oai:arXiv.org:2508.03584v2
arXiv:2508.03584v2 Announce Type: replace-cross
Abstract: Smart agriculture applications, integrating technologies like the Internet of Things and machine learning/artificial intelligence (ML/AI) into agriculture, hold promise to address modern challenges of rising food demand, environmental pollution, and water scarcity. Alongside the concept of the phytobiome, which defines the area including the plant, its environment, and associated organisms, and the recent emergence of molecular communication (MC), there exists an important opportunity to advance agricultural science and practice using communication theory. In this article, we motivate to use the communication engineering perspective for developing a holistic understanding of the phytobiome communication and bridge the gap between the phytobiome communication and smart agriculture. Firstly, an overview of phytobiome communication via molecular and electrophysiological signals is presented and a multi-scale framework modeling the phytobiome as a communication network is conceptualized. Then, how this framework is used to model electrophysiological signals is demonstrated with plant experiments. Furthermore, possible smart agriculture applications, such as smart irrigation and targeted delivery of agrochemicals, through engineering the phytobiome communication are proposed. These applications merge ML/AI methods with the Internet of Bio-Nano-Things enabled by MC and pave the way towards more efficient, sustainable, and eco-friendly agricultural production. Finally, the implementation challenges, open research issues, and industrial outlook for these applications are discussed.
Active Force Dynamics in Red Blood Cells Under Non-Invasive Optical Tweezers
oai:arXiv.org:2512.01417v2
arXiv:2512.01417v2 Announce Type: replace-cross
Abstract: Red blood cells (RBCs) sustain mechanical stresses associated with microcirculatory flow through ATP-driven plasma membrane flickering. This is an active phenomenon driven by motor proteins that regulate interactions between the spectrin cytoskeleton and the lipid bilayer; it is manifested in RBC shape fluctuations reflecting the cell's mechanical and metabolic state. Yet, direct quantification of the forces and energetic costs underlying this non-equilibrium behavior remains challenging due to the invasiveness of existing techniques. Here, a minimally invasive method that combines bead-free, low-power optical tweezers with high-speed video microscopy was employed to track local membrane forces and displacements in single RBCs during the same time window. This independent dual-channel measurement enabled the construction of a mechano-dynamic phase space for RBCs under different chemical treatments, that allowed for differentiating between metabolic and structural states based on their fluctuation-force signatures. Quantification of mechanical work during flickering demonstrated that membrane softening enhanced fluctuations while elevating energy dissipation. The proposed optical tweezers methodology provides a robust framework for mapping the active mechanics of living cells, enabling precise probing of cellular physiology and detection of biomechanical dysfunction in diseases.