Multiple major disease-associated clones of Legionella pneumophila have emerged recently and independently

Amoebae as training grounds for microbial pathogens

Grazing of amoebae on microorganisms represents one of the oldest predator-prey dynamic relationships in nature. It represents a genetic "melting pot" for an ancient and continuous multi-directional inter- and intra-kingdom horizontal gene transfer between amoebae and its preys, intracellular microbial residents, endosymbionts, and giant viruses, which has shaped the evolution, selection, and adaptation of microbes that evade degradation by predatory amoeba. Unicellular phagocytic amoebae are thought to be the ancient ancestors of macrophages with highly conserved eukaryotic processes. Selection and evolution of microbes within amoeba through their evolution to target highly conserved eukaryotic processes have facilitated the expansion of their host range to mammals, causing various infectious diseases. Legionella and environmental Chlamydia harbor an immense number of eukaryotic-like proteins that are involved in ubiquitin-related processes or are tandem repeats-containing proteins involved in protein-protein and protein-chromatin interactions. Some of these eukaryotic-like proteins exhibit novel domain architecture and novel enzymatic functions absent in mammalian cells, such as ubiquitin ligases, likely acquired from amoebae. Mammalian cells and amoebae may respond similarly to microbial factors that target highly conserved eukaryotic processes, but mammalian cells may undergo an accidental response to amoeba-adapted microbial factors. We discuss specific examples of microbes that have evolved to evade amoeba predation, including the bacterial pathogens- Legionella, Chlamydia, Coxiella, Rickettssia, Francisella, Mycobacteria, Salmonella, Bartonella, Rhodococcus, Pseudomonas, Vibrio, Helicobacter, Campylobacter, and Aliarcobacter. We also discuss the fungi Cryptococcus, and Asperigillus, as well as amoebae mimiviruses/giant viruses. We propose that amoeba-microbe interactions will continue to be a major "training ground" for the evolution, selection, adaptation, and emergence of microbial pathogens equipped with unique pathogenic tools to infect mammalian hosts. However, our progress will continue to be highly dependent on additional genomic, biochemical, and cellular data of unicellular eukaryotes.

Recurrence, Microevolution, and Spatiotemporal Dynamics of Legionella pneumophila Sequence Type 1905, Portugal, 2014-2022

We investigated molecular evolution and spatiotemporal dynamics of atypical Legionella pneumophila serogroup 1 sequence type 1905 and determined its long-term persistence and linkage to human disease in dispersed locations, far beyond the large 2014 outbreak epicenter in Portugal. Our finding highlights the need for public health interventions to prevent further disease spread.

From Investigating a Case of Cellulitis to Exploring Nosocomial Infection Control of ST1 Legionella pneumophila Using Genomic Approaches

Legionella pneumophila can cause a large panel of symptoms besides the classic pneumonia presentation. Here we present a case of fatal nosocomial cellulitis in an immunocompromised patient followed, a year later, by a second case of Legionnaires' disease in the same ward. While the first case was easily assumed as nosocomial based on the date of symptom onset, the second case required clear typing results to be assigned either as nosocomial and related to the same environmental source as the first case, or community acquired. To untangle this specific question, we applied core-genome multilocus typing (MLST), whole-genome single nucleotide polymorphism and whole-genome MLST methods to a collection of 36 Belgian and 41 international sequence-type 1 (ST1) isolates using both thresholds recommended in the literature and tailored threshold based on local epidemiological data. Based on the thresholds applied to cluster isolates together, the three methods gave different results and no firm conclusion about the nosocomial setting of the second case could been drawn. Our data highlight that despite promising results in the study of outbreaks and for large-scale epidemiological investigations, next-generation sequencing typing methods applied to ST1 outbreak investigation still need standardization regarding both wet-lab protocols and bioinformatics. A deeper evaluation of the L. pneumophila evolutionary clock is also required to increase our understanding of genomic differences between isolates sampled during a clinical infection and in the environment.

Severe Legionnaires' disease

BACKGROUND

Legionnaires' disease (LD) is a common but under-diagnosed cause of community-acquired pneumonia (CAP), although rapid detection of urine antigen testing (UAT) and advances in molecular testing have improved the diagnosis. LD entails intensive care unit (ICU) admission in almost one-third of cases, and the mortality rate ranges from 4% to 40%. This review aims to discuss recent advances in the study of this condition and to provide an update on the diagnosis, pathogenesis and management of severe LD.

RESULTS

The overall incidence of LD has increased worldwide in recent years due to the higher number of patients with risk factors, especially immunosuppression, and to improvements in diagnostic methods. Although LD is responsible for only around 5% of all-cause CAP, it is one of the three most common causes of CAP requiring ICU admission. Mortality in ICU patients, immunocompromised patients or patients with a nosocomial source of LD can reach 40% despite appropriate antimicrobial therapy. Regarding pathogenesis, no Legionella-specific virulence factors have been associated with severity; however, recent reports have found high pulmonary Legionella DNA loads, and impairments in immune response and lung microbiome in the most severe cases. The clinical picture includes severe lung injury requiring respiratory and/or hemodynamic support, extrapulmonary symptoms and non-specific laboratory findings. LD diagnostic methods have improved due to the broad use of UAT and the development of molecular methods allowing the detection of all Lp serogroups. Therapy is currently based on macrolides, quinolones, or a combination of the two, with prolonged treatment in severe cases.

CONCLUSIONS

Numerous factors influence the mortality rate of LD, such as ICU admission, the underlying immune status, and the nosocomial source of the infection. The host immune response (hyperinflammation and/or immunoparalysis) may also be associated with increased severity. Given that the incidence of LD is rising, studies on specific biomarkers of severity may be of great interest. Further assessments comparing different regimens and/or evaluating host-directed therapies are nowadays needed.

High performance Legionella pneumophila source attribution using genomics-based machine learning classification

Fundamental to effective Legionnaires' disease outbreak control is the ability to rapidly identify the environmental source(s) of the causative agent, Legionella pneumophila. Genomics has revolutionized pathogen surveillance, but L. pneumophila has a complex ecology and population structure that can limit source inference based on standard core genome phylogenetics. Here, we present a powerful machine learning approach that assigns the geographical source of Legionnaires' disease outbreaks more accurately than current core genome comparisons. Models were developed upon 534 L. pneumophila genome sequences, including 149 genomes linked to 20 previously reported Legionnaires' disease outbreaks through detailed case investigations. Our classification models were developed in a cross-validation framework using only environmental L. pneumophila genomes. Assignments of clinical isolate geographic origins demonstrated high predictive sensitivity and specificity of the models, with no false positives or false negatives for 13 out of 20 outbreak groups, despite the presence of within-outbreak polyclonal population structure. Analysis of the same 534-genome panel with a conventional phylogenomic tree and a core genome multi-locus sequence type allelic distance-based classification approach revealed that our machine learning method had the highest overall classification performance-agreement with epidemiological information. Our multivariate statistical learning approach maximizes the use of genomic variation data and is thus well-suited for supporting Legionnaires' disease outbreak investigations.IMPORTANCEIdentifying the sources of Legionnaires' disease outbreaks is crucial for effective control. Current genomic methods, while useful, often fall short due to the complex ecology and population structure of Legionella pneumophila, the causative agent. Our study introduces a high-performing machine learning approach for more accurate geographical source attribution of Legionnaires' disease outbreaks. Developed using cross-validation on environmental L. pneumophila genomes, our models demonstrate excellent predictive sensitivity and specificity. Importantly, this new approach outperforms traditional methods like phylogenomic trees and core genome multi-locus sequence typing, proving more efficient at leveraging genomic variation data to infer outbreak sources. Our machine learning algorithms, harnessing both core and accessory genomic variation, offer significant promise in public health settings. By enabling rapid and precise source identification in Legionnaires' disease outbreaks, such approaches have the potential to expedite intervention efforts and curtail disease transmission.

Ongoing diversification of the global fish pathogen Piscirickettsia salmonis through genetic isolation and transposition bursts

The management of bacterial pathogens remains a key challenge of aquaculture. The marine gammaproteobacterium Piscirickettsia salmonis is the etiological agent of piscirickettsiosis and causes multi-systemic infections in different salmon species, resulting in considerable mortality and substantial commercial losses. Here, we elucidate its global diversity, evolution, and selection during human interventions. Our comprehensive analysis of 73 closed, high quality genome sequences covered strains from major outbreaks and was supplemented by an analysis of all P. salmonis 16S rRNA gene sequences and metagenomic reads available in public databases. Genome comparison showed that Piscirickettsia comprises at least three distinct, genetically isolated species of which two showed evidence for continuing speciation. However, at least twice the number of species exist in marine fish or seawater. A hallmark of Piscirickettsia diversification is the unprecedented amount and diversity of transposases which are particularly active in subgroups undergoing rapid speciation and are key to the acquisition of novel genes and to pseudogenization. Several group-specific genes are involved in surface antigen synthesis and may explain the differences in virulence between strains. However, the frequent failure of antibiotic treatment of piscirickettsiosis outbreaks cannot be explained by horizontal acquisition of resistance genes which so far occurred only very rarely. Besides revealing a dynamic diversification of an important pathogen, our study also provides the data for improving its surveillance, predicting the emergence of novel lineages, and adapting aquaculture management, and thereby contributes towards the sustainability of salmon farming.

