Insights to Genetic Characterization Tools for Epidemiological Tracking of Francisella tularensis in Sweden

Tularaemia, caused by the bacterium Francisella tularensis, is endemic in Sweden and is poorly understood. The aim of this study was to evaluate the effectiveness of three different genetic typing systems to link a genetic type to the source and place of tularemia infection in Sweden. Canonical single nucleotide polymorphisms (canSNPs), MLVA including five variable number of tandem repeat loci and PmeI-PFGE were tested on 127 F. tularensis positive specimens collected from Swedish case-patients. All three typing methods identified two major genetic groups with near-perfect agreement. Higher genetic resolution was obtained with canSNP and MLVA compared to PFGE; F. tularensis samples were first assigned into ten phylogroups based on canSNPs followed by 33 unique MLVA types. Phylogroups were geographically analysed to reveal complex phylogeographic patterns in Sweden. The extensive phylogenetic diversity found within individual counties posed a challenge to linking specific genetic types with specific geographic locations. Despite this, a single phylogroup (B.22), defined by a SNP marker specific to a lone Swedish sequenced strain, did link genetic type with a likely geographic place. This result suggests that SNP markers, highly specific to a particular reference genome, may be found most frequently among samples recovered from the same location where the reference genome originated. This insight compels us to consider whole-genome sequencing (WGS) as the appropriate tool for effectively linking specific genetic type to geography. Comparing the WGS of an unknown sample to WGS databases of archived Swedish strains maximizes the likelihood of revealing those rare geographically informative SNPs.     

Effect of a glucose impulse on the CcpA regulon in Staphylococcus aureus

Background

A low genetic diversity in Francisella tularensis has been documented. Current DNA based genotyping methods for typing F. tularensis offer a limited and varying degree of subspecies, clade and strain level discrimination power. Whole genome sequencing is the most accurate and reliable method to identify, type and determine phylogenetic relationships among strains of a species. However, lower cost typing schemes are necessary in order to enable typing of hundreds or even thousands of isolates.

Results

We have generated a high-resolution phylogenetic tree from 40 Francisella isolates, including 13 F. tularensis subspecies holarctica (type B) strains, 26 F. tularensis subsp. tularensis (type A) strains and a single F. novicida strain. The tree was generated from global multi-strain single nucleotide polymorphism (SNP) data collected using a set of six Affymetrix GeneChip^® resequencing arrays with the non-repetitive portion of LVS (type B) as the reference sequence complemented with unique sequences of SCHU S4 (type A). Global SNP based phylogenetic clustering was able to resolve all non-related strains. The phylogenetic tree was used to guide the selection of informative SNPs specific to major nodes in the tree for development of a genotyping assay for identification of F. tularensis subspecies and clades. We designed and validated an assay that uses these SNPs to accurately genotype 39 additional F. tularensis strains as type A (A1, A2, A1a or A1b) or type B (B1 or B2).

Conclusion

Whole-genome SNP based clustering was shown to accurately identify SNPs for differentiation of F. tularensis subspecies and clades, emphasizing the potential power and utility of this methodology for selecting SNPs for typing of F. tularensis to the strain level. Additionally, whole genome sequence based SNP information gained from a representative population of strains may be used to perform evolutionary or phylogenetic comparisons of strains, or selection of unique strains for whole-genome sequencing projects.     

Lab on a chip genotyping for Brucella spp. based on 15-loci multi locus VNTR analysis

Background  

Brucellosis is an important zoonosis caused by the genus Brucella. In addition Brucella represents potential biological warfare agents due to the high contagious rates for humans and animals. Therefore, the strain typing epidemiological tool may be crucial for tracing back source of infection in outbreaks and discriminating naturally occurring outbreaks versus bioterroristic event. A Multiple Locus Variable-number tandem repeats (VNTR) Analysis (MLVA) assay based on 15 polymorphic markers was previously described. The obtained MLVA band profiles may be resolved by techniques ranging from low cost manual agarose gels to the more expensive capillary electrophoresis sequencing. In this paper a rapid, accurate and reproducible system, based on the Lab on a chip technology was set up for Brucella spp. genotyping.     

