Virulence factor rtx in Legionella pneumophila, evidence suggesting it is a modular multifunctional protein

A rapid and reliable method for early Legionella pneumophila identification and characterization in support of the epidemiology study

Introduction

Legionnaires' disease is a severe pneumonia predominantly caused by Legionella pneumophila (Lp), whose major reservoirs are artificial water systems. As most human infections are caused by L. pneumophila serogroup 1 (Lp1), a reliable method for Lp distinction can be crucial for bacterial spread prevention. As the ability to withstand in environments and to cause the waterborne disease is strongly related to specific genes, the identification of virulent strains can be of great relevance to implement water environmental monitoring and to contain harmful outbreaks to public health. We aimed to test an assay for Lp identification among different Legionella species, and to determine the serogroups. Additionally, we investigated the carriage of virulence and antimicrobial resistance genes.

Methods

A total of 90 Legionella spp. isolates identified by phenotypic tests were subjected to the designed quantitative PCR assay targeting specific mip for Lp, wzm for Lp1, pvcA and ahpD for biofilm production. Eleven serogroups were investigated in all our isolates tested positive for mip gene, subsequently analyzed for 12 virulence and 8 antimicrobial resistance genes.

Results

Only the 70 Lp isolates were positive for mip. Out of 27 Lp isolates belonging to serogroup 1 based on agglutination test, 23 (85.2%) carried wzm. The presence of ahpD and pvcA was found in 94.3 and 98.6% of Lp isolates, respectively. By multiplex PCR, all 23 wzm-positive strains were confirmed as serogroup 1 that was the most predominant (33%). At least one virulence gene was detected in all Lp isolates. The most frequent gene was ispE (100%), followed by issD (96%), icmK and enhC (93%), cpxA (91%), rtxA2 (74%), lvhB8-B9 (61%), and prpA (54%). The other genes were less diffused in Lp strains (rtxA1, 44%; lvhB3-B4, 47%; pvcB, 27%; lvrE, 24%). Of the macrolide resistance genes, the ereA was found in 84% of Lp strains, while only 14 (20%) harbored the lpeAB among the efflux pump genes.

Conclusion

The assays validated in this study enable the simultaneous Lp and Lp1 detection. The differentiation of Lp strains according to their virulence properties could be useful to predict the bacterial ability to survive and to cause the disease.

Novel induction of broad-spectrum antibiotics by the human pathogen Legionella

The majority of antibiotics are natural products, with microorganism-generated molecules and their derivatives being the most prevalent source of drugs to treat infections. Thus, identifying natural products remains the most valuable resource for novel therapeutics. Here, we report the discovery of a series of dormant bacteria in honey that have bactericidal activity toward Legionella, a bacterial pathogen that causes respiratory disease in humans. We show that, in response to bacterial products secreted by Legionella, the honey bacteria release diffusible antimicrobial molecules. Remarkably, the honey bacteria only produce these molecules in response to Legionella spp., when compared to a panel of 24 bacterial pathogens from different genera. However, the molecules induced by Legionella have broad activity against several clinically important pathogens, including many high-priority pathogens. Thus, Legionella spp. are potent drivers of antimicrobial molecule production by uncharacterized bacteria isolated from honey, providing access to new antimicrobial products and an unprecedented strategy for discovering novel antibiotics.

IMPORTANCE

Natural products generated by microorganisms remain the most viable and abundant source of new antibiotics. However, their discovery depends on the ability to isolate and culture the producing organisms and to identify conditions that promote antibiotic production. Here, we identify a series of previously undescribed bacteria isolated from raw honey and specific culture conditions that induce the production of antimicrobial molecules that are active against a wide variety of pathogenic bacteria.

Legionella pneumophila cell surface RtxA release by LapD/LapG and its role in virulence

BACKGROUND

Legionella pneumophila is a Gram-negative intracellular bacillus and is the causative agent of a severe form of pneumonia called Legionnaires' disease which accounts for 2-9% of cases of community acquired pneumonia. It produces an extremely large protein belonging to the RTX (Repeats in ToXin) family, called RtxA, and we previously reported that RtxA is transported by a dedicated type 1 secretion system (T1SS) to the cell surface. RTX proteins have been shown to participate in the virulence or biofilm formation of various bacteria, the most studied models being the pore forming hemolysin A (HlyA) of Escherichia coli and the biofilm associated protein LapA of P. fluorescens. LapA localization depends on the enzymatic release by LapD/LapG complex activity. This study aimed to elucidate the dual localization (cell surface associated or released state) of L. pneumophila RTX protein (RtxA) and whether this released versus sequestered state of RtxA plays a role in L. pneumophila virulence.

