Development of a reference data set for assigning Streptococcus and Enterococcus species based on next generation sequencing of the 16S-23S rRNA region.

Background: Many members of Streptococcus and Enterococcus genera are clinically relevant opportunistic pathogens warranting accurate and rapid identification for targeted therapy. Currently, the developed method based on next generation sequencing (NGS) of the 16S-23S rRNA region proved to be a rapid, reliable and precise approach for species identification directly from polymicrobial and challenging clinical samples. The introduction of this new method to routine diagnostics is hindered by a lack of the reference sequences for the 16S-23S rRNA region for many bacterial species. The aim of this study was to develop a careful assignment for streptococcal and enterococcal species based on NGS of the 16S-23S rRNA region. Methods: Thirty two strains recovered from clinical samples and 19 reference strains representing 42 streptococcal species and nine enterococcal species were subjected to bacterial identification by four Sanger-based sequencing methods targeting the genes encoding (i) 16S rRNA, (ii) sodA, (iii) tuf and (iv) rpoB; and NGS of the 16S-23S rRNA region. Results: This study allowed obtainment and deposition of reference sequences of the 16S-23S rRNA region for 15 streptococcal and 3 enterococcal species followed by enrichment for 27 and 6 species, respectively, for which reference sequences were available in the databases. For Streptococcus, NGS of the 16S-23S rRNA region was as discriminative as Sanger sequencing of the tuf and rpoB genes allowing for an unambiguous identification of 93% of analyzed species. For Enterococcus, sodA, tuf and rpoB genes sequencing allowed for identification of all species, while the NGS-based method did not allow for identification of only one enterococcal species. For both genera, the sequence analysis of the 16S rRNA gene was endowed with a low identification potential and was inferior to that of other tested identification methods. Moreover, in case of phylogenetically related species the sequence analysis of only the intergenic spacer region was not sufficient enough to precisely identify Streptococcus strains at the species level. Conclusions: Based on the developed reference dataset, clinically relevant streptococcal and enterococcal species can now be reliably identified by 16S-23S rRNA sequences in samples. This study will be useful for introduction of a novel diagnostic tool, NGS of the 16S-23S rRNA region, which undoubtedly is an improvement for reliable culture-independent species identification directly from polymicrobially constituted clinical samples.

Development and Validation of a Reference Data Set for Assigning Staphylococcus Species Based on Next-Generation Sequencing of the 16S-23S rRNA Region.

Many members of the Staphylococcus genus are clinically relevant opportunistic pathogens that warrant accurate and rapid identification for targeted therapy. The aim of this study was to develop a careful assignment scheme for staphylococcal species based on next-generation sequencing (NGS) of the 16S-23S rRNA region. All reference staphylococcal strains were identified at the species level using Sanger sequencing of the 16S rRNA, sodA, tuf, and rpoB genes and NGS of the 16S-23S rRNA region. To broaden the database, an additional 100 staphylococcal strains, including 29 species, were identified by routine diagnostic methods, 16S rRNA Sanger sequencing and NGS of the 16S-23S rRNA region. The results enabled development of reference sequences encompassing the 16S-23S rRNA region for 50 species (including one newly proposed species) and 6 subspecies of the Staphylococcus genus. This study showed sodA and rpoB targets were the most discriminative but NGS of the 16S-23S rRNA region was more discriminative than tuf gene sequencing and much more discriminative than 16S rRNA gene sequencing. Almost all Staphylococcus species could be distinguished when the max score was 99.0% or higher and the sequence similarity between the best and second best species was equal to or >0.2% (min. 9 nucleotides). This study allowed development of reference sequences for 21 staphylococcal species and enrichment for 29 species for which sequences were publicly available. We confirmed the usefulness of NGS of the 16S-23S rRNA region by identifying the whole species content in 45 clinical samples and comparing the results to those obtained using routine diagnostic methods. Based on the developed reference database, all staphylococcal species can be reliably detected based on the 16S-23S rRNA sequences in samples composed of both single species and more complex polymicrobial communities. This study will be useful for introduction of a novel diagnostic tool, which undoubtedly is an improvement for reliable species identification in polymicrobial samples. The introduction of this new method is hindered by a lack of reference sequences for the 16S-23S rRNA region for many bacterial species. The results will allow identification of all Staphylococcus species, which are clinically relevant pathogens.