Rapid and high resolution genotyping of all Escherichia coli serotypes using 10 genomic repeat-containing loci

Our laboratory has previously published two multiple-locus variable-number tandem-repeats analysis (MLVA) methods for rapid genotyping of Escherichia coli (E. coli), which are now in routine use for surveillance and outbreak detection. The first assay developed was specific for E. coli O157:H7; however this assay was not suitable for genotyping other E. coli serotypes. A new generic MLVA-assay was then developed with the capability of genotyping all E. coli serotypes. This generic E. coli MLVA (GECM7) was based on polymorphism in seven variable number of tandem repeats (VNTR) loci. GECM7 worked well with the majority of E. coli serotypes; however we wanted to increase the resolution for this method based in part of comparison with PFGE typing of E. coli O26:H11, where PFGE appeared to display higher resolution. The GECM7 method was improved by adding three new repeat-loci to a total of ten (GECM10), and a considerable increase in resolution was observed (from 296 to 507 genotypes on the same set of strains).     

Pathogenic potential and horizontal gene transfer in ovine gastrointestinal Escherichia coli

Aims: To demonstrate that a thorough characterization and virulotyping of Escherichia coli strains isolated from sheep over time leads to new insights into ovine E. coli potentially becoming human pathogens through horizontal gene transfer. Methods and Results: One hundred and fifty E. coli isolates from two sheep, sampled over 3 weeks, were characterized by serotyping, virulotyping, genotyping using multiple locus variable number tandem repeats analysis (MLVA) and susceptibility to phage infection in vitro. The 35 MLVA profiles and the serotype and virulotypes of the strains were closely associated. Many MLVA profiles differed in one locus independent of serotypes. Escherichia coli isolates of the same serotype or virulotype had identical or very similar MLVA profiles. No transductants that incorporated the bacteriophages were found in vivo, but six E. coli isolates were susceptible to the phage infection in vitro. Changes in MLVA profiles were seen after acquisition of Stx phages in vitro only. Conclusions: The sheep carried Stx phage susceptible E. coli that possessed virulence markers associated with human pathogenicity. Changes in bacterial genomes by phage transfer may complicate outbreak source investigations. Serotype has to be taken into account when evaluating strain relationships by MLVA. Significance and Impact of the Study: Sheep carry E. coli that encode for virulence markers and belong to serogroups known to be human pathogens. In addition, a selection of isolates was found to be susceptible to horizontal transfer of Shiga toxin genes by means of bacteriophages in vitro, and the transfer resulted in a discernible change of the MLVA patterns of E. coli.     

In Vitro vs In Silico Detected SNPs for the Development of a Genotyping Array: What Can We Learn from a Non-Model Species?

Background

There is considerable interest in the high-throughput discovery and genotyping of single nucleotide polymorphisms (SNPs) to accelerate genetic mapping and enable association studies. This study provides an assessment of EST-derived and resequencing-derived SNP quality in maritime pine (Pinus pinaster Ait.), a conifer characterized by a huge genome size (∼23.8 Gb/C).

Methodology/Principal Findings

A 384-SNPs GoldenGate genotyping array was built from i/ 184 SNPs originally detected in a set of 40 re-sequenced candidate genes (in vitro SNPs), chosen on the basis of functionality scores, presence of neighboring polymorphisms, minor allele frequencies and linkage disequilibrium and ii/ 200 SNPs screened from ESTs (in silico SNPs) selected based on the number of ESTs used for SNP detection, the SNP minor allele frequency and the quality of SNP flanking sequences. The global success rate of the assay was 66.9%, and a conversion rate (considering only polymorphic SNPs) of 51% was achieved. In vitro SNPs showed significantly higher genotyping-success and conversion rates than in silico SNPs (+11.5% and +18.5%, respectively). The reproducibility was 100%, and the genotyping error rate very low (0.54%, dropping down to 0.06% when removing four SNPs showing elevated error rates).

