Continued improvement in the development of the SARS-CoV-2 whole genome sequencing proficiency testing program.

Application of whole genome sequencing (WGS) has allowed monitoring of the emergence of variants of concern (VOC) of severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) globally. Genomic investigation of emerging variants and surveillance of clinical progress has reduced the public health impact of infection during the COVID-19 pandemic. These steps required developing and implementing a proficiency testing program (PTP), as WGS has been incorporated into routine reference laboratory practice. In this study, we describe how the PTP evaluated the capacity and capability of one New Zealand and 14 Australian public health laboratories to perform WGS of SARS-CoV-2 in 2022. The participants' performances in characterising a specimen panel of known SARS-CoV-2 isolates in the PTP were assessed based on: (1) genome coverage, (2) Pango lineage, and (3) sequence quality, with the choice of assessment metrics refined based on a previously reported assessment conducted in 2021. The participants' performances in 2021 and 2022 were also compared after reassessing the 2021 results using the more stringent metrics adopted in 2022. We found that more participants would have failed the 2021 assessment for all survey samples and a significantly higher fail rate per sample in 2021 compared to 2022. This study highlights the importance of choosing appropriate performance metrics to reflect better the laboratories' capacity to perform SARS-CoV-2 WGS, as was done in the 2022 PTP. It also displays the need for a PTP for WGS of SARS-CoV-2 to be available to public health laboratories ongoing, with continuous refinements in the design and provision of the PTP to account for the dynamic nature of the COVID-19 pandemic as SARS-CoV-2 continues to evolve.

Clinical, microbiological and genomic characterization of Gram-negative bacteria with dual carbapenemases as identified by rapid molecular testing.

Objective: Dual carbapenemase-producing organisms (DCPOs) are an emerging threat that expands the spectrum of antimicrobial resistance. There is limited literature on the clinical and genetic epidemiology of DCPOs. Methods: DCPO isolates were identified by Xpert® Carba-R PCR testing of routine diagnostic cultures performed from 2018 to 2021 at a New York City health system. WGS was performed by Illumina and/or PacBio. Medical records of patients were reviewed for clinical and epidemiological data. Results: Twenty-six DCPO isolates were obtained from 13 patients. Klebsiella pneumoniae (n = 22) was most frequent, followed by Pseudomonas aeruginosa (n = 2), Escherichia coli (n = 1) and Enterobacter cloacae (n = 1). The most common DCPO combination was blaNDM/blaOXA-48-like (n = 16). Notably, 1.05% (24/2290) of carbapenem-resistant Enterobacterales isolates were identified as DCPOs. The susceptibility profiles matched the identified resistance genes, except for a K. pneumoniae (blaKPC/blaOXA-48-like) isolate that was phenotypically susceptible to meropenem. Eleven patients were hospitalized within the year prior to admission, and received antibiotic(s) 1 month prior. Seven patients were originally from outside the USA. Hypertension, kidney disease and diabetes were frequent comorbidities. Death in two cases was attributed to DCPO infection. WGS of eight isolates showed that carbapenemases were located on distinct plasmids, except for one K. pneumoniae isolate where NDM and KPC carbapenemases were located on a single IncC-type plasmid backbone. Conclusions: Here we characterized a series of DCPOs from New York City. Foreign travel, prior hospitalization, antibiotic usage and comorbidities were common among DCPO cases. All carbapenemases were encoded on plasmids, which may facilitate horizontal transfer.

SarS Is a Repressor of Staphylococcus aureus Bicomponent Pore-Forming Leukocidins.

Staphylococcus aureus is a successful pathogen that produces a wide range of virulence factors that it uses to subvert and suppress the immune system. These include the bicomponent pore-forming leukocidins. How the expression of these toxins is regulated is not completely understood. Here, we describe a screen to identify transcription factors involved in the regulation of leukocidins. The most prominent discovery from this screen is that SarS, a known transcription factor which had previously been described as a repressor of alpha-toxin expression, was found to be a potent repressor of leukocidins LukED and LukSF-PV. We found that inactivating sarS resulted in increased virulence both in an ex vivo model using primary human neutrophils and in an in vivo infection model in mice. Further experimentation revealed that SarS represses leukocidins by serving as an activator of Rot, a critical repressor of toxins, as well as by directly binding and repressing the leukocidin promoters. By studying contemporary clinical isolates, we identified naturally occurring mutations in the sarS promoter that resulted in overexpression of sarS and increased repression of leukocidins in USA300 bloodstream clinical isolates. Overall, these data establish SarS as an important repressor of leukocidins and expand our understanding of how these virulence factors are being regulated in vitro and in vivo by S. aureus.

