Skip to main content

A python package for extracting tevolutionary trajectories from genotype-phenotype-maps

Project description

gpvolve

A python package for extracting tevolutionary trajectories from genotype-phenotype-maps

A Python API for the simulation and analysis of evolution in genotype-phenotype space. You can use this library to:

  1. Build a markov state model from a genotype-phenotype-map.
  2. Find clusters of genotypes that represent metastable states of the system, using PCCA+.
  3. Compute fluxes and pathways between pairs of genotypes and/or clusters of interest, using Transition Path Theory.
  4. Visualize the outputs of all of the above.

Basic Example

Build a Markov model from a genotype-phenotype map.

# Import base class, Transition Path Theory class and functions for building Markov Model.
from gpvolve import GenotypePhenotypeMSM, TransitionPathTheory, linear_skew, mccandlish, find_max

# Import visualization tool.
from gpvolve.visualization import plot_network

# Import GenotypePhenotypeMap class for handling genotype-phenotype data.
from gpmap import GenotypePhenotypeMap

# Helper functions.
from scipy.sparse import dok_matrix

# Genotype-phenotype map data.
wildtype = "AAA"
genotypes = ["AAA", "AAT", "ATA", "TAA", "ATT", "TAT", "TTA", "TTT"]
phenotypes = [0.8, 0.81, 0.88, 0.89, 0.82, 0.82, 0.95, 1.0]

# Instantiate Markov model class.
gpm = GenotypePhenotypeMap(wildtype=wildtype,
                            genotypes=genotypes,
                            phenotypes=phenotypes)


# Instantiate a evolutionary Markov State Model from the genotype-phenotype map.
gpmsm = GenotypePhenotypeMSM(gpm)

Apply an evolutionary model to describe transitions between genotypes.

# Map fitnesses to phenotypes.
gpmsm.apply_selection(fitness_function=linear_skew, selection_gradient=1)

# Build Markov State Model based on 'mccandlish' fixation probability function.
gpmsm.build_transition_matrix(fixation_model=mccandlish, population_size=100)

# Find global fitness peak.
fitness_peak = find_max(gpmsm=gpmsm, attribute='fitness')

Calculate and plot the fluxes between wildtype and triple mutant.

# Compute fluxes from wildtype to fitness peak.
fluxes = TransitionPathTheory(gpmsm, source=[0], target=[fitness_peak])

# Normalize flux.
norm_fluxes = fluxes.net_flux/fluxes.total_flux

# Plot the network and the fluxes
fig, ax = plot_network(gpmsm,
                        flux=dok_matrix(norm_fluxes),
                        edge_labels=True,
                        colorbar=True)

Install

To install from PyPI:

pip install gpvolve

To install a development version:

git clone https://github.com/harmslab/gpvolve
cd gpvolve
pip install  -e .

Project details


Release history Release notifications

Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Files for gpvolve, version 0.0.1
Filename, size File type Python version Upload date Hashes
Filename, size gpvolve-0.0.1-py2.py3-none-any.whl (24.3 kB) File type Wheel Python version py2.py3 Upload date Hashes View hashes
Filename, size gpvolve-0.0.1.tar.gz (1.4 MB) File type Source Python version None Upload date Hashes View hashes

Supported by

Elastic Elastic Search Pingdom Pingdom Monitoring Google Google BigQuery Sentry Sentry Error logging AWS AWS Cloud computing DataDog DataDog Monitoring Fastly Fastly CDN DigiCert DigiCert EV certificate StatusPage StatusPage Status page