A novel invasive Streptococcus pyogenes variant sublineage derived through recombinational replacement of the emm12 genomic region.

Group A streptococcal strains potentially acquire new M protein gene types through genetic recombination (emm switching). To detect such variants, we screened 12,596 invasive GAS genomes for strains of differing emm types that shared the same multilocus sequence type (ST). Through this screening we detected a variant consisting of 16 serum opacity factor (SOF)-positive, emm pattern E, emm82 isolates that were ST36, previously only associated with SOF-negative, emm pattern A, emm12. The 16 emm82/ST36 isolates were closely interrelated (pairwise SNP distance of 0-43), and shared the same emm82-containing recombinational fragment. emm82/ST36 isolates carried the sof12 structural gene, however the sof12 indel characteristic of emm12 strains was corrected to confer the SOF-positive phenotype. Five independent emm82/ST36 invasive case isolates comprised two sets of genetically indistinguishable strains. The emm82/ST36 isolates were primarily macrolide resistant (12/16 isolates), displayed at least 4 different core genomic arrangements, and carried 11 different combinations of virulence and resistance determinants. Phylogenetic analysis revealed that emm82/ST36 was within a minor (non-clade 1) portion of ST36 that featured almost all ST36 antibiotic resistance. This work documents emergence of a rapidly diversifying variant that is the first confirmed example of an emm pattern A strain switched to a pattern E strain.

Variation in pneumococcal invasiveness metrics is driven by serotype carriage duration and initial risk of disease.

Streptococcus pneumoniae is an opportunistic pathogen that, while usually carried asymptomatically, can cause severe invasive diseases like meningitis and bacteremic pneumonia. A central goal in S. pneumoniae public health management is to identify which serotypes (immunologically distinct strains) pose the most risk of invasive disease. The most common invasiveness metrics use cross-sectional data (i.e., invasive odds ratios (IOR)), or longitudinal data (i.e., attack rates (AR)). To assess the reliability of these metrics we developed an epidemiological model of carriage and invasive disease. Our mathematical analyses illustrate qualitative failures with the IOR metric (e.g., IOR can decline with increasing invasiveness parameters). Fitting the model to both longitudinal and cross-sectional data, our analysis supports previous work indicating that invasion risk is maximal at or near time of colonization. This pattern of early invasive disease risk leads to substantial (up to 5-fold) biases when estimating underlying differences in invasiveness from IOR metrics, due to the impact of carriage duration on IOR. Together, these results raise serious concerns with the IOR metric as a basis for public health decision-making and lend support for multiple alternate metrics including AR.

Expansion of Invasive Group A Streptococcus M1UK Lineage in Active Bacterial Core Surveillance, United States, 2019‒2021.

From 2015-2018 to 2019‒2021, hypertoxigenic M1UK lineage among invasive group A Streptococcus increased in the United States (1.7%, 21/1,230 to 11%, 65/603; p<0.001). M1UK was observed in 9 of 10 states, concentrated in Georgia (n = 41), Tennessee (n = 13), and New York (n = 13). Genomic cluster analysis indicated recent expansions.

Continued Increase of Erythromycin Nonsusceptibility and Clindamycin Nonsusceptibility Among Invasive Group A Streptococci Driven by Genomic Clusters, United States, 2018-2019.

We analyzed 9630 invasive GAS surveillance isolates in the USA. From 2015-2017 to 2018-2019, significant increases in erythromycin-nonsusceptibility (18% vs 25%) and clindamycin-nonsusceptibility (17% vs 24%) occurred, driven by rapid expansions of genomic subclones. Prevention and control of clustered infections appear key to containing antimicrobial resistance.

Benchmark datasets for SARS-CoV-2 surveillance bioinformatics.