RESULTS

The hereby work reveals that, in vitro, LapG periplasmic protease cleaves RtxA N-terminus in the middle of a di-alanine motif (position 108-109). Consistently, a strain lacking LapG protease maintains RtxA on the cell surface, whereas a strain lacking the c-di-GMP receptor LapD does not exhibit cell surface RtxA because of its continuous cleavage and release, as in the LapA-D-G model of Pseudomonas fluorescens. Interestingly, our data point out a key role of RtxA in enhancing the infection process of amoeba cells, regardless of its location (embedded or released); therefore, this may be the result of a secondary role of this surface protein.

CONCLUSIONS

This is the first experimental identification of the cleavage site within the RTX protein family. The primary role of RtxA in Legionella is still questionable as in many other bacterial species, hence it sounds reasonable to propose a major function in biofilm formation, promoting cell aggregation when RtxA is embedded in the outer membrane and facilitating biofilm dispersion in case of RtxA release. The role of RtxA in enhancing the infection process may be a result of its action on host cells (i.e., PDI interaction or pore-formation), and independently of its status (embedded or released).

Virulence Genes and Biofilm Formation Among Legionella pneumophila Isolates Collected from Hospital Water Sources

Legionella pneumophila can be transmitted to people, especially immunocompromised patients, via hospital water pipe systems and cause severe pneumonia. The aim of our study was to investigate the presence of major virulence factor genes, ability of biofilms formation, and correlation between presence of Legionella isolates and temperature, pH, and residual chlorine of water. Hundred water samples were collected from nine hospitals in Tehran, Iran. Temperature, pH, and residual chlorine were determined during sampling. Different virulence genes and the ability to form biofilms were subsequently analyzed among the L. pneumophila isolates. Results showed that 12 (12%) samples were positive in culture method and all of the isolates were positive as L. pneumophila species (mip). A correlation was found between Legionella culture positivity and temperature and pH of water, but there was no significant correlation between residual chlorine of water samples and the presence of Legionella. The isolation of Legionella rate in summer and spring was higher than winter and autumn. Twelve (100%) isolates were positive for mip genes, 9 (75%) for dot genes, 8 (66.66%) for hsp, 6 (50%) for lvh, and 4 (33.33%) for rtx. All of the isolates displayed strong ability for biofilm production every three days. Two of these isolates (16.6%) displayed weak ability to form biofilm on the first day of incubation. This study revealed that water sources in hospitals were colonized by virulent Legionella and should be continuously monitored to avoid elevated concentrations of Legionella with visible biofilm formation.

Evaluation of Genetic Diversity and Virulence Potential of Legionella pneumophila Isolated from Water Supply Systems of Residential Buildings in Latvia

Legionella is an opportunistic pathogen with a biphasic life cycle that occasionally infects humans. The aim of the study was to assess the distribution of virulence genes and genetic diversity among L. pneumophila isolated from water supply systems of residential buildings in Latvia. In total, 492 water samples from 200 residential buildings were collected. Identification of Legionella spp. was performed according to ISO 11731, and 58 isolates were subjected to whole-genome sequencing. At least one Legionella-positive sample was found in 112 out of 200 apartment buildings (56.0%). The study revealed extensive sequence-type diversity, where 58 L. pneumophila isolates fell into 36 different sequence types. A total of 420 virulence genes were identified, of which 260 genes were found in all sequenced L. pneumophila isolates. The virulence genes enhC, htpB, omp28, and mip were detected in all isolates, suggesting that adhesion, attachment, and entry into host cells are enabled for all isolates. The relative frequency of virulence genes among L. pneumophila isolates was high. The high prevalence, extensive genetic diversity, and the wide range of virulence genes indicated that the virulence potential of environmental Legionella is high, and proper risk management is of key importance to public health.

The inside scoop: Comparative genomics of two intranuclear bacteria, "Candidatus Berkiella cookevillensis" and "Candidatus Berkiella aquae"

