
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

52120
10.1038/s41467-024-52120-z
Article
Population genomics of Streptococcus mitis in UK and Ireland bloodstream infection and infective endocarditis cases
Kalizang’oma Akuzike akuzike.kalizang'oma.18@ucl.ac.uk

123
Richard Damien 4
http://orcid.org/0000-0002-1202-8540
Kwambana-Adams Brenda 1235
http://orcid.org/0000-0002-8972-7673
Coelho Juliana 6
Broughton Karen 6
Pichon Bruno 6
Hopkins Katie L. 6
Chalker Victoria 7
Beleza Sandra 8
http://orcid.org/0000-0001-8094-3751
Bentley Stephen D. 9
http://orcid.org/0000-0002-2108-1757
Chaguza Chrispin 19101112
http://orcid.org/0000-0003-4573-449X
Heyderman Robert S. r.heyderman@ucl.ac.uk

1
1 https://ror.org/02jx3x895 grid.83440.3b 0000 0001 2190 1201 NIHR Global Health Research Unit on Mucosal Pathogens, Division of Infection & Immunity, University College London, London, UK
2 Malawi Liverpool Wellcome Programme, Blantyre, Malawi
3 grid.517969.5 Department of Pathology, School of Medicine and Oral Health, Kamuzu University of Health Sciences, Blantyre, Malawi
4 https://ror.org/02jx3x895 grid.83440.3b 0000 0001 2190 1201 UCL Genetics Institute, University College London, London, UK
5 https://ror.org/03svjbs84 grid.48004.38 0000 0004 1936 9764 Department of Clinical Sciences, Liverpool School of Tropical Medicine, Liverpool, UK
6 https://ror.org/018h10037 Public Health Microbiology Division, UK Health Security Agency, Colindale, London, UK
7 https://ror.org/0227qpa16 grid.436365.1 0000 0000 8685 6563 NHS Blood and Transplant, London, UK
8 https://ror.org/04h699437 grid.9918.9 0000 0004 1936 8411 University of Leicester, Department of Genetics and Genome Biology, Leicester, UK
9 https://ror.org/05cy4wa09 grid.10306.34 0000 0004 0606 5382 Parasites and Microbes, Wellcome Sanger Institute, Hinxton, UK
10 https://ror.org/03v76x132 grid.47100.32 0000 0004 1936 8710 Department of Epidemiology of Microbial Diseases, Yale School of Public Health, Yale University, New Haven, CT USA
11 https://ror.org/03v76x132 grid.47100.32 0000 0004 1936 8710 Yale Institute for Global Health, Yale University, New Haven, CT USA
12 https://ror.org/04xs57h96 grid.10025.36 0000 0004 1936 8470 Department of Clinical Infection, Microbiology and Immunology, University of Liverpool, Liverpool, UK
6 9 2024
6 9 2024
2024
15 781229 11 2023
27 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Streptococcus mitis is a leading cause of infective endocarditis (IE). However, our understanding of the genomic epidemiology and pathogenicity of IE-associated S. mitis is hampered by low IE incidence. Here we use whole genome sequencing of 129 S. mitis bloodstream infection (BSI) isolates collected between 2001–2016 from clinically diagnosed IE cases in the UK to investigate genetic diversity, antimicrobial resistance, and pathogenicity. We show high genetic diversity of IE-associated S. mitis with virtually all isolates belonging to distinct lineages indicating no predominance of specific lineages. Additionally, we find a highly variable distribution of known pneumococcal virulence genes among the isolates, some of which are overrepresented in disease when compared to carriage strains. Our findings suggest that S. mitis in patients with clinically diagnosed IE is not primarily caused by specific hypervirulent or antimicrobial resistant lineages, highlighting the accidental pathogenic nature of S. mitis in patients with clinically diagnosed IE.

In this genomic analysis, authors observe high genetic diversity among Streptococcus mitis isolates obtained from infective endocarditis cases over 16 years in the United Kingdom and Ireland.

Subject terms

Bacterial infection
Bacterial genetics
Antibiotics
Bacterial genomics
https://doi.org/10.13039/501100000272 DH | National Institute for Health Research (NIHR) 16/136/46 NIHR200652 Kalizang’oma Akuzike issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Infective endocarditis (IE) is a life-threatening microbial infection of the interior surface lining of the heart1, and is associated with serious multi-system complications and high mortality, which approaches 25–30% at 1 year despite optimal treatment1–3. Although IE is rare, with incidence estimated between 1.5 and 11.6 cases per 100,000 people per year4, hospitalisation rates in the USA have increased by 10.6% between 2000 and 20115. In Europe, the incidence has doubled over the last two decades6. In England, incidence of IE admissions rose by 86% from 26.9 cases/million in 2009–2010 to 50.0 cases/million in 2018–197. Streptococcus mitis, is a leading cause of IE8, although it is widely considered as a typical oral commensal9.

In his Gulstonian Lectures in 1885, William Olser described micrococci within the heart valve vegetations that characterise IE10. In 1931, prior to the antibiotic era, Streptococcus pneumoniae accounted for up to ~10% of IE cases11,12, frequently associated with infection elsewhere, such as pneumonia and meningitis, and was usually fatal11,13,14. However, since the introduction of penicillin and pneumococcal conjugate vaccines, S. pneumoniae now accounts for less than 1% of IE cases11. In a multivariable analysis using S. pneumoniae as a reference, S. mitis or S. oralis has been associated with a higher IE risk with an odds ratio of 31.6 (95% CI, 19.8–50.5)15. Despite the increasing clinical importance of S. mitis, little is known about the genomic epidemiology and pathogenicity of this pathogen.

S. mitis and S. pneumoniae are closely related genetically and belong to the same species complex, yet have strikingly different pathogenic potential16. S. pneumoniae invasiveness is strongly associated with the presence of the polysaccharide capsule, pneumolysin and other virulence factors17, however, it is unclear whether pneumococcal virulence factors shared with S. mitis also facilitate S. mitis disease18,19. It is noteworthy that although largely commensal, S. mitis genomes may possess a range of pneumococcal virulence genes, including those encoding the capsule, IgA1 protease, pneumolysin, and autolysin16. The polysaccharide capsule facilitates adherence, and diffusion of molecules through to the cell surface, provides resistance to specific and non-specific host immune responses, and prevents desiccation20. S. mitis strains have been demonstrated to have capsular polysaccharide synthesis (cps) loci typically resembling that found in the pneumococcus, which sometimes results in the expression of identical capsules such as those corresponding to capsule serotypes 1, 5, and 19A21,22. It has also been suggested that S. mitis isolates harbouring the lytA gene, encoding an autolysin protein, are more likely to be associated with invasive disease as this enzyme facilitates the release of the potent pore-forming exotoxin pneumolysin and inflammatory peptidoglycan and teichoic acids from lysed bacterial cells23–25. However, it is unclear how widely these pneumococcal virulence factors are distributed among invasive S. mitis isolates. Killian and colleagues proposed that S. pneumoniae, S. pseudopneumoniae and S. mitis lineages have evolved from a pathogenic pneumococcus-like progenitor26,27, and that commensal streptococci, including S. mitis, have subsequently lost the majority of their disease-causing virulence genes as they adapt to their upper respiratory tract niche. We and others have shown that S. pneumoniae and S. mitis continue to exchange genetic material, particularly antimicrobial resistance (AMR) genes28. It remains unknown whether the presence or absence of specific virulence-associated genes and genetic backgrounds or lineages increases the pathogenicity of certain S. mitis strains.

We therefore hypothesised that S. mitis bloodstream infection (BSI) in patients with clinically diagnosed IE is associated with specific independently acquired pathogenicity loci, or lineages which have acquired particular virulence factors or AMR profiles. However, testing of this hypothesis has so far been problematic for several reasons. Firstly, the Viridans group streptococci (VGS), which includes S. mitis, are highly heterogeneous, consisting of over 50 species that are often difficult to differentiate using traditional microbiology, molecular, and biochemical techniques29, frequently leading to misidentification30. Second, and probably the most important reason, is the limited availability of publicly available S. mitis whole-genome sequencing (WGS) data in nucleotide sequence repositories. This is partly driven by the lack of focus on S. mitis compared to the other more virulent members of the Streptococcus species, for example, S. pneumoniae, which causes life-threatening infections, including pneumonia, sepsis, and meningitis25,31.

Here, we sought to expand our understanding of the genomic epidemiology of S. mitis invasive disease by undertaking WGS of a large and unique collection of 129 well-characterised IE-associated S. mitis BSI isolates collected between 2001–2016 from national bacteraemia surveillance programmes led by the British Society of Antimicrobial Chemotherapy (BSAC) and UK Health Security Agency (UKHSA). We exploited our recently developed molecular typing schemes32 and phylogenetic analysis to investigate the population structure, and the distribution of virulence and AMR genes among the IE-associated S. mitis isolates. Our genomic dataset expands the availability of S. mitis WGS data in publicly available sequence repositories by over two-fold, including the provision of nearly all the genomes from IE-associated invasive disease isolates. We show that S. mitis BSI in patients with clinically diagnosed IE is not dominated by specific hypervirulent or AMR lineages, but the pathogenicity and virulence of these strains may be enhanced by the acquisition of virulence-promoting genes through horizontal gene transfer (HGT). Our findings provide further evidence supporting the hypothesis that all S. mitis lineages may similarly cause BSIs highlighting the opportunistic nature of S. mitis infections.

Results

Characteristics of the BSAC and UKHSA IE BSI S. mitis isolates

To characterise the genetic diversity, AMR gene profiles, and identify virulence genes associated with invasive S. mitis infection, we obtained and performed WGS of 217 presumed S. mitis isolates from patients with BSI and clinically diagnosed IE from the UK and Ireland between 2001 and 2016 (Fig. 1 and Supplementary Fig. 1). These isolates were collected as part of BSAC’s Resistance Surveillance Project (n = 172), and UKHSA’s voluntary identification service (n = 45). While the retrospective nature of our analysis has preluded the rigorous application of the modified Duke/European Society of Cardiology (ESC) 2023 diagnostic criteria for IE33, the isolates were all from patients where IE had been clinically diagnosed (and given the BSI, likely fulfilling the “definite” or “possible” category), and were referred to the reference laboratory for species confirmation and further antibiotic sensitivity testing. Two of 172 BSAC and three of 45 UKHSA isolates were not viable after bacterial culture of transport swabs and were, therefore, not sequenced. One UKHSA isolate subsequently failed WGS quality control due to low DNA. Accurate species determination among the VGS has been a challenge using conventional and molecular approaches34, therefore, we used a bioinformatic approach for species confirmation among the presumed S. mitis isolates obtained from BSAC and UKHSA.Fig. 1 Schematic of the study design and analysis workflow.