Background: Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the cause of coronavirus disease 2019 (COVID-19), has spread globally and is being surveilled with an international genome sequencing effort. Surveillance consists of sample acquisition, library preparation, and whole genome sequencing. This has necessitated a classification scheme detailing Variants of Concern (VOC) and Variants of Interest (VOI), and the rapid expansion of bioinformatics tools for sequence analysis. These bioinformatic tools are means for major actionable results: maintaining quality assurance and checks, defining population structure, performing genomic epidemiology, and inferring lineage to allow reliable and actionable identification and classification. Additionally, the pandemic has required public health laboratories to reach high throughput proficiency in sequencing library preparation and downstream data analysis rapidly. However, both processes can be limited by a lack of a standardized sequence dataset. Methods: We identified six SARS-CoV-2 sequence datasets from recent publications, public databases and internal resources. In addition, we created a method to mine public databases to identify representative genomes for these datasets. Using this novel method, we identified several genomes as either VOI/VOC representatives or non-VOI/VOC representatives. To describe each dataset, we utilized a previously published datasets format, which describes accession information and whole dataset information. Additionally, a script from the same publication has been enhanced to download and verify all data from this study. Results: The benchmark datasets focus on the two most widely used sequencing platforms: long read sequencing data from the Oxford Nanopore Technologies platform and short read sequencing data from the Illumina platform. There are six datasets: three were derived from recent publications; two were derived from data mining public databases to answer common questions not covered by published datasets; one unique dataset representing common sequence failures was obtained by rigorously scrutinizing data that did not pass quality checks. The dataset summary table, data mining script and quality control (QC) values for all sequence data are publicly available on GitHub: https://github.com/CDCgov/datasets-sars-cov-2. Discussion: The datasets presented here were generated to help public health laboratories build sequencing and bioinformatics capacity, benchmark different workflows and pipelines, and calibrate QC thresholds to ensure sequencing quality. Together, improvements in these areas support accurate and timely outbreak investigation and surveillance, providing actionable data for pandemic management. Furthermore, these publicly available and standardized benchmark data will facilitate the development and adjudication of new pipelines.

Invasive Pneumococcal Disease Clusters Disproportionally Impact Persons Experiencing Homelessness, Injecting Drug Users, and the Western United States.

BACKGROUND: Invasive pneumococcal disease (IPD) isolates forming genomic clusters can reflect rapid disease transmission between vulnerable individuals. METHODS: We performed whole genome sequencing of 2820 IPD isolates recovered during 2019 through Centers for Disease Control and Prevention's Active Bacterial Core surveillance to provide strain information (serotypes, resistance, genotypes), and 2778 of these genomes were analyzed to detect highly related genomic clusters. RESULTS: Isolates from persons experiencing homelessness (PEH) were more often within genomic clusters than those from persons not experiencing homelessness (PNEH) (105/198 [53.0%] vs 592/2551 [23.2%]; P < .001). The 4 western sites accounted for 33.4% (929/2778) of isolates subjected to cluster analysis yet accounted for 48.7% (343/705) of clustering isolates (P < .001) and 75.8% (150/198) of isolates recovered from PEH (P < .001). Serotypes most frequent among PEH were (in rank order) 12F, 4, 3, 9N, 8, 20, and 22F, all of which were among the 10 serotypes exhibiting the highest proportions of clustering isolates among all cases. These serotypes accounted for 44.9% (1265/2820) of all IPD cases and are included within available vaccines. CONCLUSIONS: We identified serotype-specific and geographic differences in IPD transmission. We show the vulnerability of PEH within different regions to rapidly spreading IPD transmission networks representing several pneumococcal serotypes included in available vaccines.

Genomic Characterization of Group A Streptococci Causing Pharyngitis and Invasive Disease in Colorado, USA, June 2016- April 2017.

BACKGROUND: The genomic features and transmission link of circulating Group A Streptococcus (GAS) strains causing different disease types, such as pharyngitis and invasive disease, are not well understood. METHODS: We used whole-genome sequencing to characterize GAS isolates recovered from persons with pharyngitis and invasive disease in the Denver metropolitan area from June 2016 to April 2017. RESULTS: The GAS isolates were cultured from 236 invasive and 417 pharyngitis infections. Whole-genome sequencing identified 34 emm types. Compared with pharyngitis isolates, invasive isolates were more likely to carry the erm family genes (23% vs 7.4%, P<.001), which confer resistance to erythromycin and clindamycin (including inducible resistance), and covS gene inactivation (7% vs 0.5%, P<.001). Whole-genome sequencing identified 97 genomic clusters (433 isolates; 2-65 isolates per cluster) that consisted of genomically closely related isolates (median single-nucleotide polymorphism=3 [interquartile range, 1-4] within cluster). Thirty genomic clusters (200 isolates; 31% of all isolates) contained both pharyngitis and invasive isolates and were found in 11 emm types. CONCLUSIONS: In the Denver metropolitan population, mixed disease types were commonly seen in clusters of closely related isolates, indicative of overlapping transmission networks. Antibiotic-resistance and covS inactivation was disproportionally associated with invasive disease.

