Pan-Enterovirus Characterization Reveals Cryptic Circulation of Clinically Relevant Subtypes in Arizona Wastewater.

Background: Most seasonally circulating enteroviruses result in asymptomatic or mildly symptomatic infections. In rare cases, however, infection with some subtypes can result in paralysis or death. Of the 300 subtypes known, only poliovirus is reportable, limiting our understanding of the distribution of other enteroviruses that can cause clinical disease. Objective: The overarching objectives of this study were to: 1) describe the distribution of enteroviruses in Arizona during the late summer and fall of 2022, the time of year when they are thought to be most abundant, and 2) demonstrate the utility of viral pan-assay approaches for semi-agnostic discovery that can be followed up by more targeted assays and phylogenomics. Methods: This study utilizes pooled nasal samples collected from school-aged children and long-term care facility residents, and wastewater from multiple locations in Arizona during July-October of 2022. We used PCR to amplify and sequence a region common to all enteroviruses, followed by species-level bioinformatic characterization using the QIIME 2 platform. For Enterovirus-D68 (EV-D68), detection was carried out using RT-qPCR, followed by confirmation using near-complete whole EV-D68 genome sequencing using a newly designed tiled amplicon approach. Results: In the late summer and early fall of 2022, multiple enterovirus species were identified in Arizona wastewater, with Coxsackievirus A6, EV-D68, and Coxsackievirus A19 composing 86% of the characterized reads sequenced. While EV-D68 was not identified in pooled human nasal samples, and the only reported acute flaccid myelitis case in Arizona did not test positive for the virus, an in-depth analysis of EV-D68 in wastewater revealed that the virus was circulating from August through mid-October. A phylogenetic analysis on this relatively limited dataset revealed just a few importations into the state, with a single clade indicating local circulation. Significance: This study further supports the utility of wastewater-based epidemiology to identify potential public health threats. Our further investigations into EV-D68 shows how these data might help inform healthcare diagnoses for children presenting with concerning neurological symptoms.

Bacterial Genome Wide Association Studies (bGWAS) and Transcriptomics Identifies Cryptic Antimicrobial Resistance Mechanisms in Acinetobacter baumannii.

Antimicrobial resistance (AMR) in the nosocomial pathogen, Acinetobacter baumannii, is becoming a serious public health threat. While some mechanisms of AMR have been reported, understanding novel mechanisms of resistance is critical for identifying emerging resistance. One of the first steps in identifying novel AMR mechanisms is performing genotype/phenotype association studies; however, performing these studies is complicated by the plastic nature of the A. baumannii pan-genome. In this study, we compared the antibiograms of 12 antimicrobials associated with multiple drug families for 84 A. baumannii isolates, many isolated in Arizona, USA. in silico screening of these genomes for known AMR mechanisms failed to identify clear correlations for most drugs. We then performed a bacterial genome wide association study (bGWAS) looking for associations between all possible 21-mers; this approach generally failed to identify mechanisms that explained the resistance phenotype. In order to decrease the genomic noise associated with population stratification, we compared four phylogenetically-related pairs of isolates with differing susceptibility profiles. RNA-Sequencing (RNA-Seq) was performed on paired isolates and differentially-expressed genes were identified. In these isolate pairs, five different potential mechanisms were identified, highlighting the difficulty of broad AMR surveillance in this species. To verify and validate differential expression, amplicon sequencing was performed. These results suggest that a diagnostic platform based on gene expression rather than genomics alone may be beneficial in certain surveillance efforts. The implementation of such advanced diagnostics coupled with increased AMR surveillance will potentially improve A. baumannii infection treatment and patient outcomes.

Genomic analyses of Staphylococcus aureus clonal complex 45 isolates does not distinguish nasal carriage from bacteraemia.

