Development of an amplicon-based sequencing approach in response to the global emergence of mpox.

The 2022 multicountry mpox outbreak concurrent with the ongoing Coronavirus Disease 2019 (COVID-19) pandemic further highlighted the need for genomic surveillance and rapid pathogen whole-genome sequencing. While metagenomic sequencing approaches have been used to sequence many of the early mpox infections, these methods are resource intensive and require samples with high viral DNA concentrations. Given the atypical clinical presentation of cases associated with the outbreak and uncertainty regarding viral load across both the course of infection and anatomical body sites, there was an urgent need for a more sensitive and broadly applicable sequencing approach. Highly multiplexed amplicon-based sequencing (PrimalSeq) was initially developed for sequencing of Zika virus, and later adapted as the main sequencing approach for Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). Here, we used PrimalScheme to develop a primer scheme for human monkeypox virus that can be used with many sequencing and bioinformatics pipelines implemented in public health laboratories during the COVID-19 pandemic. We sequenced clinical specimens that tested presumptively positive for human monkeypox virus with amplicon-based and metagenomic sequencing approaches. We found notably higher genome coverage across the virus genome, with minimal amplicon drop-outs, in using the amplicon-based sequencing approach, particularly in higher PCR cycle threshold (Ct) (lower DNA titer) samples. Further testing demonstrated that Ct value correlated with the number of sequencing reads and influenced the percent genome coverage. To maximize genome coverage when resources are limited, we recommend selecting samples with a PCR Ct below 31 Ct and generating 1 million sequencing reads per sample. To support national and international public health genomic surveillance efforts, we sent out primer pool aliquots to 10 laboratories across the United States, United Kingdom, Brazil, and Portugal. These public health laboratories successfully implemented the human monkeypox virus primer scheme in various amplicon sequencing workflows and with different sample types across a range of Ct values. Thus, we show that amplicon-based sequencing can provide a rapidly deployable, cost-effective, and flexible approach to pathogen whole-genome sequencing in response to newly emerging pathogens. Importantly, through the implementation of our primer scheme into existing SARS-CoV-2 workflows and across a range of sample types and sequencing platforms, we further demonstrate the potential of this approach for rapid outbreak response.

Clonal Expansion of a Streptococcus pneumoniae Serotype 3 Capsule Variant Sequence Type 700 With Enhanced Vaccine Escape Potential After 13-Valent Pneumococcal Conjugate Vaccine Introduction.

BACKGROUND: Streptococcus pneumoniae serotype 3 remains a problem globally. Malawi introduced 13-valent pneumococcal conjugate vaccine (PCV13) in 2011, but there has been no direct protection against serotype 3 carriage. We explored whether vaccine escape by serotype 3 is due to clonal expansion of a lineage with a competitive advantage. METHODS: The distribution of serotype 3 Global Pneumococcal Sequence Clusters (GPSCs) and sequence types (STs) globally was assessed using sequences from the Global Pneumococcal Sequencing Project. Whole-genome sequences of 135 serotype 3 carriage isolates from Blantyre, Malawi (2015-2019) were analyzed. Comparative analysis of the capsule locus, entire genomes, antimicrobial resistance, and phylogenetic reconstructions were undertaken. Opsonophagocytosis was evaluated using serum samples from vaccinated adults and children. RESULTS: Serotype 3 GPSC10-ST700 isolates were most prominent in Malawi. Compared with the prototypical serotype 3 capsular polysaccharide locus sequence, 6 genes are absent, with retention of capsule polysaccharide biosynthesis. This lineage is characterized by increased antimicrobial resistance and lower susceptibility to opsonophagocytic killing. CONCLUSIONS: A serotype 3 variant in Malawi has genotypic and phenotypic characteristics that could enhance vaccine escape and clonal expansion after post-PCV13 introduction. Genomic surveillance among high-burden populations is essential to improve the effectiveness of next-generation pneumococcal vaccines.

DengueSeq: a pan-serotype whole genome amplicon sequencing protocol for dengue virus.

