Quantifying Nonvertical Inheritance in the Evolution of Legionella pneumophila

Antimicrobial Resistance Genes in Legionella from Artificial Water Systems: Findings from a Two-Year Study

BACKGROUND

Legionella species are the causative agent of Legionnaires' disease and, as ubiquitous waterborne bacteria, are prone to antimicrobial resistance gene (ARG) acquisition and dissemination due to the antimicrobial contamination of natural environments. Given the potential health risks associated with ARGs, it is crucial to assess their presence in the Legionella population.

METHODS

The ARGs lpeAB and tet56 were detected in 348 samples, isolates, and DNA extracts using conventional PCR. In a subset of lpeAB-positive isolates, azithromycin (AZT) MIC values were obtained using the EUCAST protocol and LpeAB activity was evaluated through an efflux pump inhibition assay.

RESULTS

The lpeAB gene was found in 19% (66/348) of samples, with higher detection rates in the L. pneumophila and L. pneumophila sg1 subgroups, at 30% and 41%, respectively. A positive association between lpeAB and L. pneumophila sg1 was found. The MIC values of the lpeAB-positive isolates ranged from 0.064 to 2 mg/L. LpeAB inhibition resulted in 2- and 4-fold MIC reductions in 10 of the 13 isolates analyzed. One sample each of L. longbeacheae and L. bozemanae was found to possess the tet56 gene.

CONCLUSIONS

The lpeAB gene is predominant in L. pneumophila sg1. A few isolates with the lpeAB gene exhibited MIC values below the EUCAST tentative highest MIC values for wild-type isolates. Expanding ARG monitoring in Legionella is essential to assess the public health risk of Legionnaires' disease.

Legionnaires' Disease in China Caused by Legionella pneumophila Corby

Legionella pneumophila is an intracellular pathogen causing pneumonia in humans. In February 2022, Legionnaires' disease caused by L. pneumophila strain Corby in a patient with lung adenocarcinoma was identified for the first time in China. This paper includes the case report and phenotypic and genomic analysis of the Corby (ICDC) strain. Its biological characteristics were evaluated by antibiotic sensitivity testing and cytology experiments, and genomic analysis was performed to understand its genetic evolution. The patient's clinical manifestations included cough, fever, pulmonary infiltration, and significantly decreased activity endurance. After empirical antimicrobial therapy, infection indicators decreased. The Corby (ICDC) strain was susceptible to nine antibiotics and exhibited strong intracellular proliferation ability. A phylogenetic tree showed that the Corby (ICDC) strain was closely related to the Corby strain, but under the pressure of a complex environment, its genome had undergone more rearrangement and inversion. The type IF CRISPR-Cas system was identified in its genome, and spacer analysis indicated that it had been invaded by several foreign plasmids, bacteria, and viruses during evolution. Legionnaires' disease caused by L. pneumophila strain Corby may be ignored in China, and it is urgent to improve long-term monitoring and investigation of aquatic environments and patients with respiratory infections to prevent a large-scale outbreak of Legionnaires' disease.

Reconstructing horizontal gene flow network to understand prokaryotic evolution

Horizontal gene transfer (HGT) is a major source of phenotypic innovation and a mechanism of niche adaptation in prokaryotes. Quantification of HGT is critical to decipher its myriad roles in microbial evolution and adaptation. Advances in genome sequencing and bioinformatics have augmented our ability to understand the microbial world, particularly the direct or indirect influence of HGT on diverse life forms. Methods for detecting HGT can be classified into phylogenetic-based and parametric or composition-based approaches. Here, we exploited the complementary strengths of both the approaches to construct a high confidence horizontal gene flow network. Our network is unique in its ability to detect the transfer of native genes of a genome to genomes from other taxa, thus establishing donor and recipient organisms (taxa), rather than through a post hoc analysis as is the practice with several other approaches. The scale-free horizontal gene flow network presented here provides new insights into modes of transfer for the exchange of genetic information and also illuminates differential gene flow across phyla.

Reconstructing horizontal gene flow network to understand prokaryotic evolution

Horizontal gene transfer (HGT) is a major source of phenotypic innovation and a mechanism of niche adaptation in prokaryotes. Quantification of HGT is critical to decipher its myriad roles in microbial evolution and adaptation. Advances in genome sequencing and bioinformatics have augmented our ability to understand the microbial world, particularly the direct or indirect influence of HGT on diverse life forms. Methods for detecting HGT can be classified into phylogenetic-based and parametric or composition-based approaches. Here, we exploited the complementary strengths of both the approaches to construct a high confidence horizontal gene flow network. Our network is unique in its ability to detect the transfer of native genes of a genome to genomes from other taxa, thus establishing donor and recipient organisms (taxa), rather than through a post hoc analysis as is the practice with several other approaches. The scale-free horizontal gene flow network presented here provides new insights into modes of transfer for the exchange of genetic information and also illuminates differential gene flow across phyla.

