Bacterial population genomics and infectious disease diagnostics

Improved metagenome assemblies through selective enrichment of bacterial genomic DNA from eukaryotic host genomic DNA using ATAC-seq

Genomics can be used to study the complex relationships between hosts and their microbiota. Many bacteria cannot be cultured in the laboratory, making it difficult to obtain adequate amounts of bacterial DNA and to limit host DNA contamination for the construction of metagenome-assembled genomes (MAGs). For example, Wolbachia is a genus of exclusively obligate intracellular bacteria that live in a wide range of arthropods and some nematodes. While Wolbachia endosymbionts are frequently described as facultative reproductive parasites in arthropods, the bacteria are obligate mutualistic endosymbionts of filarial worms. Here, we achieve 50-fold enrichment of bacterial sequences using ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) with Brugia malayi nematodes, containing Wolbachia (wBm). ATAC-seq uses the Tn5 transposase to cut and attach Illumina sequencing adapters to accessible DNA lacking histones, typically thought to be open chromatin. Bacterial and mitochondrial DNA in the lysates are also cut preferentially since they lack histones, leading to the enrichment of these sequences. The benefits of this include minimal tissue input (<1 mg of tissue), a quick protocol (<4 h), low sequencing costs, less bias, correct assembly of lateral gene transfers and no prior sequence knowledge required. We assembled the wBm genome with as few as 1 million Illumina short paired-end reads with >97% coverage of the published genome, compared to only 12% coverage with the standard gDNA libraries. We found significant bacterial sequence enrichment that facilitated genome assembly in previously published ATAC-seq data sets from human cells infected with Mycobacterium tuberculosis and C. elegans contaminated with their food source, the OP50 strain of E. coli. These results demonstrate the feasibility and benefits of using ATAC-seq to easily obtain bacterial genomes to aid in symbiosis, infectious disease, and microbiome research.

Structural genome variants of Pseudomonas aeruginosa clone C and PA14 strains

Plasticity of Pseudomonas aeruginosa chromosomes is mainly driven by an extended accessory genome that is shaped by insertion and deletion events. Further modification of the genome composition can be induced by chromosomal inversion events which lead to relocation of genes in the affected genomic DNA segments, modify the otherwise highly conserved core genome synteny and could even alter the location of the replication terminus. Although the genome of the first sequenced strain, PAO1, displayed such a large genomic inversion, knowledge on such recombination events in the P. aeruginosa population is limited. Several large inversions had been discovered in the late 1990s in cystic fibrosis isolates of the major clonal lineage C by physical genome mapping, and subsequent work on these examples led to the characterization of the DNA at the recombination breakpoints and a presumed recombination mechanism. Since then, the topic was barely addressed in spite of the compilation of thousands of P. aeruginosa genome sequences that are deposited in databases. Due to the use of second-generation sequencing, genome contig assembly had usually followed synteny blueprints provided by the existing reference genome sequences. Inversion detection was not feasible by these approaches, as the respective read lengths did not allow reliable resolution of sequence repeats that are typically found at the borders of inverted segments. In this study, we applied PacBio and MinION long-read sequencing to isolates of the mentioned clone C collection. Confirmation of inversions predicted from the physical mapping data demonstrated that unbiased sequence assembly of such read datasets allows the detection of genomic inversions and the resolution of the recombination breakpoint regions. Additional long-read sequencing of representatives of the other major clonal lineage, PA14, revealed large inversions in several isolates, from cystic fibrosis origin as well as from other sources. These findings indicated that inversion events are not restricted to strains from chronic infection background, but could be widespread in the P. aeruginosa population and contribute to genome plasticity. Moreover, the monitored examples emphasized the role of small mobile DNA units, such as IS elements or transposons, and accessory DNA elements in the inversion-related recombination processes.

Genome Sequence Analysis of Clostridium chauvoei Strains of European Origin and Evaluation of Typing Options for Outbreak Investigations