Staphylococcus aureus is a colonizing opportunistic pathogen and a leading cause of bloodstream infection with high morbidity and mortality. S. aureus carriage frequency is reportedly between 20 and 40 % among healthy adults, with S. aureus colonization considered to be a risk factor for S. aureus bacteraemia. It is unknown whether a genetic component of the bacterium is associated with S. aureus bacteraemia in comparison to nasal carriage strains. Previous association studies primarily focusing on the clinical outcome of an S. aureus infection have produced conflicting results, often limited by study design challenged by sample collections and the clonal diversity of S. aureus. To date, no study has investigated whether genomic features separate nasal carriage isolates from S. aureus bacteraemia isolates within a single clonal lineage. Here we have investigated whether genomic features, including single-nucleotide polymorphisms (SNPs), genes, or kmers, distinguish S. aureus nasal carriage isolates from bacteraemia isolates that all belong to the same clonal lineage [clonal complex 45 (CC45)] using whole-genome sequencing (WGS) and a genome-wide association (GWA) approach. From CC45, 100 isolates (50 bacteraemia and 50 nasal carriage, geographically and temporally matched) from Denmark were whole-genome sequenced and subjected to GWA analyses involving gene copy number variation, SNPs, gene content, kmers and gene combinations, while correcting for lineage effects. No statistically significant association involving SNPs, specific genes, gene variants, gene copy number variation, or a combination of genes was identified that could distinguish bacteraemia isolates from nasal carriage isolates. The presented results suggest that all S. aureus nasal CC45 isolates carry the potential to cause invasive disease, as no core or accessory genome content or variations were statistically associated with invasiveness.

Zoonotic Source Attribution of Salmonella enterica Serotype Typhimurium Using Genomic Surveillance Data, United States.

Increasingly, routine surveillance and monitoring of foodborne pathogens using whole-genome sequencing is creating opportunities to study foodborne illness epidemiology beyond routine outbreak investigations and case-control studies. Using a global phylogeny of Salmonella enterica serotype Typhimurium, we found that major livestock sources of the pathogen in the United States can be predicted through whole-genome sequencing data. Relatively steady rates of sequence divergence in livestock lineages enabled the inference of their recent origins. Elevated accumulation of lineage-specific pseudogenes after divergence from generalist populations and possible metabolic acclimation in a representative swine isolate indicates possible emergence of host adaptation. We developed and retrospectively applied a machine learning Random Forest classifier for genomic source prediction of Salmonella Typhimurium that correctly attributed 7 of 8 major zoonotic outbreaks in the United States during 1998-2013. We further identified 50 key genetic features that were sufficient for robust livestock source prediction.

Improved Subtyping of Staphylococcus aureus Clonal Complex 8 Strains Based on Whole-Genome Phylogenetic Analysis.

Strains of Staphylococcus aureus in clonal complex 8 (CC8), including USA300, USA500, and the Iberian clone, are prevalent pathogens in the United States, both inside and outside health care settings. Methods for typing CC8 strains are becoming obsolete as the strains evolve and diversify, and whole-genome sequencing has shown that some strain types fall into multiple sublineages within CC8. In this study, we attempt to clarify the strain nomenclature of CC8, classifying the major strain types based on whole-genome sequence phylogenetics using both methicillin-resistant S. aureus (MRSA) and methicillin-susceptible S. aureus (MSSA) genomes. We show that isolates of the Archaic and Iberian clones from decades ago make up the most basal clade of the main CC8 lineages and that at least one successful lineage of CC8, made up mostly of MSSA, diverged before the other well-known strain types USA500 and USA300. We also show that the USA500 type includes two clades separated by the previously described "Canadian epidemic MRSA" strain CMRSA9, that one clade containing USA500 also contains the USA300 clade, and that the USA300-0114 strain type is not a monophyletic group. Additionally, we present a rapid, simple CC8 strain-typing scheme using real-time PCR assays that target single nucleotide polymorphisms (SNPs) derived from our CC8 phylogeny and show the significant benefit of using more stable genomic markers based on evolutionary lineages over traditional S. aureus typing techniques. This more accurate and accessible S. aureus typing system may improve surveillance and better inform the epidemiology of this very important pathogen.IMPORTANCEStaphylococcus aureus is a major human pathogen worldwide in both community and health care settings. Surveillance for S. aureus strains is important to our understanding of their spread and to informing infection prevention and control. Confusion surrounding the strain nomenclature of one of the most prevalent lineages of S. aureus, clonal complex 8 (CC8), and the imprecision of current tools for typing S. aureus make surveillance and source tracing difficult and sometimes misleading. In this study, we clarify the CC8 strain designations and propose a new typing scheme for CC8 isolates that is rapid and easy to use. This typing scheme is based on relatively stable genomic markers, and we demonstrate its superiority over traditional typing techniques. This scheme has the potential to greatly improve epidemiological investigations of S. aureus.

KlebSeq, a Diagnostic Tool for Surveillance, Detection, and Monitoring of Klebsiella pneumoniae.