The S. mitis BSI isolates analysed were collected from clinically diagnosed IE cases, between 2001 and 2016, submitted to BSAC and UKHSA. The analysis involved WGS, population structure analysis, AMR genotyping and phenotyping, and identification of virulence genes using genotyping and bacterial genome-wide association analysis approaches to identify hypervirulent and dominant S. mitis lineages. The map of the UK and Ireland was generated by the authors in R software using the maps v4.0.0 package (https://cran.r-project.org/web/packages/maps/). The confirmed S. mitis isolates are described in Supplementary Data 1 and 2.

Overall, we confirmed S. mitis by WGS in 106 of 170 (62.4%) and 23 of 41 (56.1%) of the viable BSAC and UKHSA isolates, respectively (Fig. 1, Supplementary Data 1 and 2). We found that S. mitis was most frequently misidentified as S. oralis and S. infantis for 29 of 211 (13.7%) and 27 of 211 (12.8%) sequenced isolates, respectively (Supplementary Table 1). The average nucleotide identity (ANI) is a measure of nucleotide level similarity among orthologous genes shared between two genomes and it offers robust resolution between strains of the same or closely related species35. Pairwise ANI values calculated among the 129 S. mitis isolates ranged from 91.2 to 99.6% (Median = 93.4%) (Supplementary Fig. 2). Some of these ANI values fell below the 94–96% range generally been accepted to demarcate species boundaries36,37. We speculated that such low ANI values of up to 91% observed in S. mitis genomes reflected the fact that these strains form a species complex consisting of a continuum of related but genetically diverse lineages sufficient to be considered separate sister Streptococcal species38. Additionally, we constructed a phylogenetic tree of the 129 S. mitis isolates, in the context of 188 global S. mitis isolates obtained from public nucleotide sequence repositories (Supplementary Data 3) and 5 reference genomes (Supplementary Data 4), to visualise and compare the multiple speciation approaches (Supplementary Fig. 3). Our approach of taxonomic assignment based on several approaches (see methods) resolved known species inconsistencies resulting from less discriminatory phenotypic approaches and the use of different databases when using genotypic data. Altogether, this highlights the importance of WGS and complementary genotypic-based approaches for species confirmation among the VGS34. The S. mitis isolates from patients with IE were from patients of all ages (0–99 years), with the age range 50–59 years having the highest frequency of S. mitis isolates (20.2%; 26/129) (Supplementary Table 2). Of the 129 confirmed S. mitis isolates, 70 (54.3%) were collected from men and 58 (45.0%) from women (Supplementary Table 2).

Population genetic diversity of IE-associated S. mitis

S. mitis strains are known to be highly diverse genetically16,26,32,39–42, but accurately assessing S. mitis population structure has been a challenge due to limited S. mitis genome sequences and species-specific molecular genotyping tools. Because of this, it remains unclear whether specific lineages predominantly cause BSIs associated with S. mitis in localised populations and over time. To address this, we therefore quantified the genetic diversity of the IE-associated S. mitis isolates collected from the UK and Ireland, and we determined if population-level genetic diversity changed over time during the 16-year surveillance period. We calculated the ANI values and number of single nucleotide polymorphisms (SNPs) between all the pairs of isolates to quantify the genetic diversity. The number of non-ambiguous SNPs between the pairs of isolates ranged from ~45,000 to 55,000 bp out of a total of 2,146,613 nucleotide bases from mapped sequencing reads (Supplementary Fig. 4) and varied significantly over time (Kruskal–Wallis test, p < 0.001) across the 16-year period. Similarly, the ANI values ranged from 91.8–97.1% (Supplementary Fig. 5) and varied significantly over time (Kruskal–Wallis test, p < 0.001). The observed high number of SNPs distinguishing pairs of isolates sampled in the same year and the wide range of ANI values suggested a high genetic diversity of the IE-associated S. mitis strains during the surveillance period.

Next, we confirmed the high genetic diversity of the IE-associated S. mitis isolates by the identification of long internal and terminal branches separating the isolates in the constructed maximum-likelihood phylogenetic tree (Fig. 2). The long phylogenetic branches indicated the existence of several distinct lineages. To confirm the presence of multiple distinct lineages, we then used our recently developed S. mitis MLST scheme available on the PubMLST website (https://pubmlst.org/smitis) and a complementary whole-genome-based sequence clustering approach using the PopPUNK framework (https://www.bacpop.org/poppunk/)43. This analysis identified 127 (98.4%) unique sequence types (STs) based on the S. mitis MLST scheme and an equal number of the PopPUNK lineages, that we defined as Global Sequence Clusters (GSCs), among the 129 IE-associated S. mitis isolates (Fig. 2). Two isolates isolated in 2014 and 2016, which belonged to a single ST (ST30) and lineage (GSC27), differed from each other by 3024 SNPs, had a pairwise ANI value of 99.2%, and the isolated pair had the same penicillin susceptibility profile. Similarly, two ST36 isolates collected in 2007 and 2015 belonged to lineage GSC28 and differed from each other by 6413 SNPs, had a pairwise ANI value of 99.6%, but differed by their penicillin susceptibility profiles (Fig. 2). Since we did not have access to patient-identifiable data from BSAC or UKHSA, we could not exclude the possibility that S. mitis isolates of the same STs represented recurrent infection of the same patient. Nonetheless, considering that 127 out of 129 isolates (98.4%) belonged to different STs and lineages, and the large range in pairwise distances (3024–62,803 SNPs) and ANI values (91.8–97.1%), these observations suggest that BSI and IE-associated S. mitis in the UK and Ireland is not predominantly caused by a select few dominant lineages.Fig. 2 Genotypic and phenotypic characteristics of the Streptococcus mitis IE isolates.

Maximum-likelihood phylogenetic tree of 129 genotypically-confirmed IE S. mitis isolates. The phylogeny is built using 281,737 SNPs out of a total of 2,146,613 nucleotide bases from mapped sequencing reads, and is displayed next to metadata that includes GSC, ST, year of isolation, phenotypic susceptibility, and the presence or absence of acquired antibiotic resistance genes. Source data are provided as a Source Data file. The strain name in the same order as the phylogeny is shown in the source data file. The phylogenetic tree shows high genetic diversity of the IE-associated BSI S. mitis isolates, such that 98.4% of the isolates belong to different STs and GSCs. A limited number of resistance genes were explored as the analysis was limited to genes present in the curated public AMR databases89. Unlike S. pneumoniae, no genomic tools and public databases are available for exploring resistance in S. mitis, to antibiotics such as beta-lactam antibiotics, that are determined through point mutations89,92. Phenotypic and genotypic non-susceptibility to multiple antibiotic classes was observed for several isolates throughout the surveillance period and across the phylogeny. There was good phenotypic to genotypic resistance concordance, as 96.9% (31/32) of phenotypically tetracycline-resistant isolates had a tet(M) gene, and 94.9% (74/78) of phenotypically macrolide-resistant isolates had a mef(A) or erm(B) gene. Source data are provided as a Source Data file.

Temporal trends in AMR among the IE-associated S. mitis isolates

Due to the limited focus, misdiagnosis, and low incidence of S. mitis BSIs, temporal AMR trends for S. mitis associated with IE have not been well described. We, therefore, assessed the phenotypic susceptibility of the IE-associated S. mitis isolates against commonly used antibiotics to treat suspected S. mitis IE (penicillin, amoxicillin, gentamicin, and vancomycin)44. We found the distribution of penicillin and amoxicillin non-susceptible isolates across the entire phylogeny, not restricted to only specific phylogenetic branches containing closely related isolates (Fig. 2). Among isolates with phenotypic MIC data, 23.3% (30/129) and 6.2% (6/97) of the isolates were non-susceptible to penicillin and amoxicillin, respectively (Supplementary Data 1 and 2, antibiotic abbreviations used by BSAC and UKHSA are explained in Supplementary Table 3). We show that all 30 penicillin non-susceptible isolates belonged to different STs and GSCs (Fig. 2). We observed non-susceptibility to penicillin among isolates across the surveillance period (Fig. 2), which would impact the use of penicillin as a first-line antibiotic for Streptococcal IE44. All the isolates had low-level gentamicin resistance (MIC ≤ 128), which would not impact its use as a synergistic antibiotic in IE management44. Additionally, all isolates showed full susceptibility to vancomycin. Resistance to erythromycin and tetracycline, antibiotics not used for S. mitis IE treatment33, was observed throughout the surveillance period. S. mitis is a known reservoir of resistance genes45, and carriage of macrolide resistance isolates among 51% of healthy individuals, with a significant correlation with tetracycline co-resistance, has been described45. Macrolide and tetracycline resistance genes are often co-carried on mobile genetic elements46,47, self-transmissible DNA sequences that can move between and within species. Therefore, S. mitis, regardless of isolation source may be a reservoir of transmissible AMR for other Streptococcus species. Our AMR findings support the use of the current antibiotic treatment regimens for the management of IE; however, our data emphasises that continued surveillance remains critical for monitoring AMR trends, particularly for penicillin.

To further assess the temporal changes in antimicrobial susceptibility, we aggregated phenotypic MICs for several antibiotics into four 3-year intervals, namely 2001–2004, 2005–2008, 2009–2012, and 2013–2016. Due to the small number of phenotyped isolates per year, these 3-year intervals ensured the derivation of more robust estimates for the phenotypic MIC trends based on a sufficient number of isolates. Using this approach, we found no statistically significant differences in the median MICs across all four time intervals for seven out of eight antibiotics (Table 1 and Supplementary Fig. 6). Conversely, the median MICs for gentamicin showed a statistically significant decrease over time (Kruskal–Wallis test, p = 0.001). Despite the variability of the MIC changes over time due to the limited number of phenotyped isolates, our findings provide baseline data for the genomic surveillance of AMR in S. mitis-associated BSIs, including IE, in the UK and Ireland, regionally and globally.Table 1 Antibiotic median MIC for Streptococcus mitis isolates from patients with IE across four 3-year intervals

	Year	Number of isolates*	Median MIC (μg/mL)	Number of non-susceptible isolates and proportion (%)	P value**	
Amoxicillin	2001–2004	34	0.03	2 (5.9)		
n = 97	2005–2008	21	0.06	4 (19.0)		
	2009–2012	25	0.03	2 (8.0)		
	2013–2016	17	0.06	2 (11.8)	0.06	
Cefotaxime	2001–2004	38	0.06	2 (5.3)		
n = 118	2005–2008	23	0.06	1 (4.3)		
	2009–2012	29	0.06	3 (10.3)		
	2013–2016	28	0.125	0 (0)	0.403	
Gentamicin	2001–2004	38	6	0 (0)		
n = 118	2005–2008	23	4	0 (0)		
	2009–2012	29	4	0 (0)		
	2013–2016	28	2	0 (0)	0.001	
Penicillin	2001–2004	41	0.03	4 (9.8)		
n = 129	2005–2008	28	0.064	10 (35.7)		
	2009–2012	29	0.03	6 (20.7)		
	2013–2016	31	0.06	10 (32.3)	0.111	
Vancomycin	2001–2004	38	0.5	0 (0)		
n = 118	2005–2008	23	0.5	0 (0)		
	2009–2012	29	0.5	0 (0)		
	2013–2016	28	0.5	0 (0)	0.230	
*Phenotypic MIC data for each isolate is described in Supplementary Data 1 and 2.

**The P value was calculated using the Kruskal–Wallis test.

Other Streptococcus species, including S. pneumoniae, S. agalactiae or Group B Streptococcus (GBS), and S. pyogenes or Group A Streptococcus (GAS), are characterised by lineages that are more likely to cause invasive diseases in humans48–51. The identification and tracking of lineages have facilitated genomic surveillance and guided clinical interventions against these species52, an approach that could be equally valuable for monitoring the epidemiology of invasive S. mitis strains. We, therefore, analysed the 129 S. mitis isolates from patients with clinically diagnosed IE in the context of globally sampled strains to better understand the global genetic diversity and the distribution of AMR and virulence genes amongst invasive S. mitis. We compiled a total of 322 confirmed whole-genome sequenced S. mitis isolates, from the present study and publicly available genomic sequence repositories, representing 258 PopPUNK lineages and 259 STs (Fig. 3a, b). Analysis of the metadata for our sequenced isolates and the contextual publicly available sequences revealed that 158 out of 322 isolates (49.1%) were from carriage, 152 (47.20%) were from invasive disease, and 12 (3.7%) were from unknown sources (Supplementary Data 3). Of the invasive isolates, 138 (42.9%) were from patients with IE (129 of these were from this study), 13 (4.0%) were from bacteraemia, and 1 (0.3%) was from pneumonia. Overall, there were no shared lineages between the S. mitis isolates obtained from patients in the UK and Ireland with clinically diagnosed IE and other global strains from carriage or invasive disease (Fig. 3c). However, the global strains that shared STs and GSCs were part of a previous carriage study that sampled the same individuals, such that the same strain was sampled multiple times32. Therefore, S. mitis isolates from asymptomatic carriage and invasive disease are distributed across the entire phylogeny of the global isolates, indicating the potential for all, rather than a select few lineages, to cause IE. Furthermore, we found no clustering of isolates based on the AMR genes (Fig. 3c), or major virulence genes associated with pneumococcal pathogenicity. Two of the virulence genes, encoded within the cps locus region and pneumolysin (ply), were found in distinct positions on the global S. mitis phylogeny (Fig. 3c). Pneumococcal adherence and virulence protein A (pavA) and pneumococcal surface adhesin A (psaA) genes were present among all 322 isolates (100%) (Supplementary Data 5). Together, these findings demonstrate that even when viewed from the global context, invasive S. mitis strains are not predominantly associated with a single or limited number of lineages.Fig. 3 Population structure, virulence, and antimicrobial resistance gene profiles of IE Streptococcus mitis in context of global isolates.

a Frequency plot of the STs identified across the combined IE and global S. mitis dataset. The plot shows that the UK IE S. mitis belonged to unique STs, however, two isolates isolated in 2014 and 2016 belonged to ST30, and another two isolates collected in 2007 and 2015 were assigned to ST36. Multiple carriage isolates belonging to the same ST were largely from a previous study that sampled multiple isolates from individuals32. ST190 had the highest frequency of 23 isolates and was also likely due to sampling from the same individual. b Frequency plot of the GSCs identified across the combined IE and global S. mitis dataset. The plot shows that the UK IE S. mitis belonged to unique GSCs, however, two isolates isolated in 2014 and 2016 belonged to GSC27, and another two isolates collected in 2007 and 2015 were assigned to lineage GSC28. Multiple carriage isolates belonging to the same GSC were largely from a previous study that sampled multiple isolates from individuals32. GSC1 had the highest frequency of 23 isolates and was also likely due to sampling from the same individual. c Maximum-likelihood phylogeny of IE and global S. mitis is built using 473,175 SNPs out of 1,237,113 core nucleotide bases. The coloured tips of the phylogeny and the first horizontal metadata bar show the isolation condition of the S. mitis isolates. From the second to fifteenth horizontal bars, the isolate metadata shows the IE UK S. mitis isolates, the GSC lineage, and virulence genes. The virulence gene matrix has 7 capsule genes (cps4A – cps4F), autolysins (lytA and lytC), and pneumolysin (ply) genes. The phylogeny shows clustering of capsule genes among isolates in one region of the phylogeny, predominantly UK IE isolates. The last four horizontal bars show antimicrobial resistance (AMR) gene matrices. From the first to fourth matrix bar are the presence or absence of chloramphenicol (cat), macrolide (ermB and mefA), and tetracycline (tetM) genes. Source data are provided as a Source Data file.

Differential abundance of putative pathogenicity enhancing genes

Previous studies have suggested that HGT between S. mitis and other more virulent members of the Streptococcus genus, specifically S. pneumoniae, drives the spread of AMR and virulence between these species26,32,53. Additionally, S. mitis is known to harbour pneumococcal virulence genes, including those involved in the biosynthesis of serotype 1 and 5 capsules21,22 and homologues of other pathogenicity-associated genes, including Zinc metalloproteases (zmpC, zmpC, and zmpD)54, neuraminidases (nanA and nanB), pneumolysin (ply), immunoglobulin A protease (iga), and autolysins (lytA-C)18 and glucan binding protein B (gbpB or pcsB)55. However, no systematic analyses to assess the abundance of all the functionally characterised (or known) and hypothetical genes in the S. mitis pan-genome between the invasive and carriage isolates have been conducted to date, mostly due to the limited availability of genomic data. Since S. mitis is widely regarded as a source of virulence factors which enhance the pathogenicity of pneumococcal strains26, we, therefore, speculated that the S. mitis strains associated with BSIs may show a higher abundance of virulence genes compared to the isolates sampled from the asymptomatic carriage. To address this, we employed a bacterial genome-wide association study (GWAS)-type approach (Fig. 4), increasingly used to identify genomic loci associated with bacterial phenotypes56–60. Our null hypothesis was that no gene influenced the pathogenicity of S. mitis. Therefore, we expected the distribution of any gene would be similar among the BSI and carriage isolates due to the inclusion of phylogenetically similar but phenotypically distinct pairs of isolates. However, we found fifteen orthologous gene clusters, whose identifiers were arbitrarily defined with the prefix “SCLS”, for sequence cluster locus sequence, were differentially overrepresented among either IE-associated or carriage isolates (Table 2 and Supplementary Fig. 7). Among these genes were TP-binding cassette (ABC) transporters, competence-specific and Hca operon transcription regulators, phage-associated proteins, autolysin (a known pneumococcal virulence factor25,61–63), and several uncharacterised hypothetical proteins. Twelve of the genes were overrepresented in the carriage isolates when compared to invasive disease isolates, while three genes showed the opposite association.Fig. 4 Bacterial genome-wide association analysis.

The phylogeny of global confirmed S. mitis whole-genome sequences was built using 473,175 SNPs out of 1,237,113 nucleotide bases and annotated with disease status that was pruned to select for pairs of genetically closest carriage and invasive disease isolates. Source data are provided as a Source Data file. This phylogenetic-based approach provided an approximate matching of the isolates for the bacterial genome-wide association analysis. Source data are provided as a Source Data file.

Table 2 Summary of the S. mitis genes, which were differentially abundant among a subset of phylogenetically paired invasive disease and carriage isolates

Gene cluster*	Gene name	Gene presence in the paired isolates	P value**	Gene product/description	
None	Carriage	Disease	Both	
SCLS1	btuD	17	12	1	14	0.0055	ABC transporter, ATP-binding protein	
SCLS2		17	12	1	14	0.0055	ABC-2 family transporter protein	
SCLS3		31	1	9	3	0.0269	Transcriptional regulator ComX2	
SCLS4		29	9	1	5	0.0269	DNA-binding phage protein	
SCLS5		9	6	17	12	0.0371	Hypothetical protein	
SCLS6		19	10	2	13	0.0433	ComC/BlpC family leader-containing pheromone/bacteriocin***	
SCLS7		20	2	10	12	0.0433	Hypothetical protein	
SCLS8	lytA	31	10	2	1	0.0433	Autolysin	
SCLS9		34	8	1	1	0.0455	SPFH domain-containing protein***	
SCLS10		35	8	1	0	0.0455	Major Facilitator Superfamily (MFS) transporter***	
SCLS11		34	8	1	1	0.0455	Phage protein	
SCLS12	hcaR	35	8	1	0	0.0455	Hca operon transcriptional activator HcaR	
SCLS13		35	8	1	0	0.0455	YbhB/YbcL family Raf kinase inhibitor-like protein***	
SCLS14		34	8	1	1	0.0455	Hypothetical protein	
SCLS15		34	8	1	1	0.0455	Phage transcriptional regulator, Cro/CI family protein	
*We arbitrarily defined the orthologous gene clusters with the prefix “SCLS”, which stands for the sequence cluster locus sequence. Specific nucleotide sequences of the representative genes in each orthologous gene cluster inferred from pan-genome clustering analysis using Panaroo (see methods) are provided as a Source Data file.

**The two-sided P value was calculated using McNemar’s exact test based on phylogenetically paired invasive disease and carriage isolates.

***Gene description determined through the online NCBI BLAST tool, its databases, and using default parameters.

As the pneumococcus is a close relative of S. mitis, belonging to the same species complex, we screened 493 invasive pneumococcal isolates obtained from blood from patients with bacteraemia for the presence of these 15 overrepresented genes in S. mitis (Fig. 5). The pneumococcal isolates were obtained from the Global Pneumococcal Sequencing Project and were part of the Centers for Disease Control and Prevention’s (CDC) active bacterial core surveillance (Supplementary Data 6). We identified the presence of 6 out of 15 overrepresented orthologous S. mitis gene clusters that were also prevalent among the pneumococcal isolates. Among these gene clusters, the pneumococci had a high prevalence of 2 out of 3 genes that were more overrepresented in the invasive S. mitis when compared carriage isolates, a transcriptional regulator (gene cluster SCLS3) and hypothetical gene (gene cluster SCLS7). However, gene cluster SCLS5, a hypothetical gene, was absent in all the pneumococcal strains. In future work, targeted mutagenesis experiments would aid in understanding the function of these hypothetical genes in S. mitis as well as pneumococcus. Together, these findings suggest that these identified genes may potentially modulate the pathogenicity of S. mitis, facilitating a potential rare transition from a typical commensal to a pathogenic lifestyle64,65.Fig. 5 Phylogenetic distribution of orthologous gene clusters differentially overrepresented between IE and non-IE S. mitis isolates among invasive isolates belonging to the genetically related and more pathogenic sister species pneumococcus.

Maximum-likelihood phylogeny of invasive bacteraemia S. pneumoniae isolates is built using 141,880 SNPs out of 1,520,986 core nucleotide bases. The coloured tips of the phylogeny and the first metadata row shows the Global Pneumococcal Sequence Cluster (GPSC), while the second row shows the serotype. The subsequent 12 rows show the presence or absence of the orthologous gene clusters (with prefix SCLS for sequence cluster locus sequence) that were overrepresented in invasive or carriage S. mitis. The phylogeny shows there are a few gene clusters that are also prevalent among the invasive pneumococcal isolates, however, most gene clusters are absent among the pneumococcal isolates. The largely absent S. mitis lytA gene was further investigated among the pneumococci using the virulence finder database90, which identified the presence of pneumococcal lytA among 448 out of 493 (90.9%) isolates using pneumococcal specific lytA reference sequences in the virulence finder database90. Source data are provided as a Source Data file.

Discussion

Feared by clinicians for the potential for a missed diagnosis, IE has been a puzzle in medical science, including diagnostic microbiology, since before the time of William Osler10. In this study, we undertook WGS of a large and unique collection of clinically diagnosed IE-associated S. mitis isolates, expanding the number of publicly available sequenced genomes for this poorly studied opportunistic pathogen by nearly two-fold. Taking into consideration the rarity of IE4, this dataset represents one of the largest IE-associated S. mitis datasets to date. Contrary to our hypothesis, we did not find dominant hypervirulent lineages or populations characterised by the presence of unique virulence genes or AMR across the surveillance period. These data suggest that all S. mitis isolates have the potential to cause IE-associated BSI, further supporting the notion that S. mitis is likely an accidental pathogen. Our analysis suggests that the presence of pneumococcal virulence genes and phenotypic AMR has not led to a selective advantage in the carriage population whereby increased risk for causing disease is a consequence. However, our phylogeny-based GWAS suggest that there are genes that may potentially modulate the pathogenicity of S. mitis, facilitating the rare switch from a typical commensal to a pathogenic lifestyle. The identification of several S. mitis populations with different virulence profiles similar to S. pneumoniae supports the hypothesis that the commensalism of S. mitis evolved through genome reduction and loss of virulence genes to attain a commensal lifestyle16,26. Additionally, our chronologically sampled S. mitis dataset provides the foundation on which to monitor changes in the population structure and AMR of IE-associated S. mitis.

We have demonstrated that WGS and a combination of appropriate bioinformatic analytical techniques are crucial for the accurate identification of S. mitis species from related VGS strains, which have implications for understanding S. mitis disease and epidemiology. Here we show that WGS more accurately distinguishes S. mitis isolates from closely related species compared to phenotypic testing and matrix-assisted laser desorption ionization–time of flight mass spectrometry (MALDI-TOF) methods. Although MALDI-TOF is increasingly used by reference laboratories to rapidly identify bacterial species66, the method fails to accurately differentiate species of the Mitis group67. WGS, therefore, is a potential adjunct for accurately differentiating S. mitis and other VGS, which avoids false diagnosis at the patient level and incorrect ascertainment of the contribution of S. mitis to IE relative to related VGS species. Indeed, the false species ascertainment at diagnostic laboratories may have contributed to the low numbers of S. mitis isolates obtained across the 16-year period, and it is likely that other VGS-causing IE may be S. mitis. However, WGS is not routinely carried out in diagnostic laboratories, and the option for sequencing all archived VGS from UKHSA and BSAC laboratories to increase the S. mitis sample size in this study was resource-intensive and not typically feasible.

S. mitis that have been isolated from patients with clinically diagnosed IE between 2001 and 2016 in the UK and Ireland have remained susceptible to antibiotics commonly used as first-line treatment44. We did not identify high-level vancomycin or gentamicin resistance among the IE-associated isolates, suggesting that gentamicin used synergistically with penicillin or vancomycin is also likely to remain effective against S. mitis. Penicillin MICs have remained stable throughout the surveillance period, which is consistent with data from the USA between 2010 and 202068. However, as penicillin non-susceptible isolates were identified across the surveillance period, continued surveillance to monitor the AMR trends among IE-associated S. mitis isolates remains important alongside WGS, as species misidentification can skew AMR trends and distribution as species-specific differences in AMR among the VGS have been identified69.

We have shown extensive genetic diversity among IE S. mitis, marked by the association of each isolate with a unique lineage and a variable distribution of pneumococcal virulence genes. Similar to models proposed for Staphylococcus epidermidis opportunistic infection70, our findings suggest that S. mitis is a true commensal with accidental pathogenicity, such that multiple genetically divergent clones are found to cause IE-associated BSI and have virulence determinants equally distributed among isolates from carriage and disease. Therefore, as previously highlighted elsewhere19, the presence of pneumococcal virulence genes alone may not necessarily be enough to determine S. mitis invasiveness. For example, the polysaccharide capsule is a major virulence factor for the pneumococcus20, however, transformation experiments have shown that serotype 4 capsule acquisition by S. mitis does not increase resistance against early clearance in a mouse model to levels similar to S. pneumoniae serotype 4 capsule wild-type19. It is possible that a more complex combination of genes related to virulence, metabolism, and other functions, in addition to the host immunity, may determine the potential to cause IE.

PavA and PsaA adhesion genes are known to play a role in the pathogenesis of pneumococcal disease20,71, and the genes may also potentially contribute to the adhesion properties among S. mitis8,18. We established that both genes were not exclusive to the IE S. mitis isolates, but PavA and PsaA were present among all carriage and disease S. mitis isolates. It remains unclear whether there are differences in the presence of other adhesion genes among carriage and disease isolates. Targeted mutagenesis is potentially a valuable approach in identifying endovascular IE virulence factors, but there is considerable redundancy in the streptococcal virulence factors that may mediate disease. In the context of such extensive genome variation, a bacterial GWAS-type approach offers the opportunity to identify overrepresented genes associated with disease isolates72. Here we show, using a pilot phylogenetic-based GWAS-type approach, that some genes may be overrepresented in disease when compared to carriage isolates, which may potentially influence the pathogenicity of S. mitis and include a transcription regulator and uncharacterised hypothetical proteins. Although larger genome datasets are required to validate these pilot GWAS analyses, this will require the systematic collection of many hundreds to thousands of S. mitis isolates from individuals with asymptomatic carriage and disease to achieve sufficient statistical power to unravel potential associations.

Our finding of multiple populations with distinct virulence profiles is concordant with the current theory that S. mitis has transitioned from a more pathogenic species to adopt a commensal lifestyle. S. mitis and S. pneumoniae are very closely related yet have strikingly different pathogenic potentials16. Both species have been suggested to have evolved from a common ancestor with all properties associated with virulence16, however, S. mitis evolved through reductive evolution to become 15% smaller in genome size compared to the pneumococcus and adopt a more commensal lifestyle27. The mosaic pattern of virulence gene presence or absence among S. mitis lineages supports the suggestion that gene loss has occurred and may still be ongoing in the species16. Although this theory has been described as linear, it is possible that multiple selective pressures at various time points may have facilitated the emergence of several populations of S. mitis as we observed in our analysis.

S. mitis is known to be highly transformable, therefore, can acquire genetic material through HGT73. It is, therefore, also possible that S. mitis evolution through both gene loss and gain, via HGT and homologous recombination processes continuously shapes the genetic diversity of this species. We, therefore, hypothesise that S. mitis lineages with increased invasiveness could potentially emerge through acquiring a combination of virulence genes and other determinants from closely related pathogenic species such as the pneumococcus. Therefore, continued surveillance of IE-associated S. mitis infections is important for monitoring the changing population structure and virulence of S. mitis lineages.

Our study has some limitations. First, although the present study relies on 129 newly sequenced genomes, which equates to 40.1% (129/322) of the total publicly available S. mitis genomes, our sample size is still limited. Although the study was restricted to one geographical region, we utilised a valuable resource of archived IE-associated isolates with phenotypic data obtained through robust surveillance systems not readily found elsewhere. A prospective study would be advantageous, however, due to the rarity of S. mitis-associated invasive disease, such a study would be possible but challenging as it would take many years to achieve even a modest sample size. While the sample size requirement could be partly resolved by conducting a multi-site study involving several countries, such a study is likely to be costly and difficult to justify considering the overall low incidence of S. mitis invasive disease compared to diseases caused by other bacterial pathogens. In contrast, retrospective studies of already collected isolates from several countries through international collaborations, as evidenced by consortiums such as the Global Pneumococcal Sequencing (GPS) project74, may provide additional insights on S. mitis diversity and pathogenicity. However, considering the low incidence of S. mitis-associated IE, our dataset represents a unique and the largest collection of IE-associated S. mitis isolates to date, which will provide much-needed baseline genomic data for further comparative studies of S. mitis diversity and pathogenicity. Secondly, as highlighted, we were not able to retrospectively ascertain whether these S. mitis BSI in patients with clinically diagnosed IE fulfilled the modified Duke/ESC 2023 diagnostic criteria for IE75. While this potential imprecision may have affected our ability to identify rare hypervirulent lineages, given the BSI it is likely that these patients fulfilled the “definite” or “possible” IE categories, it is very unlikely that misclassification materially biased the conclusions of our genomic analysis. Thirdly, the S. mitis isolates were not collected systematically from all regions of the UK and Ireland as submission of isolates for bacterial surveillance was voluntary, possibly introducing bias in the samples submitted by the hospital laboratories to BSAC and UKHSA. Lastly, due to the retrospective nature of the study, we did not have access to S. mitis carriage isolates from the same geographical region and across a similar time frame, which would have helped to contextualise the IE-associated S. mitis isolates.

In conclusion, using a rare collection of IE-associated S. mitis isolates from the UK and Ireland, we have shown that S. mitis isolates from patients with suspected IE are highly diverse, with a wide distribution of AMR genes. We have shown that suspected S. mitis associated IE disease is not predominantly caused by a select few dominant lineages and that the presence of known virulence genes from S. pneumoniae does not noticeably influence invasiveness. However, our pilot GWAS-type approach of phylogenetically paired invasive and carriage S. mitis isolates suggest that some genes may be differentially abundant among invasive and carriage S. mitis strains and thus may likely modulate pathogenicity. These findings, therefore, provide further evidence for opportunistic and accidental pathogenicity and expand on the existing theory of the commensal lifestyle of S. mitis. While our AMR findings support the use of the current antibiotic treatment regimens for the management of IE, continued surveillance remains critical for monitoring the AMR trends particularly for penicillin. Our unique dataset expands the publicly available genomic dataset S. mitis by over two-fold and constitutes nearly all invasive IE-associated S. mitis genomic data available to date, which will therefore provide critical baseline data to inform further in-depth investigations of the epidemiology and biology of IE and other BSIs caused by this accidental pathogen.

Methods

Ethical approval

This study complies with ethical regulations applied in public health surveillance. UKHSA holds approvals to process patient-identifiable data for the purposes of infectious disease surveillance, in accordance with Section 60 of the Health and Social Care Act 2001. All isolates were anonymised to the key researcher, and patient-identifiable information was not included in the study. Isolates submitted to BSAC were collected as part of routine clinical investigations and were processed as such by the original laboratory.

Bacteraemia surveillance isolates and sample selection

The isolates included in the analysis were from BSI surveillance. The first batch of isolates were obtained from the British Society of Antimicrobial Chemotherapy (BSAC) Resistance Surveillance Project, which is a long-term study that aims to monitor AMR among bacterial isolates from lower respiratory tract infections and from BSI among patients in the UK and Ireland (https://bsac.org.uk/)76. The bacterial isolates are re-identified by the BSAC Central Testing Laboratory using matrix-assisted laser desorption ionisation-time of flight mass spectrometry (MALDI-TOF MS) to confirm the species. The second batch of isolates were obtained from the UKHSA voluntary surveillance, which includes streptococcal isolates obtained from patients with BSI. Species identification was confirmed by UKHSA using phenotypic analytical profile index (API) testing.

All available isolates that were collected, archived, and identified as S. mitis from suspected IE cases by BSAC and UKHSA from 2001–2016 were included in the study. However, well-phenotyped isolates from more recent years were not available. The diagnosis of IE assigned to the UKHSA and BSAC S. mitis isolates was made by the referring clinical teams. The modified Duke/ESC 2023 diagnostic criteria for IE33 were not available. However, in view of the BSI, and the referral of the isolates for species confirmation and antibiotic sensitivity testing for the management of IE, the patients likely fulfilled the “definite” or “possible” modified Duke/ESC 2023 diagnostic categories. The geographical locations of individual isolates were not available, however, BSAC surveillance covered the UK and Ireland76, while UKHSA surveillance covered England, Wales, and Northern Ireland. Age ranges for the patients were available from both surveillance programmes.

Bacterial culture and antimicrobial susceptibility testing

Bacterial transport swabs obtained from BSAC and UKHSA were inoculated on Columbia agar plates supplemented with 5% horse blood (CBA) (Thermo Scientific, UK). CBA plates were incubated at 37 °C in 5% CO2 for 18 hr, then a single colony was sub-cultured to obtain a plate with pure growth. Pure presumed S. mitis colonies were picked and suspended in cryovials with 1 ml of Todd Hewitt Broth with yeast extract (THY) (Merck, Germany) and 20% glycerol (Merck, Germany), then stored at −80 °C for downstream processing.

Amoxicillin, cefotaxime, clindamycin, erythromycin, gentamicin, penicillin, tetracycline, and vancomycin MICs were previously determined for the presumed S. mitis isolates by the reference laboratory using the agar dilution method. Where data was not available, we derived penicillin MICs for presumed S. mitis isolates, the first-line antibiotic option for S. mitis IE44,75, using the E-test® (bioMérieux, UK). The American Type Culture Collection (ATCC) 49619 S. pneumoniae strain was used as an internal control for all Antimicrobial Susceptibility Testing (AST). Although not consistently determined for all isolates, MICs for additional antibiotics are included in Supplementary Data 1 and 2.

Decreased phenotypic susceptibility (intermediate and resistant) was defined by MIC breakpoints established by the Clinical and Laboratory Standards Institute (CLSI) (https://clsi.org) for penicillin (≥0.25 μg/mL), cefotaxime (≥2 μg/mL), vancomycin (>1 μg/mL), erythromycin (≥0.5 μg/mL), tetracycline (≥4 μg/mL), and clindamycin (≥0.5 μg /mL). In the absence of CLSI guidelines, the European Committee on Antimicrobial Susceptibility Testing (EUCAST) breakpoints (https://www.eucast.org/) for amoxicillin (>2 μg/mL) and gentamicin were used. Gentamicin MIC of ≤128 μg/mL is associated with low-level intrinsic resistance, and isolates with gentamicin MIC of >128 μg/mL are associated with high-level resistance.

Bacterial DNA extraction, library preparation, and WGS

Genomic DNA for WGS was extracted from S. mitis colonies using the Qiagen DNeasy Blood & Tissue Kit (Qiagen, Germany) according to the manufacturer’s instructions. A pre-lysis step for Gram-positive bacteria using a solution consisting of 30 mg/ml of lysozyme (Merck, Germany) and 50 U/ml mutanolysin (Merck, Germany) dissolved in 1× TE buffer (Promega, UK) was included. The quality of extracted DNA was assessed by agarose gel electrophoresis (0.7%), and by measuring 260/280 and 260/230 ratios on a Nanodrop machine (Thermo Scientific, UK). DNA samples were stored at −20 °C prior to dispatch for WGS. DNA quantification, genomic library preparation, and WGS was done by University College London Pathogen Genomics Unit (PGU). The NEBNext Ultra II DNA Library Prep Kit for Illumina was used for library preparation (New England Biolabs, Ipswich, MA, USA), and the Illumina NextSeq platform (Illumina, San Diego, CA, USA) was used for WGS, which generated paired-end sequence reads of 150 bp in length and 50–100x coverage.

Post-sequencing quality control, genome assembly and speciation

Illumina sequencing reads of presumed IE S. mitis, from our previous work32, were checked for quality using FastQC (version 0.11.9) (https://github.com/s-andrews/FastQC), and trimmed using Trimmomatic77 (version 0.39) and a phred score of at least 33 per read was used as the minimum quality score threshold. De novo genome assembly was performed using default parameters in SPAdes78 (version 3.12), and genome quality was determined using the quality assessment tool for genome assemblies (QUAST version 5.0.2) with default parameters79 (Supplementary Data 7). Taxonomic classification of the sequenced S. mitis, that formed our curated dataset of 322 S. mitis genomes32, was firstly done using KRAKEN80 (version 1.0) against the MiniKraken DB_8GB database, and KRAKEN80 (version 2.0) against the minikraken2_v2_8GB_201904 database using default parameters (Supplementary Fig. 8). Genomes assigned as S. mitis were further screened by applying the online PathogenWatch Speciator tool (https://pathogen.watch/), where the in-house species identification tool applied MASH81 to search a curated NCBI RefSeq database82. Genomes that were not assigned as S. mitis by both KRAKEN versions and PathogenWatch were excluded. In this current analysis, we reanalysed the curated 322 S. mitis genome dataset using the speciation methods described above, and as an additional screening step, average nucleotide identity (ANI) values were calculated using fastANI35 (version 1.32). All strain pairs were tested against each other and against a list of complete S. mitis genomes using the “many to many” method and by using the “–matrix” option. Previous studies have suggested that ANI values of 94–96 % are generally accepted as a species boundary36,37, however, S. mitis has been shown to have lower ANI values of up to 91% as the group consists of a continuum of lineages38. Therefore, a relaxed approach using a 90% ANI threshold was used. Lastly, an S. mitis phylogeny was generated using the methods described below, and species assignment methods were assessed together to confirm the species.

Global S. mitis genomes were used to contextualise locally obtained isolates in a broader perspective. All publicly available S. mitis genome assemblies used in this project were downloaded from The National Center for Biotechnology Information (NCBI) genome database (https://www.ncbi.nlm.nih.gov/) and were from carriage, invasive disease, and unknown conditions (Supplementary Data 3).

Pairwise-SNP distance, phylogeny, and population structure analysis

Snippy (version 4.6.0) (https://github.com/tseemann/snippy) was used to map confirmed UK IE S. mitis sequence reads to the S. mitis B6 reference genome (GenBank Accession: GCA_000027165.1) to obtain SNPs, determine genetic diversity, and the alignment was used to construct maximum-likelihood phylogenies using fasttree83 (version 2.1.10). We used the generalised time-reversible model of nucleotide evolution to generate the phylogenies, which were visualised and annotated using the online Interactive Tree of Life (iToL) software84 (version 3.0) and microreact85 (version 240). Isolates were clustered into GSC using PopPUNK43 (version 2.4.0), and STs were defined using a novel multi-locus sequence typing (MLST) scheme (https://pubmlst.org/organisms/streptococcus-mitis)32.

To obtain a core-genome alignment using global S. mitis, genome assemblies were first annotated using Prokka86 (version 1.13.4), and a core-genome analysis was conducted using Panaroo87 (version 1.2 .9) to obtain a core-genome alignment. An alignment of SNPs was generated from the core-genome alignment using Snp-Sites88 (version 2.5.1), and phylogenies were constructed as described above. Acquired AMR and virulence genes were identified among the streptococci using Abricate (version 0.9.8) (https://github.com/tseemann/abricate). The ResFinder89 and virulence finder90 databases were used as references for AMR genes and virulence genes, respectively. Since very few S. mitis genomes have been sequenced and studied, genotypic resistance was used to determine concordance with phenotypic data.

Bacterial GWAS

We undertook a pilot bacterial genome-wide association study to identify specific genetic changes overrepresented in IE-associated S. mitis isolates when compared to those collected from nasopharyngeal carriage. Due to the high genetic diversity, and the modest dataset size, we only investigated the relative abundance of genes or gene clusters identified from the pan-genome analysis using Panaroo. Because of the extremely high within-species genetic diversity of S. mitis and the challenges of collecting matched isolates from invasive disease and carriage from the same setting and time frame, we performed a two-stage bacterial GWAS analysis. First, we generated a maximum-likelihood phylogenetic tree of recently sequenced and publicly available confirmed S. mitis whole-genome sequences. We annotated the phylogenetic tree with the disease status of the isolates based on the body isolation site, i.e., blood as an ‘IE-associate BSI’ and oropharynx or nasopharynx as ‘asymptomatic carriage’. Second, we selected pairs of genetically closest carriage and invasive disease isolates that shared the most recent ancestors regardless of their genetic divergence. We then pruned the initial phylogenetic tree of all the isolates to remain with a subtree with an equal number of invasive diseases and carriage S. mitis isolates, where each pair of carriage and invasive disease isolates formed monophyletic clades. This approach provided an approximate matching of the isolates, albeit with higher divergence than seen with similar analyses in other bacterial species, such as Staphylococcus aureus56, Staphylococcus epidermidis70, and Mycobacterium tuberculosis91, to allow for a robust assessment of the genes potentially enriched in the carriage and invasive disease isolates. Due to the phylogenetic matching or pairing of the S. mitis isolates, we used the exact McNemar’s test to identify genes or gene clusters overrepresented in IE or carriage-associated isolates. We used the function “mcnemar.test” in the stats (version 4.0.3) R package to perform the exact McNemar’s test. Genes or gene clusters with P value < 0.05 were considered to be statistically significant. Overrepresented genes identified as hypothetical genes were further checked using the online NCBI BLAST tool, its databases, and default parameters to determine any known gene functions. NCBI BLAST matches with the highest total score, sequence coverage, and sequence identity were used to assign potential gene function.

Invasive Streptococcus pneumoniae genomes used for screening of overrepresented S. mitis genes were obtained from the Global Pneumococcal Sequencing Project (Supplementary Data 6). The GPS was screened to identify S. pneumoniae genomes obtained via blood from patients with bacteremia74. The largest collection of S. pneumoniae genomes collected through bacteraemia surveillance was therefore used and is part of the CDC active bacterial core surveillance. Abricate (version 1.0.1) was used with default settings and the overrepresented S. mitis genes as the database to screen invasive pneumococcal genomes for the presence of these genes.

Statistical analysis

Statistical tests and associated diagrams were generated in R (version 2.11.1) (R Core Team 2014; https://www.R-project.org/), GraphPad Prism (version 8.0) (GraphPad Software, San Diego, California, USA), and edited in Inkscape version 1.0.0. Parametric data collected included the age group of the IE cases and were presented as frequencies. Non-parametric data, which included antibiotic minimum inhibitory concentrations (MICs) and pairwise-SNP distances, are presented as individual data points and median values. The Kruskal–Wallis test was used to compare median MICs among isolates grouped by year, and the test was also used to compare population-level genetic diversity by pairwise SNPs and ANI values across the 16-year surveillance period. Statistical significance was defined as p < 0.05.

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Supplementary information

Supplementary Information

Peer Review File

Description of Additional Supplementary Files

Supplementary Data 1

Supplementary Data 2

Supplementary Data 3

Supplementary Data 4

Supplementary Data 5

Supplementary Data 6

Supplementary Data 7

Reporting Summary

Source data

Source Data

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-52120-z.

Acknowledgements

This research was funded by the NIHR project reference 16/136/46 (R.S.H.) and NIHR200652 grant using UK aid from the UK Government to support global health research. R.S.H. is a NIHR Senior Investigator. We thank the British Society for Antimicrobial Chemotherapy (BSAC) for allowing the use of isolates from the BSAC Resistance Surveillance Programme. We are grateful to the Microbiology laboratory staff for the referral of isolates, the AMR and Healthcare-Associated Infections (AMRHAI) Reference Unit for antimicrobial susceptibility testing, and members of the Respiratory and vaccine-preventable bacteria reference unit (RVPBRU) for assistance in isolate identification and UKHSA for funding. We thank Andrea Gori for his expert input into the analysis. The funders had no role in study design, data collection and analysis, decision to publish, and preparation of the manuscript. Views expressed in this publication are those of the author(s) and not necessarily those of the NHS, the NIHR, the UK Department of Health and Social Care, UKHSA, or BSAC.

Author contributions

A.K., C.C., and R.S.H. conceived the study. J.C., K.B., B.P., K.L.H., and V.C. provided the IE isolates through the UKHSA. K.L.H. conducted phenotypic antibiotic susceptibility testing. S.B. contributed additional genomes for analysis. A.K. performed culture and DNA extractions for WGS. D.R. accessed and verified the raw data post-sequencing, while S.D.B. facilitated access to computing clusters and genomic pipelines at the Wellcome Sanger Institute. A.K. and C.C. led and conducted the analyses, with B.K., S.D.B., and R.S.H. offering guidance and suggestions on alternative approaches that were adopted. R.S.H. and C.C. supervised the study. A.K. wrote the first draft of the manuscript, and all authors reviewed and edited it. R.S.H. secured the funding. All authors approved the final manuscript and had ultimate responsibility for the decision to submit it for publication.

Peer review

Peer review information

Nature Communications thanks François Vandenesch, who co-reviewed with Coralie Bouchiat and the other anonymous reviewer(s), for their contribution to the peer review of this work. A peer review file is available.

Data availability

Genomes sequenced in this study have been deposited in the US National Center for Biotechnology Information (NCBI) database under BioProject accession code PRJEB55310. Publicly available genomes used in this project are under BioProjects PRJNA480039, PRJEB42564, PRJEB42963, and PRJEB53188. All genomes used in this study are also shared under genome assembly accessions listed in Supplementary Data 3, 4, and 6. Source data are provided with this paper Source data are provided with this paper.

Code availability

We have described all the tools and methods used for the analysis in the Material and Methods sections.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Holland TL Infective endocarditis Nat. Rev. Dis. Primer 2016 2 1 22 10.1038/nrdp.2016.59
Holland, T. L. et al. Infective endocarditis. Nat. Rev. Dis. Primer 2, 1–22 (2016).10.1038/nrdp.2016.59
2. Prendergast BD The changing face of infective endocarditis Heart 2006 92 879 885 10.1136/hrt.2005.067256 16216860
Prendergast, B. D. The changing face of infective endocarditis. Heart 92, 879–885 (2006).16216860 10.1136/hrt.2005.067256
3. Murdoch DR Clinical presentation, etiology and outcome of infective endocarditis in the 21st century: the international collaboration on endocarditis-prospective cohort study Arch. Intern. Med. 2009 169 463 473 10.1001/archinternmed.2008.603 19273776
Murdoch, D. R. et al. Clinical presentation, etiology and outcome of infective endocarditis in the 21st century: the international collaboration on endocarditis-prospective cohort study. Arch. Intern. Med. 169, 463–473 (2009).19273776 10.1001/archinternmed.2008.603
4. Bin Abdulhak AA Global and regional burden of infective endocarditis, 1990–2010: a systematic review of the literature Glob. Heart 2014 9 131 143 10.1016/j.gheart.2014.01.002 25432123
Bin Abdulhak, A. A. et al. Global and regional burden of infective endocarditis, 1990–2010: a systematic review of the literature. Glob. Heart 9, 131–143 (2014).25432123 10.1016/j.gheart.2014.01.002
5. Pant S Trends in infective endocarditis incidence, microbiology, and valve replacement in the United States from 2000 to 2011 J. Am. Coll. Cardiol. 2015 65 2070 2076 10.1016/j.jacc.2015.03.518 25975469
Pant, S. et al. Trends in infective endocarditis incidence, microbiology, and valve replacement in the United States from 2000 to 2011. J. Am. Coll. Cardiol. 65, 2070–2076 (2015).25975469 10.1016/j.jacc.2015.03.518
6. Talha KM Escalating incidence of infective endocarditis in Europe in the 21st century Open Heart 2021 8 e001846 10.1136/openhrt-2021-001846 34670832
Talha, K. M. et al. Escalating incidence of infective endocarditis in Europe in the 21st century. Open Heart 8, e001846 (2021).34670832 10.1136/openhrt-2021-001846
7. Thornhill MH An alarming rise in incidence of infective endocarditis in England since 2009: why? Lancet 2020 395 1325 1327 10.1016/S0140-6736(20)30530-4 32334690
Thornhill, M. H. et al. An alarming rise in incidence of infective endocarditis in England since 2009: why? Lancet 395, 1325–1327 (2020).32334690 10.1016/S0140-6736(20)30530-4
8. Mitchell J Streptococcus mitis: Walking the line between commensalism and pathogenesis Mol. Oral Microbiol. 2011 26 89 98 10.1111/j.2041-1014.2010.00601.x 21375700
Mitchell, J. Streptococcus mitis: Walking the line between commensalism and pathogenesis. Mol. Oral Microbiol. 26, 89–98 (2011).21375700 10.1111/j.2041-1014.2010.00601.x
9. Smith DJ Anderson JM King WF van Houte J Taubman MA Oral streptococcal colonization of infants Oral Microbiol. Immunol. 1993 8 1 4 10.1111/j.1399-302X.1993.tb00535.x 8510978
Smith, D. J., Anderson, J. M., King, W. F., van Houte, J. & Taubman, M. A. Oral streptococcal colonization of infants. Oral Microbiol. Immunol. 8, 1–4 (1993).8510978 10.1111/j.1399-302X.1993.tb00535.x
10. Osler W Culstonian lectures on malignant endocarditis Lancet 1885 125 505 508 10.1016/S0140-6736(02)17782-8
Osler, W. Culstonian lectures on malignant endocarditis. Lancet 125, 505–508 (1885).10.1016/S0140-6736(02)17782-8
11. Taylor SN Sanders CV Unusual manifestations of invasive pneumococcal infection Am. J. Med. 1999 107 12 27 10.1016/S0002-9343(99)00103-5
Taylor, S. N. & Sanders, C. V. Unusual manifestations of invasive pneumococcal infection. Am. J. Med. 107, 12–27 (1999).10.1016/S0002-9343(99)00103-5
12. Thayer WS Bacterial or infective endocarditis. The Gibson lectures for 1930: Lecture I Edinb. Med. J. 1931 38 237 29639893
Thayer, W. S. Bacterial or infective endocarditis. The Gibson lectures for 1930: Lecture I. Edinb. Med. J. 38, 237 (1931).29639893
13. Austrian R The syndrome of pneumococcal endocarditis, meningitis and rupture of the aortic valve Trans. Am. Clin. Climatol. Assoc. 1957 68 40 50
Austrian, R. The syndrome of pneumococcal endocarditis, meningitis and rupture of the aortic valve. Trans. Am. Clin. Climatol. Assoc. 68, 40–50 (1957).
14. Ruegsegger JM Pneumococcal endocarditis Am. Heart J. 1958 56 867 877 10.1016/0002-8703(58)90195-9 13594831
Ruegsegger, J. M. Pneumococcal endocarditis. Am. Heart J. 56, 867–877 (1958).13594831 10.1016/0002-8703(58)90195-9
15. Chamat-Hedemand Sandra Prevalence of infective endocarditis in Streptococcal bloodstream infections is dependent on streptococcal species Circulation 2020 142 720 730 10.1161/CIRCULATIONAHA.120.046723 32580572
Chamat-Hedemand Sandra et al. Prevalence of infective endocarditis in Streptococcal bloodstream infections is dependent on streptococcal species. Circulation 142, 720–730 (2020).32580572 10.1161/CIRCULATIONAHA.120.046723
16. Kilian M Evolution of Streptococcus pneumoniae and its close commensal relatives PloS One 2008 3 e2683 10.1371/journal.pone.0002683 18628950
Kilian, M. et al. Evolution of Streptococcus pneumoniae and its close commensal relatives. PloS One 3, e2683 (2008).18628950 10.1371/journal.pone.0002683
17. Rai P He F Kwang J Engelward BP Chow VTK Pneumococcal pneumolysin induces DNA damage and cell cycle arrest Sci. Rep. 2016 6 22972 10.1038/srep22972 27026501
Rai, P., He, F., Kwang, J., Engelward, B. P. & Chow, V. T. K. Pneumococcal pneumolysin induces DNA damage and cell cycle arrest. Sci. Rep. 6, 22972 (2016).27026501 10.1038/srep22972
18. Rasmussen LH In silico assessment of virulence factors in strains of Streptococcus oralis and Streptococcus mitis isolated from patients with Infective Endocarditis J. Med. Microbiol. 2017 66 1316 1323 10.1099/jmm.0.000573 28874232
Rasmussen, L. H. et al. In silico assessment of virulence factors in strains of Streptococcus oralis and Streptococcus mitis isolated from patients with Infective Endocarditis. J. Med. Microbiol. 66, 1316–1323 (2017).28874232 10.1099/jmm.0.000573
19. Rukke HV Protective role of the capsule and impact of serotype 4 switching on Streptococcus mitis Infect. Immun. 2014 82 3790 3801 10.1128/IAI.01840-14 24958712
Rukke, H. V. et al. Protective role of the capsule and impact of serotype 4 switching on Streptococcus mitis. Infect. Immun. 82, 3790–3801 (2014).24958712 10.1128/IAI.01840-14
20. Mitchell AM Mitchell TJ Streptococcus pneumoniae: virulence factors and variation Clin. Microbiol. Infect. 2010 16 411 418 10.1111/j.1469-0691.2010.03183.x 20132250
Mitchell, A. M. & Mitchell, T. J. Streptococcus pneumoniae: virulence factors and variation. Clin. Microbiol. Infect. 16, 411–418 (2010).20132250 10.1111/j.1469-0691.2010.03183.x
21. Pimenta F Streptococcus infantis, Streptococcus mitis, and Streptococcus oralis strains with highly similar cps5 loci and antigenic relatedness to serotype 5 pneumococci Front. Microbiol. 2018 9 3199 10.3389/fmicb.2018.03199 30671034
Pimenta, F. et al. Streptococcus infantis, Streptococcus mitis, and Streptococcus oralis strains with highly similar cps5 loci and antigenic relatedness to serotype 5 pneumococci. Front. Microbiol. 9, 3199 (2018).30671034 10.3389/fmicb.2018.03199
22. Lessa FC Streptococcus mitis expressing pneumococcal serotype 1 capsule Sci. Rep. 2018 8 17959 10.1038/s41598-018-35921-3 30568178
Lessa, F. C. et al. Streptococcus mitis expressing pneumococcal serotype 1 capsule. Sci. Rep. 8, 17959 (2018).30568178 10.1038/s41598-018-35921-3
23. Whatmore AM Genetic relationships between clinical isolates of Streptococcus pneumoniae, Streptococcus oralis, and Streptococcus mitis: Characterization of “Atypical” Pneumococci and Organisms Allied to S. mitis Harboring S. pneumoniae Virulence Factor-Encoding Genes Infect. Immun. 2000 68 1374 1382 10.1128/IAI.68.3.1374-1382.2000 10678950
Whatmore, A. M. et al. Genetic relationships between clinical isolates of Streptococcus pneumoniae, Streptococcus oralis, and Streptococcus mitis: Characterization of “Atypical” Pneumococci and Organisms Allied to S. mitis Harboring S. pneumoniae Virulence Factor-Encoding Genes. Infect. Immun. 68, 1374–1382 (2000).10678950 10.1128/IAI.68.3.1374-1382.2000
24. Johnston C Detection of large numbers of Pneumococcal virulence genes in Streptococci of the Mitis group J. Clin. Microbiol. 2010 48 2762 2769 10.1128/JCM.01746-09 20519466
Johnston, C. et al. Detection of large numbers of Pneumococcal virulence genes in Streptococci of the Mitis group. J. Clin. Microbiol. 48, 2762–2769 (2010).20519466 10.1128/JCM.01746-09
25. Kadioglu A Weiser JN Paton JC Andrew PW The role of Streptococcus pneumoniae virulence factors in host respiratory colonization and disease Nat. Rev. Microbiol. 2008 6 288 301 10.1038/nrmicro1871 18340341
Kadioglu, A., Weiser, J. N., Paton, J. C. & Andrew, P. W. The role of Streptococcus pneumoniae virulence factors in host respiratory colonization and disease. Nat. Rev. Microbiol. 6, 288–301 (2008).18340341 10.1038/nrmicro1871
26. Kilian M Riley DR Jensen A Brüggemann H Tettelin H Parallel evolution of Streptococcus pneumoniae and Streptococcus mitis to pathogenic and mutualistic lifestyles mBio 2014 5 e01490-01414 10.1128/mBio.01490-14 25053789
Kilian, M., Riley, D. R., Jensen, A., Brüggemann, H. & Tettelin, H. Parallel evolution of Streptococcus pneumoniae and Streptococcus mitis to pathogenic and mutualistic lifestyles. mBio 5, e01490-01414 (2014).25053789 10.1128/mBio.01490-14
27. Kilian M Tettelin H Identification of virulence-associated properties by comparative genome analysis of Streptococcus pneumoniae, S. pseudopneumoniae, S. mitis, Three S. oralis Subspecies, and S. infantis mBio 2019 10 e01985-19 10.1128/mBio.01985-19 31481387
Kilian, M. & Tettelin, H. Identification of virulence-associated properties by comparative genome analysis of Streptococcus pneumoniae, S. pseudopneumoniae, S. mitis, Three S. oralis Subspecies, and S. infantis. mBio 10, e01985-19 (2019).31481387 10.1128/mBio.01985-19
28. Kalizang’oma A Streptococcus pneumoniae serotypes that frequently colonise the human nasopharynx are common recipients of penicillin-binding protein gene fragments from Streptococcus mitis Microb. Genomics 2021 7 000622
Kalizang’oma, A. et al. Streptococcus pneumoniae serotypes that frequently colonise the human nasopharynx are common recipients of penicillin-binding protein gene fragments from Streptococcus mitis. Microb. Genomics 7, 000622 (2021).
29. Doern CD Burnham C-AD It’s not easy being green: the viridans group Streptococci, with a focus on pediatric clinical manifestations J. Clin. Microbiol. 2010 48 3829 3835 10.1128/JCM.01563-10 20810781
Doern, C. D. & Burnham, C.-A. D. It’s not easy being green: the viridans group Streptococci, with a focus on pediatric clinical manifestations. J. Clin. Microbiol. 48, 3829–3835 (2010).20810781 10.1128/JCM.01563-10
30. Pearce C Identification of pioneer viridans streptococci in the oral cavity of human neonates J. Med. Microbiol. 1995 42 67 72 10.1099/00222615-42-1-67 7739028
Pearce, C. et al. Identification of pioneer viridans streptococci in the oral cavity of human neonates. J. Med. Microbiol. 42, 67–72 (1995).7739028 10.1099/00222615-42-1-67
31. Weiser JN Ferreira DM Paton JC Streptococcus pneumoniae: transmission, colonization and invasion Nat. Rev. Microbiol. 2018 16 355 367 10.1038/s41579-018-0001-8 29599457
Weiser, J. N., Ferreira, D. M. & Paton, J. C. Streptococcus pneumoniae: transmission, colonization and invasion. Nat. Rev. Microbiol. 16, 355–367 (2018).29599457 10.1038/s41579-018-0001-8
32. Kalizang’oma A Novel multilocus sequence typing and global sequence clustering schemes for characterizing the population diversity of Streptococcus mitis J. Clin. Microbiol. 2022 0 e00802-22
Kalizang’oma, A. et al. Novel multilocus sequence typing and global sequence clustering schemes for characterizing the population diversity of Streptococcus mitis. J. Clin. Microbiol. 0, e00802-22 (2022).
33. Delgado V 2023 ESC Guidelines for the management of endocarditis: Developed by the task force on the management of endocarditis of the European Society of Cardiology (ESC) Endorsed by the European Association for Cardio-Thoracic Surgery (EACTS) and the European Association of Nuclear Medicine (EANM) Eur. Heart J. 2023 44 3948 4042 10.1093/eurheartj/ehad193 37622656
Delgado, V. et al. 2023 ESC Guidelines for the management of endocarditis: Developed by the task force on the management of endocarditis of the European Society of Cardiology (ESC) Endorsed by the European Association for Cardio-Thoracic Surgery (EACTS) and the European Association of Nuclear Medicine (EANM). Eur. Heart J. 44, 3948–4042 (2023).37622656 10.1093/eurheartj/ehad193
34. Sadowy E Hryniewicz W Identification of Streptococcus pneumoniae and other Mitis streptococci: importance of molecular methods Eur. J. Clin. Microbiol. Infect. Dis. 2020 39 2247 2256 10.1007/s10096-020-03991-9 32710352
Sadowy, E. & Hryniewicz, W. Identification of Streptococcus pneumoniae and other Mitis streptococci: importance of molecular methods. Eur. J. Clin. Microbiol. Infect. Dis. 39, 2247–2256 (2020).32710352 10.1007/s10096-020-03991-9
35. Jain C Rodriguez-R LM Phillippy AM Konstantinidis KT Aluru S High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries Nat. Commun. 2018 9 5114 10.1038/s41467-018-07641-9 30504855
Jain, C., Rodriguez-R, L. M., Phillippy, A. M., Konstantinidis, K. T. & Aluru, S. High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat. Commun. 9, 5114 (2018).30504855 10.1038/s41467-018-07641-9
36. Konstantinidis KT Tiedje JM Genomic insights that advance the species definition for prokaryotes Proc. Natl. Acad. Sci. 2005 102 2567 2572 10.1073/pnas.0409727102 15701695
Konstantinidis, K. T. & Tiedje, J. M. Genomic insights that advance the species definition for prokaryotes. Proc. Natl. Acad. Sci. 102, 2567–2572 (2005).15701695 10.1073/pnas.0409727102
37. Richter M Rosselló-Móra R Shifting the genomic gold standard for the prokaryotic species definition Proc. Natl. Acad. Sci. 2009 106 19126 19131 10.1073/pnas.0906412106 19855009
Richter, M. & Rosselló-Móra, R. Shifting the genomic gold standard for the prokaryotic species definition. Proc. Natl. Acad. Sci. 106, 19126–19131 (2009).19855009 10.1073/pnas.0906412106
38. Jensen A Scholz CFP Kilian M Re-evaluation of the taxonomy of the Mitis group of the genus Streptococcus based on whole genome phylogenetic analyses, and proposed reclassification of Streptococcus dentisani as Streptococcus oralis subsp. dentisani comb. nov., Streptococcus tigurinus as Streptococcus oralis subsp. tigurinus comb. nov., and Streptococcus oligofermentans as a later synonym of Streptococcus cristatus Int. J. Syst. Evol. Microbiol. 2016 66 4803 4820 10.1099/ijsem.0.001433 27534397
Jensen, A., Scholz, C. F. P. & Kilian, M. Re-evaluation of the taxonomy of the Mitis group of the genus Streptococcus based on whole genome phylogenetic analyses, and proposed reclassification of Streptococcus dentisani as Streptococcus oralis subsp. dentisani comb. nov., Streptococcus tigurinus as Streptococcus oralis subsp. tigurinus comb. nov., and Streptococcus oligofermentans as a later synonym of Streptococcus cristatus. Int. J. Syst. Evol. Microbiol. 66, 4803–4820 (2016).27534397 10.1099/ijsem.0.001433
39. Bek-Thomsen M Tettelin H Hance I Nelson KE Kilian M Population diversity and dynamics of Streptococcus mitis, Streptococcus oralis, and Streptococcus infantis in the upper respiratory tracts of adults, determined by a nonculture strategy Infect. Immun. 2008 76 1889 1896 10.1128/IAI.01511-07 18316382
Bek-Thomsen, M., Tettelin, H., Hance, I., Nelson, K. E. & Kilian, M. Population diversity and dynamics of Streptococcus mitis, Streptococcus oralis, and Streptococcus infantis in the upper respiratory tracts of adults, determined by a nonculture strategy. Infect. Immun. 76, 1889–1896 (2008).18316382 10.1128/IAI.01511-07
40. Kirchherr JL Clonal diversity and turnover of Streptococcus mitis bv. 1 on shedding and nonshedding oral surfaces of human infants during the first year of life Clin. Diagn. Lab. Immunol. 2005 12 1184 1190 16210481
Kirchherr, J. L. et al. Clonal diversity and turnover of Streptococcus mitis bv. 1 on shedding and nonshedding oral surfaces of human infants during the first year of life. Clin. Diagn. Lab. Immunol. 12, 1184–1190 (2005).16210481
41. Hohwy J Reinholdt J Kilian M Population dynamics of streptococcus mitis in its natural habitat Infect. Immun. 2001 69 6055 6063 10.1128/IAI.69.10.6055-6063.2001 11553543
Hohwy, J., Reinholdt, J. & Kilian, M. Population dynamics of streptococcus mitis in its natural habitat. Infect. Immun. 69, 6055–6063 (2001).11553543 10.1128/IAI.69.10.6055-6063.2001
42. Fitzsimmons S Clonal diversity of Streptococcus mitis biovar 1 isolates from the oral cavity of human neonates Clin. Diagn. Lab. Immunol. 1996 3 517 522 10.1128/cdli.3.5.517-522.1996 8877128
Fitzsimmons, S. et al. Clonal diversity of Streptococcus mitis biovar 1 isolates from the oral cavity of human neonates. Clin. Diagn. Lab. Immunol. 3, 517–522 (1996).8877128 10.1128/cdli.3.5.517-522.1996
43. Lees JA Fast and flexible bacterial genomic epidemiology with PopPUNK Genome Res. 2019 29 304 316 10.1101/gr.241455.118 30679308
Lees, J. A. et al. Fast and flexible bacterial genomic epidemiology with PopPUNK. Genome Res. 29, 304–316 (2019).30679308 10.1101/gr.241455.118
44. Habib G 2015 ESC Guidelines for the management of infective endocarditis: the Task Force for the Management of Infective Endocarditis of the European Society of Cardiology (ESC) Endorsed by: European Association for Cardio-Thoracic Surgery (EACTS), the European Association of Nuclear Medicine (EANM) Eur. Heart J. 2015 36 3075 3128 10.1093/eurheartj/ehv319 26320109
Habib, G. et al. 2015 ESC Guidelines for the management of infective endocarditis: the Task Force for the Management of Infective Endocarditis of the European Society of Cardiology (ESC) Endorsed by: European Association for Cardio-Thoracic Surgery (EACTS), the European Association of Nuclear Medicine (EANM). Eur. Heart J. 36, 3075–3128 (2015).26320109 10.1093/eurheartj/ehv319
45. Malhotra-Kumar S Oropharyngeal carriage of macrolide-resistant viridans group streptococci: a prevalence study among healthy adults in Belgium J. Antimicrob. Chemother. 2004 53 271 276 10.1093/jac/dkh026 14729750
Malhotra-Kumar, S. et al. Oropharyngeal carriage of macrolide-resistant viridans group streptococci: a prevalence study among healthy adults in Belgium. J. Antimicrob. Chemother. 53, 271–276 (2004).14729750 10.1093/jac/dkh026
46. Santoro F Vianna ME Roberts AP Variation on a theme; an overview of the Tn916/Tn1545 family of mobile genetic elements in the oral and nasopharyngeal streptococci Front. Microbiol. 2014 5 535 10.3389/fmicb.2014.00535 25368607
Santoro, F., Vianna, M. E. & Roberts, A. P. Variation on a theme; an overview of the Tn916/Tn1545 family of mobile genetic elements in the oral and nasopharyngeal streptococci. Front. Microbiol. 5, 535 (2014).25368607 10.3389/fmicb.2014.00535
47. Chancey ST Composite mobile genetic elements disseminating macrolide resistance in Streptococcus pneumoniae Front. Microbiol. 2015 6 26 10.3389/fmicb.2015.00026 25709602
Chancey, S. T. et al. Composite mobile genetic elements disseminating macrolide resistance in Streptococcus pneumoniae. Front. Microbiol. 6, 26 (2015).25709602 10.3389/fmicb.2015.00026
48. Turner CE The emergence of successful streptococcus pyogenes lineages through convergent pathways of capsule loss and recombination directing high toxin expression mBio 2019 10 e02521-19 10.1128/mBio.02521-19 31822586
Turner, C. E. et al. The emergence of successful streptococcus pyogenes lineages through convergent pathways of capsule loss and recombination directing high toxin expression. mBio 10, e02521-19 (2019).31822586 10.1128/mBio.02521-19
49. Zhi X Emerging invasive group A streptococcus M1UK lineage detected by allele-specific PCR, England, 20201 Emerg. Infect. Dis. 2023 29 1007 1010 10.3201/eid2905.221887 37019153
Zhi, X. et al. Emerging invasive group A streptococcus M1UK lineage detected by allele-specific PCR, England, 20201. Emerg. Infect. Dis. 29, 1007–1010 (2023).37019153 10.3201/eid2905.221887
50. Lo SW Emergence of a multidrug-resistant and virulent Streptococcus pneumoniae lineage mediates serotype replacement after PCV13: an international whole-genome sequencing study Lancet Microbe 2022 3 e735 e743 10.1016/S2666-5247(22)00158-6 35985351
Lo, S. W. et al. Emergence of a multidrug-resistant and virulent Streptococcus pneumoniae lineage mediates serotype replacement after PCV13: an international whole-genome sequencing study. Lancet Microbe 3, e735–e743 (2022).35985351 10.1016/S2666-5247(22)00158-6
51. Jamrozy D Increasing incidence of group B streptococcus neonatal infections in the Netherlands is associated with clonal expansion of CC17 and CC23 Sci. Rep. 2020 10 9539 10.1038/s41598-020-66214-3 32533007
Jamrozy, D. et al. Increasing incidence of group B streptococcus neonatal infections in the Netherlands is associated with clonal expansion of CC17 and CC23. Sci. Rep. 10, 9539 (2020).32533007 10.1038/s41598-020-66214-3
52. Lo SW Pneumococcal lineages associated with serotype replacement and antibiotic resistance in childhood invasive pneumococcal disease in the post-PCV13 era: an international whole-genome sequencing study Lancet Infect. Dis. 2019 19 759 769 10.1016/S1473-3099(19)30297-X 31196809
Lo, S. W. et al. Pneumococcal lineages associated with serotype replacement and antibiotic resistance in childhood invasive pneumococcal disease in the post-PCV13 era: an international whole-genome sequencing study. Lancet Infect. Dis. 19, 759–769 (2019).31196809 10.1016/S1473-3099(19)30297-X
53. Coffey TJ Dowson CG Daniels M Spratt BG Horizontal spread of an altered penicillin-binding protein 2B gene between Streptococcus pneumoniae and Streptococcus oralis FEMS Microbiol. Lett. 1993 110 335 339 10.1111/j.1574-6968.1993.tb06345.x 8354467
Coffey, T. J., Dowson, C. G., Daniels, M. & Spratt, B. G. Horizontal spread of an altered penicillin-binding protein 2B gene between Streptococcus pneumoniae and Streptococcus oralis. FEMS Microbiol. Lett. 110, 335–339 (1993).8354467 10.1111/j.1574-6968.1993.tb06345.x
54. Bek-Thomsen M Poulsen K Kilian M Occurrence and evolution of the paralogous zinc metalloproteases IgA1 protease, ZmpB, ZmpC, and ZmpD in streptococcus pneumoniae and related commensal species mBio 2012 3 e00303-12 10.1128/mBio.00303-12 23033471
Bek-Thomsen, M., Poulsen, K. & Kilian, M. Occurrence and evolution of the paralogous zinc metalloproteases IgA1 protease, ZmpB, ZmpC, and ZmpD in streptococcus pneumoniae and related commensal species. mBio 3, e00303-12 (2012).23033471 10.1128/mBio.00303-12
55. Harth-Chu EN PcsB expression diversity influences on streptococcus mitis phenotypes associated with host persistence and virulence Front. Microbiol. 2019 10 2567 10.3389/fmicb.2019.02567 31798545
Harth-Chu, E. N. et al. PcsB expression diversity influences on streptococcus mitis phenotypes associated with host persistence and virulence. Front. Microbiol. 10, 2567 (2019).31798545 10.3389/fmicb.2019.02567
56. Chaguza C Prophage-encoded immune evasion factors are critical for Staphylococcus aureus host infection, switching, and adaptation Cell Genomics 2022 2 100194 10.1016/j.xgen.2022.100194 36465278
Chaguza, C. et al. Prophage-encoded immune evasion factors are critical for Staphylococcus aureus host infection, switching, and adaptation. Cell Genomics 2, 100194 (2022).36465278 10.1016/j.xgen.2022.100194
57. Power RA Parkhill J de Oliveira T Microbial genome-wide association studies: lessons from human GWAS Nat. Rev. Genet. 2017 18 41 50 10.1038/nrg.2016.132 27840430
Power, R. A., Parkhill, J. & de Oliveira, T. Microbial genome-wide association studies: lessons from human GWAS. Nat. Rev. Genet. 18, 41–50 (2017).27840430 10.1038/nrg.2016.132
58. Wee BA Population analysis of Legionella pneumophila reveals a basis for resistance to complement-mediated killing Nat. Commun. 2021 12 7165 10.1038/s41467-021-27478-z 34887398
Wee, B. A. et al. Population analysis of Legionella pneumophila reveals a basis for resistance to complement-mediated killing. Nat. Commun. 12, 7165 (2021).34887398 10.1038/s41467-021-27478-z
59. Chewapreecha C Comprehensive identification of single nucleotide polymorphisms associated with beta-lactam resistance within pneumococcal mosaic genes PLoS Genet. 2014 10 e1004547 10.1371/journal.pgen.1004547 25101644
Chewapreecha, C. et al. Comprehensive identification of single nucleotide polymorphisms associated with beta-lactam resistance within pneumococcal mosaic genes. PLoS Genet. 10, e1004547 (2014).25101644 10.1371/journal.pgen.1004547
60. Farhat MR GWAS for quantitative resistance phenotypes in Mycobacterium tuberculosis reveals resistance genes and regulatory regions Nat. Commun. 2019 10 2128 10.1038/s41467-019-10110-6 31086182
Farhat, M. R. et al. GWAS for quantitative resistance phenotypes in Mycobacterium tuberculosis reveals resistance genes and regulatory regions. Nat. Commun. 10, 2128 (2019).31086182 10.1038/s41467-019-10110-6
61. Berry AM Lock RA Hansman D Paton JC Contribution of autolysin to virulence of Streptococcus pneumoniae Infect. Immun. 1989 57 2324 2330 10.1128/iai.57.8.2324-2330.1989 2568343
Berry, A. M., Lock, R. A., Hansman, D. & Paton, J. C. Contribution of autolysin to virulence of Streptococcus pneumoniae. Infect. Immun. 57, 2324–2330 (1989).2568343 10.1128/iai.57.8.2324-2330.1989
62. Whatmore AM Molecular characterization of equine isolates of streptococcus pneumoniae: natural disruption of genes encoding the virulence factors pneumolysin and autolysin Infect. Immun. 1999 67 2776 2782 10.1128/IAI.67.6.2776-2782.1999 10338480
Whatmore, A. M. et al. Molecular characterization of equine isolates of streptococcus pneumoniae: natural disruption of genes encoding the virulence factors pneumolysin and autolysin. Infect. Immun. 67, 2776–2782 (1999).10338480 10.1128/IAI.67.6.2776-2782.1999
63. Canvin JR The role of pneumolysin and autolysin in the pathology of pneumonia and septicemia in mice infected with a type 2 pneumococcus J. Infect. Dis. 1995 172 119 123 10.1093/infdis/172.1.119 7797901
Canvin, J. R. et al. The role of pneumolysin and autolysin in the pathology of pneumonia and septicemia in mice infected with a type 2 pneumococcus. J. Infect. Dis. 172, 119–123 (1995).7797901 10.1093/infdis/172.1.119
64. Sitkiewicz I How to become a killer, or is it all accidental? Virulence strategies in oral streptococci Mol. Oral Microbiol. 2018 33 1 12 10.1111/omi.12192 28727895
Sitkiewicz, I. How to become a killer, or is it all accidental? Virulence strategies in oral streptococci. Mol. Oral Microbiol. 33, 1–12 (2018).28727895 10.1111/omi.12192
65. Catto BA Jacobs MR Shlaes DM Streptococcus mitis. A cause of serious infection in adults Arch. Intern. Med. 1987 147 885 888 10.1001/archinte.1987.00370050081014 2437875
Catto, B. A., Jacobs, M. R. & Shlaes, D. M. Streptococcus mitis. A cause of serious infection in adults. Arch. Intern. Med. 147, 885–888 (1987).2437875 10.1001/archinte.1987.00370050081014
66. Welker M Van Belkum A Girard V Charrier J-P Pincus D An update on the routine application of MALDI-TOF MS in clinical microbiology Expert Rev. Proteomics 2019 16 695 710 10.1080/14789450.2019.1645603 31315000
Welker, M., Van Belkum, A., Girard, V., Charrier, J.-P. & Pincus, D. An update on the routine application of MALDI-TOF MS in clinical microbiology. Expert Rev. Proteomics 16, 695–710 (2019).31315000 10.1080/14789450.2019.1645603
67. Isaksson J Comparison of species identification of endocarditis associated viridans streptococci using rnpB genotyping and 2 MALDI-TOF systems Diagn. Microbiol. Infect. Dis. 2015 81 240 245 10.1016/j.diagmicrobio.2014.12.007 25616316
Isaksson, J. et al. Comparison of species identification of endocarditis associated viridans streptococci using rnpB genotyping and 2 MALDI-TOF systems. Diagn. Microbiol. Infect. Dis. 81, 240–245 (2015).25616316 10.1016/j.diagmicrobio.2014.12.007
68. Singh N Antibiotic susceptibility patterns of viridans group streptococci isolates in the United States from 2010 to 2020 JAC-Antimicrob. Resist. 2022 4 dlac049 10.1093/jacamr/dlac049 35599725
Singh, N. et al. Antibiotic susceptibility patterns of viridans group streptococci isolates in the United States from 2010 to 2020. JAC-Antimicrob. Resist. 4, dlac049 (2022).35599725 10.1093/jacamr/dlac049
69. Chun S Huh HJ Lee NY Species-specific difference in antimicrobial susceptibility among viridans group streptococci Ann. Lab. Med. 2015 35 205 211 10.3343/alm.2015.35.2.205 25729722
Chun, S., Huh, H. J. & Lee, N. Y. Species-specific difference in antimicrobial susceptibility among viridans group streptococci. Ann. Lab. Med. 35, 205–211 (2015).25729722 10.3343/alm.2015.35.2.205
70. Méric G Disease-associated genotypes of the commensal skin bacterium Staphylococcus epidermidis Nat. Commun. 2018 9 5034 10.1038/s41467-018-07368-7 30487573
Méric, G. et al. Disease-associated genotypes of the commensal skin bacterium Staphylococcus epidermidis. Nat. Commun. 9, 5034 (2018).30487573 10.1038/s41467-018-07368-7
71. Holmes AR The pavA gene of Streptococcus pneumoniae encodes a fibronectin-binding protein that is essential for virulence Mol. Microbiol. 2001 41 1395 1408 10.1046/j.1365-2958.2001.02610.x 11580843
Holmes, A. R. et al. The pavA gene of Streptococcus pneumoniae encodes a fibronectin-binding protein that is essential for virulence. Mol. Microbiol. 41, 1395–1408 (2001).11580843 10.1046/j.1365-2958.2001.02610.x
72. Gori A Pan-GWAS of streptococcus agalactiae highlights lineage-specific genes associated with virulence and niche adaptation mBio 2020 11 e00728-20 10.1128/mBio.00728-20 32518186
Gori, A. et al. Pan-GWAS of streptococcus agalactiae highlights lineage-specific genes associated with virulence and niche adaptation. mBio 11, e00728-20 (2020).32518186 10.1128/mBio.00728-20
73. Salvadori G Junges R Morrison DA Petersen FC Competence in Streptococcus pneumoniae and close commensal relatives: mechanisms and implications Front. Cell. Infect. Microbiol. 2019 9 94 10.3389/fcimb.2019.00094 31001492
Salvadori, G., Junges, R., Morrison, D. A. & Petersen, F. C. Competence in Streptococcus pneumoniae and close commensal relatives: mechanisms and implications. Front. Cell. Infect. Microbiol. 9, 94 (2019).31001492 10.3389/fcimb.2019.00094
74. Gladstone RA International genomic definition of pneumococcal lineages, to contextualise disease, antibiotic resistance and vaccine impact EBioMedicine 2019 43 338 346 10.1016/j.ebiom.2019.04.021 31003929
Gladstone, R. A. et al. International genomic definition of pneumococcal lineages, to contextualise disease, antibiotic resistance and vaccine impact. EBioMedicine 43, 338–346 (2019).31003929 10.1016/j.ebiom.2019.04.021
75. Delgado V 2023 ESC Guidelines for the management of endocarditis Eur. Heart J. 2023 44 3948 4042 10.1093/eurheartj/ehad193 37622656
Delgado, V. et al. 2023 ESC Guidelines for the management of endocarditis. Eur. Heart J. 44, 3948–4042 (2023).37622656 10.1093/eurheartj/ehad193
76. White AR On behalf of the BSAC Working Parties on Resistance Surveillance. The British Society for Antimicrobial Chemotherapy Resistance Surveillance Project: a successful collaborative model J. Antimicrob. Chemother. 2008 62 ii3 ii14 10.1093/jac/dkn348 18819978
White, A. R., On behalf of the BSAC Working Parties on Resistance Surveillance. The British Society for Antimicrobial Chemotherapy Resistance Surveillance Project: a successful collaborative model. J. Antimicrob. Chemother. 62, ii3–ii14 (2008).18819978 10.1093/jac/dkn348
77. Bolger AM Lohse M Usadel B Trimmomatic: a flexible trimmer for Illumina sequence data Bioinformatics 2014 30 2114 2120 10.1093/bioinformatics/btu170 24695404
Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120 (2014).24695404 10.1093/bioinformatics/btu170
78. Bankevich A SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing J. Comput. Biol. 2012 19 455 477 10.1089/cmb.2012.0021 22506599
Bankevich, A. et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 19, 455–477 (2012).22506599 10.1089/cmb.2012.0021
79. Gurevich A Saveliev V Vyahhi N Tesler G QUAST: quality assessment tool for genome assemblies Bioinformatics 2013 29 1072 1075 10.1093/bioinformatics/btt086 23422339
Gurevich, A., Saveliev, V., Vyahhi, N. & Tesler, G. QUAST: quality assessment tool for genome assemblies. Bioinformatics 29, 1072–1075 (2013).23422339 10.1093/bioinformatics/btt086
80. Wood DE Salzberg SL Kraken: ultrafast metagenomic sequence classification using exact alignments Genome Biol. 2014 15 R46 10.1186/gb-2014-15-3-r46 24580807
Wood, D. E. & Salzberg, S. L. Kraken: ultrafast metagenomic sequence classification using exact alignments. Genome Biol. 15, R46 (2014).24580807 10.1186/gb-2014-15-3-r46
81. Ondov BD Mash: fast genome and metagenome distance estimation using MinHash Genome Biol. 2016 17 132 10.1186/s13059-016-0997-x 27323842
Ondov, B. D. et al. Mash: fast genome and metagenome distance estimation using MinHash. Genome Biol. 17, 132 (2016).27323842 10.1186/s13059-016-0997-x
82. Lam MMC A genomic surveillance framework and genotyping tool for Klebsiella pneumoniae and its related species complex Nat. Commun. 2021 12 4188 10.1038/s41467-021-24448-3 34234121
Lam, M. M. C. et al. A genomic surveillance framework and genotyping tool for Klebsiella pneumoniae and its related species complex. Nat. Commun. 12, 4188 (2021).34234121 10.1038/s41467-021-24448-3
83. Price MN Dehal PS Arkin AP FastTree 2 – approximately maximum-likelihood trees for large alignments PLoS One 2010 5 e9490 10.1371/journal.pone.0009490 20224823
Price, M. N., Dehal, P. S. & Arkin, A. P. FastTree 2 – approximately maximum-likelihood trees for large alignments. PLoS One 5, e9490 (2010).20224823 10.1371/journal.pone.0009490
84. Letunic I Bork P Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees Nucleic Acids Res. 2016 44 W242 W245 10.1093/nar/gkw290 27095192
Letunic, I. & Bork, P. Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 44, W242–W245 (2016).27095192 10.1093/nar/gkw290
85. Argimón S Microreact: visualizing and sharing data for genomic epidemiology and phylogeography Microb. Genom. 2016 2 e000093 28348833
Argimón, S. et al. Microreact: visualizing and sharing data for genomic epidemiology and phylogeography. Microb. Genom. 2, e000093 (2016).28348833
86. Seemann T Prokka: rapid prokaryotic genome annotation Bioinformatics 2014 30 2068 2069 10.1093/bioinformatics/btu153 24642063
Seemann, T. Prokka: rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069 (2014).24642063 10.1093/bioinformatics/btu153
87. Tonkin-Hill G Producing polished prokaryotic pangenomes with the Panaroo pipeline Genome Biol. 2020 21 180 10.1186/s13059-020-02090-4 32698896
Tonkin-Hill, G. et al. Producing polished prokaryotic pangenomes with the Panaroo pipeline. Genome Biol. 21, 180 (2020).32698896 10.1186/s13059-020-02090-4
88. Page AJ SNP-sites: rapid efficient extraction of SNPs from multi-FASTA alignments Microb. Genomics 2016 2 e000056 10.1099/mgen.0.000056
Page, A. J. et al. SNP-sites: rapid efficient extraction of SNPs from multi-FASTA alignments. Microb. Genomics 2, e000056 (2016).10.1099/mgen.0.000056
89. Bortolaia, V. et al. ResFinder 4.0 for predictions of phenotypes from genotypes. J. Antimicrob. Chemother. 10.1093/jac/dkaa345 (2020).
90. Chen L Zheng D Liu B Yang J Jin Q VFDB 2016: hierarchical and refined dataset for big data analysis–10 years on Nucleic Acids Res. 2016 44 D694 D697 10.1093/nar/gkv1239 26578559
Chen, L., Zheng, D., Liu, B., Yang, J. & Jin, Q. VFDB 2016: hierarchical and refined dataset for big data analysis–10 years on. Nucleic Acids Res. 44, D694–D697 (2016).26578559 10.1093/nar/gkv1239
91. Farhat MR Shapiro BJ Sheppard SK Colijn C Murray M A phylogeny-based sampling strategy and power calculator informs genome-wide associations study design for microbial pathogens Genome Med. 2014 6 101 10.1186/s13073-014-0101-7 25484920
Farhat, M. R., Shapiro, B. J., Sheppard, S. K., Colijn, C. & Murray, M. A phylogeny-based sampling strategy and power calculator informs genome-wide associations study design for microbial pathogens. Genome Med. 6, 101 (2014).25484920 10.1186/s13073-014-0101-7
92. Li Y Penicillin-binding protein transpeptidase signatures for tracking and predicting β-lactam resistance levels in Streptococcus pneumoniae mBio 2016 7 e00756-16 10.1128/mBio.00756-16 27302760
Li, Y. et al. Penicillin-binding protein transpeptidase signatures for tracking and predicting β-lactam resistance levels in Streptococcus pneumoniae. mBio 7, e00756-16 (2016).27302760 10.1128/mBio.00756-16