Black quarter caused by Clostridium (C.) chauvoei is an important bacterial disease that affects cattle and sheep with high mortality. A comparative genomics analysis of 64 C. chauvoei strains, most of European origin and a few of non-European and unknown origin, was performed. The pangenome analysis showed limited new gene acquisition for the species. The accessory genome involved prophages and genomic islands, with variations in gene composition observed in a few strains. This limited accessory genome may indicate that the species replicates only in the host or that an active CRISPR/Cas system provides immunity to foreign genetic elements. All strains contained a CRISPR type I-B system and it was confirmed that the unique spacer sequences therein can be used to differentiate strains. Homologous recombination events, which may have contributed to the evolution of this pathogen, were less frequent compared to other related species from the genus. Pangenome single nucleotide polymorphism (SNP) based phylogeny and clustering indicate diverse clusters related to geographical origin. Interestingly the identified SNPs were mostly non-synonymous. The study demonstrates the possibility of the existence of polymorphic populations in one host, based on strain variability observed for strains from the same animal and strains from different animals of one outbreak. The study also demonstrates that new outbreak strains are mostly related to earlier outbreak strains from the same farm/region. This indicates the last common ancestor strain from one farm can be crucial to understand the genetic changes and epidemiology occurring at farm level. Known virulence factors for the species were highly conserved among the strains. Genetic elements involved in Nicotinamide adenine dinucleotide (NAD) precursor synthesis (via nadA, nadB, and nadC metabolic pathway) which are known as potential anti-virulence loci are completely absent in C. chauvoei compared to the partial inactivation in C. septicum. A novel core-genome MLST based typing method was compared to sequence typing based on CRISPR spacers to evaluate the usefulness of the methods for outbreak investigations.

Heuristic and Hierarchical-Based Population Mining of Salmonella enterica Lineage I Pan-Genomes as a Platform to Enhance Food Safety

The recent incorporation of bacterial whole-genome sequencing (WGS) into Public Health laboratories has enhanced foodborne outbreak detection and source attribution. As a result, large volumes of publicly available datasets can be used to study the biology of foodborne pathogen populations at an unprecedented scale. To demonstrate the application of a heuristic and agnostic hierarchical population structure guided pan-genome enrichment analysis (PANGEA), we used populations of S. enterica lineage I to achieve two main objectives: (i) show how hierarchical population inquiry at different scales of resolution can enhance ecological and epidemiological inquiries; and (ii) identify population-specific inferable traits that could provide selective advantages in food production environments. Publicly available WGS data were obtained from NCBI database for three serovars of Salmonella enterica subsp. enterica lineage I (S. Typhimurium, S. Newport, and S. Infantis). Using the hierarchical genotypic classifications (Serovar, BAPS1, ST, cgMLST), datasets from each of the three serovars showed varying degrees of clonal structuring. When the accessory genome (PANGEA) was mapped onto these hierarchical structures, accessory loci could be linked with specific genotypes. A large heavy-metal resistance mobile element was found in the Monophasic ST34 lineage of S. Typhimurium, and laboratory testing showed that Monophasic isolates have on average a higher degree of copper resistance than the Biphasic ones. In S. Newport, an extra sugE gene copy was found among most isolates of the ST45 lineage, and laboratory testing of multiple isolates confirmed that isolates of S. Newport ST45 were on average less sensitive to the disinfectant cetylpyridimium chloride than non-ST45 isolates. Lastly, data-mining of the accessory genomic content of S. Infantis revealed two cryptic Ecotypes with distinct accessory genomic content and distinct ecological patterns. Poultry appears to be the major reservoir for Ecotype 1, and temporal analysis further suggested a recent ecological succession, with Ecotype 2 apparently being displaced by Ecotype 1. Altogether, the use of a heuristic hierarchical-based population structure analysis that includes bacterial pan-genomes (core and accessory genomes) can (1) improve genomic resolution for mapping populations and accessing epidemiological patterns; and (2) define lineage-specific informative loci that may be associated with survival in the food chain.

ProkEvo: an automated, reproducible, and scalable framework for high-throughput bacterial population genomics analyses