Health care-acquired infections (HAIs) kill tens of thousands of people each year and add significantly to health care costs. Multidrug-resistant and epidemic strains are a large proportion of HAI agents, and multidrug-resistant strains of Klebsiella pneumoniae, a leading HAI agent, have caused an urgent public health crisis. In the health care environment, patient colonization by K. pneumoniae precedes infection, and transmission via colonization leads to outbreaks. Periodic patient screening for K. pneumoniae colonization has the potential to curb the number of HAIs. In this report, we describe the design and validation of KlebSeq, a highly informative screening tool that detects Klebsiella species and identifies clinically important strains and characteristics by using highly multiplexed amplicon sequencing without a live-culturing step. We demonstrate the utility of this tool on several complex specimen types, including urine, wound swabs and tissue, and several types of respiratory and fecal specimens, showing K. pneumoniae species and clonal group identification and antimicrobial resistance and virulence profiling, including capsule typing. Use of this amplicon sequencing tool to screen patients for Klebsiella carriage could inform health care staff of the risk of infection and outbreak potential. KlebSeq also serves as a model for next-generation molecular tools for public health and health care, as expansion of this tool can be used for several other HAI agents or applications.

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.

Hypervirulent emm59 Clone in Invasive Group A Streptococcus Outbreak, Southwestern United States.

The hyper-virulent emm59 genotype of invasive group A Streptococcus was identified in northern Arizona in 2015. Eighteen isolates belonging to a genomic cluster grouped most closely with recently identified isolates in New Mexico. The continued transmission of emm59 in the southwestern United States poses a public health concern.

Whole-Genome Sequencing to Determine Origin of Multinational Outbreak of Sarocladium kiliense Bloodstream Infections.

We used whole-genome sequence typing (WGST) to investigate an outbreak of Sarocladium kiliense bloodstream infections (BSI) associated with receipt of contaminated antinausea medication among oncology patients in Colombia and Chile during 2013-2014. Twenty-five outbreak isolates (18 from patients and 7 from medication vials) and 11 control isolates unrelated to this outbreak were subjected to WGST to elucidate a source of infection. All outbreak isolates were nearly indistinguishable (<5 single-nucleotide polymorphisms), and >21,000 single-nucleotide polymorphisms were identified from unrelated control isolates, suggesting a point source for this outbreak. S. kiliense has been previously implicated in healthcare-related infections; however, the lack of available typing methods has precluded the ability to substantiate point sources. WGST for outbreak investigation caused by eukaryotic pathogens without reference genomes or existing genotyping methods enables accurate source identification to guide implementation of appropriate control and prevention measures.

Whole genome SNP typing to investigate methicillin-resistant Staphylococcus aureus carriage in a health-care provider as the source of multiple surgical site infections.

BACKGROUND: Prevention of nosocomial transmission of infections is a central responsibility in the healthcare environment, and accurate identification of transmission events presents the first challenge. Phylogenetic analysis based on whole genome sequencing provides a high-resolution approach for accurately relating isolates to one another, allowing precise identification or exclusion of transmission events and sources for nearly all cases. We sequenced 24 methicillin-resistant Staphylococcus aureus (MRSA) genomes to retrospectively investigate a suspected point source of three surgical site infections (SSIs) that occurred over a one-year period. The source of transmission was believed to be a surgical team member colonized with MRSA, involved in all surgeries preceding the SSI cases, who was subsequently decolonized. Genetic relatedness among isolates was determined using whole genome single nucleotide polymorphism (SNP) data. RESULTS: Whole genome SNP typing (WGST) revealed 283 informative SNPs between the surgical team member's isolate and the closest SSI isolate. The second isolate was 286 and the third was thousands of SNPs different, indicating the nasal carriage strain from the surgical team member was not the source of the SSIs. Given the mutation rates estimated for S. aureus, none of the SSI isolates share a common ancestor within the past 16 years, further discounting any common point source for these infections. The decolonization procedures and resources spent on the point source infection control could have been prevented if WGST was performed at the time of the suspected transmission, instead of retrospectively. CONCLUSIONS: Whole genome sequence analysis is an ideal method to exclude isolates involved in transmission events and nosocomial outbreaks, and coupling this method with epidemiological data can determine if a transmission event occurred. These methods promise to direct infection control resources more appropriately.

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.

Mapping the Evolution of Hypervirulent Klebsiella pneumoniae.