The impact of genetic recombination on pathogenic Leptospira

Leptospirosis is the most common zoonosis worldwide, and is increasingly common in poor urban communities, where there is inadequate sewage disposal and abundance of domestic and peridomestic animals. There are many risk factors associated with the disease, such as contaminated water exposure, close contact with animals, floods, recreational activities related to water, wet agriculture. The symptoms of leptospirosis are common to other infectious diseases and, if not treated, it can lead to meningitis, liver failure, kidney damage and death. Leptospirosis is caused by 38 pathogenic species of Leptospira, which are divided in almost 30 serogroups and more than 300 serovars. The serological classification (serogroups and serovars) is based on the expression of distinct lipopolysaccharide (LPS) antigens. These antigens are also associated to protective immunity; antibodies against a serovar protect from any member of the same serovar. Serologic and phylogenetic analyses are not congruent probably due to genetic recombination of LPS genes among different leptospiral species. To analyze the importance of recombination in leptospiral evolution, we performed a gene-by-gene tree topology comparison on closed genomes available in public databases at two levels: among core genomes of pathogenic species (34 recombinant among 1213 core genes), and among core genomes of L. interrogans isolates (178/798). We found that most recombinant genes code for proteins involved in translation, ribosomal structure and biogenesis, but also for cell wall, membrane and envelope biogenesis. Besides, our final results showed that half of LPS genes are recombinant (18/36). This is relevant because serovar classification and vaccine development are based on these epitopes. The frequent recombination of LPS-associated genes suggests that natural selection is promoting the survival of recombinant lineages. These results may help understanding the factors that make Leptospira a successful pathogen.

One is not enough: On the effects of reference genome for the mapping and subsequent analyses of short-reads

Mapping of high-throughput sequencing (HTS) reads to a single arbitrary reference genome is a frequently used approach in microbial genomics. However, the choice of a reference may represent a source of errors that may affect subsequent analyses such as the detection of single nucleotide polymorphisms (SNPs) and phylogenetic inference. In this work, we evaluated the effect of reference choice on short-read sequence data from five clinically and epidemiologically relevant bacteria (Klebsiella pneumoniae, Legionella pneumophila, Neisseria gonorrhoeae, Pseudomonas aeruginosa and Serratia marcescens). Publicly available whole-genome assemblies encompassing the genomic diversity of these species were selected as reference sequences, and read alignment statistics, SNP calling, recombination rates, dN/dS ratios, and phylogenetic trees were evaluated depending on the mapping reference. The choice of different reference genomes proved to have an impact on almost all the parameters considered in the five species. In addition, these biases had potential epidemiological implications such as including/excluding isolates of particular clades and the estimation of genetic distances. These findings suggest that the single reference approach might introduce systematic errors during mapping that affect subsequent analyses, particularly for data sets with isolates from genetically diverse backgrounds. In any case, exploring the effects of different references on the final conclusions is highly recommended.

Mapping consists in the alignment of reads (i.e., DNA fragments) obtained through high-throughput genome sequencing to a previously assembled reference sequence. It is a common practice in genomic studies to use a single reference for mapping, usually the ‘reference genome’ of a species—a high-quality assembly. However, the selection of an optimal reference is hindered by intrinsic intra-species genetic variability, particularly in bacteria. It is known that genetic differences between the reference genome and the read sequences may produce incorrect alignments during mapping. Eventually, these errors could lead to misidentification of variants and biased reconstruction of phylogenetic trees (which reflect ancestry between different bacterial lineages). To our knowledge, this is the first work to systematically examine the effect of different references for mapping on the inference of tree topology as well as the impact on recombination and natural selection inferences. Furthermore, the novelty of this work relies on a procedure that guarantees that we are evaluating only the effect of the reference. This effect has proved to be pervasive in the five bacterial species that we have studied and, in some cases, alterations in phylogenetic trees could lead to incorrect epidemiological inferences. Hence, the use of different reference genomes may be prescriptive to assess the potential biases of mapping.

The Importance of Revisiting Legionellales Diversity