Low but Increasing Prevalence of Reduced Beta-lactam Susceptibility Among Invasive Group B Streptococcal Isolates, US Population-Based Surveillance, 1998-2018.

BACKGROUND: Invasive group B Streptococcus (iGBS) isolates with mutations in the pbp2x gene that encodes penicillin binding protein 2x can have reduced beta-lactam susceptibility (RBLS) when susceptible by Clinical and Laboratory Standards Institute (CLSI) criteria. We assessed the emergence and characteristics of RBLS strains in US iGBS isolates. METHODS: We analyzed iGBS isolates from 8 multistate population-based surveillance sites from 1998 to 2018. During 1998-2014, phenotypic antimicrobial susceptibility was determined by broth microdilution; criteria for 6 antibiotics were used to identify RBLS, followed by whole-genome sequencing (WGS). WGS for all isolates was added in 2015; we used phenotypic and genotypic results of >2000 isolates to validate phenotypic RBLS criteria and genotypic predictions. Since 2016, WGS has been used to screen for RBLS with broth microdilution confirmation of predicted RBLS isolates. RESULTS: Of 28 269 iGBS isolates, 28 (0.1%) were nonsusceptible by CLSI criteria; 137 (0.5%) met RBLS criteria. RBLS isolates were detected in all Active Bacterial Core surveillance sites. The RBLS proportion increased, especially since 2013 (odds ratio, 1.17; 95% CI, 1.03-1.32); the proportion that were nonsusceptible remained stable. CONCLUSIONS: The RBSL proportion was low but increasing among US iGBS isolates. Ongoing monitoring is needed to detect emerging threats to prevention and treatment of GBS infections.

Upsurge of Conjugate Vaccine Serotype 4 Invasive Pneumococcal Disease Clusters Among Adults Experiencing Homelessness in California, Colorado, and New Mexico.

After 7-valent pneumococcal conjugate vaccine introduction in the United States in 2000, invasive pneumococcal disease (IPD) due to serotype 4 greatly decreased in children and adults. Starting in 2013, serotype 4 IPD incidence increased among adults within 3 of 10 Active Bacterial Core surveillance sites. Of 325 serotype 4 cases among adults in 2010-2018, 36% were persons experiencing homelessness (PEH); incidence of serotype 4 IPD among PEH was 100-300 times higher than in the general population within these 3 areas. Genome sequencing for isolates recovered 2015-2018 (n = 246), revealed that increases in serotype 4 IPD were driven by lineages ST10172, ST244, and ST695. Within each lineage, clusters of near-identical isolates indicated close temporal relatedness. Increases in serotype 4 IPD were limited to Colorado, California, and New Mexico, with highest increases among PEH, who were at increased risk for exposure to and infections caused by these strains.

Invasive Pneumococcal Strain Distributions and Isolate Clusters Associated With Persons Experiencing Homelessness During 2018.

