genetorch 1.1.1

To deal with large amount of WGS data in suppressor screening and to predict suppressor gene

Gene Torch 1.1.1

What is new in version 1.1.1

1. You can get rid of some annoying background variation by just removing them when 'filter'

genetorch.finder.filter(taglist,lengthlimit,threshold,rid = ['ttn-1','cla-1'])

by setting rid = [] , the genes in the list will no longer exist in your results.

2. You can use the intersection function to find multiple candidate suppressors.

genetorch.finder.plot(result,size,c='Reds',sub_c='RdPu',intersection_factor = 0.03)

c = color map used in the major graph

sub_c = color map used in other graphs

intersection_num = 0.03

For example: If you think gene A is a suppressor, it means that in samples with other mutations,there is less opportunity for them to carry a gene A mutation. intersection_num means the least percent (default 3%) of the sample number you allow to carry more than one suppressor.

You can now find the intersection by click on the bubble of the gene on the graph

Other changes

genetorch.finder.intersection(result,samplelist)

You can use this function to find genes not in these samples. Return a Dataframe

genetorch.finder.gsamplename(result,gene)

You can get a list of your sample carry this gene.

genetorch.reader.get_impact(taglist)

Remove the lines in your data which do not have an amino_acid change

1.0.3 update:

bug FIXING........

filter is more accurate and faster

solve the problem that file cannot be decoded by 'utf-8' in MacOS

pip install genetorch

PLEASE INSTALL: pandas, matplotlib, seaborn BEFORE USE

This package can help find candidate genes in high throughput mutagenesis and suppressor screening experiments without mapping. Please call variants via freebayes and annotate with snpeff before using this package. Any advise is welcomed, please contact

[e-mail]:[e-mail]

genetorch.reader

genetorch.reader.readfile(filepath)

multiple renamed vcf files must be included in the filepath. examples:

filepath/1.vcf

filepath/2.vcf

filepath/3.vcf

genetorch.reader.getfile(filepath,filename)

multiple renamed folders must be included in the filepath, and a vcf file with filename must be included in the folders, the name of the folder must be splited with '_' to divide the folder name into strain name and WGS order name: examples:

filepath/cas113_20221011jxskaosdosh/filename.vcf

filepath/cas114_20221011jxskaosdosh/filename.vcf

filepath/cas115_20221011jxskaosdosh/filename.vcf

filepath/cas116_20221011jxskaosdosh/filename.vcf

filepath/cas117_20221011jxskaosdosh/filename.vcf

filepath/cas118_20221011jxskaosdosh/filename.vcf

a temp folder will be automatically created in the filepath including renamed vcf files:

filepath/temp/cas113.vcf

filepath/temp/cas114.vcf

filepath/temp/cas115.vcf

filepath/temp/cas116.vcf

filepath/temp/cas117.vcf

filepath/temp/cas118.vcf

a = genetorch.readfile() a = genetorch.getfile()

a.taglist : a list of Dataframes which included columns: 'gene', 'ID', 'type', 'base', 'protein'，'tag' column 'tag' will be filled with strain name

example: a.taglist[1]:

gene	ID	type	base	protein	tag
ttn-1	WBGenexxxx	missense	C<G	Asp666Asn	cas113
cla-1	WBGenexxxx	missense	C<G	Asp223Asn	cas114

genetorch.finder

genetorch.finder.find(taglist)

return a pandas Dataframe item. example:

sample	sample size	gene	variation	variation_number
cas113,cas114	2	ttn-1	Asp666Asn,1 Glu374Gln,2	2

ttn-1 mutation is found in 2 input files, there are three amino acid variation events, 2 * Glu374Gln, 1 * Asp666Asn. There are two kinds of variations in this gene.

genetorch.finder.filter(taglist,lengthlimit = 0.6,threshold = 1,rid = ['ttn-1','cla-1'])

return a pandas Dataframe item same as find().

provide filtered data:

lengthlimit: if n(amino acid variation) > lengthlimit * n(input files), exclude the variation

threshold: if n(variation_number) > threshold * n(input files) , exclude the gene

if the gene has no variation left after filtered, the gene will be deleted from the result Dataframe

genetorch.finder.plot(result,size,intersection_num)

result: result Dataframe from find() or filter() size: use the first n lines to show the plot.

x: gene name

y:variation number

size: size

color: variation number/size

click to see annotation of each bubble

genetorch.stocker

Module for data upload to LIM database

genetorch.stocker.stock(path,filename,outpath)

path & filename

multiple renamed folders must be included in the filepath, and a vcf file with filename must be included in the folders, the name of the folder must be splited with '_' to divide the folder name into strain name and WGS order name: examples:

filepath/cas113_20221011jxskaosdosh/filename.vcf

filepath/cas114_20221011jxskaosdosh/filename.vcf

filepath/cas115_20221011jxskaosdosh/filename.vcf

filepath/cas116_20221011jxskaosdosh/filename.vcf

filepath/cas117_20221011jxskaosdosh/filename.vcf

filepath/cas118_20221011jxskaosdosh/filename.vcf

a temp folder will be automatically created in the filepath including renamed vcf files:

filepath/temp/cas113.vcf

filepath/temp/cas114.vcf

filepath/temp/cas115.vcf

filepath/temp/cas116.vcf

filepath/temp/cas117.vcf

filepath/temp/cas118.vcf

outpath:

outpath is an empty folder. output csv files in this folder which can be uploaded to LIM database.

Examples

import genetorch as gt
a = gt.reader.readfile(path) # read vcf files from path
b = gt.reader.getfile(path,filename) # read vcf files with the filename in the individual folders in the path 
                                             # a has the same properties as b
c = gt.reader.get_impact(a.taglist) # get impact results from a
d = gt.finder.find(c) # get the result table
e = gt.finder.filter(c,lengthlimit = 0.6,threshold = 1, rid = ['ttn-1','cla-1']) # get filtered result table
gt.finder.plot(e,1000,1) # draw a picture using the first 1000 lines in the result table,
# allowing 1 sample have more than one suppressor gene
gt.stocker.stockfile(path,filename,outpath) #get stock file