Whole Genome Sequence (WGS) data from bacterial species is used for a variety of applications ranging from basic microbiological research, diagnostics, and epidemiological surveillance. The availability of WGS data from hundreds of thousands of individual isolates of individual microbial species poses a tremendous opportunity for discovery and hypothesis-generating research into ecology and evolution of these microorganisms. Flexibility, scalability, and user-friendliness of existing pipelines for population-scale inquiry, however, limit applications of systematic, population-scale approaches. Here, we present ProkEvo, an automated, scalable, reproducible, and open-source framework for bacterial population genomics analyses using WGS data. ProkEvo was specifically developed to achieve the following goals: (1) Automation and scaling of complex combinations of computational analyses for many thousands of bacterial genomes from inputs of raw Illumina paired-end sequence reads; (2) Use of workflow management systems (WMS) such as Pegasus WMS to ensure reproducibility, scalability, modularity, fault-tolerance, and robust file management throughout the process; (3) Use of high-performance and high-throughput computational platforms; (4) Generation of hierarchical-based population structure analysis based on combinations of multi-locus and Bayesian statistical approaches for classification for ecological and epidemiological inquiries; (5) Association of antimicrobial resistance (AMR) genes, putative virulence factors, and plasmids from curated databases with the hierarchically-related genotypic classifications; and (6) Production of pan-genome annotations and data compilation that can be utilized for downstream analysis such as identification of population-specific genomic signatures. The scalability of ProkEvo was measured with two datasets comprising significantly different numbers of input genomes (one with ~2,400 genomes, and the second with ~23,000 genomes). Depending on the dataset and the computational platform used, the running time of ProkEvo varied from ~3-26 days. ProkEvo can be used with virtually any bacterial species, and the Pegasus WMS uniquely facilitates addition or removal of programs from the workflow or modification of options within them. To demonstrate versatility of the ProkEvo platform, we performed a hierarchical-based population structure analyses from available genomes of three distinct pathogenic bacterial species as individual case studies. The specific case studies illustrate how hierarchical analyses of population structures, genotype frequencies, and distribution of specific gene functions can be integrated into an analysis. Collectively, our study shows that ProkEvo presents a practical viable option for scalable, automated analyses of bacterial populations with direct applications for basic microbiology research, clinical microbiological diagnostics, and epidemiological surveillance.

Multilocus Sequence Typing of Staphylococcus aureus

Multilocus sequence typing (MLST) has been successfully used to differentiate and trace the bacterial species and pathogens that cause outbreaks or epidemics of infectious diseases. MLST provides a powerful solution for molecular epidemiological characterization of bacterial strains, including Staphylococcus aureus, by using the sequences of the internal region of seven housekeeping genes. In previous studies, we utilized MLST to analyze the genotypes of S. aureus isolates from pediatric patients with cystic fibrosis and revealed three prevalent ST types, including ST5, ST30, and ST8 in these isolates. In this chapter, we describe a detailed procedure of MLST for genotyping S. aureus.

Interest of bacterial pangenome analyses in clinical microbiology

Thanks to the progress and decreasing costs in genome sequencing technologies, more than 250,000 bacterial genomes are currently available in public databases, covering most, if not all, of the major human-associated phylogenetic groups of these microorganisms, pathogenic or not. In addition, for many of them, sequences from several strains of a given species are available, thus enabling to evaluate their genetic diversity and study their evolution. In addition, the significant cost reduction of bacterial whole genome sequencing as well as the rapid increase in the number of available bacterial genomes have prompted the development of pangenomic software tools. The study of bacterial pangenome has many applications in clinical microbiology. It can unveil the pathogenic potential and ability of bacteria to resist antimicrobials as well identify specific sequences and predict antigenic epitopes that allow molecular or serologic assays and vaccines to be designed. Bacterial pangenome constitutes a powerful method for understanding the history of human bacteria and relating these findings to diagnosis in clinical microbiology laboratories in order to optimize patient management.

Surveillance-embedded genomic outbreak resolution of methicillin-susceptible Staphylococcus aureus in a neonatal intensive care unit