BACKGROUND: The increasing burden of dengue virus on public health due to more explosive and frequent outbreaks highlights the need for improved surveillance and control. Genomic surveillance of dengue virus not only provides important insights into the emergence and spread of genetically diverse serotypes and genotypes, but it is also critical to monitor the effectiveness of newly implemented control strategies. Here, we present DengueSeq, an amplicon sequencing protocol, which enables whole-genome sequencing of all four dengue virus serotypes. RESULTS: We developed primer schemes for the four dengue virus serotypes, which can be combined into a pan-serotype approach. We validated both approaches using genetically diverse virus stocks and clinical specimens that contained a range of virus copies. High genome coverage (>95%) was achieved for all genotypes, except DENV2 (genotype VI) and DENV 4 (genotype IV) sylvatics, with similar performance of the serotype-specific and pan-serotype approaches. The limit of detection to reach 70% coverage was 10-100 RNA copies/µL for all four serotypes, which is similar to other commonly used primer schemes. DengueSeq facilitates the sequencing of samples without known serotypes, allows the detection of multiple serotypes in the same sample, and can be used with a variety of library prep kits and sequencing instruments. CONCLUSIONS: DengueSeq was systematically evaluated with virus stocks and clinical specimens spanning the genetic diversity within each of the four dengue virus serotypes. The primer schemes can be plugged into existing amplicon sequencing workflows to facilitate the global need for expanded dengue virus genomic surveillance.

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.

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.

Rotavirus Genotypes in Hospitalized Children With Acute Gastroenteritis Before and After Rotavirus Vaccine Introduction in Blantyre, Malawi, 1997-2019.

BACKGROUND: Rotavirus vaccine (Rotarix [RV1]) has reduced diarrhea-associated hospitalizations and deaths in Malawi. We examined the trends in circulating rotavirus genotypes in Malawi over a 22-year period to assess the impact of RV1 introduction on strain distribution. METHODS: Data on rotavirus-positive stool specimens among children aged <5 years hospitalized with diarrhea in Blantyre, Malawi before (July 1997-October 2012, n = 1765) and after (November 2012-October 2019, n = 934) RV1 introduction were analyzed. Rotavirus G and P genotypes were assigned using reverse-transcription polymerase chain reaction. RESULTS: A rich rotavirus strain diversity circulated throughout the 22-year period; Shannon (H') and Simpson diversity (D') indices did not differ between the pre- and postvaccine periods (H' P < .149; D' P < .287). Overall, G1 (n = 268/924 [28.7%]), G2 (n = 308/924 [33.0%]), G3 (n = 72/924 [7.7%]), and G12 (n = 109/924 [11.8%]) were the most prevalent genotypes identified following RV1 introduction. The prevalence of G1P[8] and G2P[4] genotypes declined each successive year following RV1 introduction, and were not detected after 2018. Genotype G3 reemerged and became the predominant genotype from 2017 onward. No evidence of genotype selection was observed 7 years post-RV1 introduction. CONCLUSIONS: Rotavirus strain diversity and genotype variation in Malawi are likely driven by natural mechanisms rather than vaccine pressure.

Daily Rapid Antigen Exit Testing to Tailor University COVID-19 Isolation Policy.

We evaluated daily rapid antigen test (RAT) data from 323 COVID-19-positive university students in Connecticut, USA, during an Omicron-dominant period. Day 5 positivity was 47% for twice-weekly screeners and 26%-28% for less-frequent screeners, approximately halving each subsequent day. Testing negative >10 days before diagnosis (event time ratio (ETR) 0.85 [95% CI 0.75-0.96]) and prior infection >90 days (ETR 0.50 [95% CI 0.33-0.76]) were significantly associated with shorter RAT positivity duration. Symptoms before or at diagnosis (ETR 1.13 [95% CI 1.02-1.25]) and receipt of 3 vaccine doses (ETR 1.20 [95% CI 1.04-1.39]) were significantly associated with prolonged positivity. Exit RATs enabled 53%-74% of students to leave isolation early when they began isolation at the time of the first positive test, but 15%-22% remained positive beyond the recommended isolation period. Factors associated with RAT positivity duration should be further explored to determine relationships with infection duration.

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.

