Automated pneumococcal MLST using liquid-handling robotics and a capillary DNA sequencer

Multilocus sequence typing (MLST) is used by the Scottish Meningococcus and Pneumococcus Reference Laboratory (SMPRL) as a routine method for the characterization of certain bacterial pathogens. The SMPRL recently started performing MLST on strains of Streptococcus pneumoniae, and here we describe a fully automated method for MLST using a 96-well-format liquid-handling robot and a 96-capillary automated DNA sequencer.     

Short read sequence typing (SRST): multi-locus sequence types from short reads

ABSTRACT: BACKGROUND: Multi-locus sequence typing (MLST) has become the gold standard for population analyses of bacterial pathogens. This method focuses on the sequences of a small number of loci (usually seven) to divide the population and is simple, robust and facilitates comparison of results between laboratories and over time. Over the last decade, researchers and population health specialists have invested substantial effort in building up public MLST databases for nearly 100 different bacterial species, and these databases contain a wealth of important information linked to MLST sequence types such as time and place of isolation of isolation, host or niche, serotype and even clinical or drug resistance profiles. Recent advances in sequencing technology mean it is increasingly feasible to perform bacterial population analysis at the whole genome level. This offers massive gains in resolving power and genetic profiling compared to MLST, and will eventually replace MLST for bacterial typing and population analysis. However given the wealth of data currently available in MLST databases, it is crucial to maintain backwards compatibility with MLST schemes so that new genome analyses can be understood in their proper historical context. RESULTS: We present a software tool, SRST, for quick and accurate retrieval of sequence types from short read sets, using inputs easily downloaded from public databases. SRST assigns alleles using read mapping and an allele assignment score incorporating sequence coverage and variability, to determine the most likely allele. Analysis of over 3,500 loci in more than 500 publicly accessible Illumina read sets showed SRST to be highly accurate at allele assignment. SRST output is compatible with common analysis tools such as eBURST, Clonal Frame or PhyloViz, allowing easy comparison between novel genome data and MLST data. Alignment, fastq and pileup files can also be generated for novel alleles. CONCLUSIONS: SRST is a novel software tool for accurate assignment of sequence types using short read data. Several uses for the tool are demonstrated, including quality control for high-throughput sequencing projects, plasmid MLST and analysis of genomic data during outbreak investigation. SRST is open-source, requires Python, BWA and SamTools, and is available from http://srst.sourceforge.net.     

Automation of MLST using third-generation liquid-handling technology

The molecular characterization of bacterial pathogens of clinical significance is increasingly important. Methods, such as multilocus sequence typing (MLST), allow bacterial strains to be characterized during case clusters, for antibiotic-resistant strains to be monitored, and for the impact of new vaccines to be assessed. Our laboratory performs MLST on Neisseria meningitidis, Streptococcus pneumonaie, Haemophilus influenzae, and Staphylococcus aureus. We have developed high-throughput automated methods to allow MLST to be performed in a time scale useful in a clinical setting. Here we describe the automation of MLST on a third-generation liquid-handling robot.     

eBURST: inferring patterns of evolutionary descent among clusters of related bacterial genotypes from multilocus sequence typing data

The introduction of multilocus sequence typing (MLST) for the precise characterization of isolates of bacterial pathogens has had a marked impact on both routine epidemiological surveillance and microbial population biology. In both fields, a key prerequisite for exploiting this resource is the ability to discern the relatedness and patterns of evolutionary descent among isolates with similar genotypes. Traditional clustering techniques, such as dendrograms, provide a very poor representation of recent evolutionary events, as they attempt to reconstruct relationships in the absence of a realistic model of the way in which bacterial clones emerge and diversify to form clonal complexes. An increasingly popular approach, called BURST, has been used as an alternative, but present implementations are unable to cope with very large data sets and offer crude graphical outputs. Here we present a new implementation of this algorithm, eBURST, which divides an MLST data set of any size into groups of related isolates and clonal complexes, predicts the founding (ancestral) genotype of each clonal complex, and computes the bootstrap support for the assignment. The most parsimonious patterns of descent of all isolates in each clonal complex from the predicted founder(s) are then displayed. The advantages of eBURST for exploring patterns of evolutionary descent are demonstrated with a number of examples, including the simple Spain(23F)-1 clonal complex of Streptococcus pneumoniae, "population snapshots" of the entire S. pneumoniae and Staphylococcus aureus MLST databases, and the more complicated clonal complexes observed for Campylobacter jejuni and Neisseria meningitidis.     

