GrapeTree: visualization of core genomic relationships among 100,000 bacterial pathogens.

Current methods struggle to reconstruct and visualize the genomic relationships of large numbers of bacterial genomes. GrapeTree facilitates the analyses of large numbers of allelic profiles by a static "GrapeTree Layout" algorithm that supports interactive visualizations of large trees within a web browser window. GrapeTree also implements a novel minimum spanning tree algorithm (MSTree V2) to reconstruct genetic relationships despite high levels of missing data. GrapeTree is a stand-alone package for investigating phylogenetic trees plus associated metadata and is also integrated into EnteroBase to facilitate cutting edge navigation of genomic relationships among bacterial pathogens.

Heterogeneity of penicillin-non-susceptible group B streptococci isolated from a single patient in Germany.

OBJECTIVES: Streptococcus agalactiae [group B streptococci (GBS)] have been considered uniformly susceptible to penicillin. However, increasing reports from Asia and North America are documenting penicillin-non-susceptible GBS (PRGBS) with mutations in pbp genes. Here we report, to the best of our knowledge, the first two PRGBS isolates recovered in Europe (AC-13238-1 and AC-13238-2), isolated from the same patient. METHODS: Two different colony morphologies of GBS were noted from a surgical abscess drainage sample. Both were serotyped and antimicrobial susceptibility testing was performed by different methodologies. High-throughput sequencing was done to compare the isolates at the genomic level, to identify their capsular type and ST, to evaluate mutations in the pbp genes and to compare the isolates with the genomes of other PRGBS isolates sharing the same serotype and ST. RESULTS: Isolates AC-13238-1 and AC-13238-2 presented MICs above the EUCAST and CLSI breakpoints for penicillin susceptibility. Both shared the capsular type Ia operon and ST23. Genomic analysis uncovered differences between the two isolates in seven genes, including altered pbp genes. Deduced amino acid sequences revealed critical substitutions in PBP2X in both isolates. Comparison with serotype Ia clonal complex 23 PRGBS from the USA reinforced the similarity between AC-13238-1 and AC-13238-2, and their divergence from the US strains. CONCLUSIONS: Our results support the in-host evolution of ß-lactam-resistant GBS, with two PRGBS variants being isolated from one patient.

Lactobacillus mulieris sp. nov., a new species of Lactobacillus delbrueckii group.

One Gram-stain-positive, non-motile, non-spore-forming, catalase-negative, and coccobacilli-shaped strain, designated c10Ua161MT, was isolated from a urine sample from a reproductive-age healthy woman. Comparative 16S rRNA gene sequence analysis indicated that strain c10Ua161MT belonged to the genus Lactobacillus. Phylogenetic analysis based on pheS and rpoA gene sequences strongly supported a clade encompassing strains c10Ua161MT and eight other strains from public databases, distinct from currently recognized species of the genus Lactobacillus. In silico Average Nucleotide Identity (ANI) and Genome-to-Genome Distance Calculator (GGDC), showed 87.9 and 34.3 % identity to the closest relative Lactobacillus jensenii, respectively. The major fatty acids of strain c10Ua161MT were C18 : 1ω9c (65.0%), C16 : 0 (17.8%), and summed feature 8 (10.2 %; comprising C18 : 1ω7c, and/or C18 : 1ω6c). The DNA G+C content of the strains is 34.2 mol%. On the basis of data presented here, strain c10Ua161MT represents a novel species of the genus Lactobacillus, for which the name Lactobacillus mulieris sp. nov. is proposed. The type strain is c10Ua161MT (=CECT 9755T=DSM 108704T).

Streptococcus pneumoniae Serotype 3 in Mexico (1994 to 2017): Decrease of the Unusual Clonal Complex 4909 Lineage following PCV13 Introduction.

