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.     

Investigating Burkholderia cepacia complex populations recovered from Italian maize rhizosphere by multilocus sequence typing

The Burkholderia cepacia complex (BCC) comprises at least nine closely related species of abundant environmental microorganisms. Some of these species are highly spread in the rhizosphere of several crop plants, particularly of maize; additionally, as opportunistic pathogens, strains of the BCC are capable of colonizing humans. We have developed and validated a multilocus sequence typing (MLST) scheme for the BCC. Although widely applied to understand the epidemiology of bacterial pathogens, MLST has seen limited application to the population analysis of species residing in the natural environment; we describe its novel application to BCC populations within maize rhizospheres. 115 BCC isolates were recovered from the roots of different maize cultivars from three different Italian regions over a 9-year period (1994-2002). A total of 44 sequence types (STs) were found of which 41 were novel when compared with existing MLST data which encompassed a global database of 1000 clinical and environmental strains representing nearly 400 STs. In this study of rhizosphere isolates approximately 2.5 isolates per ST was found, comparable to that found for the whole BCC population. Multilocus sequence typing also resolved inaccuracies associated with previous identification of the maize isolates based on recA gene restriction fragment length polymorphims and species-specific polymerase chain reaction. The 115 maize isolates comprised the following BCC species groups, B. ambifaria (39%), BCC6 (29%), BCC5 (10%), B. pyrrocinia (8%), B. cenocepacia IIIB (7%) and B. cepacia (6%), with BCC5 and BCC6 potentially constituting novel species groups within the complex. Closely related clonal complexes of strains were identified within B. cepacia, B. cenocepacia IIIB, BCC5 and BCC6, with one of the BCC5 clonal complexes being distributed across all three sampling sites. Overall, our analysis demonstrates that the maize rhizosphere harbours a massive diversity of novel BCC STs, so that their addition to our global MLST database increased the ST diversity by 10%.     

Evidence of land‐sea transfer of the zoonotic pathogen Campylobacter to a wildlife marine sentinel species

Environmental pollution often accompanies the expansion and urbanization of human populations where sewage and wastewaters commonly have an impact on the marine environments. Here, we explored the potential for faecal bacterial pathogens, of anthropic origin, to spread to marine wildlife in coastal areas. The common zoonotic bacterium Campylobacter was isolated from grey seals (Halichoerus grypus), an important sentinel species for environmental pollution, and compared to isolates from wild birds, agricultural sources and clinical samples to characterize possible transmission routes. Campylobacter jejuni was present in half of all grey seal pups sampled (24/50 dead and 46/90 live pups) in the breeding colony on the Isle of May (Scotland), where it was frequently associated with histological evidence of disease. Returning yearling animals (19/19) were negative for C. jejuni suggesting clearance of infection while away from the localized colony infection source. The genomes of 90 isolates from seals were sequenced and characterized using a whole‐genome multilocus sequence typing (MLST) approach and compared to 192 published genomes from multiple sources using population genetic approaches and a probabilistic genetic attribution model to infer the source of infection from MLST data. The strong genotype‐host association has enabled the application of source attribution models in epidemiological studies of human campylobacteriosis, and here assignment analyses consistently grouped seal isolates with those from human clinical samples. These findings are consistent with either a common infection source or direct transmission of human campylobacter to grey seals, raising concerns about the spread of human pathogens to wildlife marine sentinel species in coastal areas.     

Multilocus sequence typing reveals high genetic diversity and epidemic population structure for the fish pathogen Yersinia ruckeri