"Candidatus Berkiella cookevillensis" (strain CC99) and "Candidatus Berkiella aquae" (strain HT99), belonging to the Coxiellaceae family, are gram-negative bacteria isolated from amoebae in biofilms present in human-constructed water systems. Both bacteria are obligately intracellular, requiring host cells for growth and replication. The intracellular bacteria-containing vacuoles of both bacteria closely associate with or enter the nuclei of their host cells. In this study, we analyzed the genome sequences of CC99 and HT99 to better understand their biology and intracellular lifestyles. The CC99 genome has a size of 2.9Mb (37.9% GC) and contains 2,651 protein-encoding genes (PEGs) while the HT99 genome has a size of 3.6Mb (39.4% GC) and contains 3,238 PEGs. Both bacteria encode high proportions of hypothetical proteins (CC99: 46.5%; HT99: 51.3%). The central metabolic pathways of both bacteria appear largely intact. Genes for enzymes involved in the glycolytic pathway, the non-oxidative branch of the phosphate pathway, the tricarboxylic acid pathway, and the respiratory chain were present. Both bacteria, however, are missing genes for the synthesis of several amino acids, suggesting reliance on their host for amino acids and intermediates. Genes for type I and type IV (dot/icm) secretion systems as well as type IV pili were identified in both bacteria. Moreover, both bacteria contain genes encoding large numbers of putative effector proteins, including several with eukaryotic-like domains such as, ankyrin repeats, tetratricopeptide repeats, and leucine-rich repeats, characteristic of other intracellular bacteria.

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.

Overview of the Clinical and Molecular Features of Legionella Pneumophila: Focus on Novel Surveillance and Diagnostic Strategies

Legionella pneumophila (L. pneumophila) is one of the most threatening nosocomial pathogens. The implementation of novel and more effective surveillance and diagnostic strategies is mandatory to prevent the occurrence of legionellosis outbreaks in hospital environments. On these bases, the present review is aimed to describe the main clinical and molecular features of L. pneumophila focusing attention on the latest findings on drug resistance mechanisms. In addition, a detailed description of the current guidelines for the disinfection and surveillance of the water systems is also provided. Finally, the diagnostic strategies available for the detection of Legionella spp. were critically reviewed, paying the attention to the description of the culture, serological and molecular methods as well as on the novel high-sensitive nucleic acid amplification systems, such as droplet digital PCR.

Computational prediction of secreted proteins in gram-negative bacteria

Gram-negative bacteria harness multiple protein secretion systems and secrete a large proportion of the proteome. Proteins can be exported to periplasmic space, integrated into membrane, transported into extracellular milieu, or translocated into cytoplasm of contacting cells. It is important for accurate, genome-wide annotation of the secreted proteins and their secretion pathways. In this review, we systematically classified the secreted proteins according to the types of secretion systems in Gram-negative bacteria, summarized the known features of these proteins, and reviewed the algorithms and tools for their prediction.

Growth by Insertion: The Family of Bacterial DDxP Proteins

We have identified a variety of proteins in species of the Legionella, Aeromonas, Pseudomonas, Vibrio, Nitrosomonas, Nitrosospira, Variovorax, Halomonas, and Rhizobia genera, which feature repetitive modules of different length and composition, invariably ending at the COOH side with Asp-Asp-x-Pro (DDxP) motifs. DDxP proteins range in size from 900 to 6200 aa (amino acids), and contain 1 to 5 different module types, present in one or multiple copies. We hypothesize that DDxP proteins were modeled by the action of specific endonucleases inserting DNA segments into genes encoding DDxP motifs. Target site duplications (TSDs) formed upon repair of staggered ends generated by endonuclease cleavage would explain the DDxP motifs at repeat ends. TSDs acted eventually as targets for the insertion of more modules of the same or different types. Repeat clusters plausibly resulted from amplification of both repeat and flanking TSDs. The proposed growth shown by the insertion model is supported by the identification of homologous proteins lacking repeats in Pseudomonas and Rhizobium. The 85 DDxP repeats identified in this work vary in length, and can be sorted into short (136-215 aa) and long (243-304 aa) types. Conserved Asp-Gly-Asp-Gly-Asp motifs are located 11-19 aa from the terminal DDxP motifs in all repeats, and far upstream in most long repeats.

Distribution of Virulence Genes and Sequence-Based Types Among Legionella pneumophila Isolated From the Water Systems of a Tertiary Care Hospital in India

Background: Legionnaires' disease (LD) is a potentially fatal pneumonia predominantly caused by infection due to Legionella pneumophila although more than 50 other Legionella species are described. Water systems contaminated with Legionella spp. are the implicated sources of Legionnaires' disease. In this study, we aimed to assess Legionella contamination in the water sources of a tertiary care hospital and to determine the virulence properties and molecular characteristics of L. pneumophila environmental isolates.

