Comparative genome analysis reveals a complex population structure of Legionella pneumophila subspecies

The majority of Legionnaires' disease (LD) cases are caused by Legionella pneumophila, a genetically heterogeneous species composed of at least 17 serogroups. Previously, it was demonstrated that L. pneumophila consists of three subspecies: pneumophila, fraseri and pascullei. During an LD outbreak investigation in 2012, we detected that representatives of both subspecies fraseri and pascullei colonized the same water system and that the outbreak-causing strain was a new member of the least represented subspecies pascullei. We used partial sequence based typing consensus patterns to mine an international database for additional representatives of fraseri and pascullei subspecies. As a result, we identified 46 sequence types (STs) belonging to subspecies fraseri and two STs belonging to subspecies pascullei. Moreover, a recent retrospective whole genome sequencing analysis of isolates from New York State LD clusters revealed the presence of a fourth L. pneumophila subspecies that we have termed raphaeli. This subspecies consists of 15 STs. Comparative analysis was conducted using the genomes of multiple members of all four L. pneumophila subspecies. Whereas each subspecies forms a distinct phylogenetic clade within the L. pneumophila species, they share more average nucleotide identity with each other than with other Legionella species. Unique genes for each subspecies were identified and could be used for rapid subspecies detection. Improved taxonomic classification of L. pneumophila strains may help identify environmental niches and virulence attributes associated with these genetically distinct subspecies.

Streptococcus Mutans Adhesin Biotypes that Match and Predict Individual Caries Development

Dental caries, which affects billions of people, is a chronic infectious disease that involves Streptococcus mutans, which is nevertheless a poor predictor of individual caries development. We therefore investigated if adhesin types of S.mutans with sucrose-independent adhesion to host DMBT1 (i.e. SpaP A, B or C) and collagen (i.e. Cnm, Cbm) match and predict individual differences in caries development. The adhesin types were measured in whole saliva by qPCR in 452 12-year-old Swedish children and related to caries at baseline and prospectively at a 5-year follow-up. Strains isolated from the children were explored for genetic and phenotypic properties. The presence of SpaP B and Cnm subtypes coincided with increased 5-year caries increment, and their binding to DMBT1 and saliva correlated with individual caries scores. The SpaP B subtypes are enriched in amino acid substitutions that coincided with caries and binding and specify biotypes of S. mutans with increased acid tolerance. The findings reveal adhesin subtypes of S. mutans that match and predict individual differences in caries development and provide a rationale for individualized oral care.

•

Adhesin subtypes of Streptococcus mutans match and predict individual caries development.

•

Adhesin binding to salivary DMBT1 correlates with individual caries scores.

•

The adhesin types coincide with distinct biotypes of S. mutans.

Dental caries, which affects billions of people, involves the bacterium Streptococcus mutans, which is nevertheless a poor predictor of caries development. The present findings provide the first evidence that S. mutans adhesin subtypes match and predict individual 5-year caries development in Swedish children. The binding strength of the adhesin subtypes correlates with individual caries scores, and the adhesin subtypes specify biotypes of S. mutans that also differ in acid tolerance. The present findings provide a rationale for individualized oral care and improved systemic health because chronic caries infection and carrying high-virulence strains pose a systemic disease risk.

Is Predominant Clonal Evolution a Common Evolutionary Adaptation to Parasitism in Pathogenic Parasitic Protozoa, Fungi, Bacteria, and Viruses?

We propose that predominant clonal evolution (PCE) in microbial pathogens be defined as restrained recombination on an evolutionary scale, with genetic exchange scarce enough to not break the prevalent pattern of clonal population structure. The main features of PCE are (1) strong linkage disequilibrium, (2) the widespread occurrence of stable genetic clusters blurred by occasional bouts of genetic exchange ('near-clades'), (3) the existence of a "clonality threshold", beyond which recombination is efficiently countered by PCE, and near-clades irreversibly diverge. We hypothesize that the PCE features are not mainly due to natural selection but also chiefly originate from in-built genetic properties of pathogens. We show that the PCE model obtains even in microbes that have been considered as 'highly recombining', such as Neisseria meningitidis, and that some clonality features are observed even in Plasmodium, which has been long described as panmictic. Lastly, we provide evidence that PCE features are also observed in viruses, taking into account their extremely fast genetic turnover. The PCE model provides a convenient population genetic framework for any kind of micropathogen. It makes it possible to describe convenient units of analysis (clones and near-clades) for all applied studies. Due to PCE features, these units of analysis are stable in space and time, and clearly delimited. The PCE model opens up the possibility of revisiting the problem of species definition in these organisms. We hypothesize that PCE constitutes a major evolutionary strategy for protozoa, fungi, bacteria, and viruses to adapt to parasitism.

