Genetic variability in environmental isolates of Legionella pneumophila from Comunidad Valenciana (Spain)

Prevalence of Legionella species isolated from shower water in public bath facilities in Toyama Prefecture, Japan

AIMS

We investigated the prevalence of Legionella spp. isolated from shower water in public bath facilities in Toyama Prefecture, Japan. In addition, we analyzed the genetic diversity among Legionella pneumophila isolates from shower water as well as the genetic relationship between isolates from shower water and from stock strains previously analyzed from sputum specimens.

METHODS

The isolates were characterized using serogrouping, 16S rRNA gene sequencing, and sequence-based typing.

RESULTS

Legionella spp. were isolated from 31/91 (34.1%) samples derived from 17/37 (45.9%) bath facilities. Isolates from shower water and bath water in each public bath facility were serologically or genetically different, indicating that we need to isolate several L. pneumophila colonies from both bath and shower water to identify public bath facilities as sources of legionellosis. The 61 L. pneumophila isolates from shower water were classified into 39 sequence types (STs) (index of discrimination = 0.974), including 19 new STs. Among the 39 STs, 12 STs match clinical isolates in the European Working Group for Legionella Infections database. Notably, ST505 L. pneumophila SG 1, a strain frequently isolated from patients with legionellosis and from bath water in this area, was isolated from shower water.

CONCLUSIONS

Pathogenic L. pneumophila strains including ST505 strain were widely distributed in shower water in public bath facilities, with genetic diversity showing several different origins. This study highlights the need to isolate several L. pneumophila colonies from both bath water and shower water to identify public bath facilities as infection sources in legionellosis cases.

Geographical and Temporal Structures of Legionella pneumophila Sequence Types in Comunitat Valenciana (Spain), 1998 to 2013

Legionella pneumophila is an accidental human pathogen associated with aerosol formation in water-related sources. High recombination rates make Legionella populations genetically diverse, and nearly 2,000 different sequence types (STs) have been described to date for this environmental pathogen. The spatial distribution of STs is extremely heterogeneous, with some variants being present worldwide and others being detected at only a local scale. Similarly, some STs have been associated with disease outbreaks, such as ST578 or ST23. Spain is among the European countries with the highest incidences of reported legionellosis cases, and specifically, Comunitat Valenciana (CV) is the second most affected area in the country. In this work, we aimed at studying the overall diversity of Legionella pneumophila populations found in the period from 1998 to 2013 in 79 localities encompassing 23 regions within CV. To do so, we performed sequence-based typing (SBT) on 1,088 L. pneumophila strains detected in the area from both environmental and clinical sources. A comparison with the genetic structuring detected in a global data set that included 20 European and 7 non-European countries was performed. Our results reveal a level of diversity in CV that can be considered representative of the diversity found in other countries worldwide.

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.

Tollip or not Tollip: what are the evolving questions behind it?

Tollip plays an important role in the interleukin-1 receptor IL-1R and Toll pathways. As a modulator of the immune pathway, it indirectly controls the amount of antimicrobial peptides. This could indicate a vital step in maintaining animal immune systems and preventing infection. Evolutionary questions are crucial to understanding the conservation and functioning of the biochemical pathways like the Tollip-mediated one. Through an analysis of 36 sequences of the Tollip protein from different animal taxa, downloaded from Kyoto Encyclopedia of Genes and Genomes (KEGG) databank, we inferred diverse evolutionary parameters, such as molecular selection and structure conservation, by analyzing residue by residue, beyond the canonical parameters to this type of study, as maximum likelihood trees. We found that Tollip presented different trends in its evolving history. In primates, the protein is becoming more unstable, just the opposite is observed in the arthropod group. The most interesting finding was the concentration of positively selected residues at amino terminal ends. Some observed topological incongruences in maximum likelihood trees of complete and curated Tollip data sets could be explained through horizontal transfers, evidenced by recombination detection. These results suggest that there is more to be researched and understood about this protein.

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.

Genetic Characterization of Legionella pneumophila Isolated from a Common Watershed in Comunidad Valenciana, Spain