MRSA lineage USA300 isolated from bloodstream infections exhibit altered virulence regulation.

The epidemic community-associated methicillin-resistant Staphylococcus aureus (CA-MRSA) USA300 lineage has recently become a leading cause of hospital-associated bloodstream infections (BSIs). Here, we leveraged this recent introduction into hospitals and the limited genetic variation across USA300 isolates to identify mutations that contribute to its success in a new environment. We found that USA300 BSI isolates exhibit altered virulence regulation. Using comparative genomics to delineate the genes involved in this phenotype, we discovered repeated and independent mutations in the transcriptional regulator sarZ. Mutations in sarZ resulted in increased virulence of USA300 BSI isolates in a murine model of BSI. The sarZ mutations derepressed the expression and production of the surface protein ClfB, which was critical for the pathogenesis of USA300 BSI isolates. Altogether, these findings highlight ongoing evolution of a major MRSA lineage and suggest USA300 strains can optimize their fitness through altered regulation of virulence.

Everybody nose: molecular and clinical characteristics of nasal colonization during active methicillin-resistant Staphylococcus aureus bloodstream infection.

BACKGROUND: Healthcare-associated infections pose a potentially fatal threat to patients worldwide and Staphylococcus aureus is one of the most common causes of healthcare-associated infections. S. aureus is a common commensal pathogen and a frequent cause of bacteremia, with studies demonstrating that nasal and blood isolates from single patients match more than 80% of the time. Here we report on a contemporary collection of colonizing isolates from those with methicillin-resistant S. aureus (MRSA) bloodstream infections to evaluate the diversity within hosts, and detail the clinical features associated with concomitant nasal colonization. METHODS: Swabs of the bilateral anterior nares were obtained from patients diagnosed with MRSA bacteremia. A single colony culture from the blood and an average of 6 colonies from the nares were evaluated for MRSA growth. For the nares cultures, we typed multiple isolates for staphylococcal protein A (spa) and derived the clonal complexes. Demographic and clinical data were obtained retrospectively from the electronic medical record system and analysed using univariate and multivariable regression models. RESULTS: Over an 11-month period, 68 patients were diagnosed with MRSA bloodstream infection, 53 were swabbed, and 37 (70%) were colonized with MRSA in the anterior nares. We performed molecular typing on 213 nasal colonies. Spa types and clonal complexes found in the blood were also detected in the nares in 95% of the cases. We also found that 11% of patients carried more than one clone of MRSA in the nares. Male sex and history of prior hospitalization within the past 90 days increased odds for MRSA colonization. CONCLUSION: The molecular epidemiological landscape of colonization in the setting of invasive disease is diverse and defining the interplay between colonization and invasive disease is critical to combating invasive MRSA disease.

Modified methicillin-resistant Staphylococcus aureus detected in neonatal intensive care patients.

OBJECTIVES: As part of an active MRSA surveillance programme in our neonatal ICU, we identified nares surveillance cultures from two infants that displayed heterogeneity in methicillin resistance between isolated subclones that lacked mecA and mecC. METHODS: The underlying mechanism for the modified Staphylococcus aureus (MODSA) methicillin-resistance phenotype was investigated by WGS. RESULTS: Comparison of finished-quality genomes of four MODSA and four MSSA subclones demonstrated that the resistance changes were associated with unique truncating mutations in the gene encoding the cyclic diadenosine monophosphate phosphodiesterase enzyme GdpP or a non-synonymous substitution in the gene encoding PBP2. CONCLUSIONS: These two cases highlight the difficulty in identifying non-mecA, non-mecC-mediated MRSA isolates in the clinical microbiology laboratory, which leads to difficulties in implementing appropriate therapy and infection control measures.

