kSNP3.0: SNP detection and phylogenetic analysis of genomes without genome alignment or reference genome

Assessment of the local clonal spread of Streptococcus pneumoniae serotype 12F caused invasive pneumococcal diseases among children and adults

BACKGROUND

We conducted active surveillance to elucidate the distribution of Streptococcus pneumoniae serotypes causing invasive pneumococcal disease (IPD) and clarified the genetic relatedness among the isolates in Kobe City, Japan.

METHODS

Forty-five IPD-causing S. pneumoniae strains were analyzed from March 2016 to May 2018 through active surveillance in Kobe City, Hyogo, Japan. Serotypes were determined by multiplex serotyping PCR and the Quellung reaction with pneumococcal antisera. Fourteen Sp12F strains were subjected to whole-genome sequencing (WGS).

RESULTS

Among 45 isolates, the most frequent serotypes were 12F (n=14, 31%), 24F (n=5, 11%), and 10A (n=4, 9%). Multilocus sequence typing (MLST) analysis of 14 isolates of Sp12F divided them into ST4846 (n=4) and ST6495 (n=10). WGS showed clonality of the 10 isolates of ST6495, with only 13 single nucleotide polymorphisms in the genomes. Meanwhile, ST4846 strains in Kobe differed from only the outbreak strains of Sp12F ST4846 in Tsuruoka, Japan, reported on 2018.

CONCLUSIONS

Serotype monitoring showed Sp12F to be the predominant serotype in Kobe, and WGS revealed the clonal spread of Sp12F ST6495 in this city. Thus, the spread of Sp12F could become a serious public health problem in Japan, warranting thorough monitoring in future.

Increase in Invasive Group A Streptococcal Disease and Emergence of Mucoid Strains in a Pediatric Population: February-June 2017

Background

Infection with group A Streptococcus (GAS) can cause severe systemic and locally invasive disease. Invasive group A streptococcal (iGAS) disease incidence varies both seasonally and year-to-year, and it may exhibit clustered outbreaks. We observed an upswing in iGAS cases at a tertiary care Children’s Hospital, prompting further characterization of local iGAS disease.

Methods

Cases of iGAS disease were abstracted from the medical record by manual chart review of all positive screening tests and cultures for GAS over a 4-year span. Incidence rates per 1000 hospital admissions and per 100 positive GAS tests were calculated and compared. Selected isolates were further characterized by whole-genome sequencing.

Results

Significant year-to-year differences in per-admission iGAS incidence rate were observed in February and June, although per-positive test incidence rates were not significantly different. Whole-genome sequencing revealed 2 dominant serotypes—emm3 and emm6—with high rates of mucoid phenotype and systemic bacteremia.

Conclusions

We document a significant but transient increase in iGAS disease incidence in 2 months of 2017. Genome sequencing revealed 2 dominant serotypes associated with mucoid phenotypes and severe disease, highlighting the dynamic nature of iGAS disease pattern.

The use of next generation sequencing for improving food safety: Translation into practice

Next Generation Sequencing (NGS) combined with powerful bioinformatic approaches are revolutionising food microbiology. Whole genome sequencing (WGS) of single isolates allows the most detailed comparison possible hitherto of individual strains. The two principle approaches for strain discrimination, single nucleotide polymorphism (SNP) analysis and genomic multi-locus sequence typing (MLST) are showing concordant results for phylogenetic clustering and are complementary to each other. Metabarcoding and metagenomics, applied to total DNA isolated from either food materials or the production environment, allows the identification of complete microbial populations. Metagenomics identifies the entire gene content and when coupled to transcriptomics or proteomics, allows the identification of functional capacity and biochemical activity of microbial populations. The focus of this review is on the recent use and future potential of NGS in food microbiology and on current challenges. Guidance is provided for new users, such as public health departments and the food industry, on the implementation of NGS and how to critically interpret results and place them in a broader context. The review aims to promote the broader application of NGS technologies within the food industry as well as highlight knowledge gaps and novel applications of NGS with the aim of driving future research and increasing food safety outputs from its wider use.

Assembly and Characterization of a Pathogen Strain Collection for Produce Safety Applications: Pre-growth Conditions Have a Larger Effect on Peroxyacetic Acid Tolerance Than Strain Diversity

Effective control of foodborne pathogens on produce requires science-based validation of interventions and control strategies, which typically involves challenge studies with a set of bacterial strains representing the target pathogens or appropriate surrogates. In order to facilitate these types of studies, a produce-relevant strain collection was assembled to represent strains from produce outbreaks or pre-harvest environments, including Listeria monocytogenes (n = 11), Salmonella enterica (n = 23), shiga-toxin producing Escherichia coli (STEC) (n = 13), and possible surrogate organisms (n = 8); all strains were characterized by whole genome sequencing (WGS). Strain diversity was assured by including the 10 most common S. enterica serotypes, L. monocytogenes lineages I-IV, and E. coli O157 as well as selected "non-O157" STEC serotypes. As it has previously been shown that strains and genetic lineages of a pathogen may differ in their ability to survive different stress conditions, a subset of representative strains for each "pathogen group" (e.g., Salmonella, STEC) was selected and assessed for survival of exposure to peroxyacetic acid (PAA) using strains pre-grown under different conditions including (i) low pH, (ii) high salt, (iii) reduced water activity, (iv) different growth phases, (v) minimal medium, and (vi) different temperatures (21°C, 37°C). The results showed that across the three pathogen groups pre-growth conditions had a larger effect on bacterial reduction after PAA exposure as compared to strain diversity. Interestingly, bacteria exposed to salt stress (4.5% NaCl) consistently showed the least reduction after exposure to PAA; however, for STEC, strains pre-grown at 21°C were as tolerant to PAA exposure as strains pre-grown under salt stress. Overall, our data suggests that challenge studies conducted with multi-strain cocktails (pre-grown under a single specific condition) may not necessarily reflect the relevant phenotypic range needed to appropriately assess different intervention strategies.

Catabolism of Nucleic Acids by a Cystic Fibrosis Pseudomonas aeruginosa Isolate: An Adaptive Pathway to Cystic Fibrosis Sputum Environment

Pseudomonas aeruginosa is a major cause of morbidity and mortality in patients with cystic fibrosis (CF). We undertook Biolog Phenotype Microarray testing of P. aeruginosa CF isolates to investigate their catabolic capabilities compared to P. aeruginosa laboratory strains PAO1 and PA14. One strain, PASS4, displayed an unusual phenotype, only showing strong respiration on adenosine and inosine. Further testing indicated that PASS4 could grow on DNA as a sole carbon source, with a higher biomass production than PAO1. This suggested that PASS4 was specifically adapted to metabolize extracellular DNA, a substrate present at high concentrations in the CF lung. Transcriptomic and proteomic profiling of PASS4 and PAO1 when grown with DNA as a sole carbon source identified a set of upregulated genes, including virulence and host-adaptation genes. PASS4 was unable to utilize N-Acetyl-D-glucosamine, and when we selected PASS4 mutants able to grow on this carbon source, they also displayed a gain in ability to catabolize a broad range of other carbon sources. Genome sequencing of the mutants revealed they all contained mutations within the purK gene, encoding a key protein in the de novo purine biosynthesis pathway. This suggested that PASS4 was a purine auxotroph. Growth assays in the presence of 2 mM adenosine and the complementation of PASS4 with an intact purK gene confirmed this conclusion. Purine auxotrophy may represent a viable microbial strategy for adaptation to DNA-rich environments such as the CF lung.

Mycobacterium caprae Infection in Captive Borneo Elephant, Japan

In 2016, disseminated tuberculosis caused by Mycobacterium caprae was diagnosed in a captive Borneo elephant in Japan. The bacterium was initially identified from clinical isolates. An isolate collected during a relapse showed isoniazid monoresistance and a codon 315 katG mutation.

Analysis of Genetic Characterization and Clonality of Legionella pneumophila Isolated from Cooling Towers in Japan

We investigated the genetic characteristics of 161 Legionella pneumophila strains isolated over a period of 10 years from cooling towers in Japan. Minimum spanning tree analysis based on the sequence-based typing (SBT) of them identified three clonal complexes (CCs); CC1 (105/161, 65.2%), CC2 (22 /161, 13.7%), and CC3 (20/161, 12.4%). CC1 was formed by serogroup (SG) 1 and SG7, whereas CC2 was mainly formed by SG1. All of the CC3 isolates except two strains were SG13. The major sequence types (STs) in CC1 and CC2 were ST1 (88/105, 83.8%) and ST154 (15/22, 68.2%), respectively. These STs are known as typical types of L. pneumophila SG1 in Japanese cooling tower. Additionally, we identified 15 strains of ST2603 as the major type in CC3. This ST has not been reported in Japanese cooling tower. Whole genome sequencing (WGS) analysis of the representative strains in the three CCs, which were isolated from various cooling towers over the 10 years, elucidated high clonal population of L. pneumophila in Japanese cooling tower. Moreover, it revealed that the strains of CC2 are phylogenetically distant compared to those of CC1 and CC3, and belonged to L. pneumophila subsp. fraseri.

