Frank Møller Aarestrup

Identification of acquired antimicrobial resistance genes

Objectives Identification of antimicrobial resistance genes is important for understanding the underlying mechanisms and the epidemiology of antimicrobial resistance. As the costs of whole-genome sequencing (WGS) continue to decline, it becomes increasingly available in routine diagnostic laboratories and is anticipated to substitute traditional methods for resistance gene identification. Thus, the current challenge is to extract the relevant information from the large amount of generated data. Methods We developed a web-based method, ResFinder that uses BLAST for identification of acquired antimicrobial resistance genes in whole-genome data. As input, the method can use both pre-assembled, complete or partial genomes, and short sequence reads from four different sequencing platforms. The method was evaluated on 1862 GenBank files containing 1411 different resistance genes, as well as on 23 de-novo-sequenced isolates. Results When testing the 1862 GenBank files, the method identified the resistance genes with an ID = 100% (100% identity) to the genes in ResFinder. Agreement between in silico predictions and phenotypic testing was found when the method was further tested on 23 isolates of five different bacterial species, with available phenotypes. Furthermore, ResFinder was evaluated on WGS chromosomes and plasmids of 30 isolates. Seven of these isolates were annotated to have antimicrobial resistance, and in all cases, annotations were compatible with the ResFinder results. Conclusions A web server providing a convenient way of identifying acquired antimicrobial resistance genes in completely sequenced isolates was created. ResFinder can be accessed at www.genomicepidemiology.org. ResFinder will continuously be updated as new resistance genes are identified.

In Silico Detection and Typing of Plasmids using PlasmidFinder and Plasmid Multilocus Sequence Typing

ABSTRACT In the work presented here, we designed and developed two easy-to-use Web tools for in silico detection and characterization of whole-genome sequence (WGS) and whole-plasmid sequence data from members of the family Enterobacteriaceae. These tools will facilitate bacterial typing based on draft genomes of multidrug-resistant Enterobacteriaceae species by the rapid detection of known plasmid types. Replicon sequences from 559 fully sequenced plasmids associated with the family Enterobacteriaceae in the NCBI nucleotide database were collected to build a consensus database for integration into a Web tool called PlasmidFinder that can be used for replicon sequence analysis of raw, contig group, or completely assembled and closed plasmid sequencing data. The PlasmidFinder database currently consists of 116 replicon sequences that match with at least at 80% nucleotide identity all replicon sequences identified in the 559 fully sequenced plasmids. For plasmid multilocus sequence typing (pMLST) analysis, a database that is updated weekly was generated from www.pubmlst.org and integrated into a Web tool called pMLST. Both databases were evaluated using draft genomes from a collection of Salmonella enterica serovar Typhimurium isolates. PlasmidFinder identified a total of 103 replicons and between zero and five different plasmid replicons within each of 49 S. Typhimurium draft genomes tested. The pMLST Web tool was able to subtype genomic sequencing data of plasmids, revealing both known plasmid sequence types (STs) and new alleles and ST variants. In conclusion, testing of the two Web tools using both fully assembled plasmid sequences and WGS-generated draft genomes showed them to be able to detect a broad variety of plasmids that are often associated with antimicrobial resistance in clinically relevant bacterial pathogens.

Staphylococcus aureus CC398: Host Adaptation and Emergence of Methicillin Resistance in Livestock

ABSTRACT Since its discovery in the early 2000s, methicillin-resistant Staphylococcus aureus (MRSA) clonal complex 398 (CC398) has become a rapidly emerging cause of human infections, most often associated with livestock exposure. We applied whole-genome sequence typing to characterize a diverse collection of CC398 isolates (n = 89), including MRSA and methicillin-susceptible S. aureus (MSSA) from animals and humans spanning 19 countries and four continents. We identified 4,238 single nucleotide polymorphisms (SNPs) among the 89 core genomes. Minimal homoplasy (consistency index = 0.9591) was detected among parsimony-informative SNPs, allowing for the generation of a highly accurate phylogenetic reconstruction of the CC398 clonal lineage. Phylogenetic analyses revealed that MSSA from humans formed the most ancestral clades. The most derived lineages were composed predominantly of livestock-associated MRSA possessing three different staphylococcal cassette chromosome mec element (SCCmec) types (IV, V, and VII-like) including nine subtypes. The human-associated isolates from the basal clades carried phages encoding human innate immune modulators that were largely missing among the livestock-associated isolates. Our results strongly suggest that livestock-associated MRSA CC398 originated in humans as MSSA. The lineage appears to have undergone a rapid radiation in conjunction with the jump from humans to livestock, where it subsequently acquired tetracycline and methicillin resistance. Further analyses are required to estimate the number of independent genetic events leading to the methicillin-resistant sublineages, but the diversity of SCCmec subtypes is suggestive of strong and diverse antimicrobial selection associated with food animal production. IMPORTANCE Modern food animal production is characterized by densely concentrated animals and routine antibiotic use, which may facilitate the emergence of novel antibiotic-resistant zoonotic pathogens. Our findings strongly support the idea that livestock-associated MRSA CC398 originated as MSSA in humans. The jump of CC398 from humans to livestock was accompanied by the loss of phage-carried human virulence genes, which likely attenuated its zoonotic potential, but it was also accompanied by the acquisition of tetracycline and methicillin resistance. Our findings exemplify a bidirectional zoonotic exchange and underscore the potential public health risks of widespread antibiotic use in food animal production. Modern food animal production is characterized by densely concentrated animals and routine antibiotic use, which may facilitate the emergence of novel antibiotic-resistant zoonotic pathogens. Our findings strongly support the idea that livestock-associated MRSA CC398 originated as MSSA in humans. The jump of CC398 from humans to livestock was accompanied by the loss of phage-carried human virulence genes, which likely attenuated its zoonotic potential, but it was also accompanied by the acquisition of tetracycline and methicillin resistance. Our findings exemplify a bidirectional zoonotic exchange and underscore the potential public health risks of widespread antibiotic use in food animal production.