High-throughput multilocus sequence typing: bringing molecular typing to the next level

Multilocus sequence typing (MLST) is a widely used system for typing microorganisms by sequence analysis of housekeeping genes. The main advantage of MLST in comparison to other typing techniques is the unambiguity and transferability of sequence data. However, a main disadvantage is the high cost of DNA sequencing. Here we introduce a high-throughput MLST (HiMLST) method that employs next-generation sequencing (NGS) technology (Roche 454), to generate large quantities of high-quality MLST data at low costs. The HiMLST protocol consists of two steps. In the first step MLST target genes are amplified by PCR in multi-well plates. During this PCR the amplicons of each bacterial isolate are provided with a unique DNA barcode, the multiplex identifier (MID). In the second step all amplicons are pooled and sequenced in a single NGS-run. The MLST profile of each individual isolate can be retrieved easily using its unique MID. With HiMLST we have profiled 575 isolates of Legionella pneumophila, Staphylococcus aureus, Pseudomonas aeruginosa and Streptococcus pneumoniae in mixed species HiMLST experiments. In conclusion, the introduction of HiMLST paves the way for a broad employment of the MLST as a high-quality and cost-effective method for typing microbial species.     

mlstdbNet – distributed multi-locus sequence typing (MLST) databases

Background  

Multi-locus sequence typing (MLST) is a method of typing that facilitates the discrimination of microbial isolates by comparing the sequences of housekeeping gene fragments. The mlstdbNet software enables the implementation of distributed web-accessible MLST databases that can be linked widely over the Internet.     

Multilocus sequencing--a new method of genotyping bacteria and first results of its use

Comparative characterization (molecular typing) of isolates within a bacterial species is one of the major problems in microbiology and epidemiology. However, it is rather difficult to correlate data obtained in various laboratories, because traditional, including molecular, methods employed in typing pathogenic microorganisms cannot be standardized. In 1998, Maiden et al. proposed multilocus sequence typing (MLST); through which alleles of several housekeeping genes are directly assessed by nucleotide sequencing, each unique allele combination determining a sequence type of a strain. The advantages of this approach are that the culturing of pathogenic microorganisms is avoided, as their gene fragments are amplified directly from biological samples, and that the sequencing data are unambiguous, easy to standardize, and electronically portable. The latter makes it possible to generate an expandable global database for each species at an Internet site, in order to use it for the purposes of genotyping pathogenic bacteria (and other infectious agents). MLST protocols have been elaborated for Neisseria meningitidis, Streptococcus pneumoniae, and Helicobacter pylori; those for Streptococcus pyogenes, Staphylococcus aureus, and Haemophilus influenzae are now being developed. Basic principles and the first results of MLST have been reviewed, including data on the distribution and microevolution of N. meningitidis clones causing epidemic meningococcal infection, the relative recombination and mutation rates in the N. meningitidis genome, the identification of antibiotic-resistant S. pneumoniae clones causing severe generalized infection, the grouping of H. pylori isolates from various geographic regions, etc.     

Rapid pulsed-field gel electrophoresis protocol for subtyping of Streptococcus suis serotype 2

A rapid pulsed-field gel electrophoresis (PFGE) protocol for subtyping of Streptococcus suis serotype 2 was developed and evaluated using 27 clinical isolates from 22 epidemiologically unrelated patients. Results were matched against antibiogram, virulence genotyping and multi locus sequence typing (MLST). PFGE appeared to be the most discriminatory with numerical index of discrimination (D) equal to 0.87.     

儿童感染青霉素耐药肺炎链球菌的 多位点序列分型

目的对儿童感染的青霉素耐药肺炎链球菌进行多位点序列分型,了解厦门地区肺炎链球菌青霉素耐药菌株遗传背景。方法采用多位点序列分型法对2012年1月至2014年12月期间分离的60株青霉素耐药肺炎链球菌进行分子分型。结果 60株青霉素耐药肺炎链球菌MLST法共检出24个ST型,其中发现6个新的ST型,分别被命名为ST10004、ST10005、ST10006、ST10007、ST10008和ST10009。存在一个优势型别ST271,占31.7%(19/60),发现了4个克隆群和20种单一克隆,其中主要克隆群为国际流行耐药克隆群Taiwan19F-14,占41.7%(25/60)。结论本地区分离的儿童青霉素耐药肺炎链球菌主要以ST271型为主,属国际流行耐药克隆群Taiwan19F-14,是引起儿童呼吸道感染肺炎链球菌多重耐药的主要原因。     