Evaluation of Dried Blood and Cerebrospinal Fluid Filter Paper Spots for Storing and Transporting Clinical Material for the Molecular Diagnosis of Invasive Meningococcal Disease.

To improve the storage and transport of clinical specimens for the diagnosis of Neisseria meningitidis (Nm) infections in resource-limited settings, we have evaluated the performance of dried blood spot (DBS) and dried cerebrospinal fluid spot (DCS) assays. DBS and DCS were prepared on filter paper from liquid specimens previously tested for Nm in the United Kingdom. Nm was detected and genogrouped by real-time PCR performed on crude genomic DNA extracted from the DBS (n = 226) and DCS (n = 226) specimens. Targeted whole-genome sequencing was performed on a subset of specimens, DBS (n = 4) and DCS (n = 6). The overall agreement between the analysis of liquid and dried specimens was (94.2%; 95% CI 90.8−96.7) for blood and (96.4%; 95% CI 93.5−98.0) for cerebrospinal fluid. Relative to liquid specimens as the reference, the DBS and DCS assays had sensitivities of (89.1%; 95% CI 82.7−93.8) and (94.2%; 95% CI 88.9−97.5), respectively, and both assays had specificities above 98%. A genogroup was identified by dried specimen analysis for 81.9% of the confirmed meningococcal infections. Near full-length Nm genome sequences (>86%) were obtained for all ten specimens tested which allowed determination of the sequence type, clonal complex, presence of antimicrobial resistance and other meningococcal genotyping. Dried blood and CSF filter spot assays offer a practical alternative to liquid specimens for the molecular and genomic characterisation of invasive meningococcal diseases in low-resource settings.

Whole genomic comparative analysis of Streptococcus pneumoniae serotype 1 isolates causing invasive and non-invasive infections among children under 5 years in Casablanca, Morocco.

BACKGROUND: Streptococcus pneumoniae serotype 1 remains a leading cause of invasive pneumococcal diseases, even in countries with PCV-10/PCV-13 vaccine implementation. The main objective of this study, which is part of the Pneumococcal African Genome project (PAGe), was to determine the phylogenetic relationships of serotype 1 isolates recovered from children patients in Casablanca (Morocco), compared to these from other African countries; and to investigate the contribution of accessory genes and recombination events to the genetic diversity of this serotype. RESULTS: The genome average size of the six-pneumococcus serotype 1 from Casablanca was 2,227,119 bp, and the average content of coding sequences was 2113, ranging from 2041 to 2161. Pangenome analysis of the 80 genomes used in this study revealed 1685 core genes and 1805 accessory genes. The phylogenetic tree based on core genes and the hierarchical bayesian clustering analysis revealed five sublineages with a phylogeographic structure by country. The Moroccan strains cluster in two different lineages, the five invasive strains clusters altogether in a divergent clade distantly related to the non-invasive strain, that cluster with all the serotype 1 genomes from Africa. CONCLUSIONS: The whole genome sequencing provides increased resolution analysis of the highly virulent serotype 1 in Casablanca, Morocco. Our results are concordant with previous works, showing that the phylogeography of S. pneumoniae serotype 1 is structured by country, and despite the small size (six isolates) of the Moroccan sample, our analysis shows the genetic cohesion of the Moroccan invasive isolates.

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.

Genomic analysis of Klebsiella pneumoniae isolates from Malawi reveals acquisition of multiple ESBL determinants across diverse lineages.