Bacteria of the order Legionellales, such as Legionella pneumophila, the agent of Legionnaires' disease, and Coxiella burnetii, the agent of Q fever, are widely recognized as human pathogens. While our view of the Legionellales is often limited to clinical isolates, ecological surveys are continually uncovering new members of the Legionellales that do not fall into the recognized pathogenic species. Here we emphasize that most of these Legionellales are nonpathogenic forms that have evolved symbiotic lifestyles with nonvertebrate hosts. The diversity of nonpathogenic forms remains, however, largely underexplored. We conjecture that its characterization, once contrasted with the data on pathogenic species, will reveal novel highlights on the mechanisms underlying lifestyle transitions of intracellular bacteria, including the emergence of pathogenesis and mutualism, transmission routes, and host specificity.

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.

Population Genomics of Legionella longbeachae and Hidden Complexities of Infection Source Attribution

Legionella longbeachae is the primary cause of legionellosis in Australasia and Southeast Asia and an emerging pathogen in Europe and the United States; however, our understanding of the population diversity of L. longbeachae from patient and environmental sources is limited. We analyzed the genomes of 64 L. longbeachae isolates, of which 29 were from a cluster of legionellosis cases linked to commercial growing media in Scotland in 2013 and 35 were non-outbreak-associated isolates from Scotland and other countries. We identified extensive genetic diversity across the L. longbeachae species, associated with intraspecies and interspecies gene flow, and a wide geographic distribution of closely related genotypes. Of note, we observed a highly diverse pool of L. longbeachae genotypes within compost samples that precluded the genetic establishment of an infection source. These data represent a view of the genomic diversity of L. longbeachae that will inform strategies for investigating future outbreaks.

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.

Genomic Analysis Reveals Novel Diversity among the 1976 Philadelphia Legionnaires' Disease Outbreak Isolates and Additional ST36 Strains

Legionella pneumophila was first recognized as a cause of severe and potentially fatal pneumonia during a large-scale outbreak of Legionnaires’ disease (LD) at a Pennsylvania veterans’ convention in Philadelphia, 1976. The ensuing investigation and recovery of four clinical isolates launched the fields of Legionella epidemiology and scientific research. Only one of the original isolates, “Philadelphia-1”, has been widely distributed or extensively studied. Here we describe the whole-genome sequencing (WGS), complete assembly, and comparative analysis of all Philadelphia LD strains recovered from that investigation, along with L. pneumophila isolates sharing the Philadelphia sequence type (ST36). Analyses revealed that the 1976 outbreak was due to multiple serogroup 1 strains within the same genetic lineage, differentiated by an actively mobilized, self-replicating episome that is shared with L. pneumophila str. Paris, and two large, horizontally-transferred genomic loci, among other polymorphisms. We also found a completely unassociated ST36 strain that displayed remarkable genetic similarity to the historical Philadelphia isolates. This similar strain implies the presence of a potential clonal population, and suggests important implications may exist for considering epidemiological context when interpreting phylogenetic relationships among outbreak-associated isolates. Additional extensive archival research identified the Philadelphia isolate associated with a non-Legionnaire case of “Broad Street pneumonia”, and provided new historical and genetic insights into the 1976 epidemic. This retrospective analysis has underscored the utility of fully-assembled WGS data for Legionella outbreak investigations, highlighting the increased resolution that comes from long-read sequencing and a sequence type-matched genomic data set.

Dynamics of genome change among Legionella species

Legionella species inhabit freshwater and soil ecosystems where they parasitize protozoa. L. pneumonphila (LP) serogroup-1 (Lp1) is the major cause of Legionnaires' Disease (LD), a life-threatening pulmonary infection that can spread systemically. The increased global frequency of LD caused by Lp and non-Lp species underscores the need to expand our knowledge of evolutionary forces underlying disease pathogenesis. Whole genome analyses of 43 strains, including all known Lp serogroups 1-17 and 17 emergent LD-causing Legionella species (of which 33 were sequenced in this study) in addition to 10 publicly available genomes, resolved the strains into four phylogenetic clades along host virulence demarcations. Clade-specific genes were distinct for genetic exchange and signal-transduction, indicating adaptation to specific cellular and/or environmental niches. CRISPR spacer comparisons hinted at larger pools of accessory DNA sequences in Lp than predicted by the pan-genome analyses. While recombination within Lp was frequent and has been reported previously, population structure analysis identified surprisingly few DNA admixture events between species. In summary, diverse Legionella LD-causing species share a conserved core-genome, are genetically isolated from each other, and selectively acquire genes with potential for enhanced virulence.

Legionella pneumophila strain associated with the first evidence of person-to-person transmission of Legionnaires' disease: a unique mosaic genetic backbone

