Intragenic tandem repeat variation between Legionella pneumophila strains

The Legionella collagen-like protein employs a distinct binding mechanism for the recognition of host glycosaminoglycans

Bacterial adhesion is a fundamental process which enables colonisation of niche environments and is key for infection. However, in Legionella pneumophila, the causative agent of Legionnaires' disease, these processes are not well understood. The Legionella collagen-like protein (Lcl) is an extracellular peripheral membrane protein that recognises sulphated glycosaminoglycans on the surface of eukaryotic cells, but also stimulates bacterial aggregation in response to divalent cations. Here we report the crystal structure of the Lcl C-terminal domain (Lcl-CTD) and present a model for intact Lcl. Our data reveal that Lcl-CTD forms an unusual trimer arrangement with a positively charged external surface and negatively charged solvent exposed internal cavity. Through molecular dynamics simulations, we show how the glycosaminoglycan chondroitin-4-sulphate associates with the Lcl-CTD surface via distinct binding modes. Our findings show that Lcl homologs are present across both the Pseudomonadota and Fibrobacterota-Chlorobiota-Bacteroidota phyla and suggest that Lcl may represent a versatile carbohydrate-binding mechanism.

The Legionella collagen-like protein employs a unique binding mechanism for the recognition of host glycosaminoglycans

Bacterial adhesion is a fundamental process which enables colonisation of niche environments and is key for infection. However, in Legionella pneumophila, the causative agent of Legionnaires' disease, these processes are not well understood. The Legionella collagen-like protein (Lcl) is an extracellular peripheral membrane protein that recognises sulphated glycosaminoglycans (GAGs) on the surface of eukaryotic cells, but also stimulates bacterial aggregation in response to divalent cations. Here we report the crystal structure of the Lcl C-terminal domain (Lcl-CTD) and present a model for intact Lcl. Our data reveal that Lcl-CTD forms an unusual dynamic trimer arrangement with a positively charged external surface and a negatively charged solvent exposed internal cavity. Through Molecular Dynamics (MD) simulations, we show how the GAG chondroitin-4-sulphate associates with the Lcl-CTD surface via unique binding modes. Our findings show that Lcl homologs are present across both the Pseudomonadota and Fibrobacterota-Chlorobiota-Bacteroidota phyla and suggest that Lcl may represent a versatile carbohydrate binding mechanism.

Genetic markers of Chlamydia pecorum virulence in ruminants support short term host-pathogen evolutionary relationships in the koala, Phascolarctos cinereus

In ruminants infected with Chlamydia pecorum, shorter lengths of coding tandem repeats (CTR) within two genes, the inclusion membrane protein (incA) and Type III secretor protein (ORF663), have been previously associated with pathogenic outcomes. In other chlamydial species, the presence of a chlamydial plasmid has been linked to heightened virulence, and the plasmid is not ubiquitous in C. pecorum across the koala's range. We therefore investigated these three markers: incA, ORF663 and C. pecorum plasmid, as potential indicators of virulence in two koala populations in New South Wales with differing expression of urogenital chlamydiosis; the Liverpool Plains and one across the Southern Highlands and South-west Sydney (SHSWS). We also investigated the diversity of these loci within strains characterised by the national multi-locus sequence typing (MLST) scheme. Although CTR lengths of incA and ORF663 varied across the populations, they occurred only within previously described pathogenic ranges for ruminants. This suggests a relatively short-term host-pathogen co-evolution within koalas and limits the utility of CTR lengths for incA and ORF663 as virulence markers in the species. However, in contrast to reports of evolution of C. pecorum towards lower virulence, as indicated by longer CTR lengths in ruminants and swine, CTR lengths for ORF663 appeared to be diverging towards less common shorter CTR lengths within strains recently introduced to koalas in the Liverpool Plains. We detected the plasmid across 90% and 92% of samples in the Liverpool Plains and SHSWS respectively, limiting its utility as an indicator of virulence. It would be valuable to examine the CTR lengths of these loci across koala populations nationally. Investigation of other hypervariable loci may elucidate the evolutionary trajectory of virulence in C. pecorum induced disease in koalas. Profiling of virulent strains will be important in risk assessments for strain movement to naïve or susceptible populations through translocations and wildlife corridor construction.

The evolution and role of eukaryotic-like domains in environmental intracellular bacteria: the battle with a eukaryotic cell