OBJECTIVES: ESBL-producing Klebsiella pneumoniae (KPN) pose a major threat to human health globally. We carried out a WGS study to understand the genetic background of ESBL-producing KPN in Malawi and place them in the context of other global isolates. METHODS: We sequenced genomes of 72 invasive and carriage KPN isolates collected from patients admitted to Queen Elizabeth Central Hospital, Blantyre, Malawi. We performed phylogenetic and population structure analyses on these and previously published genomes from Kenya (n = 66) and from outside sub-Saharan Africa (n = 67). We screened for presence of antimicrobial resistance (AMR) genetic determinants and carried out association analyses by genomic sequence cluster, AMR phenotype and time. RESULTS: Malawian isolates fit within the global population structure of KPN, clustering into the major lineages of KpI, KpII and KpIII. KpI isolates from Malawi were more related to those from Kenya, with both collections exhibiting more clonality than isolates from the rest of the world. We identified multiple ESBL genes, including blaCTX-M-15, several blaSHV, blaTEM-63 and blaOXA-10, and other AMR genes, across diverse lineages of the KPN isolates from Malawi. No carbapenem resistance genes were detected; however, we detected IncFII and IncFIB plasmids that were similar to the carbapenem resistance-associated plasmid pNDM-mar. CONCLUSIONS: There are multiple ESBL genes across diverse KPN lineages in Malawi and plasmids in circulation that are capable of carrying carbapenem resistance. Unless appropriate interventions are rapidly put in place, these may lead to a high burden of locally untreatable infection in vulnerable populations.

Emergence of Double- and Triple-Gene Reassortant G1P[8] Rotaviruses Possessing a DS-1-Like Backbone after Rotavirus Vaccine Introduction in Malawi.

To combat the high burden of rotavirus gastroenteritis, multiple African countries have introduced rotavirus vaccines into their childhood immunization programs. Malawi incorporated a G1P[8] rotavirus vaccine (Rotarix) into its immunization schedule in 2012. Utilizing a surveillance platform of hospitalized rotavirus gastroenteritis cases, we examined the phylodynamics of G1P[8] rotavirus strains that circulated in Malawi before (1998 to 2012) and after (2013 to 2014) vaccine introduction. Analysis of whole genomes obtained through next-generation sequencing revealed that all randomly selected prevaccine G1P[8] strains sequenced (n = 32) possessed a Wa-like genetic constellation, whereas postvaccine G1P[8] strains (n = 18) had a DS-1-like constellation. Phylodynamic analyses indicated that postvaccine G1P[8] strains emerged through reassortment events between human Wa- and DS-1-like rotaviruses that circulated in Malawi from the 1990s and hence were classified as atypical DS-1-like reassortants. The time to the most recent common ancestor for G1P[8] strains was from 1981 to 1994; their evolutionary rates ranged from 9.7 × 10-4 to 4.1 × 10-3 nucleotide substitutions/site/year. Three distinct G1P[8] lineages chronologically replaced each other between 1998 and 2014. Genetic drift was the likely driver for lineage turnover in 2005, whereas replacement in 2013 was due to reassortment. Amino acid substitution within the outer glycoprotein VP7 of G1P[8] strains had no impact on the structural conformation of the antigenic regions, suggesting that it is unlikely that they would affect recognition by vaccine-induced neutralizing antibodies. While the emergence of DS-1-like G1P[8] rotavirus reassortants in Malawi was therefore likely due to natural genotype variation, vaccine effectiveness against such strains needs careful evaluation.IMPORTANCE The error-prone RNA-dependent RNA polymerase and the segmented RNA genome predispose rotaviruses to genetic mutation and genome reassortment, respectively. These evolutionary mechanisms generate novel strains and have the potential to lead to the emergence of vaccine escape mutants. While multiple African countries have introduced a rotavirus vaccine, there are few data describing the evolution of rotaviruses that circulated before and after vaccine introduction. We report the emergence of atypical DS-1-like G1P[8] strains during the postvaccine era in Malawi. Three distinct G1P[8] lineages circulated chronologically from 1998 to 2014; mutation and reassortment drove lineage turnover in 2005 and 2013, respectively. Amino acid substitutions within the outer capsid VP7 glycoprotein did not affect the structural conformation of mapped antigenic sites, suggesting a limited effect on the recognition of G1-specific vaccine-derived antibodies. The genes that constitute the remaining genetic backbone may play important roles in immune evasion, and vaccine effectiveness against such atypical strains needs careful evaluation.

