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

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.

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.

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.