Streptococcus pneumoniae expressing serotype 3 has a high virulence and a high case fatality ratio. Most studies of serotype 3 pneumococci have focused on a single lineage, the widespread sequence type 180 (ST180). To evaluate the serotype 3 lineages causing infections in Mexico, we characterized 196 isolates recovered from 1994 to 2017. The isolates were mostly susceptible to all antimicrobials tested. A single meningitis isolate was resistant to penicillin, and the resistance to erythromycin was 5.2%. The isolates represented the widely disseminated clonal complex 180 (CC180; n = 140), the unusual CC4909 (n = 42), CC260 (n = 11), and a few singletons (n = 3). CC260 was less frequent among pneumococcal invasive disease isolates than CC180 and CC4909 (P = 0.015). There was a decrease of CC4909 (P < 0.001) following PCV13 introduction (2012 to 2017). The CC4909 isolates were represented mostly by ST1119 (n = 40), seemingly having a restricted geographic origin, with isolates in the PubMLST database having been recovered only in Mexico, the United States, and Germany. A genomic analysis of publicly available genomes showed that ST1119 isolates have less than 32% similarity with ST180 isolates, indicating that these lineages are more separated than revealed by traditional multilocus sequence typing. Considering the suggestions of a lower efficacy of the 13-valent pneumococcal conjugate vaccine against serotype 3, the different dynamics of the two major serotype 3 lineages in Mexico following the introduction of PCV13 should be closely monitored.

PHYLOViZ Online: web-based tool for visualization, phylogenetic inference, analysis and sharing of minimum spanning trees.

High-throughput sequencing methods generated allele and single nucleotide polymorphism information for thousands of bacterial strains that are publicly available in online repositories and created the possibility of generating similar information for hundreds to thousands of strains more in a single study. Minimum spanning tree analysis of allelic data offers a scalable and reproducible methodological alternative to traditional phylogenetic inference approaches, useful in epidemiological investigations and population studies of bacterial pathogens. PHYLOViZ Online was developed to allow users to do these analyses without software installation and to enable easy accessing and sharing of data and analyses results from any Internet enabled computer. PHYLOViZ Online also offers a RESTful API for programmatic access to data and algorithms, allowing it to be seamlessly integrated into any third party web service or software. PHYLOViZ Online is freely available at https://online.phyloviz.net.

Increasing macrolide resistance among Streptococcus agalactiae causing invasive disease in non-pregnant adults was driven by a single capsular-transformed lineage, Portugal, 2009 to 2015.

We characterised Lancefield group B streptococcal (GBS) isolates causing invasive disease among non-pregnant adults in Portugal between 2009 and 2015. All isolates (n = 555) were serotyped, assigned to clonal complexes (CCs) by multilocus sequence typing and characterised by surface protein and pilus island gene profiling. Antimicrobial susceptibility was tested by disk diffusion and resistance genotypes identified by PCR. Overall, serotype Ia was most frequent in the population (31%), followed by serotypes Ib (24%) and V (18%). Serotype Ib increased significantly throughout the study period (p < 0.001) to become the most frequent serotype after 2013. More than 40% of isolates clustered in the CC1/alp3/PI-1+PI-2a genetic lineage, including most isolates of serotypes Ib (n = 110) and V (n = 65). Erythromycin and clindamycin resistance rates were 35% and 34%, respectively, both increasing from 2009 to 2015 (p < 0.010) and associated with CC1 and serotype Ib (p < 0.001). The Ib/CC1 lineage probably resulted from acquisition of the type Ib capsular operon in a single recombination event by a representative of the V/CC1 macrolide-resistant lineage. Expansion of the new serotype Ib/CC1 lineage resulted in increased macrolide resistance in GBS, causing invasive disease among adults in Portugal. The presence of this clone elsewhere may predict more widespread increase in resistance.

A widespread family of polymorphic toxins encoded by temperate phages.