Comparison of the Diversilab® system with multi-locus sequence typing and pulsed-field gel electrophoresis for the characterization of Streptococcus agalactiae invasive strains

Streptococcus agalactiae (or group B streptococcus; GBS) is a leading cause of neonatal morbidity and mortality in the developed countries. Several epidemiological typing tools have been developed for GBS to investigate the association between genotype and disease and to assess genetic variations within genogroups. This study compared the semi-automated repetitive sequence-based PCR Diversilab® system (DL) with MLST and pulsed field gel electrophoresis (PFGE) for determining the relatedness of invasive GBS strains. We analysed 179 unrelated GBS strains isolated from adult (n = 108) and neonatal (n = 71) invasive infections. By MLST, strains were resolved into 6 clonal complexes (CCs) including 23 sequence-types (STs), and 4 unique STs, whereas DL differentiated these isolates into 12 rep-PCR clusters (rPCs) and 9 unique rep-PCR types. The discriminatory power of both methods was similar, with Simpson's diversity indexes of 71.9% and 70.6%, respectively. However, their clustering concordance was low with Wallace concordance coefficients inferior to 0.4. PFGE was performed on 31 isolates representative of the most relevant DLrPCs clustered within the 3 major MLST CCs (CC-17, CC-23 and CC-1). As already observed with MLST, the concordance of DL with PFGE was low for all three CCs (Wallace coefficient < 0.5), PFGE being more discriminative than rep-PCR. In summary, this work suggests that DL is less appropriate than MLST or PFGE to study GBS population genetic diversity.     

Comparison of two DNA sequence-based typing schemes for the Fusarium solani Species Complex and proposal of a new consensus method

Multilocus sequence typing (MLST) is a widely used approach for differentiating microbial isolates presenting many advantages such as easy access through online databases and straightforward interpretation. For the Fusarium solani species complex (FSSC), three gene regions have been widely used to investigate phylogenetic relationships at the interspecific level (ITS-nuLSU, EF1a, RPB2) and a nomenclature system has been proposed for the different known haplotypes. More recently, a MLST scheme was proposed for this species complex based on the polymorphisms of five housekeeping genes (ACC, ICL, GDP, MDP, SOD). Here, we compare the phylogenetic resolution and sequence discriminatory powers of these two sets of loci on 50 epidemiologically unrelated FSSC strains. Although the widely used gene set offers better phylogenetic resolution, the newly developed gene set is slightly better at discriminating isolates using a MLST method. A consensus scheme of eight loci is proposed for typing FSSC strains combining the advantages of the two previous gene sets and offering the best typing efficiency.     

Comparative analysis of whole genome structure of <Emphasis Type="Italic">Streptococcus suis</Emphasis> using whole genome PCR scanning

An outbreak associated with Streptococcus suis infection in humans emerged in Sichuan province, China in 2005. The outbreak is atypical for the apparent large number of human cases, high fatality rate and geographical spread. To determine whether the bacterium has changed, we compared both human and animal isolates from the Sichuan outbreak with those collected previously within China and in other countries using whole genome PCR scanning (WGPScaning) comparative sequencing of several known virulence factor genes and multilocus sequence typing (MLST) analysis. WGPScanning analysis showed that all primer pairs yielded PCR products of the expected sizes in all four strains tested. The nucleotide sequences of all the detected virulence factor genes are identical in the four strains and MLST results showed that the four isolates studied and reference strain all belonged to the ST1 complex. No new genetic changes were found in the genome structure of the isolates from this Sichuan outbreak.     

Minim typing--a rapid and low cost MLST based typing tool for Klebsiella pneumoniae

Here we report a single nucleotide polymorphism (SNP) based genotyping method for Klebsiella pneumoniae utilising high-resolution melting (HRM) analysis of fragments within the multilocus sequence typing (MLST) loci. The approach is termed mini-MLST or Minim typing and it has previously been applied to Streptococcus pyogenes, Staphylococcus aureus and Enterococcus faecium. Six SNPs were derived from concatenated MLST sequences on the basis of maximisation of the Simpsons Index of Diversity (D). DNA fragments incorporating these SNPs and predicted to be suitable for HRM analysis were designed. Using the assumption that HRM alleles are defined by G+C content, Minim typing using six fragments was predicted to provide a D = 0.979 against known STs. The method was tested against 202 K. pneumoniae using a blinded approach in which the MLST analyses were performed after the HRM analyses. The HRM-based alleles were indeed in accordance with G+C content, and the Minim typing identified known STs and flagged new STs. The tonB MLST locus was determined to be very diverse, and the two Minim fragments located herein contribute greatly to the resolving power. However these fragments are refractory to amplification in a minority of isolates. Therefore, we assessed the performance of two additional formats: one using only the four fragments located outside the tonB gene (D = 0.929), and the other using HRM data from these four fragments in conjunction with sequencing of the tonB MLST fragment (D = 0.995). The HRM assays were developed on the Rotorgene 6000, and the method was shown to also be robust on the LightCycler 480, allowing a 384-well high through-put format. The assay provides rapid, robust and low-cost typing with fully portable results that can directly be related to current MLST data. Minim typing in combination with molecular screening for antibiotic resistance markers can be a powerful surveillance tool kit.     

