Abstract 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 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.

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.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.

Metagenomic sequencing analysis is central to investigating microbial communities in clinical and environmental studies. Short read sequencing remains the primary data type for metagenomic research. However, long read sequencing promises the advantages of improved metagenomic assembly and resolved taxonomic identification. To assess the comparative performance of short and long read sequencing data for metagenomic analysis, we simulated short and long read datasets using increasingly complex metagenomes of 10, 20, and 50 microbial taxa. In addition, an empirical dataset of paired short and long read data from mouse fecal pellets was generated to assess feasibility. We compared metagenomic assembly quality, taxonomic classification capabilities, and metagenome-assembled genome recovery rates for both simulated and real metagenomic sequence data. We show that long read sequencing data significantly outperforms short read sequencing. For simulated long read datasets, metagenomic assemblies were completer and more contiguous with higher rates of metagenome-assembled genome recovery. This resulted in more precise taxonomic classifications and increased accuracy of relative abundance estimates. Analysis of empirical data demonstrated that sequencing technology directly affects compositional results. Overall, we highlight strengths of long read sequencing for metagenomic research of microbial communities over traditional short read approaches.

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 then assessed intrahost diversity, population structure, evolutionary history, and signatures of adaptation in the context of patient age, antibiotic treatment, and co-colonizing microbes. We conducted tests of gene-wide diversifying selection and identified instances of switched intergenic regions (IGRs), which may be associated with gene expression. 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. Coalescent analysis revealed rates of molecular evolution ranging from 2.21 to 8.64 nucleotide polymorphisms per genome per year. Estimated lineage ages were consistent with acquisition of colonization in early childhood and subsequent persistence of multiple sub-populations. Selection analysis focused on 1,622 core genes shared by all four clonal populations (n=79). Eleven genes were variable in three subjects - most notable, 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 ecology of microbes in the airways of CF patients is undeniably complex, posing 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.Contribution to the Field 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. We found that each patient possessed a single predominant strain in addition to transient non-clonal isolates. Considerable genomic diversity was observed among intrahost populations including variation in antibiotic resistance, even among isolates collected on the same day. This finding has implications for treatment. Estimated lineage ages were consistent with acquisition of colonization in early childhood and subsequent persistence of multiple sub-populations. The scarcity of non-synonymous substitutions suggests chronic S. aureus carriage in CF patients provides a sufficient opportunity for purifying selection to act. Eleven genes were found to be variable in three patients and present possible candidates for loci experiencing gene-wide diversifying selection. Overall, we 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.

Background: Streptococcus pneumoniae is a human commensal and the causative agent of pneumococcal disease. Pneumococci are naturally competent able to uptake exogenous DNA from the environment and incorporate it into their genome through homologous and non-homologous recombination. Recombination has significantly shaped the evolutionary history of S. pneumoniae, as it allows pneumococci to rapidly adapt. Recombination frequencies vary considerably across pneumococcal populations, yet underlying mechanisms for these variations are not well understood. Entry and exit into competence, a state in which the cell can uptake DNA, is tightly regulated through transcriptional changes. Methods: To elucidate differences in transformation frequency among strains as well as the underlying genetic mechanisms, we carried out in-vitro competence assays and measured gene expression changes during the competent state using RNA sequencing of representative strains belonging to Serotype 3 clonal complex (CC) 180 and a non-CC180 comparison. Results: We observed consistent differences in transformation frequencies among groups, which correlated with variation in differentially expressed genes during competence. While all strains exhibited a similar response to competence stimulating peptide (CSP) for early competence genes, we observed variation in expression of late competence genes, which encode the DNA uptake apparatus, DNA repair and recombination proteins. We also observed differences in expression of genes linked to bacteriocin production, which may partially explain observed population-level differences. Further genomic analysis suggests variation in promoter sequences governing late competence genes may be slowing transition from early to late components of the competence pathway. Conclusions: We show that there is considerable variation in competence even among closely related strains belonging to the same CC and that this variation may be the result of subtle genomic differences. Additional studies are needed to assess the phenotypic impact of genomic variations.

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 United States, 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 identify 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 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 13-valent pneumococcal conjugate vaccine (PCV-13) was introduced in the United States in 2010. Using a large pediatric 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 carriages isolates from children in Massachusetts between 2001 and 2014 were analyzed. 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% 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 geographic area, which may limit the generalizability 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.