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.

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.

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.

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.

Unexpected Predictors of Antibiotic Resistance in Housekeeping Genes of Staphylococcus Aureus

Methicillin-resistant Staphylococcus aureus (MRSA) is currently the most commonly identified antibiotic-resistant pathogen in US hospitals. Resistance to methicillin is carried by SCCmec genetic elements. Multilocus sequence typing (MLST) covers internal fragments of seven housekeeping genes of S. aureus. In conjunction with mec typing, MLST has been used to create an international nomenclature for S. aureus. MLST sequence types with a single nucleotide polymorphism (SNP) considered distinct. In this work, relationships among MLST SNPs and methicillin/oxacillin resistance or susceptibility were studied, using a public data base, by means of cross-tabulation tests, multivariable (phylogenetic) logistic regression (LR), decision trees, rule bases, and random forests (RF). Model performances were assessed through multiple cross-validation. Hierarchical clustering of SNPs was also employed to analyze mutational covariation. The number of instances with a known methicillin (oxacillin) antibiogram result was 1526 (649), where 63% (54%) was resistant to methicillin (oxacillin). In univariable analysis, several MLST SNPs were found strongly associated with antibiotic resistance/susceptibility. A RF model predicted correctly the resistance/susceptibility to methicillin and oxacillin in 75% and 63% of cases (cross-validated). Results were similar for LR. Hierarchical clustering of the aforementioned SNPs yielded a high level of covariation both within the same and different genes; this suggests strong genetic linkage between SNPs of housekeeping genes and antibiotic resistant associated genes. This finding provides a basis for rapid identification of antibiotic resistant S. arues lineages using a small number of genomic markers. The number of sites could subsequently be increased moderately to increase the sensitivity and specificity of genotypic tests for resistance that do not rely on the direct detection of the resistance marker itself.

Long-Term Intrahost Evolution of Methicillin Resistant Staphylococcus aureus Among Cystic Fibrosis Patients With Respiratory Carriage

Staphylococcus aureus is the most commonly identified airway colonizer of cystic fibrosis (CF) patients, and infections with methicillin-resistant S. aureus (MRSA) are associated with poor outcomes. Yet, little is known about the intrahost evolution of S. aureus among CF patients. We investigated convergent evolution and adaptation of MRSA among four CF patients with long-term respiratory carriage. For each patient, we performed whole-genome sequencing on an average of 21 isolates (range: 19–23) carried for a mean of 1,403 days (range: 903–1,679), including 25 pairs of isolates collected on the same day. We assessed intrahost diversity, population structure, evolutionary history, evidence of switched intergenic regions (IGRs), and signatures of adaptation in the context of patient age, antibiotic treatment, and co-colonizing microbes. Phylogenetic analysis delineated distinct multilocus sequence type ST5 (n = 3) and ST72 (n = 1) clonal populations in addition to sporadic, non-clonal isolates, and uncovered a putative transmission event. Variation in antibiotic resistance was observed within clonal populations, even among isolates collected on the same day. Rates of molecular evolution ranged from 2.21 to 8.64 nucleotide polymorphisms per year, and lineage ages were consistent with acquisition of colonization in early childhood followed by subsequent persistence of multiple sub-populations. Selection analysis of 1,622 core genes present in all four clonal populations (n = 79) found 11 genes variable in three subjects – most notably, ATP-dependent protease clpX, 2-oxoglutarate dehydrogenase odhA, fmtC, and transcription-repair coupling factor mfd. Only one gene, staphylococcal protein A (spa), was found to have evidence of gene-wide diversifying selection. We identified three instances of intrahost IGR switching events, two of which flanked genes related to quorum sensing. The complex microbial ecology of the CF airway poses challenges for management. We illustrate appreciable intrahost diversity as well as persistence of a dominant lineage. We also show that intrahost adaptation is a continual process, despite purifying selective pressure, and provide targets that should be investigated further for their function in CF adaptation.

Population genomics of pneumococcal carriage in Massachusetts children following introduction of PCV-13