The relative contributions of recombination and mutation to the divergence of clones of Neisseria meningitidis.

Multilocus sequence typing (MLST) is a recently developed nucleotide sequence-based method for the definitive assignment of isolates within bacterial populations to specific clones. MLST uses the same principles as multilocus enzyme electrophoresis and provides data that can be used to investigate aspects of the population genetics and evolution of bacterial species. We used an MLST data set consisting of the sequences of approximately 450-bp fragments from seven housekeeping loci from a large strain collection of Neisseria meningitidis to estimate the relative impact of recombination compared with point mutation in the diversification of N. meningitidis clonal complexes. 126 meningococcal isolates were assigned to 10 clonal complexes, 9 of which contained minor clonal variants. The allelic variation within each complex was classified as a recombinational exchange or a putative point mutation through a comparison of the sequences of each variant allele with that of the allele typically found in the clonal complex. The nine clonal complexes contained a total of 23 allelic variants, and analysis of the sequences of these variant alleles revealed that a single nucleotide site in a meningococcal housekeeping gene is at least 80-fold more likely to change as a result of recombination than as a result of mutation. This value is estimated to be 10-50-fold for Escherichia coli and approximately 50-fold for Streptococcus pneumoniae.     

Estimating recombinational parameters in Streptococcus pneumoniae from multilocus sequence typing data

Multilocus sequence typing (MLST) is a highly discriminatory molecular typing method that defines isolates of bacterial pathogens using the sequences of approximately 450-bp internal fragments of seven housekeeping genes. This technique has been applied to 575 isolates of Streptococcus pneumoniae and identifies a number of discrete clonal complexes. These clonal complexes are typically represented by a single group of isolates sharing identical alleles at all seven loci, plus single-locus variants that differ from this group at only one out of the seven loci. As MLST is highly discriminatory, the members of each clonal complex can be assumed to have a recent common ancestor, and the molecular events that give rise to the single-locus variants can be used to estimate the relative contributions of recombination and mutation to clonal divergence. By comparing the sequences of the variant alleles within each clonal complex with the allele typically found within that clonal complex, we estimate that recombination has generated new alleles at a frequency approximately 10-fold higher than mutation, and that a single nucleotide site is approximately 50 times more likely to change through recombination than mutation. We also demonstrate how to estimate the average length of recombinational replacements from MLST data.     

The utility of high-resolution melting analysis of SNP nucleated PCR amplicons--an MLST based Staphylococcus aureus typing scheme

High resolution melting (HRM) analysis is gaining prominence as a method for discriminating DNA sequence variants. Its advantage is that it is performed in a real-time PCR device, and the PCR amplification and HRM analysis are closed tube, and effectively single step. We have developed an HRM-based method for Staphylococcus aureus genotyping. Eight single nucleotide polymorphisms (SNPs) were derived from the S. aureus multi-locus sequence typing (MLST) database on the basis of maximized Simpson's Index of Diversity. Only G?A, G?T, C?A, C?T SNPs were considered for inclusion, to facilitate allele discrimination by HRM. In silico experiments revealed that DNA fragments incorporating the SNPs give much higher resolving power than randomly selected fragments. It was shown that the predicted optimum fragment size for HRM analysis was 200 bp, and that other SNPs within the fragments contribute to the resolving power. Six DNA fragments ranging from 83 bp to 219 bp, incorporating the resolution optimized SNPs were designed. HRM analysis of these fragments using 94 diverse S. aureus isolates of known sequence type or clonal complex (CC) revealed that sequence variants are resolved largely in accordance with G+C content. A combination of experimental results and in silico prediction indicates that HRM analysis resolves S. aureus into 268 "melt types" (MelTs), and provides a Simpson's Index of Diversity of 0.978 with respect to MLST. There is a high concordance between HRM analysis and the MLST defined CCs. We have generated a Microsoft Excel key which facilitates data interpretation and translation between MelT and MLST data. The potential of this approach for genotyping other bacterial pathogens was investigated using a computerized approach to estimate the densities of SNPs with unlinked allelic states. The MLST databases for all species tested contained abundant unlinked SNPs, thus suggesting that high resolving power is not dependent upon large numbers of SNPs.     