Water quality influences Legionella pneumophila determination

Legionellosis is a respiratory disease of public health concern. The bacterium Legionella pneumophila is the etiologic agent responsible for >90% of legionellosis cases in the United States. Legionellosis transmission primarily occurs through the inhalation or aspiration of contaminated water aerosols or droplets. Therefore, a thorough understanding of L. pneumophila detection methods and their performance in various water quality conditions is needed to develop preventive measures. Two hundred and nine potable water samples were collected from taps in buildings across the United States. L. pneumophila was determined using three culture methods: Buffered Charcoal Yeast Extract (BCYE) culture with Matrix-assisted Laser Desorption/Ionization Mass Spectrometry (MALDI-MS) identification, Legiolert® 10- and 100-mL tests, and one molecular method: quantitative Polymerase Chain Reaction (qPCR) assay. Culture and molecular positive results were confirmed by secondary testing including MALDI-MS. Eight water quality variables were studied, including source water type, secondary disinfectant, total chlorine residual, heterotrophic bacteria, total organic carbon (TOC), pH, water hardness, cold- and hot-water lines. The eight water quality variables were segmented into 28 categories, based on scale and ranges, and method performance was evaluated in each of these categories. Additionally, a Legionella genus qPCR assay was used to determine the water quality variables that promote or hinder Legionella spp. occurrence. L. pneumophila detection frequency ranged from 2 to 22% across the methods tested. Method performance parameters of sensitivity, specificity, positive and negative predictive values, and accuracy were >94% for the qPCR method but ranged from 9 to 100% for the culture methods. Water quality influenced L. pneumophila determination by culture and qPCR methods. L. pneumophila qPCR detection frequencies positively correlated with TOC and heterotrophic bacterial counts. The source water-disinfectant combination influenced the proportion of Legionella spp. that is L. pneumophila. Water quality influences L. pneumophila determination. To accurately detect L. pneumophila, method selection should consider the water quality in addition to the purpose of testing (general environmental monitoring versus disease-associated investigations).

Rapid adaptations of Legionella pneumophila to the human host

Legionella pneumophila are host-adapted bacteria that infect and reproduce primarily in amoeboid protists. Using similar infection mechanisms, they infect human macrophages, and cause Legionnaires' disease, an atypical pneumonia, and the milder Pontiac fever. We hypothesized that, despite the similarities in infection mechanisms, the hosts are different enough that there exist high-selective value mutations that would dramatically increase the fitness of Legionella inside the human host. By comparing a large number of isolates from independent infections, we identified two genes, mutated in three unrelated patients, despite the short duration of the incubation period (2-14 days). One is a gene coding for an outer membrane protein (OMP) belonging to the OmpP1/FadL family. The other is a gene coding for an EAL-domain-containing protein involved in cyclic-di-GMP regulation, which in turn modulates flagellar activity. The clinical strain, carrying the mutated EAL-domain-containing homologue, grows faster in macrophages than the wild-type strain, and thus appears to be better adapted to the human host. As human-to-human transmission is very rare, fixation of these mutations into the population and spread into the environment is unlikely. Therefore, parallel evolution - here mutations in the same genes observed in independent human infections - could point to adaptations to the accidental human host. These results suggest that despite the ability of L. pneumophila to infect, replicate in and exit from macrophages, its human-specific adaptations are unlikely to be fixed in the population.

Comparative Genomics of Legionella pneumophila Isolates from the West Bank and Germany Support Molecular Epidemiology of Legionnaires' Disease

Legionella pneumophila is an environmental bacterium and clinical pathogen that causes many life-threating outbreaks of an atypical pneumonia called Legionnaires' disease (LD). Studies of this pathogen have focused mainly on Europe and the United States. A shortage in L. pneumophila data is clearly observed for developing countries. To reduce this knowledge gap, L. pneumophila isolates were studied in two widely different geographical areas, i.e., the West Bank and Germany. For this study, we sequenced and compared the whole genome of 38 clinical and environmental isolates of L. pneumophila covering different MLVA-8(12) genotypes in the two areas. Sequencing was conducted using the Illumina HiSeq 2500 platform. In addition, two isolates (A194 and H3) were sequenced using a Pacific Biosciences (PacBio) RSII platform to generate complete reference genomes from each of the geographical areas. Genome sequences from 55 L. pneumophila strains, including 17 reference strains, were aligned with the genome sequence of the closest strain (L. pneumophila strain Alcoy). A whole genome phylogeny based on single nucleotide polymorphisms (SNPs) was created using the ParSNP software v 1.0. The reference genomes obtained for isolates A194 and H3 consisted of circular chromosomes of 3,467,904 bp and 3,691,263 bp, respectively. An average of 36,418 SNPs (min. 8569, max. 70,708 SNPs) against our reference strain L. pneumophila str. Alcoy, and 2367 core-genes were identified among the fifty-five strains. An analysis of the genomic population structure by SNP comparison divided the fifty-five L. pneumophila strains into six branches. Individual isolates in sub-lineages in these branches differed by less than 120 SNPs if they had the same MLVA genotype and were isolated from the same location. A bioinformatics analysis identified the genomic islands (GIs) for horizontal gene transfer and mobile genetic elements, demonstrating that L. pneumophila showed high genome plasticity. Four L. pneumophila isolates (H3, A29, A129 and L10-091) contained well-defined plasmids. On average, only about half of the plasmid genes could be matched to proteins in databases. In silico phage findings suggested that 43 strains contained at least one phage. However, none of them were found to be complete. BLASTp analysis of proteins from the type IV secretion Dot/Icm system showed those proteins highly conserved, with less than 25% structural differences in the new L. pneumophila isolates. Overall, we demonstrated that whole genome sequencing provides a molecular surveillance tool for L. pneumophila at the highest conceivable discriminatory level, i.e., two to eight SNPs were observed for isolates from the same location but several years apart.

Whole-genome sequencing of the clinical isolate of Legionella pneumophila ALAW1 from the West Bank allows high-resolution typing and determination of pathogenicity mechanisms

Background

Legionella pneumophila is water-based bacterium causing Legionnaires' disease (LD). We describe the first documented case of nosocomial LD caused by L. pneumophila sequence type (ST) 461 and serogroup 6. The etiology of LD was confirmed by culturing the bronchoalveolar lavage sample retrieving L. pneumophila strain ALAW1. A 7-days treatment of the LD patient with Azithromycin and Levofloxacin allowed complete recovery.

Methods

In details, we sequenced the whole genome of the L. pneumophila ALAW1 using Illumina HiSeq platform. The sequence of ALAW1 was aligned with the genome sequence from the closely related reference strain Alcoy 2300/99 and a whole-genome phylogeny based on single nucleotide polymorphisms (SNPs) was created using Parsnp software. Also, the TYGS web-server was used in order to compare the genome with type strain.

Results

An analysis of the population structure by SNP and TYGS comparison clustered ALAW1 with the reference genome Alcoy 2300/99. Blastp analysis of the type IV secretion Dot/Icm system genes showed that these genes were highly conserved with (≤25%) structural differences at the protein level.

Conclusions

Overall, this study provides insights into detailed genome structure and demonstrated the value of whole-genome sequencing as the ultimate typing tool for Legionella.

Epidemiological analysis of Legionnaires' disease in Scotland: a genomic study

BACKGROUND

Legionella pneumophila is the main cause of a severe pneumonic illness known as Legionnaires' disease and is a global public health threat. Whole-genome sequencing (WGS) can be applied to trace environmental origins of L pneumophila infections, providing information to guide appropriate interventions. We aim to explore the evolutionary and epidemiological relationships in a 36-year Scottish L pneumophila reference isolate collection.

METHODS

We investigated the genomic epidemiology of Legionnaires' disease over 36 years in Scotland, comparing genome sequences for all clinical L pneumophila isolates (1984-2020) with a sequence dataset of 3211 local and globally representative isolates. We used a stratified clustering approach to capture epidemiological relationships by core genome Multi-locus Sequence Typing, followed by high-resolution phylogenetic analysis of clusters to measure diversity and evolutionary relatedness in context with epidemiological metadata.

FINDINGS

Clustering analysis showed that 111 (57·5 %) of 193 of L pneumophila infections in Scotland were caused by ten endemic lineages with a wide temporal and geographical distribution. Phylogenetic analysis of L pneumophila identified hospital-associated sublineages that had been detected in the hospital environment up to 19 years. Furthermore, 12 (30·0%) of 40 community-associated infections (excluding a single, large outbreak) that occurred over a 13 year period (from 2000 to 2013) were caused by a single widely distributed endemic clone (ST37), consistent with enhanced human pathogenicity. Finally, our analysis revealed clusters linked by national or international travel to distinct geographical regions, indicating several previously unrecognised travel links between closely related isolates (fewer than five single nucleotide polymorphisms) connected by geography.

INTERPRETATION

Our analysis reveals the existence of previously undetected endemic clones of L pneumophila that existed for many years in hospital, community, and travel-associated environments. In light of these findings, we propose that cluster and outbreak definitions should be reconsidered, and propose WGS-based surveillance as a critical public health tool for real-time identification and mitigation of clinically important endemic clones.