A first strong evidence of person-to-person transmission of Legionnaires’ Disease (LD) was recently reported. Here, we characterize the genetic backbone of this case-related Legionella pneumophila strain (“PtVFX/2014”), which also caused a large outbreak of LD. PtVFX/2014 is phylogenetically divergent from the most worldwide studied outbreak-associated L. pneumophila subspecies pneumophila serogroup 1 strains. In fact, this strain is also from serogroup 1, but belongs to the L. pneumophila subspecies fraseri. Its genomic mosaic backbone reveals eight horizontally transferred regions encompassing genes, for instance, involved in lipopolysaccharide biosynthesis or encoding virulence-associated Dot/Icm type IVB secretion system (T4BSS) substrates. PtVFX/2014 also inherited a rare ~65 kb pathogenicity island carrying virulence factors and detoxifying enzymes believed to contribute to the emergence of best-fitted strains in water reservoirs and in human macrophages, as well as a inter-species transferred (from L. oakridgensis) ~37.5 kb genomic island (harboring a lvh/lvr T4ASS cluster) that had never been found intact within L. pneumophila species. PtVFX/2014 encodes another lvh/lvr cluster near to CRISPR-associated genes, which may boost L. pneumophila transition from an environmental bacterium to a human pathogen. Overall, this unique genomic make-up may impact PtVFX/2014 ability to adapt to diverse environments, and, ultimately, to be transmitted and cause human disease.

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.

Population structure and minimum core genome typing of Legionella pneumophila

Legionella pneumophila is an important human pathogen causing Legionnaires’ disease. In this study, whole genome sequencing (WGS) was used to study the characteristics and population structure of L. pneumophila strains. We sequenced and compared 53 isolates of L. pneumophila covering different serogroups and sequence-based typing (SBT) types (STs). We found that 1,896 single-copy orthologous genes were shared by all isolates and were defined as the minimum core genome (MCG) of L. pneumophila. A total of 323,224 single-nucleotide polymorphisms (SNPs) were identified among the 53 strains. After excluding 314,059 SNPs which were likely to be results of recombination, the remaining 9,165 SNPs were referred to as MCG SNPs. Population Structure analysis based on MCG divided the 53 L. pneumophila into nine MCG groups. The within-group distances were much smaller than the between-group distances, indicating considerable divergence between MCG groups. MCG groups were also supplied by phylogenetic analysis and may be considered as robust taxonomic units within L. pneumophila. Among the nine MCG groups, eight showed high intracellular growth ability while one showed low intracellular growth ability. Furthermore, MCG typing also showed high resolution in subtyping ST1 strains. The results obtained in this study provided significant insights into the evolution, population structure and pathogenicity of L. pneumophila.

Legionella effectors reflect strength in diversity

The Legionella genus includes opportunistic human pathogenic species that invade human cells using effector proteins that evolved during association with their natural amoeba hosts. A new study compares the genomes of 41 Legionella species to identify nearly 6,000 effectors, providing insight into these species' evolution and pathogenic lifestyles.

Comparative analyses of Legionella species identifies genetic features of strains causing Legionnaires' disease

BACKGROUND

The genus Legionella comprises over 60 species. However, L. pneumophila and L. longbeachae alone cause over 95% of Legionnaires’ disease. To identify the genetic bases underlying the different capacities to cause disease we sequenced and compared the genomes of L. micdadei, L. hackeliae and L. fallonii (LLAP10), which are all rarely isolated from humans.

RESULTS

We show that these Legionella species possess different virulence capacities in amoeba and macrophages, correlating with their occurrence in humans. Our comparative analysis of 11 Legionella genomes belonging to five species reveals highly heterogeneous genome content with over 60% representing species-specific genes; these comprise a complete prophage in L. micdadei, the first ever identified in a Legionella genome. Mobile elements are abundant in Legionella genomes; many encode type IV secretion systems for conjugative transfer, pointing to their importance for adaptation of the genus. The Dot/Icm secretion system is conserved, although the core set of substrates is small, as only 24 out of over 300 described Dot/Icm effector genes are present in all Legionella species. We also identified new eukaryotic motifs including thaumatin, synaptobrevin or clathrin/coatomer adaptine like domains.

CONCLUSIONS

Legionella genomes are highly dynamic due to a large mobilome mainly comprising type IV secretion systems, while a minority of core substrates is shared among the diverse species. Eukaryotic like proteins and motifs remain a hallmark of the genus Legionella. Key factors such as proteins involved in oxygen binding, iron storage, host membrane transport and certain Dot/Icm substrates are specific features of disease-related strains.

Comparative analyses of Legionella species identifies genetic features of strains causing Legionnaires’ disease

Background