The life stage-specific pathometabolism of Legionella pneumophila

The genus Legionella belongs to Gram-negative bacteria found ubiquitously in aquatic habitats, where it grows in natural biofilms and replicates intracellularly in various protozoa (amoebae, ciliates). L. pneumophila is known as the causative agent of Legionnaires' disease, since it is also able to replicate in human alveolar macrophages, finally leading to inflammation of the lung and pneumonia. To withstand the degradation by its host cells, a Legionella-containing vacuole (LCV) is established for intracellular replication, and numerous effector proteins are secreted into the host cytosol using a type four B secretion system (T4BSS). During intracellular replication, Legionella has a biphasic developmental cycle that alternates between a replicative and a transmissive form. New knowledge about the host-adapted and life stage-dependent metabolism of intracellular L. pneumophila revealed a bipartite metabolic network with life stage-specific usages of amino acids (e.g. serine), carbohydrates (e.g. glucose) and glycerol as major substrates. These metabolic features are associated with the differentiation of the intracellular bacteria, and thus have an important impact on the virulence of L. pneumophila.

Molecular epidemiology, phylogeny and evolution of Legionella

Legionella are opportunistic pathogens that develop in aquatic environments where they multiply in protozoa. When infected aerosols reach the human respiratory tract they may accidentally infect the alveolar macrophages leading to a severe pneumonia called Legionnaires' disease (LD). The ability of Legionella to survive within host-cells is strictly dependent on the Dot/Icm Type 4 Secretion System that translocates a large repertoire of effectors into the host cell cytosol. Although Legionella is a large genus comprising nearly 60 species that are worldwide distributed, only about half of them have been involved in LD cases. Strikingly, the species Legionella pneumophila alone is responsible for 90% of all LD cases. The present review summarizes the molecular approaches that are used for L. pneumophila genotyping with a major focus on the contribution of whole genome sequencing (WGS) to the investigation of local L. pneumophila outbreaks and global epidemiology studies. We report the newest knowledge regarding the phylogeny and the evolution of Legionella and then focus on virulence evolution of those Legionella species that are known to have the capacity to infect humans. Finally, we discuss the evolutionary forces and adaptation mechanisms acting on the Dot/Icm system itself as well as the role of mobile genetic elements (MGE) encoding T4ASSs and of gene duplications in the evolution of Legionella and its adaptation to different hosts and lifestyles.

Genetic diversity and evolutionary relationships among Legionella pneumophila clinical isolates, Portugal, 1987 to 2012

The genetic diversity of 89 clinical Legionella isolates, collected between 1987 and 2012, in 22 hospitals from the five regions of Portugal, was analysed in this study using monoclonal antibodies (MAbs) of the Dresden panel and the sequence-based typing (SBT) protocol. The eBURST algorithm was used to infer levels of relatedness between isolates. All isolates collected were Legionella pneumophila, which were further characterised into four subgroups by MAbs, and 30 sequence types (STs) by SBT. Twelve of the STs were unique to Portugal; one of them (ST100) was represented by 32 epidemiologically related isolates. The ST44 was the profile with the highest number of epidemiologically unrelated isolates. The eBURST analyses indicate that, within the group formed by the 30 STs identified in this study, 17 STs were genetically close to at least another ST in the group. The comparison between the eBURST diagrams obtained with the STs from this study and the entire SBT database of the European Working Group for Legionella, showed that 24 (seven of them unique to Portugal) of our 30 STs were related with STs identified in others countries. These results suggest that the population of L. pneumophila clinical strains in Portugal includes both worldwide and local strains.

Intragenic recombination has a critical role on the evolution of Legionella pneumophila virulence-related effector sidJ

