Geographical migration and fitness dynamics of Streptococcus pneumoniae.

Streptococcus pneumoniae is a leading cause of pneumonia and meningitis worldwide. Many different serotypes co-circulate endemically in any one location1,2. The extent and mechanisms of spread and vaccine-driven changes in fitness and antimicrobial resistance remain largely unquantified. Here using geolocated genome sequences from South Africa (n = 6,910, collected from 2000 to 2014), we developed models to reconstruct spread, pairing detailed human mobility data and genomic data. Separately, we estimated the population-level changes in fitness of strains that are included (vaccine type (VT)) and not included (non-vaccine type (NVT)) in pneumococcal conjugate vaccines, first implemented in South Africa in 2009. Differences in strain fitness between those that are and are not resistant to penicillin were also evaluated. We found that pneumococci only become homogenously mixed across South Africa after 50 years of transmission, with the slow spread driven by the focal nature of human mobility. Furthermore, in the years following vaccine implementation, the relative fitness of NVT compared with VT strains increased (relative risk of 1.68; 95% confidence interval of 1.59-1.77), with an increasing proportion of these NVT strains becoming resistant to penicillin. Our findings point to highly entrenched, slow transmission and indicate that initial vaccine-linked decreases in antimicrobial resistance may be transient.

Pangenome-spanning epistasis and coselection analysis via de Bruijn graphs.

Studies of bacterial adaptation and evolution are hampered by the difficulty of measuring traits such as virulence, drug resistance, and transmissibility in large populations. In contrast, it is now feasible to obtain high-quality complete assemblies of many bacterial genomes thanks to scalable high-accuracy long-read sequencing technologies. To exploit this opportunity, we introduce a phenotype- and alignment-free method for discovering coselected and epistatically interacting genomic variation from genome assemblies covering both core and accessory parts of genomes. Our approach uses a compact colored de Bruijn graph to approximate the intragenome distances between pairs of loci for a collection of bacterial genomes to account for the impacts of linkage disequilibrium (LD). We demonstrate the versatility of our approach to efficiently identify associations between loci linked with drug resistance and adaptation to the hospital niche in the major human bacterial pathogens Streptococcus pneumoniae and Enterococcus faecalis.

Robust analysis of prokaryotic pangenome gene gain and loss rates with Panstripe.

Horizontal gene transfer (HGT) plays a critical role in the evolution and diversification of many microbial species. The resulting dynamics of gene gain and loss can have important implications for the development of antibiotic resistance and the design of vaccine and drug interventions. Methods for the analysis of gene presence/absence patterns typically do not account for errors introduced in the automated annotation and clustering of gene sequences. In particular, methods adapted from ecological studies, including the pangenome gene accumulation curve, can be misleading as they may reflect the underlying diversity in the temporal sampling of genomes rather than a difference in the dynamics of HGT. Here, we introduce Panstripe, a method based on generalized linear regression that is robust to population structure, sampling bias, and errors in the predicted presence/absence of genes. We show using simulations that Panstripe can effectively identify differences in the rate and number of genes involved in HGT events, and illustrate its capability by analyzing several diverse bacterial genome data sets representing major human pathogens.

Genomic and functional determinants of host spectrum in Group B Streptococcus.

Group B Streptococcus (GBS) is a major human and animal pathogen that threatens public health and food security. Spill-over and spill-back between host species is possible due to adaptation and amplification of GBS in new niches but the evolutionary and functional mechanisms underpinning those phenomena are poorly known. Based on analysis of 1,254 curated genomes from all major GBS host species and six continents, we found that the global GBS population comprises host-generalist, host-adapted and host-restricted sublineages, which are found across host groups, preferentially within one host group, or exclusively within one host group, respectively, and show distinct levels of recombination. Strikingly, the association of GBS genomes with the three major host groups (humans, cattle, fish) is driven by a single accessory gene cluster per host, regardless of sublineage or the breadth of host spectrum. Moreover, those gene clusters are shared with other streptococcal species occupying the same niche and are functionally relevant for host tropism. Our findings demonstrate (1) the heterogeneity of genome plasticity within a bacterial species of public health importance, enabling the identification of high-risk clones; (2) the contribution of inter-species gene transmission to the evolution of GBS; and (3) the importance of considering the role of animal hosts, and the accessory gene pool associated with their microbiota, in the evolution of multi-host bacterial pathogens. Collectively, these phenomena may explain the adaptation and clonal expansion of GBS in animal reservoirs and the risk of spill-over and spill-back between animals and humans.