FUNDING

Chief Scientist Office, Biotechnology and Biological Sciences Research Council (UK), Medical Research Council Precision Medicine Doctoral Training Programme, Wellcome Trust, and Medical Research Council (UK).

Clinical challenges treating Stenotrophomonas maltophilia infections: an update

Stenotrophomonas maltophilia is a non-fermenting, Gram-negative bacillus that has emerged as an opportunistic nosocomial pathogen. Its intrinsic multidrug resistance makes treating infections caused by S. maltophilia a great clinical challenge. Clinical management is further complicated by its molecular heterogeneity that is reflected in the uneven distribution of antibiotic resistance and virulence determinants among different strains, the shortcomings of available antimicrobial susceptibility tests and the lack of standardized breakpoints for the handful of antibiotics with in vitro activity against this microorganism. Herein, we provide an update on the most recent literature concerning these issues, emphasizing the impact they have on clinical management of S. maltophilia infections.

Extensive epigenetic modification with large-scale chromosomal and plasmid recombination characterise the Legionella longbeachae serogroup 1 genome

Legionella longbeachae is an environmental bacterium that is the most clinically significant Legionella species in New Zealand (NZ), causing around two-thirds of all notified cases of Legionnaires’ disease. Here we report the sequencing and analysis of the geo-temporal genetic diversity of 54 L. longbeachae serogroup 1 (sg1) clinical isolates, derived from cases from around NZ over a 22-year period, including one complete genome and its associated methylome. The 54 sg1 isolates belonged to two main clades that last shared a common ancestor between 95 BCE and 1694 CE. There was diversity at the genome-structural level, with large-scale arrangements occurring in some regions of the chromosome and evidence of extensive chromosomal and plasmid recombination. This includes the presence of plasmids derived from recombination and horizontal gene transfer between various Legionella species, indicating there has been both intra- and inter-species gene flow. However, because similar plasmids were found among isolates within each clade, plasmid recombination events may pre-empt the emergence of new L. longbeachae strains. Our complete NZ reference genome consisted of a 4.1 Mb chromosome and a 108 kb plasmid. The genome was highly methylated with two known epigenetic modifications, m⁴C and m⁶A, occurring in particular sequence motifs within the genome.

Genome analysis of Legionella pneumophila ST23 from various countries reveals highly similar strains

ST23 isolated in Italy are analysed by cgMLST and SNP approaches and they are also compared with ST23 from other countries. They are found to be phylogenetically related independently on year, town, or country of isolation.

Legionella pneumophila serogroup 1 (Lp1) sequence type (ST) 23 is one of the most commonly detected STs in Italy where it currently causes all investigated outbreaks. ST23 has caused both epidemic and sporadic cases between 1995 and 2018 and was analysed at genomic level and compared with ST23 isolated in other countries to determine possible similarities and differences. A core genome multi-locus sequence typing (cgMLST), based on a previously described set of 1,521 core genes, and single-nucleotide polymorphisms (SNPs) approaches were applied to an ST23 collection including genomes from Italy, France, Denmark and Scotland. DNAs were automatically extracted, libraries prepared using NextEra library kit and MiSeq sequencing performed. Overall, 63 among clinical and environmental Italian Lp1 isolates and a further seven and 11 ST23 from Denmark and Scotland, respectively, were sequenced, and pangenome analysed. Both cgMLST and SNPs analyses showed very few loci and SNP variations in ST23 genomes. All the ST23 causing outbreaks and sporadic cases in Italy and elsewhere, were phylogenetically related independent of year, town or country of isolation. Distances among the ST23s were further shortened when SNPs due to horizontal gene transfers were removed. The Lp1 ST23 isolated in Italy have kept their monophyletic origin, but they are phylogenetically close also to ST23 from other countries. The ST23 are quite widespread in Italy, and a thorough epidemiological investigation is compelled to determine sources of infection when this ST is identified in both LD sporadic cases and outbreaks.

The molecular basis of FimT-mediated DNA uptake during bacterial natural transformation

Naturally competent bacteria encode sophisticated protein machinery for the uptake and translocation of exogenous DNA into the cell. If this DNA is integrated into the bacterial genome, the bacterium is said to be naturally transformed. Most competent bacterial species utilise type IV pili for the initial DNA uptake step. These proteinaceous cell-surface structures are composed of thousands of pilus subunits (pilins), designated as major or minor according to their relative abundance in the pilus. Here, we show that the minor pilin FimT plays an important role in the natural transformation of Legionella pneumophila. We use NMR spectroscopy, in vitro DNA binding assays and in vivo transformation assays to understand the molecular basis of FimT's role in this process. FimT binds to DNA via an electropositive patch, rich in arginines, several of which are well-conserved and located in a conformationally flexible C-terminal tail. FimT orthologues from other Gammaproteobacteria share the ability to bind to DNA. Our results suggest that FimT plays an important role in DNA uptake in a wide range of competent species.

A Tale of Four Danish Cities: Legionella pneumophila Diversity in Domestic Hot Water and Spatial Variations in Disease Incidence

Denmark has one of the highest Legionnaires' disease notification rates within Europe, averaging 4.7 cases per 100,000 population annually (2017 to 2020). The relatively high incidence of disease is not uniform across the country, and approximately 70% of all domestically acquired cases in Denmark are caused by Legionella pneumophila (LP) strains that are considered less virulent. The aim of this study was to investigate if colonization rates, levels of colonization, and/or types of LP present in hot water systems were associated with geographic differences in Legionnaires' disease incidence. Domestic water systems from four cities in Denmark were analyzed via culture and qPCR. Serogrouping and sequence typing was performed on randomly selected isolates. Single nucleotide polymorphism was used to identify clonal relationship among isolates from the four cities. The results revealed a high LP colonization rate from 68% to 87.5% among systems, composed primarily of non-serogroup 1. LP serogroup 1 reacting with the monoclonal antibody (MAb) 3/1 was not identified in any of the systems tested, while MAb 3/1 negative serogroup 1 strains were isolated from 10 systems (9.6%). We hypothesize that a combination of factors influences the incidence rate of LD in each city, including sequence type and serogroup distribution, colonization rate, concentration of Legionella in Pre-flush and Flush samples, and potentially building characteristics such as water temperature measured at the point of use.

TNF-α response in macrophages depends on clinical Legionella pneumophila isolates genotypes

Legionnaires' Disease (LD) is a severe pneumonia mainly caused in Europe by Legionella pneumophila serogroup 1 (Lp1). Sequence-based typing methods reveal that some sequence types (ST) are overrepresented in clinical samples such as ST1 and ST47, suggesting that some strains are more fit for infection than others. In the present study, a collection of 108 Lp1 clinical isolates were used to evaluate the strain-dependent immune responses from human macrophages. Clinical Lp1 isolates induced differential TNFα secretion from macrophages. ST1 isolates induced a significantly higher TNF-α secretion than non-ST1, whereas ST47 isolates induced a significantly lower TNF-α secretion than non-ST47 isolates. ST1 isolates induced a significantly higher cell death than ST47 isolates evaluated by lactate dehydrogenase activity (cytotoxicity) and caspase-3 activity (apoptosis). Treatment of macrophages with anti-TNF-α antibodies significantly reduced the cell death in macrophages infected with ST1 or ST47 strains. The TNF-α secretion was neither explained by a differential bacterial replication nor by the number or type (bystander or infected) of TNF-α producing cells following infection but by a differential response from macrophages. The Paris ST1 reference strain elicited a significantly higher TNF-α gene transcription and a higher induction of NF-κB signaling pathway than the Lorraine ST47 reference strain.Clinical Lp1 isolates induce a diverse immune response and cell death, which could be related to the genotype. The two predominant sequence-types ST1 and ST47 trigger opposite inflammatory response that could be related to the host susceptibility.

Population analysis of Legionella pneumophila reveals a basis for resistance to complement-mediated killing

Legionella pneumophila is the most common cause of the severe respiratory infection known as Legionnaires' disease. However, the microorganism is typically a symbiont of free-living amoeba, and our understanding of the bacterial factors that determine human pathogenicity is limited. Here we carried out a population genomic study of 902 L. pneumophila isolates from human clinical and environmental samples to examine their genetic diversity, global distribution and the basis for human pathogenicity. We find that the capacity for human disease is representative of the breadth of species diversity although some clones are more commonly associated with clinical infections. We identified a single gene (lag-1) to be most strongly associated with clinical isolates. lag-1, which encodes an O-acetyltransferase for lipopolysaccharide modification, has been distributed horizontally across all major phylogenetic clades of L. pneumophila by frequent recent recombination events. The gene confers resistance to complement-mediated killing in human serum by inhibiting deposition of classical pathway molecules on the bacterial surface. Furthermore, acquisition of lag-1 inhibits complement-dependent phagocytosis by human neutrophils, and promoted survival in a mouse model of pulmonary legionellosis. Thus, our results reveal L. pneumophila genetic traits linked to disease and provide a molecular basis for resistance to complement-mediated killing.

The global population structure and beta-lactamase repertoire of the opportunistic pathogen Serratia marcescens

