Julia for biologists

Reversible molecular simulation for training classical and machine-learning force fields

The next generation of force fields for molecular dynamics will be developed using a wealth of data. Training systematically with experimental data remains a challenge, however, especially for machine-learning potentials. Differentiable molecular simulation calculates gradients of observables with respect to parameters through molecular dynamics trajectories. Here, we improve this approach by explicitly calculating gradients using a reverse-time simulation with effectively constant memory cost and a computation count similar to the forward simulation. The method is applied to learn all-atom water and gas diffusion models with different functional forms and to train a machine-learning potential for diamond from scratch. Comparison to ensemble reweighting indicates that reversible simulation can provide more accurate gradients and train to match time-dependent observables.

Sainsc: A Computational Tool for Segmentation-Free Analysis of In Situ Capture Data

Spatially resolved transcriptomics (SRT) has become the method of choice for characterising the complexity of biomedical tissue samples. Until recently, scientists were restricted to SRT methods that can profile a limited set of target genes at high spatial resolution or transcriptome-wide but at a low spatial resolution. Through recent developments, there are now methods that offer both subcellular spatial resolution and full transcriptome coverage. However, utilising these new methods' high spatial resolution and gene resolution remains elusive due to several factors, including low detection efficiency and high computational costs. Here, we present Sainsc (Segmentation-free analysis of in situ capture data), which combines a cell-segmentation-free approach with efficient data processing of transcriptome-wide nanometre-resolution spatial data. Sainsc can generate cell-type maps with accurate cell-type assignment at the nanometre scale, together with corresponding maps of the assignment scores that facilitate interpretation of the local confidence of cell-type assignment. We demonstrate its utility and accuracy for different tissues and technologies. Compared to other methods, Sainsc requires lower computational resources and has scalable performance, enabling interactive data exploration. Sainsc is compatible with common data analysis frameworks and is available as open-source software in multiple programming languages.

The Evolution of Systems Biology and Systems Medicine: From Mechanistic Models to Uncertainty Quantification

Understanding interaction mechanisms within cells, tissues, and organisms is crucial for driving developments across biology and medicine. Mathematical modeling is an essential tool for simulating such biological systems. Building on experiments, mechanistic models are widely used to describe small-scale intracellular networks. The development of sequencing techniques and computational tools has recently enabled multiscale models. Combining such larger scale network modeling with mechanistic modeling provides us with an opportunity to reveal previously unknown disease mechanisms and pharmacological interventions. Here, we review systems biology models from mechanistic models to multiscale models that integrate multiple layers of cellular networks and discuss how they can be used to shed light on disease states and even wellness-related states. Additionally, we introduce several methods that increase the certainty and accuracy of model predictions. Thus, combining mechanistic models with emerging mathematical and computational techniques can provide us with increasingly powerful tools to understand disease states and inspire drug discoveries.

Enhancing medical imaging education: integrating computing technologies, digital image processing and artificial intelligence

The rapid advancement of technology has brought significant changes to various fields, including medical imaging (MI). This discussion paper explores the integration of computing technologies (e.g. Python and MATLAB), digital image processing (e.g. image enhancement, segmentation and three-dimensional reconstruction) and artificial intelligence (AI) into the undergraduate MI curriculum. By examining current educational practices, gaps and limitations that hinder the development of future-ready MI professionals are identified. A comprehensive curriculum framework is proposed, incorporating essential computational skills, advanced image processing techniques and state-of-the-art AI tools, such as large language models like ChatGPT. The proposed curriculum framework aims to improve the quality of MI education significantly and better equip students for future professional practice and challenges while enhancing diagnostic accuracy, improving workflow efficiency and preparing students for the evolving demands of the MI field.

MetaRange.jl: A Dynamic and Metabolic Species Range Model for Plant Species

Process-based models for range dynamics are urgently needed due to increasing intensity of human-induced biodiversity change. Despite a few existing models that focus on demographic processes, their use remains limited compared to the widespread application of correlative approaches. This slow adoption is largely due to the challenges in calibrating biological parameters and the high computational demands for large-scale applications. Moreover, increasing the number of simulated processes (i.e., mechanistic complexity) may further exacerbate those reasons of delay. Therefore, balancing mechanistic complexity and computational effectiveness of process-based models is a key area for improvement. A promising research direction is to expand demographically explicit metapopulation models by integrating metabolic constraints. We translated and expanded a previously developed R metapopulation model to Julia language and published it as a Julia module. The model integrates species-specific parameters such as preferred environmental conditions, biomass and dispersal ability with demographic rates (e.g., reproductive and mortality rates) derived from local temperature and biomass via the metabolic theory of ecology. We provide a simple application example for the model in which we illustrate a typical use case by predicting the future occurrence of Orchis militaris in Bavaria under different climate change scenarios. Our results show that climate change reduces habitat suitability overall, but some regions like the Franconian Forest and the Alps see increased suitability and abundance, confirming their role as refugia. Simulating metapopulation dynamics reveals that local population dynamics and dispersal are crucial for accurate predictions. For instance, increasing dispersal distance reduces overall abundance loss but also lessens population growth in refugia. This highlights the importance of measuring traits like dispersal ability to improve climate change forecasts.