SidJ is a Dot/Icm effector involved in the trafficking or retention of ER-derived vesicles to Legionella pneumophila vacuoles whose mutation causes an observable growth defect, both in macrophage and amoeba hosts. Given the crucial role of this effector in L. pneumophila virulence we investigated the mechanisms shaping its molecular evolution. The alignment of SidJ sequences revealed several alleles with amino acid variations that may influence the protein properties. The identification of HGT events and the detection of balancing selection operating on sidJ evolution emerge as a clear result. Evidence suggests that intragenic recombination is an important strategy in the evolutionary adaptive process playing an active role on sidJ genetic plasticity. This pattern of evolution is in accordance with the life style of L. pneumophila as a broad host-range pathogen by preventing host-specialization and contributing to the resilience of the species.

Real-time investigation of a Legionella pneumophila outbreak using whole genome sequencing

Legionella pneumophila is the main pathogen responsible for outbreaks of Legionnaires' disease, which can be related to contaminated water supplies such as cooling towers or water pipes. We combined conventional molecular methods and whole genome sequence (WGS) analysis to investigate an outbreak of L. pneumophila in a large Australian hospital. Typing of these isolates using sequence-based typing and virulence gene profiling, was unable to discriminate between outbreak and non-outbreak isolates. WGS analysis was performed on isolates during the outbreak, as well as on unlinked isolates from the Public Health Microbiology reference collection. The more powerful resolution provided by analysis of whole genome sequences allowed outbreak isolates to be distinguished from isolates that were temporally and spatially unassociated with the outbreak, demonstrating that this technology can be used in real-time to investigate L. pneumophila outbreaks.

Comparison of the Legionella pneumophila population structure as determined by sequence-based typing and whole genome sequencing

Background

Legionella pneumophila is an opportunistic pathogen of humans where the source of infection is usually from contaminated man-made water systems. When an outbreak of Legionnaires’ disease caused by L. pneumophila occurs, it is necessary to discover the source of infection. A seven allele sequence-based typing scheme (SBT) has been very successful in providing the means to attribute outbreaks of L. pneumophila to a particular source or sources. Particular sequence types described by this scheme are known to exhibit specific phenotypes. For instance some types are seen often in clinical cases but are rarely isolated from the environment and vice versa. Of those causing human disease some types are thought to be more likely to cause more severe disease. It is possible that the genetic basis for these differences are vertically inherited and associated with particular genetic lineages within the population. In order to provide a framework within which to test this hypothesis and others relating to the population biology of L. pneumophila, a set of genomes covering the known diversity of the organism is required.

Results

Firstly, this study describes a means to group L. pneumophila strains into pragmatic clusters, using a methodology that takes into consideration the genetic forces operating on the population. These clusters can be used as a standardised nomenclature, so those wishing to describe a group of strains can do so. Secondly, the clusters generated from the first part of the study were used to select strains rationally for whole genome sequencing (WGS). The data generated was used to compare phylogenies derived from SBT and WGS. In general the SBT sequence type (ST) accurately reflects the whole genome-based genotype. Where there are exceptions and recombination has resulted in the ST no longer reflecting the genetic lineage described by the whole genome sequence, the clustering technique employed detects these sequence types as being admixed, indicating their mixed inheritance.

Conclusions

We conclude that SBT is usually a good proxy for the genetic lineage described by the whole genome, and therefore utility of SBT is still suitable until the technology and economics of high throughput sequencing reach the point where routine WGS of L. pneumophila isolates for outbreak investigation is feasible.

Reproductive clonality of pathogens: a perspective on pathogenic viruses, bacteria, fungi, and parasitic protozoa

We propose that clonal evolution in micropathogens be defined as restrained recombination on an evolutionary scale, with genetic exchange scarce enough to not break the prevalent pattern of clonal population structure, a definition already widely used for all kinds of pathogens, although not clearly formulated by many scientists and rejected by others. The two main manifestations of clonal evolution are strong linkage disequilibrium (LD) and widespread genetic clustering ("near-clading"). We hypothesize that this pattern is not mainly due to natural selection, but originates chiefly from in-built genetic properties of pathogens, which could be ancestral and could function as alternative allelic systems to recombination genes ("clonality/sexuality machinery") to escape recombinational load. The clonal framework of species of pathogens should be ascertained before any analysis of biomedical phenotypes (phylogenetic character mapping). In our opinion, this model provides a conceptual framework for the population genetics of any micropathogen.

