An extended multi-locus molecular typing schema for Streptococcus pneumoniae demonstrates that a limited number of capsular switch events is responsible for serotype heterogeneity of closely related strains from different countries

MLSTar: automatic multilocus sequence typing of bacterial genomes in R

Multilocus sequence typing (MLST) is a standard tool in population genetics and bacterial epidemiology that assesses the genetic variation present in a reduced number of housekeeping genes (typically seven) along the genome. This methodology assigns arbitrary integer identifiers to genetic variations at these loci which allows us to efficiently compare bacterial isolates using allele-based methods. Now, the increasing availability of whole-genome sequences for hundreds to thousands of strains from the same bacterial species has allowed us to apply and extend MLST schemes by automatic extraction of allele information from the genomes. The PubMLST database is the most comprehensive resource of described schemes available for a wide variety of species. Here we present MLSTar as the first R package that allows us to (i) connect with the PubMLST database to select a target scheme, (ii) screen a desired set of genomes to assign alleles and sequence types, and (iii) interact with other widely used R packages to analyze and produce graphical representations of the data. We applied MLSTar to analyze more than 2,500 bacterial genomes from different species, showing great accuracy, and comparable performance with previously published command-line tools. MLSTar can be freely downloaded from http://github.com/iferres/MLSTar.

Efficient Inference of Recent and Ancestral Recombination within Bacterial Populations

Prokaryotic evolution is affected by horizontal transfer of genetic material through recombination. Inference of an evolutionary tree of bacteria thus relies on accurate identification of the population genetic structure and recombination-derived mosaicism. Rapidly growing databases represent a challenge for computational methods to detect recombinations in bacterial genomes. We introduce a novel algorithm called fastGEAR which identifies lineages in diverse microbial alignments, and recombinations between them and from external origins. The algorithm detects both recent recombinations (affecting a few isolates) and ancestral recombinations between detected lineages (affecting entire lineages), thus providing insight into recombinations affecting deep branches of the phylogenetic tree. In simulations, fastGEAR had comparable power to detect recent recombinations and outstanding power to detect the ancestral ones, compared with state-of-the-art methods, often with a fraction of computational cost. We demonstrate the utility of the method by analyzing a collection of 616 whole-genomes of a recombinogenic pathogen Streptococcus pneumoniae, for which the method provided a high-resolution view of recombination across the genome. We examined in detail the penicillin-binding genes across the Streptococcus genus, demonstrating previously undetected genetic exchanges between different species at these three loci. Hence, fastGEAR can be readily applied to investigate mosaicism in bacterial genes across multiple species. Finally, fastGEAR correctly identified many known recombination hotspots and pointed to potential new ones. Matlab code and Linux/Windows executables are available at https://users.ics.aalto.fi/~pemartti/fastGEAR/ (last accessed February 6, 2017).

Long Persistence of a Streptococcus pneumoniae 23F Clone in a Cystic Fibrosis Patient

Streptococcus pneumoniae is a common resident in the human nasopharynx. However, carriage can result in severe diseases due to a unique repertoire of pathogenicity factors that are rare in closely related commensal streptococci. We investigated a penicillin-resistant S. pneumoniae clone of serotype 23F isolated from a cystic fibrosis patient on multiple occasions over an unusually long period of over 3 years that was present without causing disease. Genome comparisons revealed an apparent nonfunctional pneumococcus-specific gene encoding a hyaluronidase, supporting the view that this enzyme adds to the virulence potential of the bacterium. The 23F clone harbored unique mosaic genes encoding penicillin resistance determinants, the product of horizontal gene transfer involving the commensal S. mitis as donor species. Sequences identical to one such mosaic gene were identified in an S. mitis strain from the same patient, suggesting that in this case S. pneumoniae played the role of donor.

Streptococcus pneumoniae isolates of serotype 23F with intermediate penicillin resistance were recovered on seven occasions over a period of 37 months from a cystic fibrosis patient in Berlin. All isolates expressed the same multilocus sequence type (ST), ST10523. The genome sequences of the first and last isolates, D122 and D141, revealed the absence of two phage-related gene clusters compared to the genome of another ST10523 strain, D219, isolated earlier at a different place in Germany. Genomes of all three strains carried the same novel mosaic penicillin-binding protein (PBP) genes, pbp2x, pbp2b, and pbp1a; these genes were distinct from those of other penicillin-resistant S. pneumoniae strains except for pbp1a of a Romanian S. pneumoniae isolate. All PBPs contained mutations that have been associated with the penicillin resistance phenotype. Most interestingly, a mosaic block identical to an internal pbp2x sequence of ST10523 was present in pbp2x of Streptococcus mitis strain B93-4, which was isolated from the same patient. This suggests interspecies gene transfer from S. pneumoniae to S. mitis within the host. Nearly all genes expressing surface proteins, which represent major virulence factors of S. pneumoniae and are typical for this species, were present in the genome of ST10523. One exception was the hyaluronidase gene hlyA, which contained a 12-nucleotide deletion within the promoter region and an internal stop codon. The lack of a functional hyaluronidase might contribute to the ability to persist in the host for an unusually long period of time.

