Detecting co-selection through excess linkage disequilibrium in bacterial genomes.

Population genomics has revolutionized our ability to study bacterial evolution by enabling data-driven discovery of the genetic architecture of trait variation. Genome-wide association studies (GWAS) have more recently become accompanied by genome-wide epistasis and co-selection (GWES) analysis, which offers a phenotype-free approach to generating hypotheses about selective processes that simultaneously impact multiple loci across the genome. However, existing GWES methods only consider associations between distant pairs of loci within the genome due to the strong impact of linkage-disequilibrium (LD) over short distances. Based on the general functional organisation of genomes it is nevertheless expected that majority of co-selection and epistasis will act within relatively short genomic proximity, on co-variation occurring within genes and their promoter regions, and within operons. Here, we introduce LDWeaver, which enables an exhaustive GWES across both short- and long-range LD, to disentangle likely neutral co-variation from selection. We demonstrate the ability of LDWeaver to efficiently generate hypotheses about co-selection using large genomic surveys of multiple major human bacterial pathogen species and validate several findings using functional annotation and phenotypic measurements. Our approach will facilitate the study of bacterial evolution in the light of rapidly expanding population genomic data.

Modulation of multidrug-resistant clone success in Escherichia coli populations: a longitudinal, multi-country, genomic and antibiotic usage cohort study.

BACKGROUND: The effect of antibiotic usage on the success of multidrug-resistant (MDR) clones in a population remains unclear. With this genomics-based molecular epidemiology study, we aimed to investigate the contribution of antibiotic use to Escherichia coli clone success, relative to intra-strain competition for colonisation and infection. METHODS: We sequenced all the available E coli bloodstream infection isolates provided by the British Society for Antimicrobial Chemotherapy (BSAC) from 2012 to 2017 (n=718) and combined these with published data from the UK (2001-11; n=1090) and Norway (2002-17; n=3254). Defined daily dose (DDD) data from the European Centre for Disease Prevention and Control (retrieved on Sept 21, 2021) for major antibiotic classes (ß-lactam, tetracycline, macrolide, sulfonamide, quinolone, and non-penicillin ß-lactam) were used together with sequence typing, resistance profiling, regression analysis, and non-neutral Wright-Fisher simulation-based modelling to enable systematic comparison of resistance levels, clone success, and antibiotic usage between the UK and Norway. FINDINGS: Sequence type (ST)73, ST131, ST95, and ST69 accounted for 892 (49·3%) of 1808 isolates in the BSAC collection. In the UK, the proportion of ST69 increased between 2001-10 and 2011-17 (p=0·0004), whereas the proportions of ST73 and ST95 did not vary between periods. ST131 expanded quickly after its emergence in 2003 and its prevalence remained consistent throughout the study period (apart from a brief decrease in 2009-10). The extended-spectrum ß-lactamase (ESBL)-carrying, globally disseminated MDR clone ST131-C2 showed overall greater success in the UK (154 [56·8%] of 271 isolates in 2003-17) compared with Norway (51 [18·3%] of 278 isolates in 2002-17; p<0·0001). DDD data indicated higher total use of antimicrobials in the UK, driven mainly by the class of non-penicillin ß-lactams, which were used between 2·7-times and 5·1-times more in the UK per annum (ratio mean 3·7 [SD 0·8]). This difference was associated with the higher success of the MDR clone ST131-C2 (pseudo-R2 69·1%). A non-neutral Wright-Fisher model replicated the observed expansion of non-MDR and MDR sequence types under higher DDD regimes. INTERPRETATION: Our study indicates that resistance profiles of contemporaneously successful clones can vary substantially, warranting caution in the interpretation of correlations between aggregate measures of resistance and antibiotic usage. Our study further suggests that in countries with low-to-moderate use of antibiotics, such as the UK and Norway, the extent of non-penicillin ß-lactam use modulates rather than determines the success of widely disseminated MDR ESBL-carrying E coli clones. Detailed understanding of underlying causal drivers of success is important for improved control of resistant pathogens. FUNDING: Trond Mohn Foundation, Marie Sklodowska-Curie Actions, European Research Council, Royal Society, and Wellcome Trust.