The genus Legionella comprises over 60 species. However, L. pneumophila and L. longbeachae alone cause over 95% of Legionnaires’ disease. To identify the genetic bases underlying the different capacities to cause disease we sequenced and compared the genomes of L. micdadei, L. hackeliae and L. fallonii (LLAP10), which are all rarely isolated from humans.

Results

We show that these Legionella species possess different virulence capacities in amoeba and macrophages, correlating with their occurrence in humans. Our comparative analysis of 11 Legionella genomes belonging to five species reveals highly heterogeneous genome content with over 60% representing species-specific genes; these comprise a complete prophage in L. micdadei, the first ever identified in a Legionella genome. Mobile elements are abundant in Legionella genomes; many encode type IV secretion systems for conjugative transfer, pointing to their importance for adaptation of the genus. The Dot/Icm secretion system is conserved, although the core set of substrates is small, as only 24 out of over 300 described Dot/Icm effector genes are present in all Legionella species. We also identified new eukaryotic motifs including thaumatin, synaptobrevin or clathrin/coatomer adaptine like domains.

Conclusions

Legionella genomes are highly dynamic due to a large mobilome mainly comprising type IV secretion systems, while a minority of core substrates is shared among the diverse species. Eukaryotic like proteins and motifs remain a hallmark of the genus Legionella. Key factors such as proteins involved in oxygen binding, iron storage, host membrane transport and certain Dot/Icm substrates are specific features of disease-related strains.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0505-0) contains supplementary material, which is available to authorized users.

Recombination drives genome evolution in outbreak-related Legionella pneumophila isolates

Legionella pneumophila is a strictly environmental pathogen and the etiological agent of legionellosis. It is known that non-vertical processes have a major role in the short-term evolution of pathogens, but little is known about the relevance of these and other processes in environmental bacteria. We report the whole-genome sequencing of 69 L. pneumophila strains linked to recurrent outbreaks in a single location (Alcoy, Spain) over 11 years. We found some examples where the genome sequences of isolates of the same sequence type and outbreak did not cluster together and were more closely related to sequences from different outbreaks. Our analyses identify 16 recombination events responsible for almost 98% of the SNPs detected in the core genome and an apparent acceleration in the evolutionary rate. These results have profound implications for the understanding of microbial populations and for public health interventions in Legionella outbreak investigations.

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.

Classification of Streptomyces phylogroup pratensis (Doroghazi and Buckley, 2010) based on genetic and phenotypic evidence, and proposal of Streptomyces pratensis sp. nov

The Streptomyces phylogroup pratensis (Doroghazi and Buckley, 2010) contains isolates obtained from grassy fields, as well as Streptomyces flavogriseus ATCC 33331 and strain CGMCC 4.1868. This latter strain was received as Streptomyces griseoplanus but was subsequently found to be mislabeled, and S. flavogriseus ATCC 33331 (=IAF-45-CD) was shown to be clearly distinct from the type strain S. flavogriseus ATCC 25452(T) (=CGMCC 4.1884(T)). In order to evaluate the taxonomic position of phylogroup pratensis further, sequences of the 16S rRNA gene and five protein-coding housekeeping genes (atpD, gyrB, recA, rpoB and trpB) were determined for six strains of the phylogroup and type strains of 19 related species, which were selected by a BLAST search based on the sequences of the phylogroup. The 16S rRNA gene sequences for the phylogroup were identical to those of eight species belonging to cluster I of the S. griseus clade. However, in all the individual protein-coding gene and MLSA phylogenies, the phylogroup strains without exception formed an obviously distinct cluster that could be equated with a new species status. The phylogenetic evidence for the new species assignment was also supported by corresponding DNA-DNA hybridization values and by phenotypic characteristics. It is therefore proposed that the phylogroup should be classified as Streptomyces pratensis sp. nov., and the type strain is ch24(T) (=CGMCC 4.6829(T)=NRRL B-24916(T)).

Quantifying homologous replacement of loci between haloarchaeal species

In vitro studies of the haloarchaeal genus Haloferax have demonstrated their ability to frequently exchange DNA between species, whereas rates of homologous recombination estimated from natural populations in the genus Halorubrum are high enough to maintain random association of alleles between five loci. To quantify the effects of gene transfer and recombination of commonly held (relaxed core) genes during the evolution of the class Halobacteria (haloarchaea), we reconstructed the history of 21 genomes representing all major groups. Using a novel algorithm and a concatenated ribosomal protein phylogeny as a reference, we created a directed horizontal genetic transfer (HGT) network of contemporary and ancestral genomes. Gene order analysis revealed that 90% of testable HGTs were by direct homologous replacement, rather than nonhomologous integration followed by a loss. Network analysis revealed an inverse log-linear relationship between HGT frequency and ribosomal protein evolutionary distance that is maintained across the deepest divergences in Halobacteria. We use this mathematical relationship to estimate the total transfers and amino acid substitutions delivered by HGTs in each genome, providing a measure of chimerism. For the relaxed core genes of each genome, we conservatively estimate that 11–20% of their evolution occurred in other haloarchaea. Our findings are unexpected, because the transfer and homologous recombination of relaxed core genes between members of the class Halobacteria disrupts the coevolution of genes; however, the generation of new combinations of divergent but functionally related genes may lead to adaptive phenotypes not available through cumulative mutations and recombination within a single population.

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.

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.