Intracellular pathogens that are able to thrive in different environments, such as Legionella spp. which preferentially live in protozoa in aquatic environments or environmental Chlamydiae which replicate either within protozoa or a range of animals, possess a plethora of cellular biology tools to influence their eukaryotic host. The host manipulation tools that evolved in the interaction with protozoa, confer these bacteria the capacity to also infect phylogenetically distinct eukaryotic cells, such as macrophages and thus they can also be human pathogens. To manipulate the host cell, bacteria use protein secretion systems and molecular effectors. Although these molecular effectors are encoded in bacteria, they are expressed and function in a eukaryotic context often mimicking or inhibiting eukaryotic proteins. Indeed, many of these effectors have eukaryotic-like domains. In this review we propose that the main pathways environmental intracellular bacteria need to subvert in order to establish the host eukaryotic cell as a replication niche are chromatin remodelling, ubiquitination signalling, and modulation of protein-protein interactions via tandem repeat domains. We then provide mechanistic insight into how these proteins might have evolved as molecular weapons. Finally, we highlight that in environmental intracellular bacteria the number of eukaryotic-like domains and proteins is considerably higher than in intracellular bacteria specialised to an isolated niche, such as obligate intracellular human pathogens. As mimics of eukaryotic proteins are critical components of host pathogen interactions, this distribution of eukaryotic-like domains suggests that the environment has selected them.

Characterisation of the novel clinical isolate X-4 containing a new tp0548 sequence-type

OBJECTIVES

A novel tp0548 sequence-type was identified in one clinical isolate (X-4) from a patient diagnosed with primary syphilis in Xiamen, China. To precisely define and characterise a new clinical isolate, we performed further genome-scale molecular analysis.

METHODS

The pooled segment genome sequencing method followed by Illumina sequencing was performed.

RESULTS

This novel sequence-type contained a unique nucleotide substitution 'T' at position 167 and belonged to the SS14-like clade of TPA strains, as determined by phylogenetic analysis. Multi-locus sequence analysis of nine chromosomal loci demonstrated that the X-4 isolate was clustered within a monophyletic group of TPA strains. Whole-genome phylogenetic analysis subsequently corroborated the TPA strain classification of the X-4 isolate and revealed that the isolate was closely related to the SS14 strain, with 42 single-nucleotide variations and 12 insertions/deletions. In addition, high intrastrain heterogeneity in the length of the poly G/C tracts was found in the TPAChi_0347 locus, which might indicate that this gene of the X-4 isolate is likely involved in phase variation events. The length heterogeneity of the poly A/T tracts was lower than the genetic variability of the poly G/C tracts, and all the observed intrastrain variations fell within coding regions.

CONCLUSION

The novel tp0548 sequence-type was determined to belong to a new TPA isolate, X-4. The identification of variable length in homopolymetic tracts (G/C and A/T) could provide a snapshot of the genes that potentially involved in genotype-phenotype variations. These findings provide an unequivocal characterisation for better understanding the molecular variation of this emerging isolate.

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.

Cyclic AMP-Independent Control of Twitching Motility in Pseudomonas aeruginosa

FimV is a Pseudomonas aeruginosa inner membrane hub protein that modulates levels of the second messenger, cyclic AMP (cAMP), through the activation of adenylate cyclase CyaB. Although type IVa pilus (T4aP)-dependent twitching motility is modulated by cAMP levels, mutants lacking FimV are twitching impaired, even when exogenous cAMP is provided. Here we further define FimV's cAMP-dependent and -independent regulation of twitching. We confirmed that the response regulator of the T4aP-associated Chp chemotaxis system, PilG, requires both FimV and the CyaB regulator, FimL, to activate CyaB. However, in cAMP-replete backgrounds-lacking the cAMP phosphodiesterase CpdA or the CheY-like protein PilH or expressing constitutively active CyaB-pilG and fimV mutants failed to twitch. Both cytoplasmic and periplasmic domains of FimV were important for its cAMP-dependent and -independent roles, while its septal peptidoglycan-targeting LysM motif was required only for twitching motility. Polar localization of the sensor kinase PilS, a key regulator of transcription of the major pilin, was FimV dependent. However, unlike its homologues in other species that localize flagellar system components, FimV was not required for swimming motility. These data provide further evidence to support FimV's role as a key hub protein that coordinates the polar localization and function of multiple structural and regulatory proteins involved in P. aeruginosa twitching motility.IMPORTANCEPseudomonas aeruginosa is a serious opportunistic pathogen. Type IVa pili (T4aP) are important for its virulence, because they mediate dissemination and invasion via twitching motility and are involved in surface sensing, which modulates pathogenicity via changes in cAMP levels. Here we show that the hub protein FimV and the response regulator of the Chp system, PilG, regulate twitching independently of their roles in the modulation of cAMP synthesis. These functions do not require the putative scaffold protein FimL, proposed to link PilG with FimV. PilG may regulate asymmetric functioning of the T4aP system to allow for directional movement, while FimV appears to localize both structural and regulatory elements-including the PilSR two-component system-to cell poles for optimal function.

