Comparison of whole genome sequencing typing tools for the typing of Belgian Legionella pneumophila outbreaks isolates

Whole genome sequencing (WGS) marks a turning point for outbreak investigations for microorganisms related to public health matters, like Legionella pneumophila (Lp). Here, we evaluated the available Lp WGS typing tools for isolates of previously documented Belgian outbreaks, as well as small groups of related and non-related isolates. One reference strain and 77 clinical and environmental isolates were evaluated. Seven isolates belong to a Sequence Type (ST) 36 outbreak in 1999 and sixteen (ten clinical, two matching environmental and four non-related controls) belong to another ST1 outbreak in 1985-1987. The remaining isolates belong to small groups of related and non-related isolates of diverse ST's. WGS was performed and data were analysed using whole genome (wg) and core genome (cg) multilocus sequence typing (MLST) with "Ridom SeqSphere + " (cgMLST), "Applied Maths-Bionumerics" (wgMLST) and the 50 loci cgMLST (CDC/ESGLI_ESCMID). Results of the three tools were concordant with the traditional Sequence Based Typing (SBT). The known outbreaks and small clusters could be detected and clear discrimination of ST1 non-related isolates was obtained. In addition, the 50 loci cgMLST allowed to classify the isolates into subtypes because almost all the 50 genes could be called in all the analysed isolates, which was not achieved by the other tools. This is a big advantage in terms of standardisation and comparison between laboratories for future epidemiological investigations. WGS allowed to analyse a large volume of samples and generated more accurate conclusions for outbreak investigations compared to other typing methods due to its higher discriminatory power and throughput.