Methods: During May 2015 through August 2018, a total of 201 hospital water samples were tested for L. pneumophila by standardized culture procedures; environmental isolates were examined for the presence of two virulence genes: Legionella vir homolog (lvh) and repeats in structural toxin (rtxA) by PCR. The genotyping of isolates was performed by sequence-based typing (SBT) according to the protocol of the European Study Group for Legionella Infections (ESGLI).

Results:L. pneumophila was isolated from 38/201 (18.9%) water samples; among the 46 isolates, the lvh locus was present in 45 (97.8%), the rtxA locus was found in 45 (97.8%), and both loci were found in 44 (95.7%) isolates. A total of 23 sequence types (STs) were identified among the 44 isolates (index of discrimination [IOD] of 0.929), and 11/23 (47.8%) STs were new to the ESGLI database.

Conclusions: The study results showed genetic diversity in L. pneumophila isolates from the hospital environment along with a high percentage of pathogenicity loci. Besides, certain STs may have an increased ability to cause legionellosis, thus requires specific infection control and prevention strategies whenever identified.

Essential role of calcium in extending RTX adhesins to their target

•

Elongated beta-sandwich repeats are a major part of bacterial RTX adhesins.

•

The repeats are arranged in tandem to extend away from the bacterial surface.

•

Calcium ions are coordinated in the linkers between repeats to stiffen the protein.

•

Rigidification of the tandem repeats further helps extension of the adhesin.

•

The repeats differ greatly between species, but all have Ca²⁺ in their linkers.

RTX adhesins are long, multi-domain proteins present on the outer membrane of many Gram-negative bacteria. From this vantage point, adhesins use their distal ligand-binding domains for surface attachment leading to biofilm formation. To expand the reach of the ligand-binding domains, RTX adhesins maintain a central extender region of multiple tandem repeats, which makes up most of the proteins’ large molecular weight. Alignments of the 10-15-kDa extender domains show low sequence identity between adhesins. Here we have produced and structurally characterized protein constructs of four tandem repeats (tetra-tandemers) from two different RTX adhesins. In comparing the tetra-tandemers to each other and already solved structures from Marinomonas primoryensis and Salmonella enterica, the extender domains fold as diverse beta-sandwich structures with widely differing calcium contents. However, all the tetra-tandemers have at least one calcium ion coordinated in the linker region between beta-sandwich domains whose role appears to be the rigidification of the extender region to help the adhesin extend its reach.

From Input to Output: The Lap/c-di-GMP Biofilm Regulatory Circuit

Biofilms are the dominant bacterial lifestyle. The regulation of the formation and dispersal of bacterial biofilms has been the subject of study in many organisms. Over the last two decades, the mechanisms of Pseudomonas fluorescens biofilm formation and regulation have emerged as among the best understood of any bacterial biofilm system. Biofilm formation by P. fluorescens occurs through the localization of an adhesin, LapA, to the outer membrane via a variant of the classical type I secretion system. The decision between biofilm formation and dispersal is mediated by LapD, a c-di-GMP receptor, and LapG, a periplasmic protease, which together control whether LapA is retained or released from the cell surface. LapA localization is also controlled by a complex network of c-di-GMP-metabolizing enzymes. This review describes the current understanding of LapA-mediated biofilm formation by P. fluorescens and discusses several emerging models for the regulation and function of this adhesin.

Type 1 Does the Two-Step: Type 1 Secretion Substrates with a Functional Periplasmic Intermediate

Bacteria have evolved several secretion strategies for polling and responding to environmental flux and insult. Of these, the type 1 secretion system (T1SS) is known to secrete an array of biologically diverse proteins-from small, <10-kDa bacteriocins to gigantic adhesins with a mass >1 MDa. For the last several decades, T1SSs have been characterized as a one-step translocation strategy whereby the secreted substrate is transported directly into the extracellular environment from the cytoplasm with no periplasmic intermediate. Recent phylogenetic, biochemical, and genetic evidences point to a distinct subgroup of T1SS machinery linked with a bacterial transglutaminase-like cysteine proteinase (BTLCP), which uses a two-step secretion mechanism. BTLCP-linked T1SSs transport a class of repeats-in-toxin (RTX) adhesins that are critical for biofilm formation. The prototype of this RTX adhesin group, LapA of Pseudomonas fluorescens Pf0-1, uses a novel N-terminal retention module to anchor the adhesin at the cell surface as a secretion intermediate threaded through the outer membrane-localized TolC-like protein LapE. This secretion intermediate is posttranslationally cleaved by the BTLCP family LapG protein to release LapA from its cognate T1SS pore. Thus, the secretion of LapA and related RTX adhesins into the extracellular environment appears to be a T1SS-mediated two-step process that involves a periplasmic intermediate. In this review, we contrast the T1SS machinery and substrates of the BLTCP-linked two-step secretion process with those of the classical one-step T1SS to better understand the newly recognized and expanded role of this secretion machinery.