Mge-cluster: a reference-free approach for typing bacterial plasmids.

Extrachromosomal elements of bacterial cells such as plasmids are notorious for their importance in evolution and adaptation to changing ecology. However, high-resolution population-wide analysis of plasmids has only become accessible recently with the advent of scalable long-read sequencing technology. Current typing methods for the classification of plasmids remain limited in their scope which motivated us to develop a computationally efficient approach to simultaneously recognize novel types and classify plasmids into previously identified groups. Here, we introduce mge-cluster that can easily handle thousands of input sequences which are compressed using a unitig representation in a de Bruijn graph. Our approach offers a faster runtime than existing algorithms, with moderate memory usage, and enables an intuitive visualization, classification and clustering scheme that users can explore interactively within a single framework. Mge-cluster platform for plasmid analysis can be easily distributed and replicated, enabling a consistent labelling of plasmids across past, present, and future sequence collections. We underscore the advantages of our approach by analysing a population-wide plasmid data set obtained from the opportunistic pathogen Escherichia coli, studying the prevalence of the colistin resistance gene mcr-1.1 within the plasmid population, and describing an instance of resistance plasmid transmission within a hospital environment.

Evolutionary and functional history of the Escherichia coli K1 capsule.

Escherichia coli is a leading cause of invasive bacterial infections in humans. Capsule polysaccharide has an important role in bacterial pathogenesis, and the K1 capsule has been firmly established as one of the most potent capsule types in E. coli through its association with severe infections. However, little is known about its distribution, evolution and functions across the E. coli phylogeny, which is fundamental to elucidating its role in the expansion of successful lineages. Using systematic surveys of invasive E. coli isolates, we show that the K1-cps locus is present in a quarter of bloodstream infection isolates and has emerged in at least four different extraintestinal pathogenic E. coli (ExPEC) phylogroups independently in the last 500 years. Phenotypic assessment demonstrates that K1 capsule synthesis enhances E. coli survival in human serum independent of genetic background, and that therapeutic targeting of the K1 capsule re-sensitizes E. coli from distinct genetic backgrounds to human serum. Our study highlights that assessing the evolutionary and functional properties of bacterial virulence factors at population levels is important to better monitor and predict the emergence of virulent clones, and to also inform therapies and preventive medicine to effectively control bacterial infections whilst significantly lowering antibiotic usage.

Inferring the heritability of bacterial traits in the era of machine learning.

Quantification of heritability is a fundamental desideratum in genetics, which allows an assessment of the contribution of additive genetic variation to the variability of a trait of interest. The traditional computational approaches for assessing the heritability of a trait have been developed in the field of quantitative genetics. However, the rise of modern population genomics with large sample sizes has led to the development of several new machine learning-based approaches to inferring heritability. In this article, we systematically summarize recent advances in machine learning which can be used to infer heritability. We focus on an application of these methods to bacterial genomes, where heritability plays a key role in understanding phenotypes such as antibiotic resistance and virulence, which are particularly important due to the rising frequency of antimicrobial resistance. By designing a heritability model incorporating realistic patterns of genome-wide linkage disequilibrium for a frequently recombining bacterial pathogen, we test the performance of a wide spectrum of different inference methods, including also GCTA. In addition to the synthetic data benchmark, we present a comparison of the methods for antibiotic resistance traits for multiple bacterial pathogens. Insights from the benchmarking and real data analyses indicate a highly variable performance of the different methods and suggest that heritability inference would likely benefit from tailoring of the methods to the specific genetic architecture of the target organism. Availability and implementation: The R codes and data used in the numerical experiments are available at: https://github.com/tienmt/her_MLs.

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.

Strong pathogen competition in neonatal gut colonisation.