Serratia marcescens is a global spread nosocomial pathogen. This rod-shaped bacterium displays a broad host range and worldwide geographical distribution. Here we analyze an international collection of this multidrug-resistant, opportunistic pathogen from 35 countries to infer its population structure. We show that S. marcescens comprises 12 lineages; Sm1, Sm4, and Sm10 harbor 78.3% of the known environmental strains. Sm5, Sm6, and Sm7 comprise only human-associated strains which harbor smallest pangenomes, genomic fluidity and lowest levels of core recombination, indicating niche specialization. Sm7 and Sm9 lineages exhibit the most concerning resistome; bla_KPC-2 plasmid is widespread in Sm7, whereas Sm9, also an anthropogenic-exclusive lineage, presents highest plasmid/lineage size ratio and plasmid-diversity encoding metallo-beta-lactamases comprising bla_NDM-1. The heterogeneity of resistance patterns of S. marcescens lineages elucidated herein highlights the relevance of surveillance programs, using whole-genome sequencing, to provide insights into the molecular epidemiology of carbapenemase producing strains of this species.

Legionella pneumophila-Epidemiology and Characterization of Clinical Isolates, Slovenia, 2006-2020

Legionella pneumophila is the causative agent of severe Legionnaires’ disease (LD). Although an increasing number of LD cases have been observed, published data from Slovenia are very limited and data on molecular epidemiology are even scarcer. The present retrospective study (2006–2020) reports the results of the microbiological diagnosis of LD, as well as the epidemiology and characterization of the Legionella clinical isolates. We tested urine samples from 15,540 patients with pneumonia symptoms for L. pneumophila infection by urine antigen test, of which 717 (4.6%) tested positive. Isolation of L. pneumophila was successfully performed from 88 clinical specimens, with 82 (93.2%) being identified as L. pneumophila sg 1 and six (6.8%) as L. pneumophila sg 2–14. Sequence-based typing (SBT) identified 33 different sequence types (STs), the most frequent being ST1 and ST23. Sequence type 1 mainly comprised isolates belonging to the Philadelphia subgroup, and ST23 mostly to Allentown/France. The standard SBT scheme, as well as Dresden phenotyping for L. pneumophila, presented a high diversity among isolates.

Transposon Insertion Sequencing in a Clinical Isolate of Legionella pneumophila Identifies Essential Genes and Determinants of Natural Transformation

Legionella pneumophila is the etiologic agent of a severe form of nosocomial and community-acquired pneumonia in humans. The environmental life traits of L. pneumophila are essential to its ability to accidentally infect humans.

Legionella pneumophila is a Gram-negative bacterium ubiquitous in freshwater environments which, if inhaled, can cause a severe pneumonia in humans. The emergence of L. pneumophila is linked to several traits selected in the environment, the acquisition of some of which involved intra- and interkingdom horizontal gene transfer events. Transposon insertion sequencing (TIS) is a powerful method to identify the genetic basis of selectable traits as well as to identify fitness determinants and essential genes, which are possible antibiotic targets. TIS has not yet been used to its full power in L. pneumophila, possibly because of the difficulty of obtaining a high-saturation transposon insertion library. Indeed, we found that isolates of sequence type 1 (ST1), which includes the commonly used laboratory strains, are poorly permissive to saturating mutagenesis by conjugation-mediated transposon delivery. In contrast, we obtained high-saturation libraries in non-ST1 clinical isolates, offering the prospect of using TIS on unaltered L. pneumophila strains. Focusing on one of them, we then used TIS to identify essential genes in L. pneumophila. We also revealed that TIS could be used to identify genes controlling vertical transmission of mobile genetic elements. We then applied TIS to identify all the genes required for L. pneumophila to develop competence and undergo natural transformation, defining the set of major and minor type IV pilins that are engaged in DNA uptake. This work paves the way for the functional exploration of the L. pneumophila genome by TIS and the identification of the genetic basis of other life traits of this species.

IMPORTANCELegionella pneumophila is the etiologic agent of a severe form of nosocomial and community-acquired pneumonia in humans. The environmental life traits of L. pneumophila are essential to its ability to accidentally infect humans. A comprehensive identification of their genetic basis could be obtained through the use of transposon insertion sequencing. However, this powerful approach had not been fully implemented in L. pneumophila. Here, we describe the successful implementation of the transposon-sequencing approach in a clinical isolate of L. pneumophila. We identify essential genes, potential drug targets, and genes required for horizontal gene transfer by natural transformation. This work represents an important step toward identifying the genetic basis of the many life traits of this environmental and pathogenic species.

Slowly or Nonresolving Legionnaires' Disease: Case Series and Literature Review

BACKGROUND

Rarely, Legionnaires' disease (LD) can progress into a slowly or nonresolving form.

METHODS

A nationwide retrospective study was conducted by the French National Reference Center for Legionella (2013-2017) including cases of slowly or nonresolving LD defined as persistent clinical symptoms, computed tomography (CT) scan abnormalities, and Legionella detection in lower respiratory tract specimens by culture and/or real-time (RT) polymerase chain reaction (PCR) >30 days after symptom onset.

RESULTS

Twelve cases of community-acquired slowly or nonresolving LD were identified among 1686 cases of culture-positive LD. Median (interquartile range [IQR]) age was 63 (29-82) years. Ten (83.3%) patients had ≥1 immunosuppressive factor. Clinically, 9 patients transiently recovered before further deterioration (median [IQR] symptom-free interval, 30 [18-55] days), 3 patients had uniformly persistent symptoms (median [IQR] time, 48 [41.5-54] days). Two patients had >2 recurrences. CT scan imagery found lung abscess in 5 (41.6%) cases. Slowly or nonresolving LD was diagnosed on positive Legionella cultures (n = 10, 83.3%) at 49.5 (IQR, 33.7-79) days. Two cases were documented through positive Legionella RT PCR at 52 and 53 days (cycle threshold detection of 21.5 and 33.7, respectively). No genomic microevolution and no Legionella resistance to antibiotics were detected. The median (IQR) duration of treatment was 46.5 (21-92.5) days. Two empyema cases required thoracic surgery. At a median (IQR) follow-up of 26 (14-41.5) months, LD-attributable mortality was 16.6% (n = 2).

CONCLUSIONS

Slowly or nonresolving LD may occur in immunocompromised patients, possibly leading to lung abscess and empyema.

Guidelines for the use and interpretation of assays for monitoring autophagy (4th edition)1

In 2008, we published the first set of guidelines for standardizing research in autophagy. Since then, this topic has received increasing attention, and many scientists have entered the field. Our knowledge base and relevant new technologies have also been expanding. Thus, it is important to formulate on a regular basis updated guidelines for monitoring autophagy in different organisms. Despite numerous reviews, there continues to be confusion regarding acceptable methods to evaluate autophagy, especially in multicellular eukaryotes. Here, we present a set of guidelines for investigators to select and interpret methods to examine autophagy and related processes, and for reviewers to provide realistic and reasonable critiques of reports that are focused on these processes. These guidelines are not meant to be a dogmatic set of rules, because the appropriateness of any assay largely depends on the question being asked and the system being used. Moreover, no individual assay is perfect for every situation, calling for the use of multiple techniques to properly monitor autophagy in each experimental setting. Finally, several core components of the autophagy machinery have been implicated in distinct autophagic processes (canonical and noncanonical autophagy), implying that genetic approaches to block autophagy should rely on targeting two or more autophagy-related genes that ideally participate in distinct steps of the pathway. Along similar lines, because multiple proteins involved in autophagy also regulate other cellular pathways including apoptosis, not all of them can be used as a specific marker for bona fide autophagic responses. Here, we critically discuss current methods of assessing autophagy and the information they can, or cannot, provide. Our ultimate goal is to encourage intellectual and technical innovation in the field.

Biogeography and Environmental Drivers of Legionella pneumophila Abundance and Genotype Composition Across the West Bank: Relevance of a Genotype-Based Ecology for Understanding Legionella Occurrence

The West Bank can be considered as a high-risk area for Legionella prevalence in drinking water due to high ambient temperature, intermittent water supply, frequent pressure loss, and storage of drinking water in roof containers. To assess occurrence of Legionella species, especially L. pneumophila, in the drinking water of the West Bank, the drinking water distribution systems of eight hospitals were sampled over a period of 2.3 years covering the seasonal cycle and the major geographic regions. To gain insight into potential environmental drivers, a set of physico-chemical and microbiological parameters was recorded. Sampling included drinking water and biofilm analyzed by culture and PCR-based methods. Cultivation led to the isolation of 180 strains of L. pneumophila that were genotyped by Multi-Locus Variable Number of Tandem Repeat Analysis (MLVA). Surprisingly, the abundance of culturable L. pneumophila was low in drinking water of the sampling sites, with only three out of eight sites where Legionella was observed at all (range: 30–500 CFU/Liter). By contrast, biofilm and PCR-based analyses showed a higher prevalence. Statistical analyses with physico-chemical parameters revealed a decrease of L. pneumophila abundance for water and biofilm with increasing magnesium concentrations (>30 mg/L). MLVA-genotype analysis of the L. pneumophila isolates and their spatial distribution indicated three niches characterized by distinct physico-chemical parameters and inhabited by specific consortia of genotypes. This study provides novel insights into mechanisms shaping L. pneumophila populations and triggering their abundance leading to an understanding of their genotype-specific niches and ecology in support of improved prevention measures.