The emergence and diversification of a zoonotic pathogen from within the microbiota of intensively farmed pigs.

The expansion and intensification of livestock production is predicted to promote the emergence of pathogens. As pathogens sometimes jump between species, this can affect the health of humans as well as livestock. Here, we investigate how livestock microbiota can act as a source of these emerging pathogens through analysis of Streptococcus suis, a ubiquitous component of the respiratory microbiota of pigs that is also a major cause of disease on pig farms and an important zoonotic pathogen. Combining molecular dating, phylogeography, and comparative genomic analyses of a large collection of isolates, we find that several pathogenic lineages of S. suis emerged in the 19th and 20th centuries, during an early period of growth in pig farming. These lineages have since spread between countries and continents, mirroring trade in live pigs. They are distinguished by the presence of three genomic islands with putative roles in metabolism and cell adhesion, and an ongoing reduction in genome size, which may reflect their recent shift to a more pathogenic ecology. Reconstructions of the evolutionary histories of these islands reveal constraints on pathogen emergence that could inform control strategies, with pathogenic lineages consistently emerging from one subpopulation of S. suis and acquiring genes through horizontal transfer from other pathogenic lineages. These results shed light on the capacity of the microbiota to rapidly evolve to exploit changes in their host population and suggest that the impact of changes in farming on the pathogenicity and zoonotic potential of S. suis is yet to be fully realized.

Diverse regulatory pathways modulate bet hedging of competence induction in epigenetically-differentiated phase variants of Streptococcus pneumoniae.

Despite enabling Streptococcus pneumoniae to acquire antibiotic resistance and evade vaccine-induced immunity, transformation occurs at variable rates across pneumococci. Phase variants of isolate RMV7, distinguished by altered methylation patterns driven by the translocating variable restriction-modification (tvr) locus, differed significantly in their transformation efficiencies and biofilm thicknesses. These differences were replicated when the corresponding tvr alleles were introduced into an RMV7 derivative lacking the locus. RNA-seq identified differential expression of the type 1 pilus, causing the variation in biofilm formation, and inhibition of competence induction in the less transformable variant, RMV7domi. This was partly attributable to RMV7domi's lower expression of ManLMN, which promoted competence induction through importing N-acetylglucosamine. This effect was potentiated by analogues of some proteobacterial competence regulatory machinery. Additionally, one of RMV7domi's phage-related chromosomal island was relatively active, which inhibited transformation by increasing expression of the stress response proteins ClpP and HrcA. However, HrcA increased competence induction in the other variant, with its effects depending on Ca2+ supplementation and heat shock. Hence the heterogeneity in transformation efficiency likely reflects the diverse signalling pathways by which it is affected. This regulatory complexity will modulate population-wide responses to synchronising quorum sensing signals to produce co-ordinated yet stochastic bet hedging behaviour.

Novel pneumococcal capsule type 33E results from the inactivation of glycosyltransferase WciE in vaccine type 33F.

The polysaccharide (PS) capsule is essential for immune evasion and virulence of Streptococcus pneumoniae. Existing pneumococcal vaccines are designed to elicit anticapsule antibodies; however, the effectiveness of these vaccines is being challenged by the emergence of new capsule types or variants. Herein, we characterize a newly discovered capsule type, 33E, that appears to have repeatedly emerged from vaccine type 33F via an inactivation mutation in the capsule glycosyltransferase gene, wciE. Structural analysis demonstrated that 33E and 33F share an identical repeat unit backbone [→5)-ß-D-Galf2Ac-(1→3)-ß-D-Galp-(1→3)-α-D-Galp-(1→3)-ß-D-Galf-(1→3)-ß-D-Glcp-(1→], except that a galactose (α-D-Galp) branch is present in 33F but not in 33E. Though the two capsule types were indistinguishable using conventional typing methods, the monoclonal antibody Hyp33FM1 selectively bound 33F but not 33E pneumococci. Further, we confirmed that wciE encodes a glycosyltransferase that catalyzes the addition of the branching α-D-Galp and that its inactivation in 33F strains results in the expression of the 33E capsule type. Though 33F and 33E share a structural and antigenic similarity, our pilot study suggested that immunization with a 23-valent pneumococcal PS vaccine containing 33F PS did not significantly elicit cross-opsonic antibodies to 33E. New conjugate vaccines that target capsule type 33F may not necessarily protect against 33E. Therefore, studies of new conjugate vaccines require knowledge of the newly identified capsule type 33E and reliable pneumococcal typing methods capable of distinguishing it from 33F.

