Roles of bacteriophages, plasmids and CRISPR immunity in microbial community dynamics revealed using time-series integrated meta-omics.

Viruses and plasmids (invasive mobile genetic elements (iMGEs)) have important roles in shaping microbial communities, but their dynamic interactions with CRISPR-based immunity remain unresolved. We analysed generation-resolved iMGE-host dynamics spanning one and a half years in a microbial consortium from a biological wastewater treatment plant using integrated meta-omics. We identified 31 bacterial metagenome-assembled genomes encoding complete CRISPR-Cas systems and their corresponding iMGEs. CRISPR-targeted plasmids outnumbered their bacteriophage counterparts by at least fivefold, highlighting the importance of CRISPR-mediated defence against plasmids. Linear modelling of our time-series data revealed that the variation in plasmid abundance over time explained more of the observed community dynamics than phages. Community-scale CRISPR-based plasmid-host and phage-host interaction networks revealed an increase in CRISPR-mediated interactions coinciding with a decrease in the dominant 'Candidatus Microthrix parvicella' population. Protospacers were enriched in sequences targeting genes involved in the transmission of iMGEs. Understanding the factors shaping the fitness of specific populations is necessary to devise control strategies for undesirable species and to predict or explain community-wide phenotypes.

Integration of time-series meta-omics data reveals how microbial ecosystems respond to disturbance.

The development of reliable, mixed-culture biotechnological processes hinges on understanding how microbial ecosystems respond to disturbances. Here we reveal extensive phenotypic plasticity and niche complementarity in oleaginous microbial populations from a biological wastewater treatment plant. We perform meta-omics analyses (metagenomics, metatranscriptomics, metaproteomics and metabolomics) on in situ samples over 14 months at weekly intervals. Based on 1,364 de novo metagenome-assembled genomes, we uncover four distinct fundamental niche types. Throughout the time-series, we observe a major, transient shift in community structure, coinciding with substrate availability changes. Functional omics data reveals extensive variation in gene expression and substrate usage amongst community members. Ex situ bioreactor experiments confirm that responses occur within five hours of a pulse disturbance, demonstrating rapid adaptation by specific populations. Our results show that community resistance and resilience are a function of phenotypic plasticity and niche complementarity, and set the foundation for future ecological engineering efforts.

First draft genome sequence of a strain belonging to the Zoogloea genus and its gene expression in situ.

The Gram-negative beta-proteobacterium Zoogloea sp. LCSB751 (LMG 29444) was newly isolated from foaming activated sludge of a municipal wastewater treatment plant. Here, we describe its draft genome sequence and annotation together with a general physiological and genomic analysis, as the first sequenced representative of the Zoogloea genus. Moreover, Zoogloea sp. gene expression in its environment is described using metatranscriptomic data obtained from the same treatment plant. The presented genomic and transcriptomic information demonstrate a pronounced capacity of this genus to synthesize poly-ß-hydroxyalkanoate within wastewater.

Phylogeographic, genomic, and meropenem susceptibility analysis of Burkholderia ubonensis.

The bacterium Burkholderia ubonensis is commonly co-isolated from environmental specimens harbouring the melioidosis pathogen, Burkholderia pseudomallei. B. ubonensis has been reported in northern Australia and Thailand but not North America, suggesting similar geographic distribution to B. pseudomallei. Unlike most other Burkholderia cepacia complex (Bcc) species, B. ubonensis is considered non-pathogenic, although its virulence potential has not been tested. Antibiotic resistance in B. ubonensis, particularly towards drugs used to treat the most severe B. pseudomallei infections, has also been poorly characterised. This study examined the population biology of B. ubonensis, and includes the first reported isolates from the Caribbean. Phylogenomic analysis of 264 B. ubonensis genomes identified distinct clades that corresponded with geographic origin, similar to B. pseudomallei. A small proportion (4%) of strains lacked the 920kb chromosome III replicon, with discordance of presence/absence amongst genetically highly related strains, demonstrating that the third chromosome of B. ubonensis, like other Bcc species, probably encodes for a nonessential pC3 megaplasmid. Multilocus sequence typing using the B. pseudomallei scheme revealed that one-third of strains lack the "housekeeping" narK locus. In comparison, all strains could be genotyped using the Bcc scheme. Several strains possessed high-level meropenem resistance (≥32 µg/mL), a concern due to potential transmission of this phenotype to B. pseudomallei. In silico analysis uncovered a high degree of heterogeneity among the lipopolysaccharide O-antigen cluster loci, with at least 35 different variants identified. Finally, we show that Asian B. ubonensis isolate RF23-BP41 is avirulent in the BALB/c mouse model via a subcutaneous route of infection. Our results provide several new insights into the biology of this understudied species.

