Package: genetorch
Gene Torch 1.3.5REALLY REALLY A LOT OF UPDATESinclude a protein folder which contains all alphafold model of C.elegans protein
with .apf(Aminoacid position file)You can read it using module pcluster.What is new in version 1.2.0genetorch.finder.plotpro(readfile,size,c='Reds',sub_c='RdPu',origin = a)A more professional version to distinguish gene via binominal distribution instead of overlap. The calculation is a bit slow but more accurate.
The genes you tried will be savedin a.candidateand all the sample names of the candidate gene will be listed ina.suppressor_group.
You should also notice that in the previous version,genetorch.finder.plot(result,size,c,sub_c)you shall input a result file. In this function, you may need to input the whole readfile() object or getfile() object instead.genetorch.reader.get_impact(readfile)This function can now annotate splice_donor and splice_acceptor variation with 'X_donor' and 'X_acceptor' because we noticed that splice variation do have an impact on protein function.genetorch.finder.get_p_sup(gene,readfile)return a Dataframe including all the p values the gene you chose vs others.genetorch.finder.get_p_between(readfile,genea,geneb)return the p value between genea and genebgenetorch.finder.find(readfile)/filter(readfile,lengthlimit,rid)We canceled the threshold in 'filter' because we noticed that it is not very useful and can make our users confused. These two functions will automatically do agenetorch.reader.get_impactnow.
Also, in some cases you may have to reset the  candidate and suppressor_group list. You only need to do a filter or find to do this.genetorch.simulator.false_positive(readfile,gene,[0.05,0.04,0.03,0.02,0.01])Give the estimated value of the false positive cases using the gene you chose to compare with others. This step will cost about 5 mins for each value.What is new in version 1.1.11. 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 settingrid = [], 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 graphsub_c = color map used in other graphsintersection_num = 0.03For 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 graphOther changesgenetorch.finder.intersection(result,samplelist)You can use this function to find genes not in these samples. Return a Dataframegenetorch.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 change1.0.3 update:bug FIXING........filter is more accurate and fastersolve the problem that file cannot be decoded by 'utf-8' in MacOSpip install genetorchPLEASE INSTALL: pandas, matplotlib, seaborn BEFORE USEThis 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]:guozhengyang980525@yahoo.co.jpgenetorch.readergenetorch.reader.readfile(filepath)multiple renamed vcf files must be included in the filepath.filepath
|---1.vcf
|---2.vcf
|---3.vcf
|---4.vcf
|---5.vcf
|---6.vcfgenetorch.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
|---cas114_20221011jxskaosdosh---filename.vcf
|---cas115_20221011jxskaosdosh---filename.vcf
|---cas116_20221011jxskaosdosh---filename.vcf
|---cas117_20221011jxskaosdosh---filename.vcf
|---cas118_20221011jxskaosdosh---filename.vcfa temp folder will be automatically created in the filepath including renamed vcf files:filepath
|----temp
|     |---cas113.vcf
|     |---cas114.vcf
|     |---cas115.vcf
|     |---cas116.vcf
|     |---cas117.vcf
|     |---cas118.vcg
|---cas113_20221011jxskaosdosh---filename.vcf
|---cas114_20221011jxskaosdosh---filename.vcf
|---cas115_20221011jxskaosdosh---filename.vcf
|---cas116_20221011jxskaosdosh---filename.vcf
|---cas117_20221011jxskaosdosh---filename.vcf
|---cas118_20221011jxskaosdosh---filename.vcfa = 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]:geneIDtypebaseproteintagttn-1WBGenexxxxmissenseC<GAsp666Asncas113cla-1WBGenexxxxmissenseC<GAsp223Asncas114genetorch.findergenetorch.finder.find(readfile())return a pandas Dataframe item.example:samplesample sizegenevariationvariation_number--------------------------------------cas113,cas1142ttn-1Asp666Asn,1 Glu374Gln,22ttn-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(readfile(),lengthlimit = 0.6,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 variationif the gene has no variation left after filtered, the gene will be deleted from the result Dataframegenetorch.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
automatically show annotation of each bubblegenetorch.stockerModule for data upload to LIM databasegenetorch.stocker.stock(path,filename,outpath)path & filenamemultiple 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
|---cas114_20221011jxskaosdosh---filename.vcf
|---cas115_20221011jxskaosdosh---filename.vcf
|---cas116_20221011jxskaosdosh---filename.vcf
|---cas117_20221011jxskaosdosh---filename.vcf
|---cas118_20221011jxskaosdosh---filename.vcfa temp folder will be automatically created in the filepath including renamed vcf files:filepath
|----temp
|     |---cas113.vcf
|     |---cas114.vcf
|     |---cas115.vcf
|     |---cas116.vcf
|     |---cas117.vcf
|     |---cas118.vcg
|---cas113_20221011jxskaosdosh---filename.vcf
|---cas114_20221011jxskaosdosh---filename.vcf
|---cas115_20221011jxskaosdosh---filename.vcf
|---cas116_20221011jxskaosdosh---filename.vcf
|---cas117_20221011jxskaosdosh---filename.vcf
|---cas118_20221011jxskaosdosh---filename.vcfoutpath:outpath is an empty folder.
output csv files in this folder which can be uploaded to LIM database.Examplesimport 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
gt.reader.get_impact(a) # get impact results from a
d = gt.finder.find(a) # get the result table
e = gt.finder.filter(a,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