Variable number of tandem repeats of 9 Plasmodium vivax genes among Southeast Asian isolates

The variable number of tandem repeats (VNTRs) provides valuable information about both the functional and evolutionary aspects of genetic diversity. Comparative analysis of 3 Plasmodium falciparum genomes has shown that more than 9% of its open reading frames (ORFs) harbor VNTRs. Although microsatellites and VNTR genes of P. vivax were reported, the VNTR polymorphism of genes has not been examined widely. In this study, 230 P. vivax genes were analyzed for VNTRs by SERV, and 33 kinds of TR deletions or insertions from 29 P. vivax genes (12.6%) were found. Of these, 9 VNTR fragments from 8 P. vivax genes were used for PCR amplification and sequence analysis to examine the genetic diversity among 134 isolates from four Southeast Asian countries (China, Republic of Korea, Thailand, and Myanmar) with different malaria endemicity. We confirmed the existence of extensive polymorphism of VNTR fragments in field isolates. This detection provides several suitable markers for analysis of the molecular epidemiology of P. vivax field isolates.

Temperature-Dependent Growth Modeling of Environmental and Clinical Legionella pneumophila Multilocus Variable-Number Tandem-Repeat Analysis (MLVA) Genotypes

Legionella pneumophila causes waterborne infections resulting in severe pneumonia. High-resolution genotyping of L. pneumophila isolates can be achieved by multiple-locus variable-number tandem-repeat analysis (MLVA). Recently, we found that different MLVA genotypes of L. pneumophila dominated different sites in a small drinking-water network, with a genotype-related temperature and abundance regime. The present study focuses on understanding the temperature-dependent growth kinetics of the genotypes that dominated the water network. Our aim was to model mathematically the influence of temperature on the growth kinetics of different environmental and clinical L. pneumophila genotypes and to compare it with the influence of their ecological niches. Environmental strains showed a distinct temperature preference, with significant differences among the growth kinetics of the three studied genotypes (Gt4, Gt6, and Gt15). Gt4 strains exhibited superior growth at lower temperatures (25 and 30°C), while Gt15 strains appeared to be best adapted to relatively higher temperatures (42 and 45°C). The temperature-dependent growth traits of the environmental genotypes were consistent with their distribution and temperature preferences in the water network. Clinical isolates exhibited significantly higher growth rates and reached higher maximal cell densities at 37°C and 42°C than the environmental strains. Further research on the growth preferences of L. pneumophila clinical and environmental genotypes will result in a better understanding of their ecological niches in drinking-water systems as well as in the human body.IMPORTANCELegionella pneumophila is a waterborne pathogen that threatens humans in developed countries. The bacteria inhabit natural and man-made freshwater environments. Here we demonstrate that different environmental L. pneumophila genotypes have different temperature-dependent growth kinetics. Moreover, Legionella strains that belong to the same species but were isolated from environmental and clinical sources possess adaptations for growth at different temperatures. These growth preferences may influence the bacterial colonization at specific ecological niches within the drinking-water network. Adaptations for growth at human body temperatures may facilitate the abilities of some L. pneumophila strains to infect and cause illness in humans. Our findings may be used as a tool to improve Legionella monitoring in drinking-water networks. Risk assessment models for predicting the risk of legionellosis should take into account not only Legionella concentrations but also the temperature-dependent growth kinetics of the isolates.

The Conserved Tetratricopeptide Repeat-Containing C-Terminal Domain of Pseudomonas aeruginosa FimV Is Required for Its Cyclic AMP-Dependent and -Independent Functions

