Transforming microbial genotyping: a robotic pipeline for genotyping bacterial strains

Physical Laboratory Automation in Synthetic Biology

Synthetic Biology has overcome many of the early challenges facing the field and is entering a systems era characterized by adoption of Design-Build-Test-Learn (DBTL) approaches. The need for automation and standardization to enable reproducible, scalable, and translatable research has become increasingly accepted in recent years, and many of the hardware and software tools needed to address these challenges are now in place or under development. However, the lack of connectivity between DBTL modules and barriers to access and adoption remain significant challenges to realizing the full potential of lab automation. In this review, we characterize and classify the state of automation in synthetic biology with a focus on the physical automation of experimental workflows. Though fully autonomous scientific discovery is likely a long way off, impressive progress has been made toward automating critical elements of experimentation by combining intelligent hardware and software tools. It is worth questioning whether total automation that removes humans entirely from the loop should be the ultimate goal, and considerations for appropriate automation versus total automation are discussed in this light while emphasizing areas where further development is needed in both contexts.

SCRAMBLER: A Tool for De Novo CRISPR Array Reconstruction and Its Application for Analysis of the Structure of Prokaryotic Populations

CRISPR arrays are prokaryotic genomic loci consisting of repeat sequences alternating with unique spacers acquired from foreign nucleic acids. As one of the fastest-evolving parts of the genome, CRISPR arrays can be used to differentiate closely related prokaryotic lineages and track individual strains in prokaryotic communities. However, the assembly of full-length CRISPR arrays sequences remains a problem. Here, we developed SCRAMBLER, a tool that includes several pipelines for assembling CRISPR arrays from high-throughput short-read sequencing data. We assessed its performance with model data sets (Escherichia coli strains containing different CRISPR arrays and imitating prokaryotic communities of different complexities) and intestinal microbiomes of extant and extinct pachyderms. Evaluation of SCRAMBLER's performance using model data sets demonstrated its ability to assemble CRISPR arrays correctly from reads containing pairs of spacers, yielding a precision rate of >80% and a recall rate of 60-85% when checked against ground-truth data. Likewise, SCRAMBLER successfully assembled CRISPR arrays from the environmental samples, as attested by their matching with database entries. SCRAMBLER, an open-source software (github.com/biolab-tools/SCRAMBLER), can facilitate analysis of the composition and dynamics of CRISPR arrays in complex communities.

Genomic population structure associated with repeated escape of Salmonella enterica ATCC14028s from the laboratory into nature

Salmonella enterica serovar Typhimurium strain ATCC14028s is commercially available from multiple national type culture collections, and has been widely used since 1960 for quality control of growth media and experiments on fitness (“laboratory evolution”). ATCC14028s has been implicated in multiple cross-contaminations in the laboratory, and has also caused multiple laboratory infections and one known attempt at bioterrorism. According to hierarchical clustering of 3002 core gene sequences, ATCC14028s belongs to HierCC cluster HC20_373 in which most internal branch lengths are only one to three SNPs long. Many natural Typhimurium isolates from humans, domesticated animals and the environment also belong to HC20_373, and their core genomes are almost indistinguishable from those of laboratory strains. These natural isolates have infected humans in Ireland and Taiwan for decades, and are common in the British Isles as well as the Americas. The isolation history of some of the natural isolates confirms the conclusion that they do not represent recent contamination by the laboratory strain, and 10% carry plasmids or bacteriophages which have been acquired in nature by HGT from unrelated bacteria. We propose that ATCC14028s has repeatedly escaped from the laboratory environment into nature via laboratory accidents or infections, but the escaped micro-lineages have only a limited life span. As a result, there is a genetic gap separating HC20_373 from its closest natural relatives due to a divergence between them in the late 19^th century followed by repeated extinction events of escaped HC20_373.

