Microbial sequence typing in the genomic era

Oligella spp.: A systematic review on an uncommon urinary pathogen

BACKGROUND

Oligella is an uncommon Gram-negative coccobacillus that was first thought to belong to the urogenital tract. The genus Oligella comprises two species that were recovered from various samples worldwide.

METHODS

We perform a systematic review focusing on Oligella microbiological characteristics, habitat, role in Human microbiome and infection, and antimicrobial susceptibility.

RESULTS

In humans, Oligella is mainly found as part of the microbiome of individuals with predisposing conditions. Oligella were also associated with invasive infections in patients with underlying diseases. Nevertheless, their prevalence remains to determine. Oligella culture requires up to 48 h on agar media in vitro, while urinary samples are usually incubated for 24 h. Consequently, microbiologists should be prompt to prolong the incubation of agar media when the direct examination showed Gram-negative coccobacilli. Oligella is accurately identified using MALDI-TOF mass spectrometry, but biochemical methods often provided inconsistent results. Specific guidelines for antimicrobial susceptibility testing of Oligella lack but the incubation could require up to 48 h of incubation. In contrast to O. urethralis, which is susceptible to third-generation cephalosporin, O. ureolytica is likely resistant to numerous antimicrobials. Genectic determinants of resistance were identified for beta-lactams and aminoglycosides.

CONCLUSION

Oligella is an uncommon pathogen that can be underrecognized. Microbiologists should be prompt to prolong the incubation of agar media plated with urines when the direct examination showed Gram-negative coccobacilli. Carbapenems should probably be given for the empirical treatment.

The correlation between intestinal colonization and infection of carbapenem-resistant Klebsiella pneumoniae: a systemactic review

As a widely spread Gram-negative bacteria, klebsiella pneumoniae mainly causes acquired infections in hospitals, such as lung infections, urinary tract infections, bloodstream infections, etc. In recent years, the number of multidrug-resistant K. pneumoniae strains has increased dramatically, posing a great threat to human health. Carbapenem-resistant Klebsiella pneumoniae (CRKP) can be colonized in human body, especially in gastrointestinal tract, and some colonized patients can be infected during hospitalization, among which invasive operation, underlying disease, admission to intensive care unit, antibiotic use, severity of the primary disease, advanced age, operation, coma and renal failure are common risk factors for secondary infection. Active screening and preventive measures can effectively prevent the occurrence of CRKP infection. Based on the epidemiological status, this study aims to discuss the correlation between colonization and secondary infection induced by carbapenem-resistant Klebsiella pneumoniae and risk factors for their happening, and provide some reference for nosocomial infection prevention and control.

Molecular epidemiological and antimicrobial-resistant mechanisms analysis of prolonged Neisseria gonorrhoeae collection between 1971 and 2005 in Japan

Objectives

As antimicrobial-resistant (AMR) Neisseria gonorrhoeae strains have emerged, humans have adjusted the antimicrobials used to treat infections. We identified shifts in the N. gonorrhoeae population and the determinants of AMR strains isolated during the recurring emergence of resistant strains and changes in antimicrobial therapies.

Methods

We examined 243 N. gonorrhoeae strains corrected at the Kanagawa Prefectural Institute of Public Health, Kanagawa, Japan, these isolated in 1971-2005. We performed multilocus sequence typing and AMR determinants (penA, mtrR, porB, ponA, 23S rRNA, gyrA and parC) mainly using high-throughput genotyping methods together with draft whole-genome sequencing on the MiSeq (Illumina) platform.

Results

All 243 strains were divided into 83 STs. ST1901 (n = 17) was predominant and first identified after 2001. Forty-two STs were isolated in the 1970s, 34 in the 1980s, 22 in the 1990s and 13 in the 2000s, indicating a decline in ST diversity over these decades. Among the 29 strains isolated after 2001, 28 were highly resistant to ciprofloxacin (MIC ≥ 8 mg/L) with two or more amino-acid substitutions in quinolone-resistance-determining regions. Seven strains belonging to ST7363 (n = 3), ST1596 (n = 3) and ST1901 (n = 1) were not susceptible to cefixime, and six strains carried penA alleles with mosaic-like penicillin-binding protein 2 (PBP2; penA 10.001 and 10.016) or PBP2 substitutions A501V and A517G.

Conclusions

We observed a significant reduction in the diversity of N. gonorrhoeae over 35 years in Japan. Since 2001, ST1901, which is resistant to ciprofloxacin, has superseded previous strains, becoming the predominant ST population.

Laboratory Methods for Culture and Identification

This chapter addresses the cultivation, identification, and characterization of Streptococcus suis. Here, we describe in detail the most used methodologies and expected results.

Contamination of Plant Foods with Bacillus cereus in a Province and Analysis of Its Traceability

Bacillus cereus is an important zoonotic foodborne conditional pathogen. It is found in vegetables, dairy products, rice, and other foods, thereby greatly endangering human health. Investigations on B. cereus contamination in China primarily focus on raw milk, dairy products, meat, and others, and limited research has been conducted on plant-based foodstuffs. The rapid development of sequencing technology and the application of bioinformatics-related techniques means that analysis based on whole-genome sequencing has become an important tool for the molecular-epidemiology investigation of B. cereus. In this study, we investigated the contamination of B. cereus in six types of commercially available plant foods from eight regions of a province. The molecular epidemiology of the isolated B. cereus was analyzed by whole-genome sequencing. We aimed to provide fundamental data for the surveillance and epidemiology analysis of B. cereus in food products in China. The rapid traceability system of B. cereus established in this study can provide a basis for rapid molecular epidemiology analysis of B. cereus, as well as for the prevention and surveillance of B. cereus. Moreover, it can also be expanded to monitoring and rapid tracing of more foodborne pathogens.

Genotyping-driven diversity assessment of biocontrol potent Bacillus spp. strain collection as a potential method for the development of strain-specific biomarkers