Conclusions/Significance

This study demonstrates that ESTs provide a resource for SNP identification in non-model species, which do not require any additional bench work and little bio-informatics analysis. However, the time and cost benefits of in silico SNPs are counterbalanced by a lower conversion rate than in vitro SNPs. This drawback is acceptable for population-based experiments, but could be dramatic in experiments involving samples from narrow genetic backgrounds. In addition, we showed that both the visual inspection of genotyping clusters and the estimation of a per SNP error rate should help identify markers that are not suitable to the GoldenGate technology in species characterized by a large and complex genome.     

Analysis of Salmonella enterica Serovar Typhimurium Variable-Number Tandem-Repeat Data for Public Health Investigation Based on Measured Mutation Rates and Whole-Genome Sequence Comparisons

Variable-number tandem repeats (VNTRs) mutate rapidly and can be useful markers for genotyping. While multilocus VNTR analysis (MLVA) is increasingly used in the detection and investigation of food-borne outbreaks caused by Salmonella enterica serovar Typhimurium (S. Typhimurium) and other bacterial pathogens, MLVA data analysis usually relies on simple clustering approaches that may lead to incorrect interpretations. Here, we estimated the rates of copy number change at each of the five loci commonly used for S. Typhimurium MLVA, during in vitro and in vivo passage. We found that loci STTR5, STTR6, and STTR10 changed during passage but STTR3 and STTR9 did not. Relative rates of change were consistent across in vitro and in vivo growth and could be accurately estimated from diversity measures of natural variation observed during large outbreaks. Using a set of 203 isolates from a series of linked outbreaks and whole-genome sequencing of 12 representative isolates, we assessed the accuracy and utility of several alternative methods for analyzing and interpreting S. Typhimurium MLVA data. We show that eBURST analysis was accurate and informative. For construction of MLVA-based trees, a novel distance metric, based on the geometric model of VNTR evolution coupled with locus-specific weights, performed better than the commonly used simple or categorical distance metrics. The data suggest that, for the purpose of identifying potential transmission clusters for further investigation, isolates whose profiles differ at one of the rapidly changing STTR5, STTR6, and STTR10 loci should be collapsed into the same cluster.     

MLVA as a Tool for Public Health Surveillance of Human Salmonella Typhimurium: Prospective Study in Belgium and Evaluation of MLVA Loci Stability

Surveillance of Salmonella enterica subsp. enterica serovar Typhimurium (S. Typhimurium) is generally considered to benefit from molecular techniques like multiple-locus variable-number of tandem repeats analysis (MLVA), which allow earlier detection and confinement of outbreaks. Here, a surveillance study, including phage typing, antimicrobial susceptibility testing and the in Europe most commonly used 5-loci MLVA on 1,420 S. Typhimurium isolates collected between 2010 and 2012 in Belgium, was used to evaluate the added value of MLVA for public health surveillance. Phage types DT193, DT195, DT120, DT104, DT12 and U302 dominate the Belgian S. Typhimurium population. A combined resistance to ampicillin, streptomycin, sulphonamides and tetracycline (ASSuT) with or without additional resistances was observed for 42.5% of the isolates. 414 different MLVA profiles were detected, of which 14 frequent profiles included 44.4% of the S. Typhimurium population. During a serial passage experiment on selected isolates to investigate the in vitro stability of the 5 MLVA loci, variations over time were observed for loci STTR6, STTR10, STTR5 and STTR9. This study demonstrates that MLVA improves public health surveillance of S. Typhimurium. However, the 5-loci MLVA should be complemented with other subtyping methods for investigation of possible outbreaks with frequent MLVA profiles. Also, variability in these MLVA loci should be taken into account when investigating extended outbreaks and studying dynamics over longer periods.     