UNLABELLED: Highly invasive, community-acquired Klebsiella pneumoniae infections have recently emerged, resulting in pyogenic liver abscesses. These infections are caused by hypervirulent K. pneumoniae (hvKP) isolates primarily of capsule serotype K1 or K2. Hypervirulent K1 isolates belong to clonal complex 23 (CC23), indicating that this clonal lineage has a specific genetic background conferring hypervirulence. Here, we apply whole-genome sequencing to a collection of K. pneumoniae isolates to characterize the phylogenetic background of hvKP isolates with an emphasis on CC23. Most of the hvKP isolates belonged to CC23 and grouped into a distinct monophyletic clade, revealing that CC23 is a unique clonal lineage, clearly distinct from nonhypervirulent strains. Separate phylogenetic analyses of the CC23 isolates indicated that the CC23 lineage evolved recently by clonal expansion from a single common ancestor. Limited grouping according to geographical origin was observed, suggesting that CC23 has spread globally through multiple international transmissions. Conversely, hypervirulent K2 strains clustered in genetically unrelated groups. Strikingly, homologues of a large virulence plasmid were detected in all hvKP clonal lineages, indicating a key role in K. pneumoniae hypervirulence. The plasmid encodes two siderophores, aerobactin and salmochelin, and RmpA (regulator of the mucoid phenotype); all these factors were found to be restricted to hvKP isolates. Genomic comparisons revealed additional factors specifically associated with CC23. These included a distinct variant of a genomic island encoding yersiniabactin, colibactin, and microcin E492. Furthermore, additional novel genomic regions unique to CC23 were revealed which may also be involved in the increased virulence of this important clonal lineage. IMPORTANCE: During the last 3 decades, hypervirulent Klebsiella pneumoniae (hvKP) isolates have emerged, causing severe community-acquired infections primarily in the form of pyogenic liver abscesses. This syndrome has so far primarily been found in Southeast Asia, but increasing numbers of cases are being reported worldwide, indicating that the syndrome is turning into a globally emerging disease. We applied whole-genome sequencing to a collection of K. pneumoniae clinical isolates to reveal the phylogenetic background of hvKP and to identify genetic factors associated with the increased virulence. The hvKP isolates primarily belonged to clonal complex 23 (CC23), and this clonal lineage was revealed to be clearly distinct from nonhypervirulent strains. A specific virulence plasmid was found to be associated with hypervirulence, and novel genetic determinants uniquely associated with CC23 were identified. Our findings extend the understanding of the genetic background of the emergence of hvKP clones.

Genomic Analysis of the Emergence and Rapid Global Dissemination of the Clonal Group 258 Klebsiella pneumoniae Pandemic.

Multidrug-resistant Klebsiella pneumoniae producing the KPC carbapenemase have rapidly spread throughout the world, causing severe healthcare-associated infections with limited antimicrobial treatment options. Dissemination of KPC-producing K. pneumoniae is largely attributed to expansion of a single dominant strain, ST258. In this study, we explore phylogenetic relationships and evolution within ST258 and its clonal group, CG258, using whole genome sequence analysis of 167 isolates from 20 countries collected over 17 years. Our results show a common ST258 ancestor emerged from its diverse parental clonal group around 1995 and likely acquired blaKPC prior to dissemination. Over the past two decades, ST258 has remained highly clonal despite diversity in accessory elements and divergence in the capsule polysaccharide synthesis locus. Apart from the large recombination event that gave rise to ST258, few mutations set it apart from its clonal group. However, one mutation occurs in a global transcription regulator. Characterization of outer membrane protein sequences revealed a profile in ST258 that includes a truncated OmpK35 and modified OmpK37. Our work illuminates potential genomic contributors to the pathogenic success of ST258, helps us better understand the global dissemination of this strain, and identifies genetic markers unique to ST258.

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.

Valley fever: finding new places for an old disease: Coccidioides immitis found in Washington State soil associated with recent human infection.

We used real-time polymerase chain reaction and culture to demonstrate persistent colonization of soils by Coccidioides immitis, an agent of valley fever, in Washington State linked to recent human infections and located outside the endemic range. Whole-genome sequencing confirmed genetic identity between isolates from soil and one of the case-patients.

Comparative analysis of subtyping methods against a whole-genome-sequencing standard for Salmonella enterica serotype Enteritidis.