BACKGROUND: We aimed to characterize invasive pneumococcal disease (IPD) isolates collected from multistate surveillance in the United States during 2018 and examine within-serotype propensities of isolates to form related clusters. METHODS: We predicted strain features using whole genome sequencing obtained from 2885 IPD isolates obtained through the Center for Disease Control and Prevention's Active Bacterial Core surveillance (ABCs), which has a surveillance population of approximately 34.5 million individuals distributed among 10 states. Phylogenetic analysis was provided for serotypes accounting for ≥27 isolates. RESULTS: Thirteen-valent pneumococcal conjugate vaccine (PCV13) serotypes together with 6C accounted for 23 of 105 (21.9%) of isolates from children aged <5 years and 820 of 2780 (29.5%) isolates from those aged ≥5 years. The most common serotypes from adult IPD isolates were serotypes 3 (413/2780 [14.9%]), 22F (291/2780 [10.5%]), and 9N (191/2780 [6.9%]). Among child IPD isolates, serotypes 15BC (18/105 [17.1%]), 3 (11/105 [10.5%]), and 33F (10/105 [9.5%]) were most common. Serotypes 4, 12F, 20, and 7F had the highest proportions of isolates that formed related clusters together with the highest proportions of isolates from persons experiencing homelessness (PEH). Among 84 isolates from long-term care facilities, 2 instances of highly related isolate pairs from co-residents were identified. CONCLUSIONS: Non-PCV13 serotypes accounted for >70% of IPD in ABCs; however, PCV13 serotype 3 is the most common IPD serotype overall. Serotypes most common among PEH were more often associated with temporally related clusters identified both among PEH and among persons not reportedly experiencing homelessness.

Genomic Surveillance of Streptococcus pyogenes Strains Causing Invasive Disease, United States, 2016-2017.

BACKGROUND: Streptococcus pyogenes is a major cause of severe, invasive infections in humans. The bacterial pathogen harbors a wide array of virulence factors and exhibits high genomic diversity. Rapid changes of circulating strains in a community are common. Understanding the current prevalence and dynamics of S. pyogenes lineages could inform vaccine development and disease control strategies. METHODS: We used whole-genome sequencing (WGS) to characterize all invasive S. pyogenes isolates obtained through the United States Center for Disease Control and Prevention's Active Bacterial Core surveillance (ABCs) in 2016 and 2017. We determined the distribution of strain features, including emm type, antibiotic resistance determinants, and selected virulence factors. Changes in strain feature distribution between years 2016 and 2017 were evaluated. Phylogenetic analysis was used to identify expanding lineages within emm type. RESULTS: Seventy-one emm types were identified from 3873 isolates characterized. The emm types targeted by a 30-valent M protein-based vaccine accounted for 3230 (89%) isolates. The relative frequencies of emm types collected during the 2 years were similar. While all isolates were penicillin-susceptible, erythromycin-resistant isolates increased from 273 (16% of 2016 isolates) to 432 (23% of 2017 isolates), mainly driven by increase of the erm-positive emm types 92 and 83. The prevalence of 24 virulence factors, including 11 streptococcal pyrogenic toxins, ranged from 6 to 90%. In each of three emm types (emm 49, 82, and 92), a subgroup of isolates significantly expanded between 2016 and 2017 compared to isolates outside of the subgroup (P-values < 0.0001). Specific genomic sequence changes were associated with these expanded lineages. CONCLUSIONS: While the overall population structure of invasive S. pyogenes isolates in the United States remained stable, some lineages, including several that were antibiotic-resistant, increased between 2016 and 2017. Continued genomic surveillance can help monitor and characterize bacterial features associated with emerging strains from invasive infections.

Visualizing variation within Global Pneumococcal Sequence Clusters (GPSCs) and country population snapshots to contextualize pneumococcal isolates.