Features of Variable Number of Tandem Repeats in Yersinia pestis and the Development of a Hierarchical Genotyping Scheme

Background

Variable number of tandem repeats (VNTRs) that are widely distributed in the genome of Yersinia pestis proved to be useful markers for the genotyping and source-tracing of this notorious pathogen. In this study, we probed into the features of VNTRs in the Y. pestis genome and developed a simple hierarchical genotyping system based on optimized VNTR loci.

Methodology/Principal Findings

Capillary electrophoresis was used in this study for multi-locus VNTR analysis (MLVA) in 956 Y. pestis strains. The general features and genetic diversities of 88 VNTR loci in Y. pestis were analyzed with BioNumerics, and a “14+12” loci-based hierarchical genotyping system, which is compatible with single nucleotide polymorphism-based phylogenic analysis, was established.

Conclusions/Significance

Appropriate selection of target loci reduces the impact of homoplasies caused by the rapid mutation rates of VNTR loci. The optimized “14+12” loci are highly discriminative in genotyping and source-tracing Y. pestis for molecular epidemiological or microbial forensic investigations with less time and lower cost. An MLVA genotyping datasets of representative strains will improve future research on the source-tracing and microevolution of Y. pestis.     

Molecular Strain Typing of Brucella abortus Isolates from Italy by Two VNTR Allele Sizing Technologies

Brucellosis, one of the most important re-emerging zoonoses in many countries, is caused by bacteria belonging to the genus Brucella. Furthermore these bacteria represent potential biological warfare agents and the identification of species and biovars of field strains may be crucial for tracing back source of infection, allowing to discriminate naturally occurring outbreaks instead of bioterrorist events. In the last years, multiple-locus variable-number tandem repeat analysis (MLVA) has been proposed as complement of the classical biotyping methods and it has been applied for genotyping large collections of Brucella spp. At present, the MLVA band profiles may be resolved by automated or manual procedures. The Lab on a chip technology represents a valid alternative to standard genotyping techniques (as agarose gel electrophoresis) and it has been previously used for Brucella genotyping. Recently, a new high-throughput genotyping analysis system based on capillary gel electrophoresis, the QIAxcel, has been described. The aim of the study was to evaluate the ability of two DNA sizing equipments, the QIAxcel System and the Lab chip GX, to correctly call alleles at the sixteen loci including one frequently used MLVA assay for Brucella genotyping. The results confirmed that these technologies represent a meaningful advancement in high-throughput Brucella genotyping. Considering the accuracy required to confidently resolve loci discrimination, QIAxcel shows a better ability to measure VNTR allele sizes compared to LabChip GX.     

Molecular Epidemiology of Cases of Mycoplasma californicum Infection in Japan

Bovine mastitis due to Mycoplasma californicum is often accompanied by huge economic losses, and the disease spreads very quickly. An appropriate molecular epidemiological analysis is needed to prevent and control infectious disease, but molecular epidemiological analysis methods for M. californicum have not yet been reported. Here we developed a combination of multiple-locus variable-number tandem repeat analysis (MLVA) and pulsed-field gel electrophoresis (PFGE) methods, which are common genotyping methods for various bacteria, for M. californicum. The MLVA is based on four interspersed repeat units that were found in the M. californicum genome data. The MLVA using these repeat units showed sufficient discriminatory power for a molecular epidemiological analysis; i.e., a Hunter-Gaston diversity index (HGDI) of 0.949, against M. californicum strains in Japan and M. californicum strain ATCC 33461. The PFGE for M. californicum also showed sufficient discriminatory power, with an HGDI of 0.985. Strain ATCC 33461 showed MLVA profiles and pulsotypes that differed greatly from those of strains from Japan. These results indicate that MLVA and PFGE are good tools for identifying M. californicum transmission events more accurately. Our combined MLVA and PFGE analysis suggests the persistence of M. californicum infection among herds in a specific area for a long period of time, as well as the movement of cows and heifers accompanying the expansion of M. californicum infection. Failure to identify asymptomatic infected cows is suspected as one of the central causes of the present M. californicum infection scenario in Japan.     