A retrospective investigation was performed to evaluate whole-genome sequencing as a benchmark for comparing molecular subtyping methods for Salmonella enterica serotype Enteritidis and survey the population structure of commonly encountered S. enterica serotype Enteritidis outbreak isolates in the United States. A total of 52 S. enterica serotype Enteritidis isolates representing 16 major outbreaks and three sporadic cases collected between 2001 and 2012 were sequenced and subjected to subtyping by four different methods: (i) whole-genome single-nucleotide-polymorphism typing (WGST), (ii) multiple-locus variable-number tandem-repeat (VNTR) analysis (MLVA), (iii) clustered regularly interspaced short palindromic repeats combined with multi-virulence-locus sequence typing (CRISPR-MVLST), and (iv) pulsed-field gel electrophoresis (PFGE). WGST resolved all outbreak clusters and provided useful robust phylogenetic inference results with high epidemiological correlation. While both MLVA and CRISPR-MVLST yielded higher discriminatory power than PFGE, MLVA outperformed the other methods in delineating outbreak clusters whereas CRISPR-MVLST showed the potential to trace major lineages and ecological origins of S. enterica serotype Enteritidis. Our results suggested that whole-genome sequencing makes a viable platform for the evaluation and benchmarking of molecular subtyping methods.

MetaGeniE: characterizing human clinical samples using deep metagenomic sequencing.

With the decreasing cost of next-generation sequencing, deep sequencing of clinical samples provides unique opportunities to understand host-associated microbial communities. Among the primary challenges of clinical metagenomic sequencing is the rapid filtering of human reads to survey for pathogens with high specificity and sensitivity. Metagenomes are inherently variable due to different microbes in the samples and their relative abundance, the size and architecture of genomes, and factors such as target DNA amounts in tissue samples (i.e. human DNA versus pathogen DNA concentration). This variation in metagenomes typically manifests in sequencing datasets as low pathogen abundance, a high number of host reads, and the presence of close relatives and complex microbial communities. In addition to these challenges posed by the composition of metagenomes, high numbers of reads generated from high-throughput deep sequencing pose immense computational challenges. Accurate identification of pathogens is confounded by individual reads mapping to multiple different reference genomes due to gene similarity in different taxa present in the community or close relatives in the reference database. Available global and local sequence aligners also vary in sensitivity, specificity, and speed of detection. The efficiency of detection of pathogens in clinical samples is largely dependent on the desired taxonomic resolution of the organisms. We have developed an efficient strategy that identifies "all against all" relationships between sequencing reads and reference genomes. Our approach allows for scaling to large reference databases and then genome reconstruction by aggregating global and local alignments, thus allowing genetic characterization of pathogens at higher taxonomic resolution. These results were consistent with strain level SNP genotyping and bacterial identification from laboratory culture.

TaqMan real-time PCR assays for single-nucleotide polymorphisms which identify Francisella tularensis and its subspecies and subpopulations.

Francisella tularensis, the etiologic agent of tularemia and a Class A Select Agent, is divided into three subspecies and multiple subpopulations that differ in virulence and geographic distribution. Given these differences, there is a need to rapidly and accurately determine if a strain is F. tularensis and, if it is, assign it to subspecies and subpopulation. We designed TaqMan real-time PCR genotyping assays using eleven single nucleotide polymorphisms (SNPs) that were potentially specific to closely related groups within the genus Francisella, including numerous subpopulations within F. tularensis species. We performed extensive validation studies to test the specificity of these SNPs to particular populations by screening the assays across a set of 565 genetically and geographically diverse F. tularensis isolates and an additional 21 genetic near-neighbor (outgroup) isolates. All eleven assays correctly determined the genetic groups of all 565 F. tularensis isolates. One assay differentiates F. tularensis, F. novicida, and F. hispaniensis from the more genetically distant F. philomiragia and Francisella-like endosymbionts. Another assay differentiates F. tularensis isolates from near neighbors. The remaining nine assays classify F. tularensis-confirmed isolates into F. tularensis subspecies and subpopulations. The genotyping accuracy of these nine assays diminished when tested on outgroup isolates (i.e. non F. tularensis), therefore a hierarchical approach of assay usage is recommended wherein the F. tularensis-specific assay is used before the nine downstream assays. Among F. tularensis isolates, all eleven assays were highly sensitive, consistently amplifying very low concentrations of DNA. Altogether, these eleven TaqMan real-time PCR assays represent a highly accurate, rapid, and sensitive means of identifying the species, subspecies, and subpopulation of any F. tularensis isolate if used in a step-wise hierarchical scheme. These assays would be very useful in clinical, epidemiological, and/or forensic investigations involving F. tularensis.

Origin and evolution of European community-acquired methicillin-resistant Staphylococcus aureus.