Comparative Genomic Analysis and In Vivo Modeling of Streptococcus pneumoniae ST3081 and ST618 Isolates Reveal Key Genetic and Phenotypic Differences Contributing to Clonal Replacement of Serotype 1 in The Gambia.

Streptococcus pneumoniae serotype 1 is one of the leading causes of invasive pneumococcal disease (IPD) in West Africa, with ST618 being the dominant cause of IPD in The Gambia. Recently however, a rare example of clonal replacement was observed, where the ST3081 clone of serotype 1 replaced the predominant ST618 clone as the main cause of IPD. In the current study, we sought to find the reasons for this unusual replacement event. Using whole-genome sequence analysis and clinically relevant models of in vivo infection, we identified distinct genetic and phenotypic characteristics of the emerging ST3081 clone. We show that ST3081 is significantly more virulent than ST618 in models of invasive pneumonia, and is carried at higher densities than ST618 during nasopharyngeal carriage. We also observe sequence type-specific accessory genes and a unique sequence type-specific fixed mutation in the pneumococcal toxin pneumolysin, which is associated with increased hemolytic activity in ST3081 and may contribute to increased virulence in this clone. Our study provides evidence that, within the same serotype 1 clonal complex, biological properties differ significantly from one clone to another in terms of virulence and host invasiveness, and that these differences may be the result of key genetic differences within the genome.

Genomic Epidemiology of Penicillin-Nonsusceptible Pneumococci with Nonvaccine Serotypes Causing Invasive Disease in the United States.

Conjugate vaccination against seven pneumococcal serotypes (PCV7) reduced disease prevalence due to antibiotic-resistant strains throughout the 2000s. However, diseases caused by resistant nonvaccine type (NVT) strains increased. Some of these emerging strains were derived from vaccine types (VT) that had changed their capsule by recombination. The introduction of a vaccine targeting 13 serotypes (PCV13) in 2010 has led to concern that this scenario will repeat itself. We generated high-quality draft genomes from 265 isolates of NVT pneumococci not susceptible to penicillin (PNSP) in 2009 and compared them with the genomes of 581 isolates from 2012 to 2013 collected by the Active Bacterial Core surveillance (ABCs) of the Centers for Disease Control and Prevention (CDC). Of the seven sequence clusters (SCs) identified, three SCs fell into a single lineage associated with serogroup 23, which had an origin in 1908 as dated by coalescent analysis and included isolates with a divergent 23B capsule locus. Three other SCs represented relatively deep-branching lineages associated with serotypes 35B, 15A, and 15BC. In all cases, the resistant clones originated prior to 2010, indicating that PNSP are at present dominated by descendants of NVT clones present before vaccination. With one exception (15BC/ST3280), these SCs were related to clones identified by the Pneumococcal Molecular Epidemiology Network (PMEN). We conclude that postvaccine diversity in NVT PNSP between 2009 and 2013 was driven mainly by the persistence of preexisting strains rather than through de novo adaptation, with few cases of serotype switching. Future surveillance is essential for documenting the long-term dynamics and resistance of NVT PNSP.

Genomic landscape of extended-spectrum ß-lactamase resistance in Escherichia coli from an urban African setting.

Objectives: Efforts to treat Escherichia coli infections are increasingly being compromised by the rapid, global spread of antimicrobial resistance (AMR). Whilst AMR in E. coli has been extensively investigated in resource-rich settings, in sub-Saharan Africa molecular patterns of AMR are not well described. In this study, we have begun to explore the population structure and molecular determinants of AMR amongst E. coli isolates from Malawi. Methods: Ninety-four E. coli isolates from patients admitted to Queen's Hospital, Malawi, were whole-genome sequenced. The isolates were selected on the basis of diversity of phenotypic resistance profiles and clinical source of isolation (blood, CSF and rectal swab). Sequence data were analysed using comparative genomics and phylogenetics. Results: Our results revealed the presence of five clades, which were strongly associated with E. coli phylogroups A, B1, B2, D and F. We identified 43 multilocus STs, of which ST131 (14.9%) and ST12 (9.6%) were the most common. We identified 25 AMR genes. The most common ESBL gene was bla CTX-M-15 and it was present in all five phylogroups and 11 STs, and most commonly detected in ST391 (4/4 isolates), ST648 (3/3 isolates) and ST131 [3/14 (21.4%) isolates]. Conclusions: This study has revealed a high diversity of lineages associated with AMR, including ESBL and fluoroquinolone resistance, in Malawi. The data highlight the value of longitudinal bacteraemia surveillance coupled with detailed molecular epidemiology in all settings, including low-income settings, in describing the global epidemiology of ESBL resistance.