Rapid High Resolution Genotyping of Francisella tularensis by Whole Genome Sequence Comparison of Annotated Genes (“MLST+”)

The zoonotic disease tularemia is caused by the bacterium Francisella tularensis. This pathogen is considered as a category A select agent with potential to be misused in bioterrorism. Molecular typing based on DNA-sequence like canSNP-typing or MLVA has become the accepted standard for this organism. Due to the organism’s highly clonal nature, the current typing methods have reached their limit of discrimination for classifying closely related subpopulations within the subspecies F. tularensis ssp. holarctica. We introduce a new gene-by-gene approach, MLST⁺, based on whole genome data of 15 sequenced F. tularensis ssp. holarctica strains and apply this approach to investigate an epidemic of lethal tularemia among non-human primates in two animal facilities in Germany. Due to the high resolution of MLST⁺ we are able to demonstrate that three independent clones of this highly infectious pathogen were responsible for these spatially and temporally restricted outbreaks.     

Significance of the Identification in the Horn of Africa of an Exceptionally Deep Branching Mycobacterium tuberculosis Clade

Molecular and phylogeographic studies have led to the definition within the Mycobacterium tuberculosis complex (MTBC) of a number of geotypes and ecotypes showing a preferential geographic location or host preference. The MTBC is thought to have emerged in Africa, most likely the Horn of Africa, and to have spread worldwide with human migrations. Under this assumption, there is a possibility that unknown deep branching lineages are present in this region. We genotyped by spoligotyping and multiple locus variable number of tandem repeats (VNTR) analysis (MLVA) 435 MTBC isolates recovered from patients. Four hundred and eleven isolates were collected in the Republic of Djibouti over a 12 year period, with the other 24 isolates originating from neighbouring countries. All major M. tuberculosis lineages were identified, with only two M. africanum and one M. bovis isolates. Upon comparison with typing data of worldwide origin we observed that several isolates showed clustering characteristics compatible with new deep branching. Whole genome sequencing (WGS) of seven isolates and comparison with available WGS data from 38 genomes distributed in the different lineages confirms the identification of ancestral nodes for several clades and most importantly of one new lineage, here referred to as lineage 7. Investigation of specific deletions confirms the novelty of this lineage, and analysis of its precise phylogenetic position indicates that the other three superlineages constituting the MTBC emerged independently but within a relatively short timeframe from the Horn of Africa. The availability of such strains compared to the predominant lineages and sharing very ancient ancestry will open new avenues for identifying some of the genetic factors responsible for the success of the modern lineages. Additional deep branching lineages may be readily and efficiently identified by large-scale MLVA screening of isolates from sub-Saharan African countries followed by WGS analysis of a few selected isolates.     

新疆喀什维吾尔族结核病患者结核分枝杆菌5个VNTRs基因分型研究

目的:应用多位点数目可变串联重复序列分析(multiple loci VNTR analysis,MLVA)技术,对新疆喀什地区维吾尔族结核病患者结核分枝杆菌临床分离株进行基因分型,探讨5个数目可变串联重复序列(VNTR)基因型种类及其分布。方法:收集结核分枝杆菌,采用PCR和琼脂糖凝胶电泳技术,结合BioNumerics5.0软件,对其5个VNTR位点进行结果分析。结果:分离出58株结核分枝杆菌,分为4个基因群21个基因型,分别为Ⅰ群占19.1%,含7个基因型;Ⅱ群占3.4%,含2个基因型;Ⅲ群占67.2%,含9个基因型;Ⅳ群占10.3%,含5个基因型。结论:新疆喀什地区维吾尔族结核病患者的结核分枝杆菌存在明显的基因多态性,且存在主要流行菌群。     