The effect of recombination on the evolution of a population of Neisseria meningitidis.

Neisseria meningitidis (the meningococcus) is a major human pathogen with a history of high invasive disease burden, particularly in sub-Saharan Africa. Our current understanding of the evolution of meningococcal genomes is limited by the rarity of large-scale genomic population studies and lack of in-depth investigation of the genomic events associated with routine pathogen transmission. Here, we fill this knowledge gap by a detailed analysis of 2839 meningococcal genomes obtained through a carriage study of over 50,000 samples collected systematically in Burkina Faso, West Africa, before, during, and after the serogroup A vaccine rollout, 2009-2012. Our findings indicate that the meningococcal genome is highly dynamic, with highly recombinant loci and frequent gene sharing across deeply separated lineages in a structured population. Furthermore, our findings illustrate how population structure can correlate with genome flexibility, as some lineages in Burkina Faso are orders of magnitude more recombinant than others. We also examine the effect of selection on the population, in particular how it is correlated with recombination. We find that recombination principally acts to prevent the accumulation of deleterious mutations, although we do also find an example of recombination acting to speed the adaptation of a gene. In general, we show the importance of recombination in the evolution of a geographically expansive population with deep population structure in a short timescale. This has important consequences for our ability to both foresee the outcomes of vaccination programs and, using surveillance data, predict when lineages of the meningococcus are likely to become a public health concern.

Novel Multilocus Sequence Typing and Global Sequence Clustering Schemes for Characterizing the Population Diversity of Streptococcus mitis.