Comparative Analysis of Tools and Approaches for Source Tracking Listeria monocytogenes in a Food Facility Using Whole-Genome Sequence Data

As WGS is increasingly used by food industry to characterize pathogen isolates, users are challenged by the variety of analysis approaches available, ranging from methods that require extensive bioinformatics expertise to commercial software packages. This study aimed to assess the impact of analysis pipelines (i.e., different hqSNP pipelines, a cg/wgMLST pipeline) and the reference genome selection on analysis results (i.e., hqSNP and allelic differences as well as tree topologies) and conclusion drawn. For these comparisons, whole genome sequences were obtained for 40 Listeria monocytogenes isolates collected over 18 years from a cold-smoked salmon facility and 2 other isolates obtained from different facilities as part of academic research activities; WGS data were analyzed with three hqSNP pipelines and two MLST pipelines. After initial clustering using a k-mer based approach, hqSNP pipelines were run using two types of reference genomes: (i) closely related closed genomes (“closed references”) and (ii) high-quality de novo assemblies of the dataset isolates (“draft references”). All hqSNP pipelines identified similar hqSNP difference ranges among isolates in a given cluster; use of different reference genomes showed minimal impacts on hqSNP differences identified between isolate pairs. Allelic differences obtained by wgMLST showed similar ranges as hqSNP differences among isolates in a given cluster; cgMLST consistently showed fewer differences than wgMLST. However, phylogenetic trees and dendrograms, obtained based on hqSNP and cg/wgMLST data, did show some incongruences, typically linked to clades supported by low bootstrap values in the trees. When a hqSNP cutoff was used to classify isolates as “related” or “unrelated,” use of different pipelines yielded a considerable number of discordances; this finding supports that cut-off values are valuable to provide a starting point for an investigation, but supporting and epidemiological evidence should be used to interpret WGS data. Overall, our data suggest that cgMLST-based data analyses provide for appropriate subtype differentiation and can be used without the need for preliminary data analyses (e.g., k-mer based clustering) or external closed reference genomes, simplifying data analyses needs. hqSNP or wgMLST analyses can be performed on the isolate clusters identified by cgMLST to increase the precision on determining the genomic similarity between isolates.

A Pathovar of Xanthomonas oryzae Infecting Wild Grasses Provides Insight Into the Evolution of Pathogenicity in Rice Agroecosystems

Xanthomonas oryzae (Xo) are globally important rice pathogens. Virulent lineages from Africa and Asia and less virulent strains from the United States have been well characterized. Xanthomonas campestris pv. leersiae (Xcl), first described in 1957, causes bacterial streak on the perennial grass, Leersia hexandra, and is a close relative of Xo. L. hexandra, a member of the Poaceae, is highly similar to rice phylogenetically, is globally ubiquitous around rice paddies, and is a reservoir of pathogenic Xo. We used long read, single molecule real time (SMRT) genome sequences of five strains of Xcl from Burkina Faso, China, Mali, and Uganda to determine the genetic relatedness of this organism with Xo. Novel transcription activator-like effectors (TALEs) were discovered in all five strains of Xcl. Predicted TALE target sequences were identified in the Leersia perrieri genome and compared to rice susceptibility gene homologs. Pathogenicity screening on L. hexandra and diverse rice cultivars confirmed that Xcl are able to colonize rice and produce weak but not progressive symptoms. Overall, based on average nucleotide identity (ANI), type III (T3) effector repertoires, and disease phenotype, we propose to rename Xcl to X. oryzae pv. leersiae (Xol) and use this parallel system to improve understanding of the evolution of bacterial pathogenicity in rice agroecosystems.

Genes Associated With Psychrotolerant Bacillus cereus Group Isolates

The Bacillus cereus group comprises 18 different species, including human pathogens as well as psychrotolerant strains that are an important cause of fluid milk spoilage. To enhance our understanding of the genetic markers associated with psychrotolerance (defined here as > 1 log₁₀ increase in cfu/mL after 21 days incubation at 6°C) among dairy-associated B. cereus group isolates, we used genetic (whole genome sequencing) and phenotypic methods [growth in Skim Milk Broth (SMB) and Brain Heart Infusion (BHI) broth] to characterize 23 genetically-distinct representative isolates from a collection of 503 dairy-associated isolates. Quality threshold clustering identified three categories of psychrotolerance: (i) 14 isolates that were not psychrotolerant in BHI or SMB, (ii) 6 isolates that were psychrotolerant in BHI but not in SMB, and (iii) 2 isolates that were psychrotolerant in BHI and SMB. One isolate, which was psychrotolerant in BHI broth but was just below the cut-off of >1 log₁₀ cfu/mL increase in SMB was not assigned to a cluster. A maximum likelihood phylogeny constructed with core genome single nucleotide polymorphisms classified all psychrotolerant isolates (i.e., psychrotolerant in BHI) into clade VI (representing B. mycoides/weihenstephanensis). Analysis of correlations between gene ortholog presence or absence patterns and psychrotolerance identified 206 orthologous gene clusters that were significantly overrepresented among psychrotolerant strains, including two clusters of cold shock proteins, which were identified in 8/9 and 7/9 psychrotolerant isolates. Gene ontology analyses revealed 36 gene ontology terms that were overrepresented in psychrotolerant isolates, including putrescine catabolic processes and putrescine transmembrane transporter activity. Lastly, Hidden Markov Model searches identified three protein family motifs, including cold shock domain proteins and fatty acid hydroxylases that were significantly associated with psychrotolerance in BHI broth. Analyses of CspA sequences revealed a positive association between psychrotolerant strains and a previously identified “psychrotolerant” CspA sequence. Overall, our data highlight genetic and phenotypic differences in psychrotolerance among B. cereus group dairy-associated isolates and show that psychrotolerance is dependent on the growth medium. We also identified a number of gene targets that could be used for specific detection or control of psychrotolerant B. cereus group isolates.

Long-term Colonization by Campylobacter jejuni Within a Human Host: Evolution, Antimicrobial Resistance, and Adaptation

Background

Campylobacteriosis is inflammation of the gastrointestinal tract as a result of Campylobacter infection. Most campylobacteriosis cases are acute and self-limiting, with Campylobacter excretion ceasing a few weeks after symptoms cease. We identified a patient with fecal specimens positive for Campylobacter jejuni (ST45) intermittently during a 10-year period.

Methods

Sixteen Campylobacter isolates were collected from the patient during 2006-2016. The isolates' genomes were sequenced to determine their relatedness, and their antimicrobial susceptibility patterns and motility were measured to determine the effects of antibiotic therapy and long-term excretion on the Campylobacter population.

Results

Phylogenetic analyses estimated that the isolates shared a date of common ancestor between 1998 and 2006, coinciding with the onset of symptoms for the patient. Genomic analysis identified selection for changes in motility, and antimicrobial susceptibility testing suggested that the Campylobacter population developed resistance to several antibiotics coinciding with periods of antibiotic therapy.

Conclusions

The patient was consistently colonized with organisms from a Campylobacter population that adapted to the internal environment of the patient. Genomic and phylogenetic analyses can give insight into a patient's infection history and the effect of antimicrobial treatment on Campylobacter populations in this unusual situation of long-term colonization of an individual.

Evolutionary Model of Cluster Divergence of the Emergent Marine Pathogen Vibrio vulnificus: From Genotype to Ecotype

Vibrio vulnificus is an emergent marine pathogen and is the cause of a deadly septicemia. However, the genetic factors that differentiate its clinical and environmental strains and its several biotypes remain mostly enigmatic. In this work, we investigated the underlying genomic properties and population dynamics of the V. vulnificus species to elucidate the traits that make these strains emerge as a human pathogen. The acquisition of different ecological determinants could have allowed the development of highly divergent clusters with different lifestyles within the same environment. However, we identified strains from both clusters in the mucosa of aquaculture species, indicating that manmade niches are bringing strains from the two clusters together, posing a potential risk of recombination and of emergence of novel variants. We propose a new evolutionary model that provides a perspective that could be broadly applicable to other pathogenic vibrios and facultative bacterial pathogens to pursue strategies to prevent their infections.

Vibrio vulnificus, an opportunistic pathogen, is the causative agent of a life-threatening septicemia and a rising problem for aquaculture worldwide. The genetic factors that differentiate its clinical and environmental strains remain enigmatic. Furthermore, clinical strains have emerged from every clade of V. vulnificus. In this work, we investigated the underlying genomic properties and population dynamics of the V. vulnificus species from an evolutionary and ecological point of view. Genome comparisons and bioinformatic analyses of 113 V. vulnificus isolates indicate that the population of V. vulnificus is made up of four different clusters. We found evidence that recombination and gene flow between the two largest clusters (cluster 1 [C1] and C2) have drastically decreased to the point where they are diverging independently. Pangenome and phenotypic analyses showed two markedly different lifestyles for these two clusters, indicating commensal (C2) and bloomer (C1) ecotypes, with differences in carbohydrate utilization, defense systems, and chemotaxis, among other characteristics. Nonetheless, we identified frequent intra- and interspecies exchange of mobile genetic elements (e.g., antibiotic resistance plasmids, novel “chromids,” or two different and concurrent type VI secretion systems) that provide high levels of genetic diversity in the population. Surprisingly, we identified strains from both clusters in the mucosa of aquaculture species, indicating that manmade niches are bringing strains from the two clusters together. We propose an evolutionary model of V. vulnificus that could be broadly applicable to other pathogenic vibrios and facultative bacterial pathogens to pursue strategies to prevent their infections and emergence.

