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

Recombination in Bacterial Genomes: Evolutionary Trends

Bacterial organisms have undergone homologous recombination (HR) and horizontal gene transfer (HGT) multiple times during their history. These processes could increase fitness to new environments, cause specialization, the emergence of new species, and changes in virulence. Therefore, comprehensive knowledge of the impact and intensity of genetic exchanges and the location of recombination hotspots on the genome is necessary for understanding the dynamics of adaptation to various conditions. To this end, we aimed to characterize the functional impact and genomic context of computationally detected recombination events by analyzing genomic studies of any bacterial species, for which events have been detected in the last 30 years. Genomic loci where the transfer of DNA was detected pertained to mobile genetic elements (MGEs) housing genes that code for proteins engaged in distinct cellular processes, such as secretion systems, toxins, infection effectors, biosynthesis enzymes, etc. We found that all inferences fall into three main lifestyle categories, namely, ecological diversification, pathogenesis, and symbiosis. The latter primarily exhibits ancestral events, thus, possibly indicating that adaptation appears to be governed by similar recombination-dependent mechanisms.

Rapid adaptations of Legionella pneumophila to the human host

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

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

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

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

Background

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

Methods

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

Results

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

Conclusions

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

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.

Genomic characterization and assessment of pathogenic potential of Legionella spp. isolates from environmental monitoring

Several species in the genus Legionella are known to cause an acute pneumonia when the aerosols containing the bacteria from man-made water systems are inhaled. The disease is usually caused by Legionella pneumophila, but other species have been implicated in the infection. The disease is frequently manifested as an outbreak, which means several people are affected when exposed to the common source of Legionella contamination. Therefor environmental surveillance which includes isolation and identification of Legionella is performed routinely. However, usually no molecular or genome-based methods are employed in further characterization of the isolates during routine environmental monitoring. During several years of such monitoring, isolates from different geographical locations were collected and 39 of them were sequenced by hybrid de novo approach utilizing short and long sequencing reads. In addition, the isolates were typed by standard culture and MALDI-TOF method. The sequencing reads were assembled and annotated to produce high-quality genomes. By employing discriminatory genome typing, four potential new species in the Legionella genus were identified, which are yet to be biochemically and morphologically characterized. Moreover, functional annotations concerning virulence and antimicrobial resistance were performed on the sequenced genomes. The study contributes to the knowledge on little-known non-pneumophila species present in man-made water systems and establishes support for future genetic relatedness studies as well as understanding of their pathogenic potential.

Concept about the Virulence Factor of Legionella

Pathogenic species of Legionella can infect human alveolar macrophages through Legionella-containing aerosols to cause a disease called Legionellosis, which has two forms: a flu-like Pontiac fever and severe pneumonia named Legionnaires' disease (LD). Legionella is an opportunistic pathogen that frequently presents in aquatic environments as a biofilm or protozoa parasite. Long-term interaction and extensive co-evolution with various genera of amoebae render Legionellae pathogenic to infect humans and also generate virulence differentiation and heterogeneity. Conventionally, the proteins involved in initiating replication processes and human macrophage infections have been regarded as virulence factors and linked to pathogenicity. However, because some of the virulence factors are associated with the infection of protozoa and macrophages, it would be more accurate to classify them as survival factors rather than virulence factors. Given that the molecular basis of virulence variations among non-pathogenic, pathogenic, and highly pathogenic Legionella has not yet been elaborated from the perspective of virulence factors, a comprehensive explanation of how Legionella infects its natural hosts, protozoans, and accidental hosts, humans is essential to show a novel concept regarding the virulence factor of Legionella. In this review, we overviewed the pathogenic development of Legionella from protozoa, the function of conventional virulence factors in the infections of protozoa and macrophages, the host's innate immune system, and factors involved in regulating the host immune response, before discussing a probably new definition for the virulence factors of Legionella.

Ubiquitin and Legionella: From bench to bedside

Legionella pneumophila, a Gram-negative intracellular bacterium, is one of the major causes of Legionnaires' disease, a specific type of atypical pneumonia. Despite intensive research efforts that elucidated many relevant structural, molecular and medical insights into Legionella's pathogenicity, Legionnaires' disease continues to present an ongoing public health concern. Legionella's virulence is based on its ability to simultaneously hijack multiple molecular pathways of the host cell to ensure its fast replication and dissemination. Legionella usurps the host ubiquitin system through multiple effector proteins, using the advantage of both conventional and unconventional (phosphoribosyl-linked) ubiquitination, thus providing optimal conditions for its replication. In this review, we summarize the current understanding of L. pneumophila from medical, biochemical and molecular perspectives. We describe the clinical disease presentation, its diagnostics and treatment, as well as host-pathogen interactions, with the emphasis on the ability of Legionella to target the host ubiquitin system upon infection. Furthermore, the interdisciplinary use of innovative technologies enables better insights into the pathogenesis of Legionnaires' disease and provides new opportunities for its treatment and prevention.

Legionella pneumophila: The Journey from the Environment to the Blood

An outbreak of a potentially fatal form of pneumonia in 1976 and in the annual convention of the American Legion was the first time that Legionella spp. was identified. Thereafter, the term Legionnaires’ disease (LD) was established. The infection in humans is transmitted by the inhalation of aerosols that contain the microorganisms that belong to the Legionellaceae family and the genus Legionella. The genus Legionella contains genetically heterogeneous species and serogroups. The Legionella pneumophila serogroup 1 (Lp1) is the most often detected strain in outbreaks of LD. The pathogenesis of LD infection initiates with the attachment of the bacterial cells to the host cells, and subsequent intracellular replication. Following invasion, Legionella spp. activates its virulence mechanisms: generation of specific compartments of Legionella-containing vacuole (LCV), and expression of genes that encode a type IV secretion system (T4SS) for the translocation of proteins. The ability of L. pneumophila to transmigrate across the lung’s epithelium barrier leads to bacteremia, spread, and invasion of many organs with subsequent manifestations, complications, and septic shock. The clinical manifestations of LD depend on the bacterial load in the aerosol, the virulence factors, and the immune status of the patient. The infection has two distinct forms: the non- pneumatic form or Pontiac fever, which is a milder febrile flu-like illness, and LD, a more severe form, which includes pneumonia. In addition, the extrapulmonary involvement of LD can include heart, brain, abdomen, and joints.

Presence of Viable, Clinically Relevant Legionella Bacteria in Environmental Water and Soil Sources of China