Molecular evidence of SARS-CoV-2 in New York before the first pandemic wave.

Numerous reports document the spread of SARS-CoV-2, but there is limited information on its introduction before the identification of a local case. This may lead to incorrect assumptions when modeling viral origins and transmission. Here, we utilize a sample pooling strategy to screen for previously undetected SARS-CoV-2 in de-identified, respiratory pathogen-negative nasopharyngeal specimens from 3,040 patients across the Mount Sinai Health System in New York. The patients had been previously evaluated for respiratory symptoms or influenza-like illness during the first 10 weeks of 2020. We identify SARS-CoV-2 RNA from specimens collected as early as 25 January 2020, and complete SARS-CoV-2 genome sequences from multiple pools of samples collected between late February and early March, documenting an increase prior to the later surge. Our results provide evidence of sporadic SARS-CoV-2 infections a full month before both the first officially documented case and emergence of New York as a COVID-19 epicenter in March 2020.

PathoSPOT genomic epidemiology reveals under-the-radar nosocomial outbreaks.

BACKGROUND: Whole-genome sequencing (WGS) is increasingly used to map the spread of bacterial and viral pathogens in nosocomial settings. A limiting factor for more widespread adoption of WGS for hospital infection prevention practices is the availability of standardized tools for genomic epidemiology. METHODS: We developed the Pathogen Sequencing Phylogenomic Outbreak Toolkit (PathoSPOT) to automate integration of genomic and medical record data for rapid detection and tracing of nosocomial outbreaks. To demonstrate its capabilities, we applied PathoSPOT to complete genome surveillance data of 197 MRSA bacteremia cases from two hospitals during a 2-year period. RESULTS: PathoSPOT identified 8 clonal clusters encompassing 33 patients (16.8% of cases), none of which had been recognized by standard practices. The largest cluster corresponded to a prolonged outbreak of a hospital-associated MRSA clone among 16 adults, spanning 9 wards over a period of 21 months. Analysis of precise timeline and location data with our toolkit suggested that an initial exposure event in a single ward led to infection and long-term colonization of multiple patients, followed by transmissions to other patients during recurrent hospitalizations. CONCLUSIONS: We demonstrate that PathoSPOT genomic surveillance enables the detection of complex transmission chains that are not readily apparent from epidemiological data and that contribute significantly to morbidity and mortality, enabling more effective intervention strategies.

Introductions and early spread of SARS-CoV-2 in the New York City area.

New York City (NYC) has emerged as one of the epicenters of the current severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic. To identify the early transmission events underlying the rapid spread of the virus in the NYC metropolitan area, we sequenced the virus that causes coronavirus disease 2019 (COVID-19) in patients seeking care at the Mount Sinai Health System. Phylogenetic analysis of 84 distinct SARS-CoV-2 genomes indicates multiple, independent, but isolated introductions mainly from Europe and other parts of the United States. Moreover, we found evidence for community transmission of SARS-CoV-2 as suggested by clusters of related viruses found in patients living in different neighborhoods of the city.

Blurred Molecular Epidemiological Lines Between the Two Dominant Methicillin-Resistant Staphylococcus aureus Clones.

BACKGROUND: Methicillin-resistant Staphylococcus aureus (MRSA) causes life-threatening infections in both community and hospital settings and is a leading cause of health care-associated infections (HAIs). We sought to describe the molecular epidemiological landscape of patients with MRSA bloodstream infections (BSIs) at an urban medical center by evaluating the clinical characteristics associated with the two dominant endemic clones. METHODS: Comprehensive clinical data from the electronic health records of 227 hospitalized patients ≥18 years old with MRSA BSI over a 33-month period in New York City were collected. The descriptive epidemiology and mortality associated with the two dominant clones were compared using logistic regression. RESULTS: Molecular analysis revealed that 91% of all single-patient MRSA BSIs were due to two equally represented genotypes, clonal complex (CC) 5 (n = 117) and CC8 (n = 110). MRSA BSIs were associated with a 90-day mortality rate of 27%. CC8 caused disease more frequently in younger age groups (56 ± 17 vs 67 ± 17 years old; P < .001) and in those of nonwhite race (odds ratio [OR], 3.45; 95% confidence interval [CI], 1.51-7.87; P = .003), with few other major distinguishing features. Morbidity and mortality also did not differ significantly between the two clones. CC8 caused BSIs more frequently in the setting of peripheral intravenous catheters (OR, 5.96; 95% CI, 1.51-23.50; P = .01). CONCLUSIONS: The clinical features distinguishing dominant MRSA clones continue to converge. The association of CC8 with peripheral intravenous catheter infections underscores the importance of classical community clones causing hospital-onset infections. Ongoing monitoring and analysis of the dynamic epidemiology of this endemic pathogen are crucial to inform management and prevent disease.