Mechanisms of Resistance to Folate Pathway Inhibitors in Burkholderia pseudomallei: Deviation from the Norm.

The trimethoprim and sulfamethoxazole combination, co-trimoxazole, plays a vital role in the treatment of Burkholderia pseudomallei infections. Previous studies demonstrated that the B. pseudomallei BpeEF-OprC efflux pump confers widespread trimethoprim resistance in clinical and environmental isolates, but this is not accompanied by significant resistance to co-trimoxazole. Using the excluded select-agent strain B. pseudomallei Bp82, we now show that in vitro acquired trimethoprim versus co-trimoxazole resistance is mainly mediated by constitutive BpeEF-OprC expression due to bpeT mutations or by BpeEF-OprC overexpression due to bpeS mutations. Mutations in bpeT affect the carboxy-terminal effector-binding domain of the BpeT LysR-type activator protein. Trimethoprim resistance can also be mediated by dihydrofolate reductase (FolA) target mutations, but this occurs rarely unless BpeEF-OprC is absent. BpeS is a transcriptional regulator that is 62% identical to BpeT. Mutations affecting the BpeS DNA-binding or carboxy-terminal effector-binding domains result in constitutive BpeEF-OprC overexpression, leading to trimethoprim and sulfamethoxazole efflux and thus to co-trimoxazole resistance. The majority of laboratory-selected co-trimoxazole-resistant mutants often also contain mutations in folM, encoding a pterin reductase. Genetic analyses of these mutants established that both bpeS mutations and folM mutations contribute to co-trimoxazole resistance, although the exact role of folM remains to be determined. Mutations affecting bpeT, bpeS, and folM are common in co-trimoxazole-resistant clinical isolates, indicating that mutations affecting these genes are clinically significant. Co-trimoxazole resistance in B. pseudomallei is a complex phenomenon, which may explain why resistance to this drug is rare in this bacterium.IMPORTANCEBurkholderia pseudomallei causes melioidosis, a tropical disease that is difficult to treat. The bacterium's resistance to antibiotics limits therapeutic options. The paucity of orally available drugs further complicates therapy. The oral drug of choice is co-trimoxazole, a combination of trimethoprim and sulfamethoxazole. These antibiotics target two distinct enzymes, FolA (dihydrofolate reductase) and FolP (dihydropteroate synthase), in the bacterial tetrahydrofolate biosynthetic pathway. Although co-trimoxazole resistance is minimized due to two-target inhibition, bacterial resistance due to folA and folP mutations does occur. Co-trimoxazole resistance in B. pseudomallei is rare and has not yet been studied. Co-trimoxazole resistance in this bacterium employs a novel strategy involving differential regulation of BpeEF-OprC efflux pump expression that determines the drug resistance profile. Contributing are mutations affecting folA, but not folP, and folM, a folate pathway-associated gene whose function is not yet well understood and which has not been previously implicated in folate inhibitor resistance in clinical isolates.

Mapping and comparing bacterial microbiota in the sinonasal cavity of healthy, allergic rhinitis, and chronic rhinosinusitis subjects.