Yersinia ruckeri is the causative agent of enteric redmouth in fish and one of the major bacterial pathogens causing losses in salmonid aquaculture. Previously typing methods, including restriction enzyme analysis, pulsed-field gel electrophoresis and multilocus enzyme electrophoresis (MLEE) have indicated a clonal population structure. In this work, we describe a multilocus sequence typing (MLST) scheme for Y.ruckeri based on the internal fragment sequence of six housekeeping genes. This MLST scheme was applied to 103 Y.ruckeri strains from diverse geographic areas and hosts as well as environmental sources. Sequences obtained from this work were deposited and are available in a public database (http://publmst.org/yruckeri/). Thirty different sequence types (ST) were identified, 21 of which were represented by a single isolate, evidencing high genetic diversity. ST2 comprised more than one-third of the isolates and was most frequently observed among isolates from trout. Two major clonal complexes (CC) were identified by eBURST analysis showing a common evolutionary origin for 94 isolates forming 21 STs into CC1 and for 6 isolates of 6 STs in the CC2. It was also possible to associate some unique ST with isolates from recent outbreaks in vaccinated salmonid fish.     

Multilocus Sequence Typing (MLST) for Characterization of Enterobacter cloacae

Enterobacter cloacae is an important emerging pathogen, which sometime causes respiratory infection, surgical site infection, urinary infection, sepsis, and outbreaks at neonatal units. We have developed a multilocus sequence typing (MLST) scheme utilizing seven housekeeping genes and evaluated the performance in 101 clinical isolates. The MLST scheme yielded 83 sequence types (ST) including 78 novel STs found in the clinical isolates. These findings supported the robustness of the MLST scheme developed in this study.     

Bacterial population genetics, evolution and epidemiology.

Asexual bacterial populations inevitably consist of an assemblage of distinct clonal lineages. However, bacterial populations are not entirely asexual since recombinational exchanges occur, mobilizing small genome segments among lineages and species. The relative contribution of recombination, as opposed to de novo mutation, in the generation of new bacterial genotypes varies among bacterial populations and, as this contribution increases, the clonality of a given population decreases. In consequence, a spectrum of possible population structures exists, with few bacterial species occupying the extremes of highly clonal and completely non-clonal, most containing both clonal and non-clonal elements. The analysis of collections of bacterial isolates, which accurately represent the natural population, by nucleotide sequence determination of multiple housekeeping loci provides data that can be used both to investigate the population structure of bacterial pathogens and for the molecular characterization of bacterial isolates. Understanding the population structure of a given pathogen is important since it impacts on the questions that can be addressed by, and the methods and samples required for, effective molecular epidemiological studies.     

Comparison of a semi-automated rep-PCR system and multilocus sequence typing for differentiation of Salmonella enterica isolates

The accurate sub-typing of Salmonella enterica isolates is essential for epidemiological investigations and surveillance of Salmonella infections. Salmonella isolates are currently identified using the Kauffman-White serotyping scheme. Multilocus sequence typing (MLST) schemes have been developed for the major bacterial pathogens including Salmonella and have assisted in understanding the molecular epidemiology and population biology of these organisms. Recently, the DiversiLab rep-PCR system has been developed using micro-fluidic chips to provide standardized, semi-automated fingerprinting for pathogens including S. enterica. In the current study, 71 isolates of S. enterica, representing 21 serovars, were analyzed using MLST and the DiversiLab rep-PCR system. MLST was able to identify 31 sequence types (STs), while the DiversiLab system revealed 38 DiversiLab types (DTs). The rep-PCR distinguished isolates of different serovars and showed greater discriminatory power (0.95) than MLST typing (0.89). Rep-PCR exhibited 92% concordance with MLST and 90% with serotyping, while the concordance level of MLST typing with serotyping was 96%, representing a strong association. Comparison of rep-PCR profiles with those held in an online library database led to the accurate prediction of serovar in 63% of cases and resulted in inaccurate predictions for 10% of profiles. MLST and the rep-PCR system may provide useful additional informative techniques for the molecular identification of S. enterica. We conclude that the DiversiLab rep-PCR system may provide a rapid (less than 4 h) and standardized method for sub-typing isolates of S. enterica.     

Multi-locus sequence typing: a tool for global epidemiology

The characterization of pathogenic isolates plays a pivotal role in the epidemiology of infectious diseases, generating the information necessary for identifying, tracking, and intervening against disease outbreaks. In 1998 multi-locus sequence typing (MLST) was proposed as a nucleotide sequence-based approach that could be applied to many bacterial pathogens. It combined developments in high-throughput sequencing and bioinformatics with established population genetics techniques to provide a portable, reproducible, and scalable typing system that reflected the population and evolutionary biology of bacterial pathogens. MLST schemes have been developed for a variety of procaryotic and eucaryotic pathogens and the data generated have contributed to both epidemiological surveillance and fundamental studies of pathogen biology.     

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.     

江苏水产养殖区拟态弧菌种群结构及遗传多样性

【背景】拟态弧菌(Vibrio mimicus)是一种常见的革兰氏阴性病原菌,广泛分布于水环境和水生动物体内,可导致多种水产动物和人类感染。多位点序列分型(multilocus sequence typing, MLST)已被应用于多种病原菌的分子分型,其通过分析不同菌株之间的遗传关系,监测细菌传播的时间和地理分布,确定感染和传播途径,但目前未见有关拟态弧菌MLST的报道。【目的】开发一种基于MLST的拟态弧菌分型方法,并用于江苏水产养殖区拟态弧菌的种群结构和遗传进化分析,为拟态弧菌感染所引起的疾病防治提供理论基础。【方法】选择拟态弧菌的7个管家基因dnaE、gyrB、mdh、recA、rpoD、pntA和pyrH作为靶点,对江苏水产养殖区分离的155株拟态弧菌进行PCR扩增和测序。将测序结果分配等位基因,制作等位基因谱,分配不同的序列类型(sequence type, ST),利用软件goeBURST-1.2.1和MEGA-X对分配的ST型进行克隆复合体和遗传进化树聚类分析;此外,利用Kirby-Bauer圆盘扩散法测试155株拟态弧菌的药敏特性。【结果】155株拟态弧菌被分为56个STs,其中ST11占比最高;在双位点变异(double locus variants, DLV)水平分析发现56个STs分为3个克隆复合体和3个单体;系统发育树显示,56个STs被分为3个集群(cluster I、cluster II、cluster III)。药敏结果显示,155株拟态弧菌对红霉素类抗生素的耐药性最高(88.39%, 137/155),对氯霉素类抗生素敏感性最高(91.61%, 142/155)。【结论】本研究建立的MLST方法具有良好的分辨率,可作为拟态弧菌系统发育和未来流行病学调查有用的分子分型工具。根据抗生素耐药谱结果,提示在养殖过程中可选用氟苯尼考等国家批准使用的专用抗菌药对拟态弧菌进行防治。     

Genetic diversity of Staphylococcus equorum isolates from Saeu-jeotgal evaluated by multilocus sequence typing

Staphylococcus equorum, the predominant bacterial species detected in Saeu-jeotgal, a Korean high-salt fermented seafood, is a candidate starter bacterium for Saeu-jeotgal fermentation. A multilocus sequence typing (MLST) scheme was developed to evaluate the genetic diversity and background of S. equorum strains isolated from Saeu-jeotgal. A total of 135 strains, including 117 isolates from Saeu-jeotgal, and others from Myeolchi-jeotgal, sausage, cheese and horse skin, were subjected to MLST, and the internal fragments of seven housekeeping genes, aroE, dnaJ, glpF, gmk, hsp60, mutS, and pta, were compared. This MLST scheme produced 45 sequence types (STs) and the eBURST algorithm clustered the STs into nine clonal groups and seven singletons. Clonal group 1, the major group, consisted of 30 isolates from cheese, Saeu-jeotgal and sausages, which were classified into 12 STs. The predominant ST, ST26, comprised 25 isolates and presented as a singleton. Most of the isolates from Myeolchi-jeotgal and sausages clustered on two different branches of a phylogenetic tree generated with a cluster analysis using the maximum likelihood algorithm. This MLST scheme established the genetic backgrounds of S. equorum strains isolated from different types of food. Among the housekeeping genes used for MLST, gmk had the fewest allele types and fairly low sequence identities (74.0–90.0 %) within the Staphylococcus species. Therefore, sequence analyses of the gmk gene and 16S rRNA gene can be used for the accurate and rapid identification of S. equorum.     

Multilocus sequence typing for global surveillance of meningococcal disease

The global surveillance of bacterial pathogens is particularly important for bacteria with diverse and dynamic populations that cause periodic epidemics or pandemics. The isolate characterization methods employed for surveillance should: (1) generate unambiguous data; (2) be readily implemented in a variety of scenarios and be reproducible among laboratories; (3) be scalable and preferably available in a high throughput format; and (4) be cost effective. Multilocus sequence typing (MLST) was designed to meet these criteria and has been implemented effectively for a wide range of microorganisms. The 'Impact of meningococcal epidemiology and population biology on public health in Europe (EU-MenNet)' project had amongst its objectives: (1) to disseminate meningococcal MLST and sequence-based typing throughout Europe by establishing a centre for training and data generation, and (2) to produce a comprehensive Europe-wide picture of meningococcal disease epidemiology for the first time. Data produced from the project have shown the distribution of a relatively small number of STs, clonal complexes and PorA types that account for a large proportion of the disease-associated isolates in Europe. The project demonstrates how molecular typing can be combined with epidemiological data via the Internet for global disease surveillance.     

Diversification and Distribution of Ruminant Chlamydia abortus Clones Assessed by MLST and MLVA

Chlamydia abortus, an obligate intracellular bacterium, is the most common infectious cause of abortion in small ruminants worldwide and has zoonotic potential. We applied multilocus sequence typing (MLST) together with multiple-locus variable-number tandem repeat analysis (MLVA) to genotype 94 ruminant C. abortus strains, field isolates and samples collected from 1950 to 2011 in diverse geographic locations, with the aim of delineating C. abortus lineages and clones. MLST revealed the previously identified sequence types (STs) ST19, ST25, ST29 and ST30, plus ST86, a recently-assigned type on the Chlamydiales MLST website and ST87, a novel type harbouring the hemN_21 allele, whereas MLVA recognized seven types (MT1 to MT7). Minimum-spanning-tree analysis suggested that all STs but one (ST30) belonged to a single clonal complex, possibly reflecting the short evolutionary timescale over which the predicted ancestor (ST19) has diversified into three single-locus variants (ST86, ST87 and ST29) and further, through ST86 diversification, into one double-locus variant (ST25). ST descendants have probably arisen through a point mutation evolution mode. Interestingly, MLVA showed that in the ST19 population there was a greater genetic diversity than in other STs, most of which exhibited the same MT over time and geographical distribution. However, the evolutionary pathways of C. abortus STs seem to be diverse across geographic distances with individual STs restricted to particular geographic locations. The ST30 singleton clone displaying geographic specificity and represented by the Greek strains LLG and POS was effectively distinguished from the clonal complex lineage, supporting the notion that possibly two separate host adaptations and hence independent bottlenecks of C. abortus have occurred through time. The combination of MLST and MLVA assays provides an additional level of C. abortus discrimination and may prove useful for the investigation and surveillance of emergent C. abortus clonal populations.     

Population genetics of a transformable bacterium: The influence of horizontal genetic exchange on the biology of Neisseria meningitidis

Abstract Information on the biochemistry and genetics of bacterial species, usually obtained by the study of single isolates, is enhanced by studies of populations of bacteria. Recent advances in molecular technology, particularly polymerase chain reaction-based nucleotide sequence analysis, provide powerful for the study of population genetics. Data obtained by such techniques indicate that, while some bacterial species have a clonal population structure, others are non-clonal or panmictic. Clonal populations are a consequence of asexual reproduction by binary fission; panmictic population structures results from 'horizontal' exchange of genetic material between clones. A consequence of horizontal genetic exchange is mosaic gene structures, recognisable by comparisons of nucleotide sequences. In transformable bacteria, for example the human pathogen Neisseria meningitidis , several different genes, including the gene encoding the class 1 outer membrane protein, a major surface antigen, are mosaics. This genetic process has implications both for vaccine design and in the interpretation of epidemiological data.     

Multilocus sequence typing as an approach for population analysis of Medicago-nodulating rhizobia

Multilocus sequence typing (MLST), a sequence-based method to characterize bacterial genomes, was used to examine the genetic structure in a large collection of Medicago-nodulating rhizobial strains. This is the first study where MLST has been applied in conjunction with eBURST analysis to determine the population genetic structure of nonpathogenic bacteria recovered from the soil environment. Sequence variation was determined in 10 chromosomal loci of 231 strains that predominantly originated from southwest Asia. Genetic diversity for each locus ranged from 0.351 to 0.819, and the strains examined were allocated to 91 different allelic profiles or sequence types (STs). The genus Medicago is nodulated by at least two groups of rhizobia with divergent chromosomes that have been classified as Sinorhizobium meliloti and Sinorhizobium medicae. Evidence was obtained that the degree of genetic exchange among the chromosomes across these groups is limited. The symbiosis with Medicago polymorpha of nine strains placed in one of these groups, previously identified as S. medicae, ranged from ineffective to fully effective, indicating that there was no strong relationship between symbiotic phenotype and chromosomal genotype.     

Genetic diversity and structure of Rhizobium leguminosarum populations associated with clover plants are influenced by local environmental variables

The identification and conservation of indigenous rhizobia associated with legume plants and their application as biofertilizers is becoming an agricultural worldwide priority. However, little is known about the genetic diversity and phylogeny of rhizobia in Romania. In the present study, the genetic diversity and population composition of Rhizobium leguminosarum symbiovar trifolii isolates from 12 clover plants populations located across two regions in Romania were analyzed. Red clover isolates were phenotypically evaluated and genotyped by sequencing 16S rRNA gene, 16S-23S intergenic spacer, three chromosomal genes (atpD, glnII and recA) and two plasmid genes (nifH and nodA). Multilocus sequence typing (MLST) analysis revealed that red clover plants are nodulated by a wide genetic diversity of R. leguminosarum symbiovar trifolii sequence types (STs), highly similar to the ones previously found in white clover. Rhizobial genetic variation was found mainly within the two clover populations for both chromosomal and plasmid types. Many STs appear to be unique for this region and the genetic composition of rhizobia differs significantly among the clover populations. Furthermore, our results showed that both soil pH and altitude contributed to plasmid sequence type composition while differences in chromosomal composition were affected by the altitude and were strongly correlated with distance.     

Using multilocus sequence data to define the pneumococcus

We investigated the genetic relationships between serotypeable pneumococci and nonserotypeable presumptive pneumococci using multilocus sequence typing (MLST) and partial sequencing of the pneumolysin gene (ply). Among 121 nonserotypeable presumptive pneumococci from Finland, we identified isolates of three classes: those with sequence types (STs) identical to those of serotypeable pneumococci, suggesting authentic pneumococci in which capsular expression had been downregulated or lost; isolates that clustered among serotypeable pneumococci on a tree based on the concatenated sequences of the MLST loci but which had STs that differed from those of serotypeable pneumococci in the MLST database; and a more diverse collection of isolates that did not cluster with serotypeable pneumococci. The latter isolates typically had sequences at all seven MLST loci that were 5 to 10% divergent from those of authentic pneumococci and also had distinct and divergent ply alleles. These isolates are proposed to be distinct from pneumococci but cannot be resolved from them by optochin susceptibility, bile solubility, or the presence of the ply gene. Complete resolution of pneumococci from the related but distinct population is problematic, as recombination between them was evident, and a few isolates of each population possessed alleles at one or occasionally more MLST loci from the other population. However, a tree based on the concatenated sequences of the MLST loci in most cases unambiguously distinguished whether a nonserotypeable isolate was or was not a pneumococcus, and the sequence of the ply gene fragment was found to be useful to resolve difficult cases.     

Investigation of genotypes of Borrelia burgdorferi in Ixodes scapularis ticks collected during surveillance in Canada

The genetic diversity of Borrelia burgdorferi sensu stricto, the agent of Lyme disease in North America, has consequences for the performance of serological diagnostic tests and disease severity. To investigate B. burgdorferi diversity in Canada, where Lyme disease is emerging, bacterial DNA in 309 infected adult Ixodes scapularis ticks collected in surveillance was characterized by multilocus sequence typing (MLST) and analysis of outer surface protein C gene (ospC) alleles. Six ticks carried Borrelia miyamotoi, and one tick carried the novel species Borrelia kurtenbachii. 142 ticks carried B. burgdorferi sequence types (STs) previously described from the United States. Fifty-eight ticks carried B. burgdorferi of 1 of 19 novel or undescribed STs, which were single-, double-, or triple-locus variants of STs first described in the United States. Clonal complexes with founder STs from the United States were identified. Seventeen ospC alleles were identified in 309 B. burgdorferi-infected ticks. Positive and negative associations in the occurrence of different alleles in the same tick supported a hypothesis of multiple-niche polymorphism for B. burgdorferi in North America. Geographic analysis of STs and ospC alleles were consistent with south-to-north dispersion of infected ticks from U.S. sources on migratory birds. These observations suggest that the genetic diversity of B. burgdorferi in eastern and central Canada corresponds to that in the United States, but there was evidence for founder events skewing the diversity in emerging tick populations. Further studies are needed to investigate the significance of these observations for the performance of diagnostic tests and clinical presentation of Lyme disease in Canada.