Clades of closely related bacteria exist in nature. Individual isolates from such clades are often distinguishable by genomic sequencing because genomic sequence differences can be acquired over a few years due to neutral drift and natural selection. The evolution of laboratory strains is often largely frozen, physically due to storage conditions and genetically due to long periods of storage. Thus, laboratory strains can normally be readily distinguished from natural isolates because they show much less diversity. However, laboratory strain ATCC14028s shows modest levels of sequence diversity because it has been shipped around the world to multiple laboratories and is routinely used for analyses of laboratory evolution. Closely related natural isolates also exist, but their genetic diversity is not dramatically greater at the core genome level. Indeed, many scientists doubt that such isolates are natural, and interpret them as undetected contamination by the laboratory strain. We present data indicating that ATCC14028s has repeatedly escaped from the laboratory through inadvertent contamination of the environment, infection of technical staff and deliberate bioterrorism. The escapees survive in nature long enough that some acquire mobile genomic elements by horizontal gene transfer, but eventually they go extinct. As a result, even extensive global databases of natural isolates lack closely related isolates whose ancestors diverged from ATCC14028s within the last 100 years.

Genomic diversity of Salmonella enterica -The UoWUCC 10K genomes project

Background: Most publicly available genomes of Salmonella enterica are from human disease in the US and the UK, or from domesticated animals in the US. Methods: Here we describe a historical collection of 10,000 strains isolated between 1891-2010 in 73 different countries. They encompass a broad range of sources, ranging from rivers through reptiles to the diversity of all S. enterica isolated on the island of Ireland between 2000 and 2005. Genomic DNA was isolated, and sequenced by Illumina short read sequencing. Results: The short reads are publicly available in the Short Reads Archive. They were also uploaded to EnteroBase, which assembled and annotated draft genomes. 9769 draft genomes which passed quality control were genotyped with multiple levels of multilocus sequence typing, and used to predict serovars. Genomes were assigned to hierarchical clusters on the basis of numbers of pair-wise allelic differences in core genes, which were mapped to genetic Lineages within phylogenetic trees. Conclusions: The University of Warwick/University College Cork (UoWUCC) project greatly extends the geographic sources, dates and core genomic diversity of publicly available S. enterica genomes. We illustrate these features by an overview of core genomic Lineages within 33,000 publicly available Salmonella genomes whose strains were isolated before 2011. We also present detailed examinations of HC400, HC900 and HC2000 hierarchical clusters within exemplar Lineages, including serovars Typhimurium, Enteritidis and Mbandaka. These analyses confirm the polyphyletic nature of multiple serovars while showing that discrete clusters with geographical specificity can be reliably recognized by hierarchical clustering approaches. The results also demonstrate that the genomes sequenced here provide an important counterbalance to the sampling bias which is so dominant in current genomic sequencing.

Genomic diversity of Salmonella enterica -The UoWUCC 10K genomes project

Background: Most publicly available genomes of Salmonella enterica are from human disease in the US and the UK, or from domesticated animals in the US. Methods: Here we describe a historical collection of 10,000 strains isolated between 1891-2010 in 73 different countries. They encompass a broad range of sources, ranging from rivers through reptiles to the diversity of all S. enterica isolated on the island of Ireland between 2000 and 2005. Genomic DNA was isolated, and sequenced by Illumina short read sequencing. Results: The short reads are publicly available in the Short Reads Archive. They were also uploaded to EnteroBase, which assembled and annotated draft genomes. 9769 draft genomes which passed quality control were genotyped with multiple levels of multilocus sequence typing, and used to predict serovars. Genomes were assigned to hierarchical clusters on the basis of numbers of pair-wise allelic differences in core genes, which were mapped to genetic Lineages within phylogenetic trees. Conclusions: The University of Warwick/University College Cork (UoWUCC) project greatly extends the geographic sources, dates and core genomic diversity of publicly available S. enterica genomes. We illustrate these features by an overview of core genomic Lineages within 33,000 publicly available Salmonella genomes whose strains were isolated before 2011. We also present detailed examinations of HC400, HC900 and HC2000 hierarchical clusters within exemplar Lineages, including serovars Typhimurium, Enteritidis and Mbandaka. These analyses confirm the polyphyletic nature of multiple serovars while showing that discrete clusters with geographical specificity can be reliably recognized by hierarchical clustering approaches. The results also demonstrate that the genomes sequenced here provide an important counterbalance to the sampling bias which is so dominant in current genomic sequencing.

The EnteroBase user's guide, with case studies on Salmonella transmissions, Yersinia pestis phylogeny, and Escherichia core genomic diversity