The Xylella fastidosa RTX operons: evidence for the evolution of protein mosaics through novel genetic exchanges

BACKGROUND

Xylella fastidiosa (Xf) is a gram negative bacterium inhabiting the plant vascular system. In most species this bacterium lives as a benign symbiote, but in several agriculturally important plants (e.g. coffee, citrus, grapevine) Xf is pathogenic. Xf has four loci encoding homologues to hemolysin RTX proteins, virulence factors involved in a wide range of plant pathogen interactions.

RESULTS

We show that all four genes are expressed during pathogenesis in grapevine. The sequences from these four genes have a complex repetitive structure. At the C-termini, sequence diversity between strains is what would be expected from orthologous genes. However, within strains there is no N-terminal homology, indicating these loci encode RTXs of different functions and/or specificities. More striking is that many of the orthologous loci between strains share this extreme variation at the N-termini. Thus these RTX orthologues are most easily visualized as fusions between the orthologous C-termini and different N-termini. Further, the four genes are found in operons having a peculiar structure with an extensively duplicated module encoding a small protein with homology to the N-terminal region of the full length RTX. Surprisingly, some of these small peptides are most similar not to their corresponding full length RTX, but to the N-termini of RTXs from other Xf strains, and even other remotely related species.

CONCLUSIONS

These results demonstrate that these genes are expressed in planta during pathogenesis. Their structure suggests extensive evolutionary restructuring through horizontal gene transfers and heterologous recombination mechanisms. The sum of the evidence suggests these repetitive modules are a novel kind of mobile genetic element.

Epidemiological Investigation of Legionella pneumophila Serogroup 2 to 14 Isolates from Water Samples by Amplified Fragment Length Polymorphism and Sequence-Based Typing and Detection of Virulence Traits

The aim of this study is to explore the dispersion, clonality, and virulence of Legionella pneumophila serogroups 2 to 14 in the Greek environment. Eighty L. pneumophila serogroup 2 to 14 strains isolated from water distribution systems of hotels, hospitals, athletic venues, and ferries in Greece were tested by monoclonal antibodies (MAbs) for serogroup discrimination and molecularly by amplified fragment length polymorphism (AFLP) for genetic diversity. Fifty-six of 80 strains were also typed by the sequence-based typing (SBT) method. Αll strains were further analyzed for detection of two pathogenicity loci: Legionella vir homologue (lvh) and repeats in structural toxin (rtxA). Thirty-seven strains (46.2%) belonged to serogroup 6, 26 strains (32.5%) to serogroup 3, and 7 (8.8%) to other serogroups (4, 5, 8, and 10). Ten strains (12.5%) were nontypeable (NT) into the known serogroups. Thirty-nine different AFLP types were found among the 80 L. pneumophila serogroup 2 to 14 strains, and 24 different SBT types were found among the 56 strains tested. Among the 80 strains, the lvh locus was present in 75 (93.8%), the rtxA locus was found in 76 (95%), and both loci were found in 73 (91.3%) strains. This study showed that there is genetic variability of L. pneumophila serogroups 2 to 14 in the Greek environment as well as a high percentage of the pathogenicity loci. Ιntroducing an effective diagnostic test for L. pneumophila serogroups 2 to 14 in urine and promoting the examination of respiratory specimens from patients hospitalized for pneumonia in Greek hospitals are essential.

IMPORTANCE

In this study, the dispersion, clonality, and virulence of environmental isolates of Legionella pneumophila serogroups 2 to 14 (Lp2-14) in Greece were investigated. Genetic variability of Lp2-14 in the Greek environment was identified together with the presence of the pathogenicity loci in a high percentage of the isolates. Despite the high prevalence of Lp2-14 in the Greek environment, no clinical cases were reported, which may be due to underdiagnosis of the disease. Almost all the legionellosis cases are diagnosed in Greece by using the urine antigen test, which is specific for Lp1. There is an urgent need to improve the clinical diagnosis of legionellosis by introducing an effective diagnostic test for Lp2-14 in urine and by promoting the PCR examination of respiratory specimens from patients with compatible clinical symptoms.

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.