Developmental hematopoietic stem cell variation explains clonal hematopoiesis later in life

Clonal hematopoiesis becomes increasingly common with age, but its cause is enigmatic because driver mutations are often absent. Serial observations infer weak selection indicating variants are acquired much earlier in life with unexplained initial growth spurts. Here we use fluctuating CpG methylation as a lineage marker to track stem cell clonal dynamics of hematopoiesis. We show, via the shared prenatal circulation of monozygotic twins, that weak selection conferred by stem cell variation created before birth can reliably yield clonal hematopoiesis later in life. Theory indicates weak selection will lead to dominance given enough time and large enough population sizes. Human hematopoiesis satisfies both these conditions. Stochastic loss of weakly selected variants is naturally prevented by the expansion of stem cell lineages during development. The dominance of stem cell clones created before birth is supported by blood fluctuating CpG methylation patterns that exhibit low correlation between unrelated individuals but are highly correlated between many elderly monozygotic twins. Therefore, clonal hematopoiesis driven by weak selection in later life appears to reflect variation created before birth.

Physiological signal analysis and open science using the Julia language and associated software

In this mini review, we propose the use of the Julia programming language and its software as a strong candidate for reproducible, efficient, and sustainable physiological signal analysis. First, we highlight available software and Julia communities that provide top-of-the-class algorithms for all aspects of physiological signal processing despite the language's relatively young age. Julia can significantly accelerate both research and software development due to its high-level interactive language and high-performance code generation. It is also particularly suited for open and reproducible science. Openness is supported and welcomed because the overwhelming majority of Julia software programs are open source and developed openly on public platforms, primarily through individual contributions. Such an environment increases the likelihood that an individual not (originally) associated with a software program would still be willing to contribute their code, further promoting code sharing and reuse. On the other hand, Julia's exceptionally strong package manager and surrounding ecosystem make it easy to create self-contained, reproducible projects that can be instantly installed and run, irrespective of processor architecture or operating system.

Generating synthetic signaling networks for in silico modeling studies

Predictive models of signaling pathways have proven to be difficult to develop. Traditional approaches to developing mechanistic models rely on collecting experimental data and fitting a single model to that data. This approach works for simple systems but has proven unreliable for complex systems such as biological signaling networks. Thus, there is a need to develop new approaches to create predictive mechanistic models of complex systems. To meet this need, we developed a method for generating artificial signaling networks that were reasonably realistic and thus could be treated as ground truth models. These synthetic models could then be used to generate synthetic data for developing and testing algorithms designed to recover the underlying network topology and associated parameters. We defined the reaction degree and reaction distance to measure the topology of reaction networks, especially to consider enzymes. To determine whether our generated signaling networks displayed meaningful behavior, we compared them with signaling networks from the BioModels Database. This comparison indicated that our generated signaling networks had high topological similarities with BioModels signaling networks with respect to the reaction degree and distance distributions. In addition, our synthetic signaling networks had similar behavioral dynamics with respect to both steady states and oscillations, suggesting that our method generated synthetic signaling networks comparable with BioModels and thus could be useful for building network evaluation tools.

Hypergraph animals

Here we introduce simple structures for the analysis of complex hypergraphs, hypergraph animals. These structures are designed to describe the local node neighborhoods of nodes in hypergraphs. We establish their relationships to lattice animals and network motifs, and we develop their combinatorial properties for sparse and uncorrelated hypergraphs. We make use of the tight link of hypergraph animals to partition numbers, which opens up a vast mathematical framework for the analysis of hypergraph animals. We then study their abundances in random hypergraphs. Two transferable insights result from this analysis: (i) it establishes the importance of high-cardinality edges in ensembles of random hypergraphs that are inspired by the classical Erdös-Renyí random graphs; and (ii) there is a close connection between degree and hyperedge cardinality in random hypergraphs that shapes animal abundances and spectra profoundly. Both findings imply that hypergraph animals can have the potential to affect information flow and processing in complex systems. Our analysis also suggests that we need to spend more effort on investigating and developing suitable conditional ensembles of random hypergraphs that can capture real-world structures and their complex dependency structures.