BACKGROUND: The role of microbiota in sinonasal inflammation can be further understood by targeted sampling of healthy and diseased subjects. We compared the microbiota of the middle meatus (MM) and inferior meatus (IM) in healthy, allergic rhinitis (AR), and chronic rhinosinusitis (CRS) subjects to characterize intrasubject, intersubject, and intergroup differences. METHODS: Subjects were recruited in the office, and characterized into healthy, AR, and CRS groups. Endoscopically-guided swab samples were obtained from the MM and IM bilaterally. Bacterial microbiota were characterized by sequencing the V3-V4 region of the 16S ribosomal RNA (rRNA) gene. RESULTS: Intersubject microbiome analyses were conducted in 65 subjects: 8 healthy, 11 AR, and 46 CRS (25 CRS with nasal polyps [CRSwNP]; 21 CRS without nasal polyps [CRSsNP]). Intrasubject analyses were conducted for 48 individuals (4 controls, 11 AR, 8 CRSwNP, and 15 CRSwNP). There was considerable intersubject microbiota variability, but intrasubject profiles were similar (p = 0.001, nonparametric t test). Intrasubject bacterial diversity was significantly reduced in MM of CRSsNP subjects compared to IM samples (p = 0.022, nonparametric t test). CRSsNP MM samples were enriched in Streptococcus, Haemophilus, and Fusobacterium spp. but exhibited loss of diversity compared to healthy, CRSwNP, and AR subject-samples (p < 0.05; nonparametric t test). CRSwNP patients were enriched in Staphylococcus, Alloiococcus, and Corynebacterium spp. CONCLUSION: This study presents the sinonasal microbiome profile in one of the larger populations of non-CRS and CRS subjects, and is the first office-based cohort in the literature. In contrast to healthy, AR, and CRSwNP subjects, CRSsNP MM samples exhibited decreased microbiome diversity and anaerobic enrichment. CRSsNP MM samples had reduced diversity compared to same-subject IM samples, a novel finding.

More than 50% of Clostridium difficile Isolates from Pet Dogs in Flagstaff, USA, Carry Toxigenic Genotypes.

Nosocomial acquisition of Clostridium difficile is well documented, yet recent studies have highlighted the importance of community acquired infections and identified community associated reservoirs for this pathogen. Multiple studies have implicated companion pets and farm animals as possible sources of community acquired C. difficile infections in humans. To explore the potential role of pet dogs in human C. difficile infections we systematically collected canine fecal samples (n = 197) in Flagstaff, AZ. Additionally, nineteen fecal samples were collected at a local veterinary clinic from diarrheic dogs. We used these combined samples to investigate important questions regarding C. difficile colonization in pet canines: 1) What is the prevalence and diversity of C. difficile in this companion pet population, and 2) Do C. difficile isolates collected from canines genetically overlap with isolates that cause disease in humans? We used a two-step sequence typing approach, including multilocus sequence typing to determine the overall genetic diversity of C. difficile present in Flagstaff canines, and whole-genome sequencing to assess the fine-scale diversity patterns within identical multilocus sequence types from isolates obtained within and among multiple canine hosts. We detected C. difficile in 17% of the canine fecal samples with 10% containing toxigenic strains that are known to cause human disease. Sequencing analyses revealed similar genotypes in dogs and humans. These findings suggest that companion pets are a potential source of community acquired C. difficile infections in humans.

A Bacillus anthracis Genome Sequence from the Sverdlovsk 1979 Autopsy Specimens.

Anthrax is a zoonotic disease that occurs naturally in wild and domestic animals but has been used by both state-sponsored programs and terrorists as a biological weapon. A Soviet industrial production facility in Sverdlovsk, USSR, proved deficient in 1979 when a plume of spores was accidentally released and resulted in one of the largest known human anthrax outbreaks. In order to understand this outbreak and others, we generated a Bacillus anthracis population genetic database based upon whole-genome analysis to identify all single-nucleotide polymorphisms (SNPs) across a reference genome. Phylogenetic analysis has defined three major clades (A, B, and C), B and C being relatively rare compared to A. The A clade has numerous subclades, including a major polytomy named the trans-Eurasian (TEA) group. The TEA radiation is a dominant evolutionary feature of B. anthracis, with many contemporary populations having resulted from a large spatial dispersal of spores from a single source. Two autopsy specimens from the Sverdlovsk outbreak were deep sequenced to produce draft B. anthracis genomes. This allowed the phylogenetic placement of the Sverdlovsk strain into a clade with two Asian live vaccine strains, including the Russian Tsiankovskii strain. The genome was examined for evidence of drug resistance manipulation or other genetic engineering, but none was found. The Soviet Sverdlovsk strain genome is consistent with a wild-type strain from Russia that had no evidence of genetic manipulation during its industrial production. This work provides insights into the world's largest biological weapons program and provides an extensive B. anthracis phylogenetic reference. IMPORTANCE: The 1979 Russian anthrax outbreak resulted from an industrial accident at the Soviet anthrax spore production facility in the city of Sverdlovsk. Deep genomic sequencing of two autopsy specimens generated a draft genome and phylogenetic placement of the Soviet Sverdlovsk anthrax strain. While it is known that Soviet scientists had genetically manipulated Bacillus anthracis with the potential to evade vaccine prophylaxis and antibiotic therapeutics, there was no genomic evidence of this from the Sverdlovsk production strain genome. The whole-genome SNP genotype of the Sverdlovsk strain was used to precisely identify it and its close relatives in the context of an extensive global B. anthracis strain collection. This genomic identity can now be used for forensic tracking of this weapons material on a global scale and for future anthrax investigations.