The distribution of pathogenic Legionella in the environmental soil and water of China has not been documented yet. In this study, Legionella was detected in 129 of 575 water (22.43%) and 41 of 442 soil samples (9.28%) by culture. Twelve Legionella species were identified, of which 11 were disease-associated. Of the Legionella-positive samples, 109 of 129 (84.50%) water and 29 of 41 (70.73%) soil were positive for L. pneumophila, which accounted for about 75% of Legionella isolates in both water and soil, suggesting L. pneumophila was the most frequent species. Soil showed a higher diversity of Legionella spp. as compared with water (0.6279 versus 0.4493). In contrast, serogroup (sg) 1 was more prevalent among L. pneumophila isolates from water than from soil (26.66% versus 12.21%). Moreover, many disease-associated sequence types (STs) of L. pneumophila were found in China. Intragenic recombination was acting on L. pneumophila from both water and soil. Phylogeny, population structure, and molecular evolution analyses revealed a probable existence of L. pneumophila isolates with a special genetic background that is more adaptable to soil or water sources and a small proportion of genetic difference between water and soil isolates. The detection of viable, clinically relevant Legionella demonstrates soil as another source for harboring and dissemination of pathogenic Legionella bacteria in China. Future research should assess the implication in public health with the presence of Legionella in the soil and illustrate the genetic and pathogenicity difference of Legionella between water and soil, particularly the most prevalent L. pneumophila. IMPORTANCE Pathogenic Legionella spp. is the causative agent of Legionnaires' disease (LD), and L. pneumophila is the most common one. Most studies have focused on L. pneumophila from water and clinical samples. However, the soil is another important reservoir for this bacterium, and the distribution of Legionella spp. in water and soil sources has not been compared and documented in China yet. Discovering the distribution of Legionella spp. and L. pneumophila in the two environments may help a deep understanding of the pathogenesis and molecular evolution of the bacterium. Our research systematically uncovered the distributions of Legionella spp. in different regions and sources (e.g., water and soil) of China. Moreover, phylogeny, population structure, and molecular evolution study revealed the possible existence of L. pneumophila with a special genetic background that is more adaptable to soil or water sources, and genetic difference may exist.

Horizontal Gene Transfer of Fluoroquinolone Resistance-Conferring Genes From Commensal Neisseria to Neisseria gonorrhoeae: A Global Phylogenetic Analysis of 20,047 Isolates

Antimicrobial resistance in Neisseria gonorrhoeae is an important global health concern. The genetically related commensal Neisseria act as a reservoir of resistance genes, and horizontal gene transfer (HGT) has been shown to play an important role in the genesis of resistance to cephalosporins and macrolides in N. gonorrhoeae. In this study, we evaluated if there was evidence of HGT in the genes gyrA/gyrB and parC/parE responsible for fluoroquinolone resistance. Even though the role of gyrB and parE in quinolone resistance is unclear, the subunits gyrB and parE were included as zoliflodacin, a promising new drug to treat N. gonorrhoeae targets the gyrB subunit. We analyzed a collection of 20,047 isolates; 18,800 N. gonorrhoeae, 1,238 commensal Neisseria spp., and nine Neisseria meningitidis. Comparative genomic analyses identified HGT events in genes, gyrA, gyrB, parC, and parE. Recombination events were predicted in N. gonorrhoeae and Neisseria commensals. Neisseria lactamica, Neisseria macacae, and Neisseria mucosa were identified as likely progenitors of the HGT events in gyrA, gyrB, and parE, respectively.

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

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

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

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

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

Legionella pneumophila regulates host cell motility by targeting Phldb2 with a 14-3-3ζ-dependent protease effector

The cytoskeleton network of eukaryotic cells is essential for diverse cellular processes, including vesicle trafficking, cell motility, and immunity, thus is a common target for bacterial virulence factors. A number of effectors from the bacterial pathogen Legionella pneumophila have been shown to modulate the function of host actin cytoskeleton to construct the Legionella-containing vacuole (LCV) permissive for its intracellular replication. In this study, we found that the Dot/Icm effector Lem8 (Lpg1290) is a protease whose activity is catalyzed by a Cys-His-Asp motif known to be associated with diverse biochemical activities. Intriguingly, we found that Lem8 interacts with the host regulatory protein 14-3-3ζ, which activates its protease activity. Furthermore, Lem8 undergoes self-cleavage in a process that requires 14-3-3ζ. We identified the Pleckstrin homology-like domain-containing protein Phldb2 involved in cytoskeleton organization as a target of Lem8 and demonstrated that Lem8 plays a role in the inhibition of host cell migration by attacking Phldb2.

Translocated Legionella pneumophila small RNAs mimic eukaryotic microRNAs targeting the host immune response

Legionella pneumophila is an intracellular bacterial pathogen that can cause a severe form of pneumonia in humans, a phenotype evolved through interactions with aquatic protozoa in the environment. Here, we show that L. pneumophila uses extracellular vesicles to translocate bacterial small RNAs (sRNAs) into host cells that act on host defence signalling pathways. The bacterial sRNA RsmY binds to the UTR of ddx58 (RIG-I encoding gene) and cRel, while tRNA-Phe binds ddx58 and irak1 collectively reducing expression of RIG-I, IRAK1 and cRel, with subsequent downregulation of IFN-β. Thus, RsmY and tRNA-Phe are bacterial trans-kingdom regulatory RNAs downregulating selected sensor and regulator proteins of the host cell innate immune response. This miRNA-like regulation of the expression of key sensors and regulators of immunity is a feature of L. pneumophila host-pathogen communication and likely represents a general mechanism employed by bacteria that interact with eukaryotic hosts.

Dynamics of Legionella Community Interactions in Response to Temperature and Disinfection Treatment: 7 Years of Investigation

In man-made water distribution systems, Legionella community interactions remain unknown, due to their ability to change from sessile to planktonic states or live in viable but non-culturable forms, in response to anthropic and environmental stress. During 7 years of hospital Legionella surveillance, in 191 hot water positive samples, the interactions among the Legionella species, temperature, and disinfection treatment were evaluated. Legionella was isolated following ISO 11731:2017, and identification was performed by mip gene sequencing and sequence-based typing (SBT) for L. anisa or L. rubrilucens and L. pneumophila, respectively. The species with the higher frequency of isolation was L. pneumophila serogroup 1 (78.53%; 4865.36 ± 25,479.11 cfu/L), followed by L. anisa (54.45%; 558.79 ± 2637.41 cfu/L) and L. rubrilucens (21.99%; 307.73 ± 1574.95 cfu/L), which were sometimes present together. Spearman's rho correlation test was conducted among the species with respect to temperature and disinfectant (H₂O₂/Ag⁺). The results showed a generally positive interaction among these species sharing the same environment, except for competition between L. anisa and L. rubrilucens. High temperature (48.83 ± 2.59 °C) and disinfection treatment (11.58 ± 4.99 mg/L) affected the presence of these species. An exception was observed with L. anisa, which showed disinfection treatment resistance. For the purposes of environmental surveillance, it is fundamental to better understand the interactions and dynamic of the Legionella community in man-made water systems in order to choose the proper physical or chemical treatments. The simultaneous presence of different Legionella species could result in an increased resistance to high temperature and disinfectant treatment, leading to changes in contamination level and species diversity.

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

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