Targeted next-generation sequencing of the 16S-23S rRNA region for culture-independent bacterial identification - increased discrimination of closely related species.

The aim of this study was to develop an easy-to-use culture-free diagnostic method based on next generation sequencing (NGS) of PCR amplification products encompassing whole 16S-23S rRNA region to improve the resolution of bacterial species identification. To determine the resolution of the new method 67 isolates were subjected to four identification methods: Sanger sequencing of the 16S rRNA gene; NGS of the 16S-23S rRNA region using MiSeq (Illumina) sequencer; Microflex MS (Bruker) and VITEK MS (bioMérieux). To evaluate the performance of this new method when applied directly on clinical samples, we conducted a proof of principle study with 60 urine samples from patients suspected of urinary tract infections (UTIs), 23 BacT/ALERT (bioMérieux) positive blood culture bottles and 21 clinical orthopedic samples. The resolution power of NGS of the 16S-23S rRNA region was superior to other tested identification methods. Furthermore, the new method correctly identified pathogens established as the cause of UTIs and blood stream infections with conventional culture. NGS of the 16S-23S rRNA region also showed increased detection of bacterial microorganisms in clinical samples from orthopedic patients. Therefore, we conclude that our method has the potential to increase diagnostic yield for detection of bacterial pathogenic species compared to current methods.

Reprint of "Application of next generation sequencing in clinical microbiology and infection prevention".

Current molecular diagnostics of human pathogens provide limited information that is often not sufficient for outbreak and transmission investigation. Next generation sequencing (NGS) determines the DNA sequence of a complete bacterial genome in a single sequence run, and from these data, information on resistance and virulence, as well as information for typing is obtained, useful for outbreak investigation. The obtained genome data can be further used for the development of an outbreak-specific screening test. In this review, a general introduction to NGS is presented, including the library preparation and the major characteristics of the most common NGS platforms, such as the MiSeq (Illumina) and the Ion PGM™ (ThermoFisher). An overview of the software used for NGS data analyses used at the medical microbiology diagnostic laboratory in the University Medical Center Groningen in The Netherlands is given. Furthermore, applications of NGS in the clinical setting are described, such as outbreak management, molecular case finding, characterization and surveillance of pathogens, rapid identification of bacteria using the 16S-23S rRNA region, taxonomy, metagenomics approaches on clinical samples, and the determination of the transmission of zoonotic micro-organisms from animals to humans. Finally, we share our vision on the use of NGS in personalised microbiology in the near future, pointing out specific requirements.

Application of next generation sequencing in clinical microbiology and infection prevention.

Current molecular diagnostics of human pathogens provide limited information that is often not sufficient for outbreak and transmission investigation. Next generation sequencing (NGS) determines the DNA sequence of a complete bacterial genome in a single sequence run, and from these data, information on resistance and virulence, as well as information for typing is obtained, useful for outbreak investigation. The obtained genome data can be further used for the development of an outbreak-specific screening test. In this review, a general introduction to NGS is presented, including the library preparation and the major characteristics of the most common NGS platforms, such as the MiSeq (Illumina) and the Ion PGM™ (ThermoFisher). An overview of the software used for NGS data analyses used at the medical microbiology diagnostic laboratory in the University Medical Center Groningen in The Netherlands is given. Furthermore, applications of NGS in the clinical setting are described, such as outbreak management, molecular case finding, characterization and surveillance of pathogens, rapid identification of bacteria using the 16S-23S rRNA region, taxonomy, metagenomics approaches on clinical samples, and the determination of the transmission of zoonotic micro-organisms from animals to humans. Finally, we share our vision on the use of NGS in personalised microbiology in the near future, pointing out specific requirements.