A next-generation dynamic programming language Julia: Its features and applications in biological science

BACKGROUND

The advent of Julia as a sophisticated and dynamic programming language in 2012 represented a significant milestone in computational programming, mathematical analysis, and statistical modeling. Having reached its stable release in version 1.9.0 on May 7, 2023, Julia has developed into a powerful and versatile instrument. Despite its potential and widespread adoption across various scientific and technical domains, there exists a noticeable knowledge gap in comprehending its utilization within biological sciences.

THE AIM OF REVIEW

This comprehensive review aims to address this particular knowledge gap and offer a thorough examination of Julia's fundamental characteristics and its applications in biology.

KEY SCIENTIFIC CONCEPTS OF THE REVIEW

The review focuses on a research gap in the biological science. The review aims to equip researchers with knowledge and tools to utilize Julia's capabilities in biological science effectively and to demonstrate the gap. It paves the way for innovative solutions and discoveries in this rapidly evolving field. It encompasses an analysis of Julia's characteristics, packages, and performance compared to the other programming languages in this field. The initial part of this review discusses the key features of Julia, such as its dynamic and interactive nature, fast processing speed, ease of expression manipulation, user-friendly syntax, code readability, strong support for multiple dispatch, and advanced type system. It also explores Julia's capabilities in data analysis, visualization, machine learning, and algorithms, making it suitable for scientific applications. The next section emphasizes the importance of using Julia in biological research, highlighting its seamless integration with biological studies for data analysis, and computational biology. It also compares Julia with other programming languages commonly used in biological research through benchmarking and performance analysis. Additionally, it provides insights into future directions and potential challenges in Julia's applications in biology.

RankCompV3: a differential expression analysis algorithm based on relative expression orderings and applications in single-cell RNA transcriptomics

BACKGROUND

Effective identification of differentially expressed genes (DEGs) has been challenging for single-cell RNA sequencing (scRNA-seq) profiles. Many existing algorithms have high false positive rates (FPRs) and often fail to identify weak biological signals.

RESULTS

We present a novel method for identifying DEGs in scRNA-seq data called RankCompV3. It is based on the comparison of relative expression orderings (REOs) of gene pairs which are determined by comparing the expression levels of a pair of genes in a set of single-cell profiles. The numbers of genes with consistently higher or lower expression levels than the gene of interest are counted in two groups in comparison, respectively, and the result is tabulated in a 3 × 3 contingency table which is tested by McCullagh's method to determine if the gene is dysregulated. In both simulated and real scRNA-seq data, RankCompV3 tightly controlled the FPR and demonstrated high accuracy, outperforming 11 other common single-cell DEG detection algorithms. Analysis with either regular single-cell or synthetic pseudo-bulk profiles produced highly concordant DEGs with the ground-truth. In addition, RankCompV3 demonstrates higher sensitivity to weak biological signals than other methods. The algorithm was implemented using Julia and can be called in R. The source code is available at https://github.com/pathint/RankCompV3.jl .

CONCLUSIONS

The REOs-based algorithm is a valuable tool for analyzing single-cell RNA profiles and identifying DEGs with high accuracy and sensitivity.

FlowAtlas: an interactive tool for high-dimensional immunophenotyping analysis bridging FlowJo with computational tools in Julia

As the dimensionality, throughput and complexity of cytometry data increases, so does the demand for user-friendly, interactive analysis tools that leverage high-performance machine learning frameworks. Here we introduce FlowAtlas: an interactive web application that enables dimensionality reduction of cytometry data without down-sampling and that is compatible with datasets stained with non-identical panels. FlowAtlas bridges the user-friendly environment of FlowJo and computational tools in Julia developed by the scientific machine learning community, eliminating the need for coding and bioinformatics expertise. New population discovery and detection of rare populations in FlowAtlas is intuitive and rapid. We demonstrate the capabilities of FlowAtlas using a human multi-tissue, multi-donor immune cell dataset, highlighting key immunological findings. FlowAtlas is available at https://github.com/gszep/FlowAtlas.jl.git.

Differentiable simulation to develop molecular dynamics force fields for disordered proteins