Genomic Resolution of Outbreak-Associated Legionella pneumophila Serogroup 1 Isolates from New York State

A total of 30 Legionella pneumophila serogroup 1 isolates representing 10 separate legionellosis laboratory investigations ("outbreaks") that occurred in New York State between 2004 and 2012 were selected for evaluation of whole-genome sequencing (WGS) approaches for molecular subtyping of this organism. Clinical and environmental isolates were available for each outbreak and were initially examined by pulsed-field gel electrophoresis (PFGE). Sequence-based typing alleles were extracted from WGS data yielding complete sequence types (ST) for isolates representing 8 out of the 10 outbreaks evaluated in this study. Isolates from separate outbreaks sharing the same ST also contained the fewest differences in core genome single nucleotide polymorphisms (SNPs) and the greatest proportion of identical allele sequences in a whole-genome multilocus sequence typing (wgMLST) scheme. Both core SNP and wgMLST analyses distinguished isolates from separate outbreaks, including those from two outbreaks sharing indistinguishable PFGE profiles. Isolates from a hospital-associated outbreak spanning multiple years shared indistinguishable PFGE profiles but displayed differences in their genome sequences, suggesting the presence of multiple environmental sources. Finally, the rtx gene demonstrated differences in the repeat region sequence among ST1 isolates from different outbreaks, suggesting that variation in this gene may be useful for targeted molecular subtyping approaches for L. pneumophila This study demonstrates the utility of various genome sequence analysis approaches for L. pneumophila for environmental source attribution studies while furthering the understanding of Legionella ecology.

IMPORTANCE

We demonstrate that whole-genome sequencing helps to improve resolution of Legionella pneumophila isolated during laboratory investigations of legionellosis compared to traditional subtyping methods. These data can be important in confirming the environmental sources of legionellosis outbreaks. Moreover, we evaluated various methods to analyze genome sequence data to help resolve outbreak-related isolates.

Prediction of the origin of French Legionella pneumophila strains using a mixed-genome microarray

Background

Legionella is a water and soil bacterium that can infect humans, causing a pneumonia known as Legionnaires’ disease. The pneumonia is almost exclusively caused by the species L. pneumophila, of which serogroup 1 is responsible for 90% of patients. Within serogroup 1, large differences in prevalence in clinical isolates have been described. A recent study, using a Dutch Legionella strain collection, identified five virulence associated markers. In our study, we verify whether these five Dutch markers can predict the patient or environmental origin of a French Legionella strain collection. In addition, we identify new potential virulence markers and verify whether these can predict better. A total of 219 French patient isolates and environmental strains were compared using a mixed-genome micro-array. The micro-array data were analysed to identify predictive markers, using a Random Forest algorithm combined with a logistic regression model. The sequences of the identified markers were compared with eleven known Legionella genomes, using BlastN and BlastX; the functionality for each of the predictive markers was checked in the literature.

Results

The five Dutch markers insufficiently predicted the patient or environmental origin of the French Legionella strains. Subsequent analyses identified four predictive markers for the French collection that were used for the logistic regression model. This model showed a negative predictive value of 91%. Three of the French markers differed from the Dutch markers, one showed considerable overlap and was found in one of the Legionella genomes (Lorraine strain). This marker encodes for a structural toxin protein RtxA, described for L. pneumophila as a factor involved in virulence and entry in both human cells and amoebae.

Conclusions

The combination of a mixed-genome micro-array and statistical analysis using a Random Forest algorithm has identified virulence markers in a consistent way. The Lorraine strain and related Dutch and French Legionella strains contain a marker that encodes a RtxA protein which probably is involved in the increased prevalence in clinical isolates. The current set of predictive markers is insufficient to justify its use as a reliable test in the public health field in France. Our results suggest that genetic differences in Legionella strains exist between geographically distinct entities. It may be necessary to develop region-specific mixed-genome microarrays that are constantly adapted and updated.

Phylogenetic reconstruction of the Legionella pneumophila Philadelphia-1 laboratory strains through comparative genomics