Extensive recombination events and horizontal gene transfer shaped the Legionella pneumophila genomes

Background

Legionella pneumophila is an intracellular pathogen of environmental protozoa. When humans inhale contaminated aerosols this bacterium may cause a severe pneumonia called Legionnaires' disease. Despite the abundance of dozens of Legionella species in aquatic reservoirs, the vast majority of human disease is caused by a single serogroup (Sg) of a single species, namely L. pneumophila Sg1. To get further insights into genome dynamics and evolution of Sg1 strains, we sequenced strains Lorraine and HL 0604 1035 (Sg1) and compared them to the available sequences of Sg1 strains Paris, Lens, Corby and Philadelphia, resulting in a comprehensive multigenome analysis.

Results

We show that L. pneumophila Sg1 has a highly conserved and syntenic core genome that comprises the many eukaryotic like proteins and a conserved repertoire of over 200 Dot/Icm type IV secreted substrates. However, recombination events and horizontal gene transfer are frequent. In particular the analyses of the distribution of nucleotide polymorphisms suggests that large chromosomal fragments of over 200 kbs are exchanged between L. pneumophila strains and contribute to the genome dynamics in the natural population. The many secretion systems present might be implicated in exchange of these fragments by conjugal transfer. Plasmids also play a role in genome diversification and are exchanged among strains and circulate between different Legionella species.

Conclusion

Horizontal gene transfer among bacteria and from eukaryotes to L. pneumophila as well as recombination between strains allows different clones to evolve into predominant disease clones and others to replace them subsequently within relatively short periods of time.

High-throughput typing method to identify a non-outbreak-involved Legionella pneumophila strain colonizing the entire water supply system in the town of Rennes, France

Two legionellosis outbreaks occurred in the city of Rennes, France, during the past decade, requiring in-depth monitoring of Legionella pneumophila in the water network and the cooling towers in the city. In order to characterize the resulting large collection of isolates, an automated low-cost typing method was developed. The multiplex capillary-based variable-number tandem repeat (VNTR) (multiple-locus VNTR analysis [MLVA]) assay requiring only one PCR amplification per isolate ensures a high level of discrimination and reduces hands-on and time requirements. In less than 2 days and using one 4-capillary apparatus, 217 environmental isolates collected between 2000 and 2009 and 5 clinical isolates obtained during outbreaks in 2000 and 2006 in Rennes were analyzed, and 15 different genotypes were identified. A large cluster of isolates with closely related genotypes and representing 77% of the population was composed exclusively of environmental isolates extracted from hot water supply systems. It was not responsible for the known Rennes epidemic cases, although strains showing a similar MLVA profile have regularly been involved in European outbreaks. The clinical isolates in Rennes had the same genotype as isolates contaminating a mall's cooling tower. This study further demonstrates that unknown environmental or genetic factors contribute to the pathogenicity of some strains. This work illustrates the potential of the high-throughput MLVA typing method to investigate the origin of legionellosis cases by allowing the systematic typing of any new isolate and inclusion of data in shared databases.

Investigation of the population structure of Legionella pneumophila by analysis of tandem repeat copy number and internal sequence variation

The population structure of the species Legionella pneumophila was investigated by multilocus variable number of tandem repeats (VNTR) analysis (MLVA) and sequencing of three VNTRs (Lpms01, Lpms04 and Lpms13) in selected strains. Of 150 isolates of diverse origins, 136 (86 %) were distributed into eight large MLVA clonal complexes (VACCs) and the rest were either unique or formed small clusters of up to two MLVA genotypes. In spite of the lower degree of genome-wide linkage disequilibrium of the MLVA loci compared with sequence-based typing, the clustering achieved by the two methods was highly congruent. The detailed analysis of VNTR Lpms04 alleles showed a very complex organization, with five different repeat unit lengths and a high level of internal variation. Within each MLVA-defined VACC, Lpms04 was endowed with a common recognizable pattern with some interesting exceptions. Evidence of recombination events was suggested by analysis of internal repeat variations at the two additional VNTR loci, Lpms01 and Lpms13. Sequence analysis of L. pneumophila VNTR locus Lpms04 alone provides a first-line assay for allocation of a new isolate within the L. pneumophila population structure and for epidemiological studies.