Bacillus species are among the most researched and frequently applied biocontrol agents. To estimate their potential as environmentally friendly microbial-based products, reliable and rapid plant colonization monitoring methods are essential. We evaluated repetitive element-based (rep) and Random Amplified Polymorphic DNA (RAPD) PCR (Polymerase Chain Reaction) genotyping in a diversity assessment of 251 strains from bulk soil, straw, and manure samples across Serbia, highlighting their discriminative force and the presence of unique bands. RAPD 272, OPG 5, and (GTG)₅ primers were most potent in revealing the high diversity of a sizable environmental Bacillus spp. collection. RAPD 272 also amplified a unique band for a proven biocontrol strain, opening the possibility of Sequence Characterized Amplified Region (SCAR) marker design. That will enable colonization studies using the SCAR marker for its specific detection. This study provides a guide for primer selection for diversity and monitoring studies of environmental Bacillus spp. isolates.

Rapid and Routine Molecular Typing Using Multiplex Polymerase Chain Reaction and MinION Sequencer

Molecular typing is an essential tool that has been extensively applied in laboratories as well as in clinical settings. Next-generation sequencing technologies promise high-throughput and cost-effective molecular applications; however, the accessibility of these technologies is limited due to the high capital cost. Oxford Nanopore Technologies (ONT) offers a MinION device with the advantages of real-time data analysis, rapid library preparation, and low cost per test. However, the advantages of the MinION device are often overshadowed by its lower raw accuracy. Herein, we present a concise multilocus sequence typing protocol of Staphylococcus aureus using multiplex polymerase chain reaction and Rapid Barcoding Kit for barcoding and MinION device for sequencing. Moreover, to clarify the effects of carryover DNA on tasks that require high sequence accuracy, we used the MinION flow cell in successive runs of washing and reusing. Our results revealed that the MinION flow cell could achieve accurate typing of a total of 467 samples with 3,269 kilobase-long genes within a total of 5 runs. This thus demonstrates the effectiveness of a portable nanopore MinION sequencer in providing accurate, rapid, and routine molecular typing.

Genetic analysis of plasmid-encoded mcr-1 resistance in Enterobacteriaceae derived from poultry meat in the Netherlands

Background

Colistin is classified as the highest priority and critically important antimicrobial for human medicine by WHO as it is the last resort agent for treatment of carbapenem-resistant Enterobacteriaceae in humans. Additional research is necessary to elucidate the genetic structure of mcr-1 resistance genes, commonly found on plasmids, using WGS.

Objectives

To map and compare the genetic characteristics of 35 mcr-1-mediated colistin-resistant Enterobacteriaceae isolated from chicken meat to highlight the genetic variation of the mcr-1-containing plasmids.

Methods

Sequencing was performed using Illumina HiSeq2500, Novaseq6000 and ONT’s GridION. GridION data was locally basecalled and demultiplexed using ONT’s Albacore 2.3.4 followed by Porechop 2.3. Quality filtering was performed using Filtlong 2.0. Hybrid Assembly was performed using Unicycler 4.7. Plasmids were compared with reference sequences in plasmid-RefSeq and pATLAS.

Results

A total of 35 mcr-1 positive Enterobacteriaceae were investigated, which resulted in 34 qualitatively robust hybrid assemblies of 2 Klebsiella pneumoniae and 32 Escherichia coli. mcr-1.1 was present in 33/34 isolates. One isolate contained an mcr-1.1-like resistance gene, due to a deletion of one codon. Two mcr-1.1 genes were located on the chromosome, while the majority of the mcr-1 genes were found on IncX4 type plasmids (n = 19). Almost all plasmids identified in this study were highly similar to plasmids found in human-derived strains.

Conclusions

The mcr-1.1-containing plasmids from retail chicken show high sequence similarity to human mcr-1.1 plasmids, suggesting that this may be a contributor to the presence of colistin resistance in humans.

Antimicrobial Resistance and Genomic Characterization of Six New Sequence Types in Multidrug-Resistant Pseudomonas aeruginosa Clinical Isolates from Pakistan

Pseudomonas aeruginosa (P. aeruginosa) is a major bacterial pathogen associated with a variety of infections with high mortality rates. Most of the clinical P. aeruginosa isolates belong to a limited number of genetic subgroups characterized by multiple housekeeping genes’ sequences (usually 5–7) through the Multi-Locus Sequence Typing (MLST) scheme. The emergence and dissemination of novel multidrug-resistant (MDR) sequence types (ST) in P. aeruginosa pose serious clinical concerns. We performed whole-genome sequencing on a cohort (n = 160) of MDR P. aeruginosa isolates collected from a tertiary care hospital lab in Pakistan and found six isolates belonging to six unique MLST allelic profiles. The genomes were submitted to the PubMLST database and new ST numbers (ST3493, ST3494, ST3472, ST3489, ST3491, and ST3492) were assigned to the respective allele combinations. MLST and core-genome-based phylogenetic analysis confirmed the divergence of these isolates and positioned them in separate branches. Analysis of the resistome of the new STs isolates revealed the presence of genes blaOXA-50, blaPAO, blaPDC, blaVIM-2, aph(3′)-IIb, aac(6′)-II, aac(3)-Id, fosA, catB7, dfrB2, crpP, merP and a number of missense and frame-shift mutations in chromosomal genes conferring resistance to various antipseudomonal antibiotics. The exoS, exoT, pvdE, rhlI, rhlR, lasA, lasB, lasI, and lasR genes were the most prevalent virulence-related genes among the new ST isolates. The different genotypic features revealed the adaptation of these new clones to a variety of infections by various mutations in genes affecting antimicrobial resistance, quorum sensing and biofilm formation. Close monitoring of these antibiotic-resistant pathogens and surveillance mechanisms needs to be adopted to reduce their spread to the healthcare facilities of Pakistan. We believe that these strains can be used as reference strains for future comparative analysis of isolates belonging to the same STs.

Application of simplified MLST scheme for direct typing of clinical samples from human leptospirosis cases in a tertiary hospital in the Philippines

