Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2025 Aug 2;41(8):btaf394.
doi: 10.1093/bioinformatics/btaf394.

Dynamic gene regulatory network inference from single-cell data using optimal transport

François Lamoline  1 , Isabel Haasler  2   3 , Johan Karlsson  3 , Jorge Gonçalves  1   4 , Atte Aalto  1   5
Affiliations

Dynamic gene regulatory network inference from single-cell data using optimal transport

François Lamoline et al. Bioinformatics. .

Abstract

Motivation: Modelling gene expression is a central problem in systems biology. Single-cell technologies have revolutionized the field by enabling sequencing at the resolution of individual cells. This results in a much richer data compared to what is obtained by bulk technologies, offering new possibilities and challenges for gene regulatory network inference.

Results: In this work, we introduce GRIT (gene regulation inference by transport)-a method to fit a differential equation model and to infer gene regulatory networks from single-cell data using the theory of optimal transport. The idea consists in tracking the evolution of the cell distribution over time and finding the system whose temporal marginals minimize the transport cost with the observations. GRIT is finally used to identify genes and pathways affected by two Parkinson's disease associated mutations.

Availability and implementation: Matlab implementation of the method and code for data generation are at gitlab.com/uniluxembourg/lcsb/systems-control/grit together with a user guide. A snapshot of the code used for the results of this article is at doi: 10.5281/zenodo.15582432.

PubMed Disclaimer

Figures

Figure 1.
Figure 1.
Method pipeline: The model identification iteration is run until convergence. The iteration can be interpreted as a block coordinate descent for solving (3), alternating between minimization with respect to the transport plans Mk and the model (A,b). The results are then used in a separate variable selection step to obtain confidence scores for GRN link existence.
Figure 2.
Figure 2.
(a) The squared model error ∥[A,b]−[A¯,b¯]∥F2 with data simulated from a linear discrete-time system with varying number of timepoints and number of cells per timepoint. The plot is in logarithmic scale and it shows the mean and 80th percentiles obtained from five replicates. The sloped gridlines correspond to a decay m−1/2. (b) The AUROC and AUPR scores from 20 replicates with data simulated from linear discrete-time system (LD), linear continuous-time system (LC), nonlinear system (NL), or system with non-observed protein concentrations (Pr).
Figure 3.
Figure 3.
Results summary on the BEELINE benchmarking pipeline for (a) the synthetic dataset, (b) curated dataset, and (c) the RNA-Seq data. The bars show the average rank of each method in the different tasks, and the error bars indicate the means of the bottom 50th percentile and the top 50th percentile.
Figure 4.
Figure 4.
Histogram of perturbation target scores with high-confidence genes indicated for the LRRK2 mutation (a) and for the PINK1 mutation (b).
See this image and copyright information in PMC

References

    1. Aalto A, Gonçalves J. Linear system identification from ensemble snapshot observations. In: 2019 IEEE 58th Conference on Decision and Control (CDC), Nice, France. IEEE Control Systems Society, 2019, 7554–9. 10.1109/CDC40024.2019.9029334 - DOI
    1. Aalto A, Lamoline F, Gonçalves J. Linear system identifiability from single-cell data. Syst Control Lett 2022;165:105287. 10.1016/j.sysconle.2022.105287 - DOI
    1. Aalto A, Viitasaari L, Ilmonen P et al. Gene regulatory network inference from sparsely sampled noisy data. Nat Commun 2020;11:3493. 10.1038/s41467-020-17217-1 - DOI - PMC - PubMed
    1. Ahlmann-Eltze C, Huber W. Comparison of transformations for single-cell RNA-seq data. Nat Methods 2023;20:665–72. 10.1038/s41592-023-01814-1 - DOI - PMC - PubMed
    1. Aleksander SA, Balhoff J, Carbon S et al. ; Gene Ontology Consortium. The gene ontology knowledgebase in 2023. Genetics 2023;224:iyad031. 10.1093/genetics/iyad031 - DOI - PMC - PubMed