Improving Bacterial Metagenomic Research through Long-Read Sequencing.

Metagenomic sequencing analysis is central to investigating microbial communities in clinical and environmental studies. Short-read sequencing remains the primary approach for metagenomic research; however, long-read sequencing may offer advantages of improved metagenomic assembly and resolved taxonomic identification. To compare the relative performance for metagenomic studies, we simulated short- and long-read datasets using increasingly complex metagenomes comprising 10, 20, and 50 microbial taxa. Additionally, we used an empirical dataset of paired short- and long-read data generated from mouse fecal pellets to assess real-world performance. We compared metagenomic assembly quality, taxonomic classification, and metagenome-assembled genome (MAG) recovery rates. We show that long-read sequencing data significantly improve taxonomic classification and assembly quality. Metagenomic assemblies using simulated long reads were more complete and more contiguous with higher rates of MAG recovery. This resulted in more precise taxonomic classifications. Principal component analysis of empirical data demonstrated that sequencing technology affects compositional results as samples clustered by sequence type, not sample type. Overall, we highlight strengths of long-read metagenomic sequencing for microbiome studies, including improving the accuracy of classification and relative abundance estimates. These results will aid researchers when considering which sequencing approaches to use for metagenomic projects.

Streptococcus pneumoniae serotype 3 population structure in the era of conjugate vaccines, 2001-2018.

Background. Despite use of highly effective conjugate vaccines, invasive pneumococcal disease (IPD) remains a leading cause of morbidity and mortality and disproportionately affects Indigenous populations. Although included in the 13-valent pneumococcal conjugate vaccine (PCV13), which was introduced in 2010, serotype 3 continues to cause disease among Indigenous communities in the Southwest USA. In the Navajo Nation, serotype 3 IPD incidence increased among adults (3.8/100 000 in 2001-2009 and 6.2/100 000 in 2011-2019); in children the disease persisted although the rates dropped from 5.8/100 000 to 2.3/100 000.Methods. We analysed the genomic epidemiology of serotype 3 isolates collected from 129 adults and 63 children with pneumococcal carriage (n=61) or IPD (n=131) from 2001 to 2018 of the Navajo Nation. Using whole-genome sequencing data, we determined clade membership and assessed changes in serotype 3 population structure over time.Results. The serotype 3 population structure was characterized by three dominant subpopulations: clade II (n=90, 46.9 %) and clade Iα (n=59, 30.7 %), which fall into Clonal Complex (CC) 180, and a non-CC180 clade (n=43, 22.4 %). The proportion of clade II-associated IPD cases increased significantly from 2001 to 2010 to 2011-2018 among adults (23.1-71.8 %; P<0.001) but not in children (27.3-33.3 %; P=0.84). Over the same period, the proportion of clade II-associated carriage increased; this was statistically significant among children (23.3-52.6 %; P=0.04) but not adults (0-50.0 %, P=0.08).Conclusions. In this setting with persistent serotype 3 IPD and carriage, clade II has increased since 2010. Genomic changes may be contributing to the observed trends in serotype 3 carriage and disease over time.

Environmental surface monitoring as a noninvasive method for SARS-CoV-2 surveillance in community settings: Lessons from a university campus study.

Environmental testing of high-touch objects is a potential noninvasive approach for monitoring population-level trends of SARS-CoV-2 and other respiratory viruses within a defined setting. We aimed to determine the association between SARS-CoV-2 contamination on high-touch environmental surfaces, community level case incidence, and university student health data. Environmental swabs were collected from January 2022 to November 2022 from high-touch objects and surfaces from five locations on a large university campus in Florida, USA. RT-qPCR was used to detect and quantify viral RNA, and a subset of positive samples was analyzed by viral genome sequencing to identify circulating lineages. During the study period, we detected SARS-CoV-2 viral RNA on 90.7 % of 162 tested samples. Levels of environmental viral RNA correlated with trends in community-level activity and case reports from the student health center. A significant positive correlation was observed between the estimated viral gene copy number in environmental samples and the weekly confirmed cases at the university. Viral sequencing data from environmental samples identified lineages concurrently circulating in the local community and state based on genomic surveillance data. Further, we detected emerging variants in environmental samples prior to their identification by clinical genomic surveillance. Our results demonstrate the utility of viral monitoring on high-touch environmental surfaces for SARS-CoV-2 surveillance at a community level. In communities with delayed or limited testing facilities, immediate environmental surface testing may considerably inform epidemic dynamics.