High Leptospira Diversity in Animals and Humans Complicates the Search for Common Reservoirs of Human Disease in Rural Ecuador.

BACKGROUND: Leptospirosis is a zoonotic disease responsible for high morbidity around the world, especially in tropical and low income countries. Rats are thought to be the main vector of human leptospirosis in urban settings. However, differences between urban and low-income rural communities provide additional insights into the epidemiology of the disease. METHODOLOGY/PRINCIPAL FINDINGS: Our study was conducted in two low-income rural communities near the coast of Ecuador. We detected and characterized infectious leptospira DNA in a wide variety of samples using new real time quantitative PCR assays and amplicon sequencing. We detected infectious leptospira in a high percentage of febrile patients (14.7%). In contrast to previous studies on leptospirosis risk factors, higher positivity was not found in rats (3.0%) but rather in cows (35.8%) and pigs (21.1%). Six leptospira species were identified (L. borgpetersenii, L kirschnerii, L santarosai, L. interrogans, L noguchii, and an intermediate species within the L. licerasiae and L. wolffii clade) and no significant differences in the species of leptospira present in each animal species was detected (χ2 = 9.89, adj.p-value = 0.27). CONCLUSIONS/SIGNIFICANCE: A large portion of the world's human population lives in low-income, rural communities, however, there is limited information about leptospirosis transmission dynamics in these settings. In these areas, exposure to peridomestic livestock is particularly common and high prevalence of infectious leptospira in cows and pigs suggest that they may be the most important reservoir for human transmission. Genotyping clinical samples show that multiple species of leptospira are involved in human disease. As these genotypes were also detected in samples from a variety of animals, genotype data must be used in conjunction with epidemiological data to provide evidence of transmission and the importance of different potential leptospirosis reservoirs.

The Effects of Signal Erosion and Core Genome Reduction on the Identification of Diagnostic Markers.

UNLABELLED: Whole-genome sequence (WGS) data are commonly used to design diagnostic targets for the identification of bacterial pathogens. To do this effectively, genomics databases must be comprehensive to identify the strict core genome that is specific to the target pathogen. As additional genomes are analyzed, the core genome size is reduced and there is erosion of the target-specific regions due to commonality with related species, potentially resulting in the identification of false positives and/or false negatives. IMPORTANCE: A comparative analysis of 1,130 Burkholderia genomes identified unique markers for many named species, including the human pathogens B. pseudomallei and B. mallei Due to core genome reduction and signature erosion, only 38 targets specific to B. pseudomallei/mallei were identified. By using only public genomes, a larger number of markers were identified, due to undersampling, and this larger number represents the potential for false positives. This analysis has implications for the design of diagnostics for other species where the genomic space of the target and/or closely related species is not well defined.

VNTR diversity in Yersinia pestis isolates from an animal challenge study reveals the potential for in vitro mutations during laboratory cultivation.

Underlying mutation rates and other evolutionary forces shape the population structure of bacteria in nature. Although easily overlooked, similar forces are at work in the laboratory and may influence observed mutations. Here, we investigated tissue samples and Yersinia pestis isolates from a rodent laboratory challenge with strain CO92 using whole genome sequencing and multi-locus variable-number tandem repeat (VNTR) analysis (MLVA). We identified six VNTR mutations that were found to have occurred in vitro during laboratory cultivation rather than in vivo during the rodent challenge. In contrast, no single nucleotide polymorphism (SNP) mutations were observed, either in vivo or in vitro. These results were consistent with previously published mutation rates and the calculated number of Y. pestis generations that occurred during the in vitro versus the in vivo portions of the experiment. When genotyping disease outbreaks, the potential for in vitro mutations should be considered, particularly when highly variable genetic markers such as VNTRs are used.