The 13-valent pneumococcal conjugate vaccine (PCV-13) was introduced in the United States in 2010. Using a large paediatric carriage sample collected from shortly after the introduction of PCV-7 to several years after the introduction of PCV-13, we investigate alterations in the composition of the pneumococcal population following the introduction of PCV-13, evaluating the extent to which the post-vaccination non-vaccine type (NVT) population mirrors that from prior to vaccine introduction and the effect of PCV-13 on vaccine type lineages. Draft genome assemblies from 736 newly sequenced and 616 previously published pneumococcal carriage isolates from children in Massachusetts between 2001 and 2014 were analysed. Isolates were classified into one of 22 sequence clusters (SCs) on the basis of their core genome sequence. We calculated the SC diversity for each sampling period as the probability that any two randomly drawn isolates from that period belong to different SCs. The sampling period immediately after the introduction of PCV-13 (2011) was found to have higher diversity than preceding (2007) or subsequent (2014) sampling periods {Simpson’s D 2007: 0.915 [95 % confidence interval (CI) 0.901, 0.929]; 2011:  0.935 [0.927, 0.942]; 2014 :  0.912 [0.901, 0.923]}. Amongst NVT isolates, we found the distribution of SCs in 2011 to be significantly different from that in 2007 or 2014 (Fisher’s exact test P=0.018, 0.0078), but did not find a difference comparing 2007 to 2014 (Fisher’s exact test P=0.24), indicating greater similarity between samples separated by a longer time period than between samples from closer time periods. We also found changes in the accessory gene content of the NVT population between 2007 and 2011 to have been reduced by 2014. Amongst the new serotypes targeted by PCV-13, four were present in our sample. The proportion of our sample composed of PCV-13-only vaccine serotypes 19A, 6C and 7F decreased between 2007 and 2014, but no such reduction was seen for serotype 3. We did, however, observe differences in the genetic composition of the pre- and post-PCV-13 serotype 3 population. Our isolates were collected during discrete sampling periods from a small geographical area, which may limit the generalizability of our findings. Pneumococcal diversity increased immediately following the introduction of PCV-13, but subsequently returned to pre-vaccination levels. This is reflected in the distribution of NVT lineages, and, to a lesser extent, their accessory gene frequencies. As such, there may be a period during which the population is particularly disrupted by vaccination before returning to a more stable distribution. The persistence and shifting genetic composition of serotype 3 is a concern and warrants further investigation.

Global emergence and population dynamics of divergent serotype 3 CC180 pneumococci

Streptococcus pneumoniae serotype 3 remains a significant cause of morbidity and mortality worldwide, despite inclusion in the 13-valent pneumococcal conjugate vaccine (PCV13). Serotype 3 increased in carriage since the implementation of PCV13 in the USA, while invasive disease rates remain unchanged. We investigated the persistence of serotype 3 in carriage and disease, through genomic analyses of a global sample of 301 serotype 3 isolates of the Netherlands3-31 (PMEN31) clone CC180, combined with associated patient data and PCV utilization among countries of isolate collection. We assessed phenotypic variation between dominant clades in capsule charge (zeta potential), capsular polysaccharide shedding, and susceptibility to opsonophagocytic killing, which have previously been associated with carriage duration, invasiveness, and vaccine escape. We identified a recent shift in the CC180 population attributed to a lineage termed Clade II, which was estimated by Bayesian coalescent analysis to have first appeared in 1968 [95% HPD: 1939-1989] and increased in prevalence and effective population size thereafter. Clade II isolates are divergent from the pre-PCV13 serotype 3 population in non-capsular antigenic composition, competence, and antibiotic susceptibility, the last of which resulting from the acquisition of a Tn916-like conjugative transposon. Differences in recombination rates among clades correlated with variations in the ATP-binding subunit of Clp protease, as well as amino acid substitutions in the comCDE operon. Opsonophagocytic killing assays elucidated the low observed efficacy of PCV13 against serotype 3. Variation in PCV13 use among sampled countries was not independently correlated with the CC180 population shift; therefore, genotypic and phenotypic differences in protein antigens and, in particular, antibiotic resistance may have contributed to the increase of Clade II. Our analysis emphasizes the need for routine, representative sampling of isolates from disperse geographic regions, including historically under-sampled areas. We also highlight the value of genomics in resolving antigenic and epidemiological variations within a serotype, which may have implications for future vaccine development.

The impact of serotype-specific vaccination on phylodynamic parameters of Streptococcus pneumoniae and the pneumococcal pan-genome

In the United States, the introduction of the heptavalent pneumococcal conjugate vaccine (PCV) largely eliminated vaccine serotypes (VT); non-vaccine serotypes (NVT) subsequently increased in carriage and disease. Vaccination also disrupts the composition of the pneumococcal pangenome, which includes mobile genetic elements and polymorphic non-capsular antigens important for virulence, transmission, and pneumococcal ecology. Antigenic proteins are of interest for future vaccines; yet, little is known about how the they are affected by PCV use. To investigate the evolutionary impact of vaccination, we assessed recombination, evolution, and pathogen demographic history of 937 pneumococci collected from 1998-2012 among Navajo and White Mountain Apache Native American communities. We analyzed changes in the pneumococcal pangenome, focusing on metabolic loci and 19 polymorphic protein antigens. We found the impact of PCV on the pneumococcal population could be observed in reduced diversity, a smaller pangenome, and changing frequencies of accessory clusters of orthologous groups (COGs). Post-PCV7, diversity rebounded through clonal expansion of NVT lineages and inferred in-migration of two previously unobserved lineages. Accessory COGs frequencies trended toward pre-PCV7 values with increasing time since vaccine introduction. Contemporary frequencies of protein antigen variants are better predicted by pre-PCV7 values (1998-2000) than the preceding period (2006-2008), suggesting balancing selection may have acted in maintaining variant frequencies in this population. Overall, we present the largest genomic analysis of pneumococcal carriage in the United States to date, which includes a snapshot of a true vaccine-naïve community prior to the introduction of PCV7. These data improve our understanding of pneumococcal evolution and emphasize the need to consider pangenome composition when inferring the impact of vaccination and developing future protein-based pneumococcal vaccines.