Non-canonical activation of the ER stress sensor ATF6 by Legionella pneumophila effectors

Legionella pneumophila secretes toxins into the host cell that induce the non-canonical processing and activation of the ER stress sensor and transcription factor ATF6 via a mechanism that is distinct from the canonical pathway activated by unfolded protein buildup.

The intracellular bacterial pathogen Legionella pneumophila (L.p.) secretes ∼330 effector proteins into the host cell to sculpt an ER-derived replicative niche. We previously reported five L.p. effectors that inhibit IRE1, a key sensor of the homeostatic unfolded protein response (UPR) pathway. In this study, we discovered a subset of L.p. toxins that selectively activate the UPR sensor ATF6, resulting in its cleavage, nuclear translocation, and target gene transcription. In a deviation from the conventional model, this L.p.–dependent activation of ATF6 does not require its transport to the Golgi or its cleavage by the S1P/S2P proteases. We believe that our findings highlight the unique regulatory control that L.p. exerts upon the three UPR sensors and expand the repertoire of bacterial proteins that selectively perturb host homeostatic pathways.

Virulence effector SidJ evolution in Legionella pneumophila is driven by positive selection and intragenic recombination

Effector proteins translocated by the Dot/Icm type IV secretion system determine the virulence of Legionella pneumophila (L. pneumophila). Among these effectors, members of the SidE family (SidEs) regulate several cellular processes through a unique phosphoribosyl ubiquitination mechanism mediated by another effector, SidJ. Host-cell calmodulin (CaM) activates SidJ to glutamylate the SidEs of ubiquitin (Ub) ligases and to make a balanced Ub ligase activity. Given the central role of SidJ in this regulatory process, studying the nature of evolution of sidJ is important to understand the virulence of L. pneumophila and the interaction between the bacteria and its hosts. By studying sidJ from a large number of L. pneumophila strains and using various molecular evolution algorithms, we demonstrated that intragenic recombination drove the evolution of sidJ and contributed to sidJ diversification. Additionally, we showed that four codons of sidJ which are located in the N-terminal (NTD) (codons 58 and 200) and C-terminal (CTD) (codons 868 and 869) domains, but not in the kinase domain (KD) had been subjected to strong positive selection pressure, and variable mutation profiles of these codons were identified. Protein structural modeling of SidJ provided possible explanations for these mutations. Codons 868 and 869 mutations might engage in regulating the interactions of SidJ with CaM through hydrogen bonds and affect the CaM docking to SidJ. Mutation in codon 58 of SidJ might affect the distribution of main-chain atoms that are associated with the interaction with CaM. In contrast, mutations in codon 200 might influence the α-helix stability in the NTD. These mutations might be important to balance Ub ligase activity for different L. pneumophila hosts. This study first reported that intragenic recombination and positive Darwinian selection both shaped the genetic plasticity of sidJ, contributing to a deeper understanding of the adaptive mechanisms of this intracellular bacterium to different hosts.

Legionella pneumophila CRISPR-Cas Suggests Recurrent Encounters with One or More Phages in the Family Microviridae

Legionella pneumophila is a ubiquitous freshwater pathogen and the causative agent of Legionnaires' disease. L. pneumophila growth within protists provides a refuge from desiccation, disinfection, and other remediation strategies. One outstanding question has been whether this protection extends to phages. L. pneumophila isolates are remarkably devoid of prophages and to date no Legionella phages have been identified. Nevertheless, many L. pneumophila isolates maintain active CRISPR-Cas defenses. So far, the only known target of these systems is an episomal element that we previously named Legionella mobile element 1 (LME-1). The continued expansion of publicly available genomic data promises to further our understanding of the role of these systems. We now describe over 150 CRISPR-Cas systems across 600 isolates to establish the clearest picture yet of L. pneumophila's adaptive defenses. By searching for targets of 1,500 unique CRISPR-Cas spacers, LME-1 remains the only identified CRISPR-Cas targeted integrative element. We identified 3 additional LME-1 variants-all targeted by previously and newly identified CRISPR-Cas spacers-but no other similar elements. Notably, we also identified several spacers with significant sequence similarity to microviruses, specifically those within the subfamily Gokushovirinae. These spacers are found across several different CRISPR-Cas arrays isolated from geographically diverse isolates, indicating recurrent encounters with these phages. Our analysis of the extended Legionella CRISPR-Cas spacer catalog leads to two main conclusions: current data argue against CRISPR-Cas targeted integrative elements beyond LME-1, and the heretofore unknown L. pneumophila phages are most likely lytic gokushoviruses. IMPORTANCE Legionnaires' disease is an often-fatal pneumonia caused by Legionella pneumophila, which normally grows inside amoebae and other freshwater protists. L. pneumophila trades diminished access to nutrients for the protection and isolation provided by the host. One outstanding question is whether L. pneumophila is susceptible to phages, given the protection provided by its intracellular lifestyle. In this work, we use Legionella CRISPR spacer sequences as a record of phage infection to predict that the "missing" L. pneumophila phages belong to the microvirus subfamily Gokushovirinae. Gokushoviruses are known to infect another intracellular pathogen, Chlamydia. How do gokushoviruses access L. pneumophila (and Chlamydia) inside their "cozy niches"? Does exposure to phages happen during a transient extracellular period (during cell-to-cell spread) or is it indicative of a more complicated environmental lifestyle? One thing is clear, 100 years after their discovery, phages continue to hold important secrets about the bacteria upon which they prey.

Evidence of Horizontal Gene Transfer of 50S Ribosomal Genes rplB, rplD, and rplY in Neisseria gonorrhoeae

Horizontal gene transfer (HGT) in the penA and multidrug efflux pump genes has been shown to play a key role in the genesis of antimicrobial resistance in Neisseria gonorrhoeae. In this study, we evaluated if there was evidence of HGT in the genes coding for the ribosomal proteins in the Neisseria genus. We did this in a collection of 11,659 isolates of Neisseria, including N. gonorrhoeae and commensal Neisseria species (N. cinerea, N. elongata, N. flavescens, N. mucosa, N. polysaccharea, and N. subflava). Comparative genomic analyses identified HGT events in three genes: rplB, rplD, and rplY coding for ribosomal proteins L2, L4 and L25, respectively. Recombination events were predicted in N. gonorrhoeae and N. cinerea, N. subflava, and N. lactamica were identified as likely progenitors. In total, 2,337, 2,355, and 1,127 isolates possessed L2, L4, and L25 HGT events. Strong associations were found between HGT in L2/L4 and the C2597T 23S rRNA mutation that confers reduced susceptibility to macrolides. Whilst previous studies have found evidence of HGT of entire genes coding for ribosomal proteins in other bacterial species, this is the first study to find evidence of HGT-mediated chimerization of ribosomal proteins.