Rapid Drug Susceptibility Testing of Drug-Resistant Mycobacterium tuberculosis Isolates Directly from Clinical Samples by Use of Amplicon Sequencing: a Proof-of-Concept Study.

Increasingly complex drug-resistant tuberculosis (DR-TB) is a major global health concern and one of the primary reasons why TB is now the leading infectious cause of death worldwide. Rapid characterization of a DR-TB patient's complete drug resistance profile would facilitate individualized treatment in place of empirical treatment, improve treatment outcomes, prevent amplification of resistance, and reduce the transmission of DR-TB. The use of targeted next-generation sequencing (NGS) to obtain drug resistance profiles directly from patient sputum samples has the potential to enable comprehensive evidence-based treatment plans to be implemented quickly, rather than in weeks to months, which is currently needed for phenotypic drug susceptibility testing (DST) results. In this pilot study, we evaluated the performance of amplicon sequencing of Mycobacterium tuberculosis DNA from patient sputum samples using a tabletop NGS technology and automated data analysis to provide a rapid DST solution (the Next Gen-RDST assay). One hundred sixty-six out of 176 (94.3%) sputum samples from the Republic of Moldova yielded complete Next Gen-RDST assay profiles for 7 drugs of interest. We found a high level of concordance of our Next Gen-RDST assay results with phenotypic DST (97.0%) and pyrosequencing (97.8%) results from the same clinical samples. Our Next Gen-RDST assay was also able to estimate the proportion of resistant-to-wild-type alleles down to mixtures of ≤1%, which demonstrates the ability to detect very low levels of resistant variants not detected by pyrosequencing and possibly below the threshold for phenotypic growth methods. The assay as described here could be used as a clinical or surveillance tool.

Local Population Structure and Patterns of Western Hemisphere Dispersal for Coccidioides spp., the Fungal Cause of Valley Fever.

UNLABELLED: Coccidioidomycosis (or valley fever) is a fungal disease with high morbidity and mortality that affects tens of thousands of people each year. This infection is caused by two sibling species, Coccidioides immitis and C. posadasii, which are endemic to specific arid locales throughout the Western Hemisphere, particularly the desert southwest of the United States. Recent epidemiological and population genetic data suggest that the geographic range of coccidioidomycosis is expanding, as new endemic clusters have been identified in the state of Washington, well outside the established endemic range. The genetic mechanisms and epidemiological consequences of this expansion are unknown and require better understanding of the population structure and evolutionary history of these pathogens. Here we performed multiple phylogenetic inference and population genomics analyses of 68 new and 18 previously published genomes. The results provide evidence of substantial population structure in C. posadasii and demonstrate the presence of distinct geographic clades in central and southern Arizona as well as dispersed populations in Texas, Mexico, South America, and Central America. Although a smaller number of C. immitis strains were included in the analyses, some evidence of phylogeographic structure was also detected in this species, which has been historically limited to California and Baja, Mexico. Bayesian analyses indicated that C. posadasii is the more ancient of the two species and that Arizona contains the most diverse subpopulations. We propose a southern Arizona-northern Mexico origin for C. posadasii and describe a pathway for dispersal and distribution out of this region. IMPORTANCE: Coccidioidomycosis, or valley fever, is caused by the pathogenic fungi Coccidioides posadasii and C. immitis The fungal species and disease are primarily found in the American desert southwest, with spotted distribution throughout the Western Hemisphere. Initial molecular studies suggested a likely anthropogenic movement of C. posadasii from North America to South America. Here we comparatively analyze eighty-six genomes of the two Coccidioides species and establish local and species-wide population structures to not only clarify the earlier dispersal hypothesis but also provide evidence of likely ancestral populations and patterns of dispersal for the known subpopulations of C. posadasii.

Genetic variation at the MHC DRB1 locus is similar across Gunnison's prairie dog (Cynomys gunnisoni) colonies regardless of plague history.