Opportunistic bacterial pathogen species and their strains that colonise the human gut are generally understood to compete against both each other and the commensal species colonising this ecosystem. Currently we are lacking a population-wide quantification of strain-level colonisation dynamics and the relationship of colonisation potential to prevalence in disease, and how ecological factors might be modulating these. Here, using a combination of latest high-resolution metagenomics and strain-level genomic epidemiology methods we performed a characterisation of the competition and colonisation dynamics for a longitudinal cohort of neonatal gut microbiomes. We found strong inter- and intra-species competition dynamics in the gut colonisation process, but also a number of synergistic relationships among several species belonging to genus Klebsiella, which includes the prominent human pathogen Klebsiella pneumoniae. No evidence of preferential colonisation by hospital-adapted pathogen lineages in either vaginal or caesarean section birth groups was detected. Our analysis further enabled unbiased assessment of strain-level colonisation potential of extra-intestinal pathogenic Escherichia coli (ExPEC) in comparison with their propensity to cause bloodstream infections. Our study highlights the importance of systematic surveillance of bacterial gut pathogens, not only from disease but also from carriage state, to better inform therapies and preventive medicine in the future.

Key features of pneumococcal isolates recovered in Central and Northwestern Russia in 2011-2018 determined through whole-genome sequencing.

Invasive pneumococcal disease remains one of the leading causes of morbidity and mortality worldwide. In Russia, 13- valent pneumococcal conjugate vaccine (PCV13) was introduced into the childhood immunization programme nationwide in 2014. As part of the Global Pneumococcal Sequencing Project (GPS), we used genome data to characterize 179 pneumococcal isolates collected from Russia in 2011-2018 to investigate the circulating pneumococcal strains using a standardized genomic definition of pneumococcal lineages (global pneumococcal sequence clusters, GPSCs), prevalent serotypes and antimicrobial resistance profiles.We observed high serotype and lineage diversity among the 179 isolates recovered from cerebrospinal fluid (n=77), nasopharyngeal swabs (n=99) and other non-sterile site swabs (n=3). Overall, 60 GPSCs were identified, including 48 clonal complexes (CCs) and 14 singletons, and expressed 42 serotypes (including non-typable). Among PCV13 serotypes, 19F, 6B and 23F were the top three serotypes while 11A, 15B/C and 8 were the top three among non-PCV13 serotypes in the collection. Two lineages (GPSC6 and GPSC47) expressed both PCV13 and non-PCV13 serotypes that caused invasive disease, and were penicillin- and multidrug-resistant (MDR), highlighting their potential to adapt and continue to cause infections under vaccine and antibiotic selective pressure. PCV13 serotypes comprised 92 % (11/12) of the CSF isolates from the children aged below 5 years; however, the prevalence of PCV13 serotype isolates dropped to 53 % (31/58) among the nasopharyngeal isolates. Our analysis showed that 59 % (105/179) of the isolates were predicted to be non-susceptible to at least one class of antibiotics and 26 % (46/179) were MDR. Four MDR lineages (GPSC1, GPSC6, GPSC10 and GPSC47) accounted for 65 % (30/46) of the MDR isolates and expressed PCV13 serotypes (93 %, 28/30).This study provides evidence of high genetic and serotype diversity contributed by a mix of globally spreading and regionally circulating lineages in Russia. The observations suggest that the PCV13 vaccine could be important in reducing both invasive disease and antimicrobial resistance. We also identify potential lineages (GPSC6 and GPSC47) that may evade the vaccine.

Emergence of a multidrug-resistant and virulent Streptococcus pneumoniae lineage mediates serotype replacement after PCV13: an international whole-genome sequencing study.