Association of Pneumococcal Protein Antigen Serology With Age and Antigenic Profile of Colonizing Isolates

Background

Several Streptococcus pneumoniae proteins play a role in pathogenesis and are being investigated as vaccine targets. It is largely unknown whether naturally acquired antibodies reduce the risk of colonization with strains expressing a particular antigenic variant.

Methods

Serum immunoglobulin G (IgG) titers to 28 pneumococcal protein antigens were measured among 242 individuals aged <6 months-78 years in Native American communities between 2007 and 2009. Nasopharyngeal swabs were collected >- 30 days after serum collection, and the antigen variant in each pneumococcal isolate was determined using genomic data. We assessed the association between preexisting variant-specific antibody titers and subsequent carriage of pneumococcus expressing a particular antigen variant.

Results

Antibody titers often increased across pediatric groups before decreasing among adults. Individuals with low titers against group 3 pneumococcal surface protein C (PspC) variants were more likely to be colonized with pneumococci expressing those variants. For other antigens, variant-specific IgG titers do not predict colonization.

Conclusion

We observed an inverse association between variant-specific antibody concentration and homologous pneumococcal colonization for only 1 protein. Further assessment of antibody repertoires may elucidate the nature of antipneumococcal antibody-mediated mucosal immunity while informing vaccine development.

Non-toxigenic environmental Vibrio cholerae O1 strain from Haiti provides evidence of pre-pandemic cholera in Hispaniola

Vibrio cholerae is ubiquitous in aquatic environments, with environmental toxigenic V. cholerae O1 strains serving as a source for recurrent cholera epidemics and pandemic disease. However, a number of questions remain about long-term survival and evolution of V. cholerae strains within these aquatic environmental reservoirs. Through monitoring of the Haitian aquatic environment following the 2010 cholera epidemic, we isolated two novel non-toxigenic (ctxA/B-negative) Vibrio cholerae O1. These two isolates underwent whole-genome sequencing and were investigated through comparative genomics and Bayesian coalescent analysis. These isolates cluster in the evolutionary tree with strains responsible for clinical cholera, possessing genomic components of 6^th and 7^th pandemic lineages, and diverge from "modern" cholera strains around 1548 C.E. [95% HPD: 1532-1555]. Vibrio Pathogenicity Island (VPI)-1 was present; however, SXT/R391-family ICE and VPI-2 were absent. Rugose phenotype conversion and vibriophage resistance evidenced adaption for persistence in aquatic environments. The identification of V. cholerae O1 strains in the Haitian environment, which predate the first reported cholera pandemic in 1817, broadens our understanding of the history of pandemics. It also raises the possibility that these and similar environmental strains could acquire virulence genes from the 2010 Haitian epidemic clone, including the cholera toxin producing CTXϕ.

Genomic Epidemiology of Methicillin-Resistant Staphylococcus aureus in a Neonatal Intensive Care Unit

Despite infection prevention efforts, neonatal intensive care unit (NICU) patients remain at risk of Methicillin-resistant Staphylococcus aureus (MRSA) infection. Modes of transmission for healthcare-associated (HA) and community-associated (CA) MRSA remain poorly understood and may vary by genotype, hindering the development of effective prevention and control strategies. From 2008–2010, all patients admitted to a level III NICU were screened for MRSA colonization, and all available isolates were spa-typed. Spa-type t008, the most prevalent CA- genotype in the United States, spa-type t045, a HA- related genotype, and a convenience sample of strains isolated from 2003–2011, underwent whole-genome sequencing and phylodynamic analysis. Patient risk factors were compared between colonized and noncolonized infants, and virulence and resistance genes compared between spa-type t008 and non-t008 strains. Epidemiological and genomic data were used to estimate MRSA importations and acquisitions through transmission reconstruction. MRSA colonization was identified in 9.1% (177/1940) of hospitalized infants and associated with low gestational age and birth weight. Among colonized infants, low gestational age was more common among those colonized with t008 strains. Our data suggest that approximately 70% of colonizations were the result of transmission events within the NICU, with the remainder likely to reflect importations of “outside” strains. While risk of transmission within the NICU was not affected by spa-type, patterns of acquisition and importation differed between t008 and t045 strains. Phylodynamic analysis showed the effective population size of spa-type t008 has been exponentially increasing in both community and hospital, with spa-type t008 strains possessed virulence genes not found among t045 strains; t045 strains, in contrast, appeared to be of more recent origin, with a possible hospital source. Our data highlight the importance of both intra-NICU transmission and recurrent introductions in maintenance of MRSA colonization within the NICU environment, as well as spa-type-specific differences in epidemiology.