Knowledge of pneumococcal lineages, their geographic distribution and antibiotic resistance patterns, can give insights into global pneumococcal disease. We provide interactive bioinformatic outputs to explore such topics, aiming to increase dissemination of genomic insights to the wider community, without the need for specialist training. We prepared 12 country-specific phylogenetic snapshots, and international phylogenetic snapshots of 73 common Global Pneumococcal Sequence Clusters (GPSCs) previously defined using PopPUNK, and present them in Microreact. Gene presence and absence defined using Roary, and recombination profiles derived from Gubbins are presented in Phandango for each GPSC. Temporal phylogenetic signal was assessed for each GPSC using BactDating. We provide examples of how such resources can be used. In our example use of a country-specific phylogenetic snapshot we determined that serotype 14 was observed in nine unrelated genetic backgrounds in South Africa. The international phylogenetic snapshot of GPSC9, in which most serotype 14 isolates from South Africa were observed, highlights that there were three independent sub-clusters represented by South African serotype 14 isolates. We estimated from the GPSC9-dated tree that the sub-clusters were each established in South Africa during the 1980s. We show how recombination plots allowed the identification of a 20 kb recombination spanning the capsular polysaccharide locus within GPSC97. This was consistent with a switch from serotype 6A to 19A estimated to have occured in the 1990s from the GPSC97-dated tree. Plots of gene presence/absence of resistance genes (tet, erm, cat) across the GPSC23 phylogeny were consistent with acquisition of a composite transposon. We estimated from the GPSC23-dated tree that the acquisition occurred between 1953 and 1975. Finally, we demonstrate the assignment of GPSC31 to 17 externally generated pneumococcal serotype 1 assemblies from Utah via Pathogenwatch. Most of the Utah isolates clustered within GPSC31 in a USA-specific clade with the most recent common ancestor estimated between 1958 and 1981. The resources we have provided can be used to explore to data, test hypothesis and generate new hypotheses. The accessible assignment of GPSCs allows others to contextualize their own collections beyond the data presented here.

A mosaic tetracycline resistance gene tet(S/M) detected in an MDR pneumococcal CC230 lineage that underwent capsular switching in South Africa.

OBJECTIVES: We reported tet(S/M) in Streptococcus pneumoniae and investigated its temporal spread in relation to nationwide clinical interventions. METHODS: We whole-genome sequenced 12 254 pneumococcal isolates from 29 countries on an Illumina HiSeq sequencer. Serotype, multilocus ST and antibiotic resistance were inferred from genomes. An SNP tree was built using Gubbins. Temporal spread was reconstructed using a birth-death model. RESULTS: We identified tet(S/M) in 131 pneumococcal isolates and none carried other known tet genes. Tetracycline susceptibility testing results were available for 121 tet(S/M)-positive isolates and all were resistant. A majority (74%) of tet(S/M)-positive isolates were from South Africa and caused invasive diseases among young children (59% HIV positive, where HIV status was available). All but two tet(S/M)-positive isolates belonged to clonal complex (CC) 230. A global phylogeny of CC230 (n=389) revealed that tet(S/M)-positive isolates formed a sublineage predicted to exhibit resistance to penicillin, co-trimoxazole, erythromycin and tetracycline. The birth-death model detected an unrecognized outbreak of this sublineage in South Africa between 2000 and 2004 with expected secondary infections (effective reproductive number, R) of ∼2.5. R declined to ∼1.0 in 2005 and <1.0 in 2012. The declining epidemic could be related to improved access to ART in 2004 and introduction of pneumococcal conjugate vaccine (PCV) in 2009. Capsular switching from vaccine serotype 14 to non-vaccine serotype 23A was observed within the sublineage. CONCLUSIONS: The prevalence of tet(S/M) in pneumococci was low and its dissemination was due to an unrecognized outbreak of CC230 in South Africa. Capsular switching in this MDR sublineage highlighted its potential to continue to cause disease in the post-PCV13 era.

Emergent Invasive Group A Streptococcus dysgalactiae subsp. equisimilis, United States, 2015-2018.

The term group A Streptococcus is considered synonymous for the species Streptococcus pyogenes. We describe an emergent invasive S. dysgalactiae subspecies equisimilis lineage that obtained the group A antigen through a single ancestral recombination event between a group C S. dysgalactiae subsp. equisimilis strain and a group A S. pyogenes strain.

Genome-wide association analyses of invasive pneumococcal isolates identify a missense bacterial mutation associated with meningitis.