We observed an increase in methicillin-susceptible Staphylococcus aureus (MSSA) infections at a Dutch neonatal intensive care unit. Weekly neonatal MSSA carriage surveillance and cross-sectional screenings of health care workers (HCWs) were available for outbreak tracing. Traditional clustering of MSSA isolates by spa typing and Multiple-Locus Variable number tandem repeat Analysis (MLVA) suggested that nosocomial transmission had contributed to the infections. We investigated whether whole-genome sequencing (WGS) of MSSA surveillance would provide additional evidence for transmission. MSSA isolates from neonatal infections, carriage surveillance, and HCWs were subjected to WGS and bioinformatic analysis for identification and localization of high-quality single nucleotide polymorphisms, and in-depth analysis of subsets of isolates. By measuring the genetic diversity in background surveillance, we defined transmission-level relatedness and identified isolates that had been unjustly assigned to clusters based on MLVA, while spa typing was concordant but of insufficient resolution. Detailing particular subsets of isolates provided evidence that HCWs were involved in multiple outbreaks, yet it alleviated concerns about one particular HCW. The improved resolution and accuracy of genomic outbreak analyses substantially altered the view on outbreaks, along with apposite measures. Therefore, inclusion of the circulating background population has the potential to overcome current issues in genomic outbreak inference.

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

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

Molecular epidemiology of Staphylococcus aureus bacteremia in a single large Minnesota medical center in 2015 as assessed using MLST, core genome MLST and spa typing

Staphylococcus aureus is a leading cause of bacteremia in hospitalized patients. Whether or not S. aureus bacteremia (SAB) is associated with clonality, implicating potential nosocomial transmission, has not, however, been investigated. Herein, we examined the epidemiology of SAB using whole genome sequencing (WGS). 152 SAB isolates collected over the course of 2015 at a single large Minnesota medical center were studied. Staphylococcus protein A (spa) typing was performed by PCR/Sanger sequencing; multilocus sequence typing (MLST) and core genome MLST (cgMLST) were determined by WGS. Forty-eight isolates (32%) were methicillin–resistant S. aureus (MRSA). The isolates encompassed 66 spa types, clustered into 11 spa clonal complexes (CCs) and 10 singleton types. 88% of 48 MRSA isolates belonged to spa CC-002 or -008. Methicillin-susceptible S. aureus (MSSA) isolates were more genotypically diverse, with 61% distributed across four spa CCs (CC-002, CC-012, CC-008 and CC-084). By MLST, there was 31 sequence types (STs), including 18 divided into 6 CCs and 13 singleton STs. Amongst MSSA isolates, the common MLST clones were CC5 (23%), CC30 (19%), CC8 (15%) and CC15 (11%). Common MRSA clones were CC5 (67%) and CC8 (25%); there were no MRSA isolates in CC45 or CC30. By cgMLST analysis, there were 9 allelic differences between two isolates, with the remaining 150 isolates differing from each other by over 40 alleles. The two isolates were retroactively epidemiologically linked by medical record review. Overall, cgMLST analysis resulted in higher resolution epidemiological typing than did multilocus sequence or spa typing.

The single-species metagenome: subtyping Staphylococcus aureus core genome sequences from shotgun metagenomic data

In this study we developed a genome-based method for detecting Staphylococcus aureus subtypes from metagenome shotgun sequence data. We used a binomial mixture model and the coverage counts at >100,000 known S. aureus SNP (single nucleotide polymorphism) sites derived from prior comparative genomic analysis to estimate the proportion of 40 subtypes in metagenome samples. We were able to obtain >87% sensitivity and >94% specificity at 0.025X coverage for S. aureus. We found that 321 and 149 metagenome samples from the Human Microbiome Project and metaSUB analysis of the New York City subway, respectively, contained S. aureus at genome coverage >0.025. In both projects, CC8 and CC30 were the most common S. aureus clonal complexes encountered. We found evidence that the subtype composition at different body sites of the same individual were more similar than random sampling and more limited evidence that certain body sites were enriched for particular subtypes. One surprising finding was the apparent high frequency of CC398, a lineage often associated with livestock, in samples from the tongue dorsum. Epidemiologic analysis of the HMP subject population suggested that high BMI (body mass index) and health insurance are possibly associated with S. aureus carriage but there was limited power to identify factors linked to carriage of even the most common subtype. In the NYC subway data, we found a small signal of geographic distance affecting subtype clustering but other unknown factors influence taxonomic distribution of the species around the city.