Despite the major threat of leptospirosis to public health in the Philippines, its epidemiologic data remain scarce. Multilocus sequence typing (MLST) is a method often used for identification of circulating Leptospira species and disease surveillance. Unfortunately, molecular typing of Leptospira isolates is not routinely done in most hospital settings. A simplified MLST scheme targeting three loci (adk, lipL41, mreA) was performed for rapid direct typing of Leptospira in clinical specimens. Blood samples from suspected or clinically diagnosed cases (n = 50) were initially screened via polymerase chain reaction (PCR) targeting 23S rRNA, 16S rRNA (rrs2), and lipL32 genes. From the nine positives, seven had interpretable data from MLST. Allelic profiles identified L. interrogans in all positive samples. Six were assigned to ST12 of serovar Manilae (serogroup Pyrogenes) while one sample cannot be clearly differentiated between two serovars/serogroups, Bataviae/Losbanos (serogroup Bataviae) or Australis (serogroup Australis), indicating possibility of a new ST. Phylogenetic analysis confirmed that the application of simplified MLST scheme produces consistent results with the seven-loci genetic profile of published Leptospira MLST schemes. Reduced scheme addressed the challenges often encountered in the amplification of full MLST genetic profile of Leptospira. The approach is a potential alternative to serological tests for rapid typing of clinical specimens and can also aid in investigations on disease epidemiology specifically to monitor occurrence, pathogen transmission, host specificity and susceptibility, and other factors that could lead to potential outbreaks.

Whole-genome sequencing and ad hoc shared genome analysis of Staphylococcus aureus isolates from a New Zealand primary school

Epidemiological studies of communicable diseases increasingly use large whole-genome sequencing (WGS) datasets to explore the transmission of pathogens. It is important to obtain an initial overview of datasets and identify closely related isolates, but this can be challenging with large numbers of isolates and imperfect sequencing. We used an ad hoc whole-genome multi locus sequence typing method to summarise data from a longitudinal study of Staphylococcus aureus in a primary school in New Zealand. Each pair of isolates was compared and the number of genes where alleles differed between isolates was tallied to produce a matrix of "allelic differences". We plotted histograms of the number of allelic differences between isolates for: all isolate pairs; pairs of isolates from different individuals; and pairs of isolates from the same individual. 340 sequenced isolates were included, and the ad hoc shared genome contained 445 genes. There were between 0 and 420 allelic differences between isolate pairs and the majority of pairs had more than 260 allelic differences. We found many genetically closely related S. aureus isolates from single individuals and a smaller number of closely-related isolates from separate individuals. Multiple S. aureus isolates from the same individual were usually very closely related or identical over the ad hoc shared genome. Siblings carried genetically similar, but not identical isolates. An ad hoc shared genome approach to WGS analysis can accommodate imperfect sequencing of the included isolates, and can provide insights into relationships between isolates in epidemiological studies with large WGS datasets containing diverse isolates.

The global dissemination of bacterial infections necessitates the study of reverse genomic epidemiology