Bacterial mutations predisposing pneumococcus to causing meningitis, a more severe form of invasive pneumococcal disease (IPD), are largely unknown. Knowledge of such mutations may improve our understanding of pathogenesis and inform preventive strategies. Here we report a pneumococcal pbp1b gene mutation (pbp1bA641C causing N214T change in PBP1b transglycosylase domain) that is associated with meningitis in an exploratory cohort of IPD patients (n = 2054, p = 6.8 × 10-6), in an independent confirmatory cohort (n = 2518, p = 2.3 × 10-6), and in a combined analysis (n = 4572, p = 3.0 × 10-10). Patients infected by the pbp1b641C genotype pneumococci show 2.8-fold odds (95% CI 1.7 to 4.8) of meningitis compared to those infected by non-pbp1b641C pneumococci, after controlling for pneumococcal serotype, antibiotic resistance, and patient age. The pbp1bA641C change results in longer time needed for bacterial killing by antibiotic treatment and shows evidence of being under positive selection. Thus, a pneumococcal mutation conferring increased antibiotic tolerance is associated with meningitis among IPD patients.

A Population-Based Descriptive Atlas of Invasive Pneumococcal Strains Recovered Within the U.S. During 2015-2016.

Invasive pneumococcal disease (IPD) has greatly decreased since implementation in the U.S. of the 7 valent conjugate vaccine (PCV7) in 2000 and 13 valent conjugate vaccine (PCV13) in 2010. We used whole genome sequencing (WGS) to predict phenotypic traits (serotypes, antimicrobial phenotypes, and pilus determinants) and determine multilocus genotypes from 5334 isolates (~90% of cases) recovered during 2015-2016 through Active Bacterial Core surveillance. We identified 44 serotypes; 26 accounted for 98% of the isolates. PCV13 serotypes (inclusive of serotype 6C) accounted for 1503 (28.2%) isolates, with serotype 3 most common (657/5334, 12.3%), while serotypes 1 and 5 were undetected. Of 305 isolates from children <5 yrs, 60 (19.7%) were of PCV13 serotypes 19A, 19F, 3, 6B, and 23F (58/60 were 19A, 19F, or 3). We quantitated MLST-based lineages first detected during the post-PCV era (since 2002) that potentially arose through serotype-switching. The 7 predominant emergent post-PCV strain complexes included 23B/CC338, 15BC/CC3280, 19A/CC244, 4/CC439, 15A/CC156, 35B/CC156, and 15BC/CC156. These strains accounted for 332 isolates (6.2% of total) and were more frequently observed in children <5 yrs (17.7%; 54/305). Fifty-seven categories of recently emerged (in the post PCV7 period) putative serotype-switch variants were identified, accounting for 402 isolates. Many of these putative switch variants represented newly emerged resistant strains. Penicillin-nonsusceptibility (MICs > 0.12 µg/ml) was found among 22.4% (1193/5334) isolates, with higher penicillin MICs (2-8 µg/ml) found in 8.0% (425/5334) of isolates that were primarily (372/425, 87.5%) serotypes 35B and 19A. Most (792/1193, 66.4%) penicillin-nonsusceptible isolates were macrolide-resistant, 410 (34.4%) of which were erm gene positive and clindamycin-resistant. The proportion of macrolide-resistant isolates increased with increasing penicillin MICs; even isolates with reduced penicillin susceptibility (MIC = 0.06 µg/ml) were much more likely to be macrolide-resistant than basally penicillin-susceptible isolates (MIC < 0.03 µg/ml). The contribution of recombination to strain diversification was assessed through quantitating 35B/CC558-specific bioinformatic pipeline features among non-CC558 CCs and determining the sizes of gene replacements. Although IPD has decreased greatly and stabilized in the post-PCV13 era, the species continually generates recombinants that adapt to selective pressures exerted by vaccines and antimicrobials. These data serve as a baseline for monitoring future changes within each invasive serotype.

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.

Key features of invasive pneumococcal isolates recovered in Lima, Peru determined through whole genome sequencing.

