Overcoming Data Bottlenecks in Genomic Pathogen Surveillance.

Performing whole genome sequencing (WGS) for the surveillance of antimicrobial resistance offers the ability to determine not only the antimicrobials to which rates of resistance are increasing, but also the evolutionary mechanisms and transmission routes responsible for the increase at local, national, and global scales. To derive WGS-based outputs, a series of processes are required, beginning with sample and metadata collection, followed by nucleic acid extraction, library preparation, sequencing, and analysis. Throughout this pathway there are many data-related operations required (informatics) combined with more biologically focused procedures (bioinformatics). For a laboratory aiming to implement pathogen genomics, the informatics and bioinformatics activities can be a barrier to starting on the journey; for a laboratory that has already started, these activities may become overwhelming. Here we describe these data bottlenecks and how they have been addressed in laboratories in India, Colombia, Nigeria, and the Philippines, as part of the National Institute for Health Research Global Health Research Unit on Genomic Surveillance of Antimicrobial Resistance. The approaches taken include the use of reproducible data parsing pipelines and genome sequence analysis workflows, using technologies such as Data-flo, the Nextflow workflow manager, and containerization of software dependencies. By overcoming barriers to WGS implementation in countries where genome sampling for some species may be underrepresented, a body of evidence can be built to determine the concordance of antimicrobial sensitivity testing and genome-derived resistance, and novel high-risk clones and unknown mechanisms of resistance can be discovered.

Genomic surveillance of Salmonella spp. in the Philippines during 2013-2014.

BACKGROUND: Increasing antimicrobial resistance (AMR) in Salmonella has been observed in the Philippines. We aimed to characterise the population and AMR mechanisms of Salmonella with whole genome sequencing (WGS) and compare it with laboratory surveillance methods. METHODS: The serotype, multilocus sequence type, AMR genes and relatedness between isolates were determined from the genomes of 148 Salmonella Typhi (S. Typhi) and 65 non-typhoidal Salmonella (NTS) collected by the Antimicrobial Resistance Surveillance Program during 2013-2014. Genotypic serotypes and AMR prediction were compared with phenotypic data. RESULTS: AMR rates in S. Typhi were low, with sparse acquisition of mutations associated with reduced susceptibility to fluoroquinolones or extended-spectrum beta-lactamases (ESBL) genes. By contrast, 75% of NTS isolates were insusceptible to at least one antimicrobial, with more than half carrying mutations and/or genes linked to fluoroquinolone resistance. ESBL genes were detected in five genomes, which also carried other AMR determinants. The population of S. Typhi was dominated by likely endemic genotype 3.0, which caused a putative local outbreak. The main NTS clades were global epidemic S. Enteritidis ST11 and S. Typhimurium monophasic variant (I,4,[5],12: i: -) ST34. CONCLUSION: We provide the first genomic characterisation of Salmonella from the Philippines and evidence of WGS utility for ongoing surveillance.

Genomic surveillance of methicillin-resistant Staphylococcus aureus in the Philippines, 2013-2014.

Methicillin-resistant Staphylococcus aureus (MRSA) remains one of the leading causes of both nosocomial and community infections worldwide. In the Philippines, MRSA rates have remained above 50% since 2010, but resistance to other antibiotics, including vancomycin, is low. The MRSA burden can be partially attributed to pathogen-specific characteristics of the circulating clones, but little was known about the S. aureus clones circulating in the Philippines. We sequenced the whole genomes of 116 S. aureus isolates collected in 2013-2014 within the Antimicrobial Resistance Surveillance Program. The multilocus sequence type, spa type, SCCmec type, presence of antimicrobial resistance (AMR) determinants and virulence genes and relatedness between the isolates were all derived from the sequence data. The concordance between phenotypic and genotypic resistance was also determined. The MRSA population in the Philippines comprised a limited number of genetic clones, including several international epidemic clones, such as CC30-spa-t019-SCCmec-IV-PVL+, CC5-SCCmec-typeIV and ST239-spa-t030-SCCmec-typeIII. The CC30 genomes were related to the South-West Pacific clone but formed a distinct, diverse lineage, with evidence of global dissemination. We showed independent acquisition of resistance to sulfamethoxazole/trimethoprim in various locations and genetic clones but mostly in paediatric patients with invasive infections. The concordance between phenotypic and genotypic resistance was 99.68% overall for eight antibiotics in seven classes. We have made the first comprehensive genomic survey of S. aureus in the Philippines, which bridges the gap in genomic data from the Western Pacific Region and will constitute the genetic background for contextualizing prospective surveillance.

Genomic surveillance of Neisseria gonorrhoeae in the Philippines, 2013-2014.

