Multilocus sequence typing
Multilocus sequence typing (MLST) is a technique in molecular biology for the typing of multiple loci. The procedure characterizes isolates of microbial species using the DNA sequences of internal fragments of multiple housekeeping genes. Approximately 450-500 bp internal fragments of each gene are used, as these can be accurately sequenced on both strands using an automated DNA sequencer. For each housekeeping gene, the different sequences present within a bacterial species are assigned as distinct alleles and, for each isolate, the alleles at each of the loci define the allelic profile or sequence type (ST).
The first MLST scheme to be developed was for Neisseria meningitidis, the causative agent of meningococcal meningitis and septicaemia. Since its introduction for the research of evolutionary history, MLST has been used not only for human pathogens but also for plant pathogens.
MLST directly measures the DNA sequence variations in a set of housekeeping genes and characterizes strains by their unique allelic profiles. The principle of MLST is simple: the technique involves PCR amplification followed by DNA sequencing. Nucleotide differences between strains can be checked at a variable number of genes depending on the degree of discrimination desired.
The workflow of MLST involves: 1) data collection, 2) data analysis and 3) multilocus sequence analysis. In the data collection step, definitive identification of variation is obtained by nucleotide sequence determination of gene fragments. In the data analysis step, all unique sequences are assigned allele numbers and combined into an allelic profile and assigned a sequence type (ST). If new alleles and STs are found, they are stored in the database after verification. In the final analysis step of MLST, the relatedness of isolates are made by comparing allelic profiles. Researchers do epidemiological and phylogenetical studies by comparing STs of different clonal complexes. A huge set of data is produced during the sequencing and identification process so bioinformatic techniques are used to arrange, manage, analyze and merge all of the biological data.
To strike the balance between the acceptable identification power, time and cost for the strain typing, about seven to eight house-keeping genes are commonly used in the laboratories. Quoting Staphylococcus aureus as an example, seven housekeeping genes are used in MLST typing. These genes include carbamate kinase (arcC), shikimate dehydrogenase (aroE), glycerol kinase (glpF), guanylate kinase (gmk), phosphate acetyltransferase (pta), triosephosphate isomerase (tpi) and acetyl coenzyme A acetyltransferase (yqiL) as specified by the MLST website. However, it is not uncommon for up to ten housekeeping genes to be used. For Vibrio vulnificus, the housekeeping genes used are glucose-6-phosphate isomerase (glp), DNA gyrase, subunit B (gyrB), malate-lactate dehydrogenase (mdh), methionyl-tRNA synthetase (metG), phosphoribosylaminoimidazole synthetase (purM), threonine dehydrogenase (dtdS), diaminopimelate decarboxylase (lysA), transhydrogenase alpha subunit (pntA), dihydroorotase (pyrC) and tryptophanase (tnaA). Thus both the number and type of housekeeping genes interrogated by MLST may differ from species to species.
For each of these housekeeping genes, the different sequences are assigned as alleles and the alleles at the loci provide an allelic profile. A series of profiles can then be the identification marker for strain typing. Sequences that differ at even a single nucleotide are assigned as different alleles and no weighting is given to take into account the number of nucleotide differences between alleles, as we cannot distinguish whether differences at multiple nucleotide sites are a result of multiple point mutations or a single recombinational exchange. The large number of potential alleles at each of the loci provides the ability to distinguish billions of different allelic profiles, and a strain with the most common allele at each locus would only be expected to occur by chance approximately once in 10,000 isolates. Despite MLST providing high discriminatory power, the accumulation of nucleotide changes in housekeeping genes is a relatively slow process and the allelic profile of a bacterial isolate is sufficiently stable over time for the method to be ideal for global epidemiology.
The relatedness of isolates is displayed as a dendrogram constructed using the matrix of pairwise differences between their allelic profiles, eBURST or a minimum spanning tree (MST). The dendrogram is only a convenient way of displaying those isolates that have identical or very similar allelic profiles that can be assumed to be derived from a common ancestor; the relationships between isolates that differ at more than three out of seven loci are likely to be unreliable and should not be taken to infer their phylogeny. The MST connects all samples in such a way that the summed distance of all branches of the tree is minimal.
Alternatively, the relatedness of isolates can also be analysed with MultiLocus Sequence Analysis (MLSA). This does not use the assigned alleles, but instead concatenates the sequences of the gene fragments of the housekeeping genes and uses this concatenated sequence to determine phylogenetic relationships. In contrast to MLST, this analysis does assign a higher similarity between sequences differing only a single nucleotide and a lower similarity between sequences with multiple nucleotide differences. As a result, this analysis is more suitable for organisms with a clonal evolution and less suitable for organisms in which recombinational events occur very often. It can also be used to determine phylogenetic relationships between closely related species. The terms MLST and MLSA are very often considered interchangeable. This is however not correct as each analysis method has its distinctive features and uses. Care should be taken to use the correct term.
Comparison with other techniques
Earlier serological typing approaches had been established for differentiating bacterial isolates, but immunological typing has drawbacks such as reliance on few antigenic loci and unpredictable reactivities of antibodies with different antigenic variants. Several molecular typing schemes have been proposed to determine the relatedness of pathogens such as pulsed-field gel electrophoresis (PFGE), ribotyping, and PCR-based fingerprinting. But these DNA banding-based subtyping methods do not provide meaningful evolutionary analyses. Despite PFGE being considered by many researchers as the “gold standard”, many strains are not typable by this technique due to the degradation of the DNA during the process (gel smears).
The approach of MLST is distinct from Multi locus enzyme electrophoresis (MLEE), which is based on different electrophoretic mobilities (EM) of multiple core metabolic enzymes. The alleles at each locus define the EM of their products, as different amino acid sequences between enzymes result in different mobilities and distinct bands when run on a gel. The relatedness of isolates can then be visualized with a dendrogram generated from the matrix of pairwise differences between the electrophoretic types. This method has a lower resolution than MLST for several reasons, all arising from the fact that enzymatic phenotype diversity is merely a proxy for DNA sequence diversity. First, enzymes may have different amino acid sequences without having sufficiently different EM to give distinct bands. Second, "silent mutations" may alter the DNA sequence of a gene without altering the encoded amino acids. Thirdly, the phenotype of the enzyme can easily be altered in response to environmental conditions and badly affect the reproducibility of MLEE results - common modifications of enzymes are phosphorylation, cofactor binding and cleavage of transport sequences. This also limits comparability of MLEE data obtained by different laboratories, whereas MLST provides portable and comparable DNA sequence data and has great potential for automation and standardization.
MLST should not be confused with DNA barcoding. The latter is a taxonomic method that uses short genetic markers to recognize particular species of eukaryotes. It is based on the fact that mitochondrial DNA (mtDNA) or some parts of the ribosomal DNA cistron have relatively fast mutation rates, which give significant variation in sequences between species. mtDNA methods are only possible in eukaryotes (as prokaryotes lack mitochondria), whereas MLST, although initially developed for prokaryotes, is now finding application in eukaryotes and in principle could be applied to any kingdom.
Advantages and applications
MLST is highly unambiguous and portable. Materials required for ST determination can be exchanged between laboratories. Primer sequences and protocols can be accessed electronically. It is reproducible and scalable. MLST is automated, combines advances in high throughput sequencing and bioinformatics with established population genetics techniques. MLST data can be used to investigate evolutionary relationships among bacteria. MLST provides good discriminatory power to differentiate isolates.
The application of MLST is huge, and provides a resource for the scientific, public health, and veterinary communities as well as the food industry. The following are examples of MLST applications.
Campylobacter is the common causative agent for bacterial infectious intestinal diseases, usually arising from undercooked poultry or unpasteurised milk. However, its epidemiology is poorly understood since outbreaks are rarely detected, so that the sources and transmission routes of outbreak are not easily traced. In addition, Campylobacter genomes are genetically diverse and unstable with frequent inter- and intragenomic recombination, together with phase variation, which complicates the interpretation of data from many typing methods. Until recently, with the application of MLST technique, Campylobacter typing has achieved a great success and added onto the MLST database. As at 1 May 2008, the Campylobacter MLST database contains 3516 isolates and about 30 publications that use or mention MLST in research on Campylobacter (http://pubmlst.org/campylobacter/).
MLST has provided a more richly textured picture of bacteria within human populations and on strain variants that may be pathogenic to human, plants and animals. MLST technique was first used by Maiden et al. (1) to characterize Neisseria meningitidis using six loci. The application of MLST has clearly resolved the major meningococcal lineages known to be responsible for invasive disease around the world. To improve the level of discriminatory power between the major invasive lineages, seven loci are now being used and have been accepted by many laboratories as the method of choice for characterizing meningococcal isolates. It is a well known fact that recombinational exchanges commonly occur in N. meningitidis, leading to rapid diversification of meningococcal clones. MLST has successfully provided a reliable method for characterization of clones within other bacterial species in which the rates of clonal diversification are generally lower.
S. aureus causes a number of diseases. Methicillin-resistant S. aureus (MRSA) has generated growing concerns over its resistance to almost all antibiotics except vancomycin. However, most serious S. aureus infections in the community, and many in hospitals, are caused by methicillin-susceptible isolates (MSSA) and there have been few attempts to identify the hypervirulent MSSA clones associated with serious disease. MLST was therefore developed to provide an unambiguous method of characterizing MRSA clones and for the identification of the MSSA clones associated with serious disease.
S. pyogenes causes diseases ranging from pharyngitis to life-threatening impetigo including necrotizing fasciitis. An MLST scheme for S. pyogenes has been developed. At present, the database (mlst.net) contains the allelic profiles of isolates that represent the worldwide diversity of the organism and isolates from serious invasive disease.
C. albicans is a fungal pathogen of humans and is responsible for hospital-acquired bloodstream infections. MLST technique has used to characterize C. albicans isolates. Combination of the alleles at the different loci results in unique diploid sequence types that can be used to discriminate strains. MLST has been shown successfully applied to study the epidemiology of C. albicans in the hospital as well as the diversity of C. albicans isolates obtained from diverse ecological niches including human and animal hosts.
The genus Cronobacter is composed of 7 species. Before 2007, the single species name Enterobacter sakazakii was applied to these organisms. The Cronobacter MLST was initially applied to distinguish between C. sakazakii and C. malonaticus because 16S rDNA sequencing is not always accurate enough, and biotyping is too subjective. The Cronobacter MLST scheme uses 7 alleles; atpD, fusA, glnS, gltB, gyrB, infB and ppsA giving a concatenated sequence of 3036 bp for phylogenetic analysis (MLSA) and comparative genomics. MLST has also been used in the formal recognition of new Cronobacter species. The method has revealed a strong association between one genetic lineage, sequence type 4 (ST4), and cases of neonatal meningitis., The Cronobacter MLST site is at http://www.pubMLST.org/cronobacter.
MLST appears best in population genetic study but it is expensive. Due to the sequence conservation in housekeeping genes, MLST sometimes lacks the discriminatory power to differentiate bacterial strains, which limits its use in epidemiological investigations. To improve the discriminatory power of MLST, a multi-virulence-locus sequence typing (MVLST) approach has been developed using Listeria monocytogenes . MVLST broadens the benefits of MLST but targets virulence genes, which may be more polymorphic than housekeeping genes. Population genetics is not the only relevant factor in an epidemic. Virulence factors are also important in causing disease, and population genetic studies struggle to monitor these. This is because the genes involved are often highly recombining and mobile between strains in comparison with the population genetic framework. Thus, for example in Escherichia coli, identifying strains carrying toxin genes is more important than having a population genetics-based evaluation of prevalent strains.
The advent of second-generation sequencing technologies has made it possible to obtain sequence information across the entire bacterial genome at relatively modest cost and effort, and MLST can now be assigned from whole-genome sequence information, rather than sequencing each locus separately as was the practice when MLST was first developed. Whole-genome sequencing provides richer information for differentiating bacterial strains (MLST uses approximately 0.1% of the genomic sequence to assign type while disregarding the rest of the bacterial genome). For example, whole-genome sequencing of numerous isolates has revealed the single MLST lineage ST258 of Klebsiella pneumoniae comprises two distinct genetic clades, providing additional information about the evolution and spread of these multi-drug resistant organisms, and disproving the previous hypothesis of a single clonal origin for ST258.
MLST databases contain the reference allele sequences and sequence types for each organism, and also isolate epidemiological data. The websites contain interrogation and analysis software which allow users to query their allele sequences and sequence types. MLST is widely used as a tool for researchers and public healthcare workers.
The databases hosted at each site are different and hold the organism specific reference allele sequences and lists of STs for individual organisms.
To assist the gathering and formatting of the utilized sequences a simple and free plug-in for Firefox has been developed (link).
- Maiden, MC.; Bygraves, JA.; Feil, E.; Morelli, G.; Russell, JE.; Urwin, R.; Zhang, Q.; Zhou, J.; et al. (Mar 1998). "Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms". Proc Natl Acad Sci U S A. 95 (6): 3140–5. Bibcode:1998PNAS...95.3140M. doi:10.1073/pnas.95.6.3140. PMC 19708. PMID 9501229.
- Sarris, PF; Trantas, EA; Mpalantinaki, E; Ververidis, F; Goumas, DE (2012). "Pseudomonas viridiflava, a multi host plant pathogen with significant genetic variation at the molecular level". PLOS ONE. 7 (4): e36090. Bibcode:2012PLoSO...736090S. doi:10.1371/journal.pone.0036090. PMC 3338640. PMID 22558343.
- Spratt, Brian G (1999). "Multilocus sequence typing: molecular typing of bacterial pathogens in an era of rapid DNA sequencing and the Internet". Current Opinion in Microbiology. 2 (3): 312–316. doi:10.1016/S1369-5274(99)80054-X. PMID 10383857.
- Spratt, BG; Maiden, MC (1999). "Bacterial population genetics, evolution and epidemiology". Philos Trans R Soc Lond B Biol Sci. 354 (1384): 701–710. doi:10.1098/rstb.1999.0423. PMC 1692550. PMID 10365396.
- Jolley, Keith. "Minimum Spanning Tree". pubMLST. Archived from the original on 15 September 2007. Retrieved 10 October 2013.
- Gevers D, Cohan FM, Lawrence JG, Spratt BG, Coenye T, Feil EJ, Stackebrandt E, Van de Peer Y, Vandamme P, Thompson FL, Swings J (2005). "Opinion: Re-evaluating prokaryotic species". Nat Rev Microbiol. 3 (9): 733–739. doi:10.1038/nrmicro1236. PMID 16138101.
- E J Feil; M C Maiden; M Achtman; B G Spratt (1999). "The relative contributions of recombination and mutation to the divergence of clones of Neisseria meningitidis". Mol Biol Evol. 16 (11): 1496–1502. doi:10.1093/oxfordjournals.molbev.a026061. PMID 10555280.
- Enright MC, Spratt BG, Kalia A, Cross JH, Bessen DE (2001). "Multilocus sequence typing of Streptococcus pyogenes and the relationships between emm type and clone". Infect Immun. 69 (4): 2416–2427. doi:10.1128/iai.69.4.2416-2427.2001. PMC 98174. PMID 11254602.
- McGregor KF, Spratt BG, Kalia A, Bennett A, Bilek N, Beall B, Bessen DE (2004). "Multilocus sequence typing of Streptococcus pyogenes representing most known emm types and distinctions among subpopulation genetic structures". J Bacteriol. 186 (13): 4285–4294. doi:10.1128/jb.186.13.4285-4294.2004. PMC 421626. PMID 15205431.
- Baldwin; et al. (2009). "Multilocus sequence typing of Cronobacter sakazakii and Cronobacter malonaticus reveals stable clonal structures with clinical significance which do not correlate with biotypes". BMC Microbiology. 9: 223. doi:10.1186/1471-2180-9-223. PMC 2770063. PMID 19852808.
- Kucerova; et al. (2011). "Cronobacter: diversity and ubiquity". Quality Assurance and Safety of Foods and Crops. 3 (3): 104–122. doi:10.1111/j.1757-837X.2011.00104.x.
- Joseph; et al. (2011). "Cronobacter condimenti sp. nov., isolated from spiced meat and Cronobacter universalis sp. nov., a novel species designation for Cronobacter sp. genomospecies 1, recovered from a leg infection, water, and food ingredients". Int J Syst Evol Microbiol. 62 (Pt 6): 1277–83. doi:10.1099/ijs.0.032292-0. PMID 22661070.
- Joseph and Forsythe (2011). "Association of Cronobacter sakazakii ST4 with neonatal infections". Emerging Infectious Diseases. 17 (9): 1713–5. doi:10.3201/eid1709.110260. PMC 3322087. PMID 21888801.
- Zhang, W.; Jayarao, B. M.; Knabel, S. J. (2004). "Multi-Virulence-Locus Sequence Typing of Listeria monocytogenes". Applied and Environmental Microbiology. 70 (2): 913–920. doi:10.1128/AEM.70.2.913-920.2004. ISSN 0099-2240. PMC 348834. PMID 14766571.
- "Center for Genomic Epidemiology".
- DeLeo, F.; et al. (2014). "Molecular dissection of the evolution of carbapenem-resistant multilocus sequence type 258 Klebsiella pneumoniae". Proceedings of the National Academy of Sciences. 111 (13): 4988–93. Bibcode:2014PNAS..111.4988D. doi:10.1073/pnas.1321364111. PMC 3977278. PMID 24639510.
- Woodford N, Turton JF, Livermore DM; Turton; Livermore (2011). "Multiresistant Gram-negative bacteria: the role of high-risk clones in the dissemination of antibiotic resistance". FEMS Microbiology Reviews. 35 (5): 736–55. doi:10.1111/j.1574-6976.2011.00268.x. PMID 21303394.CS1 maint: multiple names: authors list (link)
- mlst.net - Imperial College London
- PubMLST - Oxford University
- databases hosted at University College Cork
- databases held at Pasteur Institute
- BioNumerics One universal bioinformatics solution to store and analyze all your biological data. The entire MLST workflow can be performed and results compared to other typing methods and/or metadata. Allele sequences and IDs can also be derived from whole genome sequences.
- CLC Microbial Genomics Module includes a set of tools and ready-to-use workflows for MLST in the context of meta-information such as outbreak, host, geographical location, etc., and allows easy comparison to typing results obtained using other typing methods.
- Maiden MC, Bygraves JA, Feil E, et al. (March 1998). "Multilocus sequence typing: A portable approach to the identification of clones within populations of pathogenic microorganisms". Proc. Natl. Acad. Sci. U.S.A. 95 (6): 3140–5. Bibcode:1998PNAS...95.3140M. doi:10.1073/pnas.95.6.3140. PMC 19708. PMID 9501229.
- Urwin R, Maiden MC; Maiden (October 2003). "Multi-locus sequence typing: a tool for global epidemiology". Trends Microbiol. 11 (10): 479–87. doi:10.1016/j.tim.2003.08.006. PMID 14557031.
- Foley SL, White DG, McDermott PF, et al. (October 2006). "Comparison of Subtyping Methods for Differentiating Salmonella enterica Serovar Typhimurium Isolates Obtained from Food Animal Sources". J. Clin. Microbiol. 44 (10): 3569–77. doi:10.1128/JCM.00745-06. PMC 1594788. PMID 17021084.
- Johnson JK, Arduino SM, Stine OC, Johnson JA, Harris AD; Arduino; Stine; Johnson; Harris (November 2007). "Multilocus Sequence Typing Compared to Pulsed-Field Gel Electrophoresis for Molecular Typing of Pseudomonas aeruginosa". J. Clin. Microbiol. 45 (11): 3707–12. doi:10.1128/JCM.00560-07. PMC 2168506. PMID 17881548.CS1 maint: multiple names: authors list (link)
- Nallapareddy SR, Duh RW, Singh KV, Murray BE; Duh; Singh; Murray (March 2002). "Molecular Typing of Selected Enterococcus faecalis Isolates: Pilot Study Using Multilocus Sequence Typing and Pulsed-Field Gel Electrophoresis". J. Clin. Microbiol. 40 (3): 868–76. doi:10.1128/JCM.40.3.868-876.2002. PMC 120268. PMID 11880407.CS1 maint: multiple names: authors list (link)
- Fakhr MK, Nolan LK, Logue CM; Nolan; Logue (May 2005). "Multilocus Sequence Typing Lacks the Discriminatory Ability of Pulsed-Field Gel Electrophoresis for Typing Salmonella enterica Serovar Typhimurium". J. Clin. Microbiol. 43 (5): 2215–9. doi:10.1128/JCM.43.5.2215-2219.2005. PMC 1153745. PMID 15872244.CS1 maint: multiple names: authors list (link)
- Chen Y, Zhang W, Knabel SJ; Zhang; Knabel (October 2005). "Multi-Virulence-Locus Sequence Typing Clarifies Epidemiology of Recent Listeriosis Outbreaks in the United States". J. Clin. Microbiol. 43 (10): 5291–4. doi:10.1128/JCM.43.10.5291-5294.2005. PMC 1248515. PMID 16208000.CS1 maint: multiple names: authors list (link)