Characterization of Legionella pneumophila Populations by Multilocus Variable Number of Tandem Repeats (MLVA) Genotyping from Drinking Water and Biofilm in Hospitals from Different Regions of the West Bank

The West Bank can be considered a high-risk area for Legionnaires' disease (LD) due to its hot climate, intermittent water supply and roof storage of drinking water. Legionella, mostly L. pneumophila, are responsible for LD, a severe, community-acquired and nosocomial pneumonia. To date, no extensive assessment of Legionella spp and L. pneumophila using cultivation in combination with molecular approaches in the West Bank has been published. Two years of environmental surveillance of Legionella in water and biofilms in the drinking water distribution systems (DWDS) of eight hospitals was carried out; 180 L. pneumophila strains were isolated, mostly from biofilms in DWDS. Most of the isolates were identified as serogroup (Sg) 1 (60%) and 6 (30%), while a minor fraction comprised Sg 8 and 10. Multilocus Variable number of tandem repeats Analysis using 13 loci (MLVA-8(12)) was applied as a high-resolution genotyping method and compared to the standard Sequence Based Typing (SBT). The isolates were genotyped in 27 MLVA-8(12) genotypes (Gt), comprising four MLVA clonal complexes (VACC 1; 2; 5; 11). The major fraction of isolates constituted Sequence Type (ST)1 and ST461. Most of the MLVA-genotypes were highly diverse and often unique. The MLVA-genotype composition showed substantial regional variability. In general, the applied MLVA-method made it possible to reproducibly genotype the isolates, and was consistent with SBT but showed a higher resolution. The advantage of the higher resolution was most evident for the subdivision of the large strain sets of ST1 and ST461; these STs were shown to be highly pneumonia-relevant in a former study. This shows that the resolution by MLVA is advantageous for back-tracking risk sites and for the avoidance of outbreaks of L. pneumophila. Overall, our results provide important insights into the detailed population structure of L. pneumophila, allowing for better risk assessment for DWDS.

Air-conditioner cooling towers as complex reservoirs and continuous source of Legionella pneumophila infection evidenced by a genomic analysis study in 2017, Switzerland

Introduction

Water supply and air-conditioner cooling towers (ACCT) are potential sources of Legionella pneumophila infection in people. During outbreaks, traditional typing methods cannot sufficiently segregate L. pneumophila strains to reliably trace back transmissions to these artificial water systems. Moreover, because multiple L. pneumophila strains may be present within these systems, methods to adequately distinguish strains are needed. Whole genome sequencing (WGS) and core genome multilocus sequence typing (cgMLST), with their higher resolution are helpful in this respect. In summer 2017, the health administration of the city of Basel detected an increase of L. pneumophila infections compared with previous months, signalling an outbreak.

Aim

We aimed to identify L. pneumophila strains populating suspected environmental sources of the outbreak, and to assess the relations between these strains and clinical outbreak strains.

Methods

An epidemiological and WGS-based microbiological investigation was performed, involving isolates from the local water supply and two ACCTs (n = 60), clinical outbreak and non-outbreak related isolates from 2017 (n = 8) and historic isolates from 2003–2016 (n = 26).

Results

In both ACCTs, multiple strains were found. Phylogenetic analysis of the ACCT isolates showed a diversity of a few hundred allelic differences in cgMLST. Furthermore, two isolates from one ACCT showed no allelic differences to three clinical isolates from 2017. Five clinical isolates collected in the Basel area in the last decade were also identical in cgMLST to recent isolates from the two ACCTs.

Conclusion

Current outbreak-related and historic isolates were linked to ACCTs, which form a complex environmental habitat where strains are conserved over years.

A Snapshot of the Prevalence and Molecular Diversity of Legionella pneumophila in the Water Systems of Israeli Hotels

Exposure to Legionella spp. contaminated aerosols in hotel settings confers risk for travel-associated Legionnaire’s disease (TALD). In this study, we investigated the prevalence of Legionella contamination and its molecular diversity in hotels and resorts across Israel. The study was comprised of a convenience sample of water systems from 168 hotels and resorts countrywide, routinely inspected between March 2015 and February 2017. Isolation and quantitation of Legionella were performed in a water laboratory using the ISO 11731 method. The distribution of Legionella isolates was analyzed according to geography and source. The genetic diversity of a subset of isolates was analyzed by sequence-based typing (SBT) at the National Reference Laboratory for Legionella and compared to the national database. Out of 2830 samples tested, 470 (17%) obtained from 102 different premises (60% of hotels) were positive for Legionella spp. In 230 samples (49% of all positive, 8% of total samples), accounting for 37% of hotels, Legionella spp. counts exceeded the regulatory threshold of 1000 CFU/L. The most frequently contaminated water sources were cooling towers (38%), followed by faucets, hot tubs, water lines, and storage tanks (14–17% each). Furthermore, 32% and 17% of samples obtained from cooling towers and hot tubs, respectively, exceeded the regulatory thresholds. SBT was performed on 78 strains and revealed 27 different sequence types (STs), including two novel STs. The most prevalent STs found were ST1 (26%), ST87 (10%), ST93 (6%), and ST461 and ST1516 (5% each). Several L. pneumophila STs were found to be limited to certain geographical regions. This is the first study to investigate the prevalence and diversity of Legionella in hotels and resorts in Israel during non-outbreak environmental inspections. These findings will inform risk assessment, surveillance, and control measures of TALD.

Persistent Legionnaires' Disease and Associated Antibiotic Treatment Engender a Highly Disturbed Pulmonary Microbiome Enriched in Opportunistic Microorganisms

Despite the importance of pneumonia to public health, little is known about the composition of the lung microbiome during infectious diseases, such as pneumonia, and how it evolves during antibiotic therapy. To study the possible relation of the pulmonary microbiome to the severity and outcome of this respiratory disease, we analyzed the dynamics of the pathogen and the human lung microbiome during persistent infections caused by the bacterium Legionella pneumophila and their evolution during antimicrobial treatment. We collected 10 bronchoalveolar lavage fluid samples from three patients during long-term hospitalization due to pneumonia and performed a unique longitudinal study of the interkingdom microbiome, analyzing the samples for presence of bacteria, archaea, fungi, and protozoa by high-throughput Illumina sequencing of marker genes. The lung microbiome of the patients was characterized by a strong predominance of the pathogen, a low diversity of the bacterial fraction, and an increased presence of opportunistic microorganisms. The fungal fraction was more stable than the bacterial fraction. During long-term treatment, no genomic changes or antibiotic resistance-associated mutations that could explain the persistent infection occurred, according to whole-genome sequencing analyses of the pathogen. After antibiotic treatment, the microbiome did not recover rapidly but was mainly constituted of antibiotic-resistant species and enriched in bacteria, archaea, fungi, or protozoa associated with pathogenicity. The lung microbiome seems to contribute to nonresolving Legionella pneumonia, as it is strongly disturbed during infection and enriched in opportunistic and/or antibiotic-resistant bacteria and microorganisms, including fungi, archaea, and protozoa that are often associated with infections.IMPORTANCE The composition and dynamics of the lung microbiome during pneumonia are not known, although the lung microbiome might influence the severity and outcome of this infectious disease, similar to what was shown for the microbiome at other body sites. Here we report the findings of a comprehensive analysis of the lung microbiome composition of three patients with long-term pneumonia due to L. pneumophila and its evolution during antibiotic treatment. This work adds to our understanding of how the microbiome changes during disease and antibiotic treatment and points to microorganisms and their interactions that might be beneficial. In addition to bacteria and fungi, our analyses included archaea and eukaryotes (protozoa), showing that both are present in the pulmonary microbiota and that they might also play a role in the response to the microbiome disturbance.

The phylogenetic landscape and nosocomial spread of the multidrug-resistant opportunist Stenotrophomonas maltophilia

Recent studies portend a rising global spread and adaptation of human- or healthcare-associated pathogens. Here, we analyse an international collection of the emerging, multidrug-resistant, opportunistic pathogen Stenotrophomonas maltophilia from 22 countries to infer population structure and clonality at a global level. We show that the S. maltophilia complex is divided into 23 monophyletic lineages, most of which harbour strains of all degrees of human virulence. Lineage Sm6 comprises the highest rate of human-associated strains, linked to key virulence and resistance genes. Transmission analysis identifies potential outbreak events of genetically closely related strains isolated within days or weeks in the same hospitals.

Multidrug resistance of the opportunistic pathogen Stenotrophomonas maltophilia is an increasing problem. Here, analyzing strains from 22 countries, the authors show that the S. maltophilia complex is divided into 23 monophyletic lineages and find evidence for intra-hospital transmission.

Legionnaires’ Disease: State of the Art Knowledge of Pathogenesis Mechanisms of Legionella