Characterization of Emetic and Diarrheal Bacillus cereus Strains From a 2016 Foodborne Outbreak Using Whole-Genome Sequencing: Addressing the Microbiological, Epidemiological, and Bioinformatic Challenges

The Bacillus cereus group comprises multiple species capable of causing emetic or diarrheal foodborne illness. Despite being responsible for tens of thousands of illnesses each year in the U.S. alone, whole-genome sequencing (WGS) is not yet routinely employed to characterize B. cereus group isolates from foodborne outbreaks. Here, we describe the first WGS-based characterization of isolates linked to an outbreak caused by members of the B. cereus group. In conjunction with a 2016 outbreak traced to a supplier of refried beans served by a fast food restaurant chain in upstate New York, a total of 33 B. cereus group isolates were obtained from human cases (n = 7) and food samples (n = 26). Emetic (n = 30) and diarrheal (n = 3) isolates were most closely related to B. paranthracis (group III) and B. cereus sensu stricto (group IV), respectively. WGS indicated that the 30 emetic isolates (24 and 6 from food and humans, respectively) were closely related and formed a well-supported clade distinct from publicly available emetic group III genomes with an identical sequence type (ST 26). The 30 emetic group III isolates from this outbreak differed from each other by a mean of 8.3 to 11.9 core single nucleotide polymorphisms (SNPs), while differing from publicly available emetic group III ST 26 B. cereus group genomes by a mean of 301.7-528.0 core SNPs, depending on the SNP calling methodology used. Using a WST-1 cell proliferation assay, the strains isolated from this outbreak had only mild detrimental effects on HeLa cell metabolic activity compared to reference diarrheal strain B. cereus ATCC 14579. We hypothesize that the outbreak was a single source outbreak caused by emetic group III B. cereus belonging to the B. paranthracis species, although food samples were not tested for presence of the emetic toxin cereulide. In addition to showcasing how WGS can be used to characterize B. cereus group strains linked to a foodborne outbreak, we also discuss potential microbiological and epidemiological challenges presented by B. cereus group outbreaks, and we offer recommendations for analyzing WGS data from the isolates associated with them.

Risk Assessment After a Severe Hospital-Acquired Infection Associated With Carbapenemase-Producing Pseudomonas aeruginosa

IMPORTANCE

Resistance of gram-negative bacilli to carbapenems is rapidly emerging worldwide. In 2016, the World Health Organization defined the hospital-built environment as a core component of infection prevention and control programs. The hospital-built environment has recently been reported as a source for outbreaks and sporadic transmission events of carbapenemase-producing gram-negative bacilli from the environment to patients.

OBJECTIVE

To assess risk after the identification of an unexpected, severe, and lethal hospital-acquired infection caused by carbapenemase-producing Pseudomonas aeruginosa in a carbapenemase-low endemic setting.

DESIGN, SETTINGS, AND PARTICIPANTS

A case series study in which a risk assessment was performed on all 11 patients admitted to the combined cardiothoracic surgery and pulmonary diseases ward and the hospital-built environment in the Radboud University Medical Center, the Netherlands, in February 2018.

EXPOSURES

Water and aerosols containing carbapenemase-producing (Verona integron-mediated metallo-β-lactamase [VIM]) P aeruginosa.

MAIN OUTCOMES AND MEASURES

Colonization and/or infection of patients and/or contamination of the environment after the detection of 1 patient infected with carbapenemase-producing (VIM) P aeruginosa.

RESULTS

A total of 5 men (age range, 60-84 years) and 6 women (age range, 55-74 years) were admitted to the combined cardiothoracic surgery and pulmonary diseases ward. The risk assessment was performed after carbapenemase-producing (VIM) P aeruginosa was unexpectedly detected in a man in his early 60s, who had undergone a left-sided pneumonectomy and adjuvant radiotherapy. No additional cases (colonization or infection) of carbapenemase-producing (VIM) P aeruginosa were detected. Plausible transmission of carbapenemase-producing P aeruginosa from the hospital environment to the patient via the air was confirmed by whole-genome sequencing, which proved the relation of Pseudomonas strains from the patient, the shower drains in 8 patient rooms, 1 sink, and an air sample.

CONCLUSIONS AND RELEVANCE

This study suggests that rethinking the hospital-built environment, including shower drains and the sewage system, will be crucial for the prevention of severe and potential lethal hospital-acquired infections.

Phenotypic Prediction: Linking in vitro Virulence to the Genomics of 59 Salmonella enterica Strains

The increased availability of whole-genome-sequencing techniques generates a wealth of DNA data on numerous organisms, including foodborne pathogens such as Salmonella. However, how these data can be used to improve microbial risk assessment and understanding of Salmonella epidemiology remains a challenge. The aim of this study was to assess variability in in vitro virulence and genetic characteristics between and within different serovars. The phenotypic behavior of 59 strains of 32 different Salmonella enterica serovars from animal, human and food origin was assessed in an in vitro gastro-intestinal tract (GIT) system and they were analyzed for the presence of 233 putative virulence genes as markers for phenotypic prediction. The probability of in vitro infection, P(inf), defined as the fraction of infectious cells passing from inoculation to host cell invasion at the last stage of the GIT system, was interpreted as the in vitro virulence. Results showed that the (average) P(inf) of Salmonella serovars ranged from 5.3E-05 (S. Kedougou) to 5.2E-01 (S. Typhimurium). In general, a higher P(inf) on serovar level corresponded to higher reported human incidence from epidemiological reporting data. Of the 233 virulence genes investigated, only 101 showed variability in presence/absence among the strains. In vitro P(inf) was found to be positively associated with the presence of specific plasmid related virulence genes (mig-5, pef, rck, and spv). However, not all serovars with a relatively high P(inf), > 1E-02, could be linked with these specific genes. Moreover, some outbreak related strains (S. Heidelberg and S. Thompson) did not reveal this association with P(inf). No clear association with in vitro virulence P(inf) was identified when grouping serovars with the same virulence gene profile (virulence plasmid, Typhoid toxin, peg operon and stk operon). This study shows that the in vitro P(inf) variation among individual strains from the same serovar is larger than that found between serovars. Therefore, ranking P(inf) of S. enterica on serovar level alone, or in combination with a serovar specific virulence gene profile, cannot be recommended. The attribution of single biological phenomena to individual strains or serovars is not sufficient to improve the hazard characterization for S. enterica. Future microbial risk assessments, including virulence gene profiles, require a systematic approach linked to epidemiological studies rather than revealing differences in characteristics on serovar level alone.

Streptococcus infantis, Streptococcus mitis, and Streptococcus oralis Strains With Highly Similar cps5 Loci and Antigenic Relatedness to Serotype 5 Pneumococci

Streptococcus pneumoniae is a highly impactful bacterial pathogen on a global scale. The principal pneumococcal virulence factor and target of effective vaccines is its polysaccharide capsule, of which there are many structurally distinct forms. Here, we describe four distinct strains of three Mitis group commensal species (Streptococcus infantis, Streptococcus mitis, and Streptococcus oralis) recovered from upper respiratory tract specimens from adults in Kenya and the United States that were PCR-positive for the pneumococcal serotype 5 specific gene, wzy5. For each of the four strains, the 15 genes comprising the capsular polysaccharide biosynthetic gene cluster (cps5) shared the same order found in serotype 5 pneumococci, and each of the serotype 5-specific genes from the serotype 5 pneumococcal reference strain shared 76-99% sequence identity with the non-pneumococcal counterparts. Double-diffusion experiments demonstrated specific reactivity of the non-pneumococcal strains with pneumococcal serotype 5 typing sera. Antiserum raised against S. mitis strain KE67013 specifically reacted with serotype 5 pneumococci for a positive Quellung reaction and stimulated serotype 5 specific opsonophagocytic killing of pneumococci. Four additional commensal strains, identified using PCR serotyping assays on pharyngeal specimens, revealed loci highly homologous to those of pneumococci of serotypes 12F, 15A, 18C, and 33F. These data, in particular the species and strain diversity shown for serotype 5, highlight the existence of a broad non-pneumococcal species reservoir in the upper respiratory tract for the expression of capsular polysaccharides that are structurally related or identical to those corresponding to epidemiologically significant serotypes. Very little is known about the genetic and antigenic capsular diversity among the vast array of commensal streptococcal strains that represent multiple diverse species. The discovery of serotype 5 strains within three different commensal species suggests that extensive capsular serologic overlap exists between pneumococci and other members of the diverse Mitis group. These findings may have implications for our current understanding of naturally acquired immunity to S. pneumoniae and pneumococcal serotype distributions in different global regions. Further characterization of commensal strains carrying homologs of serotype-specific genes previously thought to be specific for pneumococci of known serotypes may shed light on the evolution of these important loci.