A Complete Genome Screening Program of Clinical Methicillin-Resistant Staphylococcus aureus Isolates Identifies the Origin and Progression of a Neonatal Intensive Care Unit Outbreak.

Whole-genome sequencing (WGS) of Staphylococcus aureus is increasingly used as part of infection prevention practices. In this study, we established a long-read technology-based WGS screening program of all first-episode methicillin-resistant Staphylococcus aureus (MRSA) blood infections at a major urban hospital. A survey of 132 MRSA genomes assembled from long reads enabled detailed characterization of an outbreak lasting several months of a CC5/ST105/USA100 clone among 18 infants in a neonatal intensive care unit (NICU). Available hospital-wide genome surveillance data traced the origins of the outbreak to three patients admitted to adult wards during a 4-month period preceding the NICU outbreak. The pattern of changes among complete outbreak genomes provided full spatiotemporal resolution of its progression, which was characterized by multiple subtransmissions and likely precipitated by equipment sharing between adults and infants. Compared to other hospital strains, the outbreak strain carried distinct mutations and accessory genetic elements that impacted genes with roles in metabolism, resistance, and persistence. This included a DNA recognition domain recombination in the hsdS gene of a type I restriction modification system that altered DNA methylation. Transcriptome sequencing (RNA-Seq) profiling showed that the (epi)genetic changes in the outbreak clone attenuated agr gene expression and upregulated genes involved in stress response and biofilm formation. Overall, our findings demonstrate the utility of long-read sequencing for hospital surveillance and for characterizing accessory genomic elements that may impact MRSA virulence and persistence.

The purine biosynthesis regulator PurR moonlights as a virulence regulator in Staphylococcus aureus.

The pathogen Staphylococcus aureus colonizes and infects a variety of different sites within the human body. To adapt to these different environments, S. aureus relies on a complex and finely tuned regulatory network. While some of these networks have been well-elucidated, the functions of more than 50% of the transcriptional regulators in S. aureus remain unexplored. Here, we assess the contribution of the LacI family of metabolic regulators to staphylococcal virulence. We found that inactivating the purine biosynthesis regulator purR resulted in a strain that was acutely virulent in bloodstream infection models in mice and in ex vivo models using primary human neutrophils. Remarkably, these enhanced pathogenic traits are independent of purine biosynthesis, as the purR mutant was still highly virulent in the presence of mutations that disrupt PurR's canonical role. Through the use of transcriptomics coupled with proteomics, we revealed that a number of virulence factors are differentially regulated in the absence of purR Indeed, we demonstrate that PurR directly binds to the promoters of genes encoding virulence factors and to master regulators of virulence. These results guided us into further ex vivo and in vivo studies, where we discovered that S. aureus toxins drive the death of human phagocytes and mice, whereas the surface adhesin FnbA contributes to the increased bacterial burden observed in the purR mutant. Thus, S. aureus repurposes a metabolic regulator to directly control the expression of virulence factors, and by doing so, tempers its pathogenesis.

A physical and genetic map of Cannabis sativa identifies extensive rearrangements at the THC/CBD acid synthase loci.

