Abstract:Understanding the regulation of gene expression is one of the key problems in current biology. A promising method for that purpose is the determination of the temporal dynamics between known initial and ending network states, by using simple acting rules. The huge amount of rule combinations and the nonlinear inherent nature of the problem make genetic algorithms an excellent candidate for finding optimal solutions. As this is a computationally intensive problem that needs long runtimes in conventional architectures for realistic network sizes, it is fundamental to accelerate this task. In this article, we study how to develop efficient parallel implementations of this method for the fine-grained parallel architecture of graphics processing units (GPUs) using the compute unified device architecture (CUDA) platform. An exhaustive and methodical study of various parallel genetic algorithm schemes-master-slave, island, cellular, and hybrid models, and various individual selection methods (roulette, elitist)-is carried out for this problem. Several procedures that optimize the use of the GPU's resources are presented. We conclude that the implementation that produces better results (both from the performance and the genetic algorithm fitness perspectives) is simulating a few thousands of individuals grouped in a few islands using elitist selection. This model comprises 2 mighty factors for discovering the best solutions: finding good individuals in a short number of generations, and introducing genetic diversity via a relatively frequent and numerous migration. As a result, we have even found the optimal solution for the analyzed gene regulatory network (GRN). In addition, a comparative study of the performance obtained by the different parallel implementations on GPU versus a sequential application on CPU is carried out. In our tests, a multifold speedup was obtained for our optimized parallel implementation of the method on medium class GPU over an equivalent sequential single-core implementation running on a recent Intel i7 CPU. This work can provide useful guidance to researchers in biology, medicine, or bioinformatics in how to take advantage of the parallelization on massively parallel devices and GPUs to apply novel metaheuristic algorithms powered by nature for real-world applications (like the method to solve the temporal dynamics of GRNs).

TurboGP: A flexible and advanced python based GP library

Genetic Programming Based on Non-Dominated Sorting

GP+: A Python Library for Kernel-based learning via Gaussian Processes

General Purpose Optimization Library (GPOL): A Flexible and Efficient Multi-Purpose Optimization Library in Python

Exploiting Tournament Selection for Efficient Parallel Genetic Programming

TensorGP -- Genetic Programming Engine in TensorFlow

Imperative Genetic Programming

GPfit: An R package for Gaussian Process Model Fitting using a New Optimization Algorithm

PyGAD: an intuitive genetic algorithm Python library

GSGP-CUDA -- a CUDA framework for Geometric Semantic Genetic Programming

GPflowOpt: A Bayesian Optimization Library using TensorFlow

An Empirical Study of Progressive Insular Cooperative GP

Graphics Processing Unit–Enhanced Genetic Algorithms for Solving the Temporal Dynamics of Gene Regulatory Networks

GA-lapagos, an open-source c framework including a python-based system for data analysis

A segregated genetic programming for bioprocess modelling with outliers

CGP++ : A Modern C++ Implementation of Cartesian Genetic Programming

MOGPTK: The multi-output Gaussian process toolkit

HOTGP -- Higher-Order Typed Genetic Programming

Genetic Programming to Optimize 3D Trajectories

Code Building Genetic Programming