Platforms and Analytical Tools Used in Nucleic Acid Sequence-Based Microbial Genotyping Procedures

In the decade and a half since the introduction of next-generation sequencing (NGS), the technical feasibility, cost, and overall utility of sequencing have changed dramatically, including applications for infectious disease epidemiology. Massively parallel sequencing technologies have decreased the cost of sequencing by more than 6 orders or magnitude over this time, with a corresponding increase in data generation and complexity. This review provides an overview of the basic principles, chemistry, and operational mechanics of current sequencing technologies, including both conventional Sanger and NGS approaches. As the generation of large amounts of sequence data becomes increasingly routine, the role of bioinformatics in data analysis and reporting becomes all the more critical, and the successful deployment of NGS in public health settings requires careful consideration of changing information technology, bioinformatics, workforce, and regulatory requirements. While there remain important challenges to the sustainable implementation of NGS in public health, in terms of both laboratory and bioinformatics capacity, the impact of these technologies on infectious disease surveillance and outbreak investigations has been nothing short of revolutionary. Understanding the important role that NGS plays in modern public health laboratory practice is critical, as is the need to ensure appropriate workforce, infrastructure, facilities, and funding consideration for routine NGS applications, future innovation, and rapidly scaling NGS-based infectious disease surveillance and outbreak response activities. *This article is part of a curated collection.

Draft genomes of Cronobacter sakazakii strains isolated from dried spices bring unique insights into the diversity of plant-associated strains

Cronobacter sakazakii is a Gram-negative opportunistic pathogen that causes life- threatening infantile infections, such as meningitis, septicemia, and necrotizing enterocolitis, as well as pneumonia, septicemia, and urinary tract and wound infections in adults. Here, we report 26 draft genome sequences of C. sakazakii, which were obtained from dried spices from the USA, the Middle East, China, and the Republic of Korea. The average genome size of the C. sakazakii genomes was 4393 kb, with an average of 4055 protein coding genes, and an average genome G + C content of 56.9%. The genomes contained genes related to carbohydrate transport and metabolism, amino acid transport and metabolism, and cell wall/membrane biogenesis. In addition, we identified genes encoding proteins involved in osmotic responses such as DnaJ, Aquaproin Z, ProQ, and TreF, as well as virulence-related and heat shock-related proteins. Interestingly, a metabolic island comprised of a variably-sized xylose utilization operon was found within the spice-associated C. sakazakii genomes, which supports the hypothesis that plants may serve as transmission vectors or alternative hosts for Cronobacter species. The presence of the genes identified in this study can support the remarkable phenotypic traits of C. sakazakii such as the organism's capabilities of adaptation and survival in response to adverse growth environmental conditions (e.g. osmotic and desiccative stresses). Accordingly, the genome analyses provided insights into many aspects of physiology and evolutionary history of this important foodborne pathogen.

Enterobacter cloacae Complex Sequence Type 171 Isolates Expressing KPC-4 Carbapenemase Recovered from Canine Patients in Ohio

Companion animals are likely relevant in the transmission of antimicrobial-resistant bacteria. Enterobacter xiangfangensis sequence type 171 (ST171), a clone that has been implicated in clusters of infections in humans, was isolated from two dogs with clinical disease in Ohio. The canine isolates contained IncHI2 plasmids encoding bla _KPC-4 Whole-genome sequencing was used to put the canine isolates in phylogenetic context with available human ST171 sequences, as well as to characterize their bla _KPC-4 plasmids.

Complete genome sequencing of three human clinical isolates of Staphylococcus caprae reveals virulence factors similar to those of S. epidermidis and S. capitis

Background

Staphylococcus caprae is an animal-associated bacterium regarded as part of goats’ microflora. Recently, S. caprae has been reported to cause human nosocomial infections such as bacteremia and bone and joint infections. However, the mechanisms responsible for the development of nosocomial infections remain largely unknown. Moreover, the complete genome sequence of S. caprae has not been determined.

Results

We determined the complete genome sequences of three methicillin-resistant S. caprae strains isolated from humans and compared these sequences with the genomes of S. epidermidis and S. capitis, both of which are closely related to S. caprae and are inhabitants of human skin capable of causing opportunistic infections. The genomes showed that S. caprae JMUB145, JMUB590, and JMUB898 strains contained circular chromosomes of 2,618,380, 2,629,173, and 2,598,513 bp, respectively. JMUB145 carried type V SCCmec, while JMUB590 and JMUB898 had type IVa SCCmec. A genome-wide phylogenetic SNP tree constructed using 83 complete genome sequences of 24 Staphylococcus species and 2 S. caprae draft genome sequences confirmed that S. caprae is most closely related to S. epidermidis and S. capitis. Comparative complete genome analysis of eight S. epidermidis, three S. capitis and three S. caprae strains revealed that they shared similar virulence factors represented by biofilm formation genes. These factors include wall teichoic acid synthesis genes, poly-gamma-DL-glutamic acid capsule synthesis genes, and other genes encoding nonproteinaceous adhesins. The 17 proteinases/adhesins and extracellular proteins known to be associated with biofilm formation in S. epidermidis were also conserved in these three species, and their biofilm formation could be detected in vitro. Moreover, two virulence-associated gene clusters, the type VII secretion system and capsular polysaccharide biosynthesis gene clusters, identified in S. aureus were present in S. caprae but not in S. epidermidis and S. capitis genomes.

Conclusion

The complete genome sequences of three methicillin-resistant S. caprae isolates from humans were determined for the first time. Comparative genome analysis revealed that S. caprae is closely related to S. epidermidis and S. capitis at the species level, especially in the ability to form biofilms, which may lead to increased virulence during the development of S. caprae infections.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5185-9) contains supplementary material, which is available to authorized users.

Analysis of enterotoxigenic Bacillus cereus strains from dried foods using whole genome sequencing, multi-locus sequence analysis and toxin gene prevalence and distribution using endpoint PCR analysis

Bacillus cereus strains were isolated from dried foods, which included international brands of spices from South East Asia, Mexico and India purchased from several retail stores, samples of powdered infant formula (PIF), medicated fish feed and dietary supplements. The genetic diversity of 64 strains from spices and PIF was determined using a multiplex endpoint PCR assay designed to identify hemolysin BL, nonhemolytic enterotoxin, cytotoxin K, and enterotoxin FM toxin genes. Thirteen different B. cereus toxigenic gene patterns or profiles were identified among the strains. Randomly selected B. cereus strains were sequenced and compared with reference Genomic Groups from National Center Biotechnology Information using bioinformatics tools. A comprehensive multi-loci sequence analysis (MLSA) was designed using alleles from 25 known MLST genes specifically tailored for use with whole genome assemblies. A cohort of representative genomes of strains from a few FDA regulated commodities like dry foods and medicated fish feed was used to demonstrate the utility of the 25-MLSA approach for rapid clustering and identification of Genome Groups. The analysis clustered the strains from medicated fish feed, dry foods, and dietary supplements into phylogenetically-related groups. 25-MLSA also pointed to a greater diversity of B. cereus strains from foods and feed than previously recognized. Our integrated approach of toxin gene PCR, and to our knowledge, whole genome sequencing (WGS) based sequence analysis, may be the first of its kind that demonstrates enterotoxigenic potential and genomic diversity in parallel.

Safety of the fungal workhorses of industrial biotechnology: update on the mycotoxin and secondary metabolite potential of Aspergillus niger, Aspergillus oryzae, and Trichoderma reesei

This review presents an update on the current knowledge of the secondary metabolite potential of the major fungal species used in industrial biotechnology, i.e., Aspergillus niger, Aspergillus oryzae, and Trichoderma reesei. These species have a long history of safe use for enzyme production. Like most microorganisms that exist in a challenging environment in nature, these fungi can produce a large variety and number of secondary metabolites. Many of these compounds present several properties that make them attractive for different industrial and medical applications. A description of all known secondary metabolites produced by these species is presented here. Mycotoxins are a very limited group of secondary metabolites that can be produced by fungi and that pose health hazards in humans and other vertebrates when ingested in small amounts. Some mycotoxins are species-specific. Here, we present scientific basis for (1) the definition of mycotoxins including an update on their toxicity and (2) the clarity on misclassification of species and their mycotoxin potential reported in literature, e.g., A. oryzae has been wrongly reported as an aflatoxin producer, due to misclassification of Aspergillus flavus strains. It is therefore of paramount importance to accurately describe the mycotoxins that can potentially be produced by a fungal species that is to be used as a production organism and to ensure that production strains are not capable of producing mycotoxins during enzyme production. This review is intended as a reference paper for authorities, companies, and researchers dealing with secondary metabolite assessment, risk evaluation for food or feed enzyme production, or considerations on the use of these species as production hosts.

Whole genome sequencing for investigations of meningococcal outbreaks in the United States: a retrospective analysis