Legionella species are environmental gram-negative bacteria able to cause a severe form of pneumonia in humans known as Legionnaires’ disease. Since the identification of Legionella pneumophila in 1977, four decades of research on Legionella biology and Legionnaires’ disease have brought important insights into the biology of the bacteria and the molecular mechanisms that these intracellular pathogens use to cause disease in humans. Nowadays, Legionella species constitute a remarkable model of bacterial adaptation, with a genus genome shaped by their close coevolution with amoebae and an ability to exploit many hosts and signaling pathways through the secretion of a myriad of effector proteins, many of which have a eukaryotic origin. This review aims to discuss current knowledge of Legionella infection mechanisms and future research directions to be taken that might answer the many remaining open questions. This research will without a doubt be a terrific scientific journey worth taking.

A real-time PCR for specific detection of the Legionella pneumophila serogroup 1 ST1 complex

OBJECTIVE

Legionella pneumophila serogroup 1 (Lp1) sequence type (ST) 1 is globally widespread in the environment and accounts for a significant proportion of Legionella infections, including nosocomial Legionnaires' disease (LD). This study aimed to design a sensitive and specific detection method for Lp ST1 that will underpin epidemiological investigations and risk assessment.

METHODS

A total of 628 Lp genomes (126 ST1s) were analyzed by comparative genomics. Interrogation of more than 900 accessory genes revealed seven candidate targets for specific ST1 detection and specific primers and hydrolysis probes were designed and evaluated. The analytical sensitivity and specificity of the seven primer and probe sets were evaluated on serially diluted DNA extracted from the reference strain CIP107629 and via qPCR applied on 200 characterized isolates. The diagnostic performance of the assay was evaluated on 142 culture-proven clinical samples from LD cases and a real-life investigation of a case cluster.

RESULTS

Of seven qPCR assays that underwent analytical validation, one PCR target (lpp1868) showed higher sensitivity and specificity for ST1 and ST1-like strains. The diagnostic performance of the assay using respiratory samples corresponded to a sensitivity of 95% (19/20) (95% CI (75.1-99.9)) and specificity of 100% (122/122) (95% CI (97-100)). The ST1 PCR assay could link two out of three culture-negative hospitalized LD cases to ST1 during a real-time investigation.

CONCLUSION

Using whole genome sequencing (WGS) data, we developed and validated a sensitive and specific qPCR assay for the detection of Lp1 belonging to the ST1 clonal complex by amplification of the lpp1868 gene. The ST1 qPCR is expected to deliver an added value for Lp control and prevention, in conjunction with other recently developed molecular assays.

Diverse conjugative elements silence natural transformation in Legionella species

Natural transformation is a major mechanism of horizontal gene transfer. Although the genes required for natural transformation are nearly ubiquitous in bacteria, it is commonly reported that some isolates of transformable species fail to transform. In Legionella pneumophila, we show that the inability of multiple clinical isolates to transform is caused by a conjugative element that shuts down expression of genes required for transformation. Diverse conjugative elements in the Legionella genus have adopted the same inhibition strategy. We propose that inhibition of natural transformation by episomal and integrated conjugative elements can explain the lack of transformability of isolates and also the apparent lack of natural transformation in some species.

Natural transformation (i.e., the uptake of DNA and its stable integration in the chromosome) is a major mechanism of horizontal gene transfer in bacteria. Although the vast majority of bacterial genomes carry the genes involved in natural transformation, close relatives of naturally transformable species often appear not competent for natural transformation. In addition, unexplained extensive variations in the natural transformation phenotype have been reported in several species. Here, we addressed this phenomenon by conducting a genome-wide association study (GWAS) on a panel of isolates of the opportunistic pathogen Legionella pneumophila. GWAS revealed that the absence of the transformation phenotype is associated with the conjugative plasmid pLPL. The plasmid inhibits transformation by simultaneously silencing the genes required for DNA uptake and recombination. We identified a small RNA (sRNA), RocRp, as the sole plasmid-encoded factor responsible for the silencing of natural transformation. RocRp is homologous to the highly conserved and chromosome-encoded sRNA RocR which controls the transient expression of the DNA uptake system. Assisted by the ProQ/FinO-domain RNA chaperone RocC, RocRp acts as a substitute of RocR, ensuring that the bacterial host of the conjugative plasmid does not become naturally transformable. Distinct homologs of this plasmid-encoded sRNA are found in diverse conjugative elements in other Legionella species. Their low to high prevalence may result in the lack of transformability of some isolates up to the apparent absence of natural transformation in the species. Generally, our work suggests that conjugative elements obscure the widespread occurrence of natural transformability in bacteria.

Legionella: a reemerging pathogen

PURPOSE OF REVIEW

The present review summarizes new knowledge about Legionella epidemiology, clinical characteristics, community-associated and hospital-based outbreaks, molecular typing and molecular epidemiology, prevention, and detection in environmental and clinical specimens.

RECENT FINDINGS

The incidence of Legionnaire's disease is rising and the mortality rate remains high, particularly for immunocompromised patients. Extracorporeal membrane oxygenation may help support patients with severe respiratory failure. Fluoroquinolones and macrolides appear to be equally efficacious for treating Legionnaires' disease. Whole genome sequencing is an important tool for determining the source for Legionella infections and for understanding routes of transmission and mechanisms by which new pathogenic clones emerge. Real-time quantitative polymerase chain reaction testing of respiratory specimens may improve our ability to diagnose Legionnaire's disease. The frequency of viable but nonculturable organisms is quite high in some water systems but their role in causing clinical disease has not been defined.

SUMMARY

Legionellosis remains an important public health threat. To prevent these infections, staff of municipalities and large buildings must implement effective water system management programs that reduce Legionella growth and transmission and all Medicare-certified healthcare facilities must have water management policies. In addition, we need better methods for detecting Legionella in water systems and in clinical specimens to improve prevention strategies and clinical diagnosis.

Legionella pneumophila levels and sequence-type distribution in hospital hot water samples from faucets to connecting pipes

Recent studies have reported increased levels of Legionella pneumophila (Lp) at points of use compared to levels in primary and secondary components of hot water systems, suggesting possible selection by environmental conditions. In this study, concentrations of Lp in a hospital hot water system were evaluated by profile sampling, collecting successive water samples to determine the prevalence at the faucet (distal) and upstream piping before and after a system intervention to increase temperature. Lp strain diversity was compared between different points of use and different areas of the hot water system (i.e., tap, intermediate piping and main upflow piping). In total, 47 isolates were recovered from 32 positive hot water samples collected from designated taps, showers and recirculation loops; these isolates were subsequently analyzed by sequence-based typing (SBT). Lp levels were comparable between first draw (500 mL) and flushed (2 and 5 min) samples, whereas a decrease was observed in the amount of culturable cells (1 log). Two sequence types (STs) were identified throughout the system. ST378 (sg4/10) was present in 91% of samples, while ST154-like (sg1) was present in 41%; both STs were simultaneously recovered in 34% of samples. Isolated STs displayed comparable tolerance to copper (0.8-5 mg/L) and temperature (55 °C, 1 h) exposure. The ability to replicate within THP1 cells and Acanthamoeba castellanii was similar between the two STs and a comparative environmental outbreak strain. The low Lp diversity and the detection of both Lp sequence types in repeated subsequent samples collected from positive faucets in a hospital wing suggest a minimal impact of the distal conditions on strain selection for the sampled points, as well as a possible adaptation to stressors present in the system, leading to the predominance of a few strains.

Evaluation of SNP calling methods for closely related bacterial isolates and a novel high-accuracy pipeline: BactSNP

Bacteria are highly diverse, even within a species; thus, there have been many studies which classify a single species into multiple types and analyze the genetic differences between them. Recently, the use of whole-genome sequencing (WGS) has been popular for these analyses, and the identification of single-nucleotide polymorphisms (SNPs) between isolates is the most basic analysis performed following WGS. The performance of SNP-calling methods therefore has a significant effect on the accuracy of downstream analyses, such as phylogenetic tree inference. In particular, when closely related isolates are analyzed, e.g. in outbreak investigations, some SNP callers tend to detect a high number of false-positive SNPs compared with the limited number of true SNPs among isolates. However, the performances of various SNP callers in such a situation have not been validated sufficiently. Here, we show the results of realistic benchmarks of commonly used SNP callers, revealing that some of them exhibit markedly low accuracy when target isolates are closely related. As an alternative, we developed a novel pipeline BactSNP, which utilizes both assembly and mapping information and is capable of highly accurate and sensitive SNP calling in a single step. BactSNP is also able to call SNPs among isolates when the reference genome is a draft one or even when the user does not input the reference genome. BactSNP is available at https://github.com/IEkAdN/BactSNP.

Genomic heterogeneity differentiates clinical and environmental subgroups of Legionella pneumophila sequence type 1