Rhizome of life, catastrophes, sequence exchanges, gene creations, and giant viruses: how microbial genomics challenges Darwin

Darwin's theory about the evolution of species has been the object of considerable dispute. In this review, we have described seven key principles in Darwin's book The Origin of Species and tried to present how genomics challenge each of these concepts and improve our knowledge about evolution. Darwin believed that species evolution consists on a positive directional selection ensuring the “survival of the fittest.” The most developed state of the species is characterized by increasing complexity. Darwin proposed the theory of “descent with modification” according to which all species evolve from a single common ancestor through a gradual process of small modification of their vertical inheritance. Finally, the process of evolution can be depicted in the form of a tree. However, microbial genomics showed that evolution is better described as the “biological changes over time.” The mode of change is not unidirectional and does not necessarily favors advantageous mutations to increase fitness it is rather subject to random selection as a result of catastrophic stochastic processes. Complexity is not necessarily the completion of development: several complex organisms have gone extinct and many microbes including bacteria with intracellular lifestyle have streamlined highly effective genomes. Genomes evolve through large events of gene deletions, duplications, insertions, and genomes rearrangements rather than a gradual adaptative process. Genomes are dynamic and chimeric entities with gene repertoires that result from vertical and horizontal acquisitions as well as de novo gene creation. The chimeric character of microbial genomes excludes the possibility of finding a single common ancestor for all the genes recorded currently. Genomes are collections of genes with different evolutionary histories that cannot be represented by a single tree of life (TOL). A forest, a network or a rhizome of life may be more accurate to represent evolutionary relationships among species.

Lateral gene exchanges shape the genomes of amoeba-resisting microorganisms

Based on Darwin's concept of the tree of life, vertical inheritance was thought to be dominant, and mutations, deletions, and duplication were streaming the genomes of living organisms. In the current genomic era, increasing data indicated that both vertical and lateral gene inheritance interact in space and time to trigger genome evolution, particularly among microorganisms sharing a given ecological niche. As a paradigm to their diversity and their survival in a variety of cell types, intracellular microorganisms, and notably intracellular bacteria, were considered as less prone to lateral genetic exchanges. Such specialized microorganisms generally have a smaller gene repertoire because they do rely on their host's factors for some basic regulatory and metabolic functions. Here we review events of lateral gene transfer (LGT) that illustrate the genetic exchanges among intra-amoebal microorganisms or between the microorganism and its amoebal host. We tentatively investigate the functions of laterally transferred genes in the light of the interaction with their host as they should confer a selective advantage and success to the amoeba-resisting microorganisms (ARMs).

How microbiology helps define the rhizome of life

In contrast to the tree of life (TOF) theory, species are mosaics of gene sequences with different origins. Observations of the extensive lateral sequence transfers in all organisms have demonstrated that the genomes of all life forms are collections of genes with different evolutionary histories that cannot be represented by a single TOF. Moreover, genes themselves commonly have several origins due to recombination. The human genome is not free from recombination events, so it is a mosaic like other organisms' genomes. Recent studies have demonstrated evidence for the integration of parasitic DNA into the human genome. Lateral transfer events have been accepted as major contributors of genome evolution in free-living bacteria. Furthermore, the accumulation of genomic sequence data provides evidence for extended genetic exchanges in intracellular bacteria and suggests that such events constitute an agent that promotes and maintains all bacterial species. Archaea and viruses also form chimeras containing primarily bacterial but also eukaryotic sequences. In addition to lateral transfers, orphan genes are indicative of the fact that gene creation is a permanent and unsettled phenomenon. Currently, a rhizome may more adequately represent the multiplicity and de novo creation of a genome. We wanted to confirm that the term “rhizome” in evolutionary biology applies to the entire cellular life history. This view of evolution should resemble a clump of roots representing the multiple origins of the repertoires of the genes of each species.

Phylogenomics of T4 cyanophages: lateral gene transfer in the 'core' and origins of host genes