Virulence, Antimicrobial Resistance Properties and Phylogenetic Background of Non-H7 Enteropathogenic Escherichia coli O157.

Escherichia coli (E.coli) O157 that do not produce Shiga toxin and do not possess flagellar antigen H7 are of diverse H serotypes. In this study, the antibiotic resistance properties, genotype of a set of virulence associated genes and the phylogenetic background of E. coli O157:non-H7 groups were compared. Whole genome sequencing was performed on fourteen O157:non-H7 isolates collected in the STEC-ID-net study. The genomes were compared with E. coli O157 genomes and a typical Enteropathogenic E. coli (tEPEC) genome downloaded from NCBI. Twenty-six (86%) of the analyzed genomes had the intimin encoding gene eae but of different types mostly correlating with their H types, e.g., H16, H26, H39, and H45 carried intimin type ε, ß, κ, and α, respectively. They belonged to several E. coli phylogenetic groups, i.e., to phylogenetic group A, B1, B2, and D. Seven (50%) of our collected O157:non-H7 isolates were resistant to two or more antibiotics. Several mobile genetic elements, such as plasmids, insertion elements, and pathogenicity islands, carrying a set of virulence and resistance genes were found in the E. coli O157:non-H7 isolates. Core genome phylogenetic analysis showed that O157:non-H7 isolates probably evolved from different phylogenetic lineages and were distantly related to the E. coli O157:H7 lineage. We hypothesize that independent acquisition of mobile genetic elements by isolates of different lineages have contributed to the different molecular features of the O157:non-H7 strains. Although distantly related to the STEC O157, E. coli O157:non-H7 isolates from multiple genetic background could be considered as pathogen of concern for their diverse virulence and antibiotic resistance properties.

Comprehensive Characterization of Escherichia coli O104:H4 Isolated from Patients in the Netherlands.

In 2011, a Shiga toxin-producing Enteroaggregative Escherichia coli (EAEC Stx2a+) O104:H4 strain caused a serious outbreak of acute gastroenteritis and hemolytic-uremic syndrome (HUS) in Germany. In 2013, E. coli O104:H4 isolates were obtained from a patient with HUS and her friend showing only gastrointestinal complaints. The antimicrobial resistance and virulence profiles of these isolates together with three EAEC Stx2a+ O104:H4 isolates from 2011 were determined and compared. Whole-genome sequencing (WGS) was performed for detailed characterization and to determine genetic relationship of the isolates. Four additional genomes of EAEC Stx2a+ O104:H4 isolates of 2009 and 2011 available on NCBI were included in the virulence and phylogenetic analysis. All E. coli O104:H4 isolates tested were positive for stx2a, aatA, and terD but were negative for escV. All, except one 2011 isolate, were positive for aggR and were therefore considered EAEC. The EAEC Stx2a+ O104:H4 isolates of 2013 belonged to sequence type (ST) ST678 as the 2011 isolates and showed slightly different resistance and virulence patterns compared to the 2011 isolates. Core-genome phylogenetic analysis showed that the isolates of 2013 formed a separate cluster from the isolates of 2011 and 2009 by 27 and 20 different alleles, respectively. In addition, only a one-allele difference was found between the isolate of the HUS-patient and that of her friend. Our study shows that EAEC Stx2a+ O104:H4 strains highly similar to the 2011 outbreak clone in their core genome are still circulating necessitating proper surveillance to prevent further outbreaks with these potentially pathogenic strains. In addition, WGS not only provided a detailed characterization of the isolates but its high discriminatory power also enabled us to discriminate the 2013 isolates from the isolates of 2009 and 2011 expediting the use of WGS in public health services to rapidly apply proper infection control strategies.