Multilocus Sequence Typing of Total Genome Sequenced Bacteria

ABSTRACT Accurate strain identification is essential for anyone working with bacteria. For many species, multilocus sequence typing (MLST) is considered the “gold standard” of typing, but it is traditionally performed in an expensive and time-consuming manner. As the costs of whole-genome sequencing (WGS) continue to decline, it becomes increasingly available to scientists and routine diagnostic laboratories. Currently, the cost is below that of traditional MLST. The new challenges will be how to extract the relevant information from the large amount of data so as to allow for comparison over time and between laboratories. Ideally, this information should also allow for comparison to historical data. We developed a Web-based method for MLST of 66 bacterial species based on WGS data. As input, the method uses short sequence reads from four sequencing platforms or preassembled genomes. Updates from the MLST databases are downloaded monthly, and the best-matching MLST alleles of the specified MLST scheme are found using a BLAST-based ranking method. The sequence type is then determined by the combination of alleles identified. The method was tested on preassembled genomes from 336 isolates covering 56 MLST schemes, on short sequence reads from 387 isolates covering 10 schemes, and on a small test set of short sequence reads from 29 isolates for which the sequence type had been determined by traditional methods. The method presented here enables investigators to determine the sequence types of their isolates on the basis of WGS data. This method is publicly available at www.cbs.dtu.dk/services/MLST.

Effect of Abolishment of the Use of Antimicrobial Agents for Growth Promotion on Occurrence of Antimicrobial Resistance in Fecal Enterococci from Food Animals in Denmark

ABSTRACT From 1995 to 2000, a total of 673 Enterococcus faeciumand 1,088 Enterococcus faecalis isolates from pigs together with 856 E. faecium isolates from broilers were isolated and tested for susceptibility to four classes of antimicrobial agents used for growth promotion as part of the Danish program of monitoring for antimicrobial resistance. The four antimicrobials were avilamycin, erythromycin, vancomycin, and virginiamycin. Major changes in the use of antimicrobial agents for growth promotion have occurred during the last 6 years in Denmark. The government banned the use of avoparcin in 1995 and of virginiamycin in 1998. Furthermore, the producers have voluntarily stopped all use beginning in 1999. The avoparcin ban in 1995 was followed by a decrease in the occurrence of glycopeptide-resistant E. faecium (GRE) in broilers, from 72.7% in 1995 to 5.8% in 2000. The occurrence of glycopeptide resistance among isolates from pigs remained constant at around 20% from 1995 to 1997. It was shown that, in GRE from pigs, the genes encoding macrolide and glycopeptide resistance were genetically linked and that, following the decrease in the use of tylosin during 1998 and 1999, the occurrence of GRE in pigs decreased to 6.0% in 2000. From 1995 to 1997 the occurrence of erythromycin resistance among E. faecium and E. faecalis isolates from pigs was almost 90%. Use of tylosin decreased considerably during 1998 and 1999, and this decrease was followed by decreases in the occurrence of resistance to 46.7 and 28.1% among E. faecium and E. faecalis isolates from pigs, respectively. Erythromycin resistance among E. faecium isolates from broilers reached a maximum of 76.3% in 1997 but decreased to 12.7% in 2000 concomitantly with more limited use of virginiamycin. Use of virginiamycin increased from 1995 to 1997 and was followed by an increased occurrence of virginiamycin resistance among E. faecium isolates in broilers, from 27.3% in 1995 to 66.2% in 1997. In January 1998 the use of virginiamycin was banned in Denmark, and the occurrence of virginiamycin resistance decreased to 33.9% in 2000. Use of avilamycin increased from 1995 to 1996 and was followed by an increase in avilamycin resistance among E. faeciumisolates from broilers, from 63.6% in 1995 to 77.4% in 1996. Since 1996 avilamycin usage has decreased, followed by a decrease in resistance to 4.8% in 2000. Our observations show that it is possible to reduce the occurrence of antimicrobial resistance in a national population of food animals when the selective pressure is removed. Cases in which resistance to vancomycin was linked to resistance to erythromycin were exceptions. In such cases resistance did not decrease until the use of both avoparcin and tylosin was limited.