BACKGROUND: Serotype 24F is one of the emerging pneumococcal serotypes after the introduction of pneumococcal conjugate vaccine (PCV). We aimed to identify lineages driving the increase of serotype 24F in France and place these findings into a global context. METHODS: Whole-genome sequencing was performed on a collection of serotype 24F pneumococci from asymptomatic colonisation (n=229) and invasive disease (n=190) isolates among individuals younger than 18 years in France, from 2003 to 2018. To provide a global context, we included an additional collection of 24F isolates in the Global Pneumococcal Sequencing (GPS) project database for analysis. A Global Pneumococcal Sequence Cluster (GPSC) and a clonal complex (CC) were assigned to each genome. Phylogenetic, evolutionary, and spatiotemporal analysis were conducted using the same 24F collection and supplemented with a global collection of genomes belonging to the lineage of interest from the GPS project database (n=25 590). FINDINGS: Serotype 24F was identified in numerous countries mainly due to the clonal spread of three lineages: GPSC10 (CC230), GPSC16 (CC156), and GPSC206 (CC7701). GPSC10 was the only multidrug-resistant lineage. GPSC10 drove the increase in 24F in France and had high invasive disease potential. The international dataset of GPSC10 (n=888) revealed that this lineage expressed 16 other serotypes, with only six included in 13-valent PCV (PCV13). All serotype 24F isolates were clustered in a single clade within the GPSC10 phylogeny and long-range transmissions were detected from Europe to other continents. Spatiotemporal analysis showed GPSC10-24F took 3-5 years to spread across France and a rapid change of serotype composition from PCV13 serotype 19A to 24F during the introduction of PCV13 was observed in neighbouring country Spain. INTERPRETATION: Our work reveals that GPSC10 alone is a challenge for serotype-based vaccine strategy. More systematic investigation to identify lineages like GPSC10 will better inform and improve next-generation preventive strategies against pneumococcal diseases. FUNDING: Bill & Melinda Gates Foundation, Wellcome Sanger Institute, and the US Centers for Disease Control and Prevention.

Changes in serotype prevalence of Streptococcus pneumoniae in Southampton, UK between 2006 and 2018.

Streptococcus pneumoniae continues to cause significant disease burden. Whilst pneumococcal conjugate vaccines (PCV) have substantially reduced this burden, serotype replacement partially negates this success due to increased disease associated with non-vaccine serotypes (NVTs). Continued surveillance is therefore essential to provide crucial epidemiological data. Annual cross-sectional surveillance of paediatric pneumococcal carriage was started in Southampton, UK following PCV7 roll-out in 2006. Nasopharyngeal swabs were collected from children < 5 years old each winter (October to March) from 2006/07 and for each consecutive year until 2017/18. Pneumococcal serotype was inferred from whole genome sequencing data. A total of 1429 (32.5%) pneumococci were isolated from 4093 children. Carriage ranged from 27.8% (95%CI 23.7-32.7) in 2008/09 to 37.9% (95%CI 32.8-43.2) in 2014/15. Analyses showed that carriage increased in children aged 24-35 months (p < 0.001) and 47-60 months (p < 0.05). Carriage of PCV serotypes decreased markedly following PCV7 and/or PCV13 introduction, apart from serotype 3 where the relative frequency was slightly lower post-PCV13 (pre-PCV13 n = 7, 1.67%; post-PCV13 n = 13, 1.27%). Prevalence of NVTs implicated in increased disease was low with 24F (n = 19, 1.4%) being the most common followed by 9N (n = 11, 0.8%), 8 (n = 7, 0.5%) and 12F (n = 3, 0.2%).

Comparative Genomics of Disease and Carriage Serotype 1 Pneumococci.

The isolation of Streptococcus pneumoniae serotypes in systemic tissues of patients with invasive disease versus the nasopharynx of healthy individuals with asymptomatic carriage varies widely. Some serotypes are hyper-invasive, particularly serotype 1, but the underlying genetics remain poorly understood due to the rarity of carriage isolates, reducing the power of comparison with invasive isolates. Here, we use a well-controlled genome-wide association study to search for genetic variation associated with invasiveness of serotype 1 pneumococci from a serotype 1 endemic setting in Africa. We found no consensus evidence that certain genomic variation is overrepresented among isolates from patients with invasive disease than asymptomatic carriage. Overall, the genomic variation explained negligible phenotypic variability, suggesting a minimal effect on the disease status. Furthermore, changes in lineage distribution were seen with lineages replacing each other over time, highlighting the importance of continued pathogen surveillance. Our findings suggest that the hyper-invasiveness is an intrinsic property of the serotype 1 strains, not specific for a "disease-associated" subpopulation disproportionately harboring unique genomic variation.

A Streptococcus pneumoniae lineage usually associated with pneumococcal conjugate vaccine (PCV) serotypes is the most common cause of serotype 35B invasive disease in South Africa, following routine use of PCV.