Whole genome protein microarrays for serum profiling of immunodominant antigens of Bacillus anthracis

A commercial Bacillus anthracis (Anthrax) whole genome protein microarray has been used to identify immunogenic Anthrax proteins (IAP) using sera from groups of donors with (a) confirmed B. anthracis naturally acquired cutaneous infection, (b) confirmed B. anthracis intravenous drug use-acquired infection, (c) occupational exposure in a wool-sorters factory, (d) humans and rabbits vaccinated with the UK Anthrax protein vaccine and compared to naïve unexposed controls. Anti-IAP responses were observed for both IgG and IgA in the challenged groups; however the anti-IAP IgG response was more evident in the vaccinated group and the anti-IAP IgA response more evident in the B. anthracis-infected groups. Infected individuals appeared somewhat suppressed for their general IgG response, compared with other challenged groups. Immunogenic protein antigens were identified in all groups, some of which were shared between groups whilst others were specific for individual groups. The toxin proteins were immunodominant in all vaccinated, infected or other challenged groups. However, a number of other chromosomally-located and plasmid encoded open reading frame proteins were also recognized by infected or exposed groups in comparison to controls. Some of these antigens e.g., BA4182 are not recognized by vaccinated individuals, suggesting that there are proteins more specifically expressed by live Anthrax spores in vivo that are not currently found in the UK licensed Anthrax Vaccine (AVP). These may perhaps be preferentially expressed during infection and represent expression of alternative pathways in the B. anthracis “infectome.” These may make highly attractive candidates for diagnostic and vaccine biomarker development as they may be more specifically associated with the infectious phase of the pathogen. A number of B. anthracis small hypothetical protein targets have been synthesized, tested in mouse immunogenicity studies and validated in parallel using human sera from the same study.

SRST2: Rapid genomic surveillance for public health and hospital microbiology labs

Rapid molecular typing of bacterial pathogens is critical for public health epidemiology, surveillance and infection control, yet routine use of whole genome sequencing (WGS) for these purposes poses significant challenges. Here we present SRST2, a read mapping-based tool for fast and accurate detection of genes, alleles and multi-locus sequence types (MLST) from WGS data. Using >900 genomes from common pathogens, we show SRST2 is highly accurate and outperforms assembly-based methods in terms of both gene detection and allele assignment. We include validation of SRST2 within a public health laboratory, and demonstrate its use for microbial genome surveillance in the hospital setting. In the face of rising threats of antimicrobial resistance and emerging virulence among bacterial pathogens, SRST2 represents a powerful tool for rapidly extracting clinically useful information from raw WGS data. Source code is available from http://katholt.github.io/srst2/.

SRST2: Rapid genomic surveillance for public health and hospital microbiology labs

Rapid molecular typing of bacterial pathogens is critical for public health epidemiology, surveillance and infection control, yet routine use of whole genome sequencing (WGS) for these purposes poses significant challenges. Here we present SRST2, a read mapping-based tool for fast and accurate detection of genes, alleles and multi-locus sequence types (MLST) from WGS data. Using >900 genomes from common pathogens, we show SRST2 is highly accurate and outperforms assembly-based methods in terms of both gene detection and allele assignment. We include validation of SRST2 within a public health laboratory, and demonstrate its use for microbial genome surveillance in the hospital setting. In the face of rising threats of antimicrobial resistance and emerging virulence among bacterial pathogens, SRST2 represents a powerful tool for rapidly extracting clinically useful information from raw WGS data. Source code is available from http://katholt.github.io/srst2/.