Comparison of Minisatellite Polymorphisms in the Bacillus cereus Complex: a Simple Assay for Large-Scale Screening and Identification of Strains Most Closely Related to Bacillus anthracis

Polymorphism of five tandem repeats that are monomorphic in Bacillus anthracis was investigated in 230 isolates of the B. cereus group and in 5 sequenced B. cereus genomes in search for markers allowing identification of B. cereus and B. thuringiensis strains most closely related to B. anthracis. Using this multiple-locus variable number of tandem repeat analysis (MLVA), a cluster of 30 strains was selected for further characterization. Eventually, six of these were characterized by multilocus sequence type analysis. One of the strains is only six point mutations (of almost 3,000 bp) away from B. anthracis and was also proposed to be closest to B. anthracis by MLVA analysis. However, this strain remains separated from B. anthracis by a number of significant genetic events observed in B. anthracis, including the loss of the hemolysin activity, the presence of four prophages, and the presence of the two virulence plasmids, pXO1 and pXO2. One particular minisatellite marker provides an efficient assay to identify the subset of B. cereus and B. thuringiensis strains closely related to B. anthracis. Based on these results, a very simple assay is proposed that allows the screening of hundreds of strains from the B. cereus complex, with modest equipment and at a low cost, to eventually fill the gap with B. anthracis and better understand the origin and making of this dangerous pathogen.     

Development of Multiple-Locus Variable-Number Tandem-Repeat Analysis for the Molecular Subtyping of Enterobacter sakazakii

The genomic content of Enterobacter sakazakii strain ATCC BAA-894 was analyzed for variable-number tandem repeats (VNTRs). In this study we report the development of a multiple-locus VNTR analysis (MLVA) strategy for the subtyping of E. sakazakii. The method is based on a GeneScan analysis of four VNTR loci labeled with multiple fluorescent dyes. This approach was applied to a collection of 112 isolates representing all 16 of the currently defined E. sakazakii biogroups. MLVA successfully discriminated among these isolates and compared favorably with pulsed-field gel electrophoresis. The method was relatively fast and easy to perform. The potential value of MLVA as an epidemiological tool is discussed.     

Identification of polymorphic tandem repeats by direct comparison of genome sequence from different bacterial strains : a web-based resource

Background

Polymorphic tandem repeat typing is a new generic technology which has been proved to be very efficient for bacterial pathogens such as B. anthracis, M. tuberculosis, P. aeruginosa, L. pneumophila, Y. pestis. The previously developed tandem repeats database takes advantage of the release of genome sequence data for a growing number of bacteria to facilitate the identification of tandem repeats. The development of an assay then requires the evaluation of tandem repeat polymorphism on well-selected sets of isolates. In the case of major human pathogens, such as S. aureus, more than one strain is being sequenced, so that tandem repeats most likely to be polymorphic can now be selected in silico based on genome sequence comparison.

Results

In addition to the previously described general Tandem Repeats Database, we have developed a tool to automatically identify tandem repeats of a different length in the genome sequence of two (or more) closely related bacterial strains. Genome comparisons are pre-computed. The results of the comparisons are parsed in a database, which can be conveniently queried over the internet according to criteria of practical value, including repeat unit length, predicted size difference, etc. Comparisons are available for 16 bacterial species, and the orthopox viruses, including the variola virus and three of its close neighbors.

Conclusions

We are presenting an internet-based resource to help develop and perform tandem repeats based bacterial strain typing. The tools accessible at http://minisatellites.u-psud.fr now comprise four parts. The Tandem Repeats Database enables the identification of tandem repeats across entire genomes. The Strain Comparison Page identifies tandem repeats differing between different genome sequences from the same species. The "Blast in the Tandem Repeats Database" facilitates the search for a known tandem repeat and the prediction of amplification product sizes. The "Bacterial Genotyping Page" is a service for strain identification at the subspecies level.

     

Phylogeography of Francisella tularensis ssp. holarctica in France