Cannabis sativa is widely cultivated for medicinal, food, industrial, and recreational use, but much remains unknown regarding its genetics, including the molecular determinants of cannabinoid content. Here, we describe a combined physical and genetic map derived from a cross between the drug-type strain Purple Kush and the hemp variety "Finola." The map reveals that cannabinoid biosynthesis genes are generally unlinked but that aromatic prenyltransferase (AP), which produces the substrate for THCA and CBDA synthases (THCAS and CBDAS), is tightly linked to a known marker for total cannabinoid content. We further identify the gene encoding CBCA synthase (CBCAS) and characterize its catalytic activity, providing insight into how cannabinoid diversity arises in cannabis. THCAS and CBDAS (which determine the drug vs. hemp chemotype) are contained within large (>250 kb) retrotransposon-rich regions that are highly nonhomologous between drug- and hemp-type alleles and are furthermore embedded within ∼40 Mb of minimally recombining repetitive DNA. The chromosome structures are similar to those in grains such as wheat, with recombination focused in gene-rich, repeat-depleted regions near chromosome ends. The physical and genetic map should facilitate further dissection of genetic and molecular mechanisms in this commercially and medically important plant.

Genome Plasticity of agr-Defective Staphylococcus aureus during Clinical Infection.

Therapy for bacteremia caused by Staphylococcus aureus is often ineffective, even when treatment conditions are optimal according to experimental protocols. Adapted subclones, such as those bearing mutations that attenuate agr-mediated virulence activation, are associated with persistent infection and patient mortality. To identify additional alterations in agr-defective mutants, we sequenced and assembled the complete genomes of clone pairs from colonizing and infected sites of several patients in whom S. aureus demonstrated a within-host loss of agr function. We report that events associated with agr inactivation result in agr-defective blood and nares strain pairs that are enriched in mutations compared to pairs from wild-type controls. The random distribution of mutations between colonizing and infecting strains from the same patient, and between strains from different patients, suggests that much of the genetic complexity of agr-defective strains results from prolonged infection or therapy-induced stress. However, in one of the agr-defective infecting strains, multiple genetic changes resulted in increased virulence in a murine model of bloodstream infection, bypassing the mutation of agr and raising the possibility that some changes were selected. Expression profiling correlated the elevated virulence of this agr-defective mutant to restored expression of the agr-regulated ESAT6-like type VII secretion system, a known virulence factor. Thus, additional mutations outside the agr locus can contribute to diversification and adaptation during infection by S. aureus agr mutants associated with poor patient outcomes.

Genetic Basis of Emerging Vancomycin, Linezolid, and Daptomycin Heteroresistance in a Case of Persistent Enterococcus faecium Bacteremia.

Whole-genome sequencing was used to examine a persistent Enterococcus faecium bacteremia that acquired heteroresistance to three antibiotics in response to prolonged multidrug therapy. A comparison of the complete genomes before and after each change revealed the emergence of known resistance determinants for vancomycin and linezolid and suggested that a novel mutation in fabF, encoding a fatty acid synthase, was responsible for daptomycin nonsusceptibility. Plasmid recombination contributed to the progressive loss of vancomycin resistance after withdrawal of the drug.

HapFlow: visualizing haplotypes in sequencing data.