Draft genome sequence of a highly proteolytic Staphylococcus aureus USA300 isolate from human urine.

The secreted proteases of Staphylococcus aureus have been shown to be critical during infection. Here, we present the draft genome sequence of S. aureus TGH337, a hyper-proteolytic USA300 strain isolated from human urine.

Joining Forces against Antibiotic Resistance: The One Health Solution.

Antibiotic resistance is a significant global health concern that affects both human and animal populations. The One Health approach acknowledges the interconnectedness of human health, animal health, and the environment. It emphasizes the importance of collaboration and coordination across these sectors to tackle complex health challenges such as antibiotic resistance. In the context of One Health, antibiotic resistance refers to the ability of bacteria to withstand the efficacy of antibiotics, rendering them less effective or completely ineffective in treating infections. The emergence and spread of antibiotic-resistant bacteria pose a threat to human and animal health, as well as to the effectiveness of medical treatments and veterinary interventions. In particular, One Health recognizes that antibiotic use in human medicine, animal agriculture, and the environment are interconnected factors contributing to the development and spread of antibiotic resistance. For example, the misuse and overuse of antibiotics in human healthcare, including inappropriate prescribing and patient non-compliance, can contribute to the selection and spread of resistant bacteria. Similarly, the use of antibiotics in livestock production for growth promotion and disease prevention can contribute to the development of antibiotic resistance in animals and subsequent transmission to humans through the food chain. Addressing antibiotic resistance requires a collaborative One Health approach that involves multiple participants, including healthcare professionals, veterinarians, researchers, and policymakers.

A collaborative approach to improving representation in viral genomic surveillance.

The lack of routine viral genomic surveillance delayed the initial detection of SARS-CoV-2, allowing the virus to spread unfettered at the outset of the U.S. epidemic. Over subsequent months, poor surveillance enabled variants to emerge unnoticed. Against this backdrop, long-standing social and racial inequities have contributed to a greater burden of cases and deaths among minority groups. To begin to address these problems, we developed a new variant surveillance model geared toward building 'next generation' genome sequencing capacity at universities in or near rural areas and engaging the participation of their local communities. The resulting genomic surveillance network has generated more than 1,000 SARS-CoV-2 genomes to date, including the first confirmed case in northeast Louisiana of Omicron, and the first and sixth confirmed cases in Georgia of the emergent BA.2.75 and BQ.1.1 variants, respectively. In agreement with other studies, significantly higher viral gene copy numbers were observed in Delta variant samples compared to those from Omicron BA.1 variant infections, and lower copy numbers were seen in asymptomatic infections relative to symptomatic ones. Collectively, the results and outcomes from our collaborative work demonstrate that establishing genomic surveillance capacity at smaller academic institutions in rural areas and fostering relationships between academic teams and local health clinics represent a robust pathway to improve pandemic readiness.

Antigen test swabs are comparable to nasopharyngeal swabs for sequencing of SARS-CoV-2.

Viral genomic surveillance has been integral in the global response to the SARS-CoV-2 pandemic. Surveillance efforts rely on the availability of representative clinical specimens from ongoing testing activities. However, testing practices have recently shifted due to the widespread availability and use of rapid antigen tests, which could lead to gaps in future monitoring efforts. As such, genomic surveillance strategies must adapt to include laboratory workflows that are robust to sample type. To that end, we compare the results of RT-qPCR and viral genome sequencing using samples from positive BinaxNOW COVID-19 Antigen Card swabs (N = 555) to those obtained from nasopharyngeal (NP) swabs used for nucleic acid amplification testing (N = 135). We show that swabs obtained from antigen cards are comparable in performance to samples from NP swabs, providing a viable alternative and allowing for the potential expansion of viral genomic surveillance to outpatient clinic as well as other settings where rapid antigen tests are often used.

Favipiravir Suppresses Zika Virus (ZIKV) through Activity as a Mutagen.