FimV is a Pseudomonas aeruginosa inner membrane protein that regulates intracellular cyclic AMP (cAMP) levels-and thus type IV pilus (T4P)-mediated twitching motility and type II secretion (T2S)-by activating the adenylate cyclase CyaB. Its cytoplasmic domain contains three predicted tetratricopeptide repeat (TPR) motifs separated by an unstructured region: two proximal to the inner membrane and one within the "FimV C-terminal domain," which is highly conserved across diverse homologs. Here, we present the crystal structure of the FimV C terminus, FimV861-919, containing a TPR motif decorated with solvent-exposed, charged side chains, plus a C-terminal capping helix. FimV689, a truncated form lacking this C-terminal motif, did not restore wild-type levels of twitching or surface piliation compared to the full-length protein. FimV689 failed to restore wild-type levels of the T4P motor ATPase PilU or T2S, suggesting that it was unable to activate cAMP synthesis. Bacterial two-hybrid analysis showed that TPR3 interacts directly with the CyaB activator, FimL. However, FimV689 failed to restore wild-type motility in a fimV mutant expressing a constitutively active CyaB (fimV cyaB-R456L), suggesting that the C-terminal motif is also involved in cAMP-independent functions of FimV. The data show that the highly conserved TPR-containing C-terminal domain of FimV is critical for its cAMP-dependent and -independent functions.

IMPORTANCE

FimV is important for twitching motility and cAMP-dependent virulence gene expression in P. aeruginosa FimV homologs have been identified in several human pathogens, and their functions are not limited to T4P expression. The C terminus of FimV is remarkably conserved among otherwise very diverse family members, but its role is unknown. We provide here biological evidence for the importance of the C-terminal domain in both cAMP-dependent (through FimL) and -independent functions of FimV. We present X-ray crystal structures of the conserved C-terminal domain and identify a consensus sequence for the C-terminal TPR within the conserved domain. Our data extend our knowledge of FimV's functionally important domains, and the structures and consensus sequences provide a foundation for studies of FimV and its homologs.

Comparative genome analysis of Weissella ceti, an emerging pathogen of farm-raised rainbow trout

Background

The genus Weissella belongs to the lactic acid bacteria and includes 18 currently identified species, predominantly isolated from fermented food but rarely from cases of bacteremia in animals. Recently, a new species, designated Weissella ceti, has been correlated with hemorrhagic illness in farm-raised rainbow trout in China, Brazil, and the USA, with high transmission and mortality rates during outbreaks. Although W. ceti is an important emerging veterinary pathogen, little is known about its genomic features or virulence mechanisms. To better understand these and to characterize the species, we have previously sequenced the genomes of W. ceti strains WS08, WS74, and WS105, isolated from different rainbow trout farms in Brazil and displaying different pulsed-field gel electrophoresis patterns. Here, we present a comparative analysis of the three previously sequenced genomes of W. ceti strains from Brazil along with W. ceti NC36 from the USA and those of other Weissella species.

Results

Phylogenomic and orthology-based analyses both showed a high-similarity in the genetic structure of these W. ceti strains. This structure is corroborated by the highly syntenic order of their genes and the neutral evolution inferred from Tajima’s D. A whole-genome multilocus sequence typing analysis distinguished strains WS08 and NC36 from strains WS74 and WS105. We predicted 10 putative genomic islands (GEI), among which PAIs 3a and 3b are phage sequences that occur only in WS105 and WS74, respectively, whereas PAI 1 is species specific.

Conclusions

We identified several genes putatively involved in the basic processes of bacterial physiology and pathogenesis, including survival in aquatic environment, adherence in the host, spread inside the host, resistance to immune-system-mediated stresses, and antibiotic resistance. These data provide new insights in the molecular epidemiology and host adaptation for this emerging pathogen in aquaculture.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2324-4) contains supplementary material, which is available to authorized users.

Hybrid de novo tandem repeat detection using short and long reads

Background

As one of the most studied genome rearrangements, tandem repeats have a considerable impact on genetic backgrounds of inherited diseases. Many methods designed for tandem repeat detection on reference sequences obtain high quality results. However, in the case of a de novo context, where no reference sequence is available, tandem repeat detection remains a difficult problem. The short reads obtained with the second-generation sequencing methods are not long enough to span regions that contain long repeats. This length limitation was tackled by the long reads obtained with the third-generation sequencing platforms such as Pacific Biosciences technologies. Nevertheless, the gain on the read length came with a significant increase of the error rate. The main objective of nowadays studies on long reads is to handle the high error rate up to 16%.

Methods