Molecular evolution of Legionella pneumophila dotA gene, the contribution of natural environmental strains

Given the role of DotA protein in establishing successful infections and the diversity of host cells interacting with Legionella pneumophila in nature, it is possible that this gene product is a target for adaptive evolution. We investigated the influence of L. pneumophila isolates from natural environments with the molecular evolution of this crucial virulence-related gene. The population genetic structure of L. pneumophila was inferred from the partial sequences of rpoB and dotA of 303 worldwide strains. The topology of the two inferred trees was not congruent and in the inferred dotA tree the vast majority of the natural environmental isolates were clustered in a discrete group. The Ka/Ks ratio demonstrated that this group, contrary to all others, has been under strong diversifying selection. The alignment of all DotA sequences allowed the identification of several alleles and the amino acid variations were not randomly distributed. Moreover, from these results we can conclude that dotA from L. pneumophila clinical and man-made environmental strains belong to a sub-set of all genotypes existing in nature. A split graph analysis showed evidence of a network-like organization and several intergenic recombination events were detected within L. pneumophila strains resulting in mosaic genes in which different gene segments exhibited different evolutionary histories. We have determined that the allelic diversity of dotA is predominantly found in L. pneumophila isolates from natural environments, suggesting that niche-specific selection pressures have been operating on this gene. Indeed, the high level of dotA allelic diversity may reflect fitness variation in the persistence of those strains in distinct environmental niches and/or tropism to various protozoan hosts.

Molecular pathogenesis of infections caused by Legionella pneumophila

The genus Legionella contains more than 50 species, of which at least 24 have been associated with human infection. The best-characterized member of the genus, Legionella pneumophila, is the major causative agent of Legionnaires' disease, a severe form of acute pneumonia. L. pneumophila is an intracellular pathogen, and as part of its pathogenesis, the bacteria avoid phagolysosome fusion and replicate within alveolar macrophages and epithelial cells in a vacuole that exhibits many characteristics of the endoplasmic reticulum (ER). The formation of the unusual L. pneumophila vacuole is a feature of its interaction with the host, yet the mechanisms by which the bacteria avoid classical endosome fusion and recruit markers of the ER are incompletely understood. Here we review the factors that contribute to the ability of L. pneumophila to infect and replicate in human cells and amoebae with an emphasis on proteins that are secreted by the bacteria into the Legionella vacuole and/or the host cell. Many of these factors undermine eukaryotic trafficking and signaling pathways by acting as functional and, in some cases, structural mimics of eukaryotic proteins. We discuss the consequences of this mimicry for the biology of the infected cell and also for immune responses to L. pneumophila infection.

New endemic Legionella pneumophila serogroup I clones, Ontario, Canada

Identifying geographic distribution can improve surveillance and clinical testing procedures.

The water-borne pathogen Legionella pneumophila serogroup 1 (Lp1) is the most commonly reported etiologic agent of legionellosis. To examine the genetic diversity, the long-term epidemiology, and the molecular evolution of Lp1 clinical isolates, we conducted sequence-based typing on a collection of clinical isolates representing 3 decades of culture-confirmed legionellosis in Ontario, Canada. Analysis showed that the population of Lp1 in Ontario is highly diverse and combines lineages identified worldwide with local strains. Identical types were identified in sporadic and outbreak-associated strains. In the past 15 years, the incidence of some lineages distributed worldwide has tended to decrease, and local endemic clones and lineages have emerged. Comparative geographic distribution analysis suggests that some lineages are specific to eastern North America. These findings have general clinical implications for the study of Lp1 molecular evolution and for the identification of Lp1 circulating strains in North America.

Generation of diversity in Streptococcus mutans genes demonstrated by MLST