Environmental Management of Legionella in Domestic Water Systems: Consolidated and Innovative Approaches for Disinfection Methods and Risk Assessment

Legionella is able to remain in water as free-living planktonic bacteria or to grow within biofilms that adhere to the pipes. It is also able to enter amoebas or to switch into a viable but not culturable (VBNC) state, which contributes to its resistance to harsh conditions and hinders its detection in water. Factors regulating Legionella growth, such as environmental conditions, type and concentration of available organic and inorganic nutrients, presence of protozoa, spatial location of microorganisms, metal plumbing components, and associated corrosion products are important for Legionella survival and growth. Finally, water treatment and distribution conditions may affect each of these factors. A deeper comprehension of Legionella interactions in water distribution systems with the environmental conditions is needed for better control of the colonization. To this purpose, the implementation of water management plans is the main prevention measure against Legionella. A water management program requires coordination among building managers, health care providers, and Public Health professionals. The review reports a comprehensive view of the state of the art and the promising perspectives of both monitoring and disinfection methods against Legionella in water, focusing on the main current challenges concerning the Public Health sector.

'Candidatus Sarmatiella mevalonica' endosymbiont of the ciliate Paramecium provides insights on evolutionary plasticity among Rickettsiales

Members of the bacterial order Rickettsiales are obligatorily associated with a wide range of eukaryotic hosts. Their evolutionary trajectories, in particular concerning the origin of shared or differential traits among distant sub-lineages, are still poorly understood. Here, we characterized a novel Rickettsiales bacterium associated with the ciliate Paramecium tredecaurelia and phylogenetically related to the Rickettsia genus. Its genome encodes significant lineage-specific features, chiefly the mevalonate pathway gene repertoire, involved in isoprenoid precursor biosynthesis. Not only this pathway has never been described in Rickettsiales, it also is very rare among bacteria, though typical in eukaryotes, thus likely representing a horizontally acquired trait. The presence of these genes could enable an efficient exploitation of host-derived intermediates for isoprenoid synthesis. Moreover, we hypothesize the reversed reactions could have replaced canonical pathways for producing acetyl-CoA, essential for phospholipid biosynthesis. Additionally, we detected phylogenetically unrelated mevalonate pathway genes in metagenome-derived Rickettsiales sequences, likely indicating evolutionary convergent effects of independent horizontal gene transfer events. Accordingly, convergence, involving both gene acquisitions and losses, is highlighted as a relevant evolutionary phenomenon in Rickettsiales, possibly favoured by plasticity and comparable lifestyles, representing a potentially hidden origin of other more nuanced similarities among sub-lineages.

Ecogenetics of antibiotic resistance in Listeria monocytogenes

The acquisition process of antibiotic resistance in an otherwise susceptible organism is shaped by the ecology of the species. Unlike other relevant human pathogens, Listeria monocytogenes has maintained a high rate of susceptibility to the antibiotics used for decades to treat human and animal infections. However, L. monocytogenes can acquire antibiotic resistance genes from other organisms' plasmids and conjugative transposons. Ecological factors could account for its susceptibility. L. monocytogenes is ubiquitous in nature, most frequently including reservoirs unexposed to antibiotics, including intracellular sanctuaries. L. monocytogenes has a remarkably closed genome, reflecting limited community interactions, small population sizes and high niche specialization. The L. monocytogenes species is divided into variants that are specialized in small specific niches, which reduces the possibility of coexistence with potential donors of antibiotic resistance. Interactions with potential donors are also hampered by interspecies antagonism. However, occasional increases in population sizes (and thus the possibility of acquiring antibiotic resistance) can derive from selection of the species based on intrinsic or acquired resistance to antibiotics, biocides, heavy metals or by a natural tolerance to extreme conditions. High-quality surveillance of the emergence of resistance to the key drugs used in primary therapy is mandatory.

Horizontal genetic exchange of chromosomally encoded markers between Campylobacter jejuni cells

Many epidemiological studies provide us with the evidence of horizontal gene transfer (HGT) contributing to the bacterial genomic diversity that benefits the bacterial populations with increased ability to adapt to the dynamic environments. Campylobacter jejuni, a major cause of acute enteritis in the U.S., often linked with severe post-infection neuropathies, has been reported to exhibit a non-clonal population structure and comparatively higher strain-level genetic variation. In this study, we provide evidence of the HGT of chromosomally encoded genetic markers between C. jejuni cells in the biphasic MH medium. We used two C. jejuni NCTC-11168 mutants harbouring distinct antibiotic-resistance genes [chloramphenicol (Cm) and kanamycin (Km)] present at two different neutral genomic loci. Cultures of both marker strains were mixed together and incubated for 5 hrs, then plated on MH agar plates supplemented with both antibiotics. The recombinant cells with double antibiotic markers were generated at the frequency of 0.02811 ± 0.0035% of the parental strains. PCR assays using locus-specific primers confirmed that transfer of the antibiotic-resistance genes was through homologous recombination. Also, the addition of chicken cecal content increased the recombination efficiency approximately up to 10-fold as compared to the biphasic MH medium (control) at P < 0.05. Furthermore, treating the co-culture with DNase I decreased the available DNA, which in turn significantly reduced recombination efficiency by 99.92% (P < 0.05). We used the cell-free supernatant of 16 hrs-culture of Wild-type C. jejuni as a template for PCR and found DNA sequences from six different genomic regions were easily amplified, indicating the presence of released chromosomal DNA in the culture supernatant. Our findings suggest that HGT in C. jejuni is facilitated in the chicken gut environment contributing to in vivo genomic diversity. Additionally, C. jejuni might have an active mechanism to release its chromosomal DNA into the extracellular environment, further expediting HGT in C. jejuni populations.

Evolution and function of bacterial RCC1 repeat effectors

Intracellular bacterial pathogens harbour genes, the closest homologues of which are found in eukaryotes. Regulator of chromosome condensation 1 (RCC1) repeat proteins are phylogenetically widespread and implicated in protein-protein interactions, such as the activation of the small GTPase Ran by its cognate guanine nucleotide exchange factor, RCC1. Legionella pneumophila and Coxiella burnetii, the causative agents of Legionnaires' disease and Q fever, respectively, harbour RCC1 repeat coding genes. Legionella pneumophila secretes the RCC1 repeat 'effector' proteins LegG1, PpgA and PieG into eukaryotic host cells, where they promote the activation of the pleiotropic small GTPase Ran, microtubule stabilisation, pathogen vacuole motility and intracellular bacterial growth as well as host cell migration. The RCC1 repeat effectors localise to the pathogen vacuole or the host plasma membrane and target distinct components of the Ran GTPase cycle, including Ran modulators and the small GTPase itself. Coxiella burnetii translocates the RCC1 repeat effector NopA into host cells, where the protein localises to nucleoli. NopA binds to Ran GTPase and promotes the nuclear accumulation of Ran(GTP), thus pertubing the import of the transcription factor NF-κB and innate immune signalling. Hence, divergent evolution of bacterial RCC1 repeat effectors defines the range of Ran GTPase cycle targets and likely allows fine-tuning of Ran GTPase activation by the pathogens at different cellular sites.