Standing of nucleic acid testing strategies in veterinary diagnosis laboratories to uncover Mycobacterium tuberculosis complex members

Nucleic acid testing (NAT) designate any molecular approach used for the detection, identification, and characterization of pathogenic microorganisms, enabling the rapid, specific, and sensitive diagnostic of infectious diseases, such as tuberculosis. These assays have been widely used since the 90s of the last century in human clinical laboratories and, subsequently, also in veterinary diagnostics. Most NAT strategies are based in the polymerase chain reaction (PCR) and its several enhancements and variations. From the conventional PCR, real-time PCR and its combinations, isothermal DNA amplification, to the nanotechnologies, here we review how the NAT assays have been applied to decipher if and which member of the Mycobacterium tuberculosis complex is present in a clinical sample. Recent advances in DNA sequencing also brought new challenges and have made possible to generate rapidly and at a low cost, large amounts of sequence data. This revolution with the high-throughput sequencing (HTS) technologies makes whole genome sequencing (WGS) and metagenomics the trendiest NAT strategies, today. The ranking of NAT techniques in the field of clinical diagnostics is rising, and we provide a SWOT (Strengths, Weaknesses, Opportunities, and Threats) analysis with our view of the use of molecular diagnostics for detecting tuberculosis in veterinary laboratories, notwithstanding the gold standard being still the classical culture of the agent. The complementary use of both classical and molecular diagnostics approaches is recommended to speed the diagnostic, enabling a fast decision by competent authorities and rapid tackling of the disease.

Nucleic acid detection technologies and marker molecules in bacterial diagnostics

There is a growing need for quick and reliable methods for microorganism detection and identification worldwide. Although traditional culture-based technologies are trustworthy and accurate at a relatively low cost, they are also time- and labor-consuming and are limited to culturable bacteria. Those weaknesses have created a necessity for alternative technologies that are capable for faster and more precise bacterial identification from medical, food or environmental samples. The most common current approach is to analyze the nucleic acid component of analyte solution and determine the bacterial composition according to the specific nucleic acid profiles that are present. This review aims to give an up-to-date overview of different nucleic acid target sequences and respective analytical technologies.

Four genotyping schemes for phylogenetic analysis of Pseudomonas aeruginosa: comparison of their congruence with multi-locus sequence typing

Several molecular typing schemes have been proposed to differentiate among isolates and clonal groups, and hence establish epidemiological or phylogenetic links. It has been widely accepted that multi-locus sequence typing (MLST) is the gold standard for phylogenetic typing/long-term epidemiological surveillance, but other recently described methods may be easier to carry out, especially in settings with limited access to DNA sequencing. Comparing the performance of such techniques to MLST is therefore of relevance. A study was therefore carried out with a collection of P. aeruginosa strains (n = 133) typed by four typing schemes: MLST, multiple-locus variable number tandem repeat analysis (MLVA), pulsed-field gel electrophoresis (PFGE) and the commercial DiversiLab microbial typing system (DL). The aim of this study was to compare the results of each typing method with MLST. The Simpson's indices of diversity were 0.989, 0.980, 0.961 and 0.906 respectively for PFGE, MLVA, DL and MLST. The congruence between techniques was measured by the adjusted Wallace index (W): this coefficient indicates the probability that a pair of isolates which is assigned to the same type by one typing method is also typed as identical by the other. In this context, the congruence between techniques was recorded as follow: MLVA-type to predict MLST-type (93%), PFGE to MLST (92%), DL to MLST (64.2%), PFGE to MLVA (63.5%) and PFGE to DL (61.7%). Conversely, for all above combinations, prediction was very poor. The congruence was increased at the clonal complex (CC) level. MLST is regarded the gold standard for phylogenetic classification of bacteria, but is rather laborious to carry out in many settings. Our data suggest that MLVA can predict the MLST-type with high accuracy, and even higher when studying the clonal complex level. Of the studied three techniques MLVA was therefore the best surrogate method to predict MLST.