Whole genome sequencing (WGS) has revolutionized the genotyping of bacterial pathogens and is expected to become the new gold standard for tracing the transmissions of bacterial infectious diseases for public health purposes. Traditional genomic epidemiology often uses WGS as a verification tool, namely, when a common source or epidemiological link is suspected, the collected isolates are sequenced for the determination of clonal relationships. However, increasingly frequent international travel and food transportation, and the associated potential for the cross-border transmission of bacterial pathogens, often lead to an absence of information on bacterial transmission routes. Here we introduce the concept of 'reverse genomic epidemiology', i.e. when isolates are inspected by genome comparisons to be sufficiently similar to one another, they are assumed to be a consequence of infection from a common source. Through BacWGSTdb (http://bacdb.org/BacWGSTdb/), a database we have developed for bacterial genome typing and source tracking, we have found that almost the entire analyzed 20 bacterial species exhibit the phenomenon of cross-border clonal dissemination. Five networks were further identified in which isolates sharing nearly identical genomes were collected from at least five different countries. Three of these have been documented as real infectious disease outbreaks, therefore demonstrating the feasibility and authority of reverse genomic epidemiology. Our survey and proposed strategy would be of potential value in establishing a global surveillance system for tracing bacterial transmissions and outbreaks; the related database and techniques require urgent standardization.

The Application of High-Throughput Technologies for the Study of Microbiome and Cancer

Human gut microbiome research, especially gut microbiome, has been developing at a considerable pace over the last decades, driven by a rapid technological advancement. The emergence of high-throughput technologies, such as genomics, transcriptomics, and others, has afforded the generation of large volumes of data, and in relation to specific pathologies such as different cancer types. The current review identifies high-throughput technologies as they have been implemented in the study of microbiome and cancer. Four main thematic areas have emerged: the characterization of microbial diversity and composition, microbial functional analyses, biomarker prediction, and, lastly, potential therapeutic applications. The majority of studies identified focus on the microbiome diversity characterization, which is reaching technological maturity, while the remaining three thematic areas could be described as emerging.

Antibiotic resistance: Time of synthesis in a post-genomic age

Antibiotic resistance has been highlighted by international organizations, including World Health Organization, World Bank and United Nations, as one of the most relevant global health problems. Classical approaches to study this problem have focused in infected humans, mainly at hospitals. Nevertheless, antibiotic resistance can expand through different ecosystems and geographical allocations, hence constituting a One-Health, Global-Health problem, requiring specific integrative analytic tools. Antibiotic resistance evolution and transmission are multilayer, hierarchically organized processes with several elements (from genes to the whole microbiome) involved. However, their study has been traditionally gene-centric, each element independently studied. The development of robust-economically affordable whole genome sequencing approaches, as well as other -omic techniques as transcriptomics and proteomics, is changing this panorama. These technologies allow the description of a system, either a cell or a microbiome as a whole, overcoming the problems associated with gene-centric approaches. We are currently at the time of combining the information derived from -omic studies to have a more holistic view of the evolution and spread of antibiotic resistance. This synthesis process requires the accurate integration of -omic information into computational models that serve to analyse the causes and the consequences of acquiring AR, fed by curated databases capable of identifying the elements involved in the acquisition of resistance. In this review, we analyse the capacities and drawbacks of the tools that are currently in use for the global analysis of AR, aiming to identify the more useful targets for effective corrective interventions.

Characterization of genetic diversity and population structure within Staphylococcus chromogenes by multilocus sequence typing

Staphylococcus chromogenes is a common skin commensal in cattle and has been identified as a frequent cause of bovine mastitis and intramammary infections. We have developed a seven locus Multilocus Sequence Typing (MLST) scheme for typing S. chromogenes. Sequence-based typing systems, such as MLST, have application in studies of genetic diversity, population structure, and epidemiology, including studies of strain variation as a factor in pathogenicity or host adaptation. The S. chromogenes scheme was tested on 120 isolates collected from three geographic locations, Vermont and Washington State in the United States and Belgium. A total of 46 sequence types (STs) were identified with most of the STs being location specific. The utility of the typing scheme is indicated by a discrimination power of 95.6% for all isolates and greater than 90% for isolates from each of the three locations. Phylogenetic analysis placed 39 of the 46 STs into single core group consistent with a common genetic lineage; the STs in this group differ by less than 0.5% at the nucleotide sequence level. Most of the diversification in this lineage group can be attributed to mutation; recombination plays a limited role. This lineage group includes two clusters of single nucleotide variants in starburst configurations indicative of recent clonal expansion; nearly 50% of the isolates sampled in this study are in these two clusters. The remaining seven STs were set apart from the core group by having alleles with highly variable sequences at one or more loci. Recombination had a higher impact than mutation in the diversification of these outlier STs. Alleles with hypervariable sequences were detected at five of the seven loci used in the MLST scheme; the average sequence distances between the hypervariable alleles and the common core alleles ranged from 12 to 34 nucleotides. The extent of these sequence differences suggests the hypervariable alleles may be remnants of an ancestral genotype.

Amplicon-Based Next Generation Sequencing for Rapid Identification of Rickettsia and Ectoparasite Species from Entomological Surveillance in Thailand

BACKGROUND

Next generation sequencing (NGS) technology has been used for a wide range of epidemiological and surveillance studies. Here, we used amplicon-based NGS to species identify Rickettsia and their arthropod hosts from entomological surveillance.

METHODS

During 2015-2016, we screened 1825 samples of rodents and ectoparasites collected from rodents and domestic mammals (dog, cat, and cattle) across Thailand for Rickettsia. The citrate synthase gene was amplified to identify Rickettsia to species, while the Cytochrome Oxidase subunit I (COI) and subunit II (COII) genes were used as target genes for ectoparasite identification. All target gene amplicons were pooled for library preparation and sequenced with Illumina MiSeq platform.

RESULT

The highest percentage of Rickettsia DNA was observed in fleas collected from domestic animals (56%) predominantly dogs. Only a few samples of ticks from domestic animals, rodent fleas, and rodent tissue were positive for Rickettisia DNA. NGS based characterization of Rickettsia by host identified Rickettsia asembonensis as the most common bacteria in positive fleas collected from dogs (83.2%) while "Candidatus Rickettsia senegalensis" was detected in only 16.8% of Rickettsia positive dog fleas. Sequence analysis of COI and COII revealed that almost all fleas collected from dogs were Ctenocephalides felis orientis. Other Rickettsia species were detected by NGS including Rickettsia heilongjiangensis from two Haemaphysalis hystricis ticks, and Rickettsia typhi in two rodent tissue samples.

CONCLUSION

This study demonstrates the utility of NGS for high-throughput sequencing in the species characterization/identification of bacteria and ectoparasite for entomological surveillance of rickettsiae. A high percentage of C. f. orientis are positive for R. asembonensis. In addition, our findings indicate there is a risk of tick-borne Spotted Fever Group rickettsiosis, and flea-borne murine typhus transmission in Tak and Phangnga provinces of Thailand.

One is not enough: On the effects of reference genome for the mapping and subsequent analyses of short-reads

Mapping of high-throughput sequencing (HTS) reads to a single arbitrary reference genome is a frequently used approach in microbial genomics. However, the choice of a reference may represent a source of errors that may affect subsequent analyses such as the detection of single nucleotide polymorphisms (SNPs) and phylogenetic inference. In this work, we evaluated the effect of reference choice on short-read sequence data from five clinically and epidemiologically relevant bacteria (Klebsiella pneumoniae, Legionella pneumophila, Neisseria gonorrhoeae, Pseudomonas aeruginosa and Serratia marcescens). Publicly available whole-genome assemblies encompassing the genomic diversity of these species were selected as reference sequences, and read alignment statistics, SNP calling, recombination rates, dN/dS ratios, and phylogenetic trees were evaluated depending on the mapping reference. The choice of different reference genomes proved to have an impact on almost all the parameters considered in the five species. In addition, these biases had potential epidemiological implications such as including/excluding isolates of particular clades and the estimation of genetic distances. These findings suggest that the single reference approach might introduce systematic errors during mapping that affect subsequent analyses, particularly for data sets with isolates from genetically diverse backgrounds. In any case, exploring the effects of different references on the final conclusions is highly recommended.

Mapping consists in the alignment of reads (i.e., DNA fragments) obtained through high-throughput genome sequencing to a previously assembled reference sequence. It is a common practice in genomic studies to use a single reference for mapping, usually the ‘reference genome’ of a species—a high-quality assembly. However, the selection of an optimal reference is hindered by intrinsic intra-species genetic variability, particularly in bacteria. It is known that genetic differences between the reference genome and the read sequences may produce incorrect alignments during mapping. Eventually, these errors could lead to misidentification of variants and biased reconstruction of phylogenetic trees (which reflect ancestry between different bacterial lineages). To our knowledge, this is the first work to systematically examine the effect of different references for mapping on the inference of tree topology as well as the impact on recombination and natural selection inferences. Furthermore, the novelty of this work relies on a procedure that guarantees that we are evaluating only the effect of the reference. This effect has proved to be pervasive in the five bacterial species that we have studied and, in some cases, alterations in phylogenetic trees could lead to incorrect epidemiological inferences. Hence, the use of different reference genomes may be prescriptive to assess the potential biases of mapping.

Genome-Based Analysis of a Sequence Type 1049 Hypervirulent Klebsiella pneumoniae Causing Bacteremic Neck Abscess

Hypervirulent Klebsiella pneumoniae (hvKP) has raised grave concerns in recent years and can cause severe infections with diverse anatomic locations including liver abscess, meningitis, and endophthalmitis. However, there is limited data about neck abscess caused by hvKP. A K. pneumoniae strain Kp_whw was isolated from neck abscess. We characterized the genetic background, virulence determinates of the strain by genomic analysis and dertermined the virulence level by serum resistance assay. Kp_whw belonged to sequence type (ST) 1049 K locus (KL) 5. Kp_whw showed hypermucoviscosity phenotype and was resistant to ampicillin but susceptible to the majority of the other antimicrobial agents. A pLVPK-like virulence plasmid and a chromosomal ICEKp5-like mobile genetic element were carried by Kp_whw, resulting in the risk of dissemination of hypervirulence. The strain exhibited relative higher level of core genome allelic diversity than accessory genome profile, in comparison to hvKP of K1/K2 serotype. Kp_whw was finally demonstrated as virulent as the ST23 K1 serotype hvKP strain NTUH-K2044 in vitro. In conclusion, this work elaborates the genetic background of a clinical hvKP strain with an uncommon ST, reinforcing our understanding of virulence mechanisms of hvKP.

Multilevel genome typing: genomics-guided scalable resolution typing of microbial pathogens

Background

Both long- and short-term epidemiology are fundamental to disease control and require accurate bacterial typing. Genomic data resulting from implementation of whole genome sequencing in many public health laboratories can potentially provide highly sensitive and accurate descriptions of strain relatedness. Previous typing efforts using these data have mainly focussed on outbreak detection.

Aim

We aimed to develop multilevel genome typing (MGT), using consecutive multilocus sequence typing (MLST) schemes of increasing sizes, stepping up from seven-gene MLST to core genome MLST, to allow examination of genetic relatedness at multiple resolution levels.

Methods

The system was applied to Salmonellaenterica serovar Typhimurium. The MLST scheme used at each step (MGT level), defined a given MGT-level specific sequence type (ST). The list of STs generated from all of these increasing MGT levels, was named a genome type (GT). Using MGT, we typed 9,096 previously characterised isolates with publicly available data.

Results

Our approach could identify previously described S. Typhimurium populations, such as the DT104 multidrug resistance lineage (GT 19-2-11) and two invasive lineages of African isolates (GT 313-2-3 and 313-2-752). Further, we showed that MGT-derived clusters can accurately distinguish five outbreaks from each other and five background isolates.

Conclusion

MGT provides a universal and stable nomenclature at multiple resolutions for S. Typhimurium strains and could be implemented as an internationally standardised strain identification system. While established so far only for S. Typhimurium, the results here suggest that MGT could form the basis for typing systems in other similar microorganisms.

nanoMLST: accurate multilocus sequence typing using Oxford Nanopore Technologies MinION with a dual-barcode approach to multiplex large numbers of samples

Multilocus sequence typing (MLST) is one of the most commonly used methods for studying microbial lineage worldwide. However, the traditional MLST process using Sanger sequencing is time-consuming and expensive. We have designed a workflow that simultaneously sequenced seven full-length housekeeping genes of 96 meticillin-resistant Staphylococcus aureus isolates with dual-barcode multiplexing using just a single flow cell of an Oxford Nanopore Technologies MinION system, and then we performed bioinformatic analysis for strain typing. Fifty-one of the isolates comprising 34 sequence types had been characterized using Sanger sequencing. We demonstrate that the allele assignments obtained by our nanopore workflow (nanoMLST, available at https://github.com/jade-nhri/nanoMLST) were identical to those obtained by Sanger sequencing (359/359, with 100 % agreement rate). In addition, we estimate that our multiplex system is able to perform MLST for up to 1000 samples simultaneously; thus, providing a rapid and cost-effective solution for molecular typing.

The phylogenomics of CRISPR-Cas system and revelation of its features in Salmonella

Salmonellae display intricate evolutionary patterns comprising over 2500 serovars having diverse pathogenic profiles. The acquisition and/or exchange of various virulence factors influences the evolutionary framework. To gain insights into evolution of Salmonella in association with the CRISPR-Cas genes we performed phylogenetic surveillance across strains of 22 Salmonella serovars. The strains differed in their CRISPR1-leader and cas operon features assorting into two main clades, CRISPR1-STY/cas-STY and CRISPR1-STM/cas-STM, comprising majorly typhoidal and non-typhoidal Salmonella serovars respectively. Serovars of these two clades displayed better relatedness, concerning CRISPR1-leader and cas operon, across genera than between themselves. This signifies the acquisition of CRISPR1/Cas region could be through a horizontal gene transfer event owing to the presence of mobile genetic elements flanking CRISPR1 array. Comparison of CRISPR and cas phenograms with that of multilocus sequence typing (MLST) suggests differential evolution of CRISPR/Cas system. As opposed to broad-host-range, the host-specific serovars harbor fewer spacers. Mapping of protospacer sources suggested a partial correlation of spacer content with habitat diversity of the serovars. Some serovars like serovar Enteritidis and Typhimurium that inhabit similar environment/infect similar hosts hardly shared their protospacer sources.

Mini Review: Clinical Routine Microbiology in the Era of Automation and Digital Health

Clinical microbiology laboratories are the first line to combat and handle infectious diseases and antibiotic resistance, including newly emerging ones. Although most clinical laboratories still rely on conventional methods, a cascade of technological changes, driven by digital imaging and high-throughput sequencing, will revolutionize the management of clinical diagnostics for direct detection of bacteria and swift antimicrobial susceptibility testing. Importantly, such technological advancements occur in the golden age of machine learning where computers are no longer acting passively in data mining, but once trained, can also help physicians in making decisions for diagnostics and optimal treatment administration. The further potential of physically integrating new technologies in an automation chain, combined to machine-learning-based software for data analyses, is seducing and would indeed lead to a faster management in infectious diseases. However, if, from one side, technological advancement would achieve a better performance than conventional methods, on the other side, this evolution challenges clinicians in terms of data interpretation and impacts the entire hospital personnel organization and management. In this mini review, we discuss such technological achievements offering practical examples of their operability but also their limitations and potential issues that their implementation could rise in clinical microbiology laboratories.

Multilocus sequence typing analysis of Leishmania clinical isolates from cutaneous leishmaniasis patients of Iran

Cutaneous leishmaniasis (CL) is mainly caused by L. major and L. tropica in Old World and might be represented as typical skin lesion(s) or sometimes as a spectrum of atypical manifestations. We applied multilocus sequence typing (MLST) to explore genetic variations of Leishmania strains isolated from atypical vs. typical CL patients from Iran. A PCR-sequencing was performed for seven housekeeping genes (g6pd, mpi, asat, icd, 6pgd, fh, and trys) and genetic diversity indices and phylogenetic relationships were analyzed. A total of 41 isolates of L. major (28/41) and L. tropica (13/41) from 21 (51.2%) atypical CL and 20 (48.8%) typical CL cases were included. A set of additional sequences of 41 strains of 17 species of Leishmania were retrieved from databases. Different SNP variations were detected and the highest rate of heterozygous sites was found in g6pd and 6pgd genes (6 sites) for L. tropica and in asat and 6pgd genes (7 sites) for L. major strains. All strains were clustered into 58 unique sequence types (STs) including 17 STs related to 41 strains of Leishmania of this study. Concatenated tree clustered all strains in 6 main clades (A to F) including L. major (clade D) and L. tropica (clade B) strains. Two strains of L. major (codes 28 and 42) with highest nucleotide variations were more close to L. tropica and were grouped in Clade B. All of the STs were related in clonal complexes by using eBURST with the prediction of founder genotypes. A high rate of genetic variations and heterozygocity was evident in L. tropica and L. major strains; nevertheless, there was no significant difference in the diversity of Leishmania strains between typical CL and atypical CL groups. This study represents the first successful application of MLST approach to L. tropica and L. major strains in Iran.

Genomic Variation and Arsenic Tolerance Emerged as Niche Specific Adaptations by Different Exiguobacterium Strains Isolated From the Extreme Salar de Huasco Environment in Chilean - Altiplano

Polyextremophilic bacteria can thrive in environments with multiple stressors such as the Salar de Huasco (SH). Microbial communities in SH are exposed to low atmospheric pressure, high UV radiation, wide temperature ranges, salinity gradient and the presence of toxic compounds such as arsenic (As). In this work we focus on arsenic stress as one of the main adverse factors in SH and bacteria that belong to the Exiguobacterium genus due to their plasticity and ubiquity. Therefore, our aim was to shed light on the effect of niche conditions pressure (particularly arsenic), on the adaptation and divergence (at genotypic and phenotypic levels) of Exiguobacterium strains from five different SH sites. Also, to capture greater diversity in this genus, we use as outgroup five As(III) sensitive strains isolated from Easter Island (Chile) and The Great Salt Lake (United States). For this, samples were obtained from five different SH sites under an arsenic gradient (9 to 321 mg/kg: sediment) and isolated and sequenced the genomes of 14 Exiguobacterium strains, which had different arsenic tolerance levels. Then, we used comparative genomic analysis to assess the genomic divergence of these strains and their association with phenotypic differences such as arsenic tolerance levels and the ability to resist poly-stress. Phylogenetic analysis showed that SH strains share a common ancestor. Consequently, populations were separated and structured in different SH microenvironments, giving rise to multiple coexisting lineages. Hence, this genotypic variability is also evidenced by the COG (Clusters of Orthologous Groups) composition and the size of their accessory genomes. Interestingly, these observations correlate with physiological traits such as growth patterns, gene expression, and enzyme activity related to arsenic response and/or tolerance. Therefore, Exiguobacterium strains from SH are adapted to physiologically overcome the contrasting environmental conditions, like the arsenic present in their habitat.

Estimates of the reproductive numbers and demographic reconstruction of outbreak associated with C:P1.5-1,10-8:F3-6:ST-11(cc11) Neisseria meningitidis strains

BACKGROUND

Neisseria meningitidis can cause sporadic cases and outbreaks of invasive disease, including meningitis and sepsis. The meningococcal serogroup C (MenC) is the second most common serogroup in Italy after MenB. In this study we have estimated the reproductive numbers and the demographic reconstruction on the genomes of invasive N. meningitidis C:P1.5-1,10-8:F3-6:ST-11(cc11) strains isolated in Italy in 2012 - 2017, a period that includes the outbreak in Tuscany.

METHODS

The genomes of N. meningitidis were sequenced using the Illumina MiSeq platform, through the whole genome sequencing (WGS) method and were analyzed by the core genome MLST (cgMLST) approach, using the BIGSdb Genome Comparator tool implemented on the PubMLST website. A Bayesian method was applied to study population dynamics across the entire N. meningitidis dataset. The basic reproduction number R_0, which indicates the average number of secondary cases generated by a single primary case, was calculated using a Bayesian method, on the dataset and on the two subsets. The effective reproduction number R(t), defined as the average number of secondary cases per infectious case in a population, made up of susceptible and non-susceptible hosts was studied on the Tuscany dataset, with a Bayesian method.

RESULTS

An increase in the effective number of the N. meningitidis infections was observed between 2013 and 2016. The estimated R₀ parameter was 1.31 (95% HPD: 1.03 - 1.64), 1.22 (95% HPD: 0.90 - 1.64) and 1.4 (95% HPD: 0.91 - 1.9) for the entire dataset, first and second subset, respectively. The BDSKY estimated an initial R(t) of about 2.0 (95% HPD: 0.04 - 5.0), which showed a growing trend at the end of 2014, reaching an average value of 3.22 in 2015, and then declining below 1 from the year 2016.

CONCLUSION

Monitoring the effective reproduction number can help to inform future vaccination strategies. The increase in the reproductive number for the Tuscany dataset, was consistent with the amplification event that led to the Tuscany outbreak. Subsequently, the intervention that led to the decline of the cases was followed, suggesting a high effectiveness of the vaccination campaign.

Genomics of Serratia marcescens Isolates Causing Outbreaks in the Same Pediatric Unit 47 Years Apart: Position in an Updated Phylogeny of the Species

The first documented nosocomial outbreak caused by Serratia marcescens in Spain occurred in 1969 at the neonatal intensive care unit (NICU) of the tertiary La Paz Children's Hospital in Madrid, Spain, and based on the available phenotyping techniques at this time, it was considered as a monoclonal outbreak. Only 47 years later, another S. marcescens outbreak of an equivalent dimension occurred at the same NICU. The aim of the present study was to study isolates from these historical and contemporary outbreaks by phenotypic analysis and whole-genome sequencing techniques and to position these strains along with 444 publicly available S. marcescens genomes, separately comparing core genome and accessory genome contents. Clades inferred by both approaches showed high correlation, indicating that core and accessory genomes seem to evolve in the same manner for S. marcescens. Nine S. marcescens clusters were identified, and isolates were grouped in two of them according to sampling year. One exception was isolate 13F-69, the most genetically distant strain, located in a different cluster. Categorical functions in the annotated accessory genes of both collections were preserved among all isolates. No significant differences in frequency of insertion sequences in historical (0.18-0.20)-excluding the outlier strain-versus contemporary isolates (0.11-0.19) were found despite the expected resting effect. The most dissimilar isolate, 13F-69, contains a highly preserved plasmid previously described in Bordetella bronchiseptica. This strain exhibited a few antibiotic resistance genes not resulting in a resistant phenotype, suggesting the value of gene down expression in adaptation to long-term starvation.

Molecular typing of Cyclospora cayetanensis in produce and clinical samples using targeted enrichment of complete mitochondrial genomes and next-generation sequencing

BACKGROUND

Outbreaks of cyclosporiasis, a diarrheal illness caused by Cyclospora cayetanensis, have been a public health issue in the USA since the mid 1990's. In 2018, 2299 domestically acquired cases of cyclosporiasis were reported in the USA as a result of multiple large outbreaks linked to different fresh produce commodities. Outbreak investigations are hindered by the absence of standardized molecular epidemiological tools for C. cayetanensis. For other apicomplexan coccidian parasites, multicopy organellar DNA such as mitochondrial genomes have been used for detection and molecular typing.

METHODS

We developed a workflow to obtain complete mitochondrial genome sequences from cilantro samples and clinical samples for typing of C. cayetanensis isolates. The 6.3 kb long C. cayetanensis mitochondrial genome was amplified by PCR in four overlapping amplicons from genomic DNA extracted from cilantro, seeded with oocysts, and from stool samples positive for C. cayetanensis by diagnostic methods. DNA sequence libraries of pooled amplicons were prepared and sequenced via next-generation sequencing (NGS). Sequence reads were assembled using a custom bioinformatics pipeline.

RESULTS

This approach allowed us to sequence complete mitochondrial genomes from the samples studied. Sequence alterations, such as single nucleotide polymorphism (SNP) profiles and insertion and deletions (InDels), in mitochondrial genomes of 24 stool samples from patients with cyclosporiasis diagnosed in 2014, exhibited discriminatory power. The cluster dendrogram that was created based on distance matrices of the complete mitochondrial genome sequences, indicated distinct strain-level diversity among the 2014 C. cayetanensis outbreak isolates analyzed in this study.

CONCLUSIONS

Our results suggest that genomic analyses of mitochondrial genome sequences may help to link outbreak cases to the source.

Leptospira spp., a genus in the stage of diversity and genomic data expansion

Leptospirosis is a widespread global zoonotic bacterial disease with a noteworthy human-animal-ecosystem interface. The disease presents different clinical manifestations and a high mortality and morbidity rates in humans and animals throughout the world. Characterization and correct classification of Leptospira isolates is essential for a better understanding the epidemiological properties of the disease. In the last ten years, molecular typing tools have been developed and applied to this field. These methods together with the availability of hundreds of new whole genome sequences that belong to known and new described species are shaping the understanding and structure of the entire genus.

A cross-sectional study to characterize local HIV-1 dynamics in Washington, DC using next-generation sequencing

Washington, DC continues to experience a generalized HIV-1 epidemic. We characterized the local phylodynamics of HIV-1 in DC using next-generation sequencing (NGS) data. Viral samples from 68 participants from 2016 through 2017 were sequenced and paired with epidemiological data. Phylogenetic and network inferences, drug resistant mutations (DRMs), subtypes and HIV-1 diversity estimations were completed. Haplotypes were reconstructed to infer transmission clusters. Phylodynamic inferences based on the HIV-1 polymerase (pol) and envelope genes (env) were compared. Higher HIV-1 diversity (n.s.) was seen in men who have sex with men, heterosexual, and male participants in DC. 54.0% of the participants contained at least one DRM. The 40-49 year-olds showed the highest prevalence of DRMs (22.9%). Phylogenetic analysis of pol and env sequences grouped 31.9-33.8% of the participants into clusters. HIV-TRACE grouped 2.9-12.8% of participants when using consensus sequences and 9.0-64.2% when using haplotypes. NGS allowed us to characterize the local phylodynamics of HIV-1 in DC more broadly and accurately, given a better representation of its diversity and dynamics. Reconstructed haplotypes provided novel and deeper phylodynamic insights, which led to networks linking a higher number of participants. Our understanding of the HIV-1 epidemic was expanded with the powerful coupling of HIV-1 NGS data with epidemiological data.

Escherichia coli strains producing a novel Shiga toxin 2 subtype circulate in China

Shiga toxin (Stx) is the key virulence factor in Shiga toxin producing Escherichia coli (STEC), which can cause diarrhea and hemorrhagic colitis with life-threatening complications. Stx comprises two toxin types, Stx1 and Stx2. Several Stx1/Stx2 subtypes have been identified in E. coli, which are variable in sequences, toxicity and host specificity. Here, we report the identification of a novel Stx2 subtype, designated Stx2k, in E. coli strains widely detected from diarrheal patients, animals, and raw meats in China over time. Stx2k exhibits varied cytotoxicity in vitro among individual strains. The Stx2k converting prophages displayed considerable heterogeneity in terms of insertion site, genetic content and structure. Whole genome analysis revealed that the stx_2k-containing strains were genetically heterogeneous with diverse serotypes, sequence types, and virulence gene profiles. The nine stx_2k-containing strains formed two major phylogenetic clusters closely with strains belonging to STEC, enterotoxigenic E. coli (ETEC), and STEC/ETEC hybrid. One stx_2k-containing strain harbored one plasmid-encoded heat-stable enterotoxin sta gene and two identical copies of chromosome-encoded stb gene, exhibiting STEC/ETEC hybrid pathotype. Our finding enlarges the pool of Stx2 subtypes and highlights the extraordinary genomic plasticity of STEC strains. Given the wide distribution of the Stx2k-producing strains in diverse sources and their pathogenic potential, Stx2k should be taken into account in epidemiological surveillance of STEC infections and clinical diagnosis.

Application of mini-MLST and whole genome sequencing in low diversity hospital extended-spectrum beta-lactamase producing Klebsiella pneumoniae population

Studying bacterial population diversity is important to understand healthcare associated infections’ epidemiology and has a significant impact on dealing with multidrug resistant bacterial outbreaks. We characterised the extended-spectrum beta-lactamase producing K. pneumoniae (ESBLp KPN) population in our hospital using mini-MLST. Then we used whole genome sequencing (WGS) to compare selected isolates belonging to the most prevalent melting types (MelTs) and the colonization/infection pair isolates collected from one patient to study the ESBLp KPN population’s genetic diversity. A total of 922 ESBLp KPN isolates collected between 7/2016 and 5/2018 were divided into 38 MelTs using mini-MLST with only 6 MelTs forming 82.8% of all isolates. For WGS, 14 isolates from the most prominent MelTs collected in the monitored period and 10 isolates belonging to the same MelTs collected in our hospital in 2014 were randomly selected. Resistome, virulome and ST were MelT specific and stable over time. A maximum of 23 SNV per core genome and 58 SNV per core and accessory genome were found. To determine the SNV relatedness cut-off values, 22 isolates representing colonization/infection pair samples obtained from 11 different patients were analysed by WGS with a maximum of 22 SNV in the core genome and 40 SNV in the core and accessory genome within pairs. The mini-MLST showed its potential for real-time epidemiology in clinical practice. However, for outbreak evaluation in a low diversity bacterial population, mini-MLST should be combined with more sensitive methods like WGS. Our findings showed there were only minimal differences within the core and accessory genome in the low diversity hospital population and gene based SNV analysis does not have enough discriminatory power to differentiate isolate relatedness. Thus, intergenic regions and mobile elements should be incorporated into the analysis scheme to increase discriminatory power.

Update on the epidemiological typing methods for Acinetobacter baumannii

The outstanding ability of Acinetobacter baumannii to cause outbreaks and acquire multidrug resistance motivated the development of a plethora of typing techniques, which can help infection preventionists and hospital epidemiologists to more efficiently implement intervention controls. Nowadays, the world is witnessing a gradual transition from traditional typing methodology to whole genome sequencing-based approaches. Such approaches are opening new prospects and applications never achieved by existing typing methods. Herein, we provide the reader with an updated review on A. baumannii typing methods recapping the added value of well-established techniques previously applied for A. baumannii and detailing new ones (as clustered regularly interspaced short palindromic repeats-based typing) with a special focus on whole genome sequencing.

Operational models and criteria for incorporating microbial whole genome sequencing in hospital microbiology - A systematic literature review

BACKGROUND

Microbial whole genome sequencing (WGS) has many advantages over standard microbiological methods. However, it is not yet widely implemented in routine hospital diagnostics due to notable challenges.

OBJECTIVES

The aim was to extract managerial, financial and clinical criteria supporting the decision to implement WGS in routine diagnostic microbiology, across different operational models of implementation in the hospital setting.

METHODS

This was a systematic review of literature identified through PubMed and Web of Science. English literature studies discussing the applications of microbial WGS without limitation on publication date were eligible. A narrative approach for categorization and synthesis of the sources identified was adopted.

RESULTS

A total of 98 sources were included. Four main alternative operational models for incorporating WGS in clinical microbiology laboratories were identified: full in-house sequencing and analysis, full outsourcing of sequencing and analysis and two hybrid models combining in-house/outsourcing of the sequencing and analysis components. Six main criteria (and multiple related sub-criteria) for WGS implementation emerged from our review and included cost (e.g. the availability of resources for capital and operational investment); manpower (e.g. the ability to provide training programmes or recruit trained personnel), laboratory infrastructure (e.g. the availability of supplies and consumables or sequencing platforms), bioinformatics requirements (e.g. the availability of valid analysis tools); computational infrastructure (e.g. the availability of storage space or data safety arrangements); and quality control (e.g. the existence of standardized procedures).

CONCLUSIONS

The decision to incorporate WGS in routine diagnostics involves multiple, sometimes competing, criteria and sub-criteria. Mapping these criteria systematically is an essential stage in developing policies for adoption of this technology, e.g. using a multicriteria decision tool. Future research that will prioritize criteria and sub-criteria that were identified in our review in the context of operational models will inform decision-making at clinical and managerial levels with respect to effective implementation of WGS for routine use. Beyond WGS, similar decision-making challenges are expected with respect to future integration of clinical metagenomics.

Enterohaemorrhagic and other Shiga toxin-producing Escherichia coli (STEC): Where are we now regarding diagnostics and control strategies?

Escherichia coli comprises a highly diverse group of Gram-negative bacteria and is a common member of the intestinal microflora of humans and animals. Generally, such colonization is asymptomatic; however, some E. coli strains have evolved to become pathogenic and thus cause clinical disease in susceptible hosts. One pathotype, the Shiga toxigenic E. coli (STEC) comprising strains expressing a Shiga-like toxin is an important foodborne pathogen. A subset of STEC are the enterohaemorrhagic E. coli (EHEC), which can cause serious human disease, including haemolytic uraemic syndrome (HUS). The diagnosis of EHEC infections and the surveillance of STEC in the food chain and the environment require accurate, cost-effective and timely tests. In this review, we describe and evaluate tests now in routine use, as well as upcoming test technologies for pathogen detection, including loop-mediated isothermal amplification (LAMP) and whole-genome sequencing (WGS). We have considered the need for improved diagnostic tools in current strategies for the control and prevention of these pathogens in humans, the food chain and the environment. We conclude that although significant progress has been made, STEC still remains an important zoonotic issue worldwide. Substantial reductions in the public health burden due to this infection will require a multipronged approach, including ongoing surveillance with high-resolution diagnostic techniques currently being developed and integrated into the routine investigations of public health laboratories. However, additional research requirements may be needed before such high-resolution diagnostic tools can be used to enable the development of appropriate interventions, such as vaccines and decontamination strategies.