Legionella Diversity and Spatiotemporal Variation in The Occurrence of Opportunistic Pathogens within a Large Building Water System

Understanding Legionella survival mechanisms within building water systems (BWSs) is challenging due to varying engineering, operational, and water quality characteristics unique to each system. This study aimed to evaluate Legionella, mycobacteria, and free-living amoebae occurrence within a BWS over 18-28 months at six locations differing in plumbing material and potable water age, quality, and usage. A total of 114 bulk water and 57 biofilm samples were analyzed. Legionella culturability fluctuated seasonally with most culture-positive samples being collected during the winter compared to the spring, summer, and fall months. Positive and negative correlations between Legionella and L. pneumophila occurrence and other physiochemical and microbial water quality parameters varied between location and sample types. Whole genome sequencing of 19 presumptive Legionella isolates, from four locations across three time points, identified nine isolates as L. pneumophila serogroup (sg) 1 sequence-type (ST) 1; three as L. pneumophila sg5 ST1950 and ST2037; six as L. feeleii; and one as Ochrobactrum. Results showed the presence of a diverse Legionella population with consistent and sporadic occurrence at four and two locations, respectively. Viewed collectively with similar studies, this information will enable a better understanding of the engineering, operational, and water quality parameters supporting Legionella growth within BWSs.

Type I-F CRISPR-Cas Distribution and Array Dynamics in Legionella pneumophila

In bacteria and archaea, several distinct types of CRISPR-Cas systems provide adaptive immunity through broadly similar mechanisms: short nucleic acid sequences derived from foreign DNA, known as spacers, engage in complementary base pairing with invasive genetic elements setting the stage for nucleases to degrade the target DNA. A hallmark of type I CRISPR-Cas systems is their ability to acquire spacers in response to both new and previously encountered invaders (naïve and primed acquisition, respectively). Our phylogenetic analyses of 43 L. pneumophila type I-F CRISPR-Cas systems and their resident genomes suggest that many of these systems have been horizontally acquired. These systems are frequently encoded on plasmids and can co-occur with nearly identical chromosomal loci. We show that two such co-occurring systems are highly protective and undergo efficient primed acquisition in the lab. Furthermore, we observe that targeting by one system’s array can prime spacer acquisition in the other. Lastly, we provide experimental and genomic evidence for a model in which primed acquisition can efficiently replenish a depleted type I CRISPR array following a mass spacer deletion event.

Gene Transmission in the One Health Microbiosphere and the Channels of Antimicrobial Resistance

Antibiotic resistance is a field in which the concept of One Health can best be illustrated. One Health is based on the definition of communication spaces among diverse environments. Antibiotic resistance is encoded by genes, however, these genes are propagated in mobile genetic elements (MGEs), circulating among bacterial species and clones that are integrated into the multiple microbiotas of humans, animals, food, sewage, soil, and water environments, the One Health microbiosphere. The dynamics and evolution of antibiotic resistance depend on the communication networks linking all these ecological, biological, and genetic entities. These communications occur by environmental overlapping and merging, a critical issue in countries with poor sanitation, but also favored by the homogenizing power of globalization. The overwhelming increase in the population of highly uniform food animals has contributed to the parallel increase in the absolute size of their microbiotas, consequently enhancing the possibility of microbiome merging between humans and animals. Microbial communities coalescence might lead to shared microbiomes in which the spread of antibiotic resistance (of human, animal, or environmental origin) is facilitated. Intermicrobiome communication is exerted by shuttle bacterial species (or clones within species) belonging to generalist taxa, able to multiply in the microbiomes of various hosts, including humans, animals, and plants. Their integration into local genetic exchange communities fosters antibiotic resistance gene flow, following the channels of accessory genome exchange among bacterial species. These channels delineate a topology of gene circulation, including dense clusters of species with frequent historical and recent exchanges. The ecological compatibility of these species, sharing the same niches and environments, determines the exchange possibilities. In summary, the fertility of the One Health approach to antibiotic resistance depends on the progress of understanding multihierarchical systems, encompassing communications among environments (macro/microaggregates), among microbiotas (communities), among bacterial species (clones), and communications among MGEs.

Diverse conjugative elements silence natural transformation in Legionella species

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

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

Intracellular parasitism, the driving force of evolution of Legionella pneumophila and the genus Legionella

Legionella pneumophila is an intracellular pathogen that causes a severe pneumonia called Legionnaires' disease that is often fatal when not promptly diagnosed and treated. Legionella parasitize aquatic protozoa with which it co-evolved over an evolutionary long time. The close relationship between hosts and pathogens, their co-evolution, led to molecular interactions such as the exchange of genetic material through horizontal gene transfer (HGT). Genome sequencing of L. pneumophila and of the entire genus Legionella that comprises over 60 species revealed that Legionellae have co-opted genes and thus cellular functions from their eukaryotic hosts to a surprisingly high extent. Acquisition and loss of these eukaryotic-like genes and domains is an on-going process underlining the highly dynamic nature of the Legionella genomes. Although the large amount and diversity of HGT in Legionella seems to be unique in the prokaryotic world the analyses of more and more genomes from environmental organisms and symbionts of amoeba revealed that such genetic exchanges occur among all amoeba associated bacteria and also among the different microorganisms that infect amoeba. This dynamic reshuffling and gene-acquisition has led to the emergence of Legionella as human pathogen and may lead to the emergence of new human pathogens from the environment.

Intracellular parasitism, the driving force of evolution of Legionella pneumophila and the genus Legionella

Legionella pneumophila is an intracellular pathogen that causes a severe pneumonia called Legionnaires' disease that is often fatal when not promptly diagnosed and treated. However, L. pneumophila is mainly an environmental pathogen of protozoa. This bacterium parasitizes free-living amoeba and other aquatic protozoa with which it co-evolved over an evolutionary long time. Due to the close relationship between hosts and pathogens, their co-evolution leads to molecular interactions such as the exchange of genetic material through horizontal gene transfer (HGT). Those genes that confer an advantage to the bacteria were fixed in their genomes and help these pathogens to subvert host functions to their advantage. Genome sequencing of L. pneumophila and recently of the entire genus Legionella that comprises over 60 species revealed that Legionellae have co-opted genes and thus cellular functions from their eukaryotic hosts to a surprisingly high extent never observed before for an prokaryotic organism. Acquisition and loss of these eukaryotic-like genes and eukaryotic domains is an ongoing process underlining the highly dynamic nature of the Legionella genomes. Although the large amount and diversity of HGT that occurred between Legionella and their protozoan hosts seems to be unique in the prokaryotic world, the analyses of more and more genomes from environmental organisms and symbionts of amoeba revealed that such genetic exchanges occur among all amoeba-associated bacteria and also among the different microorganisms that infect amoeba such as viruses. This dynamic reshuffling and gene-acquisition has led to the emergence of major human pathogens such as Legionella and may lead to the emergence of new human pathogens from the environment.