SUMMARY: HapFlow is a python application for visualizing haplotypes present in sequencing data. It identifies variant profiles present and reads and creates an abstract visual representation of these profiles to make haplotypes easier to identify. AVAILABILITY AND IMPLEMENTATION: HapFlow is freely available (under a GPL license) for download (for Mac OS X, Unix and Microsoft Windows) from github (http://mjsull.github.io/HapFlow). CONTACT: apolking@usc.edu.au.

Culture-independent genome sequencing of clinical samples reveals an unexpected heterogeneity of infections by Chlamydia pecorum.

Chlamydia pecorum is an important global pathogen of livestock, and it is also a significant threat to the long-term survival of Australia's koala populations. This study employed a culture-independent DNA capture approach to sequence C. pecorum genomes directly from clinical swab samples collected from koalas with chlamydial disease as well as from sheep with arthritis and conjunctivitis. Investigations into single-nucleotide polymorphisms within each of the swab samples revealed that a portion of the reads in each sample belonged to separate C. pecorum strains, suggesting that all of the clinical samples analyzed contained mixed populations of genetically distinct C. pecorum isolates. This observation was independent of the anatomical site sampled and the host species. Using the genomes of strains identified in each of these samples, whole-genome phylogenetic analysis revealed that a clade containing a bovine and a koala isolate is distinct from other clades comprised of livestock or koala C. pecorum strains. Providing additional evidence to support exposure of koalas to Australian livestock strains, two minor strains assembled from the koala swab samples clustered with livestock strains rather than koala strains. Culture-independent probe-based genome capture and sequencing of clinical samples provides the strongest evidence yet to suggest that naturally occurring chlamydial infections are comprised of multiple genetically distinct strains.

Molecular analysis of asymptomatic bacteriuria Escherichia coli strain VR50 reveals adaptation to the urinary tract by gene acquisition.

Urinary tract infections (UTIs) are among the most common infectious diseases of humans, with Escherichia coli responsible for >80% of all cases. One extreme of UTI is asymptomatic bacteriuria (ABU), which occurs as an asymptomatic carrier state that resembles commensalism. To understand the evolution and molecular mechanisms that underpin ABU, the genome of the ABU E. coli strain VR50 was sequenced. Analysis of the complete genome indicated that it most resembles E. coli K-12, with the addition of a 94-kb genomic island (GI-VR50-pheV), eight prophages, and multiple plasmids. GI-VR50-pheV has a mosaic structure and contains genes encoding a number of UTI-associated virulence factors, namely, Afa (afimbrial adhesin), two autotransporter proteins (Ag43 and Sat), and aerobactin. We demonstrated that the presence of this island in VR50 confers its ability to colonize the murine bladder, as a VR50 mutant with GI-VR50-pheV deleted was attenuated in a mouse model of UTI in vivo. We established that Afa is the island-encoded factor responsible for this phenotype using two independent deletion (Afa operon and AfaE adhesin) mutants. E. coli VR50afa and VR50afaE displayed significantly decreased ability to adhere to human bladder epithelial cells. In the mouse model of UTI, VR50afa and VR50afaE displayed reduced bladder colonization compared to wild-type VR50, similar to the colonization level of the GI-VR50-pheV mutant. Our study suggests that E. coli VR50 is a commensal-like strain that has acquired fitness factors that facilitate colonization of the human bladder.

Genomic evolution of the pathogenic Wolbachia strain, wMelPop.

Most strains of the widespread endosymbiotic bacterium Wolbachia pipientis are benign or behave as reproductive parasites. The pathogenic strain wMelPop is a striking exception, however: it overreplicates in its insect hosts and causes severe life shortening. The mechanism of this pathogenesis is currently unknown. We have sequenced the genomes of three variants of wMelPop and of the closely related nonpathogenic strain wMelCS. We show that the genomes of wMelCS and wMelPop appear to be identical in the nonrepeat regions of the genome and differ detectably only by the triplication of a 19-kb region that is unlikely to be associated with life shortening, demonstrating that dramatic differences in the host phenotype caused by this endosymbiont may be the result of only minor genetic changes. We also compare the genomes of the original wMelPop strain from Drosophila melanogaster and two sequential derivatives, wMelPop-CLA and wMelPop-PGYP. To develop wMelPop as a novel biocontrol agent, it was first transinfected into and passaged in mosquito cell lines for approximately 3.5 years, generating wMelPop-CLA. This cell line-passaged strain was then transinfected into Aedes aegypti mosquitoes, creating wMelPop-PGYP, which was sequenced after 4 years in the insect host. We observe a rapid burst of genomic changes during cell line passaging, but no further mutations were detected after transinfection into mosquitoes, indicating either that host preadaptation had occurred in cell lines, that cell lines are a more selectively permissive environment than animal hosts, or both. Our results provide valuable data on the rates of genomic and phenotypic change in Wolbachia associated with host shifts over short time scales.

Easyfig: a genome comparison visualizer.

UNLABELLED: Easyfig is a Python application for creating linear comparison figures of multiple genomic loci with an easy-to-use graphical user interface. BLAST comparisons between multiple genomic regions, ranging from single genes to whole prokaryote chromosomes, can be generated, visualized and interactively coloured, enabling a rapid transition between analysis and the preparation of publication quality figures. AVAILABILITY: Easyfig is freely available (under a GPL license) for download (for Mac OS X, Unix and Microsoft Windows) from the SourceForge web site: http://easyfig.sourceforge.net/.