Over 20 years ago, two groups independently domesticated Legionella pneumophila from a clinical isolate of bacteria collected during the first recognized outbreak of Legionnaires' disease (at the 1976 American Legion's convention in Philadelphia). These two laboratory strains, JR32 and Lp01, along with their derivatives, have been disseminated to a number of laboratories around the world and form the cornerstone of much of the research conducted on this important pathogen to date. Nevertheless, no exhaustive examination of the genetic distance between these strains and their clinical progenitor has been performed thus far. Such information is of paramount importance for making sense of several phenotypic differences observed between these strains. As environmental replication of L. pneumophila is thought to exclusively occur within natural protozoan hosts, retrospective analysis of the domestication and axenic culture of the Philadelphia-1 progenitor strain by two independent groups also provides an excellent opportunity to uncover evidence of adaptation to the laboratory environment. To reconstruct the phylogenetic relationships between the common laboratory strains of L. pneumophila Philadelphia-1 and their clinical ancestor, we performed whole-genome Illumina resequencing of the two founders of each laboratory lineage: JR32 and Lp01. As expected from earlier, targeted studies, Lp01 and JR32 contain large deletions in the lvh and tra regions, respectively. By sequencing additional strains derived from Lp01 (Lp02 and Lp03), we retraced the phylogeny of these strains relative to their reported ancestor, thereby reconstructing the evolutionary dynamics of each laboratory lineage from genomic data.

Evolution at two levels in fire ants: the relationship between patterns of gene expression and protein sequence evolution

Variation in protein sequence and gene expression each contribute to phenotypic diversity, and may be subject to similar selective pressures. Eusocial insects are particularly useful for investigating the evolutionary link between protein sequence and condition-dependent patterns of gene expression because gene expression plays a central role in determining differences between eusocial insect sexes and castes. We investigated the relationship between protein coding sequence evolution and gene expression patterns in the fire ants Solenopsis invicta, S. richteri, and their hybrids to gain greater insight into how selection jointly operates on gene expression and coding sequence. We found that genes with high expression variability within castes and sexes were frequently differentially expressed between castes and sexes, as well as between species and hybrids. These results indicate that genes showing high variation in expression in one context also tend to show high variation in expression in other contexts. Our analyses further revealed that variation in both intra- and interspecific gene expression was positively associated with rate of protein sequence evolution in Solenopsis. This suggests that selective constraints on a gene operate both at the level of protein sequence and at the level of gene expression regulation. Overall, our study provides one of the strongest demonstrations that selective constraints mediate both protein sequence evolution and gene expression variability across different biological contexts and timescales.

Predictive sequence analysis of the Candidatus Liberibacter asiaticus proteome

Candidatus Liberibacter asiaticus (Ca. L. asiaticus) is a parasitic Gram-negative bacterium that is closely associated with Huanglongbing (HLB), a worldwide citrus disease. Given the difficulty in culturing the bacterium and thus in its experimental characterization, computational analyses of the whole Ca. L. asiaticus proteome can provide much needed insights into the mechanisms of the disease and guide the development of treatment strategies. In this study, we applied state-of-the-art sequence analysis tools to every Ca. L. asiaticus protein. Our results are available as a public website at http://prodata.swmed.edu/liberibacter_asiaticus/. In particular, we manually curated the results to predict the subcellular localization, spatial structure and function of all Ca. L. asiaticus proteins (http://prodata.swmed.edu/liberibacter_asiaticus/curated/). This extensive information should facilitate the study of Ca. L. asiaticus proteome function and its relationship to disease. Pilot studies based on the information from our website have revealed several potential virulence factors, discussed herein.

Quantifying protein modularity and evolvability: a comparison of different techniques

Modularity increases evolvability by reducing constraints on adaptation and by allowing preexisting parts to function in new contexts for novel uses. Protein evolution provides an excellent context to study the causes and consequences of biological modularity. In order to address such questions, however, an index for protein modularity is necessary. This paper proposes a simple index for protein modularity-"module density"-which is the number of evolutionarily independent modules that compose a protein divided by the number of amino acids in the protein. The decomposition of proteins into constituent modules can be accomplished by either of two classes of methods. The first class of methods relies on "suppositional" criteria to assign amino acids to modules, whereas the second class of methods relies on "coevolutionary" criteria for this task. One simple and practical method from the first class consists of approximating the number of modules in a protein as the number of regular secondary structure elements (i.e., helices and sheets). Methods based on coevolutionary criteria require more elaborate data, but they have the advantage of being able to specify modules without prior assumptions about why they exist. Given the increasing availability of datasets sampling protein mutational spectra (e.g., from comparative genomics, experimental evolution, and computational prediction), methods based on coevolutionary criteria will likely become more promising in the near future. The ability to meaningfully quantify protein modularity via simple indices has the potential to aid future efforts to understand protein evolutionary rate determinants, improve molecular evolution models and engineer novel proteins.

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.

