Comparison of qPCR and metagenomic sequencing methods for quantifying antibiotic resistance genes in wastewater.

Surveillance methods of circulating antibiotic resistance genes (ARGs) are of utmost importance in order to tackle what has been described as one of the greatest threats to humanity in the 21st century. In order to be effective, these methods have to be accurate, quickly deployable, and scalable. In this study, we compare metagenomic shotgun sequencing (TruSeq DNA sequencing) of wastewater samples with a state-of-the-art PCR-based method (Resistomap HT-qPCR) on four wastewater samples that were taken from hospital, industrial, urban and rural areas. ARGs that confer resistance to 11 antibiotic classes have been identified in these wastewater samples using both methods, with the most abundant observed classes of ARGs conferring resistance to aminoglycoside, multidrug-resistance (MDR), macrolide-lincosamide-streptogramin B (MLSB), tetracycline and beta-lactams. In comparing the methods, we observed a strong correlation of relative abundance of ARGs obtained by the two tested methods for the majority of antibiotic classes. Finally, we investigated the source of discrepancies in the results obtained by the two methods. This analysis revealed that false negatives were more likely to occur in qPCR due to mutated primer target sites, whereas ARGs with incomplete or low coverage were not detected by the sequencing method due to the parameters set in the bioinformatics pipeline. Indeed, despite the good correlation between the methods, each has its advantages and disadvantages which are also discussed here. By using both methods together, a more robust ARG surveillance program can be established. Overall, the work described here can aid wastewater treatment plants that plan on implementing an ARG surveillance program.

Epidemiology of vancomycin-resistant enterococci in the United Arab Emirates: a retrospective analysis of 12 years of national AMR surveillance data.

Introduction: Enterococci are usually low pathogenic, but can cause invasive disease under certain circumstances, including urinary tract infections, bacteremia, endocarditis, and meningitis, and are associated with peritonitis and intra-abdominal abscesses. Increasing resistance of enterococci to glycopeptides and fluoroquinolones, and high-level resistance to aminoglycosides is a concern. National antimicrobial resistance (AMR) surveillance data for enterococci from the Middle East and North Africa (MENA) and the Gulf region is scarce. Methods: A retrospective 12-year analysis of N = 37,909 non-duplicate diagnostic Enterococcus spp. isolates from the United Arab Emirates (UAE) was conducted. Data was generated by routine patient care during 2010-2021, collected by trained personnel and reported by participating surveillance sites to the UAE National AMR Surveillance program. Data analysis was conducted with WHONET. Results: Enterococcus faecalis was the most commonly reported species (81.5%), followed by Enterococcus faecium (8.5%), and other enterococci species (4.8%). Phenotypically vancomycin-resistant enterococci (VRE) were found in 1.8% of Enterococcus spp. isolates. Prevalence of VRE (%VRE) was highest for E. faecium (8.1%), followed by E. faecalis (0.9%). A significant level of resistance to glycopeptides (%VRE) for these two species has been observed in the majority of observed years [E. faecalis (0-2.2%), 2010: 0%, 2021: 0.6%] and E. faecium (0-14.2%, 2010: 0%, 2021: 5.8%). Resistance to fluoroquinolones was between 17 and 29% (E. faecalis) and was higher for E. faecium (between 42 and 83%). VRE were associated with higher patient mortality (RR: 2.97), admission to intensive care units (RR: 2.25), and increased length of stay (six excess inpatient days per VRE case), as compared to vancomycin-susceptible Enterococcus spp. Discussion: Published data on Enterococcus infections, in particular VRE-infections, in the UAE and MENA region is scarce. Our data demonstrates that VRE-enterococci are relatively rare in the UAE, however showing an increasing resistance trend for several clinically important antibiotic classes, causing a concern for the treatment of serious infections caused by enterococci. This study also demonstrates that VRE were associated with higher mortality, increased intensive care unit admission rates, and longer hospitalization, thus poorer clinical outcome and higher associated costs in the UAE. We recommend the expansion of current surveillance techniques (e.g., local VRE screening), stricter infection prevention and control strategies, and better stewardship interventions. Further studies on the molecular epidemiology of enterococci are needed.

Methicillin resistant Staphylococcus aureus in the United Arab Emirates: a 12-year retrospective analysis of evolving trends.

