Long-read sequencing for reliably calling the mompS allele in Legionella pneumophila sequence-based typing.

Sequence-based typing (SBT) of Legionella pneumophila is a valuable tool in epidemiological studies and outbreak investigations of Legionnaires' disease. In the L. pneumophila SBT scheme, mompS2 is one of seven genes that determine the sequence type (ST). The Legionella genome typically contains two copies of mompS (mompS1 and mompS2). When they are non-identical it can be challenging to determine the mompS2 allele, and subsequently the ST, from Illumina short-reads. In our collection of 233 L. pneumophila genomes, there were 62 STs, 18 of which carried non-identical mompS copies. Using short-reads, the mompS2 allele was misassembled or untypeable in several STs. Genomes belonging to ST154 and ST574, which carried mompS1 allele 7 and mompS2 allele 15, were assigned an incorrect mompS2 allele and/or mompS gene copy number when short-read assembled. For other isolates, mainly those carrying non-identical mompS copies, short-read assemblers occasionally failed to resolve the structure of the mompS-region, also resulting in untypeability from the short-read data. In this study, we wanted to understand the challenges we observed with calling the mompS2 allele from short-reads, assess if other short-read methods were able to resolve the mompS-region, and investigate the possibility of using long-reads to obtain the mompS alleles, and thereby perform L. pneumophila SBT from long-reads only. We found that the choice of short-read assembler had a major impact on resolving the mompS-region and thus SBT from short-reads, but no method consistently solved the mompS2 allele. By using Oxford Nanopore Technology (ONT) sequencing together with Trycycler and Medaka for long-read assembly and polishing we were able to resolve the mompS copies and correctly identify the mompS2 allele, in accordance with Sanger sequencing/EQA results for all tested isolates (n=35). The remaining six genes of the SBT profile could also be determined from the ONT-only reads. The STs called from ONT-only assemblies were also consistent with hybrid-assemblies of Illumina and ONT reads. We therefore propose ONT sequencing as an alternative method to perform L. pneumophila SBT to overcome the mompS challenge observed with short-reads. To facilitate this, we have developed ONTmompS (https://github.com/marithetland/ONTmompS), an in silico approach to determine L. pneumophila ST from long-read or hybrid assemblies.

Legionella pneumophila in Municipal Shower Systems in Stavanger, Norway; A Longitudinal Surveillance Study Using Whole Genome Sequencing in Risk Management.

Following an incidence of Legionnaires disease (LD) in 2007, where a municipal shower system was the likely source of infection, Stavanger municipality initiated a surveillance program for Legionella as part of establishing internal risk evaluation and prevention routines. More than 250 shower systems were examined for cultivatable Legionella pneumophila. The prevalence and diversity of serogroups (sg) and sequence types (STs) of L. pneumophila were mapped using available typing techniques over a period of more than 10 years (2010-2021). The surveillance showed an overall reduction in the L. pneumophila colonisation rate in municipal systems from 11 to 4.5% following prevention measures during the period, with the highest colonisation rate in complex systems (e.g., larger nursing homes and sports complexes). Further, an approximately even distribution between sg1 and 2-14 was seen. Whole genome sequencing (WGS) revealed that only a limited number of STs were detected, and they were consistent at specific locations over time. This study showed that environmental surveillance data in combination with available typing techniques and WGS can give the municipality a better tool for risk management and an overview of ST distributions that can be a valuable asset in future source investigations.