In this paper we present MixTaR, the first de novo method for tandem repeat detection that combines the high-quality of short reads and the large length of long reads. Our hybrid algorithm uses the set of short reads for tandem repeat pattern detection based on a de Bruijn graph. These patterns are then validated using the long reads, and the tandem repeat sequences are constructed using local greedy assemblies.

Results

MixTaR is tested with both simulated and real reads from complex organisms. For a complete analysis of its robustness to errors, we use short and long reads with different error rates. The results are then analysed in terms of number of tandem repeats detected and the length of their patterns.

Conclusions

Our method shows high precision and sensitivity. With low false positive rates even for highly erroneous reads, MixTaR is able to detect accurate tandem repeats with pattern lengths varying within a significant interval.

The discovery, function and development of the variable number tandem repeats in different Mycobacterium species

The method of genotyping by variable number tandem repeats (VNTRs) facilitates the epidemiological studies of different Mycobacterium species worldwide. Until now, the VNTR method is not fully understood, for example, its discovery, function and classification. The inconsistent nomenclature and terminology of VNTR is especially confusing. In this review, we first describe in detail the VNTRs in Mycobacterium tuberculosis (M. tuberculosis), as this pathogen resulted in more deaths than any other microbial pathogen as well as for which extensive studies of VNTRs were carried out, and then we outline the recent progress of the VNTR-related epidemiological research in several other Mycobacterium species, such as M. abscessus, M. africanum, M. avium, M. bovis, M. canettii, M. caprae, M. intracellulare, M. leprae, M. marinum, M. microti, M. pinnipedii and M. ulcerans from different countries and regions. This article is aimed mainly at the practical notes of VNTR to help the scientists in better understanding and performing this method.

Host adaptation of Chlamydia pecorum towards low virulence evident in co-evolution of the ompA, incA, and ORF663 Loci

Chlamydia (C.) pecorum, an obligate intracellular bacterium, may cause severe diseases in ruminants, swine and koalas, although asymptomatic infections are the norm. Recently, we identified genetic polymorphisms in the ompA, incA and ORF663 genes that potentially differentiate between high-virulence C. pecorum isolates from diseased animals and low-virulence isolates from asymptomatic animals. Here, we expand these findings by including additional ruminant, swine, and koala strains. Coding tandem repeats (CTRs) at the incA locus encoded a variable number of repeats of APA or AGA amino acid motifs. Addition of any non-APA/AGA repeat motif, such as APEVPA, APAVPA, APE, or APAPE, associated with low virulence (P<10⁻⁴), as did a high number of amino acids in all incA CTRs (P = 0.0028). In ORF663, high numbers of 15-mer CTRs correlated with low virulence (P = 0.0001). Correction for ompA phylogram position in ORF663 and incA abolished the correlation between genetic changes and virulence, demonstrating co-evolution of ompA, incA, and ORF663 towards low virulence. Pairwise divergence of ompA, incA, and ORF663 among isolates from healthy animals was significantly higher than among strains isolated from diseased animals (P≤10⁻⁵), confirming the longer evolutionary path traversed by low-virulence strains. All three markers combined identified 43 unique strains and 4 pairs of identical strains among all 57 isolates tested, demonstrating the suitability of these markers for epidemiological investigations.

Deep conservation of human protein tandem repeats within the eukaryotes

Tandem repeats (TRs) are a major element of protein sequences in all domains of life. They are particularly abundant in mammals, where by conservative estimates one in three proteins contain a TR. High generation-scale duplication and deletion rates were reported for nucleic TR units. However, it is not known whether protein TR units can also be frequently lost or gained providing a source of variation for rapid adaptation of protein function, or alternatively, tend to have conserved TR unit configurations over long evolutionary times. To obtain a systematic picture, we performed a proteome-wide analysis of the mode of evolution for human protein TRs. For this purpose, we propose a novel method for the detection of orthologous TRs based on circular profile hidden Markov models. For all detected TRs, we reconstructed bispecies TR unit phylogenies across 61 eukaryotes ranging from human to yeast. Moreover, we performed additional analyses to correlate functional and structural annotations of human TRs with their mode of evolution. Surprisingly, we find that the vast majority of human TRs are ancient, with TR unit number and order preserved intact since distant speciation events. For example, ≥61% of all human TRs have been strongly conserved at least since the root of all mammals, approximately 300 Ma. Further, we find no human protein TR that shows evidence for strong recent duplications and deletions. The results are in contrast to the high generation-scale mutability of nucleic TRs. Presumably, most protein TRs fold into stable and conserved structures that are indispensable for the function of the TR-containing protein. All of our data and results are available for download from http://www.atgc-montpellier.fr/TRE.