Before PCV7 introduction, invasive pneumococcal disease (IPD) was responsible for approximately 12,000-18,000 deaths annually among children <5years in Latin America. In Peru, PCV7 was introduced in 2009. We used whole genome sequencing to deduce key features of invasive strains collected in Lima, Peru from 2006 to 2011. We sequenced 212 IPD isolates from 16 hospitals in Lima pre (2006-2009; n=133) and post (2010-2011; n=79) PCV7 introduction; 130 (61.3%) isolates were from children≤5years old. CDC's Streptococcus lab bioinformatics pipeline revealed serotypes, sequence types (STs), pilus genes, PBP types and other resistance determinants. During the pre-PCV7 period, serotype 14 was the most common serotype (24.8%), followed by 6B (20.3%), 19F (10.5%), and 23F (6.8%). Post-PCV7, the proportion of PCV7 serotype 6B decreased significantly (to 6.3%), while 19F (16.3%), 14 (15.0%), 23F (7.5%), and 19A (7.5%) were the most common serotypes; only serotypes 3 and 10A increased significantly. Overall, 82% (n=173) of all isolates carried at least one resistance determinant, including 72 (34%) isolates that carried resistance determinants against 3 or more antimicrobial classes; of these 72 isolates, 56 (78%) belonged to a PCV7 serotype. Eighty-two STs were identified, with 53 of them organized in 14 clonal complexes. ST frequencies were distributed differently pre and post-PCV7 introduction, with only 18 of the 57 STs identified in years 2006-2009 isolates also observed in years 2010-2011 isolates. The apparent expansion of a 19F/ST1421 lineage with predicted ß-lactam resistance (PBP type 13:16:20) and carrying resistance determinants against four additional antimicrobial classes was observed.

Validation of ß-lactam minimum inhibitory concentration predictions for pneumococcal isolates with newly encountered penicillin binding protein (PBP) sequences.

BACKGROUND: Genomic sequence-based deduction of antibiotic minimum inhibitory concentration (MIC) has great potential to enhance the speed and sensitivity of antimicrobial susceptibility testing. We previously developed a penicillin-binding protein (PBP) typing system and two methods (Random Forest (RF) and Mode MIC (MM)) that accurately predicted ß-lactam MICs for pneumococcal isolates carrying a characterized PBP sequence type (phenotypic ß-lactam MICs known for at least one isolate of this PBP type). This study evaluates the prediction performance for previously uncharacterized (new) PBP types and the probability of encountering new PBP types, both of which impact the overall prediction accuracy. RESULTS: The MM and RF methods were used to predict MICs of 4309 previously reported pneumococcal isolates in 2 datasets and the results were compared to the known broth microdilution MICs to 6 ß-lactams. Based on a method that specifically evaluated predictions for new PBP types, the RF results were more accurate than MM results for new PBP types and showed percent essential agreement (MICs agree within ±1 dilution) >97%, percent category agreement (interpretive results agree) >93%, major discrepancy (sensitive isolate predicted as resistant) rate < 1.2%, and very major discrepancy (resistant isolate predicted as sensitive) rate < 1.4% for all 6 ß-lactams. The identification of new PBP types over time was well approximated by a diminishingly increasing curve (Pearson's r = 0.99) and minimally impacted overall MIC prediction performance. CONCLUSIONS: MIC prediction using the RF method could be an accurate alternative of phenotypic susceptibility testing even in the presence of previously uncharacterized PBP types.

Invasive Serotype 35B Pneumococci Including an Expanding Serotype Switch Lineage, United States, 2015-2016.

We used whole-genome sequencing to characterize 199 nonvaccine serotype 35B pneumococcal strains that caused invasive pneumococcal disease (IPD) in the United States during 2015-2016 and related these findings to previous serotype 35B IPD data obtained by Active Bacterial Core surveillance. Penicillin-nonsusceptible 35B IPD increased during post-pneumococcal 7-valent conjugate vaccine years (2001-2009) and increased further after implementation of pneumococcal 13-valent conjugate vaccine in 2010. This increase was caused primarily by the 35B/sequence type (ST) 558 lineage. 35B/ST558 and vaccine serotype 9V/ST156 lineages were implicated as cps35B donor and recipient, respectively, for a single capsular switch event that generated emergent 35B/ST156 progeny in 6 states during 2015-2016. Three additional capsular switch 35B variants were identified, 2 of which also involved 35B/ST558 as cps35B donor. Spread of 35B/ST156 is of concern in view of past global predominance of pathogenic ST156 vaccine serotype strains. Protection against serotype 35B should be considered in next-generation pneumococcal vaccines.