Pneumococcal serotype 35B is an important non-conjugate vaccine (non-PCV) serotype. Its continued emergence, post-PCV7 in the USA, was associated with expansion of a pre-existing 35B clone (clonal complex [CC] 558) along with post-PCV13 emergence of a non-35B clone previously associated with PCV serotypes (CC156). This study describes lineages circulating among 35B isolates in South Africa before and after PCV introduction. We also compared 35B isolates belonging to a predominant 35B lineage in South Africa (GPSC5), with isolates belonging to the same lineage in other parts of the world. Serotype 35B isolates that caused invasive pneumococcal disease in South Africa in 2005-2014 were characterized by whole-genome sequencing (WGS). Multi-locus sequence types and global pneumococcal sequence clusters (GPSCs) were derived from WGS data of 63 35B isolates obtained in 2005-2014. A total of 262 isolates that belong to GPSC5 (115 isolates from South Africa and 147 from other countries) that were sequenced as part of the global pneumococcal sequencing (GPS) project were included for comparison. Serotype 35B isolates from South Africa were differentiated into seven GPSCs and GPSC5 was most common (49 %, 31/63). While 35B was the most common serotype among GPSC5/CC172 isolates in South Africa during the PCV13 period (66 %, 29/44), 23F was the most common serotype during both the pre-PCV (80 %, 37/46) and PCV7 period (32 %, 8/25). Serotype 35B represented 15 % (40/262) of GPSC5 isolates within the global GPS database and 75 % (31/40) were from South Africa. The predominance of the GPSC5 lineage within non-vaccine serotype 35B, is possibly unique to South Africa and warrants further molecular surveillance of pneumococci.

Widespread sharing of pneumococcal strains in a rural African setting: proximate villages are more likely to share similar strains that are carried at multiple timepoints.

The transmission dynamics of Streptococcus pneumoniae in sub-Saharan Africa are poorly understood due to a lack of adequate epidemiological and genomic data. Here we leverage a longitudinal cohort from 21 neighbouring villages in rural Africa to study how closely related strains of S. pneumoniae are shared among infants. We analysed 1074 pneumococcal genomes isolated from 102 infants from 21 villages. Strains were designated for unique serotype and sequence-type combinations, and we arbitrarily defined strain sharing where the pairwise genetic distance between strains could be accounted for by the mean within host intra-strain diversity. We used non-parametric statistical tests to assess the role of spatial distance and prolonged carriage on strain sharing using a logistic regression model. We recorded 458 carriage episodes including 318 (69.4 %) where the carried strain was shared with at least one other infant. The odds of strain sharing varied significantly across villages (χ2=47.5, df=21, P-value <0.001). Infants in close proximity to each other were more likely to be involved in strain sharing, but we also show a considerable amount of strain sharing across longer distances. Close geographic proximity (<5 km) between shared strains was associated with a significantly lower pairwise SNP distance compared to strains shared over longer distances (P-value <0.005). Sustained carriage of a shared strain among the infants was significantly more likely to occur if they resided in villages within a 5 km radius of each other (P-value <0.005, OR 3.7). Conversely, where both infants were transiently colonized by the shared strain, they were more likely to reside in villages separated by over 15 km (P-value <0.05, OR 1.5). PCV7 serotypes were rare (13.5 %) and were significantly less likely to be shared (P-value <0.001, OR -1.07). Strain sharing was more likely to occur over short geographical distances, especially where accompanied by sustained colonization. Our results show that strain sharing is a useful proxy for studying transmission dynamics in an under-sampled population with limited genomic data. This article contains data hosted by Microreact.

A high-throughput multiplexing and selection strategy to complete bacterial genomes.