Comparison of antimicrobial resistance phenotypes and resistance genes in Enterococcus faecalis and Enterococcus faecium from humans in the community, broilers, and pigs in Denmark

Enterococcus faecalis and E. faecium isolated from humans in the community (98 and 65 isolates), broilers (126 and 122), and pigs (102 and 88) during 1998 were tested for susceptibility to 12 different antimicrobial agents and for the presence of selected genes encoding resistance using PCR. Furthermore, the presence of vancomycin resistant enterococci was examined in 38 human stool samples using selective enrichment. Widespread resistance to chloramphenicol, macrolides, kanamycin, streptomycin, and tetracycline was found among isolates from all three sources. All E. faecium isolates from humans and pigs were susceptible to avilamycin, whereas 35% of isolates from broilers were resistant. All E. faecium isolates from humans were susceptible to vancomycin, whereas 10% and 17% of isolates from broilers and pigs, respectively, were resistant. A vancomycin resistant E. faecium isolate was found in one of the 38 human fecal samples examined using selective enrichment. All vancomycin resistant isolates contained the vanA gene, all chloramphenicol resistant isolates the cat(pIP501) gene, and all five gentamicin resistant isolates the aac6-aph2 gene. Sixty-one (85%) of 72 erythromycin resistant E. faecalis examined and 57 (90%) of 63 erythromycin resistant E. faecium isolates examined contained ermB. Forty (91%) of the kanamycin resistant E. faecalis and 18 (72%) of the kanamycin resistant E. faecium isolates contained aphA3. The tet(M) gene was found in 95% of the tetracycline resistant E. faecalis and E. faecium isolates of human and animal origin, examined. tet(K) was not observed, whereas tet(L) was detected in 17% of tetracycline resistant E. faecalis isolates and in 16% of the E. faecium isolates. tet(O) was not detected in any of the isolates from pigs, but was observed in 38% of E. faecalis isolates from broilers, in two E. faecalis isolates from humans and in three E. faecium isolates from broilers. tet(S) was not detected among isolates from animals, but was observed in 31% of E. faecalis and one E. faecium isolate from humans. This study showed a frequent occurrence of antimicrobial resistance and the presence of selected resistance genes in E. faecalis and E. faecium isolated from humans, broilers and pigs. Differences in the occurrence of resistance and tetracycline resistance genes were observed among isolates from the different sources. However, similar resistance patterns and resistance genes were detected frequently indicating that transmission of resistant enterococci or resistance genes takes place between humans, broilers, and pigs.

qnrD, a Novel Gene Conferring Transferable Quinolone Resistance in Salmonella enterica Serovar Kentucky and Bovismorbificans Strains of Human Origin

ABSTRACT In a previous study, four Salmonella isolates from humans in the Henan province of China showed reduced susceptibility to ciprofloxacin (MIC, 0.125 to 0.25 μg/ml) but were susceptible to nalidixic acid (MIC, 4 to 8 μg/ml). All isolates were negative for known qnr genes (A, B, and S), aac(6′)Ib-cr, and mutations in gyrA and parC. Plasmid DNA was extracted from all four isolates and transformed into Escherichia coli TG1 and DH10B cells by electroporation, and transformants were selected on 0.06 μg/ml ciprofloxacin containing brain heart infusion agar plates. Resistance to ciprofloxacin could be transferred by electroporation, and a similar 4,270-bp plasmid was found in all transformants. By sequence analysis, the plasmid was found to carry an open reading frame that had similarities to other qnr genes and that encoded a 214-amino-acid pentapeptide repeat protein. This gene, designated qnrD, showed 48% similarity to qnrA1, 61% similarity to qnrB1, and 41% similarity to qnrS1. Further subcloning of the qnrD coding region into the constitutively expressed tetA gene of vector pBR322 showed that the gene conferred an increase in the MIC of ciprofloxacin by a factor of 32 (from an MIC of 0.002 to an MIC of 0.06 μg/ml). For comparison, qnrA1 and qnrS1 were also subcloned into pBR322 and transformed into DH10B cells, conferring MICs of 0.125 and 0.5 μg/ml, respectively. A phylogenetic analysis of all known qnr sequences was performed and showed that qnrD was more closely related to the qnrB variants but formed an independent cluster. To our knowledge, this is the first description of this qnrD gene.