EnteroBase is an integrated software environment that supports the identification of global population structures within several bacterial genera that include pathogens. Here, we provide an overview of how EnteroBase works, what it can do, and its future prospects. EnteroBase has currently assembled more than 300,000 genomes from Illumina short reads from Salmonella, Escherichia, Yersinia, Clostridioides, Helicobacter, Vibrio, and Moraxella and genotyped those assemblies by core genome multilocus sequence typing (cgMLST). Hierarchical clustering of cgMLST sequence types allows mapping a new bacterial strain to predefined population structures at multiple levels of resolution within a few hours after uploading its short reads. Case Study 1 illustrates this process for local transmissions of Salmonella enterica serovar Agama between neighboring social groups of badgers and humans. EnteroBase also supports single nucleotide polymorphism (SNP) calls from both genomic assemblies and after extraction from metagenomic sequences, as illustrated by Case Study 2 which summarizes the microevolution of Yersinia pestis over the last 5000 years of pandemic plague. EnteroBase can also provide a global overview of the genomic diversity within an entire genus, as illustrated by Case Study 3, which presents a novel, global overview of the population structure of all of the species, subspecies, and clades within Escherichia.

WGS based study of the population structure of Salmonella enterica serovar Infantis

BACKGROUND

Salmonella Infantis (S. Infantis) is one of the most frequent Salmonella serovars isolated from human cases of salmonellosis and the most detected serovar from animal and food sources in Europe. The serovar is commonly associated with poultry and there is increasing concern over multidrug resistant clones spreading worldwide, as the dominating clones are characterized by presence of large plasmids carrying multiple resistance genes. Increasing the knowledge of the S. Infantis population and evolution is important for understanding and preventing further spread. In this study, we analysed a collection of strains representing different decades, sources and geographic locations. We analysed the population structure and the accessory genome, in particular we identified prophages with a view to understand the role of prophages in relation to the evolution of this serovar.

RESULTS

We sequenced a global collection of 100 S. Infantis strains. A core-genome SNP analysis separated five strains in e-Burst Group (eBG) 297 with a long branch. The remaining strains, all in eBG31, were divided into three lineages that were estimated to have separated approximately 150 years ago. One lineage contained the vast majority of strains. In five of six clusters, no obvious correlation with source or geographical locations was seen. However, one cluster contained mostly strains from human and avian sources, indicating a clone with preference for these sources. The majority of strains within this cluster harboured a pESI-like plasmid with multiple resistance genes. Another lineage contained three genetic clusters with more rarely isolated strains of mainly animal origin, possibly less sampled or less infectious clones. Conserved prophages were identified in all strains, likely representing bacteriophages which integrated into the chromosome of a common ancestor to S. Infantis. We also saw that some prophages were specific to clusters and were probably introduced when the clusters were formed.

CONCLUSIONS

This study analysed a global S. Infantis population and described its genetic structure. We hypothesize that the population has evolved in three separate lineages, with one more successfully emerging lineage. We furthermore detected conserved prophages present in the entire population and cluster specific prophages, which probably shaped the population structure.

A genomic overview of the population structure of Salmonella

For many decades, Salmonella enterica has been subdivided by serological properties into serovars or further subdivided for epidemiological tracing by a variety of diagnostic tests with higher resolution. Recently, it has been proposed that so-called eBurst groups (eBGs) based on the alleles of seven housekeeping genes (legacy multilocus sequence typing [MLST]) corresponded to natural populations and could replace serotyping. However, this approach lacks the resolution needed for epidemiological tracing and the existence of natural populations had not been independently validated by independent criteria. Here, we describe EnteroBase, a web-based platform that assembles draft genomes from Illumina short reads in the public domain or that are uploaded by users. EnteroBase implements legacy MLST as well as ribosomal gene MLST (rMLST), core genome MLST (cgMLST), and whole genome MLST (wgMLST) and currently contains over 100,000 assembled genomes from Salmonella. It also provides graphical tools for visual interrogation of these genotypes and those based on core single nucleotide polymorphisms (SNPs). eBGs based on legacy MLST are largely consistent with eBGs based on rMLST, thus demonstrating that these correspond to natural populations. rMLST also facilitated the selection of representative genotypes for SNP analyses of the entire breadth of diversity within Salmonella. In contrast, cgMLST provides the resolution needed for epidemiological investigations. These observations show that genomic genotyping, with the assistance of EnteroBase, can be applied at all levels of diversity within the Salmonella genus.