Is Shiga Toxin-Negative Escherichia coli O157:H7 Enteropathogenic or Enterohemorrhagic Escherichia coli? Comprehensive Molecular Analysis Using Whole-Genome Sequencing.

The ability of Escherichia coli O157:H7 to induce cellular damage leading to disease in humans is related to numerous virulence factors, most notably the stx gene, encoding Shiga toxin (Stx) and carried by a bacteriophage. Loss of the Stx-encoding bacteriophage may occur during infection or culturing of the strain. Here, we collected stx-positive and stx-negative variants of E. coli O157:H7/NM (nonmotile) isolates from patients with gastrointestinal complaints. Isolates were characterized by whole-genome sequencing (WGS), and their virulence properties and phylogenetic relationship were determined. Because of the presence of the eae gene but lack of the bfpA gene, the stx-negative isolates were considered atypical enteropathogenic E. coli (aEPEC). However, they had phenotypic characteristics similar to those of the Shiga toxin-producing E. coli (STEC) isolates and belonged to the same sequence type, ST11. Furthermore, EPEC and STEC isolates shared similar virulence genes, the locus of enterocyte effacement region, and plasmids. Core genome phylogenetic analysis using a gene-by-gene typing approach showed that the sorbitol-fermenting (SF) stx-negative isolates clustered together with an SF STEC isolate and that one non-sorbitol-fermenting (NSF) stx-negative isolate clustered together with NSF STEC isolates. Therefore, these stx-negative isolates were thought either to have lost the Stx phage or to be a progenitor of STEC O157:H7/NM. As detection of STEC infections is often based solely on the identification of the presence of stx genes, these may be misdiagnosed in routine laboratories. Therefore, an improved diagnostic approach is required to manage identification, strategies for treatment, and prevention of transmission of these potentially pathogenic strains.

Assessing the public health risk of Shiga toxin-producing Escherichia coli by use of a rapid diagnostic screening algorithm.

Shiga toxin-producing Escherichia coli (STEC) is an enteropathogen of public health concern because of its ability to cause serious illness and outbreaks. In this prospective study, a diagnostic screening algorithm to categorize STEC infections into risk groups was evaluated. The algorithm consists of prescreening stool specimens with real-time PCR (qPCR) for the presence of stx genes. The qPCR-positive stool samples were cultured in enrichment broth and again screened for stx genes and additional virulence factors (escV, aggR, aat, bfpA) and O serogroups (O26, O103, O104, O111, O121, O145, O157). Also, PCR-guided culture was performed with sorbitol MacConkey agar (SMAC) and CHROMagar STEC medium. The presence of virulence factors and O serogroups was used for presumptive pathotype (PT) categorization in four PT groups. The potential risk for severe disease was categorized from high risk for PT group I to low risk for PT group III, whereas PT group IV consists of unconfirmed stx qPCR-positive samples. In total, 5,022 stool samples of patients with gastrointestinal symptoms were included. The qPCR detected stx genes in 1.8% of samples. Extensive screening for virulence factors and O serogroups was performed on 73 samples. After enrichment, the presence of stx genes was confirmed in 65 samples (89%). By culture on selective media, STEC was isolated in 36% (26/73 samples). Threshold cycle (CT) values for stx genes were significantly lower after enrichment compared to direct qPCR (P < 0.001). In total, 11 (15%), 19 (26%), 35 (48%), and 8 (11%) samples were categorized into PT groups I, II, III, and IV, respectively. Several virulence factors (stx2, stx2a, stx2f, toxB, eae, efa1, cif, espA, tccP, espP, nleA and/or nleB, tir cluster) were associated with PT groups I and II, while others (stx1, eaaA, mch cluster, ireA) were associated with PT group III. Furthermore, the number of virulence factors differed between PT groups (analysis of variance, P < 0.0001). In conclusion, a diagnostic algorithm enables fast discrimination of STEC infections associated with a high to moderate risk for severe disease (PT groups I and II) from less-virulent STEC (PT group III).