Introduction: Methicillin resistant Staphylococcus aureus (MRSA) is a major contributor to the global burden of antimicrobial resistance (AMR). As MRSA continues to evolve, the need for continued surveillance to evaluate trends remains crucial. This study was carried out to assess MRSA trends in the United Arab Emirates (UAE) based on analysis of data from the national AMR surveillance program. Methods: We carried out a 12-year (2010-2021) retrospective analysis of MRSA demographic and microbiological data collected as part of the UAE national AMR surveillance program. Participating centers from across the country routinely submit AMR surveillance data collected by trained personnel to the National AMR Surveillance Committee, where data is analyzed using a unified WHONET platform. Data on non-duplicate isolates associated with clinical infections were obtained and included in the analysis. Results: A total of 29,414 non-duplicate MRSA isolates associated with clinical infections were reported between 2010 and 2021 (2010: n = 259; 2021: n = 4,996). MRSA represented 26.4% of all S. aureus (n = 111,623) isolates identified during the study period. In 2010, among the S. aureus isolates with reported oxacillin testing, 21.9% (n/N = 259/1,181) were identified as MRSA and this showed an increase to 33.5% (n/N = 4,996/14,925) in 2021. Although there was variation in the distribution of MRSA across the seven emirates of the country, most had an upward trend. Patient demographics reflected a male preponderance, with most being adults and from the outpatient setting. Isolates were mostly from skin and soft tissue infection specimens (72.5%; n/N = 21,335/29,414). Among the inpatients (N = 8,282), a total of 3,313 MRSA isolates were from specimens obtained ≤ 48 h after admission indicative of community acquired infection. Increasing resistance trends were observed for most antibiotics including ciprofloxacin, levofloxacin, moxifloxacin, erythromycin, gentamicin, trimethoprim-sulfamethoxazole, and quinupristin/dalfopristin. Low levels of resistance (0.0-0.8%) were sustained for linezolid except for 2015, 2016, and 2017 with 2.5, 2.6, and 2.9%, respectively. No confirmed vancomycin resistance was reported. Conclusion: The increasing trend of MRSA isolates associated with clinical infections in the hospital and community settings is a concern. Continued monitoring including incorporation of genomic surveillance and infection control measures are recommended to stem the dissemination.

Evolving trends among Pseudomonas aeruginosa: a 12-year retrospective study from the United Arab Emirates.

Introduction: Pseudomonas is a group of ubiquitous non-fermenting Gram-negative bacteria (NFGNB). Of the several species associated with humans, Pseudomonas aeruginosa (PA) can acclimate to diverse environments. The global frequency of PA infections is rising and is complicated by this organism's high intrinsic and acquired resistance to several clinically relevant antibiotics. Data on the epidemiology, levels, and trends of antimicrobial resistance of PA in clinical settings in the MENA/GCC region is scarce. Methods: A retrospective 12-year analysis of 56,618 non-duplicate diagnostic Pseudomonas spp. from the United Arab Emirates (UAE) was conducted. Data was generated at 317 surveillance sites by routine patient care during 2010-2021, collected by trained personnel and reported by participating surveillance sites to the UAE National antimicrobial resistance (AMR) Surveillance program. Data analysis was conducted with WHONET (https://whonet.org/). Results: Among the total isolates (N = 56,618), the majority were PA (95.6%). Data on nationality revealed 44.1% were UAE nationals. Most isolates were from soft tissue (55.7%), followed by respiratory tract (26.7%). PA was more commonly found among inpatients than among outpatients, followed by ICUs. PA showed a horizontal trend for resistance to fluoroquinolones, 3rd- and 4th-generation cephalosporins, and decreasing trends of resistance for aminoglycosides and meropenem. The highest percentage of multidrug resistant (MDR) isolates was reported in 2011 at 35.6%. As an overall trend, the percentage of MDR, extensively drug-resistant (XDR), and possible pandrug-resistant (PDR) isolates generally declined over the study period. Carbapenem-resistant PA (CRPA) were associated with a higher mortality (RR: 2.7), increased admission to ICU (RR: 2.3), and increased length of stay (LOS) (12 excess inpatient days per case), as compared to carbapenem-susceptible PA (CSPA). Conclusion: The resistance trends in Pseudomonas species in the UAE indicated a decline in AMR and in percentages of Pseudomonas isolates with MDR and XDR profiles. The sustained Pseudomonas spp. circulation particularly in the hospital settings highlights the importance of surveillance techniques, infection control strategies, and stewardship to limit the continued dissemination. This data also shows that CRPA are associated with higher mortality, increased ICU admission rates, and a longer hospitalization, thus higher costs due to increased number of in-hospital and ICU days.

Carbapenem resistant Enterobacterales in the United Arab Emirates: a retrospective analysis from 2010 to 2021.