Characterization of Invasive Salmonella Serogroup C1 Infections in Mali

Nontyphoidal Salmonella (NTS) are the leading cause of foodborne infections worldwide and a major cause of bloodstream infections in infants and HIV-infected adults in sub-Saharan Africa (SSA). Salmonella Typhimurium (serogroup B) and Salmonella Enteritidis (serogroup D) are the most common serovars in this region. However, data describing rarer invasive NTS serovars, particularly those belonging to serogroups C1 and C2, circulating in SSA are lacking. We previously conducted systematic blood culture surveillance on pediatric patients in Bamako, Mali, from 2002 to 2014, and the results showed that serovars Typhimurium and Enteritidis accounted for 32% and 36% of isolates, respectively. Here, we present data on 27 Salmonella serogroup C1 strains that were isolated during this previous study. The strains were typed by serum agglutination and multilocus sequence typing (MLST). Sixteen strains were Salmonella Paratyphi C, four were Salmonella Colindale, and two were Salmonella Virchow. Interestingly, five strains were identified as the very rare Salmonella Brazzaville using a combination of serum agglutination and flagellin gene typing. Phenotypic characterization showed that Salmonella Brazzaville produced biofilm and exhibited catalase activity, which were not statistically different from the gastroenteritis-associated Salmonella Typhimurium sequence type (ST) 19. All tested Salmonella Paratyphi C strains were poor biofilm producers and showed significantly less catalase activity than Salmonella Typhimurium ST19. Overall, our study provides insight into the Salmonella serogroup C1 serovars that cause invasive disease in infants in Mali. In addition, we show that MLST and flagellin gene sequencing, in association with traditional serum agglutination, are invaluable tools to help identify rare Salmonella serovars.

Incidence, Antimicrobial Resistance, and Molecular Characteristics of Nontyphoidal Salmonella Including Extended-Spectrum β-Lactamase Producers in Retail Chicken Meat

The present study was undertaken to determine the prevalence of Salmonella in 100 chicken carcass samples from five integrated broiler operation brands in Korea. Serotypes, antibiotic resistance patterns, extended-spectrum β-lactamase (ESBL) genotype, and clonal divergence using multilocus sequence typing of the isolated strains were analyzed. A total of 42 chicken samples were contaminated with nontyphoidal Salmonella (NTS) isolates: 16 isolates (38%) were Salmonella Virchow, 9 (21%) were Salmonella Bareilly, and 8 (19%) were Salmonella Infantis. A multidrug resistance (MDR; resistant to more than three classes of antibiotics) phenotype was observed in 29% of the isolates, which were resistant to five or more classes of antibiotics. The dominant MDR type was resistance to classes of penicillin, cephalosporins, aminoglycosides, quinolones, and tetracyclines. All the MDR isolates were positive for ESBL producers, and all but one (with the CTX-M-1 genotype) had the CTX-M-15 genotype. Multilocus sequence typing of the isolates revealed ST16 as the dominant sequence type; Salmonella Virchow, Salmonella Infantis, and Salmonella Richmond were all ST16, indicating a close genetic relationship between these serovars. This is the first study in Korea showing the CTX-M-1 type of NTS and the prevalence of ESBL-producing strains among NTS isolated from retail chicken meat. Our findings suggest that MDR Salmonella contamination is widely prevalent in retail chicken meat, and consumption of inadequately cooked products could lead to dissemination of NTS, which is hazardous to human health.

Mammalian Genotyping Using Acoustic Droplet Ejection for Enhanced Data Reproducibility, Superior Throughput, and Minimized Cross-Contamination