Implicit solvent force fields are computationally efficient but can be unsuitable for running molecular dynamics on disordered proteins. Here I improve the a99SB-disp force field and the GBNeck2 implicit solvent model to better describe disordered proteins. Differentiable molecular simulations with 5 ns trajectories are used to jointly optimise 108 parameters to better match explicit solvent trajectories. Simulations with the improved force field better reproduce the radius of gyration and secondary structure content seen in experiments, whilst showing slightly degraded performance on folded proteins and protein complexes. The force field, called GB99dms, reproduces the results of a small molecule binding study and improves agreement with experiment for the aggregation of amyloid peptides. GB99dms, which can be used in OpenMM, is available at https://github.com/greener-group/GB99dms. This work is the first to show that gradients can be obtained directly from nanosecond-length differentiable simulations of biomolecules and highlights the effectiveness of this approach to training whole force fields to match desired properties.

Developmental hematopoietic stem cell variation explains clonal hematopoiesis later in life

Clonal hematopoiesis becomes increasingly common with age, but its cause is enigmatic because driver mutations are often absent. Serial observations infer weak selection indicating variants are acquired much earlier in life with unexplained initial growth spurts. Here we use fluctuating CpG methylation as a lineage marker to track stem cell clonal dynamics of hematopoiesis. We show, via the shared prenatal circulation of monozygotic twins, that weak selection conferred by stem cell variation created before birth can reliably yield clonal hematopoiesis later in life. Theory indicates weak selection will lead to dominance given enough time and large enough population sizes. Human hematopoiesis satisfies both these conditions. Stochastic loss of weakly selected variants is naturally prevented by the expansion of stem cell lineages during development. The dominance of stem cell clones created before birth is supported by blood fluctuating CpG methylation patterns that exhibit low correlation between unrelated individuals but are highly correlated between many elderly monozygotic twins. Therefore, clonal hematopoiesis driven by weak selection in later life appears to reflect variation created before birth.

SBMLToolkit.jl: a Julia package for importing SBML into the SciML ecosystem

Julia is a general purpose programming language that was designed for simplifying and accelerating numerical analysis and computational science. In particular the Scientific Machine Learning (SciML) ecosystem of Julia packages includes frameworks for high-performance symbolic-numeric computations. It allows users to automatically enhance high-level descriptions of their models with symbolic preprocessing and automatic sparsification and parallelization of computations. This enables performant solution of differential equations, efficient parameter estimation and methodologies for automated model discovery with neural differential equations and sparse identification of nonlinear dynamics. To give the systems biology community easy access to SciML, we developed SBMLToolkit.jl. SBMLToolkit.jl imports dynamic SBML models into the SciML ecosystem to accelerate model simulation and fitting of kinetic parameters. By providing computational systems biologists with easy access to the open-source Julia ecosystevnm, we hope to catalyze the development of further Julia tools in this domain and the growth of the Julia bioscience community. SBMLToolkit.jl is freely available under the MIT license. The source code is available at https://github.com/SciML/SBMLToolkit.jl.

Mathematical Model of Intrinsic Drug Resistance in Lung Cancer

Drug resistance is a bottleneck in cancer treatment. Commonly, a molecular treatment for cancer leads to the emergence of drug resistance in the long term. Thus, some drugs, despite their initial excellent response, are withdrawn from the market. Lung cancer is one of the most mutated cancers, leading to dozens of targeted therapeutics available against it. Here, we have developed a mechanistic mathematical model describing sensitization to nine groups of targeted therapeutics and the emergence of intrinsic drug resistance. As we focus only on intrinsic drug resistance, we perform the computer simulations of the model only until clinical diagnosis. We have utilized, for model calibration, the whole-exome sequencing data combined with clinical information from over 1000 non-small-cell lung cancer patients. Next, the model has been applied to find an answer to the following questions: When does intrinsic drug resistance emerge? And how long does it take for early-stage lung cancer to grow to an advanced stage? The results show that drug resistance is inevitable at diagnosis but not always detectable and that the time interval between early and advanced-stage tumors depends on the selection advantage of cancer cells.

BondGraphs.jl: composable energy-based modelling in systems biology

SUMMARY

BondGraphs.jl is a Julia implementation of bond graphs. Bond graphs provide a modelling framework that describes energy flow through a physical system and by construction enforce thermodynamic constraints. The framework is widely used in engineering and has recently been shown to be a powerful approach for modelling biology. Models are mutable, hierarchical, multiscale, and multiphysics, and BondGraphs.jl is compatible with the Julia modelling ecosystem.

AVAILABILITY AND IMPLEMENTATION

BondGraphs.jl is freely available under the MIT license. Source code and documentation can be found at https://github.com/jedforrest/BondGraphs.jl.