In a companion paper, we demonstrated that the nucleoside analogue favipiravir (FAV) suppressed Zika virus (ZIKV) replication in three human-derived cell lines-HeLa, SK-N-MC, and HUH-7. Our results revealed that FAV's effect was most pronounced in HeLa cells. In this work, we aimed to explain variation in FAV activity, investigating its mechanism of action and characterizing host cell factors relevant to tissue-specific differences in drug effect. Using viral genome sequencing, we show that FAV therapy was associated with an increase in the number of mutations and promoted the production of defective viral particles in all three cell lines. Our findings demonstrate that defective viral particles made up a larger portion of the viral population released from HeLa cells both at increasing FAV concentrations and at increasing exposure times. Taken together, our companion papers show that FAV acts via lethal mutagenesis against ZIKV and highlight the host cell's influence on the activation and antiviral activity of nucleoside analogues. Furthermore, the information gleaned from these companion papers can be applied to gain a more comprehensive understanding of the activity of nucleoside analogues and the impact of host cell factors against other viral infections for which we currently have no approved antiviral therapies.

A collaborative approach to improve representation in viral genomic surveillance.

The lack of routine viral genomic surveillance delayed the initial detection of SARS-CoV-2, allowing the virus to spread unfettered at the outset of the U.S. epidemic. Over subsequent months, poor surveillance enabled variants to emerge unnoticed. Against this backdrop, long-standing social and racial inequities have contributed to a greater burden of cases and deaths among minority groups. To begin to address these problems, we developed a new variant surveillance model geared toward building microbial genome sequencing capacity at universities in or near rural areas and engaging the participation of their local communities. The resulting genomic surveillance network has generated more than 1,000 SARS-CoV-2 genomes to date, including the first confirmed case in northeast Louisiana of Omicron, and the first and sixth confirmed cases in Georgia of the emergent BA.2.75 and BQ.1.1 variants, respectively. In agreement with other studies, significantly higher viral gene copy numbers were observed in Delta variant samples compared to those from Omicron BA.1 variant infections, and lower copy numbers were seen in asymptomatic infections relative to symptomatic ones. Collectively, the results and outcomes from our collaborative work demonstrate that establishing genomic surveillance capacity at smaller academic institutions in rural areas and fostering relationships between academic teams and local health clinics represent a robust pathway to improve pandemic readiness. Author summary: Genomic surveillance involves decoding a pathogen’s genetic code to track its spread and evolution. During the pandemic, genomic surveillance programs around the world provided valuable data to scientists, doctors, and public health officials. Knowing the complete SARS-CoV-2 genome has helped detect the emergence of new variants, including ones that are more transmissible or cause more severe disease, and has supported the development of diagnostics, vaccines, and therapeutics. The impact of genomic surveillance on public health depends on representative sampling that accurately reflects the diversity and distribution of populations, as well as rapid turnaround time from sampling to data sharing. After a slow start, SARS-CoV-2 genomic surveillance in the United States grew exponentially. Despite this, many rural regions and ethnic minorities remain poorly represented, leaving significant gaps in the data that informs public health responses. To address this problem, we formed a network of universities and clinics in Louisiana, Georgia, and Mississippi with the goal of increasing SARS-CoV-2 sequencing volume, representation, and equity. Our results demonstrate the advantages of rapidly sequencing pathogens in the same communities where the cases occur and present a model that leverages existing academic and clinical infrastructure for a powerful decentralized genomic surveillance system.

Complete Genome Sequences of Nine Streptococcus pneumoniae Serotype 3 Clonal Complex 180 Strains.

We announce the complete genomes of nine Streptococcus pneumoniae strains belonging to serotype 3 clonal complex 180 (CC180). The genomes consist of a single circularized contig with an average length of 2.033 Mbp. Pangenome analysis identified 1,762 core genes and 412 accessory genes. These genomes are the basis for future population genomic studies.

Carriage prevalence and genomic epidemiology of Staphylococcus aureus among Native American children and adults in the Southwestern USA.

Native American individuals in the Southwestern USA experience a higher burden of invasive Staphylococcus aureus disease than the general population. However, little is known about S. aureus carriage in these communities. A cross-sectional study was conducted to determine the carriage prevalence, risk factors and genomic epidemiology of S. aureus among Native American children (<5 years, n=121) and adults (≥18 years, n=167) in the Southwestern USA. Short- and long-read sequencing data were generated using Illumina and Oxford Nanopore Technology platforms to produce high-quality hybrid assemblies, and antibiotic-resistance, virulence and pangenome analyses were performed. S. aureus carriage prevalence was 20.7 % among children, 30.2 % among adults 18-64 years and 16.7 % among adults ≥65 years. Risk factors among adults included recent surgery, prior S. aureus infection among household members, and recent use of gyms or locker rooms by household members. No risk factors were identified among children. The bacterial population structure was dominated by clonal complex 1 (CC1) (21.1 %), CC5 (22.2 %) and CC8 (22.2 %). Isolates from children and adults were intermixed throughout the phylogeny. While the S. aureus population was diverse, the carriage prevalence was comparable to that in the general USA population. Genomic and risk-factor data suggest household, community and healthcare transmission are important components of the local epidemiology.