Identification of MazF Homologue in Legionella pneumophila Which Cleaves RNA at the AACU Sequence

MazF is a sequence-specific endoribonuclease that is widely conserved in bacteria and archaea. Here, we found an MazF homologue (MazF-lp; LPO-p0114) in Legionella pneumophila. The mazF-lp gene overlaps 14 base pairs with the upstream gene mazE-lp (MazE-lp; LPO-p0115). The induction of mazF-lp caused cell growth arrest, while mazE-lp co-induction recovered cell growth in Escherichia coli. In vivo and in vitro primer extension experiments showed that MazF-lp is a sequence-specific endoribonuclease cleaving RNA at AACU. The endoribonuclease activity of purified MazF-lp was inhibited by purified MazE-lp. We found that MazE-lp and the MazEF-lp complex specifically bind to the palindromic sequence present in the 5'-untranslated region of the mazEF-lp operon. MazE-lp and MazEF-lp both likely function as a repressor for the mazEF-lp operon and for other genes, including icmR, whose gene product functions as a secretion chaperone for the IcmQ pore-forming protein, by specifically binding to the palindromic sequence in 5'-UTR of these genes.

The pleiotropic Legionella transcription factor LvbR links the Lqs and c-di-GMP regulatory networks to control biofilm architecture and virulence

The causative agent of Legionnaires' disease, Legionella pneumophila, colonizes amoebae and biofilms in the environment. The opportunistic pathogen employs the Lqs (Legionella quorum sensing) system and the signalling molecule LAI-1 (Legionella autoinducer-1) to regulate virulence, motility, natural competence and expression of a 133 kb genomic "fitness island", including a putative novel regulator. Here, we show that the regulator termed LvbR is an LqsS-regulated transcription factor that binds to the promoter of lpg1056/hnox1 (encoding an inhibitor of the diguanylate cyclase Lpg1057), and thus, regulates proteins involved in c-di-GMP metabolism. LvbR determines biofilm architecture, since L. pneumophila lacking lvbR accumulates less sessile biomass and forms homogeneous mat-like structures, while the parental strain develops more compact bacterial aggregates. Comparative transcriptomics of sessile and planktonic ΔlvbR or ΔlqsR mutant strains revealed concerted (virulence, fitness island, metabolism) and reciprocally (motility) regulated genes in biofilm and broth respectively. Moreover, ΔlvbR is hyper-competent for DNA uptake, defective for phagocyte infection, outcompeted by the parental strain in amoebae co-infections and impaired for cell migration inhibition. Taken together, our results indicate that L. pneumophila LvbR is a novel pleiotropic transcription factor, which links the Lqs and c-di-GMP regulatory networks to control biofilm architecture and pathogen-host cell interactions.

The Legionella effector LtpM is a new type of phosphoinositide-activated glucosyltransferase

Legionella pneumophila causes Legionnaires' disease, a severe form of pneumonia. L. pneumophila translocates more than 300 effectors into host cells via its Dot/Icm (Defective in organelle trafficking/Intracellular multiplication) type IV secretion system to enable its replication in target cells. Here, we studied the effector LtpM, which is encoded in a recombination hot spot in L. pneumophila Paris. We show that a C-terminal phosphoinositol 3-phosphate (PI3P)-binding domain, also found in otherwise unrelated effectors, targets LtpM to the Legionella-containing vacuole and to early and late endosomes. LtpM expression in yeast caused cytotoxicity. Sequence comparison and structural homology modeling of the N-terminal domain of LtpM uncovered a remote similarity to the glycosyltransferase (GT) toxin PaTox from the bacterium Photorhabdus asymbiotica; however, instead of the canonical DxD motif of GT-A type glycosyltransferases, essential for enzyme activity and divalent cation coordination, we found that a DxN motif is present in LtpM. Using UDP-glucose as sugar donor, we show that purified LtpM nevertheless exhibits glucohydrolase and autoglucosylation activity in vitro and demonstrate that PI3P binding activates LtpM's glucosyltransferase activity toward protein substrates. Substitution of the aspartate or the asparagine in the DxN motif abolished the activity of LtpM. Moreover, whereas all glycosyltransferase toxins and effectors identified so far depend on the presence of divalent cations, LtpM is active in their absence. Proteins containing LtpM-like GT domains are encoded in the genomes of other L. pneumophila isolates and species, suggesting that LtpM is the first member of a novel family of glycosyltransferase effectors employed to subvert hosts.

More than 18,000 effectors in the Legionella genus genome provide multiple, independent combinations for replication in human cells

The genus Legionella comprises 65 species, among which Legionella pneumophila is a human pathogen causing severe pneumonia. To understand the evolution of an environmental to an accidental human pathogen, we have functionally analyzed 80 Legionella genomes spanning 58 species. Uniquely, an immense repository of 18,000 secreted proteins encoding 137 different eukaryotic-like domains and over 200 eukaryotic-like proteins is paired with a highly conserved type IV secretion system (T4SS). Specifically, we show that eukaryotic Rho- and Rab-GTPase domains are found nearly exclusively in eukaryotes and Legionella Translocation assays for selected Rab-GTPase proteins revealed that they are indeed T4SS secreted substrates. Furthermore, F-box, U-box, and SET domains were present in >70% of all species, suggesting that manipulation of host signal transduction, protein turnover, and chromatin modification pathways are fundamental intracellular replication strategies for legionellae. In contrast, the Sec-7 domain was restricted to L. pneumophila and seven other species, indicating effector repertoire tailoring within different amoebae. Functional screening of 47 species revealed 60% were competent for intracellular replication in THP-1 cells, but interestingly, this phenotype was associated with diverse effector assemblages. These data, combined with evolutionary analysis, indicate that the capacity to infect eukaryotic cells has been acquired independently many times within the genus and that a highly conserved yet versatile T4SS secretes an exceptional number of different proteins shaped by interdomain gene transfer. Furthermore, we revealed the surprising extent to which legionellae have coopted genes and thus cellular functions from their eukaryotic hosts, providing an understanding of how dynamic reshuffling and gene acquisition have led to the emergence of major human pathogens.

Evasion of phagotrophic predation by protist hosts and innate immunity of metazoan hosts by Legionella pneumophila