Streptococcus mitis is a common oral commensal and an opportunistic pathogen that causes bacteremia and infective endocarditis; however, the species has received little attention compared to other pathogenic streptococcal species. Effective and easy-to-use molecular typing tools are essential for understanding bacterial population diversity and biology, but schemes specific for S. mitis are not currently available. We therefore developed a multilocus sequence typing (MLST) scheme and defined sequence clusters or lineages of S. mitis using a comprehensive global data set of 322 genomes (148 publicly available and 174 newly sequenced). We used internal 450-bp sequence fragments of seven housekeeping genes (accA, gki, hom, oppC, patB, rlmN, and tsf) to define the MLST scheme and derived the global S. mitis sequence clusters using the PopPUNK clustering algorithm. We identified an initial set of 259 sequence types (STs) and 258 global sequence clusters. The schemes showed high concordance (100%), capturing extensive S. mitis diversity with strains assigned to multiple unique STs and global sequence clusters. The tools also identified extensive within- and between-host S. mitis genetic diversity among isolates sampled from a cohort of healthy individuals, together with potential transmission events, supported by both phylogeny and pairwise single nucleotide polymorphism (SNP) distances. Our novel molecular typing and strain clustering schemes for S. mitis allow for the integration of new strain data, are electronically portable at the PubMLST database (https://pubmlst.org/smitis), and offer a standardized approach to understanding the population structure of S. mitis. These robust tools will enable new insights into the epidemiology of S. mitis colonization, disease and transmission.

Discovery and Characterization of Pneumococcal Serogroup 36 Capsule Subtypes, Serotypes 36A and 36B.

Streptococcus pneumoniae can produce a wide breadth of antigenically diverse capsule types, a fact that poses a looming threat to the success of vaccines that target pneumococcal polysaccharide (PS) capsule. Yet, many pneumococcal capsule types remain undiscovered and/or uncharacterized. Prior sequence analysis of pneumococcal capsule synthesis (cps) loci suggested the existence of capsule subtypes among isolates identified as "serotype 36" according to conventional capsule typing methods. We discovered these subtypes represent two antigenically similar but distinguishable pneumococcal capsule serotypes, 36A and 36B. Biochemical analysis of their capsule PS structure reveals that both have the shared repeat unit backbone [→5)-α-d-Galf-(1→1)-d-Rib-ol-(5→P→6)-ß-d-ManpNAc-(1→4)-ß-d-Glcp-(1→] with two branching structures. Both serotypes have a ß-d-Galp branch to Ribitol. Serotypes 36A and 36B differ by the presence of a α-d-Glcp-(1→3)-ß-d-ManpNAc or α-d-Galp-(1→3)-ß-d-ManpNAc branch, respectively. Comparison of the phylogenetically distant serogroup 9 and 36 cps loci, which all encode this distinguishing glycosidic bond, revealed that the incorporation of Glcp (in types 9N and 36A) versus Galp (in types 9A, 9V, 9L, and 36B) is associated with the identity of four amino acids in the cps-encoded glycosyltransferase WcjA. Identifying functional determinants of cps-encoded enzymes and their impact on capsule PS structure is key to improving the resolution and reliability of sequencing-based capsule typing methods and discovering novel capsule variants indistinguishable by conventional serotyping methods.

RCandy: an R package for visualizing homologous recombinations in bacterial genomes.

SUMMARY: Homologous recombination is an important evolutionary process in bacteria and other prokaryotes, which increases genomic sequence diversity and can facilitate adaptation. Several methods and tools have been developed to detect genomic regions recently affected by recombination. Exploration and visualization of such recombination events can reveal valuable biological insights, but it remains challenging. Here, we present RCandy, a platform-independent R package for rapid, simple and flexible visualization of recombination events in bacterial genomes. AVAILABILITY AND IMPLEMENTATION: RCandy is an R package freely available for use under the MIT license. It is platform-independent and has been tested on Windows, Linux and MacOSX. The source code comes together with a detailed vignette available on GitHub at https://github.com/ChrispinChaguza/RCandy. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Pneumococcal Vaccines: Host Interactions, Population Dynamics, and Design Principles.

Streptococcus pneumoniae (the pneumococcus) is a nasopharyngeal commensal and respiratory pathogen. Most isolates express a capsule, the species-wide diversity of which has been immunologically classified into ∼100 serotypes. Capsule polysaccharides have been combined into multivalent vaccines widely used in adults, but the T cell independence of the antibody response means they are not protective in infants. Polysaccharide conjugate vaccines (PCVs) trigger a T cell-dependent response through attaching a carrier protein to capsular polysaccharides. The immune response stimulated by PCVs in infants inhibits carriage of vaccine serotypes (VTs), resulting in population-wide herd immunity. These were replaced in carriage by non-VTs. Nevertheless, PCVs drove reductions in infant pneumococcal disease, due to the lower mean invasiveness of the postvaccination bacterial population; age-varying serotype invasiveness resulted in a smaller reduction in adult disease. Alternative vaccines being tested in trials are designed to provide species-wide protection through stimulating innate and cellular immune responses, alongside antibodies to conserved antigens.

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.

Pneumococcal Colonization and Virulence Factors Identified Via Experimental Evolution in Infection Models.

Streptococcus pneumoniae is a commensal of the human nasopharynx and a major cause of respiratory and invasive disease. We examined adaptation and evolution of pneumococcus, within nasopharynx and lungs, in an experimental system where the selective pressures associated with transmission were removed. This was achieved by serial passage of pneumococci, separately, in mouse models of nasopharyngeal carriage or pneumonia. Passaged pneumococci became more effective colonizers of the respiratory tract and we observed several examples of potential parallel evolution. The cell wall-modifying glycosyltransferase LafA was under strong selection during lung passage, whereas the surface expressed pneumococcal vaccine antigen gene pvaA and the glycerol-3-phosphate dehydrogenase gene gpsA were frequent targets of mutation in nasopharynx-passaged pneumococci. These mutations were not identified in pneumococci that were separately evolved by serial passage on laboratory agar. We focused on gpsA, in which the same single nucleotide polymorphism arose in two independently evolved nasopharynx-passaged lineages. We describe a new role for this gene in nasopharyngeal carriage and show that the identified single nucleotide change confers resistance to oxidative stress and enhanced nasopharyngeal colonization potential. We demonstrate that polymorphisms in gpsA arise and are retained during human colonization. These findings highlight how within-host environmental conditions can determine trajectories of bacterial evolution. Relative invasiveness or attack rate of pneumococcal lineages may be defined by genes that make niche-specific contributions to bacterial fitness. Experimental evolution in animal infection models is a powerful tool to investigate the relative roles played by pathogen virulence and colonization factors within different host niches.

Fast and flexible bacterial genomic epidemiology with PopPUNK.

The routine use of genomics for disease surveillance provides the opportunity for high-resolution bacterial epidemiology. Current whole-genome clustering and multilocus typing approaches do not fully exploit core and accessory genomic variation, and they cannot both automatically identify, and subsequently expand, clusters of significantly similar isolates in large data sets spanning entire species. Here, we describe PopPUNK (Population Partitioning Using Nucleotide K -mers), a software implementing scalable and expandable annotation- and alignment-free methods for population analysis and clustering. Variable-length k-mer comparisons are used to distinguish isolates' divergence in shared sequence and gene content, which we demonstrate to be accurate over multiple orders of magnitude using data from both simulations and genomic collections representing 10 taxonomically widespread species. Connections between closely related isolates of the same strain are robustly identified, despite interspecies variation in the pairwise distance distributions that reflects species' diverse evolutionary patterns. PopPUNK can process 103-104 genomes in a single batch, with minimal memory use and runtimes up to 200-fold faster than existing model-based methods. Clusters of strains remain consistent as new batches of genomes are added, which is achieved without needing to reanalyze all genomes de novo. This facilitates real-time surveillance with consistent cluster naming between studies and allows for outbreak detection using hundreds of genomes in minutes. Interactive visualization and online publication is streamlined through the automatic output of results to multiple platforms. PopPUNK has been designed as a flexible platform that addresses important issues with currently used whole-genome clustering and typing methods, and has potential uses across bacterial genetics and public health research.

Persistent and Emerging Pneumococcal Carriage Serotypes in a Rural Gambian Community After 10 Years of Pneumococcal Conjugate Vaccine Pressure.

BACKGROUND: The continuing impact of pneumococcal conjugate vaccines (PCVs) in regions with high pneumococcal transmission is threatened by the persistence of vaccine serotypes (VTs) and the emergence of nonvaccine serotypes (NVTs). METHODS: In 2016, we conducted a cross-sectional carriage survey (CSS5) in a community where PCV7 was first introduced in 2006 during a cluster-randomized trial conducted before nationwide introduction of PCV7 (2009) and PCV13 (2011). We estimated prevalence of PCV13 VT and NVT by age and compared these with earlier surveys before (CSS0), during (CSS1-3), and after the trial but before PCV13 (CSS4). Genomic analysis was conducted for the nontypeable pneumococci. RESULTS: Prevalence of PCV13 VT carriage decreased during the 10 years between CSS0 and CSS5 across all age groups (67.6% to 13.5%, P < .001; 59.8% to 14.4%, P < .001; 43.1% to 17.9%, P < .001; and 24.0% to 5.1%, P < .001, in <2, 2-4, 5-14, and ≥15 years, respectively). However, there was no difference between CSS4 and CSS5 in children ≥2 years and adults (children <2 years, no data). The prevalence of PCV13 NVT increased between CSS0 and CSS5 for children <2 years but decreased in older children and adults. In CSS5, serotypes 3, 6A, and 19F were the most common VT and nontypeable isolates were the most common NVT. Among nontypeable isolates, 73.0% lost the ability to express a capsule. Of these, 70.8% were from a VT background. CONCLUSIONS: The decrease in PCV13 VT that has occurred since the introduction of PCV13 appears to have plateaued. Significant carriage of these serotypes remains in all age groups.

A Nationwide Outbreak of Invasive Pneumococcal Disease in Israel Caused by Streptococcus pneumoniae Serotype 2.

BACKGROUND: Invasive pneumococcal disease (IPD) caused by Streptococcus pneumoniae serotype 2 (Sp2) is infrequent. Large-scale outbreaks were not been reported following pneumococcal conjugate vaccine (PCV) implementation. We describe a Sp2 IPD outbreak in Israel, in the PCV13 era, with focus on Sp2 population structure and evolutionary dynamics. METHODS: The data were derived from a population-based, nationwide active surveillance of IPD since 2009. PCV7/PCV13 vaccines were introduced in July 2009 and November 2010, respectively. Sp2 isolates were tested for antimicrobial susceptibility, multilocus sequence typing, and whole-genome sequencing (WGS) analysis. RESULTS: Overall, 170 Sp2 IPD cases were identified during 2009-2019; Sp2 increased in 2015 and caused 6% of IPD during 2015-2019, a 7-fold increase compared with 2009-2014. The outbreak was caused by a previously unreported molecular type (ST-13578), initially observed in Israel in 2014. This clone caused 88% of Sp2 during 2015-2019. ST-13578 is a single-locus variant of ST-1504, previously reported globally including in Israel. WGS analysis confirmed clonality among the ST-13578 population. Single-nucleotide polymorphism-dense regions support a hypothesis that the ST-13578 outbreak clone evolved from ST-1504 by recombination. All tested strains were penicillin-susceptible (minimum inhibitory concentration <0.06 µg/mL). The ST-13578 clone was identified almost exclusively (99%) in the Jewish population and was mainly distributed in 3 of 7 Israeli districts. The outbreak is still ongoing, although it began declining in 2017. CONCLUSIONS: To the best of our knowledge, this is the first widespread Sp2 outbreak since PCV13 introduction worldwide, caused by the emerging ST-13578 clone.

Fast hierarchical Bayesian analysis of population structure.

We present fastbaps, a fast solution to the genetic clustering problem. Fastbaps rapidly identifies an approximate fit to a Dirichlet process mixture model (DPM) for clustering multilocus genotype data. Our efficient model-based clustering approach is able to cluster datasets 10-100 times larger than the existing model-based methods, which we demonstrate by analyzing an alignment of over 110 000 sequences of HIV-1 pol genes. We also provide a method for rapidly partitioning an existing hierarchy in order to maximize the DPM model marginal likelihood, allowing us to split phylogenetic trees into clades and subclades using a population genomic model. Extensive tests on simulated data as well as a diverse set of real bacterial and viral datasets show that fastbaps provides comparable or improved solutions to previous model-based methods, while being significantly faster. The method is made freely available under an open source MIT licence as an easy to use R package at https://github.com/gtonkinhill/fastbaps.

Genome-wide epistasis and co-selection study using mutual information.

Covariance-based discovery of polymorphisms under co-selective pressure or epistasis has received considerable recent attention in population genomics. Both statistical modeling of the population level covariation of alleles across the chromosome and model-free testing of dependencies between pairs of polymorphisms have been shown to successfully uncover patterns of selection in bacterial populations. Here we introduce a model-free method, SpydrPick, whose computational efficiency enables analysis at the scale of pan-genomes of many bacteria. SpydrPick incorporates an efficient correction for population structure, which adjusts for the phylogenetic signal in the data without requiring an explicit phylogenetic tree. We also introduce a new type of visualization of the results similar to the Manhattan plots used in genome-wide association studies, which enables rapid exploration of the identified signals of co-evolution. Simulations demonstrate the usefulness of our method and give some insight to when this type of analysis is most likely to be successful. Application of the method to large population genomic datasets of two major human pathogens, Streptococcus pneumoniae and Neisseria meningitidis, revealed both previously identified and novel putative targets of co-selection related to virulence and antibiotic resistance, highlighting the potential of this approach to drive molecular discoveries, even in the absence of phenotypic data.

Early Signals of Vaccine-driven Perturbation Seen in Pneumococcal Carriage Population Genomic Data.

BACKGROUND: Pneumococcal conjugate vaccines (PCVs) have reduced pneumococcal diseases globally. Pneumococcal genomic surveys elucidate PCV effects on population structure but are rarely conducted in low-income settings despite the high disease burden. METHODS: We undertook whole-genome sequencing (WGS) of 660 pneumococcal isolates collected through surveys from healthy carriers 2 years from 13-valent PCV (PCV13) introduction and 1 year after rollout in northern Malawi. We investigated changes in population structure, within-lineage serotype dynamics, serotype diversity, and frequency of antibiotic resistance (ABR) and accessory genes. RESULTS: In children <5 years of age, frequency and diversity of vaccine serotypes (VTs) decreased significantly post-PCV, but no significant changes occurred in persons ≥5 years of age. Clearance of VT serotypes was consistent across different genetic backgrounds (lineages). There was an increase of nonvaccine serotypes (NVTs)-namely 7C, 15B/C, and 23A-in children <5 years of age, but 28F increased in both age groups. While carriage rates have been recently shown to remain stable post-PCV due to replacement serotypes, there was no change in diversity of NVTs. Additionally, frequency of intermediate-penicillin-resistant lineages decreased post-PCV. Although frequency of ABR genes remained stable, other accessory genes, especially those associated with mobile genetic element and bacteriocins, showed changes in frequency post-PCV. CONCLUSIONS: We demonstrate evidence of significant population restructuring post-PCV driven by decreasing frequency of vaccine serotypes and increasing frequency of few NVTs mainly in children under 5. Continued surveillance with WGS remains crucial to fully understand dynamics of the residual VTs and replacement NVT serotypes post-PCV.