Validation of 'variable number of tandem repeat'-based approach for examination of 'Candidatus Liberibacter asiaticus' diversity and its applications for the analysis of the pathogen populations in the areas of recent introduction

Citrus greening (Huanglongbing, HLB) is one of the most destructive diseases of citrus worldwide. In South Asia HLB has been known for more than a century, while in Americas the disease was found relatively recently. HLB is associated with three species of 'Candidatus Liberibacter' among which 'Ca. Liberibacter asiaticus' (CLas) has most wide distribution. Recently, a number of studies identified different regions in the CLas genome with variable number of tandem repeats (VNTRs) that could be used for examination of CLas diversity. One of the objectives of the work presented here was to further validate the VNTR analysis-based approach by assessing the stability of these repeats upon multiplication of the pathogen in a host over an extended period of time and upon its passaging from a host to a host using CLas populations from Florida. Our results showed that the numbers of tandem repeats in the four loci tested display very distinguishable "signature profiles" for the two Florida-type CLas haplotype groups. Remarkably, the profiles do not change upon passage of the pathogen in citrus and psyllid hosts as well as after its presence within a host over a period of five years, suggesting that VNTR analysis-based approach represents a valid methodology for examination of the pathogen populations in various geographical regions. Interestingly, an extended analysis of CLas populations in different locations throughout Florida and in several countries in the Caribbean and Central America regions and in Mexico where the pathogen has been introduced recently demonstrated the dispersion of the same haplotypes of CLas. On the other hand, these CLas populations appeared to differ significantly from those obtained from locations where the disease has been present for a much longer time.

A pilot study of rapid whole-genome sequencing for the investigation of a Legionella outbreak

OBJECTIVES

Epidemiological investigations of Legionnaires' disease outbreaks rely on the rapid identification and typing of clinical and environmental Legionella isolates in order to identify and control the source of infection. Rapid bacterial whole-genome sequencing (WGS) is an emerging technology that has the potential to rapidly discriminate outbreak from non-outbreak isolates in a clinically relevant time frame.

METHODS

We performed a pilot study to determine the feasibility of using bacterial WGS to differentiate outbreak from non-outbreak isolates collected during an outbreak of Legionnaires' disease. Seven Legionella isolates (three clinical and four environmental) were obtained from the reference laboratory and sequenced using the Illumina MiSeq platform at Addenbrooke's Hospital, Cambridge. Bioinformatic analysis was performed blinded to the epidemiological data at the Wellcome Trust Sanger Institute.

RESULTS

We were able to distinguish outbreak from non-outbreak isolates using bacterial WGS, and to confirm the probable environmental source. Our analysis also highlighted constraints, which were the small number of Legionella pneumophila isolates available for sequencing, and the limited number of published genomes for comparison.

CONCLUSIONS

We have demonstrated the feasibility of using rapid WGS to investigate an outbreak of Legionnaires' disease. Future work includes building larger genomic databases of L pneumophila from both clinical and environmental sources, developing automated data interpretation software, and conducting a cost-benefit analysis of WGS versus current typing methods.

Whole genome sequences of three Treponema pallidum ssp. pertenue strains: yaws and syphilis treponemes differ in less than 0.2% of the genome sequence

Background

The yaws treponemes, Treponema pallidum ssp. pertenue (TPE) strains, are closely related to syphilis causing strains of Treponema pallidum ssp. pallidum (TPA). Both yaws and syphilis are distinguished on the basis of epidemiological characteristics, clinical symptoms, and several genetic signatures of the corresponding causative agents.

Methodology/Principal Findings

To precisely define genetic differences between TPA and TPE, high-quality whole genome sequences of three TPE strains (Samoa D, CDC-2, Gauthier) were determined using next-generation sequencing techniques. TPE genome sequences were compared to four genomes of TPA strains (Nichols, DAL-1, SS14, Chicago). The genome structure was identical in all three TPE strains with similar length ranging between 1,139,330 bp and 1,139,744 bp. No major genome rearrangements were found when compared to the four TPA genomes. The whole genome nucleotide divergence (d_A) between TPA and TPE subspecies was 4.7 and 4.8 times higher than the observed nucleotide diversity (π) among TPA and TPE strains, respectively, corresponding to 99.8% identity between TPA and TPE genomes. A set of 97 (9.9%) TPE genes encoded proteins containing two or more amino acid replacements or other major sequence changes. The TPE divergent genes were mostly from the group encoding potential virulence factors and genes encoding proteins with unknown function.