A TaqMan-based multiplex real-time PCR assay for the simultaneous detection of Wuchereria bancrofti and Brugia malayi

With the Global Program for the Elimination of Lymphatic Filariasis continuing to make strides towards disease eradication, many locations endemic for the causative parasites of lymphatic filariasis are realizing a substantial decrease in levels of infection and rates of disease transmission. However, with measures of disease continuing to decline, the need for time-saving and economical molecular diagnostic assays capable of detecting low levels of parasite presence is increasing. This need is greatest in locations co-endemic for both Wuchereria bancrofti and Brugia parasites because testing for both causative agents individually results in significant increases in labor and reagent costs. Here we describe a multiplex, TaqMan-based, real-time PCR assay capable of simultaneously detecting W. bancrofti and Brugia malayi DNA extracted from human bloodspots or vector mosquito pools. With comparable sensitivity to established singleplex assays, this assay provides significant cost and labor savings for disease monitoring efforts in co-endemic locations.

Software for selecting the most informative sets of genomic loci for multi-target microbial typing

Background

High-throughput sequencing can identify numerous potential genomic targets for microbial strain typing, but identification of the most informative combinations requires the use of computational screening tools. This paper describes novel software – Automated Selection of Typing Target Subsets (AuSeTTS) - that allows intelligent selection of optimal targets for pathogen strain typing. The objective of this software is to maximise both discriminatory power, using Simpson’s index of diversity (D), and concordance with existing typing methods, using the adjusted Wallace coefficient (AW). The program interrogates molecular typing results for panels of isolates, based on large target sets, and iteratively examines each target, one-by-one, to determine the most informative subset.

Results