Genetically engineered animal models are major tools of a drug discovery pipeline because they facilitate understanding of the molecular and biochemical basis of disease. These highly complex models of human disease often require increasingly convoluted genetic analysis. With growing needs for throughput and consistency, we find that traditional aspiration-and-dispense liquid-handling robots no longer have the required speed, quality, or reproducibility.We present an adaptation and installation of an acoustic droplet ejection (ADE) liquid-handling system for ultra-high-throughput screening of genetically engineered models. An ADE system is fully integrated with existing laboratory processes and platforms to facilitate execution of PCR and quantitative PCR (qPCR) reactions. Such a configuration permits interrogation of highly complex genetic models in a variety of backgrounds. Our findings demonstrate that a single ADE system replaces 8-10 traditional liquid-handling robots while increasing quality and reproducibility.We demonstrate significant improvements achieved by transitioning to an ADE device: extremely low detectable cross-contamination in PCR and qPCR despite extensive use, greatly increased data reproducibility (large increases in data quality and Cq consistency), lowered reaction volumes for large cost savings, and nearly a magnitude increase in speed per instrument. We show several comparisons between traditional- and ADE-based pipetting for a qPCR-based workflow.

Clinical microbiology informatics

The clinical microbiology laboratory has responsibilities ranging from characterizing the causative agent in a patient's infection to helping detect global disease outbreaks. All of these processes are increasingly becoming partnered more intimately with informatics. Effective application of informatics tools can increase the accuracy, timeliness, and completeness of microbiology testing while decreasing the laboratory workload, which can lead to optimized laboratory workflow and decreased costs. Informatics is poised to be increasingly relevant in clinical microbiology, with the advent of total laboratory automation, complex instrument interfaces, electronic health records, clinical decision support tools, and the clinical implementation of microbial genome sequencing. This review discusses the diverse informatics aspects that are relevant to the clinical microbiology laboratory, including the following: the microbiology laboratory information system, decision support tools, expert systems, instrument interfaces, total laboratory automation, telemicrobiology, automated image analysis, nucleic acid sequence databases, electronic reporting of infectious agents to public health agencies, and disease outbreak surveillance. The breadth and utility of informatics tools used in clinical microbiology have made them indispensable to contemporary clinical and laboratory practice. Continued advances in technology and development of these informatics tools will further improve patient and public health care in the future.

CanSNPer: a hierarchical genotype classifier of clonal pathogens

SUMMARY

Advances in typing methodologies have recently reformed the field of molecular epidemiology of pathogens. The falling cost of sequencing technologies is creating a deluge of whole genome sequencing data that burdens bioinformatics resources and tool development. In particular, single nucleotide polymorphisms in core genomes of pathogens are recognized as the most important markers for inferring genetic relationships because they are evolutionarily stable and amenable to high-throughput detection methods. Sequence data will provide an excellent opportunity to extend our understanding of infectious disease when the challenge of extracting knowledge from available sequence resources is met. Here, we present an efficient and user-friendly genotype classification pipeline, CanSNPer, based on an easily expandable database of predefined canonical single nucleotide polymorphisms.

AVAILABILITY AND IMPLEMENTATION

All documentation and Python-based source code for the CanSNPer are freely available at http://github.com/adrlar/CanSNPer.

Automated extraction of typing information for bacterial pathogens from whole genome sequence data: Neisseria meningitidis as an exemplar

Whole genome sequence (WGS) data are increasingly used to characterise bacterial pathogens. These data provide detailed information on the genotypes and likely phenotypes of aetiological agents, enabling the relationships of samples from potential disease outbreaks to be established precisely. However, the generation of increasing quantities of sequence data does not, in itself, resolve the problems that many microbiological typing methods have addressed over the last 100 years or so; indeed, providing large volumes of unstructured data can confuse rather than resolve these issues. Here we review the nascent field of storage of WGS data for clinical application and show how curated sequence-based typing schemes on websites have generated an infrastructure that can exploit WGS for bacterial typing efficiently. We review the tools that have been implemented within the PubMLST website to extract clinically useful, strain-characterisation information that can be provided to physicians and public health professionals in a timely, concise and understandable way. These data can be used to inform medical decisions such as how to treat a patient, whether to instigate public health action, and what action might be appropriate. The information is compatible both with previous sequence-based typing data and also with data obtained in the absence of WGS, providing a flexible infrastructure for WGS-based clinical microbiology.