Streptococcus mutans, consisting of serotypes c, e, f and k, is an oral aciduric organism associated with the initiation and progression of dental caries. A total of 135 independent Streptococcus mutans strains from caries-free and caries-active subjects isolated from various geographical locations were examined in two versions of an MLST scheme consisting of either 6 housekeeping genes [accC (acetyl-CoA carboxylase biotin carboxylase subunit), gki (glucokinase), lepA (GTP-binding protein), recP (transketolase), sodA (superoxide dismutase), and tyrS (tyrosyl-tRNA synthetase)] or the housekeeping genes supplemented with 2 extracellular putative virulence genes [gtfB (glucosyltransferase B) and spaP (surface protein antigen I/II)] to increase sequence type diversity. The number of alleles found varied between 20 (lepA) and 37 (spaP). Overall, 121 sequence types (STs) were defined using the housekeeping genes alone and 122 with all genes. However pi, nucleotide diversity per site, was low for all loci being in the range 0.019-0.007. The virulence genes exhibited the greatest nucleotide diversity and the recombination/mutation ratio was 0.67 [95% confidence interval 0.3-1.15] compared to 8.3 [95% confidence interval 5.0-14.5] for the 6 concatenated housekeeping genes alone. The ML trees generated for individual MLST loci were significantly incongruent and not significantly different from random trees. Analysis using ClonalFrame indicated that the majority of isolates were singletons and no evidence for a clonal structure or evidence to support serotype c strains as the ancestral S. mutans strain was apparent. There was also no evidence of a geographical distribution of individual isolates or that particular isolate clusters were associated with caries. The overall low sequence diversity suggests that S. mutans is a newly emerged species which has not accumulated large numbers of mutations but those that have occurred have been shuffled as a consequence of intra-species recombination generating genotypes which can be readily distinguished by sequence analysis.

Distribution of lag-1 alleles and sequence-based types among Legionella pneumophila serogroup 1 clinical and environmental isolates in the United States

Approximately 84% of legionellosis cases are due to Legionella pneumophila serogroup 1. Moreover, a majority of L. pneumophila serogroup 1 clinical isolates react positively with monoclonal antibody 2 (MAb2) of the international standard panel. Over 94% of the legionellosis outbreaks investigated by the Centers for Disease Control and Prevention are due to this subset of L. pneumophila serogroup 1. To date, there is no complete explanation for the enhanced ability of these strains to cause disease. To better characterize these organisms, we subtyped 100 clinical L. pneumophila serogroup 1 isolates and 50 environmental L. pneumophila serogroup 1 isolates from the United States by (i) reactivity with MAb2, (ii) presence of a lag-1 gene required for the MAb2 epitope, and (iii) sequence-based typing analysis. Our results showed that the MAb2 epitope and lag-1 gene are overrepresented in clinical L. pneumophila serogroup 1 isolates. MAb2 recognized 75% of clinical isolates but only 6% of environmental isolates. Similarly, 75% of clinical isolates but only 8% of environmental isolates harbored lag-1. We identified three distinct lag-1 alleles, referred to as Philadelphia, Arizona, and Lens alleles, among 79 isolates carrying this gene. The Arizona allele is described for the first time in this study. We identified 59 different sequence types (STs), and 34 STs (58%) were unique to the United States. Our results support the hypothesis that a select group of STs may have an enhanced ability to cause legionellosis. Combining sequence typing and lag-1 analysis shows that STs tend to associate with a single lag-1 allele type, suggesting a hierarchy of virulence genotypes. Further analysis of ST and lag-1 profiles may identify genotypes of L. pneumophila serogroup 1 that warrant immediate intervention.

Direct sequencing of Legionella pneumophila from respiratory samples for sequence-based typing analysis

We have developed a procedure to test the efficiency and reliability of sequencing of Legionella pneumophila genes directly from respiratory samples and have compared the results with those derived from cultured isolates. We tried to obtain the nucleotide sequences of six protein-coding loci included in the sequence-based typing scheme for Legionella pneumophila and three intergenic regions from 132 samples corresponding to 106 patients positive for urine antigen. A seminested PCR approach was used to amplify and sequence these nine loci directly from respiratory samples. Nucleotide sequences were directly obtained for 23 Legionella isolates and also for 66 respiratory secretions from a total of 69 patients. The efficiency of sequencing from respiratory secretions was higher than that of sequencing after the isolation of the Legionella isolates. Moreover, the perfect match between the sequences obtained by both approaches when respiratory samples and cultured isolates from the same patient were available corroborates the suitability of the direct sequencing approach for the identification of Legionella species and molecular epidemiology studies with Legionella species.