Background: Carbapenem-resistant Enterobacterales (CRE) are spreading in the United Arab Emirates (UAE) where their dissemination is facilitated by international travel, trade, and tourism. The objective of this study is to describe the longitudinal changes of CRE as reported by the national AMR surveillance system of the UAE. Methods: In this study, we retrospectively describe CRE isolated from 317 surveillance sites, including 87 hospitals and 230 centers/clinics from 2010 to 2021. The associated clinical, demographic, and microbiological characteristics are presented by relying on the UAE national AMR surveillance program. Data was analyzed using WHONET microbiology laboratory database software (http://www.whonet.org). Results: A total of 14,593 carbapenem resistant Enterobacterales were analyzed, of which 48.1% were carbapenem resistant Klebsiella pneumoniae (CRKp), 25.1% carbapenem resistant Escherichia coli (CREc), and 26.8% represented 72 other carbapenem resistant species. Carbapenem resistant strains were mostly associated with adults and isolated from urine samples (36.9% of CRKp and 66.6% of CREc) followed by respiratory samples (26.95% for CRKp) and soft tissue samples (19.5% for CRKp). Over the studied period carbapenem resistance rates remained high, especially in K. pneumoniae, and in 2021 were equivalent to 67.6% for imipenem, 76.2% for meropenem, and 91.6% for ertapenem. Nevertheless, there was a statistically significant decreasing trend for imipenem and meropenem resistance in Klebsiella species (p < 0.01) while the decrease in ertapenem resistance was non-significant. Concerning E. coli, there was a statistically significant decreasing trend for meropenem and imipenem resistance over the 12 years, while ertapenem resistance increased significantly with 83.8% of E. coli exhibiting ertapenem resistance in 2021. Resistance rates to ceftazidime and cefotaxime remained higher than 90% (in 2021) for CRKp and cefotaxime rates increased to 90.5% in 2021 for CREc. Starting 2014, resistance to colistin and tigecycline was observed in carbapenem resistant Enterobacterales. CRE were associated with a higher mortality (RR: 6.3), admission to ICU (RR 3.9), and increased length of stay (LOS; 10 excess inpatient days per CRE case). Conclusion: This study supports the need to monitor CRE in the UAE and draws attention to the significant increase of ertapenem resistance in E. coli. Future surveillance analysis should include a genetic description of carbapenem resistance to provide new strategies.

Epidemiology and antimicrobial resistance trends of Acinetobacter species in the United Arab Emirates: a retrospective analysis of 12 years of national AMR surveillance data.

Introduction: Acinetobacter spp., in particular A. baumannii, are opportunistic pathogens linked to nosocomial pneumonia (particularly ventilator-associated pneumonia), central-line catheter-associated blood stream infections, meningitis, urinary tract infections, surgical-site infections, and other types of wound infections. A. baumannii is able to acquire or upregulate various resistance determinants, making it frequently multidrug-resistant, and contributing to increased mortality and morbidity. Data on the epidemiology, levels, and trends of antimicrobial resistance of Acinetobacter spp. in clinical settings is scarce in the Gulf Cooperation Council (GCC) and Middle East and North Africa (MENA) regions. Methods: A retrospective 12-year analysis of 17,564 non-duplicate diagnostic Acinetobacter spp. isolates from the United Arab Emirates (UAE) was conducted. Data was generated at 317 surveillance sites by routine patient care during 2010-2021, collected by trained personnel and reported by participating surveillance sites to the UAE National AMR Surveillance program. Data analysis was conducted with WHONET. Results: Species belonging to the A. calcoaceticus-baumannii complex were mostly reported (86.7%). They were most commonly isolated from urine (32.9%), sputum (29.0%), and soft tissue (25.1%). Resistance trends to antibiotics from different classes during the surveillance period showed a decreasing trend. Specifically, there was a significant decrease in resistance to imipenem, meropenem, and amikacin. Resistance was lowest among Acinetobacter species to both colistin and tigecycline. The percentages of multidrug-resistant (MDR) and possibly extensively drug-resistant (XDR) isolates was reduced by almost half between the beginning of the study in 2010 and its culmination in 2021. Carbapenem-resistant Acinetobacter spp. (CRAB) was associated with a higher mortality (RR: 5.7), a higher admission to ICU (RR 3.3), and an increased length of stay (LOS; 13 excess inpatient days per CRAB case), as compared to Carbapenem-susceptible Acinetobacter spp. Conclusion: Carbapenem-resistant Acinetobacter spp. are associated with poorer clinical outcomes, and higher associated costs, as compared to carbapenem-susceptible Acinetobacter spp. A decreasing trend of MDR Acinetobacter spp., as well as resistance to all antibiotic classes under surveillance was observed during 2010 to 2021. Further studies are needed to explore the reasons and underlying factors leading to this remarkable decrease of resistance over time.

Emergence of highly resistant Candida auris in the United Arab Emirates: a retrospective analysis of evolving national trends.

Introduction: The Centers for Disease Prevention and Control lists Candida auris, given its global emergence, multidrug resistance, high mortality, and persistent transmissions in health care settings as one of five urgent threats. As a new threat, the need for surveillance of C. auris is critical. This is particularly important for a cosmopolitan setting and global hub such as the United Arab Emirates (UAE) where continued introduction and emergence of resistant variant strains is a major concern. Methods: The United Arab Emirates has carried out a 12 years of antimicrobial resistance surveillance (2010-2021) across the country, spanning all seven Emirates. A retrospective analysis of C. auris emergence from 2018-2021 was undertaken, utilising the demographic and microbiological data collected via a unified WHONET platform for AMR surveillance. Results: Nine hundred eight non-duplicate C. auris isolates were reported from 2018-2021. An exponential upward trend of cases was found. Most isolates were isolated from urine, blood, skin and soft tissue, and the respiratory tract. UAE nationals nationals comprised 29% (n = 186 of 632) of all patients; the remainder were from 34 other nations. Almost all isolates were from inpatient settings (89.0%, n = 809). The cases show widespread distribution across all reporting sites in the country. C. auris resistance levels remained consistently high across all classes of antifungals used. C. auris in this population remains highly resistant to azoles (fluconazole, 72.6% in 2021) and amphotericin. Echinocandin resistance has now emerged and is increasing annually. There was no statistically significant difference in mortality between Candida auris and Candida spp. (non-auris) patients (p-value: 0.8179), however Candida auris patients had a higher intensive care unit (ICU) admission rate (p-value <0.0001) and longer hospital stay (p < 0.0001) compared to Candida spp. (non-auris) patients. Conclusion: The increasing trend of C. auris detection and associated multidrug resistant phenotypes in the UAE is alarming. Continued C. auris circulation in hospitals requires enhanced infection control measures to prevent continued dissemination.

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

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

Comparative Genomics of Disease and Carriage Serotype 1 Pneumococci.

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

Impact and effectiveness of 13-valent pneumococcal conjugate vaccine on population incidence of vaccine and non-vaccine serotype invasive pneumococcal disease in Blantyre, Malawi, 2006-18: prospective observational time-series and case-control studies.

BACKGROUND: The population impact of pneumococcal conjugate vaccines (PCVs) depends on direct and indirect protection. Following Malawi's introduction of the 13-valent PCV (PCV13) in 2011, we examined its impact on vaccine and non-vaccine serotype invasive pneumococcal disease among vaccine-eligible-age and vaccine-ineligible-age children and adults. METHODS: We did a prospective observational time-series analysis and a case-control study. We used data from between Jan 1, 2006, and Dec 31, 2018, from laboratory-based surveillance at a government hospital in Malawi. This period included 6 years before and 7 years after introduction of PCV13. By use of negative-binomial regression, we evaluated secular trend-adjusted incidence rate ratio (IRR) in vaccine serotype and non-vaccine serotype invasive pneumococcal disease before and after introduction of PCV. We compared predicted counterfactual incidence in hypothetical absence of vaccine with empirically observed incidence following vaccine introduction. The case-control study assessed vaccine effectiveness, comparing PCV uptake among cases of vaccine-eligible-age invasive pneumococcal disease versus matched community controls. FINDINGS: Surveillance covered 10 281 476 person-years of observation, with 140 498 blood and 63 291 cerebrospinal fluid cultures. A reduction in total (vaccine serotype plus non-vaccine serotype) invasive pneumococcal disease incidence preceded introduction of PCV: 19% (IRR 0·81, 95% CI 0·74 to 0·88, p<0·0001) among infants (<1 year old), 14% (0·86, 0·80 to 0·93, p<0·0001) among children aged 1-4 years, and 8% (0·92, 0·83 to 1·01, p=0·084) among adolescents and adults (≥15 years old). Among children aged 5-14 years there was a 2% increase in total invasive pneumococcal disease (1·02, 0·93 to 1·11, p=0·72). Compared with the counterfactually predicted incidence, incidence of post-PCV13 vaccine serotype invasive pneumococcal disease was 74% (95% CI 70 to 78) lower among children aged 1-4 years and 79% (76 to 83) lower among children aged 5-14 years, but only 38% (37 to 40) lower among infants and 47% (44 to 51) lower among adolescents and adults. Although non-vaccine serotype invasive pneumococcal disease has increased in incidence since 2015, observed incidence remains low. The case-control study (19 cases and 76 controls) showed vaccine effectiveness against vaccine serotype invasive pneumococcal disease of 80·7% (-73·7 to 97·9). INTERPRETATION: In a high-mortality, high-HIV-prevalence setting in Africa, there were significant pre-vaccine reductions in the incidence of invasive pneumococcal disease. 7 years after PCV introduction, although vaccine-attributable impact among vaccine-eligible-age children was significant, indirect effects benefitting unvaccinated infants and adults were not. Policy decisions should consider multiple alternative strategies for reducing disease burden, including targeted vaccination outside infant Expanded Programme of Immunization to benefit vulnerable populations. FUNDING: Bill & Melinda Gates Foundation, Wellcome Trust, and National Institute for Health Research.

Characterization of DNA methylation in Malawian Mycobacterium tuberculosis clinical isolates.

BACKGROUND: Although Mycobacterium tuberculosis (Mtb) strains exhibit genomic homology of >99%, there is considerable variation in the phenotype. The underlying mechanisms of phenotypic heterogeneity in Mtb are not well understood but epigenetic variation is thought to contribute. At present the methylome of Mtb has not been completely characterized. METHODS: We completed methylomes of 18 Mycobacterium tuberculosis (Mtb) clinical isolates from Malawi representing the largest number of Mtb genomes to be completed in a single study using Single Molecule Real Time (SMRT) sequencing to date. RESULTS: We replicate and confirm four methylation disrupting mutations in 4 lineages of Mtb. For the first time we report complete loss of methylation courtesy of C758T (S253L) mutation in the MamB gene of Indo-oceanic lineage of Mtb. Additionally, we report a novel missense mutation G454A (G152S) in the MamA gene of the Euro-American lineage which could potentially be attributed to total disruption of methylation in the CCCAG motif but partial loss in a partner motif. Through a genomic and methylome comparative analysis with a global sample of sixteen, we report previously unknown mutations affecting the pks15/1 locus in L6 isolates. We confirm that methylation in Mtb is lineage specific although some unresolved issues still remain.

Lower Density and Shorter Duration of Nasopharyngeal Carriage by Pneumococcal Serotype 1 (ST217) May Explain Its Increased Invasiveness over Other Serotypes.

Streptococcus pneumoniae is a frequent colonizer of the human nasopharynx and a major cause of life-threating invasive infections such as pneumonia, meningitis and sepsis. Over 1 million people die every year due to invasive pneumococcal disease (IPD), mainly in developing countries. Serotype 1 is a common cause of IPD; however, unlike other serotypes, it is rarely found in the carrier state in the nasopharynx, which is often considered a prerequisite for disease. The aim of this study was to understand this dichotomy. We used murine models of carriage and IPD to characterize the pathogenesis of African serotype 1 (sequence type 217) pneumococcal strains obtained from the Queen Elizabeth Central Hospital in Blantyre, Malawi. We found that ST217 pneumococcal strains were highly virulent in a mouse model of invasive pneumonia, but in contrast to the generally accepted assumption, can also successfully establish nasopharyngeal carriage. Interestingly, we found that cocolonizing serotypes may proliferate in the presence of serotype 1, suggesting that acquisition of serotype 1 carriage could increase the risk of developing IPD by other serotypes. RNA sequencing analysis confirmed that key virulence genes associated with inflammation and tissue invasiveness were upregulated in serotype 1. These data reveal important new insights into serotype 1 pathogenesis, with implications for carriage potential and risk of invasive disease through interactions with other cocolonizing serotypes, an often-overlooked factor in transmission and disease progression.IMPORTANCE The pneumococcus causes serious diseases such as pneumonia, sepsis, and meningitis and is a major cause of morbidity and mortality worldwide. Serotype 1 accounts for the majority of invasive pneumococcal disease cases in sub-Saharan Africa but is rarely found during nasopharyngeal carriage. Understanding the mechanisms leading to nasopharyngeal carriage and invasive disease by this serotype can help reduce its burden on health care systems worldwide. In this study, we also uncovered the potential impact of serotype 1 on disease progression of other coinfecting serotypes, which can have important implications for vaccine efficacy. Understanding the interactions between different serotypes during nasopharyngeal carriage may lead to improved intervention methods and therapies to reduce pneumococcal invasive disease levels.

Bacterial genome-wide association study of hyper-virulent pneumococcal serotype 1 identifies genetic variation associated with neurotropism.

Hyper-virulent Streptococcus pneumoniae serotype 1 strains are endemic in Sub-Saharan Africa and frequently cause lethal meningitis outbreaks. It remains unknown whether genetic variation in serotype 1 strains modulates tropism into cerebrospinal fluid to cause central nervous system (CNS) infections, particularly meningitis. Here, we address this question through a large-scale linear mixed model genome-wide association study of 909 African pneumococcal serotype 1 isolates collected from CNS and non-CNS human samples. By controlling for host age, geography, and strain population structure, we identify genome-wide statistically significant genotype-phenotype associations in surface-exposed choline-binding (P = 5.00 × 10-08) and helicase proteins (P = 1.32 × 10-06) important for invasion, immune evasion and pneumococcal tropism to CNS. The small effect sizes and negligible heritability indicated that causation of CNS infection requires multiple genetic and other factors reflecting a complex and polygenic aetiology. Our findings suggest that certain pathogen genetic variation modulate pneumococcal survival and tropism to CNS tissue, and therefore, virulence for meningitis.

Visualizing variation within Global Pneumococcal Sequence Clusters (GPSCs) and country population snapshots to contextualize pneumococcal isolates.

Knowledge of pneumococcal lineages, their geographic distribution and antibiotic resistance patterns, can give insights into global pneumococcal disease. We provide interactive bioinformatic outputs to explore such topics, aiming to increase dissemination of genomic insights to the wider community, without the need for specialist training. We prepared 12 country-specific phylogenetic snapshots, and international phylogenetic snapshots of 73 common Global Pneumococcal Sequence Clusters (GPSCs) previously defined using PopPUNK, and present them in Microreact. Gene presence and absence defined using Roary, and recombination profiles derived from Gubbins are presented in Phandango for each GPSC. Temporal phylogenetic signal was assessed for each GPSC using BactDating. We provide examples of how such resources can be used. In our example use of a country-specific phylogenetic snapshot we determined that serotype 14 was observed in nine unrelated genetic backgrounds in South Africa. The international phylogenetic snapshot of GPSC9, in which most serotype 14 isolates from South Africa were observed, highlights that there were three independent sub-clusters represented by South African serotype 14 isolates. We estimated from the GPSC9-dated tree that the sub-clusters were each established in South Africa during the 1980s. We show how recombination plots allowed the identification of a 20 kb recombination spanning the capsular polysaccharide locus within GPSC97. This was consistent with a switch from serotype 6A to 19A estimated to have occured in the 1990s from the GPSC97-dated tree. Plots of gene presence/absence of resistance genes (tet, erm, cat) across the GPSC23 phylogeny were consistent with acquisition of a composite transposon. We estimated from the GPSC23-dated tree that the acquisition occurred between 1953 and 1975. Finally, we demonstrate the assignment of GPSC31 to 17 externally generated pneumococcal serotype 1 assemblies from Utah via Pathogenwatch. Most of the Utah isolates clustered within GPSC31 in a USA-specific clade with the most recent common ancestor estimated between 1958 and 1981. The resources we have provided can be used to explore to data, test hypothesis and generate new hypotheses. The accessible assignment of GPSCs allows others to contextualize their own collections beyond the data presented here.

Increased pathogenicity of pneumococcal serotype 1 is driven by rapid autolysis and release of pneumolysin.

Streptococcus pneumoniae serotype 1 is the predominant cause of invasive pneumococcal disease in sub-Saharan Africa, but the mechanism behind its increased invasiveness is not well understood. Here, we use mouse models of lung infection to identify virulence factors associated with severe bacteraemic pneumonia during serotype-1 (ST217) infection. We use BALB/c mice, which are highly resistant to pneumococcal pneumonia when infected with other serotypes. However, we observe 100% mortality and high levels of bacteraemia within 24 hours when BALB/c mice are intranasally infected with ST217. Serotype 1 produces large quantities of pneumolysin, which is rapidly released due to high levels of bacterial autolysis. This leads to substantial levels of cellular cytotoxicity and breakdown of tight junctions between cells, allowing a route for rapid bacterial dissemination from the respiratory tract into the blood. Thus, our results offer an explanation for the increased invasiveness of serotype 1.

Early Signals of Vaccine-driven Perturbation Seen in Pneumococcal Carriage Population Genomic Data.

BACKGROUND: Pneumococcal conjugate vaccines (PCVs) have reduced pneumococcal diseases globally. Pneumococcal genomic surveys elucidate PCV effects on population structure but are rarely conducted in low-income settings despite the high disease burden. METHODS: We undertook whole-genome sequencing (WGS) of 660 pneumococcal isolates collected through surveys from healthy carriers 2 years from 13-valent PCV (PCV13) introduction and 1 year after rollout in northern Malawi. We investigated changes in population structure, within-lineage serotype dynamics, serotype diversity, and frequency of antibiotic resistance (ABR) and accessory genes. RESULTS: In children <5 years of age, frequency and diversity of vaccine serotypes (VTs) decreased significantly post-PCV, but no significant changes occurred in persons ≥5 years of age. Clearance of VT serotypes was consistent across different genetic backgrounds (lineages). There was an increase of nonvaccine serotypes (NVTs)-namely 7C, 15B/C, and 23A-in children <5 years of age, but 28F increased in both age groups. While carriage rates have been recently shown to remain stable post-PCV due to replacement serotypes, there was no change in diversity of NVTs. Additionally, frequency of intermediate-penicillin-resistant lineages decreased post-PCV. Although frequency of ABR genes remained stable, other accessory genes, especially those associated with mobile genetic element and bacteriocins, showed changes in frequency post-PCV. CONCLUSIONS: We demonstrate evidence of significant population restructuring post-PCV driven by decreasing frequency of vaccine serotypes and increasing frequency of few NVTs mainly in children under 5. Continued surveillance with WGS remains crucial to fully understand dynamics of the residual VTs and replacement NVT serotypes post-PCV.

A mosaic tetracycline resistance gene tet(S/M) detected in an MDR pneumococcal CC230 lineage that underwent capsular switching in South Africa.

OBJECTIVES: We reported tet(S/M) in Streptococcus pneumoniae and investigated its temporal spread in relation to nationwide clinical interventions. METHODS: We whole-genome sequenced 12 254 pneumococcal isolates from 29 countries on an Illumina HiSeq sequencer. Serotype, multilocus ST and antibiotic resistance were inferred from genomes. An SNP tree was built using Gubbins. Temporal spread was reconstructed using a birth-death model. RESULTS: We identified tet(S/M) in 131 pneumococcal isolates and none carried other known tet genes. Tetracycline susceptibility testing results were available for 121 tet(S/M)-positive isolates and all were resistant. A majority (74%) of tet(S/M)-positive isolates were from South Africa and caused invasive diseases among young children (59% HIV positive, where HIV status was available). All but two tet(S/M)-positive isolates belonged to clonal complex (CC) 230. A global phylogeny of CC230 (n=389) revealed that tet(S/M)-positive isolates formed a sublineage predicted to exhibit resistance to penicillin, co-trimoxazole, erythromycin and tetracycline. The birth-death model detected an unrecognized outbreak of this sublineage in South Africa between 2000 and 2004 with expected secondary infections (effective reproductive number, R) of ∼2.5. R declined to ∼1.0 in 2005 and <1.0 in 2012. The declining epidemic could be related to improved access to ART in 2004 and introduction of pneumococcal conjugate vaccine (PCV) in 2009. Capsular switching from vaccine serotype 14 to non-vaccine serotype 23A was observed within the sublineage. CONCLUSIONS: The prevalence of tet(S/M) in pneumococci was low and its dissemination was due to an unrecognized outbreak of CC230 in South Africa. Capsular switching in this MDR sublineage highlighted its potential to continue to cause disease in the post-PCV13 era.

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.

BACKGROUND: Invasive pneumococcal disease remains an important health priority owing to increasing disease incidence caused by pneumococci expressing non-vaccine serotypes. We previously defined 621 Global Pneumococcal Sequence Clusters (GPSCs) by analysing 20 027 pneumococcal isolates collected worldwide and from previously published genomic data. In this study, we aimed to investigate the pneumococcal lineages behind the predominant serotypes, the mechanism of serotype replacement in disease, as well as the major pneumococcal lineages contributing to invasive pneumococcal disease in the post-vaccine era and their antibiotic resistant traits. METHODS: We whole-genome sequenced 3233 invasive pneumococcal disease isolates from laboratory-based surveillance programmes in Hong Kong (n=78), Israel (n=701), Malawi (n=226), South Africa (n=1351), The Gambia (n=203), and the USA (n=674). The genomes represented pneumococci from before and after pneumococcal conjugate vaccine (PCV) introductions and were from children younger than 3 years. We identified predominant serotypes by prevalence and their major contributing lineages in each country, and assessed any serotype replacement by comparing the incidence rate between the pre-PCV and PCV periods for Israel, South Africa, and the USA. We defined the status of a lineage as vaccine-type GPSC (≥50% 13-valent PCV [PCV13] serotypes) or non-vaccine-type GPSC (>50% non-PCV13 serotypes) on the basis of its initial serotype composition detected in the earliest vaccine period to measure their individual contribution toward serotype replacement in each country. Major pneumococcal lineages in the PCV period were identified by pooled incidence rate using a random effects model. FINDINGS: The five most prevalent serotypes in the PCV13 period varied between countries, with only serotypes 5, 12F, 15B/C, 19A, 33F, and 35B/D common to two or more countries. The five most prevalent serotypes in the PCV13 period varied between countries, with only serotypes 5, 12F, 15B/C, 19A, 33F, and 35B/D common to two or more countries. These serotypes were associated with more than one lineage, except for serotype 5 (GPSC8). Serotype replacement was mainly mediated by expansion of non-vaccine serotypes within vaccine-type GPSCs and, to a lesser extent, by increases in non-vaccine-type GPSCs. A globally spreading lineage, GPSC3, expressing invasive serotypes 8 in South Africa and 33F in the USA and Israel, was the most common lineage causing non-vaccine serotype invasive pneumococcal disease in the PCV13 period. We observed that same prevalent non-vaccine serotypes could be associated with distinctive lineages in different countries, which exhibited dissimilar antibiotic resistance profiles. In non-vaccine serotype isolates, we detected significant increases in the prevalence of resistance to penicillin (52 [21%] of 249 vs 169 [29%] of 575, p=0·0016) and erythromycin (three [1%] of 249 vs 65 [11%] of 575, p=0·0031) in the PCV13 period compared with the pre-PCV period. INTERPRETATION: Globally spreading lineages expressing invasive serotypes have an important role in serotype replacement, and emerging non-vaccine serotypes associated with different pneumococcal lineages in different countries might be explained by local antibiotic-selective pressures. Continued genomic surveillance of the dynamics of the pneumococcal population with increased geographical representation in the post-vaccine period will generate further knowledge for optimising future vaccine design. FUNDING: Bill & Melinda Gates Foundation, Wellcome Sanger Institute, and the US Centers for Disease Control.

International genomic definition of pneumococcal lineages, to contextualise disease, antibiotic resistance and vaccine impact.

BACKGROUND: Pneumococcal conjugate vaccines have reduced the incidence of invasive pneumococcal disease, caused by vaccine serotypes, but non-vaccine-serotypes remain a concern. We used whole genome sequencing to study pneumococcal serotype, antibiotic resistance and invasiveness, in the context of genetic background. METHODS: Our dataset of 13,454 genomes, combined with four published genomic datasets, represented Africa (40%), Asia (25%), Europe (19%), North America (12%), and South America (5%). These 20,027 pneumococcal genomes were clustered into lineages using PopPUNK, and named Global Pneumococcal Sequence Clusters (GPSCs). From our dataset, we additionally derived serotype and sequence type, and predicted antibiotic sensitivity. We then measured invasiveness using odds ratios that relating prevalence in invasive pneumococcal disease to carriage. FINDINGS: The combined collections (n = 20,027) were clustered into 621 GPSCs. Thirty-five GPSCs observed in our dataset were represented by >100 isolates, and subsequently classed as dominant-GPSCs. In 22/35 (63%) of dominant-GPSCs both non-vaccine serotypes and vaccine serotypes were observed in the years up until, and including, the first year of pneumococcal conjugate vaccine introduction. Penicillin and multidrug resistance were higher (p < .05) in a subset dominant-GPSCs (14/35, 9/35 respectively), and resistance to an increasing number of antibiotic classes was associated with increased recombination (R2 = 0.27 p < .0001). In 28/35 dominant-GPSCs, the country of isolation was a significant predictor (p < .05) of its antibiogram (mean misclassification error 0.28, SD ±â€¯0.13). We detected increased invasiveness of six genetic backgrounds, when compared to other genetic backgrounds expressing the same serotype. Up to 1.6-fold changes in invasiveness odds ratio were observed. INTERPRETATION: We define GPSCs that can be assigned to any pneumococcal genomic dataset, to aid international comparisons. Existing non-vaccine-serotypes in most GPSCs preclude the removal of these lineages by pneumococcal conjugate vaccines; leaving potential for serotype replacement. A subset of GPSCs have increased resistance, and/or serotype-independent invasiveness.

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

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