Although rare in the U.S., outbreaks due to Neisseria meningitidis do occur. Rapid, early outbreak detection is important for timely public health response. In this study, we characterized U.S. meningococcal isolates (N = 201) from 15 epidemiologically defined outbreaks (2009-2015) along with temporally and geographically matched sporadic isolates using multilocus sequence typing, pulsed-field gel electrophoresis (PFGE), and six whole genome sequencing (WGS) based methods. Recombination-corrected maximum likelihood (ML) and Bayesian phylogenies were reconstructed to identify genetically related outbreak isolates. All WGS analysis methods showed high degree of agreement and distinguished isolates with similar or indistinguishable PFGE patterns, or the same strain genotype. Ten outbreaks were caused by a single strain; 5 were due to multiple strains. Five sporadic isolates were phylogenetically related to 2 outbreaks. Analysis of 9 outbreaks using timed phylogenies identified the possible origin and estimated the approximate time that the most recent common ancestor emerged for outbreaks analyzed. U.S. meningococcal outbreaks were caused by single- or multiple-strain introduction, with organizational outbreaks mainly caused by a clonal strain and community outbreaks by divergent strains. WGS can infer linkage of meningococcal cases when epidemiological links are uncertain. Accurate identification of outbreak-associated cases requires both WGS typing and epidemiological data.

Genomic heterogeneity differentiates clinical and environmental subgroups of Legionella pneumophila sequence type 1

Legionella spp. are the cause of a severe bacterial pneumonia known as Legionnaires' disease (LD). In some cases, current genetic subtyping methods cannot resolve LD outbreaks caused by common, potentially endemic L. pneumophila (Lp) sequence types (ST), which complicates laboratory investigations and environmental source attribution. In the United States (US), ST1 is the most prevalent clinical and environmental Lp sequence type. In order to characterize the ST1 population, we sequenced 289 outbreak and non-outbreak associated clinical and environmental ST1 and ST1-variant Lp strains from the US and, together with international isolate sequences, explored their genetic and geographic diversity. The ST1 population was highly conserved at the nucleotide level; 98% of core nucleotide positions were invariant and environmental isolates unassociated with human disease (n = 99) contained ~65% more nucleotide diversity compared to clinical-sporadic (n = 139) or outbreak-associated (n = 28) ST1 subgroups. The accessory pangenome of environmental isolates was also ~30-60% larger than other subgroups and was enriched for transposition and conjugative transfer-associated elements. Up to ~10% of US ST1 genetic variation could be explained by geographic origin, but considerable genetic conservation existed among strains isolated from geographically distant states and from different decades. These findings provide new insight into the ST1 population structure and establish a foundation for interpreting genetic relationships among ST1 strains; these data may also inform future analyses for improved outbreak investigations.

Coli Surface Antigen 26 Acts as an Adherence Determinant of Enterotoxigenic Escherichia coli and Is Cross-Recognized by Anti-CS20 Antibodies

The coli surface antigen 26 (CS26) of enterotoxigenic Escherichia coli (ETEC) had been described as a putative adhesive pilus based on the partial sequence of the crsH gene, detected in isolates from children with diarrhea in Egypt. However, its production and activity as adherence determinant has not been experimentally addressed. The crsH was identified as a homolog of genes encoding structural subunits of ETEC colonization factors (CFs) CS12, CS18, and CS20. These CFs, along with the recently discovered CS30, belong to the γ₂ family of pili assembled by the chaperone-usher pathway (CU pili). Further, the complete CS26 locus, crsHBCDEFG, was described in an O141 ETEC strain (ETEC 100664) obtained from a diarrhea case in The Gambia, during the Global Enterics Multicenter Study. Here, we report that CS26 is a pilus of ∼10 nm in diameter, with the capacity to increase the cell adherence of the non-pathogenic strain E. coli DH10B. As for other related pili, production of CS26 seems to be regulated by phase variation. Deletion of crsHBCDEFG in ETEC 100664 significantly decreased its adherence capacity, which was recovered by in trans complementation. Furthermore, CrsH was cross-recognized by polyclonal antibodies directed against the major structural subunit of CS20, CsnA, as determined by Western blotting and immunogold labeling. ETEC CS26+ strains were found to harbor the heat-labile enterotoxin only, within three different sequence types of phylogroups A and B1, the latter suggesting acquisition through independent events of horizontal transfer. Overall, our results demonstrate that CS26 is an adhesive pilus of human ETEC. In addition, cross-reactivity with anti-CsnA antibodies indicate presence of common epitopes in γ₂-CFs.

Comparative genomics of Czech vaccine strains of Bordetella pertussis

Bordetella pertussis is a strictly human pathogen causing the respiratory infectious disease called whooping cough or pertussis. B. pertussis adaptation to acellular pertussis vaccine pressure has been repeatedly highlighted, but recent data indicate that adaptation of circulating strains started already in the era of the whole cell pertussis vaccine (wP) use. We sequenced the genomes of five B. pertussis wP vaccine strains isolated in the former Czechoslovakia in the pre-wP (1954-1957) and early wP (1958-1965) eras, when only limited population travel into and out of the country was possible. Four isolates exhibit a similar genome organization and form a distinct phylogenetic cluster with a geographic signature. The fifth strain is rather distinct, both in genome organization and SNP-based phylogeny. Surprisingly, despite isolation of this strain before 1966, its closest sequenced relative appears to be a recent isolate from the US. On the genome content level, the five vaccine strains contained both new and already described regions of difference. One of the new regions contains duplicated genes potentially associated with transport across the membrane. The prevalence of this region in recent isolates indicates that its spread might be associated with selective advantage leading to increased strain fitness.

Metagenomic analysis with strain-level resolution reveals fine-scale variation in the human pregnancy microbiome

Recent studies suggest that the microbiome has an impact on gestational health and outcome. However, characterization of the pregnancy-associated microbiome has largely relied on 16S rRNA gene amplicon-based surveys. Here, we describe an assembly-driven, metagenomics-based, longitudinal study of the vaginal, gut, and oral microbiomes in 292 samples from 10 subjects sampled every three weeks throughout pregnancy. Nonhuman sequences in the amount of 1.53 Gb were assembled into scaffolds, and functional genes were predicted for gene- and pathway-based analyses. Vaginal assemblies were binned into 97 draft quality genomes. Redundancy analysis (RDA) of microbial community composition at all three body sites revealed gestational age to be a significant source of variation in patterns of gene abundance. In addition, health complications were associated with variation in community functional gene composition in the mouth and gut. The diversity of Lactobacillus iners-dominated communities in the vagina, unlike most other vaginal community types, significantly increased with gestational age. The genomes of co-occurring Gardnerella vaginalis strains with predicted distinct functions were recovered in samples from two subjects. In seven subjects, gut samples contained strains of the same Lactobacillus species that dominated the vaginal community of that same subject and not other Lactobacillus species; however, these within-host strains were divergent. CRISPR spacer analysis suggested shared phage and plasmid populations across body sites and individuals. This work underscores the dynamic behavior of the microbiome during pregnancy and suggests the potential importance of understanding the sources of this behavior for fetal development and gestational outcome.

Whole-Genome-Sequencing characterization of bloodstream infection-causing hypervirulent Klebsiella pneumoniae of capsular serotype K2 and ST374

Hypervirulent K. pneumoniae variants (hvKP) have been increasingly reported worldwide, causing metastasis of severe infections such as liver abscesses and bacteremia. The capsular serotype K2 hvKP strains show diverse multi-locus sequence types (MLSTs), but with limited genetics and virulence information. In this study, we report a hypermucoviscous K. pneumoniae strain, RJF293, isolated from a human bloodstream sample in a Chinese hospital. It caused a metastatic infection and fatal septic shock in a critical patient. The microbiological features and genetic background were investigated with multiple approaches. The Strain RJF293 was determined to be multilocis sequence type (ST) 374 and serotype K2, displayed a median lethal dose (LD50) of 1.5 × 10² CFU in BALB/c mice and was as virulent as the ST23 K1 serotype hvKP strain NTUH-K2044 in a mouse lethality assay. Whole genome sequencing revealed that the RJF293 genome codes for 32 putative virulence factors and exhibits a unique presence/absence pattern in comparison to the other 105 completely sequenced K. pneumoniae genomes. Whole genome SNP-based phylogenetic analysis revealed that strain RJF293 formed a single clade, distant from those containing either ST66 or ST86 hvKP. Compared to the other sequenced hvKP chromosomes, RJF293 contains several strain-variable regions, including one prophage, one ICEKp1 family integrative and conjugative element and six large genomic islands. The sequencing of the first complete genome of an ST374 K2 hvKP clinical strain should reinforce our understanding of the epidemiology and virulence mechanisms of this bloodstream infection-causing hvKP with clinical significance.

Microbial sequence typing in the genomic era

Next-generation sequencing (NGS), also known as high-throughput sequencing, is changing the field of microbial genomics research. NGS allows for a more comprehensive analysis of the diversity, structure and composition of microbial genes and genomes compared to the traditional automated Sanger capillary sequencing at a lower cost. NGS strategies have expanded the versatility of standard and widely used typing approaches based on nucleotide variation in several hundred DNA sequences and a few gene fragments (MLST, MLVA, rMLST and cgMLST). NGS can now accommodate variation in thousands or millions of sequences from selected amplicons to full genomes (WGS, NGMLST and HiMLST). To extract signals from high-dimensional NGS data and make valid statistical inferences, novel analytic and statistical techniques are needed. In this review, we describe standard and new approaches for microbial sequence typing at gene and genome levels and guidelines for subsequent analysis, including methods and computational frameworks. We also present several applications of these approaches to some disciplines, namely genotyping, phylogenetics and molecular epidemiology.