BACKGROUND: Polymorphic toxins (PTs) are multi-domain bacterial exotoxins belonging to distinct families that share common features in terms of domain organization. PTs are found in all major bacterial clades, including many toxic effectors of type V and type VI secretion systems. PTs modulate the dynamics of microbial communities by killing or inhibiting the growth of bacterial competitors lacking protective immunity proteins. RESULTS: In this work, we identified a novel widespread family of PTs, named MuF toxins, which were exclusively encoded within temperate phages and their prophages. By analyzing the predicted proteomes of 1845 bacteriophages and 2464 bacterial genomes, we found that MuF-containing proteins were frequently part of the DNA packaging module of tailed phages. Interestingly, MuF toxins were abundant in the human gut microbiome. CONCLUSIONS: Our results uncovered the presence of the MuF toxin family in the temperate phages of Firmicutes. The MuF toxin family is likely to play an important role in the ecology of the human microbiota where pathogens and commensal species belonging to the Firmicutes are abundant. We propose that MuF toxins could be delivered by phages into host bacteria and either influence the lysogeny decision or serve as bacterial weapons by inhibiting the growth of competing bacteria.

Citrobacter portucalensis sp. nov., isolated from an aquatic sample.

A Gram-stain-negative strain, A60T, isolated from a water well sample in Portugal, was characterized phenotypically, genotypically and phylogenetically. Phylogenetic analysis based on 16S rRNA gene sequences indicated that strain A60T belonged to the genus Citrobacter, and recN gene phylogeny revealed one strongly supported clade encompassing strain A60T and 13 other strains from public databases, distinct from currently recognized species of the genus Citrobacter. Furthermore, multilocus sequence analysis (MLSA) based on concatenated partial fusA, leuS, pyrG and rpoB sequences confirmed the classification obtained with the recN sequence. In silico genomic comparisons, including average nucleotide identity (ANI) and the genome-to-genome distance calculator (GGDC), showed 94.6 % and 58.4 % identity to the closest relative Citrobacter freundii ATCC 8090T, respectively. The ability to metabolize different compounds further discriminated strain A60T from other species of the genus Citrobacter. The G+C content of strain A60T is 52.0 %. The results obtained support the description of a novel species within the genus Citrobacter, for which the name Citrobacter portucalensis sp. nov. is proposed, with the type strain A60T (=DSM 104542T=CECT 9236T).

Automatic determination of NET (neutrophil extracellular traps) coverage in fluorescent microscopy images.

MOTIVATION: Neutrophil extracellular traps (NETs) are believed to be essential in controlling several bacterial pathogens. Quantification of NETs in vitro is an important tool in studies aiming to clarify the biological and chemical factors contributing to NET production, stabilization and degradation. This estimation can be performed on the basis of fluorescent microscopy images using appropriate labelings. In this context, it is desirable to automate the analysis to eliminate both the tedious process of manual annotation and possible operator-specific biases. RESULTS: We propose a framework for the automated determination of NET content, based on visually annotated images which are used to train a supervised machine-learning method. We derive several methods in this framework. The best results are obtained by combining these into a single prediction. The overall Q(2) of the combined method is 93%. By having two experts label part of the image set, we were able to compare the performance of the algorithms to the human interoperator variability. We find that the two operators exhibited a very high correlation on their overall assessment of the NET coverage area in the images (R(2) is 97%), although there were consistent differences in labeling at pixel level (Q(2), which unlike R(2) does not correct for additive and multiplicative biases, was only 89%). AVAILABILITY AND IMPLEMENTATION: Open source software (under the MIT license) is available at https://github.com/luispedro/Coelho2015_NetsDetermination for both reproducibility and application to new data.

Comparing Matrix-Assisted Laser Desorption Ionization-Time of Flight Mass Spectrometry and Phenotypic and Molecular Methods for Identification of Species within the Streptococcus anginosus Group.