Detection of Campylobacter species and Arcobacter butzleri in stool samples by use of real-time multiplex PCR.

The presence of Campylobacter (or Campylobacter-like) species in stools from patients suspected of infectious gastroenteritis (n = 493) was investigated using real-time PCR for detection of Arcobacter butzleri (hsp60 gene), Campylobacter coli (ceuE gene), Campylobacter jejuni (mapA), five acknowledged pathogenic Campylobacter spp. (C16S_Lund assay), and the Campylobacter genus (C16S_LvI assay). In total, 71.4% of the samples were positive for Campylobacter DNA (n = 352) by a Campylobacter genus-specific (C16S_LvI) assay. A total of 23 samples (4.7%) were positive in the C16S_Lund assay, used for detection of C. jejuni, C. coli, C. lari, C. upsaliensis, and C. hyointestinalis. Subsequent identification of these samples yielded detection frequencies (DF) of 4.1% (C. jejuni), 0.4% (C. coli), and 0.4% (C. upsaliensis). The DF of A. butzleri was 0.4%. Interestingly, sequencing of a subgroup (n = 46) of C16S_LvI PCR-positive samples resulted in a considerable number of Campylobacter concisus-positive samples (n = 20). PCR-positive findings with the C16S_Lund and C. jejuni/C. coli-specific assays were associated with more serious clinical symptoms (diarrhea and blood). Threshold cycle (C(T)) values of C. jejuni/C. coli PCR-positive samples were comparable to those of the C16S_Lund PCR (P = 0.21). C(T) values for both assays were significantly lower than those of the C16S_LvI assay (P < 0.001 and P < 0.00001, respectively). In conclusion, this study demonstrated that in combination, the C. jejuni/C coli-specific assays and the C16S_Lund assay are both useful for routine screening purposes. Furthermore, the DF of the emerging pathogen C. concisus was at least similar to the DF of C. jejuni.

The KIzSS network, a sentinel surveillance system for infectious diseases in day care centers: study protocol.

BACKGROUND: Day care-associated infectious diseases are widely recognized as a public health problem but rarely studied. Insights into their dynamics and their association with the day care setting are important for effective decision making in management of infectious disease control. This paper describes the purpose, design and potential of our national multi-center, day care-based sentinel surveillance network for infectious diseases (the KIzSS network). The aim of the KIzSS network is to acquire a long-term insight into the syndromic and microbiological aspects of day care-related infectious diseases and associated disease burden and to model these aspects with day care setting characteristics. METHODS/DESIGN: The KIzSS network applies a prospective cohort design, following day care centers rather than individual children or staff members over time. Data on infectious disease symptoms and related morbidity (children and staff), medical consumption, absenteeism and circulating enteric pathogens (children) are collected on a daily, weekly or monthly basis. Every two years, a survey is performed to assess the characteristics of participating day care centers. DISCUSSION: The KIzSS network offers a unique potential to study infectious disease dynamics in the day care setting over a sustained period of time. The created (bio)databases will help us to assess day care-related disease burden of infectious diseases among attending children and staff and their relation with the day care setting. This will support the much needed development of evidence-based and pragmatic guidelines for infectious disease control in day care centers.

Five commercial DNA extraction systems tested and compared on a stool sample collection.

