to analyze nvt penicillin - nonsusceptible pneumococcus ( pnsp ) detected in patients with invasive pneumococcal disease , we used data from the active bacterial core surveillance ( abcs ) system , a population- and laboratory - based collaborative system between the centers for disease control and prevention and state health departments and academic institutions in 10 states ( california , colorado , connecticut , georgia , maryland , minnesota , new mexico , new york , oregon , and tennessee )
. we considered pnsp non - pcv13 serotypes detected in patients in all age groups from 2009 ( pre - pcv13 , n = 285 patients ) through 2012 ( post - pcv13 , 339 patients ) .
nonsusceptibility was based on the meningitis breakpoint ( mic > 0.12 g / ml ) , as recommended by the clinical and laboratory standards institute ( 10 ) .
/ c were grouped together as 15bc because of the reported reversible switching between the 2 serotypes , which makes the precise differentiation of these serotypes difficult ( 11 ) . to determine whether geographic differences in the proportions of pnsp were consistent across serotypes , we calculated the proportions of pnsp for each of the 7 most common nvt serotypes ( 15a , 15bc , 16f , 23a , 23b , 33f , and 35b ) across the 10 sites for 2009 and 2012 .
we found that serotypes with the highest proportions of pnsp in 2012 already had high resistance in 2009 ( figure 1 ) .
we calculated the spearman correlation coefficient between the proportion of pnsp for each pair of serotypes across states in 2009 ( range 0.09 to 0.66 ) and 2012 ( range 0.300.79 ) ( technical appendix ) .
we found significant overall correlation in 2009 and 2012 ( p<0.001 for both years ) , indicating that sites with high proportions of pnsp in 1 serotype typically will also have high proportions of pnsp in other serotypes .
this finding suggests that differences in selection pressure account for the geographic variation in the proportions of pnsp .
comparison of proportion of nonvaccine type serotypes with penicillin resistance , by serotype , united states , 2009 and 2012 . based on active bacterial core surveillance system data from 10 us states .
we next implemented a standardized regression approach , used previously to analyze the pneumococcal - resistance patterns pre - pcv7 ( 7 ) and post - pcv7 ( 8) ( technical appendix ) . to investigate the source of geographic variation in the proportion of pnsp , we tested the hypotheses that either geographic heterogeneity in serotype distribution ( std1 ) , or serotype - specific differences in penicillin resistance ( std2 ) were responsible for the observed variation .
these effects were quantified by regressing crude versus standardized prevalence of penicillin resistance ( figure 2 ; technical appendix ) , by which a regression slope of 1 would indicate that the factor considered had zero effect . by using 2009 data , we found that regression slopes for std1 ( 0.445 , 95% ci 0.083 to 0.972 ) and std2 ( 0.463 , 95% ci 0.013 to 0.939 ) indicate that both factors played a similar intermediate role in generating this geographic variation in penicillin resistance , with neither 95% ci containing 1 . in 2012 ,
the regression coefficient for std2 was higher ( 0.634 , 95% ci 0.141.128 ) , whereas the coefficient for std1 decreased ( 0.367 , 95% ci 0.025 to 0.758 ) .
although these changes are not statistically significant relative to 2009 , they might suggest shifting contributions to the observed variation in proportions of pnsp after the introduction of the pcv13 vaccine in 2010 , with geographic differences in serotype distribution having an increased role and differences in serotype - specific pnsp becoming less important .
crude versus standardized proportions of nonvaccine type serotypes with penicillin nonsusceptibility , by state , united states , 2009 and 2012 , based on active bacterial core surveillance system data from 10 us states .
regression slopes with 95% cis are indicated in the upper left corner of each panel .
dashed lines represent the inverse - variance weighted ( red ) and unweighted ( gray ) regression slopes .
finally , we sought to quantify the rate of change in penicillin resistance during 20092012 in each state .
we documented the proportion of pnsp by state for the pre- and post - pcv13 periods ( technical appendix table 4 ) .
new mexico , maryland , and georgia saw the highest increases in the proportion of pnsp during 20092012 , whereas a slight decline was observed for colorado , new york , and connecticut .
although the distribution of serotypes might greatly fluctuate among geographic regions immediately after vaccine introduction , the overall proportions of nvt serotypes with penicillin resistance across the country might not vary significantly between the pre- and post - vaccine periods . of potential importance
are the small increases in the proportions of pnsp serotypes not included in either vaccine that were observed between the implementation of pcv7 and pcv13 ( 12 ) , which might lay the foundation for changes post - pcv13 .
the marked variation in the proportion of penicillin resistance among states highlights the potential of local selective pressures to favor certain serotypes and resistant strains within each serotype to increase in frequency as the population returns to equilibrium ( 13 ) .
previous studies have already shown significant regional differences in antibiotic use and vaccination coverage across the united states ( 14,15 ) .
a combination of these factors , which will likely vary between and within regions , would greatly affect proportions of resistance across the country .   in our study , we observed that the dynamics of penicillin resistance continue to shift in the wake of vaccine introduction .
our postvaccine observations were recorded shortly after the introduction of the vaccine ; additional observations would be valuable to determine the stability of the postvaccine dynamics and any potential importance of the temporal changes we observed to factors contributing to variation in resistance levels .
further long - term nationwide surveillance of serotype dynamics is required to assess the multiple ecologic factors that influence antibiotic resistance in the pneumococcus in the conjugate vaccine era .
methods and additional results for a study assessing penicillin resistance of nonvaccine type pneumococcus before and after pcv13 introduction , united states .
 introduction of 13-valent pneumococcal conjugate vaccine in the united states was not associated with a significant change in prevalence of penicillin resistance in nonvaccine type serotypes because of the variable success of highly resistant serotypes . 
 differences in regional serotype distribution and serotype - specific resistance contributed to geographic heterogeneity of penicillin resistance . 