Legionella pneumophila is a ubiquitous environmental bacterium that has evolved to infect and proliferate within amoebae and other protists. It is thought that accidental inhalation of contaminated water particles by humans is what has enabled this pathogen to proliferate within alveolar macrophages and cause pneumonia. However, the highly evolved macrophages are equipped with more sophisticated innate defence mechanisms than are protists, such as the evolution of phagotrophic feeding into phagocytosis with more evolved innate defence processes. Not surprisingly, the majority of proteins involved in phagosome biogenesis (~80%) have origins in the phagotrophy stage of evolution. There are a plethora of highly evolved cellular and innate metazoan processes, not represented in protist biology, that are modulated by L. pneumophila, including TLR2 signalling, NF-κB, apoptotic and inflammatory processes, histone modification, caspases, and the NLRC-Naip5 inflammasomes. Importantly, L. pneumophila infects haemocytes of the invertebrate Galleria mellonella, kill G. mellonella larvae, and proliferate in and kill Drosophila adult flies and Caenorhabditis elegans. Although coevolution with protist hosts has provided a substantial blueprint for L. pneumophila to infect macrophages, we discuss the further evolutionary aspects of coevolution of L. pneumophila and its adaptation to modulate various highly evolved innate metazoan processes prior to becoming a human pathogen.

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

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

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.

Multiplex polymerase chain reaction of genetic markers for detection of potentially pathogenic environmental Legionella pneumophila isolates

Background & objectives:

Genomic constitution of the bacterium Legionella pneumophila plays an important role in providing them a pathogenic potential. Here, we report the standardization and application of multiplex polymerase chain reaction (PCR) for the detection of molecular markers of pathogenic potential in L. pneumophila in hospital environment.

Methods:

Culture of the standard strains of L. pneumophila was performed in buffered charcoal-yeast extract agar with L-cysteine at pH 6.9. Primers were designed for multiplex PCR, and standardization for the detection of five markers annotated to L. pneumophila plasmid pLPP (11A2), lipopolysaccharide synthesis (19H4), CMP-N-acetylneuraminic acid synthetase (10B12), conjugative coupling factor (24B1) and hypothetical protein (8D6) was done. A total of 195 water samples and 200 swabs were collected from the hospital environment. The bacterium was isolated from the hospital environment by culture and confirmed by 16S rRNA gene PCR and restriction enzyme analysis. A total of 45 L. pneumophila isolates were studied using the standardized multiplex PCR.

Results:

The PCR was sensitive to detect 0.1 ng/μl DNA and specific for the two standard strains used in the study. Of the 45 hospital isolates tested, 11 isolates had four markers, 12 isolates had three markers, 10 isolates had two markers, nine isolates had one marker and three isolates had none of the markers. None of the isolates had all the five markers.

Interpretation & conclusions:

The findings of this study showed the presence of gene markers of pathogenic potential of the bacterium L. pneumophila. However, the genomic constitution of the environmental isolates should be correlated with clinical isolates to prove their pathogenic potential. Rapid diagnostic methods such as multiplex PCR reported here, for elucidating gene markers, could help in future epidemiological studies of bacterium L. pneumophila.

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

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

Mammalian Solute Carrier (SLC)-like transporters of Legionella pneumophila

Acquisition of nutrients during intra-vacuolar growth of L. pneumophila within macrophages or amoebae is poorly understood. Since many genes of L. pneumophila are acquired by inter-kingdom horizontal gene transfer from eukaryotic hosts, we examined the presence of human solute carrier (SLC)-like transporters in the L. pneumophila genome using I-TASSER to assess structural alignments. We identified 11 SLC-like putative transporters in L. pneumophila that are structurally similar to SLCs, eight of which are amino acid transporters, and one is a tricarboxylate transporter. The two other transporters, LstA and LstB, are structurally similar to the human glucose transporter, SLC2a1/Glut1. Single mutants of lstA or lstB have decreased ability to import, while the lstA/lstB double mutant is severely defective for uptake of glucose. While lstA or lstB single mutants are not defective in intracellular proliferation within Acanthamoeba polyphaga and human monocyte-derived macrophages, the lstA/lstB double mutant is severely defective in both host cells. The two phenotypic defects of the lstA/lstB double mutant in uptake of glucose and intracellular replication are both restored upon complementation of either lstA or lstB. Our data show that the two glucose transporters, LstA and LstB, are redundant and are required for intracellular replication within human macrophages and amoebae.

A Legionella pneumophila amylase is essential for intracellular replication in human macrophages and amoebae

Legionella pneumophila invades protozoa with an "accidental" ability to cause pneumonia upon transmission to humans. To support its nutrition during intracellular residence, L. pneumophila relies on host amino acids as the main source of carbon and energy to feed the TCA cycle. Despite the apparent lack of a requirement for glucose for L. pneumophila growth in vitro and intracellularly, the organism contains multiple amylases, which hydrolyze polysaccharides into glucose monomers. Here we describe one predicted putative amylase, LamB, which is uniquely present only in L. pneumophila and L. steigerwaltii among the ~60 species of Legionella. Our data show that LamB has a strong amylase activity, which is abolished upon substitutions of amino acids that are conserved in the catalytic pocket of amylases. Loss of LamB or expression of catalytically-inactive variants of LamB results in a severe growth defect of L. pneumophila in Acanthamoeba polyphaga and human monocytes-derived macrophages. Importantly, the lamB null mutant is severely attenuated in intra-pulmonary proliferation in the mouse model and is defective in dissemination to the liver and spleen. Our data show an essential role for LamB in intracellular replication of L. pneumophila in amoeba and human macrophages and in virulence in vivo.

Population structure of Environmental and Clinical Legionella pneumophila isolates in Catalonia

Legionella is the causative agent of Legionnaires’ disease (LD). In Spain, Catalonia is the region with the highest incidence of LD cases. The characterisation of clinical and environmental isolates using molecular epidemiology techniques provides epidemiological data for a specific geographic region and makes it possible to carry out phylogenetic and population-based analyses. The aim of this study was to describe and compare environmental and clinical isolates of Legionella pneumophila in Catalonia using sequence-based typing and monoclonal antibody subgrouping. A total of 528 isolates were characterised. For data analysis, the isolates were filtered to reduce redundancies, and 266 isolates (109 clinical and 157 environmental) were finally included. Thirty-two per cent of the clinical isolates were ST23, ST37 and ST1 while 40% of the environmental isolates were ST284 and ST1. Although the index of diversity was higher in clinical than in environmental ST isolates, we observed that clinical STs were similar to those recorded in other regions but that environmental STs were more confined to particular study areas. This observation supports the idea that only certain STs trigger cases or outbreaks in humans. Therefore, comparison of the genomes of clinical and environmental isolates could provide important information about the traits that favour infection or environmental persistence.

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.

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

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

Molecular typing of Legionella pneumophila isolates from environmental water samples and clinical samples using a five-gene sequence typing and standard Sequence-Based Typing