The heterogeneity of members of the Streptococcus anginosus group (SAG) has traditionally hampered their correct identification. Recently, the group was subdivided into 6 taxa whose prevalence among human infections is poorly described. We evaluated the accuracy of the Rapid ID32 Strep test, matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS), and a PCR multiplex method to identify 212 SAG isolates recovered from human infections to the species and subspecies level by using multilocus sequence analysis (MLSA) as the gold standard. We also determined the antimicrobial susceptibilities of the isolates. Representatives of all SAG taxa were found among our collection. MALDI-TOF MS and the Rapid ID32 Strep test correctly identified 92% and 68% of the isolates to the species level, respectively, but showed poor performance at the subspecies level, and the latter was responsible for major identification errors. The multiplex PCR method results were in complete agreement with the MLSA identifications but failed to distinguish the subspecies Streptococcus constellatus subsp. pharyngis and S. constellatus subsp. viborgensis. A total of 145 MLSA sequence types were present in our collection, indicating that within each taxon a number of different lineages are capable of causing infection. Significant antibiotic resistance was observed only to tetracycline, erythromycin, and clindamycin and was present in most taxa. MALDI-TOF MS is a reliable method for routine SAG species identification, while the need for identification to the subspecies level is not clearly established.

Spontaneous Genomic Variation as a Survival Strategy of Nosocomial Staphylococcus haemolyticus.

Staphylococcus haemolyticus is one of the most important nosocomial human pathogens frequently isolated in bloodstream and medical device-related infections. However, its mechanisms of evolution and adaptation are still poorly explored. To characterize the strategies of genetic and phenotypic diversity in S. haemolyticus, we analyzed an invasive strain for genetic and phenotypic stability after serial passage in vitro in the absence and presence of beta-lactam antibiotics. We performed pulsed-field gel electrophoresis (PFGE) of the culture and analyzed five colonies at seven time points during stability assays for beta-lactam susceptibility, hemolysis, mannitol fermentation, and biofilm production. We compared their whole genomes and performed phylogenetic analysis based on core single-nucleotide polymorphisms (SNPs). We observed a high instability in the PFGE profiles at the different time points in the absence of antibiotic. Analysis of WGS data for individual colonies showed the occurrence of six large-scale genomic deletions within the oriC environ, smaller deletions in non-oriC environ regions, and nonsynonymous mutations in clinically relevant genes. The regions of deletion and point mutations included genes encoding amino acid and metal transporters, resistance to environmental stress and beta-lactams, virulence, mannitol fermentation, metabolic processes, and insertion sequence (IS) elements. Parallel variation was detected in clinically significant phenotypic traits such as mannitol fermentation, hemolysis, and biofilm formation. In the presence of oxacillin, PFGE profiles were overall stable over time and mainly corresponded to a single genomic variant. Our results suggest that S. haemolyticus populations are composed of subpopulations of genetic and phenotypic variants. The maintenance of subpopulations in different physiological states may be a strategy to adapt rapidly to stress situations imposed by the host, particularly in the hospital environment. IMPORTANCE The introduction of medical devices and antibiotics into clinical practice have substantially improved patient quality of life and contributed to extended life expectancy. One of its most cumbersome consequences was the emergence of medical device-associated infections caused by multidrug-resistant and opportunistic bacteria such as Staphylococcus haemolyticus. However, the reason for this bacterium's success is still elusive. We found that in the absence of environmental stresses, S. haemolyticus can spontaneously produce subpopulations of genomic and phenotypic variants with deletions/mutations in clinically relevant genes. However, when exposed to selective pressures, such as the presence of antibiotics, a single genomic variant will be recruited and become dominant. We suggest that the maintenance of these cell subpopulations in different physiological states is an extremely effective strategy to adapt to stresses imposed by the host or the infection environment and might contribute for S. haemolyticus survival and persistence in the hospital.

LMAS: evaluating metagenomic short de novo assembly methods through defined communities.