UNLABELLED: Community-acquired methicillin-resistant Staphylococcus aureus (CA-MRSA) was recognized in Europe and worldwide in the late 1990s. Within a decade, several genetically and geographically distinct CA-MRSA lineages carrying the small SCCmec type IV and V genetic elements and the Panton-Valentine leukocidin (PVL) emerged around the world. In Europe, the predominant CA-MRSA strain belongs to clonal complex 80 (CC80) and is resistant to kanamycin/amikacin and fusidic acid. CC80 was first reported in 1993 but was relatively rare until the late 1990s. It has since been identified throughout North Africa, the Middle East, and Europe, with recent sporadic reports in sub-Saharan Africa. While strongly associated with skin and soft tissue infections, it is rarely found among asymptomatic carriers. Methicillin-sensitive S. aureus (MSSA) CC80 strains are extremely rare except in sub-Saharan Africa. In the current study, we applied whole-genome sequencing to a global collection of both MSSA and MRSA CC80 isolates. Phylogenetic analyses strongly suggest that the European epidemic CA-MRSA lineage is derived from a PVL-positive MSSA ancestor from sub-Saharan Africa. Moreover, the tree topology suggests a single acquisition of both the SCCmec element and a plasmid encoding the fusidic acid resistance determinant. Four canonical SNPs distinguish the derived CA-MRSA lineage and include a nonsynonymous mutation in accessory gene regulator C (agrC). These changes were associated with a star-like expansion into Europe, the Middle East, and North Africa in the early 1990s, including multiple cases of cross-continent imports likely driven by human migrations. IMPORTANCE: With increasing levels of CA-MRSA reported from most parts of the Western world, there is a great interest in understanding the origin and factors associated with the emergence of these epidemic lineages. To trace the origin, evolution, and dissemination pattern of the European CA-MRSA clone (CC80), we sequenced a global collection of strains of the S. aureus CC80 lineage. Our study determined that a single descendant of a PVL-positive methicillin-sensitive ancestor circulating in sub-Saharan Africa rose to become the dominant CA-MRSA clone in Europe, the Middle East, and North Africa. In the transition from a methicillin-susceptible lineage to a successful CA-MRSA clone, it simultaneously became resistant to fusidic acid, a widely used antibiotic for skin and soft tissue infections, thus demonstrating the importance of antibiotic selection in the success of this clone. This finding furthermore highlights the significance of horizontal gene acquisitions and underscores the combined importance of these factors for the success of CA-MRSA.

Cryptococcus gattii in North American Pacific Northwest: whole-population genome analysis provides insights into species evolution and dispersal.

The emergence of distinct populations of Cryptococcus gattii in the temperate North American Pacific Northwest (PNW) was surprising, as this species was previously thought to be confined to tropical and semitropical regions. Beyond a new habitat niche, the dominant emergent population displayed increased virulence and caused primary pulmonary disease, as opposed to the predominantly neurologic disease seen previously elsewhere. Whole-genome sequencing was performed on 118 C. gattii isolates, including the PNW subtypes and the global diversity of molecular type VGII, to better ascertain the natural source and genomic adaptations leading to the emergence of infection in the PNW. Overall, the VGII population was highly diverse, demonstrating large numbers of mutational and recombinational events; however, the three dominant subtypes from the PNW were of low diversity and were completely clonal. Although strains of VGII were found on at least five continents, all genetic subpopulations were represented or were most closely related to strains from South America. The phylogenetic data are consistent with multiple dispersal events from South America to North America and elsewhere. Numerous gene content differences were identified between the emergent clones and other VGII lineages, including genes potentially related to habitat adaptation, virulence, and pathology. Evidence was also found for possible gene introgression from Cryptococcus neoformans var. grubii that is rarely seen in global C. gattii but that was present in all PNW populations. These findings provide greater understanding of C. gattii evolution in North America and support extensive evolution in, and dispersal from, South America. Importance: Cryptococcus gattii emerged in the temperate North American Pacific Northwest (PNW) in the late 1990s. Beyond a new environmental niche, these emergent populations displayed increased virulence and resulted in a different pattern of clinical disease. In particular, severe pulmonary infections predominated in contrast to presentation with neurologic disease as seen previously elsewhere. We employed population-level whole-genome sequencing and analysis to explore the genetic relationships and gene content of the PNW C. gattii populations. We provide evidence that the PNW strains originated from South America and identified numerous genes potentially related to habitat adaptation, virulence expression, and clinical presentation. Characterization of these genetic features may lead to improved diagnostics and therapies for such fungal infections. The data indicate that there were multiple recent introductions of C. gattii into the PNW. Public health vigilance is warranted for emergence in regions where C. gattii is not thought to be endemic.