BACKGROUND: Bacterial whole-genome sequencing based on short-read technologies often results in a draft assembly formed by contiguous sequences. The introduction of long-read sequencing technologies permits those contiguous sequences to be unambiguously bridged into complete genomes. However, the elevated costs associated with long-read sequencing frequently limit the number of bacterial isolates that can be long-read sequenced. Here we evaluated the recently released 96 barcoding kit from Oxford Nanopore Technologies (ONT) to generate complete genomes on a high-throughput basis. In addition, we propose an isolate selection strategy that optimizes a representative selection of isolates for long-read sequencing considering as input large-scale bacterial collections. RESULTS: Despite an uneven distribution of long reads per barcode, near-complete chromosomal sequences (assembly contiguity = 0.89) were generated for 96 Escherichia coli isolates with associated short-read sequencing data. The assembly contiguity of the plasmid replicons was even higher (0.98), which indicated the suitability of the multiplexing strategy for studies focused on resolving plasmid sequences. We benchmarked hybrid and ONT-only assemblies and showed that the combination of ONT sequencing data with short-read sequencing data is still highly desirable (i) to perform an unbiased selection of isolates for long-read sequencing, (ii) to achieve an optimal genome accuracy and completeness, and (iii) to include small plasmids underrepresented in the ONT library. CONCLUSIONS: The proposed long-read isolate selection ensures the completion of bacterial genomes that span the genome diversity inherent in large collections of bacterial isolates. We show the potential of using this multiplexing approach to close bacterial genomes on a high-throughput basis.

Genomic surveillance of invasive Streptococcus pneumoniae isolates in the period pre-PCV10 and post-PCV10 introduction in Brazil.

In 2010, Brazil introduced the 10-valent pneumococcal conjugate vaccine (PCV10) into the national children's immunization programme. This study describes the genetic characteristics of invasive Streptococcus pneumoniae isolates before and after PCV10 introduction. A subset of 466 [pre-PCV10 (2008-2009): n=232, post-PCV10 (2012-2013): n=234;<5 years old: n=310, ≥5 years old: n=156] pneumococcal isolates, collected through national laboratory surveillance, were whole-genome sequenced (WGS) to determine serotype, pilus locus, antimicrobial resistance and genetic lineages. Following PCV10 introduction, in the <5 years age group, non-vaccine serotypes (NVT) serotype 3 and serotype 19A were the most frequent, and serotypes 12F, 8 and 9 N in the ≥5 years old group. The study identified 65 Global Pneumococcal Sequence Clusters (GPSCs): 49 (88 %) were GPSCs previously described and 16 (12 %) were Brazilian clusters. In total, 36 GPSCs (55 %) were NVT lineages, 18 (28 %) vaccine serotypes (VT) and 11 (17 %) were both VT and NVT lineages. In both sampling periods, the most frequent lineage was GPSC6 (CC156, serotypes 14/9V). In the <5 years old group, a decrease in penicillin (P=0.0123) and cotrimoxazole (P<0.0001) resistance and an increase in tetracycline (P=0.019) were observed. Penicillin nonsusceptibility was predicted in 40 % of the isolates; 127 PBP combinations were identified (51 predicted MIC≥0.125 mg l-1); cotrimoxazole (folA and/or folP alterations), macrolide (mef and/or ermB) and tetracycline (tetM, tetO or tetS/M) resistance were predicted in 63, 13 and 21.6 % of pneumococci studied, respectively. The main lineages associated with multidrug resistance in the post-PCV10 period were composed of NVT, GPSC1 (CC320, serotype 19A), and GPSC47 (ST386, serotype 6C). The study provides a baseline for future comparisons and identified important NVT lineages in the post-PCV10 period in Brazil.

Population genetic structure, serotype distribution and antibiotic resistance of Streptococcus pneumoniae causing invasive disease in children in Argentina.