The last two decades have revealed that phages (viruses that infect bacteria) are abundant and play fundamental roles in the Earth System, with the T4-like myoviruses (herein T4-like phages) emerging as a dominant 'signal' in wild populations. Here we examine 27 T4-like phage genomes, with a focus on 17 that infect ocean picocyanobacteria (cyanophages), to evaluate lateral gene transfer (LGT) in this group. First, we establish a reference tree by evaluating concatenated core gene supertrees and whole genome gene content trees. Next, we evaluate what fraction of these 'core genes' shared by all 17 cyanophages appear prone to LGT. Most (47 out of 57 core genes) were vertically transferred as inferred from tree tests and genomic synteny. Of those 10 core genes that failed the tree tests, the bulk (8 of 10) remain syntenic in the genomes with only a few (3 of the 10) having identifiable signatures of mobile elements. Notably, only one of these 10 is shared not only by the 17 cyanophages, but also by all 27 T4-like phages (thymidylate synthase); its evolutionary history suggests cyanophages may be the origin of these genes to Prochlorococcus. Next, we examined intragenic recombination among the core genes and found that it did occur, even among these core genes, but that the rate was significantly higher between closely related phages, perhaps reducing any detectable LGT signal and leading to taxon cohesion. Finally, among 18 auxiliary metabolic genes (AMGs, a.k.a. 'host' genes), we found that half originated from their immediate hosts, in some cases multiple times (e.g. psbA, psbD, pstS), while the remaining have less clear evolutionary origins ranging from cyanobacteria (4 genes) or microbes (5 genes), with particular diversity among viral TalC and Hsp20 sequences. Together, these findings highlight the patterns and limits of vertical evolution, as well as the ecological and evolutionary roles of LGT in shaping T4-like phage genomes.

The rhizome of life: the sympatric Rickettsia felis paradigm demonstrates the random transfer of DNA sequences

The intracellular flea symbiont, Rickettsia felis, may meet other organisms intracellularly such as R. typhi. We used a single-gene phylogenetic approach of the 1375 R. felis genes to look for horizontal transfers that occurred as a result of the bacterial promiscuity with other organisms. Our results showed that besides genes that are linked to the Spotted Fever Group, 165 genes have a different history and are linked to other Rickettsia such as R. bellii (107 genes), R. typhi (15 genes), or to other bacteria such as Legionella sp. and Francisella sp. or to eukaryotes. Among these genes, we identified 73 individual genes and 34 spatial clusters containing 2-4 adjacent genes, a total of 79 genes, with evidence of en bloc transfer. We described 13 chimeric genes resulting from gene recombination with sympatric R. typhi. The transferred DNA sequences present different sizes and functions, suggesting that the horizontal transfer in R. felis is random and neutral within its specific host. Our study shows that the strict intracellular bacteria R. felis exhibits a mosaic genome. We therefore developed a new representation for the evolutionary history of R. felis showing its different putative ancestors in the form of a rhizome.

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.

The influence of rickettsiologists on post-modern microbiology

Many of the definitions in microbiology are currently false. We have reviewed the great denominations of microbiology and attempted to free microorganisms from the theories of the twentieth century. The presence of compartmentation and a nucleoid in Planctomycetes clearly calls into question the accuracy of the definitions of eukaryotes and prokaryotes. Archaea are viewed as prokaryotes resembling bacteria. However, the name archaea, suggesting an archaic origin of lifestyle, is inconsistent with the lifestyle of this family. Viruses are defined as small, filterable infectious agents, but giant viruses challenge the size criteria used for the definition of a virus. Pathogenicity does not require the acquisition of virulence factors (except for toxins), and in many cases, gene loss is significantly inked to the emergence of virulence. Species classification based on 16S rRNA is useless for taxonomic purposes of human pathogens, as a 2% divergence would classify all Rickettsiae within the same species and would not identify bacteria specialized for mammal infection. The use of metagenomics helps us to understand evolution and physiology by elucidating the structure, function, and interactions of the major microbial communities, but it neglects the minority populations. Finally, Darwin’s descent with modification theory, as represented by the tree of life, no longer matches our current genomic knowledge because genomics has revealed the occurrence of de novo-created genes and the mosaic structure of genomes, the Rhizome of life is therefore more appropriate.

Trends and barriers to lateral gene transfer in prokaryotes