In this study, 5 different commercial DNA extraction systems were tested on a stool sample collection containing 81 clinical stool specimens that were culture-positive for diarrheagenic Escherichia coli, Campylobacter jejuni, Salmonella enterica, or Clostridium difficile. The purified DNAs were analyzed by polymerase chain reaction (PCR) directed toward the relevant organisms. The results showed that conventional PCR combined with the extraction systems BioRobot EZ1 (Qiagen, Hilden, Germany), Bugs'n Beads (Genpoint, Oslo, Norway), ChargeSwitch (Invitrogen, Paisley, UK), QIAamp Stool Mini Kit (Qiagen), and 2 protocols (generic and Specific A) for EasyMag (BioMérieux, Marcy I'Etoile, France) were able to identify 89%, 62%, 85%, 88%, 85%, and 91%, respectively, of the pathogens originally identified by conventional culture-based methods. When TaqMan PCR was combined with the EasyMag Specific A protocol, 99% of the samples were correctly identified. The results demonstrate that the extraction efficiencies can vary significantly among different extraction systems, careful optimization may have a significant positive effect, and the use of sensitive and specific detection methods like TaqMan PCR is an ideal choice for this type of analysis.

Improved detection of five major gastrointestinal pathogens by use of a molecular screening approach.

The detection of bacterial and parasitic gastrointestinal pathogens through culture and microscopy is laborious and time-consuming. We evaluated a molecular screening approach (MSA) for the detection of five major enteric pathogens: Salmonella enterica, Campylobacter jejuni, Giardia lamblia, Shiga toxin-producing Escherichia coli (STEC), and Shigella spp./enteroinvasive E. coli (EIEC), for use in the daily practice of a clinical microbiology laboratory. The MSA consists of prescreening of stool specimens with two real-time multiplex PCR (mPCR) assays, which give results within a single working day, followed by guided culture/microscopy of the positive or mPCR-inhibited samples. In the present 2-year overview, 28,185 stool specimens were included. The MSA was applied to 13,974 stool samples (49.6%), whereas 14,211 samples were tested by conventional methods only (50.4%). The MSA significantly increased the total detection rate compared to that of conventional methods (19.2% versus 6.4%). The detection of all included pathogens, with the exception of S. enterica, significantly improved. MSA detection frequencies were as follows: C. jejuni, 8.1%; G. lamblia, 4.7%; S. enterica, 3.0%; STEC, 1.9%; and Shigella spp./EIEC, 1.4%. The guided culture/microscopy was positive in 76.8%, 58.1%, 88.9%, 16.8%, and 18.1% of mPCR-positive specimens, respectively. Of all mPCRs, only 1.8% was inhibited. Other findings were that detection of mixed infections was increased (0.9% versus 0.02%) and threshold cycle (C(T)) values for MSA guided culture/microscopy-positive samples were significantly lower than those for guided culture/microscopy-negative samples. In conclusion, an MSA for detection of gastrointestinal pathogens resulted in markedly improved detection rates and a substantial decrease in time to reporting of (preliminary) results.

Prospective study of use of PCR amplification and sequencing of 16S ribosomal DNA from cerebrospinal fluid for diagnosis of bacterial meningitis in a clinical setting.

We have evaluated the use of a broad-range PCR aimed at the 16S rRNA gene in detecting bacterial meningitis in a clinical setting. To achieve a uniform DNA extraction procedure for both gram-positive and gram-negative organisms, a combination of physical disruption (bead beating) and a silica-guanidiniumthiocyanate procedure was used for nucleic acid preparation. To diminish the risk of contamination as much as possible, we chose to amplify almost the entire 16S rRNA gene. The analytical sensitivity of the assay was approximately 1 x 10(2) to 2 x 10(2) CFU/ml of cerebrospinal fluid (CSF) for both gram-negative and gram-positive bacteria. In a prospective study of 227 CSF samples, broad-range PCR proved to be superior to conventional methods in detecting bacterial meningitis when antimicrobial therapy had already started. Overall, our assay showed a sensitivity of 86%, a specificity of 97%, a positive predictive value of 80%, and a negative predictive value of 98% compared to culture. We are currently adapting the standard procedures in our laboratory for detecting bacterial meningitis; broad-range 16S ribosomal DNA PCR detection is indicated when antimicrobial therapy has already started at time of lumbar puncture or when cultures remain negative, although the suspicion of bacterial meningitis remains.