Invasive disease caused by Streptococcus pneumoniae (IPD) is one of the leading causes of morbidity and mortality in young children worldwide. In Argentina, PCV13 was introduced into the childhood immunization programme nationwide in 2012 and PCV7 was available from 2000, but only in the private market. Since 1993 the National IPD Surveillance Programme, consisting of 150 hospitals, has conducted nationwide pneumococcal surveillance in Argentina in children under 6 years of age, as part of the SIREVA II-OPS network. A total of 1713 pneumococcal isolates characterized by serotype (Quellung) and antimicrobial resistance (agar dilution) to ten antibiotics, belonging to three study periods: pre-PCV7 era 1998-1999 (pre-PCV), before the introduction of PCV13 2010-2011 (PCV7) and after the introduction of PCV13 2012-2013 (PCV13), were available for inclusion. Fifty-four serotypes were identified in the entire collection and serotypes 14, 5 and 1 represented 50 % of the isolates. Resistance to penicillin was 34.9 %, cefotaxime 10.6 %, meropenem 4.9 %, cotrimoxazole 45 %, erythromycin 21.5 %, tetracycline 15.4 % and chloramphenicol 0.4 %. All the isolates were susceptible to levofloxacin, rifampin and vancomycin. Of 1713 isolates, 1061 (61.9 %) were non-susceptible to at least one antibiotic and 235(13.7 %) were multidrug resistant. A subset of 413 isolates was randomly selected and whole-genome sequenced as part of Global Pneumococcal Sequencing Project (GPS). The genome data was used to investigate the population structure of S. pneumoniae defining pneumococcal lineages using Global Pneumococcal Sequence Clusters (GPSCs), sequence types (STs) and clonal complexes (CCs), prevalent serotypes and their associated pneumococcal lineages and genomic inference of antimicrobial resistance. The collection showed a great diversity of strains. Among the 413 isolates, 73 known and 36 new STs were identified belonging to 38 CCs and 25 singletons, grouped into 52 GPSCs. Important changes were observed among vaccine types when pre-PCV and PCV13 periods were compared; a significant decrease in serotypes 14, 6B and 19F and a significant increase in 7F and 3. Among non-PCV13 types, serogroup 24 increased from 0 % in pre-PCV to 3.2 % in the PCV13 period. Our analysis showed that 66.1 % (273/413) of the isolates were predicted to be non-susceptible to at least one antibiotic and 11.9 % (49/413) were multidrug resistant. We found an agreement of 100 % when comparing the serotype determined by Quellung and WGS-based serotyping and 98.4 % of agreement in antimicrobial resistance. Continued surveillance of the pneumococcal population is needed to reveal the dynamics of pneumococcal isolates in Argentina in post-PCV13. This article contains data hosted by Microreact.

Streptococcus pneumoniae genomic datasets from an Indian population describing pre-vaccine evolutionary epidemiology using a whole genome sequencing approach.

Globally, India has a high burden of pneumococcal disease, and pneumococcal conjugate vaccine (PCV) has been rolled out in different phases across the country since May 2017 in the national infant immunization programme (NIP). To provide a baseline for assessing the impact of the vaccine on circulating pneumococci in India, genetic characterization of pneumococcal isolates detected prior to introduction of PCV would be helpful. Here we present a population genomic study of 480 Streptococcus pneumoniae isolates collected across India and from all age groups before vaccine introduction (2009-2017), including 294 isolates from pneumococcal disease and 186 collected through nasopharyngeal surveys. Population genetic structure, serotype and antimicrobial susceptibility profile were characterized and predicted from whole-genome sequencing data. Our findings revealed high levels of genetic diversity represented by 110 Global Pneumococcal Sequence Clusters (GPSCs) and 54 serotypes. Serotype 19F and GPSC1 (CC320) was the most common serotype and pneumococcal lineage, respectively. Coverage of PCV13 (Pfizer) and 10-valent Pneumosil (Serum Institute of India) serotypes in age groups of ≤2 and 3-5 years were 63-75 % and 60-69 %, respectively. Coverage of PPV23 (Merck) serotypes in age groups of ≥50 years was 62 % (98/158). Among the top five lineages causing disease, GPSC10 (CC230), which ranked second, is the only lineage that expressed both PCV13 (serotypes 3, 6A, 14, 19A and 19F) and non-PCV13 (7B, 13, 10A, 11A, 13, 15B/C, 22F, 24F) serotypes. It exhibited multidrug resistance and was the largest contributor (17 %, 18/103) of NVTs in the disease-causing population. Overall, 42 % (202/480) of isolates were penicillin-resistant (minimum inhibitory concentration ≥0.12 µg ml-1) and 45 % (217/480) were multidrug-resistant. Nine GPSCs (GPSC1, 6, 9, 10, 13, 16, 43, 91, 376) were penicillin-resistant and among them six were multidrug-resistant. Pneumococci expressing PCV13 serotypes had a higher prevalence of antibiotic resistance. Sequencing of pneumococcal genomes has significantly improved our understanding of the biology of these bacteria. This study, describing the pneumococcal disease and carriage epidemiology pre-PCV introduction, demonstrates that 60-75 % of pneumococcal serotypes in children ≤5 years are covered by PCV13 and Pneumosil. Vaccination against pneumococci is very likely to reduce antibiotic resistance. A multidrug-resistant pneumococcal lineage, GPSC10 (CC230), is a high-risk clone that could mediate serotype replacement.