Yersinia pestis was introduced to North America around 1900 and leads to nearly 100% mortality in prairie dog (Cynomys spp.) colonies during epizootic events, which suggests this pathogen may exert a strong selective force. We characterized genetic diversity at an MHC class II locus (DRB1) in Gunnison's prairie dog (C. gunnisoni) and quantified population genetic structure at the DRB1 versus 12 microsatellite loci in three large Arizona colonies. Two colonies, Seligman (SE) and Espee Ranch (ES), have experienced multiple plague-related die-offs in recent years, whereas plague has never been documented at Aubrey Valley (AV). We found fairly low allelic diversity at the DRB1 locus, with one allele (DRB1*01) at high frequency (0.67-0.87) in all colonies. Two other DRB1 alleles appear to be trans-species polymorphisms shared with the black-tailed prairie dog (C. ludovicianus), indicating that these alleles have been maintained across evolutionary time frames. Estimates of genetic differentiation were generally lower at the MHC locus (F ST = 0.033) than at microsatellite markers (F ST = 0.098). The reduced differentiation at DRB1 may indicate that selection has been important for shaping variation at MHC loci, regardless of the presence or absence of plague in recent decades. However, genetic drift has probably also influenced the DRB1 locus because its level of differentiation was not different from that of microsatellites in an F ST outlier analysis. We then compared specific MHC alleles to plague survivorship in 60 C. gunnisoni that had been experimentally infected with Y. pestis. We found that survival was greater in individuals that carried at least one copy of the most common allele (DRB1*01) compared to those that did not (60% vs. 20%). Although the sample sizes of these two groups were unbalanced, this result suggests the possibility that this MHC class II locus, or a nearby linked gene, could play a role in plague survival.

NASP: an accurate, rapid method for the identification of SNPs in WGS datasets that supports flexible input and output formats.

Whole-genome sequencing (WGS) of bacterial isolates has become standard practice in many laboratories. Applications for WGS analysis include phylogeography and molecular epidemiology, using single nucleotide polymorphisms (SNPs) as the unit of evolution. NASP was developed as a reproducible method that scales well with the hundreds to thousands of WGS data typically used in comparative genomics applications. In this study, we demonstrate how NASP compares with other tools in the analysis of two real bacterial genomics datasets and one simulated dataset. Our results demonstrate that NASP produces similar, and often better, results in comparison with other pipelines, but is much more flexible in terms of data input types, job management systems, diversity of supported tools and output formats. We also demonstrate differences in results based on the choice of the reference genome and choice of inferring phylogenies from concatenated SNPs or alignments including monomorphic positions. NASP represents a source-available, version-controlled, unit-tested method and can be obtained from tgennorth.github.io/NASP.

Diverse Burkholderia Species Isolated from Soils in the Southern United States with No Evidence of B. pseudomallei.

The global distribution of the soil-dwelling bacterium Burkholderia pseudomallei, causative agent of melioidosis, is poorly understood. We used established culturing methods developed for B. pseudomallei to isolate Burkholderia species from soil collected at 18 sampling sites in three states in the southern United States (Arizona (n = 4), Florida (n = 7), and Louisiana (n = 7)). Using multi-locus sequence typing (MLST) of seven genes, we identified 35 Burkholderia isolates from these soil samples. All species belonged to the B. cepacia complex (Bcc), including B. cenocepacia, B. cepacia, B. contaminans, B. diffusa, B. metallica, B. seminalis, B. vietnamiensis and two unnamed members of the Bcc. The MLST analysis provided a high level of resolution among and within these species. Despite previous clinical cases within the U.S. involving B. pseudomallei and its close phylogenetic relatives, we did not isolate any of these taxa. The Bcc contains a number of opportunistic pathogens that cause infections in cystic fibrosis patients. Interestingly, we found that B. vietnamiensis was present in soil from all three states, suggesting it may be a common component in southern U.S. soils. Most of the Burkholderia isolates collected in this study were from Florida (30/35; 86%), which may be due to the combination of relatively moist, sandy, and acidic soils found there compared to the other two states. We also investigated one MLST gene, recA, for its ability to identify species within Burkholderia. A 365bp fragment of recA recovered nearly the same species-level identification as MLST, thus demonstrating its cost effective utility when conducting environmental surveys for Burkholderia. Although we did not find B. pseudomallei, our findings document that other diverse Burkholderia species are present in soils in the southern United States.

Best practices for evaluating single nucleotide variant calling methods for microbial genomics.