Gene acquisition by lateral gene transfer (LGT) is an important mechanism for natural variation among prokaryotes. Laboratory experiments show that protein-coding genes can be laterally transferred extremely fast among microbial cells, inherited to most of their descendants, and adapt to a new regulatory regime within a short time. Recent advance in the phylogenetic analysis of microbial genomes using networks approach reveals a substantial impact of LGT during microbial genome evolution. Phylogenomic networks of LGT among prokaryotes reconstructed from completely sequenced genomes uncover barriers to LGT in multiple levels. Here we discuss the kinds of barriers to gene acquisition in nature including physical barriers for gene transfer between cells, genomic barriers for the integration of acquired DNA, and functional barriers for the acquisition of new genes.

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.

Origins of bacterial diversity through horizontal genetic transfer and adaptation to new ecological niches

Horizontal genetic transfer (HGT) has played an important role in bacterial evolution at least since the origins of the bacterial divisions, and HGT still facilitates the origins of bacterial diversity, including diversity based on antibiotic resistance. Adaptive HGT is aided by unique features of genetic exchange in bacteria such as the promiscuity of genetic exchange and the shortness of segments transferred. Genetic exchange rates are limited by the genetic and ecological similarity of organisms. Adaptive transfer of genes is limited to those that can be transferred as a functional unit, provide a niche-transcending adaptation, and are compatible with the architecture and physiology of other organisms. Horizontally transferred adaptations may bring about fitness costs, and natural selection may ameliorate these costs. The origins of ecological diversity can be analyzed by comparing the genomes of recently divergent, ecologically distinct populations, which can be discovered as sequence clusters. Such genome comparisons demonstrate the importance of HGT in ecological diversification. Newly divergent populations cannot be discovered as sequence clusters when their ecological differences are coded by plasmids, as is often the case for antibiotic resistance; the discovery of such populations requires a screen for plasmid-coded functions. This paper reviews the features of bacterial genetics that allow HGT, the similarities between organisms that foster HGT between them, the limits to the kinds of adaptations that can be transferred, and amelioration of fitness costs associated with HGT; the paper also reviews approaches to discover the origins of new, ecologically distinct bacterial populations and the role that HGT plays in their founding.

Of woods and webs: possible alternatives to the tree of life for studying genomic fluidity in E. coli

Background

We introduce several forest-based and network-based methods for exploring microbial evolution, and apply them to the study of thousands of genes from 30 strains of E. coli. This case study illustrates how additional analyses could offer fast heuristic alternatives to standard tree of life (TOL) approaches.

Results

We use gene networks to identify genes with atypical modes of evolution, and genome networks to characterize the evolution of genetic partnerships between E. coli and mobile genetic elements. We develop a novel polychromatic quartet method to capture patterns of recombination within E. coli, to update the clanistic toolkit, and to search for the impact of lateral gene transfer and of pathogenicity on gene evolution in two large forests of trees bearing E. coli. We unravel high rates of lateral gene transfer involving E. coli (about 40% of the trees under study), and show that both core genes and shell genes of E. coli are affected by non-tree-like evolutionary processes. We show that pathogenic lifestyle impacted the structure of 30% of the gene trees, and that pathogenic strains are more likely to transfer genes with one another than with non-pathogenic strains. In addition, we propose five groups of genes as candidate mobile modules of pathogenicity. We also present strong evidence for recent lateral gene transfer between E. coli and mobile genetic elements.

Conclusions

Depending on which evolutionary questions biologists want to address (i.e. the identification of modules, genetic partnerships, recombination, lateral gene transfer, or genes with atypical evolutionary modes, etc.), forest-based and network-based methods are preferable to the reconstruction of a single tree, because they provide insights and produce hypotheses about the dynamics of genome evolution, rather than the relative branching order of species and lineages. Such a methodological pluralism - the use of woods and webs - is to be encouraged to analyse the evolutionary processes at play in microbial evolution.

This manuscript was reviewed by: Ford Doolittle, Tal Pupko, Richard Burian, James McInerney, Didier Raoult, and Yan Boucher

Intracellular Rickettsiales: Insights into manipulators of eukaryotic cells

The order Rickettsiales comprises obligate intracellular bacteria that are the ancestors of modern eukaryotes. These bacteria infect various vectors and hosts, with some species being pathogenic to man. Rickettsiales have small, degraded genomes and provide a paradigm for increased pathogenicity despite gene loss; significant levels of genetic exchange occur between bacteria that infect the same host and with the eukaryotic hosts themselves. Crosstalk between host and bacteria appears to be mediated by a Type IV secretion system and proteins containing eukaryotic-like repeat motifs. Rickettsiales also manipulate host reproduction and induce host resistance to viruses. Manipulation of its host by Rickettsiales has long been misunderstood because of technical difficulties, but recent advances in understanding bacterial-eukaryotes interactions have been made and are reviewed here.