| title | tags | authors | date | bibliography | |||||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Multilocus sequence typing by blast from de novo assemblies against PubMLST |
|
|
3 Nov 2016 |
paper.bib |
Multilocus sequence typing (MLST) is a standard method for categorising genomes [@Maiden1998] based on variation in a small set of conserved house keeping genes. It allows for rapid identification of genomes into high level categories and is extremely useful for epidemiological investigations [@Urwin2003] making it a key tool for public health reference laboratories. Whilst MLST is more commonly associated with classical sequencing methods, it is possible to extract the same information from Next Generation Sequencing data, in particular from de novo assemblies which are generated routinely for bacterial sequencing data [@PAGE2016].
We provide a scalable command line tool, MLSTcheck, which can take multiple de novo assemblies and output detailed information about the sequence type of the samples. It provides access to 124 MLST databases covering all of the major human disease causing bacterial pathogens. MLSTcheck can search one or more databases at once, is parallelisable, fast and robust. When a sample contains more than one allele, it flags the contaminant since there should only be 1 copy of a house keeping gene in a well designed MLST scheme. A multiple FASTA alignment of the concatenated MLST genes is optionally produced, allowing for the creation of phylogenetic trees. This allows for rapid epidemiological outbreak investigations. Whilst other software applications can perform similar functions [@Seemann2016; @Jolley2010], this application follows more rigorous software engineering principles, including automated testing, continuous integration, object orientated code, and is installable via CPAN (a Perl package manager). In a large diverse set of 6814 publicly accessible draft assemblies, MLSTcheck was able to assign a sequence type in 99.6% of cases [@PAGE2016].