To estimate the tRNA adaptation index (tAI)
Genetic tRNA Adaptation index (gtAI)
gtAI is a new package implemented in python to effectively estimate the tRNA adaptation index (tAI).
For more information about the gtAI: Not yet published
For more information about the tAI: Mario dos Reis et. al.,.
Python >=3.7 is required.
pip install gtAI
Contributions to the software are welcome
For bugs and suggestions, the most effective way is by raising an issue on the github issue tracker. Github allows you to classify your issues so that we know if it is a bug report, feature request or feedback to the authors.
from gtAI import Run_gtAI df_tai, dict_wi, rel_values = Run_gtAI.gtai_analysis(main_fasta, GtRNA, genetic_code_number, size_pop, generation_number=50, ref_fasta= ref_fasta, bacteria=False)
main_fasta (str): A main fasta file containing the genes to be analyzed. GtRNA (dict): The tRNA genes count ref_fasta (str): Reference genes with the highest gene expression in a genome. genetic_code_number (int): default = 1, The Genetic Codes number described by NCBI (https://www.ncbi.nlm.nih.gov/Taxonomy/Utils/wprintgc.cgi) size_pop (int): A parameter for the genetic algorithm to identify the population size containing the possible solutions to optimize Sij-values. (default = 60) generation_number (int): A parameter for the genetic algorithm to identify the generation number. (default = 100) bacteria (bool): True If the tested organism is prokaryotic or archaeans, else equal to False (default = False)
Note: for ref_fasta parameter, the user is able to use a reference set of interest (in fasta format). Otherwise, the package will automatically generate a reference set based on the ENc values of the tested genome. For more information: API documentation.
Note: Population size must be an even number
df_tai (dataframe): Contains each gene id and its gtAI value. final_dict_wi (dict): Contains each codon and its absolute adaptiveness value. rel_values (dict): Contains each codon and its relative adaptiveness values.
1- Import gtAI functions.
from gtAI import Run_gtAI from gtAI import gtAI
2- In this example, we will use Saccharomyces cerevisiae S288C coding sequences.
3- Prepare the tRNA gene copy number of the tested genome.
The user has two options; a) input the tRNA gene copy number as python dictionary or, b) using GtRNAdb() function, the user can get it automatically from the GtRNA database, using the link to the tested genome (In our case Saccharomyces cerevisiae S288C). Or by tRNADB_CE() function to get the tRNA gene copy number from tRNADB_CE database using also the link to the tested genome.
In this example, the second option (b) will be used.
url_GtRNAdb = "http://gtrnadb.ucsc.edu/genomes/eukaryota/Scere3/" #### From GtRNAdb GtRNA = gtAI.GtRNAdb(url_GtRNAdb)
for more infromation about GtRNAdb() as well as tRNADB_CE(); API documentation.
4- Parameter settings for gtai_analysis() function.
genetic_code_number = 1 ref_fasta = "" bacteria = False size_pop = 60 generation_number = 100
for more information about gtai_analysis() and the parameters; API documentation.
5- Run gtAI.
df_tai , final_dict_wi, rel_values = Run_gtAI.gtai_analysis(main_fasta,GtRNA,genetic_code_number,bacteria=bacteria, size_pop=size_pop,generation_number=generation_number)
df_tai (dataframe): Contains each gene id and its gtAI value final_dict_wi (dict): Contains each codon and its absolute adaptiveness value rel_values (dict): Contains each codon and its relative adaptiveness values
6- To save the gtAI result as a CSV file.
import pandas as pd df_tai.to_csv("test.csv", header=True)
You can access the API documentation from here: gtAI Documentation
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
|Filename, size||File type||Python version||Upload date||Hashes|
|Filename, size gtAI-1.0.2.tar.gz (19.0 kB)||File type Source||Python version None||Upload date||Hashes View|