Diversity and impact of prokaryotic toxins on aquatic environments: a review

Microorganisms are ubiquitous in all habitats and are recognized by their metabolic versatility and ability to produce many bioactive compounds, including toxins. Some of the most common toxins present in water are produced by several cyanobacterial species. As a result, their blooms create major threats to animal and human health, tourism, recreation and aquaculture. Quite a few cyanobacterial toxins have been described, including hepatotoxins, neurotoxins, cytotoxins and dermatotoxins. These toxins are secondary metabolites, presenting a vast diversity of structures and variants. Most of cyanobacterial secondary metabolites are peptides or have peptidic substructures and are assumed to be synthesized by non-ribosomal peptide synthesis (NRPS), involving peptide synthetases, or NRPS/PKS, involving peptide synthetases and polyketide synthases hybrid pathways. Besides cyanobacteria, other bacteria associated with aquatic environments are recognized as significant toxin producers, representing important issues in food safety, public health, and human and animal well being. Vibrio species are one of the most representative groups of aquatic toxin producers, commonly associated with seafood-born infections. Some enterotoxins and hemolysins have been identified as fundamental for V. cholerae and V. vulnificus pathogenesis, but there is evidence for the existence of other potential toxins. Campylobacter spp. and Escherichia coli are also water contaminants and are able to produce important toxins after infecting their hosts. Other bacteria associated with aquatic environments are emerging as toxin producers, namely Legionella pneumophila and Aeromonas hydrophila, described as responsible for the synthesis of several exotoxins, enterotoxins and cytotoxins. Furthermore, several Clostridium species can produce potent neurotoxins. Although not considered aquatic microorganisms, they are ubiquitous in the environment and can easily contaminate drinking and irrigation water. Clostridium members are also spore-forming bacteria and can persist in hostile environmental conditions for long periods of time, contributing to their hazard grade. Similarly, Pseudomonas species are widespread in the environment. Since P. aeruginosa is an emergent opportunistic pathogen, its toxins may represent new hazards for humans and animals. This review presents an overview of the diversity of toxins produced by prokaryotic microorganisms associated with aquatic habitats and their impact on environment, life and health of humans and other animals. Moreover, important issues like the availability of these toxins in the environment, contamination sources and pathways, genes involved in their biosynthesis and molecular mechanisms of some representative toxins are also discussed.

Legionella pneumophila pangenome reveals strain-specific virulence factors

BACKGROUND

Legionella pneumophila subsp. pneumophila is a gram-negative gamma-Proteobacterium and the causative agent of Legionnaires' disease, a form of epidemic pneumonia. It has a water-related life cycle. In industrialized cities L. pneumophila is commonly encountered in refrigeration towers and water pipes. Infection is always via infected aerosols to humans. Although many efforts have been made to eradicate Legionella from buildings, it still contaminates the water systems. The town of Alcoy (Valencian Region, Spain) has had recurrent outbreaks since 1999. The strain "Alcoy 2300/99" is a particularly persistent and recurrent strain that was isolated during one of the most significant outbreaks between the years 1999-2000.

RESULTS

We have sequenced the genome of the particularly persistent L. pneumophila strain Alcoy 2300/99 and have compared it with four previously sequenced strains known as Philadelphia (USA), Lens (France), Paris (France) and Corby (England).Pangenome analysis facilitated the identification of strain-specific features, as well as some that are shared by two or more strains. We identified: (1) three islands related to anti-drug resistance systems; (2) a system for transport and secretion of heavy metals; (3) three systems related to DNA transfer; (4) two CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats) systems, known to provide resistance against phage infections, one similar in the Lens and Alcoy strains, and another specific to the Paris strain; and (5) seven islands of phage-related proteins, five of which seem to be strain-specific and two shared.

CONCLUSIONS

The dispensable genome disclosed by the pangenomic analysis seems to be a reservoir of new traits that have mainly been acquired by horizontal gene transfer and could confer evolutionary advantages over strains lacking them.

Intragenic tandem repeat variation between Legionella pneumophila strains

Background

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

Results

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

Conclusion

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