Understanding pneumococcal serotype 1 biology through population genomic analysis.

BACKGROUND: Pneumococcus kills over one million children annually and over 90 % of these deaths occur in low-income countries especially in Sub-Saharan Africa (SSA) where HIV exacerbates the disease burden. In SSA, serotype 1 pneumococci particularly the endemic ST217 clone, causes majority of the pneumococcal disease burden. To understand the evolution of the virulent ST217 clone, we analysed ST217 whole genomes from isolates sampled from African and Asian countries. METHODS: We analysed 226 whole genome sequences from the ST217 lineage sampled from 9 African and 4 Asian countries. We constructed a whole genome alignment and used it for phylogenetic and coalescent analyses. We also screened the genomes to determine presence of antibiotic resistance conferring genes. RESULTS: Population structure analysis grouped the ST217 isolates into five sequence clusters (SCs), which were highly associated with different geographical regions and showed limited intracontinental and intercontinental spread. The SCs showed lower than expected genomic sequence, which suggested strong purifying selection and small population sizes caused by bottlenecks. Recombination rates varied between the SCs but were lower than in other successful clones such as PMEN1. African isolates showed higher prevalence of antibiotic resistance genes than Asian isolates. Interestingly, certain West African isolates harbored a defective chloramphenicol and tetracycline resistance-conferring element (Tn5253) with a deletion in the loci encoding the chloramphenicol resistance gene (cat pC194), which caused lower chloramphenicol than tetracycline resistance. Furthermore, certain genes that promote colonisation were absent in the isolates, which may contribute to serotype 1's rarity in carriage and consequently its lower recombination rates. CONCLUSIONS: The high phylogeographic diversity of the ST217 clone shows that this clone has been in circulation globally for a long time, which allowed its diversification and adaptation in different geographical regions. Such geographic adaptation reflects local variations in selection pressures in different locales. Further studies will be required to fully understand the biological mechanisms which makes the ST217 clone highly invasive but unable to successfully colonise the human nasopharynx for long durations which results in lower recombination rates.

Comparative Genomic Analysis of Meningitis- and Bacteremia-Causing Pneumococci Identifies a Common Core Genome.

Streptococcus pneumoniae is a nasopharyngeal commensal that occasionally invades normally sterile sites to cause bloodstream infection and meningitis. Although the pneumococcal population structure and evolutionary genetics are well defined, it is not clear whether pneumococci that cause meningitis are genetically distinct from those that do not. Here, we used whole-genome sequencing of 140 isolates of S. pneumoniae recovered from bloodstream infection (n = 70) and meningitis (n = 70) to compare their genetic contents. By fitting a double-exponential decaying-function model, we show that these isolates share a core of 1,427 genes (95% confidence interval [CI], 1,425 to 1,435 genes) and that there is no difference in the core genome or accessory gene content from these disease manifestations. Gene presence/absence alone therefore does not explain the virulence behavior of pneumococci that reach the meninges. Our analysis, however, supports the requirement of a range of previously described virulence factors and vaccine candidates for both meningitis- and bacteremia-causing pneumococci. This high-resolution view suggests that, despite considerable competency for genetic exchange, all pneumococci are under considerable pressure to retain key components advantageous for colonization and transmission and that these components are essential for access to and survival in sterile sites.

High multiple carriage and emergence of Streptococcus pneumoniae vaccine serotype variants in Malawian children.