BACKGROUND: The de novo assembly of raw sequence data is key in metagenomic analysis. It allows recovering draft genomes from a pool of mixed raw reads, yielding longer sequences that offer contextual information and provide a more complete picture of the microbial community. FINDINGS: To better compare de novo assemblers for metagenomic analysis, LMAS (Last Metagenomic Assembler Standing) was developed as a flexible platform allowing users to evaluate assembler performance given known standard communities. Overall, in our test datasets, k-mer De Bruijn graph assemblers outperformed the alternative approaches but came with a greater computational cost. Furthermore, assemblers branded as metagenomic specific did not consistently outperform other genomic assemblers in metagenomic samples. Some assemblers still in use, such as ABySS, MetaHipmer2, minia, and VelvetOptimiser, perform relatively poorly and should be used with caution when assembling complex samples. Meaningful strain resolution at the single-nucleotide polymorphism level was not achieved, even by the best assemblers tested. CONCLUSIONS: The choice of a de novo assembler depends on the computational resources available, the replicon of interest, and the major goals of the analysis. No single assembler appeared an ideal choice for short-read metagenomic prokaryote replicon assembly, each showing specific strengths. The choice of metagenomic assembler should be guided by user requirements and characteristics of the sample of interest, and LMAS provides an interactive evaluation platform for this purpose. LMAS is open source, and the workflow and its documentation are available at https://github.com/B-UMMI/LMAS and https://lmas.readthedocs.io/, respectively.

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

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

Draft Genome Sequence of Escherichia coli DSM 12242, a Highly Efficient Host Strain for the Isolation of Somatic Coliphages.

Here, we report the draft genome sequence and characterization of the commercial strain Escherichia coli DSM 12242 (=ATCC 13706/60 =NZRM 3262) derived from strain C (ATCC 13706), which is suitable for the isolation of coliphages from environmental and clinical samples.

Natural Occurrence of Escherichia coli-Infecting Bacteriophages in Clinical Samples.

The interaction between bacteriophages, bacteria and the human host as a tripartite system has recently captured attention. The taxonomic diversity of bacteriophages, as a natural parasite of bacteria, still remains obscure in human body biomes, representing a so-called "viral dark matter." Here, we isolated and characterized coliphages from blood, urine and tracheal aspirates samples collected at a tertiary care hospital in Austria. Phages were more often isolated from blood, followed by urine and tracheal aspirates. Phylogenetic analysis and genome comparisons allowed the identification of phages belonging to the Tunavirinae subfamily, and to the Peduovirus and Tequintavirus genera. Tunavirinae phages cluster together and are found in samples from 14 patients, suggesting their prevalence across a variety of human samples. When compared with other phage genomes, the highest similarity level was at 87.69% average nucleotide identity (ANI), which suggests that these are in fact a newly isolated phage species. Tequintavirus phages share a 95.90% with phage 3_29, challenging the ANI threshold currently accepted to differentiate phage species. The isolated phages appear to be virulent, with the exception of the Peduovirus members, which are integrative and seem to reside as prophages in bacterial genomes.

Origin, evolution, and distribution of the molecular machinery for biosynthesis of sialylated lipooligosaccharide structures in Campylobacter coli.

Campylobacter jejuni and Campylobacter coli are the most common cause of bacterial gastroenteritis worldwide. Additionally, C. jejuni is the most common bacterial etiological agent in the autoimmune Guillain-Barré syndrome (GBS). Ganglioside mimicry by C. jejuni lipooligosaccharide (LOS) is the triggering factor of the disease. LOS-associated genes involved in the synthesis and transfer of sialic acid (glycosyltranferases belonging to family GT-42) are essential in C. jejuni to synthesize ganglioside-like LOS. Despite being isolated from GBS patients, scarce genetic evidence supports C. coli role in the disease. In this study, through data mining and bioinformatics analysis, C. coli is shown to possess a larger GT-42 glycosyltransferase repertoire than C. jejuni. Although GT-42 glycosyltransferases are widely distributed in C. coli population, only a fraction of C. coli strains (1%) are very likely able to express ganglioside mimics. Even though the activity of C. coli specific GT-42 enzymes and their role in shaping the bacterial population are yet to be explored, evidence presented herein suggest that loss of function of some LOS-associated genes occurred during agriculture niche adaptation.