Comparative analysis of whole genome structure of Streptococcus suis using whole genome PCR scanning

An outbreak associated with Streptococcus suis infection in humans emerged in Sichuan province, China in 2005. The outbreak is atypical for the apparent large number of human cases, high fatality rate and geographical spread. To determine whether the bacterium has changed, we compared both human and animal isolates from the Sichuan outbreak with those collected previously within China and in other countries using whole genome PCR scanning (WGPScaning) comparative sequencing of several known virulence factor genes and multilocus sequence typing (MLST) analysis. WGPScanning analysis showed that all primer pairs yielded PCR products of the expected sizes in all four strains tested. The nucleotide sequences of all the detected virulence factor genes are identical in the four strains and MLST results showed that the four isolates studied and reference strain all belonged to the ST1 complex. No new genetic changes were found in the genome structure of the isolates from this Sichuan outbreak. Contributed equally to this work Supported by the National Key Technologies Research and Development Program (Grant No. 2005BA711A09) from the Ministry of Science and Technology of China     

How clonal is Staphylococcus aureus?

Staphylococcus aureus is an important human pathogen and represents a growing public health burden owing to the emergence and spread of antibiotic-resistant clones, particularly within the hospital environment. Despite this, basic questions about the evolution and population biology of the species, particularly with regard to the extent and impact of homologous recombination, remain unanswered. We address these issues through an analysis of sequence data obtained from the characterization by multilocus sequence typing (MLST) of 334 isolates of S. aureus, recovered from a well-defined population, over a limited time span. We find no significant differences in the distribution of multilocus genotypes between strains isolated from carriers and those from patients with invasive disease; there is, therefore, no evidence from MLST data, which index variation within the stable "core" genome, for the existence of hypervirulent clones of this pathogen. Examination of the sequence changes at MLST loci during clonal diversification shows that point mutations give rise to new alleles at least 15-fold more frequently than does recombination. This contrasts with the naturally transformable species Neisseria meningitidis and Streptococcus pneumoniae, in which alleles change between 5- and 10-fold more frequently by recombination than by mutation. However, phylogenetic analysis suggests that homologous recombination does contribute toward the evolution of this species over the long term. Finally, we note a striking excess of nonsynonymous substitutions in comparisons between isolates belonging to the same clonal complex compared to isolates belonging to different clonal complexes, suggesting that the removal of deleterious mutations by purifying selection may be relatively slow.     

Sequence Analysis of 96 Genomic Regions Identifies Distinct Evolutionary Lineages within CC156, the Largest Streptococcus pneumoniae Clonal Complex in the MLST Database

Multi-Locus Sequence Typing (MLST) of Streptococcus pneumoniae is based on the sequence of seven housekeeping gene fragments. The analysis of MLST allelic profiles by eBURST allows the grouping of genetically related strains into Clonal Complexes (CCs) including those genotypes with a common descent from a predicted ancestor. However, the increasing use of MLST to characterize S. pneumoniae strains has led to the identification of a large number of new Sequence Types (STs) causing the merger of formerly distinct lineages into larger CCs. An example of this is the CC156, displaying a high level of complexity and including strains with allelic profiles differing in all seven of the MLST loci, capsular type and the presence of the Pilus Islet-1 (PI-1). Detailed analysis of the CC156 indicates that the identification of new STs, such as ST4945, induced the merging of formerly distinct clonal complexes. In order to discriminate the strain diversity within CC156, a recently developed typing schema, 96-MLST, was used to analyse 66 strains representative of 41 different STs. Analysis of allelic profiles by hierarchical clustering and a minimum spanning tree identified ten genetically distinct evolutionary lineages. Similar results were obtained by phylogenetic analysis on the concatenated sequences with different methods. The identified lineages are homogenous in capsular type and PI-1 presence. ST4945 strains were unequivocally assigned to one of the lineages. In conclusion, the identification of new STs through an exhaustive analysis of pneumococcal strains from various laboratories has highlighted that potentially unrelated subgroups can be grouped into a single CC by eBURST. The analysis of additional loci, such as those included in the 96-MLST schema, will be necessary to accurately discriminate the clonal evolution of the pneumococcal population.