AuSeTTS was evaluated using three target sets: 51 binary targets (13 toxin genes, 16 phage-related loci and 22 SCCmec elements), used for multilocus typing of 153 methicillin-resistant Staphylococcus aureus (MRSA) isolates; 17 MLVA loci in 502 Streptococcus pneumoniae isolates from the MLVA database (http://www.mlva.eu) and 12 MLST loci for 98 Cryptococcus spp. isolates.

The maximum D for MRSA, 0.984, was achieved with a subset of 20 targets and a D value of 0.954 with 7 targets. Twelve targets predicted MLST with a maximum AW of 0.9994. All 17 S. pneumoniae MLVA targets were required to achieve maximum D of 0.997, but 4 targets reached D of 0.990. Twelve targets predicted pneumococcal serotype with a maximum AW of 0.899 and 9 predicted MLST with maximum AW of 0.963. Eight of the 12 MLST loci were sufficient to achieve the maximum D of 0.963 for Cryptococcus spp.

Conclusions

Computerised analysis with AuSeTTS allows rapid selection of the most discriminatory targets for incorporation into typing schemes. Output of the program is presented in both tabular and graphical formats and the software is available for free download from http://www.cidmpublichealth.org/pages/ausetts.html.

Bacterial genomics in infectious disease and the clinical pathology laboratory

CONTEXT

Throughout history, technologic advancements have fueled the engine of innovation, which, in turn, has driven discovery. Accordingly, recent advancements in DNA sequencing technology are revolutionizing bacterial genomics.

OBJECTIVE

To review important developments from the literature. The current state of bacterial genomics, with an emphasis on human pathogens and the clinical pathology laboratory, will be discussed.

DATA SOURCES

A comprehensive review was performed of the relevant literature indexed in PubMed (National Library of Medicine) and referenced medical texts.

CONCLUSIONS

Many important discoveries bearing on infectious disease research and pathology laboratory practice have been achieved through whole-genome sequencing strategies. Bacterial genomics has improved our understanding of molecular pathogenesis, host-pathogen interactions, and antibiotic-resistance mechanisms. Bacterial genomics has also facilitated the study of population structures, epidemics and outbreaks, and newly identified pathogens. Many opportunities now exist for clinical pathologists to contribute to bacterial genomics, including in the design of new diagnostic tests, therapeutic agents, and vaccines.

Population structure of hyperinvasive serotype 12F, clonal complex 218 Streptococcus pneumoniae revealed by multilocus boxB sequence typing

At least four outbreaks of invasive disease caused by serotype 12F, clonal complex 218 Streptococcus pneumoniae have occurred in the United States over the past two decades. We studied the population structure of this clonal complex using a sample of 203 outbreak and surveillance isolates that were collected over 22 years from 34 US states and eight other countries. Conventional multilocus sequence typing identified five types and distinguished a single outbreak from the others. To improve typing resolution, multilocus boxB sequence typing (MLBT) was developed from 10 variable boxB minisatellite loci. MLBT identified 86 types and distinguished between each of the four outbreaks. Diversity across boxB loci tended to be positively correlated with repeat array size and, overall, best fit the infinite alleles mutation model. Multilocus linkage disequilibrium was strong, but pairwise disequilibrium decreased with the physical distance between loci and was strongest in one large region of the chromosome, indicating recent recombinations. Two major clusters were identified in the sample, and they were differentiated geographically, as western and more easterly US clusters, and temporally, as clusters that predominated before and after the licensure of pneumococcal conjugate vaccines. The diversity and linkage disequilibrium within these two clusters also differed, suggesting different population dynamics. MLBT revealed hidden aspects of the population structure of these hyperinvasive pneumococci, and it may provide a useful adjunct tool for outbreak investigations, surveillance, and population genetics studies of other pneumococcal clonal complexes.

Advances in bacterial transcriptome and transposon insertion-site profiling using second-generation sequencing

The arrival of second-generation sequencing has revolutionized the study of bacteria within a short period. The sequence information generated from these platforms has helped in our understanding of bacterial development, adaptation and diversity and how bacteria cause disease. Furthermore, these technologies have quickly been adapted for high-throughput studies that were previously performed using DNA cloning or microarray-based applications. This has facilitated a more comprehensive study of bacterial transcriptomes through RNA sequencing (RNA-Seq) and the systematic determination of gene function by 'transposon monitoring'. In this review, we provide an outline of these powerful tools and the in silico analyses used in their application, and also highlight the biological questions being addressed in these approaches.

Interplay of recombination and selection in the genomes of Chlamydia trachomatis

Background

Chlamydia trachomatis is an obligate intracellular bacterial parasite, which causes several severe and debilitating diseases in humans. This study uses comparative genomic analyses of 12 complete published C. trachomatis genomes to assess the contribution of recombination and selection in this pathogen and to understand the major evolutionary forces acting on the genome of this bacterium.

Results

The conserved core genes of C. trachomatis are a large proportion of the pan-genome: we identified 836 core genes in C. trachomatis out of a range of 874-927 total genes in each genome. The ratio of recombination events compared to mutation (ρ/θ) was 0.07 based on ancestral reconstructions using the ClonalFrame tool, but recombination had a significant effect on genetic diversification (r/m = 0.71). The distance-dependent decay of linkage disequilibrium also indicated that C. trachomatis populations behaved intermediately between sexual and clonal extremes. Fifty-five genes were identified as having a history of recombination and 92 were under positive selection based on statistical tests. Twenty-three genes showed evidence of being under both positive selection and recombination, which included genes with a known role in virulence and pathogencity (e.g., ompA, pmps, tarp). Analysis of inter-clade recombination flux indicated non-uniform currents of recombination between clades, which suggests the possibility of spatial population structure in C. trachomatis infections.

Conclusions

C. trachomatis is the archetype of a bacterial species where recombination is relatively frequent yet gene gains by horizontal gene transfer (HGT) and losses (by deletion) are rare. Gene conversion occurs at sites across the whole C. trachomatis genome but may be more often fixed in genes that are under diversifying selection. Furthermore, genome sequencing will reveal patterns of serotype specific gene exchange and selection that will generate important research questions for understanding C. trachomatis pathogenesis.

Reviewers

This article was reviewed by Dr. Jeremy Selengut, Dr. Lee S. Katz (nominated by Dr. I. King Jordan) and Dr. Arcady Mushegian.