BACKGROUND: Carriage of either single or multiple pneumococcal serotypes (multiple carriage) is a prerequisite for developing invasive pneumococcal disease. However, despite the reported high rates of pneumococcal carriage in Malawi, no data on carriage of multiple serotypes has been reported previously. Our study provides the first description of the prevalence of multiple pneumococcal carriage in Malawi. METHODS: The study was conducted in Blantyre and Karonga districts in Malawi, from 2008 to 2012. We recruited 116 children aged 0-13 years. These children were either HIV-infected (N = 44) or uninfected (N = 72). Nasopharyngeal samples were collected using sterile swabs. Pneumococcal serotypes in the samples were identified by microarray. Strains that could not be typed by microarray were sequenced to characterise possible genetic alterations within the capsular polysaccharide (CPS) locus. RESULTS: The microarray identified 179 pneumococcal strains (from 116 subjects), encompassing 43 distinct serotypes and non-typeable (NT) strains. Forty per cent (46/116) of children carried multiple serotypes. Carriage of vaccine type (VT) strains was higher (p = 0.028) in younger (0-2 years) children (71 %, 40/56) compared to older (3-13 years) children (50 %, 30/60). Genetic variations within the CPS locus of known serotypes were observed in 19 % (34/179) of the strains identified. The variants included 13-valent pneumococcal conjugate vaccine (PCV13) serotypes 6B and 19A, and the polysaccharide vaccine serotype 20. Serotype 6B variants were the most frequently isolated (47 %, 16/34). Unlike the wild type, the CPS locus of the 6B variants contained an insertion of the licD-family phosphotransferase gene. The CPS locus of 19A- and 20-variants contained an inversion in the sugar-biosynthesis (rmlD) gene and a 717 bp deletion within the transferase (whaF) gene, respectively. CONCLUSIONS: The high multiple carriage in Malawian children provides opportunities for genetic exchange through horizontal gene transfer. This may potentially lead to CPS locus variants and vaccine escape. Variants reported here occurred naturally, however, PCV13 introduction could exacerbate the CPS genetic variations. Further studies are therefore recommended to assess the invasive potential of these variants and establish whether PCV13 would offer cross-protection. We have shown that younger children (0-2 years) are a reservoir of VT serotypes, which makes them an ideal target for vaccination.

Expansion of pneumococcal serotype 23F and 14 lineages with genotypic changes in capsule polysaccharide locus and virulence gene profiles post introduction of pneumococcal conjugate vaccine in Blantyre, Malawi.

Since the introduction of the 13-valent pneumococcal conjugate vaccine (PCV13) in Malawi in 2011, there has been persistent carriage of vaccine serotype (VT) Streptococcus pneumoniae, despite high vaccine coverage. To determine if there has been a genetic change within the VT capsule polysaccharide (cps) loci since the vaccine's introduction, we compared 1022 whole-genome-sequenced VT isolates from 1998 to 2019. We identified the clonal expansion of a multidrug-resistant, penicillin non-susceptible serotype 23F GPSC14-ST2059 lineage, a serotype 14 GPSC9-ST782 lineage and a novel serotype 14 sequence type GPSC9-ST18728 lineage. Serotype 23F GPSC14-ST2059 had an I253T mutation within the capsule oligosaccharide repeat unit polymerase Wzy protein, which is predicted in silico to alter the protein pocket cavity. Moreover, serotype 23F GPSC14-ST2059 had SNPs in the DNA binding sites for the cps transcriptional repressors CspR and SpxR. Serotype 14 GPSC9-ST782 harbours a non-truncated version of the large repetitive protein (Lrp), containing a Cna protein B-type domain which is also present in proteins associated with infection and colonisation. These emergent lineages also harboured genes associated with antibiotic resistance, and the promotion of colonisation and infection which were absent in other lineages of the same serotype. Together these data suggest that in addition to serotype replacement, modifications of the capsule locus associated with changes in virulence factor expression and antibiotic resistance may promote vaccine escape. In summary, the study highlights that the persistence of vaccine serotype carriage despite high vaccine coverage in Malawi may be partly caused by expansion of VT lineages post-PCV13 rollout.