In vivo evolution of drug-resistant Mycobacterium tuberculosis in patients during long-term treatment

Background

In the current scenario, the drug-resistant tuberculosis is a significant challenge in the control of tuberculosis worldwide. In order to investigate the in vivo evolution of drug-resistant M. tuberculosis, the present study envisaged sequencing of the draft genomes of 18 serial isolates from four pre-extensively drug-resistant (pre-XDR) tuberculosis patients for continuous genetic alterations.

Results

All of the isolates harbored single nucleotide polymorphisms (SNPs) ranging from 1303 to 1309 with M. tuberculosis H37Rv as the reference. SNPs ranged from 0 to 12 within patients. The evolution rates were higher than the reported SNPs of 0.5 in the four patients. All the isolates exhibited mutations at sites of known drug targets, while some contained mutations in uncertain drug targets including folC, proZ, and pyrG. The compensatory substitutions for rescuing these deleterious mutations during evolution were only found in RpoC I491T in one patient. Many loci with microheterogeneity showed transient mutations in different isolates. Ninety three SNPs exhibited significant association with refractory pre-XDR TB isolates.

Conclusions

Our results showed evolutionary changes in the serial genetic characteristics of the pre-XDR TB patients due to accumulation of the fixed drug-resistant related mutations, and the transient mutations under continuous antibiotics pressure over several years.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5010-5) contains supplementary material, which is available to authorized users.

Whole genome shotgun sequencing revealed highly polymorphic genome regions and genes in Escherichia coli O157:H7 isolates collected from a single feedlot

Escherichia coli serotype O157:H7 continues to pose a serious health threat to human beings. Cattle, a major reservoir of the pathogen, harbor E. coli O157:H7 in their gastrointestinal tract and shed variable concentrations of E. coli O157:H7 into the environment. Genetic characterization of cattle-shed E. coli O157 strains is of interest to the livestock industry, food business, and public health community. The present study applied whole genome shotgun sequencing (WGS) and single nucleotide variant (SNV) calling to characterize 279 cattle-shed E. coli O157:H7 strains isolated from a single feedlot located in southwestern region of the US. More than 4,000 SNVs were identified among the strains and the resultant phylogenomic tree revealed three major groups. Using the Sakai strain genome as reference, more than 2,000 SNVs were annotated and a detailed SNV map generated. Results clearly revealed highly polymorphic loci along the E. coli O157:H7 genome that aligned with the prophage regions and highly variant genes involved in processing bacterial genetic information. The WGS data were further profiled against a comprehensive virulence factor database (VFDB) for virulence gene identification. Among the total 285 virulence genes identified, only 132 were present in all the strains. There were six virulence genes unique to single isolates. Our findings suggested that the genome variations of the E. coli O157:H7 were mainly attributable to dynamics of certain phages, and the bacterial strains have variable virulence gene profiles, even though they came from a single cattle population, which may explain the differences in pathogenicity, host prevalence, and transmissibility by E. coli O157:H7.

Acquisition of antimicrobial-resistant variants in repeated infections caused by Pseudomonas aeruginosa revealed by whole genome sequencing

Pseudomonas aeruginosa, responsible for serious nosocomial-acquired infections, possesses intrinsic antibiotic resistance mechanisms and commonly exhibits multidrug resistance. Here, we report the evolving resistance profiles of strains isolated from the sputum of a patient being treated for repeated P. aeruginosa infections following cancer resection. Whole genome sequencing of six isolates obtained over a 2-month period revealed two key single nucleotide polymorphisms in the mexR and gyrB genes that affected efflux pump expression and antimicrobial resistance.

Genome-driven evaluation and redesign of PCR tools for improving the detection of virulence-associated genes in aeromonads

Many virulence factors have been described for opportunistic pathogens within the genus Aeromonas. Polymerase Chain Reactions (PCRs) are commonly used in population studies of aeromonads to detect virulence-associated genes in order to better understand the epidemiology and emergence of Aeromonas from the environment to host, but their performances have never been thoroughly evaluated. We aimed to determine diagnostic sensitivity and specificity of PCR assays for the detection of virulence-associated genes in a collection of Aeromonas isolates representative for the genetic diversity in the genus. Thirty-nine Aeromonas strains belonging to 27 recognized species were screened by published PCR assays for virulence-associated genes (act, aerA, aexT, alt, ascFG, ascV, ast, lafA, lip, ser, stx1, stx2A). In parallel, homologues of the 12 putative virulence genes were searched from the genomes of the 39 strains. Of the 12 published PCR assays for virulence factors, the comparison of PCR results and genome analysis estimated diagnostic sensitivities ranging from 34% to 100% and diagnostic specificities ranged from 71% to 100% depending upon the gene. To improve the detection of virulence-associated genes in aeromonads, we have designed new primer pairs for aerA/act, ser, lafA, ascFG and ascV, which showed excellent diagnostic sensitivity and specificity. Altogether, the analysis of high quality genomic data, which are more and more easy to obtain, provides significant improvements in the genetic detection of virulence factors in bacterial strains.

Genomic characterization of extensively drug-resistant Acinetobacter baumannii strain, KAB03 belonging to ST451 from Korea

