Skip to content

Data and code to accompany our paper on identifying patterns of gene degradation associated with invasiveness in Salmonella

License

Notifications You must be signed in to change notification settings

Gardner-BinfLab/invasive_salmonella

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

75 Commits
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

invasive_salmonella

This repository contains data and code to accompany our paper on identifying patterns of gene degradation associated with invasiveness in Salmonella.

The method relies on the DeltaBS metric (https://github.com/Gardner-BinfLab/deltaBS) for identifying deleterious mutations in protein-coding genes. The DeltaBS approach uses profile hidden Markov models to take a diverse collection of sequences of the same gene from different organisms and capture patterns of sequence variation that occur commonly in nature (in this case, within the Gammaproteobacteria). These models are then used to score the protein coding genes of different strains of bacteria to produce bitscores, which indicate how well each protein conforms to modelled sequence constraints on that gene. The distribution of these bitscores can then be used by the random forest to identify genes which show a difference in adherence to sequence constraints between strains from different niches. A difference in bitscore distribution is likely to be caused by a change in selective pressures on a gene in one niche, leading to increased accumulation of mutations not commonly observed in nature. This accumulation of rare mutations is likely to result in the degradation of the gene, however these mutations could also result in a change to the function or structural stability of the protein. For a given protein family, the difference in bitscore for one strain and the median bitscore for the strain collection (DeltaBS) can be used as an indication of the likelihood that this accumulation of mutations is deleterious.

The repository is separated into three directories:

  • One for building your own model for identifying genes associated with a phenotype.
  • One containing code for the original analysis.
  • One containing code to analyse your own Salmonella strains using our model.

See the Wiki page for instructions on building your own random forest model or testing your own Salmonella isolates using our invasiveness index: https://github.com/UCanCompBio/invasive_salmonella/wiki

About

Data and code to accompany our paper on identifying patterns of gene degradation associated with invasiveness in Salmonella

Resources

License

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published