Antimicrobial-resistant Neisseria gonorrhoeae is a major threat to public health and is of particular concern in the Western Pacific Region, where the incidence of gonorrhoea is high. The Antimicrobial Resistance Surveillance Program (ARSP) has been capturing information on resistant gonorrhoea since 1996, but genomic epidemiology studies on this pathogen are lacking in the Philippines. We sequenced the whole genomes of 21 N. gonorrhoeae isolates collected in 2013-2014 by ARSP. The multilocus sequence type, multiantigen sequence type, presence of determinants of antimicrobial resistance and relatedness among the isolates were all derived from the sequence data. The concordance between phenotypic and genotypic resistance was also determined. Ten of 21 isolates were resistant to penicillin, ciprofloxacin and tetracycline, due mainly to the presence of the blaTEM gene, the S91F mutation in the gyrA gene and the tetM gene, respectively. None of the isolates was resistant to ceftriaxone or cefixime. The concordance between phenotypic and genotypic resistance was 92.38% overall for five antibiotics in four classes. Despite the small number of isolates studied, they were genetically diverse, as shown by the sequence types, the N. gonorrhoeae multiantigen sequence typing types and the tree. Comparison with global genomes placed the Philippine genomes within global lineage A and led to the identification of an international transmission route. This first genomic survey of N. gonorrhoeae isolates collected by ARSP will be used to contextualize prospective surveillance. It highlights the importance of genomic surveillance in the Western Pacific and other endemic regions for understanding the spread of drug-resistant gonorrhoea worldwide.

Genomic surveillance of Acinetobacter baumannii in the Philippines, 2013-2014.

OBJECTIVE: Acinetobacter baumannii is an opportunistic nosocomial pathogen that has increasingly become resistant to carbapenems worldwide. In the Philippines, rates of carbapenem resistance and multidrug resistance are above 50%. We undertook a genomic study of carbapenem-resistant A. baumannii in the Philippines to characterize the population diversity and antimicrobial resistance mechanisms. METHODS: We sequenced the whole genomes of 117 A. baumannii isolates recovered by 16 hospitals in the Philippines between 2013 and 2014. From the genome sequences, we determined the multilocus sequence type, presence of acquired determinants of antimicrobial resistance and relatedness between isolates. We also compared the phenotypic and genotypic resistance results. RESULTS: Carbapenem resistance was mainly explained by acquisition of the class-D ß-lactamase gene blaOXA-23. The concordance between phenotypic and genotypic resistance to imipenem was 98.15%, and it was 94.97% overall for the seven antibiotics analysed. Twenty-two different sequence types were identified, including 7 novel types. The population was dominated by the high-risk international clone 2 (i.e. clonal complex 92), in particular by ST195 and ST208 and their single locus variants. Using whole-genome sequencing, we identified local clusters representing potentially undetected nosocomial outbreaks, as well as multihospital clusters that indicated interhospital dissemination. Comparison with global genomes suggested that the establishment of carbapenem-resistant international clone 2 in the Philippines is likely the result of clonal expansion and geographical dissemination, and at least partly explained by inadequate hospital infection control and prevention. DISCUSSION: This is the first extensive genomic study of carbapenem-resistant A. baumannii in the Philippines, and it underscores the importance of hospital infection control and prevention measures to contain high-risk clones.

Genomic surveillance of Pseudomonas aeruginosa in the Philippines, 2013-2014.

Pseudomonas aeruginosa is an opportunistic pathogen that often causes nosocomial infections resistant to treatment. Rates of antimicrobial resistance (AMR) are increasing, as are rates of multidrug-resistant (MDR) and possible extensively drug-resistant (XDR) infections. Our objective was to characterize the molecular epidemiology and AMR mechanisms of this pathogen. We sequenced the whole genome for each of 176 P. aeruginosa isolates collected in the Philippines in 2013-2014; derived the multilocus sequence type (MLST), presence of AMR determinants and relatedness between isolates; and determined concordance between phenotypic and genotypic resistance. Carbapenem resistance was associated with loss of function of the OprD porin and acquisition of the metallo-ß-lactamase (MBL) gene bla VIM. Concordance between phenotypic and genotypic resistance was 93.27% overall for six antibiotics in three classes, but varied among aminoglycosides. The population of P. aeruginosa was diverse, with clonal expansions of XDR genomes belonging to MLSTs ST235, ST244, ST309 and ST773. We found evidence of persistence or reintroduction of the predominant clone ST235 in one hospital, and of transfer between hospitals. Most of the ST235 genomes formed a distinct lineage from global genomes, thus raising the possibility that they may be unique to the Philippines. In addition, long-read sequencing of one representative XDR ST235 isolate identified an integron carrying multiple resistance genes (including bla VIM-2), with differences in gene composition and synteny from the P. aeruginosa class 1 integrons described previously. The survey bridges the gap in genomic data from the Western Pacific Region and will be useful for ongoing surveillance; it also highlights the importance of curtailing the spread of ST235 within the Philippines.