IMPORTANCEStreptococcus pneumoniae is a common resident in the human nasopharynx. However, carriage can result in severe diseases due to a unique repertoire of pathogenicity factors that are rare in closely related commensal streptococci. We investigated a penicillin-resistant S. pneumoniae clone of serotype 23F isolated from a cystic fibrosis patient on multiple occasions over an unusually long period of over 3 years that was present without causing disease. Genome comparisons revealed an apparent nonfunctional pneumococcus-specific gene encoding a hyaluronidase, supporting the view that this enzyme adds to the virulence potential of the bacterium. The 23F clone harbored unique mosaic genes encoding penicillin resistance determinants, the product of horizontal gene transfer involving the commensal S. mitis as donor species. Sequences identical to one such mosaic gene were identified in an S. mitis strain from the same patient, suggesting that in this case S. pneumoniae played the role of donor.

Genotyping Streptococcus pneumoniae

ABSTRACT 

Streptococcus pneumoniae is a potentially deadly human pathogen associated with high morbidity, mortality and global economic burden. The universally used bacterial genotyping methods are multilocus sequence typing and pulsed field gel electrophoresis. However, another highly discriminatory, rapid and less expensive genotyping technique, multilocus variable number of tandem repeat analysis (MLVA), has been developed. Unfortunately, no universal MLVA protocol exists, and some MLVA protocols do not amplify certain loci for all pneumococcal serotypes, leaving genotyping profiles incomplete. A number of other genotyping or characterization methods have been developed and will be discussed. This review examines the various protocols for genotyping S. pneumoniae and highlights the current direction technology and research is heading to understand this bacterium.

High-resolution melt PCR analysis for genotyping of Ureaplasma parvum isolates directly from clinical samples

Ureaplasma sp. infection in neonates and adults underlies a variety of disease pathologies. Of the two human Ureaplasma spp., Ureaplasma parvum is clinically the most common. We have developed a high-resolution melt (HRM) PCR assay for the differentiation of the four serovars of U. parvum in a single step. Currently U. parvum strains are separated into four serovars by sequencing the promoter and coding region of the multiple-banded antigen (MBA) gene. We designed primers to conserved sequences within this region for PCR amplification and HRM analysis to generate reproducible and distinct melt profiles that distinguish clonal representatives of serovars 1, 3, 6, and 14. Furthermore, our HRM PCR assay could classify DNA extracted from 74 known (MBA-sequenced) test strains with 100% accuracy. Importantly, HRM PCR was also able to identify U. parvum serovars directly from 16 clinical swabs. HRM PCR performed with DNA consisting of mixtures of combined known serovars yielded profiles that were easily distinguished from those for single-serovar controls. These profiles mirrored clinical samples that contained mixed serovars. Unfortunately, melt curve analysis software is not yet robust enough to identify the composition of mixed serovar samples, only that more than one serovar is present. HRM PCR provides a single-step, rapid, cost-effective means to differentiate the four serovars of U. parvum that did not amplify any of the known 10 serovars of Ureaplasma urealyticum tested in parallel. Choice of reaction reagents was found to be crucial to allow sufficient sensitivity to differentiate U. parvum serovars directly from clinical swabs rather than requiring cell enrichment using microbial culture techniques.

Sequence analysis of 96 genomic regions identifies distinct evolutionary lineages within CC156, the largest Streptococcus pneumoniae clonal complex in the MLST database

Multi-Locus Sequence Typing (MLST) of Streptococcus pneumoniae is based on the sequence of seven housekeeping gene fragments. The analysis of MLST allelic profiles by eBURST allows the grouping of genetically related strains into Clonal Complexes (CCs) including those genotypes with a common descent from a predicted ancestor. However, the increasing use of MLST to characterize S. pneumoniae strains has led to the identification of a large number of new Sequence Types (STs) causing the merger of formerly distinct lineages into larger CCs. An example of this is the CC156, displaying a high level of complexity and including strains with allelic profiles differing in all seven of the MLST loci, capsular type and the presence of the Pilus Islet-1 (PI-1). Detailed analysis of the CC156 indicates that the identification of new STs, such as ST4945, induced the merging of formerly distinct clonal complexes. In order to discriminate the strain diversity within CC156, a recently developed typing schema, 96-MLST, was used to analyse 66 strains representative of 41 different STs. Analysis of allelic profiles by hierarchical clustering and a minimum spanning tree identified ten genetically distinct evolutionary lineages. Similar results were obtained by phylogenetic analysis on the concatenated sequences with different methods. The identified lineages are homogenous in capsular type and PI-1 presence. ST4945 strains were unequivocally assigned to one of the lineages. In conclusion, the identification of new STs through an exhaustive analysis of pneumococcal strains from various laboratories has highlighted that potentially unrelated subgroups can be grouped into a single CC by eBURST. The analysis of additional loci, such as those included in the 96-MLST schema, will be necessary to accurately discriminate the clonal evolution of the pneumococcal population.