chewBBACA: A complete suite for gene-by-gene schema creation and strain identification.

Gene-by-gene approaches are becoming increasingly popular in bacterial genomic epidemiology and outbreak detection. However, there is a lack of open-source scalable software for schema definition and allele calling for these methodologies. The chewBBACA suite was designed to assist users in the creation and evaluation of novel whole-genome or core-genome gene-by-gene typing schemas and subsequent allele calling in bacterial strains of interest. chewBBACA performs the schema creation and allele calls on complete or draft genomes resulting from de novo assemblers. The chewBBACA software uses Python 3.4 or higher and can run on a laptop or in high performance clusters making it useful for both small laboratories and large reference centers. ChewBBACA is available at https://github.com/B-UMMI/chewBBACA.

Epidemiological Surveillance and Typing Methods to Track Antibiotic Resistant Strains Using High Throughput Sequencing.

High-Throughput Sequencing (HTS) technologies transformed the microbial typing and molecular epidemiology field by providing the cost-effective ability for researchers to probe draft genomes, not only for epidemiological markers but also for antibiotic resistance and virulence determinants. In this chapter, we provide protocols for the analysis of HTS data for the determination of multilocus sequence typing (MLST) information and for determining presence or absence of antibiotic resistance genes.

Clonal and serotype dynamics of serogroup 6 isolates causing invasive pneumococcal disease in Portugal: 1999-2012.

Although serogroup 6 was among the first to be recognized among Streptococcus pneumoniae, several new serotypes were identified since the introduction of pneumococcal conjugate vaccines (PCVs). A decrease of the 6B-2 variant among invasive pneumococcal disease (IPD), but not 6B-1, was noted post conjugate vaccine introduction, underpinned by a decrease of CC273 isolates. Serotype 6C was associated with adult IPD and increased in this age group representing two lineages (CC315 and CC395), while the same lineages expressed other serogroup 6 serotypes in children. Taken together, these findings suggest a potential cross-protection of PCVs against serotype 6C IPD among vaccinated children but not among adults. Serotype 6A became the most important serogroup 6 serotype in children but it decreased in adult IPD. No other serogroup 6 serotypes were detected, so available phenotypic or simple genotypic assays remain adequate for distinguishing serotypes within serogroup 6 isolates.

Serum lipid alterations in GBA-associated Parkinson's disease.

INTRODUCTION: Mutations in the GBA gene, encoding for the lysosomal enzyme glucocerebrosidase, are associated with Gaucher disease. Alterations in plasma sphingolipids have been reported in Gaucher, and similarly in brain extracts in Lewy body disease. As GBA mutations are prevalent risk factors for Parkinson's disease and overlap of molecular pathways are presumable, here we assessed the lipid profiles in Parkinson's patients with and without GBA mutations. METHODS: We sequenced all GBA exons in 415 Parkinson's patients, previously genotyped for LRRK2. 64 patients (29 GBA positive vs. 35 non-GBA-carriers including 18 LRRK2 positive and 17 non-mutated) were analyzed for chitotriosidase activity and for the concentration of 40 lipid classes using HPLC-MS. RESULTS: 29/415 patients (6.9%) carried 8 different GBA mutations associated with Gaucher or Parkinson's, including one novel mutation. Chitotriosidase activity was similar across the genetic groups, while the levels of key lipids were altered in GBA mutation carriers: Monohexosylceramide, Ceramide and Sphingomyelin were elevated; while Phosphatidic acid (PA), Phosphatidylethanolamine (PE), Plasmalogen phosphatidylethanolamine (PEp) and Acyl Phosphatidylglycerol (AcylPG) were decreased. CONCLUSION: The results suggest an important role for these lipids in GBA mediated Parkinson's disease and assist in the identification of common pathways between Gaucher and Parkinson's. Ultimately, our findings may lead to the identification of novel biomarkers for individuals at increased risk of developing Parkinson's disease.