Conclusions/Significance

Hypothetical genes, with genetic differences, consistently found between TPE and TPA strains are candidates for syphilitic treponemes virulence factors. Seventeen TPE genes were predicted under positive selection, and eleven of them coded either for predicted exported proteins or membrane proteins suggesting their possible association with the cell surface. Sequence changes between TPE and TPA strains and changes specific to individual strains represent suitable targets for subspecies- and strain-specific molecular diagnostics.

Spirochete Treponema pallidum ssp. pertenue (TPE) is the causative agent of yaws while strains of Treponema pallidum ssp. pallidum (TPA) cause syphilis. Both yaws and syphilis are distinguished on the basis of epidemiological characteristics and clinical symptoms. Neither treponeme can reproduce outside the host organism, which precludes the use of standard molecular biology techniques used to study cultivable pathogens. In this study, we determined high quality whole genome sequences of TPE strains and compared them to known genetic information for T. pallidum ssp. pallidum strains. The genome structure was identical in all three TPE strains and also between TPA and TPE strains. The TPE genome length ranged between 1,139,330 bp and 1,139,744 bp. The overall sequence identity between TPA and TPE genomes was 99.8%, indicating that the two pathogens are extremely closely related. A set of 34 TPE genes (3.5%) encoded proteins containing six or more amino acid replacements or other major sequence changes. These genes more often belonged to the group of genes with predicted virulence and unknown functions suggesting their involvement in infection differences between yaws and syphilis.

Investigation of the population structure of Legionella pneumophila by analysis of tandem repeat copy number and internal sequence variation

The population structure of the species Legionella pneumophila was investigated by multilocus variable number of tandem repeats (VNTR) analysis (MLVA) and sequencing of three VNTRs (Lpms01, Lpms04 and Lpms13) in selected strains. Of 150 isolates of diverse origins, 136 (86 %) were distributed into eight large MLVA clonal complexes (VACCs) and the rest were either unique or formed small clusters of up to two MLVA genotypes. In spite of the lower degree of genome-wide linkage disequilibrium of the MLVA loci compared with sequence-based typing, the clustering achieved by the two methods was highly congruent. The detailed analysis of VNTR Lpms04 alleles showed a very complex organization, with five different repeat unit lengths and a high level of internal variation. Within each MLVA-defined VACC, Lpms04 was endowed with a common recognizable pattern with some interesting exceptions. Evidence of recombination events was suggested by analysis of internal repeat variations at the two additional VNTR loci, Lpms01 and Lpms13. Sequence analysis of L. pneumophila VNTR locus Lpms04 alone provides a first-line assay for allocation of a new isolate within the L. pneumophila population structure and for epidemiological studies.

Genome analysis of Treponema pallidum subsp. pallidum and subsp. pertenue strains: most of the genetic differences are localized in six regions

The genomes of eight treponemes including T. p. pallidum strains (Nichols, SS14, DAL-1 and Mexico A), T. p. pertenue strains (Samoa D, CDC-2 and Gauthier), and the Fribourg-Blanc isolate, were amplified in 133 overlapping amplicons, and the restriction patterns of these fragments were compared. The approximate sizes of the genomes investigated based on this whole genome fingerprinting (WGF) analysis ranged from 1139.3-1140.4 kb, with the estimated genome sequence identity of 99.57-99.98% in the homologous genome regions. Restriction target site analysis, detecting the presence of 1773 individual restriction sites found in the reference Nichols genome, revealed a high genome structure similarity of all strains. The unclassified simian Fribourg-Blanc isolate was more closely related to T. p. pertenue than to T. p. pallidum strains. Most of the genetic differences between T. p. pallidum and T. p. pertenue strains were accumulated in six genomic regions. These genome differences likely contribute to the observed differences in pathogenicity between T. p. pallidum and T. p. pertenue strains. These regions of sequence divergence could be used for the molecular detection and discrimination of syphilis and yaws strains.

The peptidoglycan-binding protein FimV promotes assembly of the Pseudomonas aeruginosa type IV pilus secretin