Extensively drug-resistant (XDR) Acinetobacter baumannii strains have emerged rapidly worldwide. The antibiotic resistance characteristics of XDR A. baumannii strains show regional differences; therefore, it is necessary to analyze both genomic and proteomic characteristics of emerging XDR A. baumannii clinical strains isolated in Korea to elucidate their multidrug resistance. Here, we isolated new sequence type of XDR A. baumannii clinical strain (KAB03) from Korean hospitals and performed comprehensive genome analyses. The strain belongs to new sequence type, ST451. Single nucleotide polymorphism (SNP) analysis with other types of A. baumannii strains revealed that KAB03 has unique SNP pattern in the regions of gyrB and gpi of MLST profiles. A. baumannii KAB03 harbours three antibiotic resistance islands (AbGRI1, 2, and 3). AbGRI1 harbours two copies of Tn2006 containing bla_OXA-23, which play an important role in antibiotic resistance. AbGRI2 possesses aminoglycoside resistant gene aph(3')-Ic and class A β-lactamase bla_TEM. AbGIR3 has macrolide resistant genes and aminoglycoside resistant gene armA. A. baumannii KAB03 harbours mutations in pmrB and pmrC, which are believed to confer colistin resistance. In addition, proteomic and transcriptional analysis of KAB03 confirmed that β-lactamases (ADC-73 and OXA-23), Ade efflux pumps (AdeIJK), outer membrane proteins (OmpA and OmpW), and colistin resistance genes (PmrCAB) were major proteins responsible for antibiotic resistance. Our proteogenomic results provide valuable information for multi-drug resistance in emerging XDR A. baumannii strains belonging to ST451.

Molecular characterization of invasive meningococcal isolates in Burkina Faso as the relative importance of serogroups X and W increases, 2008-2012

Background

Neisseria meningitidis serogroup A disease in Burkina Faso has greatly decreased following introduction of a meningococcal A conjugate vaccine in 2010, yet other serogroups continue to pose a risk of life-threatening disease. Capsule switching among epidemic-associated serogroup A N. meningitidis strains could allow these lineages to persist despite vaccination. The introduction of new strains at the national or sub-national levels could affect the epidemiology of disease.

Methods

Isolates collected from invasive meningococcal disease in Burkina Faso between 2008 and 2012 were characterized by serogrouping and molecular typing. Genome sequences from a subset of isolates were used to infer phylogenetic relationships.

Results

The ST-5 clonal complex (CC5) was identified only among serogroup A isolates, which were rare after 2010. CC181 and CC11 were the most common clonal complexes after 2010, having serogroup X and W isolates, respectively. Whole-genome phylogenetic analysis showed that the CC181 isolates collected during and after the epidemic of 2010 formed a single clade that was closely related to isolates collected in Niger during 2005 and Burkina Faso during 2007. Geographic population structure was identified among the CC181 isolates, where pairs of isolates collected from the same region of Burkina Faso within a single year had less phylogenetic diversity than the CC181 isolate collection as a whole. However, the reduction of phylogenetic diversity within a region did not extend across multiple years. Instead, CC181 isolates collected during the same year had lower than average diversity, even when collected from different regions, indicating geographic mixing of strains across years. The CC11 isolates were primarily collected during the epidemic of 2012, with sparse sampling during 2011. These isolates belong to a clade that includes previously described isolates collected in Burkina Faso, Mali, and Niger from 2011 to 2015. Similar to CC181, reduced phylogenetic diversity was observed among CC11 isolate pairs collected from the same regions during a single year.

Conclusions

The population of disease-associated N. meningitidis strains within Burkina Faso was highly dynamic between 2008 and 2012, reflecting both vaccine-imposed selection against serogroup A strains and potentially complex clonal waves of serogroup X and serogroup W strains.

Electronic supplementary material

The online version of this article (10.1186/s12879-018-3247-x) contains supplementary material, which is available to authorized users.

Characterization of a novel plasmid encoding F4-like fimbriae present in a Shiga-toxin producing enterotoxigenic Escherichia coli isolated during the investigation on a case of hemolytic-uremic syndrome

In February 2017 a case of Hemolytic-Uremic Syndrome (HUS) was reported to the National Registry of HUS in an adult living in Northern Italy. Stool specimens from the patient and his family contacts were collected and the analyses led to the isolation of a Locus of Enterocyte Effacement (LEE)-negative Shiga toxin 2 (Stx2)-producing Escherichia coli. The epidemiological investigations performed brought to collect fecal samples from the animals reared in a farm held by the case's family and a mixture of bovine and swine feces proved positive for Shiga toxin-producing E. coli (STEC) and yielded the isolation of a LEE-negative stx2-positive E. coli strain. Further characterization by whole genome sequencing led to identify the isolates as two identical O2:H27 hybrid Enterotoxigenic Shiga toxin-producing E. coli (ETEC-STEC). Sequencing of a high molecular weight plasmid present in the human isolate disclosed a peculiar plasmid harboring virulence genes characteristic for both pathotypes, including the enterohemolysin-coding gene and sta1, encoding the heat stable enterotoxin. Moreover, a complete fae locus encoding the ETEC F4 fimbriae could be identified, including a novel variant of faeG gene responsible for the production of the main structural subunit of the fimbriae. This novel faeG showed great diversity in the nucleotidic sequence when compared with the reference genes encoding the swine F4 allelic variants, whereas at the amino acid sequence level the predicted protein sequence showed some similarity with FaeG from E. coli strains of bovine origin. Further investigation on the plasmid region harboring the newly identified faeG allelic variant allowed to identify similar plasmids in NCBI sequence database, as part of the genome of other previously uncharacterized ETEC-STEC strains of bovine origin, suggesting that the novel F4-like fimbriae may play a role in bovine host specificity.

Transcriptome Analysis of Neisseria gonorrhoeae during Natural Infection Reveals Differential Expression of Antibiotic Resistance Determinants between Men and Women

Recent emergence of antimicrobial resistance of Neisseria gonorrhoeae worldwide has resulted in limited therapeutic choices for treatment of infections caused by this organism. We performed global transcriptomic analysis of N. gonorrhoeae in subjects with gonorrhea who attended a Nanjing, China, sexually transmitted infection (STI) clinic, where antimicrobial resistance of N. gonorrhoeae is high and increasing. We found that N. gonorrhoeae transcriptional responses to infection differed in genital specimens taken from men and women, particularly antibiotic resistance gene expression, which was increased in men. These sex-specific findings may provide a new approach to guide therapeutic interventions and preventive measures that are also sex specific while providing additional insight to address antimicrobial resistance of N. gonorrhoeae.

Neisseria gonorrhoeae is a bacterial pathogen responsible for the sexually transmitted infection gonorrhea. Emergence of antimicrobial resistance (AMR) of N. gonorrhoeae worldwide has resulted in limited therapeutic choices for this infection. Men who seek treatment often have symptomatic urethritis; in contrast, gonococcal cervicitis in women is usually minimally symptomatic, but may progress to pelvic inflammatory disease. Previously, we reported the first analysis of gonococcal transcriptome expression determined in secretions from women with cervical infection. Here, we defined gonococcal global transcriptional responses in urethral specimens from men with symptomatic urethritis and compared these with transcriptional responses in specimens obtained from women with cervical infections and in vitro-grown N. gonorrhoeae isolates. This is the first comprehensive comparison of gonococcal gene expression in infected men and women. RNA sequencing analysis revealed that 9.4% of gonococcal genes showed increased expression exclusively in men and included genes involved in host immune cell interactions, while 4.3% showed increased expression exclusively in women and included phage-associated genes. Infected men and women displayed comparable antibiotic-resistant genotypes and in vitro phenotypes, but a 4-fold higher expression of the Mtr efflux pump-related genes was observed in men. These results suggest that expression of AMR genes is programed genotypically and also driven by sex-specific environments. Collectively, our results indicate that distinct N. gonorrhoeae gene expression signatures are detected during genital infection in men and women. We propose that therapeutic strategies could target sex-specific differences in expression of antibiotic resistance genes.

IMPORTANCE Recent emergence of antimicrobial resistance of Neisseria gonorrhoeae worldwide has resulted in limited therapeutic choices for treatment of infections caused by this organism. We performed global transcriptomic analysis of N. gonorrhoeae in subjects with gonorrhea who attended a Nanjing, China, sexually transmitted infection (STI) clinic, where antimicrobial resistance of N. gonorrhoeae is high and increasing. We found that N. gonorrhoeae transcriptional responses to infection differed in genital specimens taken from men and women, particularly antibiotic resistance gene expression, which was increased in men. These sex-specific findings may provide a new approach to guide therapeutic interventions and preventive measures that are also sex specific while providing additional insight to address antimicrobial resistance of N. gonorrhoeae.

Absence of Curli in Soil-Persistent Escherichia coli Is Mediated by a C-di-GMP Signaling Defect and Suggests Evidence of Biofilm-Independent Niche Specialization

Escherichia coli is commonly viewed as a gastrointestinal commensal or pathogen although an increasing body of evidence suggests that it can persist in non-host environments as well. Curli are a major component of biofilm in many enteric bacteria including E. coli and are important for adherence to different biotic and abiotic surfaces. In this study we investigated curli production in a unique collection of soil-persistent E. coli isolates and examined the role of curli formation in environmental persistence. Although most soil-persistent E. coli were curli-positive, 10% of isolates were curli-negative (17 out of 170). Curli-producing E. coli (COB583, COB585, and BW25113) displayed significantly more attachment to quartz sand than the curli-negative strains. Long-term soil survival experiments indicated that curli production was not required for long-term survival in live soil (over 110 days), as a curli-negative mutant BW25113ΔcsgB had similar survival compared to wild type BW25113. Mutations in two genes associated with c-di-GMP metabolism, dgcE and pdeR, correlated with loss of curli in eight soil-persistent strains, although this did not significantly impair their survival in soil compared to curli-positive strains. Overall, the data indicate that curli-deficient and biofilm-defective strains, that also have a defect in attachment to quartz sand, are able to reside in soil for long periods of time thus pointing to the possibility that niches may exist in the soil that can support long-term survival independently of biofilm formation.

Complete Genome Sequence of Escherichia coli ME8067, an Azide-Resistant Laboratory Strain Used for Conjugation Experiments

We report here the complete genome sequence of Escherichia coli ME8067, an azide-resistant laboratory strain used for conjugation experiments. The ME8067 genome was closely related to E. coli strain K-12 substrain W3110. This genome sequence will support further genetic analysis of conjugative elements.

Codon usage revisited: Lack of correlation between codon usage and the number of tRNA genes in enterobacteria

It is widely believed that if a high number of genes are found for any tRNA in a rapidly replicating bacteria, then the cytoplasmic levels of that tRNA will be high and an open reading frame containing a higher frequency of the complementary codon will be translated faster. This idea is based on correlations between the number of tRNA genes, tRNA concentration and the frequency of codon usage observed in a limited number of strains as well as from the fact that artificially changing the number of tRNA genes alters translation efficiency and consequently the amount of properly folded protein synthesized. tRNA gene number may greatly vary in a genome due to duplications, deletions and lateral transfer which in turn would alter the levels and functionality of many proteins. Such changes are potentially deleterious for fitness and as a result it is expected that changes in tRNA gene numbers should be accompanied by a modification of the frequency of codon usage. In contrast to this model, when comparing the number of tRNA genes and the frequency of codon usage of several Salmonella enterica and Escherichia coli strains we found that changes in the number of tRNA genes are not correlated to changes in codon usage. Furthermore, these changes are not correlated with a change in the efficiency of codon translation. These results suggest that once a genome gains or loses tRNA genes, it responds by modulating the concentrations of tRNAs rather than modifying its frequency of codon usage.

Pseudomonas aeruginosa L10: A Hydrocarbon-Degrading, Biosurfactant-Producing, and Plant-Growth-Promoting Endophytic Bacterium Isolated From a Reed (Phragmites australis)

Bacterial endophytes with the capacity to degrade petroleum hydrocarbons and promote plant growth may facilitate phytoremediation for the removal of petroleum hydrocarbons from contaminated soils. A hydrocarbon-degrading, biosurfactant-producing, and plant-growth-promoting endophytic bacterium, Pseudomonas aeruginosa L10, was isolated from the roots of a reed, Phragmites australis, in the Yellow River Delta, Shandong, China. P. aeruginosa L10 efficiently degraded C₁₀-C₂₆n-alkanes from diesel oil, as well as common polycyclic aromatic hydrocarbons (PAHs) such as naphthalene, phenanthrene, and pyrene. In addition, P. aeruginosa L10 could produce biosurfactant, which was confirmed by the oil spreading method, and surface tension determination of inocula. Moreover, P. aeruginosa L10 had plant growth-stimulating attributes, including siderophore and indole-3-acetic acid (IAA) release, along with 1-aminocyclopropane-1-carboxylic (ACC) deaminase activity. To explore the mechanisms underlying the phenotypic traits of endophytic P. aeruginosa L10, we sequenced its complete genome. From the genome, we identified genes related to petroleum hydrocarbon degradation, such as putative genes encoding monooxygenase, dioxygenase, alcohol dehydrogenase, and aldehyde dehydrogenase. Genome annotation revealed that P. aeruginosa L10 contained a gene cluster involved in the biosynthesis of rhamnolipids, rhlABRI, which should be responsible for the observed biosurfactant activity. We also identified two clusters of genes involved in the biosynthesis of siderophore (pvcABCD and pchABCDREFG). The genome also harbored tryptophan biosynthetic genes (trpAB, trpDC, trpE, trpF, and trpG) that are responsible for IAA synthesis. Moreover, the genome contained the ACC deaminase gene essential for ACC deaminase activity. This study will facilitate applications of endophytic P. aeruginosa L10 to phytoremediation by advancing the understanding of hydrocarbon degradation, biosurfactant synthesis, and mutualistic interactions between endophytes and host plants.

Complete Genome Sequence of Escherichia coli J53, an Azide-Resistant Laboratory Strain Used for Conjugation Experiments

We report here the complete genome sequence of Escherichia coli J53, which is used as a recipient in conjugation experiments and is a laboratory strain derived from E. coli K-12. This genome sequence will help in the development of a comprehensive genetic analysis of conjugative elements.

Neisseria genomics: current status and future perspectives

High-throughput whole genome sequencing has unlocked a multitude of possibilities enabling members of the Neisseria genus to be examined with unprecedented detail, including the human pathogens Neisseria meningitidis and Neisseria gonorrhoeae. To maximise the potential benefit of this for public health, it is becoming increasingly important to ensure that this plethora of data are adequately stored, disseminated and made readily accessible. Investigations facilitating cross-species comparisons as well as the analysis of global datasets will allow differences among and within species and across geographic locations and different times to be identified, improving our understanding of the distinct phenotypes observed. Recent advances in high-throughput platforms that measure the transcriptome, proteome and/or epigenome are also becoming increasingly employed to explore the complexities of Neisseria biology. An integrated approach to the analysis of these is essential to fully understand the impact these may have in the Neisseria genus. This article reviews the current status of some of the tools available for next generation sequence analysis at the dawn of the ‘post-genomic’ era.

Genomic Analysis of Salmonella enterica Serovar Typhimurium DT160 Associated with a 14-Year Outbreak, New Zealand, 1998-2012

During 1998-2012, an extended outbreak of Salmonella enterica serovar Typhimurium definitive type 160 (DT160) affected >3,000 humans and killed wild birds in New Zealand. However, the relationship between DT160 within these 2 host groups and the origin of the outbreak are unknown. Whole-genome sequencing was used to compare 109 Salmonella Typhimurium DT160 isolates from sources throughout New Zealand. We provide evidence that DT160 was introduced into New Zealand around 1997 and rapidly propagated throughout the country, becoming more genetically diverse over time. The genetic heterogeneity was evenly distributed across multiple predicted functional protein groups, and we found no evidence of host group differentiation between isolates collected from human, poultry, bovid, and wild bird sources, indicating ongoing transmission between these host groups. Our findings demonstrate how a comparative genomic approach can be used to gain insight into outbreaks, disease transmission, and the evolution of a multihost pathogen after a probable point-source introduction.

Invasive Serotype 35B Pneumococci Including an Expanding Serotype Switch Lineage, United States, 2015-2016

We used whole-genome sequencing to characterize 199 nonvaccine serotype 35B pneumococcal strains that caused invasive pneumococcal disease (IPD) in the United States during 2015–2016 and related these findings to previous serotype 35B IPD data obtained by Active Bacterial Core surveillance. Penicillin-nonsusceptible 35B IPD increased during post–pneumococcal 7-valent conjugate vaccine years (2001–2009) and increased further after implementation of pneumococcal 13-valent conjugate vaccine in 2010. This increase was caused primarily by the 35B/sequence type (ST) 558 lineage. 35B/ST558 and vaccine serotype 9V/ST156 lineages were implicated as cps35B donor and recipient, respectively, for a single capsular switch event that generated emergent 35B/ST156 progeny in 6 states during 2015–2016. Three additional capsular switch 35B variants were identified, 2 of which also involved 35B/ST558 as cps35B donor. Spread of 35B/ST156 is of concern in view of past global predominance of pathogenic ST156 vaccine serotype strains. Protection against serotype 35B should be considered in next-generation pneumococcal vaccines.

Emergence of quinolone-resistant Shigella flexneri in Italy (March 2017)

In March 2017, a 45-year-old Italian man who has sex with men was admitted to the Infectious Diseases Department of Trieste Hospital (northeast Italy), because of fever, abdominal pain and dysentery. The patient had neither foreign travel history nor sexual contact with non-Italian partners. Stool cultures grew multidrug-resistant Shigella flexneri (resistant to ampicillin, chloramphenicol, streptomycin, tetracycline, trimethoprim, amoxicillin/clavulanic acid and ciprofloxacin) and whole genome sequencing detailed the resistance features. The phylogenetic analysis showed that the strain was unrelated to any previously reported strain. The patient was treated successfully with ceftriaxone. We hereby report the first case of locally-acquired, multidrug-resistant S. flexneri infection in Italy and also the emergence of a new clone.

Characteristics of Carbapenemase-Producing Enterobacteriaceae in Wastewater Revealed by Genomic Analysis

Wastewater is considered a major source of antibiotic-resistant bacteria released into the environment. Here, we characterized carbapenemase-producing Enterobacteriaceae (CPE) in wastewater by whole-genome analysis. Wastewater samples (n = 40) were collected from municipal wastewater treatment plants and hospital wastewater in Japan and Taiwan. Samples were screened for CPE using selective media, and the obtained isolates were sequenced using an Illumina MiSeq. The isolates (n = 45) included the following microorganisms: Klebsiella quasipneumoniae (n = 12), Escherichia coli (n = 10), Enterobacter cloacae complex (n = 10), Klebsiella pneumoniae (n = 8), Klebsiella variicola (n = 2), Raoultella ornithinolytica (n = 1), Citrobacter freundii (n = 1), and Citrobacter amalonaticus (n = 1). Among the 45 isolates, 38 harbored at least one carbapenemase-encoding gene. Of these, the bla_GES (bla_GES-5, bla_GES-6, and bla_GES-24) genes were found in 29 isolates. The genes were situated in novel class 1 integrons, but the integron structures were different between the Japanese (In1439 with bla_GES-24 and In1440 with bla_GES-5) and Taiwanese (In1441 with bla_GES-5 and In1442 with bla_GES-6) isolates. Other carbapenemase-encoding genes (bla_VIM-1, bla_NDM-5, bla_IMP-8, bla_IMP-19, and bla_KPC-2) were found in one to three isolates. Notably, class 1 integrons previously reported among clinical isolates obtained in the same regions as the present study, namely, In477 with bla_IMP-19 and In73 with bla_IMP-8, were found among the Japanese and Taiwanese isolates, respectively. The results indicate that CPE with various carbapenemase-encoding genes in different genetic contexts were present in biologically treated wastewater, highlighting the need to monitor for antibiotic resistance in wastewater.

Phylogeny and putative virulence gene analysis of Bartonella bovis

Bartonella bovis is a small Gram-negative bacterium recognized as an etiological agent for bacteremia and endocarditis in cattle. As few reports are available on the taxonomic position of B. bovis and its mechanism of virulence, this study aims to resolve the phylogeny of B. bovis and investigate putative virulence genes based on whole genome sequence analysis. Genome-wide comparisons based on single nucleotide polymorphisms (SNP) and orthologous genes were performed in this study for phylogenetic inference of 27 Bartonella species. Rapid Annotation using Subsystem Technology (RAST) analysis was used for annotation of putative virulence genes. The phylogenetic tree generated from the genome-wide comparison of orthologous genes exhibited a topology almost similar to that of the tree generated from SNP-based comparison, indicating a high concordance in the nucleotide and amino acid sequences of Bartonella spp. The analyses show consistent grouping of B. bovis in a cluster related to ruminant-associated species, including Bartonella australis, Bartonella melophagi and Bartonella schoenbuchensis. RAST analysis revealed genes encoding flagellar components, in corroboration with the observation of flagella-like structure of BbUM strain under negative straining. Genes associated with virulence, disease and defence, prophages, membrane transport, iron acquisition, motility and chemotaxis are annotated in B. bovis genome. The flagellin (flaA) gene of B. bovis is closely related to Bartonella bacilliformis and Bartonella clarridgeiae but distinct from other Gram-negative bacteria. The absence of type IV secretion systems, the bona fide pathogenicity factors of bartonellae, in B. bovis suggests that it may have a different mechanism of pathogenicity.