Avoparcin used as a growth promoter is associated with the occurrence of vancomycin-resistant Enterococcus faecium on Danish poultry and pig farms

We determined the association between the use of the glycopeptide antibiotic avoparcin as a growth promoter and the occurrence of Enterococcus faecium (VREF) with high-level resistance to vancomycin (MIC > or = 64 micrograms ml-1) on poultry and pig farms. The investigations were conducted as retrospective cohort studies, where groups of farms exposed or not exposed to avoparcin between September 1994 and April 1995 were compared. In poultry, the association between the use of avoparcin and the occurrence of VREF was confounded by the use of broad-spectrum antibiotics, and the adjusted relative risk was 2.9 (1.4-5.9). In pigs, the association had a similar magnitude with a non-adjusted relative risk of 3.3 (0.9-12.3). The similar findings in the two studies provide evidence in favour of a causal association between the use of avoparcin and the occurrence of VREF on farms, and suggest that food animals constitute a potential reservoir of infection for VREF in humans.

Global Monitoring of Salmonella Serovar Distribution from the World Health Organization Global Foodborne Infections Network Country Data Bank: Results of Quality Assured Laboratories from 2001 to 2007

Salmonella enterica is commonly acquired from contaminated food and is an important cause of illness worldwide. Interventions are needed to control Salmonella; subtyping Salmonella by serotyping is useful for targeting such interventions. We, therefore, analyzed the global distribution of the 15 most frequently identified serovars of Salmonella isolated from humans from 2001 to 2007 in laboratories from 37 countries that participated in World Health Organization Global Foodborne Infections Network and demonstrated serotyping proficiency in the Global Foodborne Infections Network External Quality Assurance System. In all regions throughout the study period, with the exception of the Oceania and North American regions, Salmonella serovars Enteritidis and Typhimurium ranked as the most common and second most common serovar, respectively. In the North American and Oceania (Australia and New Zealand) regions, Salmonella serovar Typhimurium was the most common serovar reported, and Salmonella serovar Enteritidis was the second most common serovar. During the study period, the proportion of Salmonella isolates reported from humans that were Salmonella serovar Enteritidis was 43.5% (range: 40.6% [2007] to 44.9% [2003]), and Salmonella serovar Typhimurium was 17.1% (range: 15% [2007] to 18.9% [2001]). Salmonella serovars Newport (mainly observed in Latin and North American and European countries), Infantis (dominating in all regions), Virchow (mainly observed in Asian, European, and Oceanic countries), Hadar (profound in European countries), and Agona (intense in Latin and North American and European countries) were also frequently isolated with an overall proportion of 3.5%, 1.8%, 1.5%, 1.5%, and 0.8%, respectively. There were large differences in the most commonly isolated serovars between regions, but lesser differences between countries within the same region. The results also highlight the complexity of the global epidemiology of Salmonella and the need and importance for improving monitoring data of those serovars of highest epidemiologic importance.

Real-Time Whole-Genome Sequencing for Routine Typing, Surveillance, and Outbreak Detection of Verotoxigenic Escherichia coli

ABSTRACT Fast and accurate identification and typing of pathogens are essential for effective surveillance and outbreak detection. The current routine procedure is based on a variety of techniques, making the procedure laborious, time-consuming, and expensive. With whole-genome sequencing (WGS) becoming cheaper, it has huge potential in both diagnostics and routine surveillance. The aim of this study was to perform a real-time evaluation of WGS for routine typing and surveillance of verocytotoxin-producing Escherichia coli (VTEC). In Denmark, the Statens Serum Institut (SSI) routinely receives all suspected VTEC isolates. During a 7-week period in the fall of 2012, all incoming isolates were concurrently subjected to WGS using IonTorrent PGM. Real-time bioinformatics analysis was performed using web-tools (www.genomicepidemiology.org) for species determination, multilocus sequence type (MLST) typing, and determination of phylogenetic relationship, and a specific VirulenceFinder for detection of E. coli virulence genes was developed as part of this study. In total, 46 suspected VTEC isolates were characterized in parallel during the study. VirulenceFinder proved successful in detecting virulence genes included in routine typing, explicitly verocytotoxin 1 (vtx1), verocytotoxin 2 (vtx2), and intimin (eae), and also detected additional virulence genes. VirulenceFinder is also a robust method for assigning verocytotoxin (vtx) subtypes. A real-time clustering of isolates in agreement with the epidemiology was established from WGS, enabling discrimination between sporadic and outbreak isolates. Overall, WGS typing produced results faster and at a lower cost than the current routine. Therefore, WGS typing is a superior alternative to conventional typing strategies. This approach may also be applied to typing and surveillance of other pathogens.