The role of interspecies recombination in the evolution of antibiotic-resistant pneumococci.

Multidrug-resistant Streptococcus pneumoniae emerge through the modification of core genome loci by interspecies homologous recombinations, and acquisition of gene cassettes. Both occurred in the otherwise contrasting histories of the antibiotic-resistant S. pneumoniae lineages PMEN3 and PMEN9. A single PMEN3 clade spread globally, evading vaccine-induced immunity through frequent serotype switching, whereas locally circulating PMEN9 clades independently gained resistance. Both lineages repeatedly integrated Tn916-type and Tn1207.1-type elements, conferring tetracycline and macrolide resistance, respectively, through homologous recombination importing sequences originating in other species. A species-wide dataset found over 100 instances of such interspecific acquisitions of resistance cassettes and flanking homologous arms. Phylodynamic analysis of the most commonly sampled Tn1207.1-type insertion in PMEN9, originating from a commensal and disrupting a competence gene, suggested its expansion across Germany was driven by a high ratio of macrolide-to-ß-lactam consumption. Hence, selection from antibiotic consumption was sufficient for these atypically large recombinations to overcome species boundaries across the pneumococcal chromosome.

Phylogeography and resistome of pneumococcal meningitis in West Africa before and after vaccine introduction.

Despite contributing to the large disease burden in West Africa, little is known about the genomic epidemiology of Streptococcus pneumoniae which cause meningitis among children under 5 years old in the region. We analysed whole-genome sequencing data from 185 S. pneumoniae isolates recovered from suspected paediatric meningitis cases as part of the World Health Organization (WHO) invasive bacterial diseases surveillance from 2010 to 2016. The phylogeny was reconstructed, accessory genome similarity was computed and antimicrobial-resistance patterns were inferred from the genome data and compared to phenotypic resistance from disc diffusion. We studied the changes in the distribution of serotypes pre- and post-pneumococcal conjugate vaccine (PCV) introduction in the Central and Western sub-regions separately. The overall distribution of non-vaccine, PCV7 (4, 6B, 9V, 14, 18C, 19F and 23F) and additional PCV13 serotypes (1, 3, 5, 6A, 19A and 7F) did not change significantly before and after PCV introduction in the Central region (Fisher's test P value 0.27) despite an increase in the proportion of non-vaccine serotypes to 40 % (n=6) in the post-PCV introduction period compared to 21.9 % (n=14). In the Western sub-region, PCV13 serotypes were more dominant among isolates from The Gambia following the introduction of PCV7, 81 % (n=17), compared to the pre-PCV period in neighbouring Senegal, 51 % (n=27). The phylogeny illustrated the diversity of strains associated with paediatric meningitis in West Africa and highlighted the existence of phylogeographical clustering, with isolates from the same sub-region clustering and sharing similar accessory genome content. Antibiotic-resistance genotypes known to confer resistance to penicillin, chloramphenicol, co-trimoxazole and tetracycline were detected across all sub-regions. However, there was no discernible trend linking the presence of resistance genotypes with the vaccine introduction period or whether the strain was a vaccine or non-vaccine serotype. Resistance genotypes appeared to be conserved within selected sub-clades of the phylogenetic tree, suggesting clonal inheritance. Our data underscore the need for continued surveillance on the emergence of non-vaccine serotypes as well as chloramphenicol and penicillin resistance, as these antibiotics are likely still being used for empirical treatment in low-resource settings. This article contains data hosted by Microreact.