Legionella spp. are the cause of a severe bacterial pneumonia known as Legionnaires' disease (LD). In some cases, current genetic subtyping methods cannot resolve LD outbreaks caused by common, potentially endemic L. pneumophila (Lp) sequence types (ST), which complicates laboratory investigations and environmental source attribution. In the United States (US), ST1 is the most prevalent clinical and environmental Lp sequence type. In order to characterize the ST1 population, we sequenced 289 outbreak and non-outbreak associated clinical and environmental ST1 and ST1-variant Lp strains from the US and, together with international isolate sequences, explored their genetic and geographic diversity. The ST1 population was highly conserved at the nucleotide level; 98% of core nucleotide positions were invariant and environmental isolates unassociated with human disease (n = 99) contained ~65% more nucleotide diversity compared to clinical-sporadic (n = 139) or outbreak-associated (n = 28) ST1 subgroups. The accessory pangenome of environmental isolates was also ~30-60% larger than other subgroups and was enriched for transposition and conjugative transfer-associated elements. Up to ~10% of US ST1 genetic variation could be explained by geographic origin, but considerable genetic conservation existed among strains isolated from geographically distant states and from different decades. These findings provide new insight into the ST1 population structure and establish a foundation for interpreting genetic relationships among ST1 strains; these data may also inform future analyses for improved outbreak investigations.

Spatial structuring of a Legionella pneumophila population within the water system of a large occupational building

The diversity of Legionella pneumophila populations within single water systems is not well understood, particularly in those unassociated with cases of Legionnaires’ disease. Here, we performed genomic analysis of 235 L. pneumophila isolates obtained from 28 water samples in 13 locations within a large occupational building. Despite regular treatment, the water system of this building is thought to have been colonized by L. pneumophila for at least 30 years without evidence of association with Legionnaires’ disease cases. All isolates belonged to one of three sequence types (STs), ST27 (n=81), ST68 (n=122) and ST87 (n=32), all three of which have been recovered from Legionnaires’ disease patients previously. Pairwise single nucleotide polymorphism differences amongst isolates of the same ST were low, ranging from 0 to 19 in ST27, from 0 to 30 in ST68 and from 0 to 7 in ST87, and no homologous recombination was observed in any lineage. However, there was evidence of horizontal transfer of a plasmid, which was found in all ST87 isolates and only one ST68 isolate. A single ST was found in 10/13 sampled locations, and isolates of each ST were also more similar to those from the same location compared with those from different locations, demonstrating spatial structuring of the population within the water system. These findings provide the first insights into the diversity and genomic evolution of a L. pneumophila population within a complex water system not associated with disease.

Legionella pneumophila and Other Legionella Species Isolated from Legionellosis Patients in Japan between 2008 and 2016

The Legionella Reference Center in Japan collected 427 Legionella clinical isolates between 2008 and 2016, including 7 representative isolates from corresponding outbreaks. The collection included 419 Legionella pneumophila isolates, of which 372 belonged to serogroup 1 (SG1) (87%) and the others belonged to SG2 to SG15 except for SG7 and SG11, and 8 isolates of other Legionella species (Legionella bozemanae, Legionella dumoffii, Legionella feeleii, Legionella longbeachae, Legionella londiniensis, and Legionella rubrilucens). L. pneumophila isolates were genotyped by sequence-based typing (SBT) and represented 187 sequence types (STs), of which 126 occurred in a single isolate (index of discrimination of 0.984). These STs were analyzed using minimum spanning tree analysis, resulting in the formation of 18 groups. The pattern of overall ST distribution among L. pneumophila isolates was diverse. In particular, some STs were frequently isolated and were suggested to be related to the infection sources. The major STs were ST23 (35 isolates), ST120 (20 isolates), and ST138 (16 isolates). ST23 was the most prevalent and most causative ST for outbreaks in Japan and Europe. ST138 has been observed only in Japan, where it has caused small-scale outbreaks; 81% of those strains (13 isolates) were suspected or confirmed to infect humans through bath water sources. On the other hand, 11 ST23 strains (31%) and 5 ST120 strains (25%) were suspected or confirmed to infect humans through bath water. These findings suggest that some ST strains frequently cause legionellosis in Japan and are found under different environmental conditions.IMPORTANCELegionella pneumophila serogroup 1 (SG1) is the most frequent cause of legionellosis. Our previous genetic analysis indicated that SG1 environmental isolates represented 8 major clonal complexes, consisting of 3 B groups, 2 C groups, and 3 S groups, which included major environmental isolates derived from bath water, cooling towers, and soil and puddles, respectively. Here, we surveyed clinical isolates collected from patients with legionellosis in Japan between 2008 and 2016. Most strains belonging to the B group were isolated from patients for whom bath water was the suspected or confirmed source of infection. Among the isolates derived from patients whose suspected infection source was soil or dust, most belonged to the S1 group and none belonged to the B or C groups. Additionally, the U group was discovered as a new group, which mainly included clinical isolates with unknown infection sources.

Clusters of Healthcare-Associated Legionnaires' Disease in Two Hospitals of Central Greece

Healthcare-associated Legionnaires' disease often leads to fatal respiratory tract infection among hospitalized patients. In this report, three cases of Legionnaires' disease among patients in two different hospitals (Hospital A and Hospital B) were investigated. After conducting an epidemiologic and environmental investigation, the water distribution systems (WDSs) were identified as the possible source of infection, as Legionella pneumophila serogroup 1 (Lp1) was isolated from both clinical and environmental samples. Patients received aerosol therapy with nebulizers during their hospitalization. Based on the results of the investigation, the hospitals' infection control committees reviewed their policies for Legionnaires' disease prevention and implemented control measures focusing on using sterile fluids for aerosol treatments.

Insights into the long-term persistence of Legionella in facilities from whole-genome sequencing

We investigated the value of whole-genome sequencing (WGS) and single nucleotide polymorphism (SNP) analyses in determining the relationships among and evolutionary rates of Legionella species with long-term persistence in three healthcare facilities. We examined retrospective clinical and environmental isolates of Legionella micdadei and Legionella pneumophila serogroup 1 isolates with identical PFGE DNA fingerprints sampled over the course of up to 18 years. WGS analyses demonstrated that heterogeneous populations of Legionella were present within each facility despite displaying the same PFGE profiles. Additionally, clustering of some clinical isolates with those from a separate but related institution exposed a source of infection not previously detected, underscoring the importance of considering phylogenetic relationships when assessing epidemiological links. The data supported an average substitution rate of 0.80 SNPs per genome per year for L. micdadei but a reliable estimate for L. pneumophila serogroup 1 could not be obtained due to complicating factors such as non-chronological links among isolates and inadequate sampling depths. While the substitution rate for L. micdadei is consistent with previous estimates for L. pneumophila, the lack of a temporal signal in our sequence data for L. pneuomphila serogroup 1 isolates suggests either insufficient change to provide an estimate or variable evolutionary rates, which could reflect the presence of both actively dividing and viable but non-culturable Legionella spp. in the built environment. This study highlights the increased discriminatory power of WGS SNP analysis as compared to PFGE, emphasizes the need for extended sampling, and provides insight into the evolution of Legionella from longitudinal investigations.

Epidemiological information is key when interpreting whole genome sequence data - lessons learned from a large Legionella pneumophila outbreak in Warstein, Germany, 2013

Whole genome sequencing (WGS) is increasingly used in Legionnaires’ disease (LD) outbreak investigations, owing to its higher resolution than sequence-based typing, the gold standard typing method for Legionella pneumophila, in the analysis of endemic strains. Recently, a gene-by-gene typing approach based on 1,521 core genes called core genome multilocus sequence typing (cgMLST) was described that enables a robust and standardised typing of L. pneumophila. Methods: We applied this cgMLST scheme to isolates obtained during the largest outbreak of LD reported so far in Germany. In this outbreak, the epidemic clone ST345 had been isolated from patients and four different environmental sources. In total 42 clinical and environmental isolates were retrospectively typed. Results: Epidemiologically unrelated ST345 isolates were clearly distinguishable from the epidemic clone. Remarkably, epidemic isolates split up into two distinct clusters, ST345-A and ST345-B, each respectively containing a mix of clinical and epidemiologically-related environmental samples. Discussion/conclusion: The outbreak was therefore likely caused by both variants of the single sequence type, which pre-existed in the environmental reservoirs. The two clusters differed by 40 alleles located in two neighbouring genomic regions of ca 42 and 26 kb. Additional analysis supported horizontal gene transfer of the two regions as responsible for the difference between the variants. Both regions comprise virulence genes and have previously been reported to be involved in recombination events. This corroborates the notion that genomic outbreak investigations should always take epidemiological information into consideration when making inferences. Overall, cgMLST proved helpful in disentangling the complex genomic epidemiology of the outbreak.

Differences in Virulence Between Legionella pneumophila Isolates From Human and Non-human Sources Determined in Galleria mellonella Infection Model

Legionella pneumophila is a ubiquitous bacterium in freshwater environments and in many man-made water systems capable of inducing pneumonia in humans. Despite its ubiquitous character most studies on L. pneumophila virulence focused on clinical strains and isolates from man-made environments, so little is known about the nature and extent of virulence variation in strains isolated from natural environments. It has been established that clinical isolates are less diverse than man-made and natural environmental strains, suggesting that only a subset of environmental isolates is specially adapted to infect humans. In this work we intended to determine if unrelated L. pneumophila strains, isolated from different environments and with distinct virulence-related genetic backgrounds, displayed differences in virulence, using the Wax Moth Galleria mellonella infection model. We found that all tested strains were pathogenic in G. mellonella, regardless of their origin. Indeed, a panoply of virulence-related phenotypes was observed sustaining the existence of significant differences on the ability of L. pneumophila strains to induce disease. Taken together our results suggest that the occurrence of human infection is not related with the increased capability of some strains to induce disease since we also found a concentration threshold above which L. pneumophila strains are equally able to cause disease. In addition, no link could be established between the sequence-type (ST) and L. pneumophila pathogenicity. We envision that in man-made water distribution systems environmental filtering selection and biotic competition acts structuring L. pneumophila populations by selecting more resilient and adapted strains that can rise to high concentration if no control measures are implemented. Therefore, public health strategies based on the sequence based typing (STB) scheme analysis should take into account that the major disease-associated clones of L. pneumophila were not related with higher virulence in G. mellonella infection model, and that potential variability of virulence-related phenotypes was found within the same ST.