Intrahost Evolution of Methicillin-Resistant Staphylococcus aureus USA300 Among Individuals With Reoccurring Skin and Soft-Tissue Infections

BACKGROUND

Methicillin-resistant Staphylococcus aureus (MRSA) USA300 is the leading cause of MRSA infections in the United States and has caused an epidemic of skin and soft-tissue infections. Recurrent infections with USA300 MRSA are common, yet intrahost evolution during persistence on an individual has not been studied. This gap hinders the ability to clinically manage recurrent infections and reconstruct transmission networks.

METHODS

To characterize bacterial intrahost evolution, we examined the clinical courses of 4 subjects with 3-6 recurrent USA300 MRSA infections, using patient clinical data, including antibiotic exposure history, and whole-genome sequencing and phylogenetic analysis of all available MRSA isolates (n = 29).

RESULTS

Among sequential isolates, we found variability in diversity, accumulation of mutations, and mobile genetic elements. Selection for antimicrobial-resistant populations was observed through both an increase in the number of plasmids conferring multidrug resistance and strain replacement by a resistant population. Two of 4 subjects had strain replacement with a genetically distinct USA300 MRSA population.

DISCUSSIONS

During a 5-year period in 4 subjects, we identified development of antimicrobial resistance, intrahost evolution, and strain replacement among isolates from patients with recurrent MRSA infections. This calls into question the efficacy of decolonization to prevent recurrent infections and highlights the adaptive potential of USA300 and the need for effective sampling.

Impact of spatial dispersion, evolution, and selection on Ebola Zaire Virus epidemic waves

Ebola virus Zaire (EBOV) has reemerged in Africa, emphasizing the global importance of this pathogen. Amidst the response to the current epidemic, several gaps in our knowledge of EBOV evolution are evident. Specifically, uncertainty has been raised regarding the potential emergence of more virulent viral variants through amino acid substitutions. Glycoprotein (GP), an essential component of the EBOV genome, is highly variable and a potential site for the occurrence of advantageous mutations. For this study, we reconstructed the evolutionary history of EBOV by analyzing 65 GP sequences from humans and great apes over diverse locations across epidemic waves between 1976 and 2014. We show that, although patterns of spatial dispersion throughout Africa varied, the evolution of the virus has largely been characterized by neutral genetic drift. Therefore, the radical emergence of more transmissible variants is unlikely, a positive finding, which is increasingly important on the verge of vaccine deployment.

Molecular epidemiology of community-associated methicillin-resistant Staphylococcus aureus in the genomic era: a cross-sectional study

Methicillin-resistant Staphylococcus aureus (MRSA) is a leading cause of healthcare-associated infections and significant contributor to healthcare cost. Community-associated-MRSA (CA-MRSA) strains have now invaded healthcare settings. A convenience sample of 97 clinical MRSA isolates was obtained from seven hospitals during a one-week period in 2010. We employed a framework integrating Staphylococcus protein A typing and full-genome next-generation sequencing. Single nucleotide polymorphisms were analyzed using phylodynamics. Twenty-six t002, 48 t008, and 23 other strains were identified. Phylodynamic analysis of 30 t008 strains showed ongoing exponential growth of the effective population size the basic reproductive number (R0) ranging from 1.24 to 1.34. No evidence of hospital clusters was identified. The lack of phylogeographic clustering suggests that community introduction is a major contributor to emergence of CA-MRSA strains within hospitals. Phylodynamic analysis provides a powerful framework to investigate MRSA transmission between the community and hospitals, an understanding of which is essential for control.

HIV-1 Subtype distribution in morocco based on national sentinel surveillance data 2004-2005

BACKGROUND

Little is known about HIV-1 subtype distribution in Morocco. Some data suggest an emergence of new HIV subtypes. We conducted phylogenetic analysis on a nationally representative sample of 60 HIV-1 viral specimens collected during 2004-2005 through the Morocco national HIV sentinel surveillance survey.

RESULTS

While subtype B is still the most prevalent, 23.3% of samples represented non-B subtypes, the majority of which were classified as CRF02_AG (15%). Molecular clock analysis confirmed that the initial introduction of HIV-1B in Morocco probably came from Europe in the early 1980s. In contrast, the CRF02_AG strain appeared to be introduced from sub-Saharan Africa in two separate events in the 1990s.

CONCLUSIONS

Subtype CRF02_AG has been emerging in Morocco since the 1990s. More information about the factors introducing HIV subtype-specific transmission will inform the prevention strategy in the region.