Population structure, epidemiology and antibiotic resistance patterns of Streptococcus pneumoniae serotype 5: prior to PCV-13 vaccine introduction in Eastern Gambia
© Ashu et al. 2016
Received: 11 May 2015
Accepted: 22 January 2016
Published: 28 January 2016
Streptococcus pneumoniae serotype 5 is among the most common serotypes causing invasive pneumococcal disease (IPD) in The Gambia. We anticipate that introduction of the 13-valent pneumococcal conjugate vaccine (PCV-13) into routine vaccination in The Gambia will reduce serotype 5 IPD. However, the emergence of new clones that have altered their genetic repertoire through capsular switching or genetic recombination after vaccination with PCV-13 poses a threat to this public health effort. In order to monitor for potential genetic changes post-PCV-13 vaccination, we established the baseline population structure, epidemiology, and antibiotic resistance patterns of serotype 5 before the introduction of PCV-13.
Fifty-five invasive S. pneumoniae serotype 5 isolates were recovered from January 2009 to August 2011 in a population-based study in the Upper River Region of The Gambia. Serotyping was done by latex agglutination and confirmed by serotype-specific Polymerase Chain Reaction (PCR). Genotyping was undertaken using Multilocus Sequence Typing (MLST). Antimicrobial sensitivity was done using disc diffusion. Contingency table analyses were conducted using Pearson’s Chi2 and Fisher’s exact test. Clustering was performed using Bionumerics version 6.5.
MLST resolved S. pneumoniae serotype 5 isolates into 3 sequence types (ST), namely ST 289(6/55), ST 3339(19/55) and ST 3404(30/55). ST 289 was identified as the major clonal complex. ST 3339, the prevalent genotype in 2009 [84.6 % (11/13)], was replaced by ST 3404 [70.4 % (19/27)] in 2010 as the dominant ST. Interestingly, ST 3404 showed lower resistance to tetracycline and oxacillin (P < 0.001), an empirical surrogate to penicillin in The Gambia.
There has been an emergence of ST 3404 in The Gambia prior to the introduction of PCV-13. Our findings provide important background data for future assessment of the impact of PCV-13 into routine immunization in developing countries, such as The Gambia.
KeywordsS. pneumoniae Serotype 5 ST 3404 Population Structure The Gambia PCV-13
Streptococcus pneumoniae is a Gram positive bacterium that colonizes the upper respiratory tract of humans and is capable of causing life-threatening infections, particularly among immunosuppressed individuals [1, 2]. S. pneumoniae has a circular genome containing between 2 million to 2.1 million base pairs . Two-thirds of its 2236 open reading frames have been assigned roles based on their predicted gene products [2, 3]. Conserved proteins with unknown functions constitute 16 % of S. pneumoniae’s open reading frames, of which, 20 % are unique to S. pneumoniae [2, 3]. S. pneumoniae has 1553 and 154 genes that are important for viability and virulence, respectively. Genetic information between strains can vary by up to 10 % [2, 4].
Based on its capsular structure, S. pneumoniae can be categorized into more than 90 different serotypes which differ in their infection ability and severity [5, 6]. S. pneumoniae causes invasive pneumococcal diseases (IPD) such as septicemia and meningitis, and is the leading cause of bacterial disease in The Gambia . Carriage of S. pneumoniae among Gambian children aged less than 1 year is 97 % and 93 % among babies aged less than 1 month . Serotype 5 is one of the most common serotypes causing IPD in The Gambia [9, 10]. Pneumococcal diseases in infants and young children can be prevented by the pneumococcal conjugate vaccine (PCV). Three PCV vaccines — PCV-7 (7-valent vaccine), PCV-10 (10-valent vaccine), and PCV-13 (13-valent vaccine) — have been used in routine immunization in many countries. Even though IPD in The Gambia is commonly caused by serotype 5 , the 13-valent vaccine containing serotype 5 was only recently introduced in July 2011.
The incidence of IPD due to serotype 5 in Eastern Gambia in 2008 was estimated to be 7 (95 % CI 0, 38) per 100 000 person-years in children aged 2–59 months. Two years later in 2010, the incidence of serotype 5 rose to 51 (95 % CI 28, 86) per 100 000 person-years [Unpublished data, Grant Mackenzie, October 2011]. Understanding whether the introduction of PCV-13 will have an impact on the population structure, epidemiology and antibiotic resistance patterns of serotype 5 is important to scientists, public health practitioners and policy makers. The introduction of PCV-13 into the routine vaccination schedule in July 2011 is expected to lead to the reduction of serotype 5 IPD. However, the possibility of the emergence of new clones through capsular switching or recombination after vaccination is a threat to the public health impact of the vaccination program. Following introduction of PCV-7, a study in Canada showed the emergence of multi-drug resistant (MDR) serotype 19A through recombination . The Canadian study painted a complex genetic picture where vaccine-selection pressure, high-level drug resistance, and S. pneumoniae mutational events created a “perfect storm” for the emergence of the multi-drug-resistant 19A serotype . However, in order to surveil for potential post-vaccine changes, it is essential to establish a reference point from which a comparison could be made. The objective of this study was therefore to determine the population structure, epidemiology, and resistance patterns of S. pneumoniae serotype 5 prior to PCV-13’s introduction in to The Gambia, making available a reference point to evaluate the impact of PCV-13’s introduction into routine immunization.
Study area and ethical statement
The Pneumococcal Surveillance Project (PSP), hosted by the Medical Research Council Unit (MRC) in collaboration with the Gambian government is aimed at evaluating the effectiveness of the introduction of PCV into The Gambia. The PSP study area covers the south bank of the Upper River Region (1111 km2), which has a population of about 160,000 inhabitants . The surveillance area is set in the Basse Health and Demographic Surveillance System, in which the population is enumerated every four months. The study base is located in the town of Basse. The study population comprised all cases of pneumonia, sepsis and meningitis among patients aged 2 weeks and older in the study area. This study was conducted in accordance with the Declaration of Helsinki and was approved by the Gambia Government/MRC Joint Ethics Committee (EC) — Scientific Coordinating Committee (SCC) number 1087. Written informed consent was obtained from all study participants.
Isolation and antimicrobial susceptibility testing
Between January 2009 and August 2011, 55 serotype 5 isolates were recovered from 19 female and 24 male in a population-based study in Eastern Gambia. Isolates were recovered from either blood, cerebrospinal fluid, or lung aspirates. A total of 13 isolates were obtained in 2009, while 27 and 15 were obtained in 2010 and 2011, respectively. More than one isolate was recovered from 11 patients. Ten patients provided 2 isolates each while 1 patient provided 3 isolates. When more than one isolate was obtained from the same patient they were obtained at separate time intervals. S. pneumoniae was identified as previously described . Oxoid antibiotic susceptibility disc (Thermofisher Scientific, UK) used to test for antimicrobial susceptibility included co-trimoxazole (25 μg), chloramphenicol (30 μg), tetracycline (30 μg), erythromycin (15 μg) and oxacillin (1 μg). Colonies of S. pneumonie from an overnight Mueller-Hinton agar plate supplemented with 5 % sheep blood (MHBA) were suspended in 2 ml of Mueller-Hinton broth (MHB). Tubes containing MHB were briefly vortexed to achieve a uniform suspension and turbidity was adjusted to a 0.5 McFarland standard (Biomerieux SA, France) equivalent to 1.5 × 10^8 cfu/ml. Sterile cotton tips were thereafter dipped into the suspensions and streaked to cover the entire surface of MHBA plates ensuring even spread. Antimicrobial discs were carefully placed on plates under antiseptic conditions. Plates were incubated overnight at 37 °C in 5 % CO2. Zones of Inhibition (ZOI) were measured and interpreted as resistant, intermediate or susceptible, using Clinical Laboratory Standard Institute (CLSI) break-point, ATCC 49619 was used as a control .
Oligonucleotide primer sequences used in this study
Primer sequence (5′- 3′)
GenBank accession no.
ATA CCT ACA CAA CTT CTG ATT ATG CCT TTG TG
GCT CGA TAA ACA TAA TCA ATA TTT GAA AAA GTA TG
GCA GTA CAG CAG TTT GTT GGA CTG ACC
GAA TAT TTT CAT TAT CAG TCC CAG TC
Genotyping was done by means of Multilocus Sequencing Typing (MLST). Seven housekeeping genes including aroE (shikimate dehydrogenase), ddl (D-alanine-D-alanine ligase), gdh (glucose-6-phosphate dehydrogenase), gki (glucose kinase), recP (transketolase), spi (signal peptidase I) and xpt (xanthine phosphoribosyltransferase) were targeted for amplification by PCR. The PCRs were performed in 25 μl volumes per reaction as previously described . Thermal cycling was performed in a Palm cycler (Corbett research) under the following conditions: 95 °C for 10 min followed by 30 amplification cycles of 95 °C for 30 s, 55 °C for 30 s, 72 °C for 1 min, and a final extension at 72 °C for 10 min. Sequencing was done by incorporation of chain-terminating dideoxynucleotides . A sequencing reaction was performed for both forward and reverse strands of each targeted gene (Big Dye Terminator Cycle Sequencing kit; Applied Biosystems, UK). PCR products were purified using Qiagen purfication kit (Qiagen Sample & Assay Technologies, UK) and submitted to the in-house core sequencing facility in Fajara (Medical Research Council Unit, The Gambia).
Editing and alignment of sequences for submission to MLST database  for allele and sequence type numbers was performed using Laser Gene DNA star, version 7.1. Data was exported into an Excel sheet (Microsoft 2007) for analyses of temporal trends. STATA version 9 (Stata Corporation, College Station, TX, USA) was used to produce summary data. The Wilcoxon signed ranked test was used to compare differences in the number of serotype 5 cases between male and female. In order to illustrate the differences between various groups with categorical data, contingency table analyses were conducted using Pearson’s Chi2 and Fisher’s exact test (two-sided). Genotype clustering was performed using Bionumerics version 6.5 (Applied Maths, Saint-Martens-Latem, Belgium).
Results and discussion
STs 3339 and 3404 differed from ST289 at two loci, ddl and spi, respectively. ddl codes for an enzyme involved in cell wall biosynthesis while spi has been associated with the processing of signal peptides for classical MHC class I molecules [18, 19]. It would be luring to infer that the changes in ddl and spi giving rise to clones ST3339 and ST3404 are of little significance because they are non-capsular (serotype) changes. For example, ddl gene diversity has been linked to the recombinational replacement of penicillin binding protein 2b (pbp2b), a gene located 783 bp upstream of ddl . This linkage suggests that changes in beta-lactam susceptibility patterns involving pbp2b could be reflected in ddl diversity. Consequently, observed disparities in ddl and even spi incidence could potentially be very important. These findings taken together emphasize the need for whole genome analysis in order to determine the relevance of the emergence of ST 3404 and the prevalence of ST3339 in The Gambia, especially because both clones originate from a recent and local differentiation event.
The median age of patients with invasive disease caused by serotype 5 was 29 months (lower quartile = 0.8 months, upper quartile = 240 months). We report no significant difference between the prevalence of ST3404 or serotype 5 — including all three clones — in male and female patients. These findings suggest that sex is not a risk factor as pertains to the acquisition of ST 3404 and of S. pneumoniae serotype 5 as a whole. However, these findings are not unique to serotype 5; the same has been reported about serotype 1 and most other pneumococcal serotypes that cause IPD in The Gambia .
The most number of serotype 5 isolates collected were recorded between the months of May to June each year. These seasonal peaks could be explained by the fact that in The Gambia this period marks the end of the dry season. During the dry seasons in the West African belt, dry and dusty winds damage the mucosa membrane of the oral cavity . Low absolute humidity and dust inhibit the mucosal immune defenses by damaging the mucosal membrane leading to an increased risk of transmission during the dry season. Transmission capacities of other bacterial pathogens such as Meningococcus have been reported to be less during the rainy seasons due to higher humidity [22, 23]. Serotype 5, like most other serotypes in The Gambia, peaks after the dry and hot season . Given that Meningococcus and Streptococcus are both transmitted by close contact with respiratory fluids, assessment of the contribution of low absolute humidity and dust as risk factors for outbreaks of S. pneumoniae is of interest in The Gambia.
All 55 S. pneumoniae serotype 5 isolates were resistant to co-trimazole. Among these isolates, 44 % (24/55) and 73 % (40/55) were susceptible to tetracycline and oxacillin, respectively. Furthermore, 98 % (54/55) of all isolates were susceptible to chloramphenicol. Observed chloramphenicol susceptibility is consistent with a previous study in The Gambia which showed that all invasive S. Pneumoniae serotype 1 isolates were susceptible to chloramphenicol. We however note that susceptibility patterns vary slightly by serotype with respect to tetracycline, co-trimazole, and oxacillin . The proportion of serotype 5 isolates susceptible to tetracycline between January 2009 to April 2010 was 0.34 (10/29) while that between May 2010 to August 2011 was 0.62 (16/26). Compared to the period between January 2009 to April 2010, we found that the number of isolates susceptible to tetracycline significantly increased between May 2010 to August 2011 (Pearson chi2 = 0.02, Fisher’s exact test = 0.03). The same was true for oxacillin although the increase in susceptibility was not statistically significant. The proportion of ST 3404 isolates susceptible to tetracycline and oxacillin was 0.83 (25/30) and 0.97 (29/30), while that for ST 3339 was 0.05 (1/19) and 0.37(7/19), respectively. There is strong evidence that the number of isolates susceptible to tetracycline and oxacillin is different for ST 3404 and ST 3339, with the emerging clone ST3404 showing greater susceptibility. All P-values for both Pearson chi2 and Fisher’s exact test were < 0.001.
Resistance to 3 or more classes of antibiotics is commonly referred to as multiple-antibiotic resistance (MDR) . About Sixty-one percent (8/13) of all isolates collected in 2009 were resistant to multiple antibiotics, while 30 % (8/27) of all isolates in 2010 and 7 % (1/15) in 2011 were resistant to multiple antibiotics, implying a decrease in the prevalence of MDR serotype 5 over time. Overall, only one ST 3404 isolate showed MDR, suggesting the drop in MDR over time could be attributed to the emergence of ST 3404. Studies carried out in Canada, Greece, Sweden and Europe as a whole reveal that antibiotic resistance patterns are dependent on local prescribing patterns [25–28]. Drug administration patterns such as dose, route of administration, and duration of antibiotic regimen also influence antibiotic resistance . However, the molecular basis of antibiotic resistance is a major factor to consider since it influences the likelihood of selection. Resistance requiring the introduction of a gene from a donor species is less likely to occur than resistance mediated by a single DNA base change, e.g. pneumococcal resistance to rifampin and fluoroquinolones . Results from our study are similar to those from a study in Canada where decline in MDR 19 F, 6B, and 23 F were documented post PCV-7 . The above study suggested that among other factors, vaccine pressure could influence changes in antibiotic susceptibility patterns.
The emergent strain in this study, ST3404, allegedly showed increased susceptibility to tetracycline, oxacillin, and consequently penicillin. In The Gambia, oxacillin is used as an empirical surrogate to penicillin; likewise it is predictive of reduced penicillin susceptibility . However, caution must be exercised when using oxacillin disc diffusion to predict penicillin susceptibility patterns as predictions can at times be limited by false susceptible and resistance interpretations . Nonetheless, by reason of the magnitude of the difference in susceptibility to oxacillin shown by STs 3404 and 3339 (P < 0.001) and the overall low interpretive error rate of oxacillin disc diffusion [29, 30], ST 3404’s alleged increased susceptibility to penicillin might be genuine. In light of these results, surveillance to determine if these patterns will change after the introduction of PCV-13 will be vital to highlight the role of immunization on susceptibility patterns of various serotype 5 clones (ST) in The Gambia.
Conclusion, implications and future perspectives
There was an emergence of ST 3404 in The Gambia prior to the introduction of the PCV-13. Interestingly, compared to ST 3339, ST 3404 showed lower resistance to tetracycline and oxacillin. Findings from this study provide important background data to assess the potential impact of PCV-13 in routine immunization in The Gambia. Finally, we recommend a comparative genome study on STs 3404, 3339 and 289 to better understand differences antimicrobial sensitivity and virulence, if there be any.
We are very grateful to Medical Research Council; we would also like to acknowledge, Aji konteh, Sambou Suso, The molecular microbiology team at Basse and Fajara, Core sequencing facility at MRC Fajara, The MLST data base hosted at Imperial College London, The pneumococcal surveillance project clinical team and field workers.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Ryan KJ. Streptococci and Enterococci. In: Sherris JC, Ryan KJ, Ray CG (eds) Medical Microbiology: An Introduction to Infectious Diseases. Fourth Edition. McGraw-Hill, USA. 2004;273–296.Google Scholar
- Van der Poll T, Opal SM. Pathogenesis, treatment, and prevention of pneumococcal pneumonia. Lancet. 2009;374(9700):1543–56.View ArticlePubMedGoogle Scholar
- Tettelin H, Nelson KE, Paulsen IT, Eisen JA, Read TD, Peterson S, et al. Complete genome sequence of a virulent isolate of Streptococcus pneumoniae. Science. 2001;293:498–506.View ArticlePubMedGoogle Scholar
- Obert C, Sublett J, Kaushal D, Hinojosa E, Barton T, Tuomanen EI, et al. Identification of a candidate Streptococcus pneumoniae core genome and regions of diversity correlated with invasive pneumococcal disease. Infect Immun. 2006;74:4766–77.View ArticlePubMedPubMed CentralGoogle Scholar
- Antonio M, Dada-Adegbola H, Biney E, Awine T, O’Callaghan J, Pfluger V, et al. Molecular epidemiology of pneumococci obtained from Gambian children aged 2–29 months with invasive pneumococcal disease during a trial of a 9-valent pneumococcal conjugate vaccine. BMC Infect Dis. 2008;8:81.View ArticlePubMedPubMed CentralGoogle Scholar
- Elberse KEM, van der Pol P, Witteveen S, van der Heide J, Schot CS, van Dijk AP, et al. Population structure of invasive streptococcus pneumoniae in the Netherlands in the Pre-vaccination Era assessed by MLVA and capsular sequence typing. PLoS One. 2011;6(5):e20390.View ArticlePubMedPubMed CentralGoogle Scholar
- Usen S, Adegbola R, Mulholland K, Jaffar S, Hilton S, Oparaugo A, et al. Epidemiology of invasive pneumococcal disease in the Western Region, The Gambia. Pediatr Infect Dis J. 1998;17(1):23–8.View ArticlePubMedGoogle Scholar
- Antonio M, Hakeem I, Sankareh K, Cheung YB, Adegbola RA. Evaluation of sequential multiplex PCR for direct detection of multiple serotypes of Streptococcus pneumoniae from nasopharyngeal secretions. J Med Microbiol. 2009;58(Pt 3):296–302.View ArticlePubMedGoogle Scholar
- Donkor ES, Bishop CJ, Gould K, Hinds J, Antonio M, Wren B, et al. High Levels of Recombination among Streptococcus pneumoniae Isolates from The Gambia. mBio. 2011;2(3):e00040–11.View ArticlePubMedPubMed CentralGoogle Scholar
- Ortqvist A, Hedlund J, Kalin M. Streptococcus pneumoniae: epidemiology, risk factors, and clinical features. Semin Respir Crit Care Med. 2005;26:563–74.View ArticlePubMedGoogle Scholar
- Pillai DR, Shahinas D, Buzina A, Pollock RA, Lau R, Khairnar K, et al. Genome-wide dissection of globally emergent multi-drug resistant serotype 19A Streptococcus pneumonia. BMC Genomics. 2009;10:642.View ArticlePubMedPubMed CentralGoogle Scholar
- Regions of Gambia. [http://www.statoids.com/ugm.html]. Accessed 5/01/2016.
- Adegbola R, Hill P, Secka O, Ikumapayi U, Lahai G, Greenwood B, et al. Serotype and antimicrobial susceptibility patterns of isolates of Streptococcus pneumoniae causing invasive disease in The Gambia 1996–2003. Trop Med Int Health. 2006;11:1128–35.View ArticlePubMedGoogle Scholar
- PCR Deduction of Pneumococcal Serotypes. [http://www.cdc.gov/streplab/downloads/triplex-pcr-africa.pdf]. Accessed 5/01/2016.
- Antonio M, Hakeem I, Awine T, Secka O, Sankareh K, Nsekpong D, et al. Seasonality and outbreak of a predominant Streptococcus pneumoniae serotype 1 clone from The Gambia: expansion of ST217 hypervirulent clonal complex in West Africa. BMC Microbiol. 2008;8:198.View ArticlePubMedPubMed CentralGoogle Scholar
- Sanger F, Nicklen S, Coulson AR. DNA sequencing with chain-terminating inhibitors. Proc Natl Acad Sci U S A. 1977;74(12):5463–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Aanensen DM, Spratt BG. The multilocus sequence typing network: mlst.net. Nucleic Acids Res. 2005;33:W728–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Roper DI, Huyton T, Vagin A, Dodson G. The molecular basis of vancomycin resistance in clinically relevant Enterococci: crystal structure of D-alanyl-D-lactate ligase (VanA). Proc Natl Acad Sci U S A. 2000;97(16):8921–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Braud VM, Allan DS, O’Callaghan CA, Söderström K, D’Andrea A, Ogg GS, et al. HLA-E binds to natural killer cell receptors CD94/NKG2A, B and C”. Nature. 1998;391(6669):795–9. 1.View ArticlePubMedGoogle Scholar
- Enright MC, Spratt BG. Extensive variation in the ddl gene of penicillin-resistant Streptococcus pneumoniae results from a hitchhiking effect driven by the penicillin-binding protein 2b gene. Mol Biol Evol. 1999;16(12):1687–95.View ArticlePubMedGoogle Scholar
- Sultan B, Labadi K, Guégan JF, Janicot S. Climate drives the meningitis epidemics onset in West Africa. PLoS Med. 2005;2(1):e6.View ArticlePubMedPubMed CentralGoogle Scholar
- Molesworth A, Cuevas LE, Morse AP, Herman JR, Thomson MC. Dust clouds and spread of infection. Lancet. 2002;359:81–2.View ArticlePubMedGoogle Scholar
- Molesworth AM, Thomson MC, Connor SJ, Cresswell MP, Morse AP, Shears P, et al. Where is the Meningitis Belt? Defining an area at risk of epidemic meningitis in Africa. Trans R Soc Med Hyg. 2002;96:242–9.View ArticleGoogle Scholar
- Klugman KP. Antibiotic selection of multiply resistant pneumococci. Clinical Infectious Disease. 2001;33:489–91.View ArticleGoogle Scholar
- Dunais B, Carsenti-Etesse H, Dellamonica P. Pneumococcal resistance patterns in Europe. Eur J Clin Microbiol Infect Dis. 1997;16:644–7.View ArticlePubMedGoogle Scholar
- Diekema DJ, Brueggemann AB, Doern GV. Antimicrobial-drug use and changes in resistance in Streptococcus pneumoniae. Emerg Infect Dis. 2000;6:552–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Melander E, Molstad S, Persson K, Hansson HB, Soderstrom M, Ekdahl K. Previous antibiotic consumption and other risk factors for carriage of penicillin-resistant Streptococcus pneumoniae in children. Eur J Clin Microbiol Infect Dis. 1998;17:834–8.View ArticlePubMedGoogle Scholar
- Arason VA, Kristinsson KG, Sigurdsson JA, Stefansdottir G, Molstad S, Gudmundsson S. Do antimicrobials increase the carriage rate of penicillin resistant pneumococci in children? Cross sectional prevalence study. BMJ. 1996;313:387–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Manninen R, Leinonen M, Huovinen P, Nissinen A. Reliability of disc diffusion susceptibility testing of Streptococcus pneumoniae and adjustment of laboratory-specific breakpoints. Finnish Study Group for Antimicrobial Resistance. J Antimicrob Chemother. 1998;41(1):19–26.View ArticlePubMedGoogle Scholar
- Jetté LP, Sinave C. Use of an oxacillin disk screening test for detection of penicillin- and ceftriaxone-resistant pneumococci. J Clin Microbiol. 1999;37(4):1178–81.PubMedPubMed CentralGoogle Scholar
- Francisco AP, Bugalho M, Ramirez M, Carriço JA. Global optimal eBURST analysis of multilocus typing data using a graphic matroid approach. BMC Bioinformatics. 2009;10:152.View ArticlePubMedPubMed CentralGoogle Scholar
- Pai R, Gertz RE, Beall B. Sequential multiplex PCR approach for determining capsular serotypes of streptococcus pneumoniae isolates. J Clin Microbiol. 2006;44(1):124–31.View ArticlePubMedPubMed CentralGoogle Scholar