Aim: To investigate the phylogeography of French Francisella tularensis ssp. holarctica isolates. Methods and Results: Canonical SNPs and MLVA were used to genotype 103 French F. tularensis ssp. holarctica isolates. We confirmed the presence of one subclade, the central and western European group (B.Br.FTNF002‐00), and identified four major MLVA genotypes with no obvious geographical differentiation. Conclusions: The lack of geographical resolution among MLVA genotypes suggests rapid dispersal, convergent evolution or a combination of the two. Significance and Impact of the Study: This study expands knowledge of the phylogeography of one of the two dominant European F. tularensis ssp. holarctica subclades and illustrates the need for additional SNP discovery within this subclade.     

Multi Locus Variable-Number Tandem Repeat (MLVA) Typing Tools Improved the Surveillance of Salmonella Enteritidis: A 6 Years Retrospective Study

Surveillance of Salmonella enterica subsp. enterica serovar Enteritidis is generally considered to benefit from molecular techniques like multiple-locus variable-number of tandem repeats analysis (MLVA), which allow early detection and confinement of outbreaks. Here, a surveillance study, including phage typing, antimicrobial susceptibility testing and MLVA on 1,535 S. Enteritidis isolates collected between 2007 and 2012, was used to evaluate the added value of MLVA for public health surveillance in Belgium. Phage types PT4, PT8, PT21, PT1, PT6, PT14b, PT28 and PT13 dominate the Belgian S. Enteritidis population. The isolates of S. Enteritidis were most frequently susceptible to all antibiotics tested. 172 different MLVA profiles were detected, of which 9 frequent profiles included 67.2% of the S. Enteritidis population. During a serial passage experiment on selected isolates to investigate the in vitro stability of the 5 MLVA loci, no variations over time were observed indicating that the MLVA profiles were stable. The MLVA profile of isolates originating from different outbreaks in the Democratic Republic of the Congo (DRC) between 2010 and 2011 were distinct from any of the MLVA profiles found in Belgian isolates throughout the six year observational period and demonstrates that MLVA improves public health surveillance of S. Enteritidis. However, MLVA should be complemented with other subtyping methods when investigating outbreaks is caused by the most common MLVA profile.     

Modified MLVA for Genotyping Queensland Invasive Streptococcus pneumoniae

Background

Globally, over 800 000 children under five die each year from infectious diseases caused by Streptococcus pneumoniae. To understand genetic relatedness between isolates, study transmission routes, assess the impact of human interventions e.g. vaccines, and determine infection sources, genotyping methods are required. The ‘gold standard’ genotyping method, Multi-Locus Sequence Typing (MLST), is useful for long-term and global studies. Another genotyping method, Multi-Locus Variable Number of Tandem Repeat Analysis (MLVA), has emerged as a more discriminatory, inexpensive and faster technique; however there is no universally accepted method and it is currently suitable for short-term and localised epidemiology studies. Currently Australia has no national MLST database, nor has it adopted any MLVA method for short-term or localised studies. This study aims to improve S. pneumoniae genotyping methods by modifying the existing MLVA techniques to be more discriminatory, faster, cheaper and technically less demanding than previously published MLVA methods and MLST.

Methods

Four different MLVA protocols, including a modified method, were applied to 317 isolates of serotyped invasive S. pneumoniae isolated from sterile body sites of Queensland children under 15 years from 2007–2012. MLST was applied to 202 isolates for comparison.

Results

The modified MLVA4 is significantly more discriminatory than the ‘gold standard’ MLST method. MLVA4 has similar discrimination compared to other MLVA techniques in this study). The failure to amplify particular loci in previous MLVA methods were minimised in MLVA4. Failure to amplify BOX-13 and Spneu19 were found to be serotype specific.

Conclusion

We have modified a highly discriminatory MLVA technique for genotyping Queensland invasive S. pneumoniae. MLVA4 has the ability to enhance our understanding of the pneumococcal epidemiology and the changing genetics of the pneumococcus in localised and short-term studies.