Rolling epidemic of Legionnaires' disease outbreaks in small geographic areas

Legionnaires’ disease (LD) is reported from many parts of the world, mostly linked to drinking water sources or cooling towers. We reviewed two unusual rolling outbreaks in Sydney and New York, each clustered in time and space. Data on these outbreaks were collected from public sources and compared to previous outbreaks in Australia and the US. While recurrent outbreaks of LD over time linked to an identified single source have been described, multiple unrelated outbreaks clustered in time and geography have not been previously described. We describe unusual geographic and temporal clustering of Legionella outbreaks in two cities, each of which experienced multiple different outbreaks within a small geographic area and within a short timeframe. The explanation for this temporal and spatial clustering of LD outbreaks in two cities is not clear, but climate variation and deteriorating water sanitation are two possible explanations. There is a need to critically analyse LD outbreaks and better understand changing trends to effectively prevent disease.

Genome Sequencing Links Persistent Outbreak of Legionellosis in Sydney (New South Wales, Australia) to an Emerging Clone of Legionella pneumophila Sequence Type 211

The city of Sydney, Australia, experienced a persistent outbreak of Legionella pneumophila serogroup 1 (Lp1) pneumonia in 2016. To elucidate the source and guide public health actions, the genomes of clinical and environmental Lp1 isolates recovered over 7 weeks were examined. A total of 48 isolates from human cases and cooling towers were sequenced and compared using single-nucleotide polymorphism (SNP)-based core-genome multilocus sequencing typing (MLST) and pangenome approaches. All three methods confirmed phylogenetic relatedness between isolates associated with outbreaks in the Central Business District (CBD) in March and May and those in suburb 1. These isolates were designated the "main cluster" and consisted of isolates from two patients from the CBD March outbreak, one patient and one tower isolate from suburb 1, and isolates from two cooling towers and three patients from the CBD May outbreak. All main cluster isolates were sequence type 211 (ST211), which previously has only been reported in Canada. Significantly, pangenome analysis identified mobile genetic elements containing a unique type IV A F-type secretion system (T4ASS), which was specific to the main cluster, and cocirculating clinical strains, suggesting a potential mechanism for increased fitness and persistence of the outbreak clone. Genome sequencing enabled linking of the geographically dispersed environmental sources of infection among the spatially and temporally coinciding cases of legionellosis in a highly populated urban setting. The discovery of a unique T4ASS emphasizes the role of genome recombination in the emergence of successful Lp1 clones.IMPORTANCE A new emerging clone has been responsible for a prolonged legionellosis outbreak in Sydney, Australia. The use of whole-genome sequencing linked two outbreaks thought to be unrelated and confirmed the outliers. These findings led to the resampling and subsequent identification of the source, guiding public health actions and bringing the outbreak to a close. Significantly, the outbreak clone was identified as sequence type 211 (ST211). Our study reports this ST in the Southern Hemisphere and presents a description of ST211 genomes from both clinical and environmental isolates. A unique mobile genetic element containing a type IV secretion system was identified in Lp1 ST211 isolates linked to the main cluster and Lp1 ST42 isolates that were cocirculating at the time of the outbreak.

Prevalence of Infection-Competent Serogroup 6 Legionella pneumophila within Premise Plumbing in Southeast Michigan

Coinciding with major changes to its municipal water system, Flint, MI, endured Legionnaires’ disease outbreaks in 2014 and 2015. By sampling premise plumbing in Flint in the fall of 2016, we found that 12% of homes harbored legionellae, a frequency similar to that in residences in neighboring areas. To evaluate the genetic diversity of Legionella pneumophila in Southeast Michigan, we determined the sequence type (ST) and serogroup (SG) of the 18 residential isolates from Flint and Detroit, MI, and the 33 clinical isolates submitted by hospitals in three area counties in 2013 to 2016. Common to one environmental and four clinical samples were strains of L. pneumophila SG1 and ST1, the most prevalent ST worldwide. Among the Flint premise plumbing isolates, 14 of 16 strains were of ST367 and ST461, two closely related SG6 strain types isolated previously from patients and corresponding environmental samples. Each of the representative SG1 clinical strains and SG6 environmental isolates from Southeast Michigan infected and survived within macrophage cultures at least as well as a virulent laboratory strain, as judged by microscopy and by enumerating CFU. Likewise, 72 h after infection, the yield of viable-cell counts increased >100-fold for each of the representative SG1 clinical isolates, Flint premise plumbing SG6 ST367 and -461 isolates, and two Detroit residential isolates. We verified by immunostaining that SG1-specific antibody does not cross-react with the SG6 L. pneumophila environmental strains. Because the widely used urinary antigen diagnostic test does not readily detect non-SG1 L. pneumophila, Legionnaires’ disease caused by SG6 L. pneumophila is likely underreported worldwide.

L. pneumophila is the leading cause of disease outbreaks associated with drinking water in the United States. Compared to what is known of the established risks of colonization within hospitals and hotels, relatively little is known about residential exposure to L. pneumophila. One year after two outbreaks of Legionnaires’ disease in Genesee County, MI, that coincided with damage to the Flint municipal water system, our multidisciplinary team launched an environmental surveillance and laboratory research campaign aimed at informing risk management strategies to provide safe public water supplies. The most prevalent L. pneumophila strains isolated from residential plumbing were closely related strains of SG6. In laboratory tests of virulence, the SG6 environmental isolates resembled SG1 clinical strains, yet they are not readily detected by the common diagnostic urinary antigen test, which is specific for SG1. Therefore, our study complements the existing epidemiological literature indicating that Legionnaires’ disease due to non-SG1 strains is underreported around the globe.

Seeding and Establishment of Legionella pneumophila in Hospitals: Implications for Genomic Investigations of Nosocomial Legionnaires' Disease

Background.

Legionnaires’ disease is an important cause of hospital-acquired pneumonia and is caused by infection with the bacterium Legionella. Because current typing methods often fail to resolve the infection source in possible nosocomial cases, we aimed to determine whether whole-genome sequencing (WGS) could be used to support or refute suspected links between cases and hospitals. We focused on cases involving a major nosocomial-associated strain, L. pneumophila sequence type (ST) 1.

Methods.

WGS data from 229 L. pneumophila ST1 isolates were analyzed, including 99 isolates from the water systems of 17 hospitals and 42 clinical isolates from patients with confirmed or suspected hospital-acquired infections, as well as isolates obtained from or associated with community-acquired sources of Legionnaires’ disease.

Results.

Phylogenetic analysis demonstrated that all hospitals from which multiple isolates were obtained have been colonized by 1 or more distinct ST1 populations. However, deep sampling of 1 hospital also revealed the existence of substantial diversity and ward-specific microevolution within the population. Across all hospitals, suspected links with cases were supported with WGS, although the degree of support was dependent on the depth of environmental sampling and available contextual information. Finally, phylogeographic analysis revealed that hospitals have been seeded with L. pneumophila via both local and international spread of ST1.

Conclusions.

WGS can be used to support or refute suspected links between hospitals and Legionnaires’ disease cases. However, deep hospital sampling is frequently required due to the potential coexistence of multiple populations, existence of substantial diversity, and similarity of hospital isolates to local populations.

Genomic investigation of a suspected outbreak of Legionella pneumophila ST82 reveals undetected heterogeneity by the present gold-standard methods, Denmark, July to November 2014

Between July and November 2014, 15 community-acquired cases of Legionnaires´ disease (LD), including four with Legionella pneumophila serogroup 1 sequence type (ST) 82, were diagnosed in Northern Zealand, Denmark. An outbreak was suspected. No ST82 isolates were found in environmental samples and no external source was established. Four putative-outbreak ST82 isolates were retrospectively subjected to whole genome sequencing (WGS) followed by phylogenetic analyses with epidemiologically unrelated ST82 sequences. The four putative-outbreak ST82 sequences fell into two clades, the two clades were separated by ca 1,700 single nt polymorphisms (SNP)s when recombination regions were included but only by 12 to 21 SNPs when these were removed. A single putative-outbreak ST82 isolate sequence segregated in the first clade. The other three clustered in the second clade, where all included sequences had < 5 SNP differences between them. Intriguingly, this clade also comprised epidemiologically unrelated isolate sequences from the UK and Denmark dating back as early as 2011. The study confirms that recombination plays a major role in L. pneumophila evolution. On the other hand, strains belonging to the same ST can have only few SNP differences despite being sampled over both large timespans and geographic distances. These are two important factors to consider in outbreak investigations.