Inadequate discriminatory power to distinguish between L. pneumophila isolates, especially those belonging to disease-related prevalent sequence types (STs) such as ST1, ST36 and ST47, is an issue of SBT scheme. In this study, we developed a multilocus sequence typing (MLST) scheme based on two non-virulence loci (trpA, cca) and three virulence loci (icmK, lspE, lssD), to genotype 110 L. pneumophila isolates from various natural and artificial water sources in Guangdong province of China, and compared with the SBT. The isolates were assigned to 33 STs of the SBT and 91 new sequence types (nSTs) of the MLST. The indices of discrimination (IODs) of SBT and MLST were 0.920 and 0.985, respectively. Maximum likelihood trees of the concatenated SBT and MLST sequences both showed distinct phylogenetic relationships between the isolates from the two environments. More intragenic recombinations were detected in nSTs than in STs, and they were both more abundant in natural water isolates. We found out the MLST had a high discriminatory ability for the disease-associated ST1 isolates: 22 ST1 isolates were assigned to 19 nSTs. Furthermore, we assayed the discrimination of the MLST for 29 reference strains (19 clinical and 10 environmental). The clinical strains were assigned to eight STs and ten nSTs. The MLST could also subtype the prevalent clinical ST36 or ST47 strains: eight ST36 strains were subtyped into three nSTs and two ST47 strains were subtyped into two nSTs. We found different distribution patterns of nSTs between the environmental and clinical ST36 isolates, and between the outbreak clinical ST36 isolates and the sporadic clinical ST36 isolates. These results together revealed the MLST scheme could be used as part of a typing scheme that increased discrimination when necessary.

Molecular evolution of virulence genes and non-virulence genes in clinical, natural and artificial environmental Legionella pneumophila isolates

Background

L. pneumophila is the main causative agent of Legionnaires’ disease. Free-living amoeba in natural aquatic environments is the reservoir and shelter for L. pneumophila. From natural water sources, L. pneumophila can colonize artificial environments such as cooling towers and hot-water systems, and then spread in aerosols, infecting the susceptible person. Therefore, molecular phylogeny and genetic variability of L. pneumophila from different sources (natural water, artificial water, and human lung tissue) might be distinct because of the selection pressure in different environments. Several studies researched genetic differences between L. pneumophila clinical isolates and environmental isolates at the nucleotide sequence level. These reports mainly focused on the analysis of virulence genes, and rarely distinguished artificial and natural isolates.

Methods

We have used 139 L. pneumophila isolates to study their genetic variability and molecular phylogeny. These isolates include 51 artificial isolates, 59 natural isolates, and 29 clinical isolates. The nucleotide sequences of two representative non-virulence (NV) genes (trpA, cca) and three representative virulence genes (icmK, lspE, lssD) were obtained using PCR and DNA sequencing and were analyzed.

Results

Levels of genetic variability including haplotypes, haplotype diversity, nucleotide diversity, nucleotide difference and the total number of mutations in the virulence loci were higher in the natural isolates. In contrast, levels of genetic variability including polymorphic sites, theta from polymorphic sites and the total number of mutations in the NV loci were higher in clinical isolates. A phylogenetic analysis of each individual gene tree showed three to six main groups, but not comprising the same L. pneumophila isolates. We detected recombination events in every virulence loci of natural isolates, but only detected them in the cca locus of clinical isolates. Neutrality tests showed that variations in the virulence genes of clinical and environmental isolates were under neutral evolution. TrpA and cca loci of clinical isolates showed significantly negative values of Tajima’s D, Fu and Li’s D* and F*, suggesting the presence of negative selection in NV genes of clinical isolates.

Discussion

Our findingsreinforced the point that the natural environments were the primary training place for L. pneumophila virulence, and intragenic recombination was an important strategy in the adaptive evolution of virulence gene. Our study also suggested the selection pressure had unevenly affected these genes and contributed to the different evolutionary patterns existed between NV genes and virulence genes. This work provides clues for future work on population-level and genetics-level questions about ecology and molecular evolution of L. pneumophila, as well as genetic differences of NV genes and virulence genes between this host-range pathogen with different lifestyles.

Serine/Threonine Protein Kinases from Bacteria, Archaea and Eukarya Share a Common Evolutionary Origin Deeply Rooted in the Tree of Life

The main family of serine/threonine/tyrosine protein kinases present in eukarya was defined and described by Hanks et al. in 1988 (Science, 241, 42-52). It was initially believed that these kinases do not exist in bacteria, but extensive genome sequencing revealed their existence in many bacteria. For historical reasons, the term "eukaryotic-type kinases" propagated in the literature to describe bacterial members of this protein family. Here, we argue that this term should be abandoned as a misnomer, and we provide several lines of evidence to support this claim. Our comprehensive phylostratigraphic analysis suggests that Hanks-type kinases present in eukarya, bacteria and archaea all share a common evolutionary origin in the lineage leading to the last universal common ancestor (LUCA). We found no evidence to suggest substantial horizontal transfer of genes encoding Hanks-type kinases from eukarya to bacteria. Moreover, our systematic structural comparison suggests that bacterial Hanks-type kinases resemble their eukaryal counterparts very closely, while their structures appear to be dissimilar from other kinase families of bacterial origin. This indicates that a convergent evolution scenario, by which bacterial kinases could have evolved a kinase domain similar to that of eukaryal Hanks-type kinases, is not very likely. Overall, our results strongly support a monophyletic origin of all Hanks-type kinases, and we therefore propose that this term should be adopted as a universal name for this protein family.

Benefits of Genomic Insights and CRISPR-Cas Signatures to Monitor Potential Pathogens across Drinking Water Production and Distribution Systems

The occurrence of pathogenic bacteria in drinking water distribution systems (DWDSs) is a major health concern, and our current understanding is mostly related to pathogenic species such as Legionella pneumophila and Mycobacterium avium but not to bacterial species closely related to them. In this study, genomic-based approaches were used to characterize pathogen-related species in relation to their abundance, diversity, potential pathogenicity, genetic exchange, and distribution across an urban drinking water system. Nine draft genomes recovered from 10 metagenomes were identified as Legionella (4 draft genomes), Mycobacterium (3 draft genomes), Parachlamydia (1 draft genome), and Leptospira (1 draft genome). The pathogenicity potential of these genomes was examined by the presence/absence of virulence machinery, including genes belonging to Type III, IV, and VII secretion systems and their effectors. Several virulence factors known to pathogenic species were detected with these retrieved draft genomes except the Leptospira-related genome. Identical clustered regularly interspaced short palindromic repeats-CRISPR-associated proteins (CRISPR-Cas) genetic signatures were observed in two draft genomes recovered at different stages of the studied system, suggesting that the spacers in CRISPR-Cas could potentially be used as a biomarker in the monitoring of Legionella related strains at an evolutionary scale of several years across different drinking water production and distribution systems. Overall, metagenomics approach was an effective and complementary tool of culturing techniques to gain insights into the pathogenic characteristics and the CRISPR-Cas signatures of pathogen-related species in DWDSs.