Innovations in sequencing technologies have allowed biologists to make incredible advances in understanding biological systems. As experience grows, researchers increasingly recognize that analyzing the wealth of data provided by these new sequencing platforms requires careful attention to detail for robust results. Thus far, much of the scientific Communit's focus for use in bacterial genomics has been on evaluating genome assembly algorithms and rigorously validating assembly program performance. Missing, however, is a focus on critical evaluation of variant callers for these genomes. Variant calling is essential for comparative genomics as it yields insights into nucleotide-level organismal differences. Variant calling is a multistep process with a host of potential error sources that may lead to incorrect variant calls. Identifying and resolving these incorrect calls is critical for bacterial genomics to advance. The goal of this review is to provide guidance on validating algorithms and pipelines used in variant calling for bacterial genomics. First, we will provide an overview of the variant calling procedures and the potential sources of error associated with the methods. We will then identify appropriate datasets for use in evaluating algorithms and describe statistical methods for evaluating algorithm performance. As variant calling moves from basic research to the applied setting, standardized methods for performance evaluation and reporting are required; it is our hope that this review provides the groundwork for the development of these standards.

Using Whole Genome Analysis to Examine Recombination across Diverse Sequence Types of Staphylococcus aureus.

Staphylococcus aureus is an important clinical pathogen worldwide and understanding this organism's phylogeny and, in particular, the role of recombination, is important both to understand the overall spread of virulent lineages and to characterize outbreaks. To further elucidate the phylogeny of S. aureus, 35 diverse strains were sequenced using whole genome sequencing. In addition, 29 publicly available whole genome sequences were included to create a single nucleotide polymorphism (SNP)-based phylogenetic tree encompassing 11 distinct lineages. All strains of a particular sequence type fell into the same clade with clear groupings of the major clonal complexes of CC8, CC5, CC30, CC45 and CC1. Using a novel analysis method, we plotted the homoplasy density and SNP density across the whole genome and found evidence of recombination throughout the entire chromosome, but when we examined individual clonal lineages we found very little recombination. However, when we analyzed three branches of multiple lineages, we saw intermediate and differing levels of recombination between them. These data demonstrate that in S. aureus, recombination occurs across major lineages that subsequently expand in a clonal manner. Estimated mutation rates for the CC8 and CC5 lineages were different from each other. While the CC8 lineage rate was similar to previous studies, the CC5 lineage was 100-fold greater. Fifty known virulence genes were screened in all genomes in silico to determine their distribution across major clades. Thirty-three genes were present variably across clades, most of which were not constrained by ancestry, indicating horizontal gene transfer or gene loss.

Phylogenetically typing bacterial strains from partial SNP genotypes observed from direct sequencing of clinical specimen metagenomic data.

We describe an approach for genotyping bacterial strains from low coverage genome datasets, including metagenomic data from complex samples. Sequence reads from unknown samples are aligned to a reference genome where the allele states of known SNPs are determined. The Whole Genome Focused Array SNP Typing (WG-FAST) pipeline can identify unknown strains with much less read data than is needed for genome assembly. To test WG-FAST, we resampled SNPs from real samples to understand the relationship between low coverage metagenomic data and accurate phylogenetic placement. WG-FAST can be downloaded from https://github.com/jasonsahl/wgfast.

Detection of Low-Level Mixed-Population Drug Resistance in Mycobacterium tuberculosis Using High Fidelity Amplicon Sequencing.

Undetected and untreated, low-levels of drug resistant (DR) subpopulations in clinical Mycobacterium tuberculosis (Mtb) infections may lead to development of DR-tuberculosis, potentially resulting in treatment failure. Current phenotypic DR susceptibility testing has a theoretical potential for 1% sensitivity, is not quantitative, and requires several weeks to complete. The use of "single molecule-overlapping reads" (SMOR) analysis with next generation DNA sequencing for determination of ultra-rare target alleles in complex mixtures provides increased sensitivity over standard DNA sequencing. Ligation free amplicon sequencing with SMOR analysis enables the detection of resistant allele subpopulations at ≥0.1% of the total Mtb population in near real-time analysis. We describe the method using standardized mixtures of DNA from resistant and susceptible Mtb isolates and the assay's performance for detecting ultra-rare DR subpopulations in DNA extracted directly from clinical sputum samples. SMOR analysis enables rapid near real-time detection and tracking of previously undetectable DR sub-populations in clinical samples allowing for the evaluation of the clinical relevance of low-level DR subpopulations. This will provide insights into interventions aimed at suppressing minor DR subpopulations before they become clinically significant.