Legionella pneumophila infects humans to produce legionellosis and Pontiac fever only from environmental sources. In order to establish control measures and study the sources of outbreaks it is essential to know extent and distribution of strain variants of this bacterium in the environment. Sporadic and outbreak-related cases of legionellosis have been historically frequent in the Comunidad Valenciana region (CV, Spain), with a high prevalence in its Southeastern-most part (BV). Environmental investigations for the detection of Legionella pneumophila are performed in this area routinely. We present a population genetics study of 87 L. pneumophila strains isolated in 13 different localities of the BV area irrigated from the same watershed and compare them to a dataset of 46 strains isolated in different points of the whole CV. Our goal was to compare environmental genetic variation at two different geographic scales, at county and regional levels. Genetic diversity, recombination and population structure were analyzed with Sequence-Based Typing data and three intergenic regions. The results obtained reveal a low, but detectable, level of genetic differentiation between both datasets, mainly, but not only, attributed to the occurrence of unusual variants of the neuA locus present in the BV populations. This differentiation is still detectable when the 10 loci considered are analyzed independently, despite the relatively high incidence of the most common genetic variant in this species, sequence type 1 (ST-1). However, when the genetic data are considered without their associated geographic information, four major groups could be inferred at the genetic level which did not show any correlation with sampling locations. The overall results indicate that the population structure of these environmental samples results from the joint action of a global, widespread ST-1 along with genetic differentiation at shorter geographic distances, which in this case are related to the common watershed for the BV localities.

Patterns of nucleotide diversity of the ldpA circadian gene in closely related species of cyanobacteria from extreme cold deserts

In the circadian system of cyanobacteria, the ldpA gene is a component of the input to the clock. We comparatively analyzed nucleotide polymorphism of this gene in populations of two closely related species of cyanobacteria (denoted as Synechococcus species S1 and S2, respectively) from extreme cold deserts in Antarctica, the Canadian Arctic, and Tibet. Although both species manifested similarly high haplotype diversities (0.990 and 0.809, respectively), the nucleotide diversity differed significantly (0.0091 in S1 and 0.0037 in S2). The populations of species S2 were more differentiated (F(ST) = 0.2242) compared to those of species S1 (F(ST) between 0.0296 and 0.1188). An analysis of positive selection with several tests yielded highly significant values (P < 0.01) for both species. On the other hand, these results may be somewhat compromised by fluctuating population sizes of the species. The apparent selection pressure coupled with the pronounced demographic factors, such as population expansion, small effective population size, and genetic drift, may thus result in the observed significant interpopulation differentiation and subsequent speciation of cyanobacteria.

Molecular evolution of key genes for type II secretion in Legionella pneumophila

Given the role of type II protein secretion system (T2S) in the ecology and pathogenesis of Legionella pneumophila, it is possible that this system is a target for adaptive evolution. The population genetic structure of L.pneumophila was inferred from the partial sequences of rpoB and from the complete sequence of three T2S structural components (lspD, lspE and pilD) and from two T2S effectors critical for intracellular infection of protozoa (proA and srnA) of 37 strains isolated from natural and man-made environments and disease-related from worldwide sources. A phylogenetic analysis was obtained for the concatenated alignment and for each individual locus. Seven main groups were identified containing the same L.pneumophila strains, suggesting an ancient divergence for each cluster and indicating that the operating selective pressures have equally affected the evolution of the five genes. Although linkage disequilibrium analysis indicate a clonal nature for population structure in this sample, our results indicate that recombination is a common phenomenon among T2S-related genes on this species, as 24 of the 37 L.pneumophila isolates contained at least one locus in which recombination was identified. Furthermore, neutral selection acting on the analysed T2S-related genes emerged as a clear result, namely on T2S effectors, ProA and SrnA, indicating that they are probably implicated in conserved virulence mechanisms through legionellae hosts.

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.

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.

Insertion sequences as highly resolutive genomic markers for sequence type 1 Legionella pneumophila Paris