The Pseudomonas aeruginosa inner membrane protein FimV is among several proteins of unknown function required for type IV pilus-mediated twitching motility, arising from extension and retraction of pili from their site of assembly in the inner membrane. The pili transit the periplasm and peptidoglycan (PG) layer, ultimately exiting the cell through the PilQ secretin. Although fimV mutants are nonmotile, they are susceptible to killing by pilus-specific bacteriophage, a hallmark of retractable surface pili. Here we show that levels of recoverable surface pili were markedly decreased in fimV pilT retraction-deficient mutants compared with levels in the pilT control, demonstrating that FimV acts at the level of pilus assembly. Levels of inner membrane assembly subcomplex proteins PilM/N/O/P were decreased in fimV mutants, but supplementation of these components in trans did not restore pilus assembly or motility. Loss of FimV dramatically reduced the levels of the PilQ secretin multimer through which pili exit the cell, in part due to decreased levels of PilQ monomers, while PilF pilotin levels were unchanged. Expression of pilQ in trans in the wild type or fimV mutants increased total PilQ monomer levels but did not alter secretin multimer levels or motility. PG pulldown assays showed that the N terminus of FimV bound PG in a LysM motif-dependent manner, and a mutant with an in-frame chromosomal deletion of the LysM motif had reduced motility, secretin levels, and surface piliation. Together, our data show that FimV's role in pilus assembly is to promote secretin formation and that this function depends upon its PG-binding domain.

Variable numbers of tandem repeats in Plasmodium falciparum genes

Genome variation studies in Plasmodium falciparum have focused on SNPs and, more recently, large-scale copy number polymorphisms and ectopic rearrangements. Here, we examine another source of variation: variable number tandem repeats (VNTRs). Interspersed low complexity features, including the well-studied P. falciparum microsatellite sequences, are commonly classified as VNTRs; however, this study is focused on longer coding VNTR polymorphisms, a small class of copy number variations. Selection against frameshift mutation is a main constraint on tandem repeats (TRs) in coding regions, while limited propagation of TRs longer than 975 nt total length is a minor restriction in coding regions. Comparative analysis of three P. falciparum genomes reveals that more than 9% of all P. falciparum ORFs harbor VNTRs, much more than has been reported for any other species. Moreover, genotyping of VNTR loci in a drug-selected line, progeny of a genetic cross, and 334 field isolates demonstrates broad variability in these sequences. Functional enrichment analysis of ORFs harboring VNTRs identifies stress and DNA damage responses along with chromatin modification activities, suggesting an influence on genome mutability and functional variation. Analysis of the repeat units and their flanking regions in both P. falciparum and Plasmodium reichenowi sequences implicates a replication slippage mechanism in the generation of TRs from an initially unrepeated sequence. VNTRs can contribute to rapid adaptation by localized sequence duplication. They also can confound SNP-typing microarrays or mapping short-sequence reads and therefore must be accounted for in such analyses.

The role of fimV and the importance of its tandem repeat copy number in twitching motility, pigment production, and morphology in Legionella pneumophila

Twitching motility, a flagella-independent type of translocation of bacteria over moist surfaces, requires an array of proteins, including FimV. To investigate the role of this protein in twitching motility in Legionella pneumophila, we have generated a knockout mutant of fimV and characterized its phenotypic effects. In addition to a major reduction in twitching motility, deletion of the fimV gene caused a number of other phenotypic effects including decreased protective pigment formation, and it also affected cell morphology. Since fimV contains a variable number of tandem repeats, which can vary according to the origin of a given strain, we have examined the importance of this variability found within the coding region of this gene. By complementing the knockout strain with constructs containing a different number of this tandem repeat, we have been able to also show that repeat copy number is important in the functioning of this gene.

Twitching motility in Legionella pneumophila

Twitching motility is a form of bacterial translocation over solid or semi-solid surfaces mediated by the extension, tethering, and subsequent retraction of type IV pili. These pili are also known to be involved in virulence, biofilm formation, formation of fruiting bodies, horizontal gene transfer, and protein secretion. We have characterized the presence of twitching motility on agar plates in Legionella pneumophila, the etiological agent of Legionnaires' disease. By examining twitching motility zones, we have demonstrated that twitching motility was dependent on agar thickness/concentration, the chemical composition of the media, the presence of charcoal and cysteine, proximity to other bacteria, and temperature. A knockout mutant of the pilus subunit, pilE, exhibited a total loss of twitching motility at 37 degrees C, but not at 27 degrees C, suggesting either the existence of a compensating pilus subunit or of another twitching motility system in this organism.