Meaningful Use of Pathogen Genomic Data.

Population genomic analysis is a powerful tool to understand the evolutionary history of pathogens and the factors contributing to the success or failure of lineages. These studies have significant implications for human health, as evident from our ongoing tracking of SARS-CoV-2. In their article, Gill et al. (J. L. Gill, J. Hedge, D. J. Wilson, and R. C. MacLean, mBio 12:e02168-21, 2021, https://doi.org/10.1128/mBio.02168-21) demonstrate the utility of pathogen genomic data by comprehensively elucidating the origin of methicillin-resistant Staphylococcus aureus ST239. To accomplish this, they leveraged newly developed tools for querying large genomic data sets. Overall, these analyses rely on the availability of representative genomic data along with their associated metadata-information about where and when samples were collected, clinical and epidemiological characteristics, and phenotypic properties. However, in many instances, these data are missing. Here, I borrow the term "meaningful use" from the Health IT field to describe the need to maximize the utility of genomic data and make suggestions for how to address the current limitations.

Early Emergence Phase of SARS-CoV-2 Delta Variant in Florida, US.

SARS-CoV-2, the causative agent of COVID-19, emerged in late 2019. The highly contagious B.1.617.2 (Delta) variant of concern (VOC) was first identified in October 2020 in India and subsequently disseminated worldwide, later becoming the dominant lineage in the US. Understanding the local transmission dynamics of early SARS-CoV-2 introductions may inform actionable mitigation efforts during subsequent pandemic waves. Yet, despite considerable genomic analysis of SARS-CoV-2 in the US, several gaps remain. Here, we explore the early emergence of the Delta variant in Florida, US using phylogenetic analysis of representative Florida and globally sampled genomes. We find multiple independent introductions into Florida primarily from North America and Europe, with a minority originating from Asia. These introductions led to three distinct clades that demonstrated varying relative rates of transmission and possessed five distinct substitutions that were 3-21 times more prevalent in the Florida sample as compared to the global sample. Our results underscore the benefits of routine viral genomic surveillance to monitor epidemic spread and support the need for more comprehensive genomic epidemiology studies of emerging variants. In addition, we provide a model of epidemic spread of newly emerging VOCs that can inform future public health responses.

Genomic Epidemiology and Global Population Structure of Exfoliative Toxin A-Producing Staphylococcus aureus Strains Associated With Staphylococcal Scalded Skin Syndrome.

Staphylococci producing exfoliative toxins are the causative agents of staphylococcal scalded skin syndrome (SSSS). Exfoliative toxin A (ETA) is encoded by eta, which is harbored on a temperate bacteriophage ΦETA. A recent increase in the incidence of SSSS in North America has been observed; yet it is largely unknown whether this is the result of host range expansion of ΦETA or migration and emergence of established lineages. Here, we detail an outbreak investigation of SSSS in a neonatal intensive care unit, for which we applied whole-genome sequencing (WGS) and phylogenetic analysis of Staphylococcus aureus isolates collected from cases and screening of healthcare workers. We identified the causative strain as a methicillin-susceptible S. aureus (MSSA) sequence type 582 (ST582) possessing ΦETA. To then elucidate the global distribution of ΦETA among staphylococci, we used a recently developed tool to query extant bacterial WGS data for biosamples containing eta, which yielded 436 genomes collected between 1994 and 2019 from 32 countries. Applying population genomic analysis, we resolved the global distribution of S. aureus with lysogenized ΦETA and assessed antibiotic resistance determinants as well as the diversity of ΦETA. The population is highly structured with eight dominant sequence clusters (SCs) that generally aligned with S. aureus ST clonal complexes. The most prevalent STs included ST109 (24.3%), ST15 (13.1%), ST121 (10.1%), and ST582 (7.1%). Among strains with available data, there was an even distribution of isolates from carriage and disease. Only the SC containing ST121 had significantly more isolates collected from disease (69%, n = 46) than carriage (31%, n = 21). Further, we identified 10.6% (46/436) of strains as methicillin-resistant S. aureus (MRSA) based on the presence of mecA and the SCCmec element. Assessment of ΦETA diversity based on nucleotide identity revealed 27 phylogroups, and prophage gene content further resolved 62 clusters. ΦETA was relatively stable within lineages, yet prophage variation is geographically structured. This suggests that the reported increase in incidence is associated with migration and expansion of existing lineages, not the movement of ΦETA to new genomic backgrounds. This revised global view reveals that ΦETA is diverse and is widely distributed on multiple genomic backgrounds whose distribution varies geographically.

SARS-CoV-2 shifting transmission dynamics and hidden reservoirs potentially limit efficacy of public health interventions in Italy.

We investigated SARS-CoV-2 transmission dynamics in Italy, one of the countries hit hardest by the pandemic, using phylodynamic analysis of viral genetic and epidemiological data. We observed the co-circulation of multiple SARS-CoV-2 lineages over time, which were linked to multiple importations and characterized by large transmission clusters concomitant with a high number of infections. Subsequent implementation of a three-phase nationwide lockdown strategy greatly reduced infection numbers and hospitalizations. Yet we present evidence of sustained viral spread among sporadic clusters acting as "hidden reservoirs" during summer 2020. Mathematical modelling shows that increased mobility among residents eventually catalyzed the coalescence of such clusters, thus driving up the number of infections and initiating a new epidemic wave. Our results suggest that the efficacy of public health interventions is, ultimately, limited by the size and structure of epidemic reservoirs, which may warrant prioritization during vaccine deployment.

Complete Genome Sequence of Exfoliative Toxin-Producing Staphylococcus aureus Strain MSSA_SSSS_01, Obtained from a Case of Staphylococcal Scalded-Skin Syndrome.

Here, we announce the complete genome sequence of an exfoliative toxin-producing strain of Staphylococcus aureus sequence type 582 (ST582), isolated from a case of staphylococcal scalded-skin syndrome. The genome consists of a single circularized unitig with a total length of 2,792,190 bp carrying 2,699 genes. The genome is the basis for future epidemiological and genomic studies.

Frequency-dependent selection can forecast evolution in Streptococcus pneumoniae.

Predicting how pathogen populations will change over time is challenging. Such has been the case with Streptococcus pneumoniae, an important human pathogen, and the pneumococcal conjugate vaccines (PCVs), which target only a fraction of the strains in the population. Here, we use the frequencies of accessory genes to predict changes in the pneumococcal population after vaccination, hypothesizing that these frequencies reflect negative frequency-dependent selection (NFDS) on the gene products. We find that the standardized predicted fitness of a strain, estimated by an NFDS-based model at the time the vaccine is introduced, enables us to predict whether the strain increases or decreases in prevalence following vaccination. Further, we are able to forecast the equilibrium post-vaccine population composition and assess the invasion capacity of emerging lineages. Overall, we provide a method for predicting the impact of an intervention on pneumococcal populations with potential application to other bacterial pathogens in which NFDS is a driving force.

Toxigenic Vibrio cholerae evolution and establishment of reservoirs in aquatic ecosystems.

The spread of cholera in the midst of an epidemic is largely driven by direct transmission from person to person, although it is well-recognized that Vibrio cholerae is also capable of growth and long-term survival in aquatic ecosystems. While prior studies have shown that aquatic reservoirs are important in the persistence of the disease on the Indian subcontinent, an epidemiological view postulating that locally evolving environmental V. cholerae contributes to outbreaks outside Asia remains debated. The single-source introduction of toxigenic V. cholerae O1 in Haiti, one of the largest outbreaks occurring this century, with 812,586 suspected cases and 9,606 deaths reported through July 2018, provided a unique opportunity to evaluate the role of aquatic reservoirs and assess bacterial transmission dynamics across environmental boundaries. To this end, we investigated the phylogeography of both clinical and aquatic toxigenic V. cholerae O1 isolates and show robust evidence of the establishment of aquatic reservoirs as well as ongoing evolution of V. cholerae isolates from aquatic sites. Novel environmental lineages emerged from sequential population bottlenecks, carrying mutations potentially involved in adaptation to the aquatic ecosystem. Based on such empirical data, we developed a mixed-transmission dynamic model of V. cholerae, where aquatic reservoirs actively contribute to genetic diversification and epidemic emergence, which underscores the complexity of transmission pathways in epidemics and endemic settings and the need for long-term investments in cholera control at both human and environmental levels.