The causative agent of legionellosis, Legionella pneumophila, colonizes all natural and human-made water networks, thus constituting the source of contaminated aerosols responsible for airborne human infections. Efficient control of infections, especially during epidemics, necessitates the fastest and most resolutive identification possible of the bacterial source for subsequent disinfection of reservoirs. We thus compared recognized typing approaches for Legionella with a method based on characterization of insertion sequence (IS) content. A total of 86 clinical or environmental isolates of L. pneumophila, including 84 Paris isolates, sampled from 25 clinical investigations in France between 2001 and 2007, were obtained from the Legionella National Reference Center. All strains were typed by monoclonal antibody subgrouping, sequence-based typing, pulsed-field gel electrophoresis, and restriction fragment length polymorphism based on the presence or absence of IS elements. We identified six different types of IS elements in L. pneumophila Paris and used them as genomic markers in hybridization experiments. One IS type, ISLpn11, revealed a high discriminatory power. Simpson's index of discrimination, calculated from the distribution of IS elements, was higher than that obtained with the other typing methods used for L. pneumophila Paris. Moreover, specific ISLpn11 copies were found only in strains isolated from particular cities. In more than half of the cases, each clinical isolate had an ISLpn11 profile that was recovered in at least one environmental isolate from the same geographical location, suggesting that our method could identify the infection source. Phylogenetic analysis suggests a clonal expansion for the L. pneumophila Paris strain.

Intragenic tandem repeat variation between Legionella pneumophila strains

Background

Bacterial genomes harbour a large number of tandem repeats, yet the possible phenotypic effects of those found within the coding region of genes are only beginning to be examined. Evidence exists from other organisms that these repeats can be involved in the evolution of new genes, gene regulation, adaptation, resistance to environmental stresses, and avoidance of the immune system.

Results

In this study, we have investigated the presence and variability in copy number of intragenic tandemly repeated sequences in the genome of Legionella pneumophila, the etiological agent of a severe pneumonia known as Legionnaires' disease. Within the genome of the Philadelphia strain, we have identified 26 intragenic tandem repeat sequences using conservative selection criteria. Of these, seven were "polymorphic" in terms of repeat copy number between a large number of L. pneumophila serogroup 1 strains. These strains were collected from a wide variety of environments and patients in several geographical regions. Within this panel of strains, all but one of these seven genes exhibited statistically different patterns in repeat copy number between samples from different origins (environmental, clinical, and hot springs).

Conclusion

These results support the hypothesis that intragenic tandem repeats could play a role in virulence and adaptation to different environments. While tandem repeats are an increasingly popular focus of molecular typing studies in prokaryotes, including in L. pneumophila, this study is the first examining the difference in tandem repeat distribution as a function of clinical or environmental origin.

Population structure and recombination in environmental isolates of Legionella pneumophila

Legionella pneumophila is a water-borne bacteria responsible for most cases of legionellosis, an emerging disease with an increasing incidence in industrialized countries. Although early analysis based on multilocus enzyme electrophoresis (MLEE) described the population structure of this species as clonal, more recent reports have suggested that recombination also contributes to shaping variation across its genome. We report here the results of analysing the nucleotide sequences of 19 loci in 31 environmental samples of L. pneumophila from a small Spanish region (near Alcoi, province of Alicante) where legionellosis has become almost endemic. We analysed the six loci currently incorporated to the sequence-based typing scheme developed by European Working Group for Legionella Infections (EWGLI) for L. pneumophila and 13 intergenic regions, for which we developed primers anchored in flanking, conserved genes. Our results show that recombination among natural isolates of this species is a common phenomenon, as 20 of the 31 isolates contained at least one locus in which recombination was revealed by at least three different methods. The mapping of the recombination events on the maximum likelihood tree of the concatenate sequence of the 19 loci indicated that at least nine independent recombination events might explain the observed distribution of recombinant loci among isolates. In consequence, we have shown that recombination in L. pneumophila is much more frequent than previously considered and that it does not seem to be restricted to already described pathogenicity islands or other genome constituents which provide it with a high plasticity.