Emergence of methicillin resistance predates the clinical use of antibiotics.

The discovery of antibiotics more than 80 years ago has led to considerable improvements in human and animal health. Although antibiotic resistance in environmental bacteria is ancient, resistance in human pathogens is thought to be a modern phenomenon that is driven by the clinical use of antibiotics1. Here we show that particular lineages of methicillin-resistant Staphylococcus aureus-a notorious human pathogen-appeared in European hedgehogs in the pre-antibiotic era. Subsequently, these lineages spread within the local hedgehog populations and between hedgehogs and secondary hosts, including livestock and humans. We also demonstrate that the hedgehog dermatophyte Trichophyton erinacei produces two ß-lactam antibiotics that provide a natural selective environment in which methicillin-resistant S. aureus isolates have an advantage over susceptible isolates. Together, these results suggest that methicillin resistance emerged in the pre-antibiotic era as a co-evolutionary adaptation of S. aureus to the colonization of dermatophyte-infected hedgehogs. The evolution of clinically relevant antibiotic-resistance genes in wild animals and the connectivity of natural, agricultural and human ecosystems demonstrate that the use of a One Health approach is critical for our understanding and management of antibiotic resistance, which is one of the biggest threats to global health, food security and development.

West Nile virus spread in Europe: Phylogeographic pattern analysis and key drivers.

BACKGROUND: West Nile virus (WNV) outbreaks in birds, humans, and livestock have occurred in multiple areas in Europe and have had a significant impact on animal and human health. The patterns of emergence and spread of WNV in Europe are very different from those in the US and understanding these are important for guiding preparedness activities. METHODS: We mapped the evolution and spread history of WNV in Europe by incorporating viral genome sequences and epidemiological data into phylodynamic models. Spatially explicit phylogeographic models were developed to explore the possible contribution of different drivers to viral dispersal direction and velocity. A "skygrid-GLM" approach was used to identify how changes in environments would predict viral genetic diversity variations over time. FINDINGS: Among the six lineages found in Europe, WNV-2a (a sub-lineage of WNV-2) has been predominant (accounting for 73% of all sequences obtained in Europe that have been shared in the public domain) and has spread to at least 14 countries. In the past two decades, WNV-2a has evolved into two major co-circulating clusters, both originating from Central Europe, but with distinct dynamic history and transmission patterns. WNV-2a spreads at a high dispersal velocity (88km/yr-215 km/yr) which is correlated to bird movements. Notably, amongst multiple drivers that could affect the spread of WNV, factors related to land use were found to strongly influence the spread of WNV. Specifically, the intensity of agricultural activities (defined by factors related to crops and livestock production, such as coverage of cropland, pasture, cultivated and managed vegetation, livestock density) were positively associated with both spread direction and velocity. In addition, WNV spread direction was associated with high coverage of wetlands and migratory bird flyways. CONCLUSION: Our results suggest that-in addition to ecological conditions favouring bird- and mosquito- presence-agricultural land use may be a significant driver of WNV emergence and spread. Our study also identified significant gaps in data and the need to strengthen virological surveillance in countries of Central Europe from where WNV outbreaks are likely seeded. Enhanced monitoring for early detection of further dispersal could be targeted to areas with high agricultural activities and habitats of migratory birds.

ARGprofiler-a pipeline for large-scale analysis of antimicrobial resistance genes and their flanking regions in metagenomic datasets.

MOTIVATION: Analyzing metagenomic data can be highly valuable for understanding the function and distribution of antimicrobial resistance genes (ARGs). However, there is a need for standardized and reproducible workflows to ensure the comparability of studies, as the current options involve various tools and reference databases, each designed with a specific purpose in mind. RESULTS: In this work, we have created the workflow ARGprofiler to process large amounts of raw sequencing reads for studying the composition, distribution, and function of ARGs. ARGprofiler tackles the challenge of deciding which reference database to use by providing the PanRes database of 14 078 unique ARGs that combines several existing collections into one. Our pipeline is designed to not only produce abundance tables of genes and microbes but also to reconstruct the flanking regions of ARGs with ARGextender. ARGextender is a bioinformatic approach combining KMA and SPAdes to recruit reads for a targeted de novo assembly. While our aim is on ARGs, the pipeline also creates Mash sketches for fast searching and comparisons of sequencing runs. AVAILABILITY AND IMPLEMENTATION: The ARGprofiler pipeline is a Snakemake workflow that supports the reuse of metagenomic sequencing data and is easily installable and maintained at https://github.com/genomicepidemiology/ARGprofiler.

A curated data resource of 214K metagenomes for characterization of the global antimicrobial resistome.

The growing threat of antimicrobial resistance (AMR) calls for new epidemiological surveillance methods, as well as a deeper understanding of how antimicrobial resistance genes (ARGs) have been transmitted around the world. The large pool of sequencing data available in public repositories provides an excellent resource for monitoring the temporal and spatial dissemination of AMR in different ecological settings. However, only a limited number of research groups globally have the computational resources to analyze such data. We retrieved 442 Tbp of sequencing reads from 214,095 metagenomic samples from the European Nucleotide Archive (ENA) and aligned them using a uniform approach against ARGs and 16S/18S rRNA genes. Here, we present the results of this extensive computational analysis and share the counts of reads aligned. Over 6.76∙108 read fragments were assigned to ARGs and 3.21∙109 to rRNA genes, where we observed distinct differences in both the abundance of ARGs and the link between microbiome and resistome compositions across various sampling types. This collection is another step towards establishing global surveillance of AMR and can serve as a resource for further research into the environmental spread and dynamic changes of ARGs.

Multiclonal human origin and global expansion of an endemic bacterial pathogen of livestock.

Most new pathogens of humans and animals arise via switching events from distinct host species. However, our understanding of the evolutionary and ecological drivers of successful host adaptation, expansion, and dissemination are limited. Staphylococcus aureus is a major bacterial pathogen of humans and a leading cause of mastitis in dairy cows worldwide. Here we trace the evolutionary history of bovine S. aureus using a global dataset of 10,254 S. aureus genomes including 1,896 bovine isolates from 32 countries in 6 continents. We identified 7 major contemporary endemic clones of S. aureus causing bovine mastitis around the world and traced them back to 4 independent host-jump events from humans that occurred up to 2,500 y ago. Individual clones emerged and underwent clonal expansion from the mid-19th to late 20th century coinciding with the commercialization and industrialization of dairy farming, and older lineages have become globally distributed via established cattle trade links. Importantly, we identified lineage-dependent differences in the frequency of host transmission events between humans and cows in both directions revealing high risk clones threatening veterinary and human health. Finally, pangenome network analysis revealed that some bovine S. aureus lineages contained distinct sets of bovine-associated genes, consistent with multiple trajectories to host adaptation via gene acquisition. Taken together, we have dissected the evolutionary history of a major endemic pathogen of livestock providing a comprehensive temporal, geographic, and gene-level perspective of its remarkable success.

Ongoing mpox outbreak in Kamituga, South Kivu province, associated with monkeypox virus of a novel Clade I sub-lineage, Democratic Republic of the Congo, 2024.

Since the beginning of 2023, the number of people with suspected monkeypox virus (MPXV) infection have sharply increased in the Democratic Republic of the Congo (DRC). We report near-to-complete MPXV genome sequences derived from six cases from the South Kivu province. Phylogenetic analyses reveal that the MPXV affecting the cases belongs to a novel Clade I sub-lineage. The outbreak strain genome lacks the target sequence of the probe and primers of a commonly used Clade I-specific real-time PCR.

Various arrangements of mobile genetic elements among CC147 subpopulations of Klebsiella pneumoniae harboring blaNDM-1: a comparative genomic analysis of carbapenem resistant strains.

BACKGROUND: Certain clonal complexes (CCs) of Klebsiella pneumoniae such as CC147 (ST147 and ST392) are major drivers of blaNDM dissemination across the world. ST147 has repeatedly reported from our geographical region, but its population dynamics and evolutionary trajectories need to be further studied. METHODS: Comparative genomic analysis of 51 carbapenem-nonsusceptible strains as well as three hypervirulent K. pneumoniae (hvKp) recovered during 16-months of surveillance was performed using various bioinformatics tools. We investigated the genetic proximity of our ST147 strains with publicly available corresponding genomes deposited globally and from neighbor countries in our geographic region. RESULTS: While IncL/M plasmid harboring blaOXA-48 was distributed among divergent clones, blaNDM-1 was circulated by twenty of the 25 CC147 dominant clone and were mostly recovered from the ICU. The NDM-1 core structure was bracketed by a single isoform of mobile genetic elements (MGEs) [ΔISKpn26-NDM-TnAs3-ΔIS3000-Tn5403] and was located on Col440I plasmid in 68.7% of ST392. However, various arrangements of MGEs including MITESen1/MITESen1 composite transposon or combination of MITESen1/ISSen4/IS903B/IS5/ISEhe3 on IncFIb (pB171) were identified in ST147. It seems that ST392 circulated blaNDM-1 in 2018 before being gradually replaced by ST147 from the middle to the end of sample collection in 2019. ST147 strains possessed the highest number of resistance markers and showed high genetic similarity with four public genomes that harbored blaNDM-1 on the same replicon type. Mainly, there was a convergence between clusters and isolated neighboring countries in the minimum-spanning tree. A conserved arrangement of resistance markers/MGEs was linked to methyltransferase armA which was embedded in class 1 integron in 8 isolates of ST147/ST48 high-risk clones. CONCLUSION: Our findings highlight the dynamic nature of blaNDM-1 transmission among K. pneumoniae in Iran that occurs both clonally and horizontally via various combinations of MGEs. This is the first analysis of Iranian ST147/NDM + clone in the global context.

Antimicrobial resistance monitoring in the Danish swine production by phenotypic methods and metagenomics from 1999 to 2018.

BackgroundIn Denmark, antimicrobial resistance (AMR) in pigs has been monitored since 1995 by phenotypic approaches using the same indicator bacteria. Emerging methodologies, such as metagenomics, may allow novel surveillance ways.AimThis study aimed to assess the relevance of indicator bacteria (Escherichia coli and Enterococcus faecalis) for AMR surveillance in pigs, and the utility of metagenomics.MethodsWe collated existing data on AMR and antimicrobial use (AMU) from the Danish surveillance programme and performed metagenomics sequencing on caecal samples that had been collected/stored through the programme during 1999-2004 and 2015-2018. We compared phenotypic and metagenomics results regarding AMR, and the correlation of both with AMU.ResultsVia the relative abundance of AMR genes, metagenomics allowed to rank these genes as well as the AMRs they contributed to, by their level of occurrence. Across the two study periods, resistance to aminoglycosides, macrolides, tetracycline, and beta-lactams appeared prominent, while resistance to fosfomycin and quinolones appeared low. In 2015-2018 sulfonamide resistance shifted from a low occurrence category to an intermediate one. Resistance to glycopeptides consistently decreased during the entire study period. Outcomes of both phenotypic and metagenomics approaches appeared to positively correlate with AMU. Metagenomics further allowed to identify multiple time-lagged correlations between AMU and AMR, the most evident being that increased macrolide use in sow/piglets or fatteners led to increased macrolide resistance with a lag of 3-6 months.ConclusionWe validated the long-term usefulness of indicator bacteria and showed that metagenomics is a promising approach for AMR surveillance.

An Overview of Antimicrobial Resistance Profiles of Publicly Available Salmonella Genomes with Sufficient Quality and Metadata.

Salmonella enterica (S. enterica) is a commensal organism or pathogen causing diseases in animals and humans, as well as widespread in the environment. Antimicrobial resistance (AMR) has increasingly affected both animal and human health and continues to raise public health concerns. A decade ago, it was estimated that the increased use of whole genome sequencing (WGS) combined with sharing of public data would drastically change and improve the surveillance and understanding of Salmonella epidemiology and AMR. This study aimed to evaluate the current usefulness of public WGS data for Salmonella surveillance and to investigate the associations between serovars, antibiotic resistance genes (ARGs), and metadata. Out of 191,306 Salmonella genomes deposited in European Nucleotide Archive and NCBI databases, 47,452 WGS with sufficient minimum metadata (country, year, and source) of S. enterica were retrieved from 116 countries and isolated between 1905 and 2020. For in silico analysis of the WGS data, KmerFinder, SISTR, and ResFinder were used for species, serovars, and AMR identification, respectively. The results showed that the five common isolation sources of S. enterica are human (29.10%), avian (22.50%), environment (11.89%), water (9.33%), and swine (6.62%). The most common ARG profiles for each class of antimicrobials are ß-lactam (blaTEM-1B; 6.78%), fluoroquinolone [(parC[T57S], qnrB19); 0.87%], folate pathway antagonist (sul2; 8.35%), macrolide [mph(A); 0.39%], phenicol (floR; 5.94%), polymyxin B (mcr-1.1; 0.09%), and tetracycline [tet(A); 12.95%]. Our study reports the first overview of ARG profiles in publicly available Salmonella genomes from online databases. All data sets from this study can be searched at Microreact.

Rice protein concentrate as a fish meal substitute in Oreochromis niloticus: Effects on immune response, intestinal cytokines, Aeromonas veronii resistance, and gut microbiota composition.

The potential of rice protein concentrate (RPC) to substitute fishmeal (FM) protein in the diet of Oreochromis niloticus was assessed in a five-month-long feeding trial. Fishmeal protein was replaced by RPC at rates of 0% (control), 25%, 50%, and 75% (RPC0, RPC25, RPC50, and RPC75, respectively). RPC25 had no significant effect on antioxidant capacity (total antioxidant capacity; superoxide dismutase, catalase, and glutathione peroxidase activities) and immune indices (lysozyme, nitric oxide, antiprotease, and bactericidal activities) after one, two, and five months of feeding, while the values for these parameters were significantly lower in the RPC75 group compared to those in the RPC0 group. The RPC25 group showed higher mRNA levels of the intestinal cytokines IL-1ß, IL-10ß, TGF-ß, and TNF-α than the control group. In fish affected by Aeromonas veronii, the highest significant cumulative mortality was recorded in the RPC75 group, followed by the RPC50, RPC25, and control groups. Gut microbiome analyses showed a reduction in microbial diversity in response to the addition of RPC, regardless of the RPC content, and the composition of the community of the RPC samples differed from that of the control. RPC-enriched diets resulted in higher relative abundances of Bacteroidetes and Fusobacteria in the gut compared to that in the gut of the control fish. In summary, RPC can be used to replace up to 25% of the FM protein in the diet of O. niloticus, while improving the antioxidant capacity, immunocompetence, and disease resistance of the fish.

Monitoring SARS-CoV-2 Circulation and Diversity through Community Wastewater Sequencing, the Netherlands and Belgium.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has rapidly become a major global health problem, and public health surveillance is crucial to monitor and prevent virus spread. Wastewater-based epidemiology has been proposed as an addition to disease-based surveillance because virus is shed in the feces of ≈40% of infected persons. We used next-generation sequencing of sewage samples to evaluate the diversity of SARS-CoV-2 at the community level in the Netherlands and Belgium. Phylogenetic analysis revealed the presence of the most prevalent clades (19A, 20A, and 20B) and clustering of sewage samples with clinical samples from the same region. We distinguished multiple clades within a single sewage sample by using low-frequency variant analysis. In addition, several novel mutations in the SARS-CoV-2 genome were detected. Our results illustrate how wastewater can be used to investigate the diversity of SARS-CoV-2 viruses circulating in a community and identify new outbreaks.

Detection of mobile genetic elements associated with antibiotic resistance in Salmonella enterica using a newly developed web tool: MobileElementFinder.

OBJECTIVES: Antimicrobial resistance (AMR) in clinically relevant bacteria is a growing threat to public health globally. In these bacteria, antimicrobial resistance genes are often associated with mobile genetic elements (MGEs), which promote their mobility, enabling them to rapidly spread throughout a bacterial community. METHODS: The tool MobileElementFinder was developed to enable rapid detection of MGEs and their genetic context in assembled sequence data. MGEs are detected based on sequence similarity to a database of 4452 known elements augmented with annotation of resistance genes, virulence factors and detection of plasmids. RESULTS: MobileElementFinder was applied to analyse the mobilome of 1725 sequenced Salmonella enterica isolates of animal origin from Denmark, Germany and the USA. We found that the MGEs were seemingly conserved according to multilocus ST and not restricted to either the host or the country of origin. Moreover, we identified putative translocatable units for specific aminoglycoside, sulphonamide and tetracycline genes. Several putative composite transposons were predicted that could mobilize, among others, AMR, metal resistance and phosphodiesterase genes associated with macrophage survivability. This is, to our knowledge, the first time the phosphodiesterase-like pdeL has been found to be potentially mobilized into S. enterica. CONCLUSIONS: MobileElementFinder is a powerful tool to study the epidemiology of MGEs in a large number of genome sequences and to determine the potential for genomic plasticity of bacteria. This web service provides a convenient method of detecting MGEs in assembled sequence data. MobileElementFinder can be accessed at https://cge.cbs.dtu.dk/services/MobileElementFinder/.

Salmonella enterica serovar Typhi H58 clone has been endemic in Zimbabwe from 2012 to 2019.

BACKGROUND: Typhoid fever, caused by S. enterica ser. Typhi, continues to be a substantial health burden in developing countries. Little is known of the genotypic diversity of S. enterica ser. Typhi in Zimbabwe, but this is key for understanding the emergence and spread of this pathogen and devising interventions for its control. OBJECTIVES: To report the molecular epidemiology of S. enterica ser. Typhi outbreak strains circulating from 2012 to 2019 in Zimbabwe, using comparative genomics. METHODS: A review of typhoid cases records from 2012 to 2019 in Zimbabwe was performed. The phylogenetic relationship of outbreak isolates from 2012 to 2019 and emergence of antibiotic resistance was investigated by whole-genome sequence analysis. RESULTS: A total 22 479 suspected typhoid cases, 760 confirmed cases were reported from 2012 to 2019 and 29 isolates were sequenced. The majority of the sequenced isolates were predicted to confer resistance to aminoglycosides, ß-lactams, phenicols, sulphonamides, tetracycline and fluoroquinolones (including qnrS detection). The qnrS1 gene was associated with an IncN (subtype PST3) plasmid in 79% of the isolates. Whole-genome SNP analysis, SNP-based haplotyping and resistance determinant analysis showed that 93% of the isolates belonged to a single clade represented by multidrug-resistant H58 lineage I (4.3.1.1), with a maximum pair-wise distance of 22 SNPs. CONCLUSIONS: This study has provided detailed genotypic characterization of the outbreak strain, identified as S. Typhi 4.3.1.1 (H58). The strain has reduced susceptibility to ciprofloxacin due to qnrS carried by an IncN (subtype PST3) plasmid resulting from ongoing evolution to full resistance.

Economic evaluation of whole genome sequencing for pathogen identification and surveillance - results of case studies in Europe and the Americas 2016 to 2019.

BackgroundWhole genome sequencing (WGS) is increasingly used for pathogen identification and surveillance.AimWe evaluated costs and benefits of routine WGS through case studies at eight reference laboratories in Europe and the Americas which conduct pathogen surveillance for avian influenza (two laboratories), human influenza (one laboratory) and food-borne pathogens (five laboratories).MethodsThe evaluation focused on the institutional perspective, i.e. the 'investment case' for implementing WGS compared with conventional methods, based on costs and benefits during a defined reference period, mostly covering at least part of 2017. A break-even analysis estimated the number of cases of illness (for the example of Salmonella surveillance) that would need to be avoided through WGS in order to 'break even' on costs.ResultsOn a per-sample basis, WGS was between 1.2 and 4.3 times more expensive than routine conventional methods. However, WGS brought major benefits for pathogen identification and surveillance, substantially changing laboratory workflows, analytical processes and outbreaks detection and control. Between 0.2% and 1.1% (on average 0.7%) of reported salmonellosis cases would need to be prevented to break even with respect to the additional costs of WGS.ConclusionsEven at cost levels documented here, WGS provides a level of additional information that more than balances the additional costs if used effectively. The substantial cost differences for WGS between reference laboratories were due to economies of scale, degree of automation, sequencing technology used and institutional discounts for equipment and consumables, as well as the extent to which sequencers are used at full capacity.

ResFinder 4.0 for predictions of phenotypes from genotypes.

OBJECTIVES: WGS-based antimicrobial susceptibility testing (AST) is as reliable as phenotypic AST for several antimicrobial/bacterial species combinations. However, routine use of WGS-based AST is hindered by the need for bioinformatics skills and knowledge of antimicrobial resistance (AMR) determinants to operate the vast majority of tools developed to date. By leveraging on ResFinder and PointFinder, two freely accessible tools that can also assist users without bioinformatics skills, we aimed at increasing their speed and providing an easily interpretable antibiogram as output. METHODS: The ResFinder code was re-written to process raw reads and use Kmer-based alignment. The existing ResFinder and PointFinder databases were revised and expanded. Additional databases were developed including a genotype-to-phenotype key associating each AMR determinant with a phenotype at the antimicrobial compound level, and species-specific panels for in silico antibiograms. ResFinder 4.0 was validated using Escherichia coli (n = 584), Salmonella spp. (n = 1081), Campylobacter jejuni (n = 239), Enterococcus faecium (n = 106), Enterococcus faecalis (n = 50) and Staphylococcus aureus (n = 163) exhibiting different AST profiles, and from different human and animal sources and geographical origins. RESULTS: Genotype-phenotype concordance was ≥95% for 46/51 and 25/32 of the antimicrobial/species combinations evaluated for Gram-negative and Gram-positive bacteria, respectively. When genotype-phenotype concordance was <95%, discrepancies were mainly linked to criteria for interpretation of phenotypic tests and suboptimal sequence quality, and not to ResFinder 4.0 performance. CONCLUSIONS: WGS-based AST using ResFinder 4.0 provides in silico antibiograms as reliable as those obtained by phenotypic AST at least for the bacterial species/antimicrobial agents of major public health relevance considered.

Proficiency Testing of Virus Diagnostics Based on Bioinformatics Analysis of Simulated In Silico High-Throughput Sequencing Data Sets.

Quality management and independent assessment of high-throughput sequencing-based virus diagnostics have not yet been established as a mandatory approach for ensuring comparable results. The sensitivity and specificity of viral high-throughput sequence data analysis are highly affected by bioinformatics processing using publicly available and custom tools and databases and thus differ widely between individuals and institutions. Here we present the results of the COMPARE [Collaborative Management Platform for Detection and Analyses of (Re-)emerging and Foodborne Outbreaks in Europe] in silico virus proficiency test. An artificial, simulated in silico data set of Illumina HiSeq sequences was provided to 13 different European institutes for bioinformatics analysis to identify viral pathogens in high-throughput sequence data. Comparison of the participants' analyses shows that the use of different tools, programs, and databases for bioinformatics analyses can impact the correct identification of viral sequences from a simple data set. The identification of slightly mutated and highly divergent virus genomes has been shown to be most challenging. Furthermore, the interpretation of the results, together with a fictitious case report, by the participants showed that in addition to the bioinformatics analysis, the virological evaluation of the results can be important in clinical settings. External quality assessment and proficiency testing should become an important part of validating high-throughput sequencing-based virus diagnostics and could improve the harmonization, comparability, and reproducibility of results. There is a need for the establishment of international proficiency testing, like that established for conventional laboratory tests such as PCR, for bioinformatics pipelines and the interpretation of such results.

Using WGS to identify antibiotic resistance genes and predict antimicrobial resistance phenotypes in MDR Acinetobacter baumannii in Tanzania.

BACKGROUND: Reliable phenotypic antimicrobial susceptibility testing can be a challenge in clinical settings in low- and middle-income countries. WGS is a promising approach to enhance current capabilities. AIM: To study diversity and resistance determinants and to predict and compare resistance patterns from WGS data of Acinetobacter baumannii with phenotypic results from classical microbiological testing at a tertiary care hospital in Tanzania. METHODS AND RESULTS: MLST using Pasteur/Oxford schemes yielded eight different STs from each scheme. Of the eight, two STs were identified to be global clones 1 (n = 4) and 2 (n = 1) as per the Pasteur scheme. Resistance testing using classical microbiology determined between 50% and 92.9% resistance across all drugs. Percentage agreement between phenotypic and genotypic prediction of resistance ranged between 57.1% and 100%, with coefficient of agreement (κ) between 0.05 and 1. Seven isolates harboured mutations at significant loci (S81L in gyrA and S84L in parC). A number of novel plasmids were detected, including pKCRI-309C-1 (219000 bp) carrying 10 resistance genes, pKCRI-43-1 (34935 bp) carrying two resistance genes and pKCRI-49-1 (11681 bp) and pKCRI-28-1 (29606 bp), each carrying three resistance genes. New ampC alleles detected included ampC-69, ampC-70 and ampC-71. Global clone 1 and 2 isolates were found to harbour ISAba1 directly upstream of the ampC gene. Finally, SNP-based phylogenetic analysis of the A. baumannii isolates revealed closely related isolates in three clusters. CONCLUSIONS: The validity of the use of WGS in the prediction of phenotypic resistance can be appreciated, but at this stage is not sufficient for it to replace conventional antimicrobial susceptibility testing in our setting.

The antimicrobial resistome in relation to antimicrobial use and biosecurity in pig farming, a metagenome-wide association study in nine European countries.

OBJECTIVES: Previous studies in food-producing animals have shown associations between antimicrobial use (AMU) and resistance (AMR) in specifically isolated bacterial species. Multi-country data are scarce and only describe between-country differences. Here we investigate associations between the pig faecal mobile resistome and characteristics at the farm-level across Europe. METHODS: A cross-sectional study was conducted among 176 conventional pig farms from nine European countries. Twenty-five faecal samples from fattening pigs were pooled per farm and acquired resistomes were determined using shotgun metagenomics and the Resfinder reference database, i.e. the full collection of horizontally acquired AMR genes (ARGs). Normalized fragments resistance genes per kilobase reference per million bacterial fragments (FPKM) were calculated. Specific farm-level data (AMU, biosecurity) were collected. Random-effects meta-analyses were performed by country, relating farm-level data to relative ARG abundances (FPKM). RESULTS: Total AMU during fattening was positively associated with total ARG (total FPKM). Positive associations were particularly observed between widely used macrolides and tetracyclines, and ARGs corresponding to the respective antimicrobial classes. Significant AMU-ARG associations were not found for ß-lactams and only few colistin ARGs were found, despite high use of these antimicrobial classes in younger pigs. Increased internal biosecurity was directly related to higher abundances of ARGs mainly encoding macrolide resistance. These effects of biosecurity were independent of AMU in mutually adjusted models. CONCLUSIONS: Using resistome data in association studies is unprecedented and adds accuracy and new insights to previously observed AMU-AMR associations. Major components of the pig resistome are positively and independently associated with on-farm AMU and biosecurity conditions.

Associations between antimicrobial use and the faecal resistome on broiler farms from nine European countries.

OBJECTIVES: To determine associations between farm- and flock-level antimicrobial usage (AMU), farm biosecurity status and the abundance of faecal antimicrobial resistance genes (ARGs) on broiler farms. METHODS: In the cross-sectional pan-European EFFORT study, conventional broiler farms were visited and faeces, AMU information and biosecurity records were collected. The resistomes of pooled faecal samples were determined by metagenomic analysis for 176 farms. A meta-analysis approach was used to relate total and class-specific ARGs (expressed as fragments per kb reference per million bacterial fragments, FPKM) to AMU (treatment incidence per DDD, TIDDDvet) per country and subsequently across all countries. In a similar way, the association between biosecurity status (Biocheck.UGent) and the resistome was explored. RESULTS: Sixty-six (38%) flocks did not report group treatments but showed a similar resistome composition and roughly similar ARG levels to antimicrobial-treated flocks. Nevertheless, we found significant positive associations between ß-lactam, tetracycline, macrolide and lincosamide, trimethoprim and aminoglycoside antimicrobial flock treatments and ARG clusters conferring resistance to the same class. Similar associations were found with purchased products. In gene-level analysis for ß-lactams and macrolides, lincosamides and streptogramins, a significant positive association was found with the most abundant gene clusters blaTEM and erm(B). Little evidence was found for associations with biosecurity. CONCLUSIONS: The faecal microbiome in European broilers contains a high diversity of ARGs, even in the absence of current antimicrobial selection pressure. Despite this, the relative abundance of genes and the composition of the resistome is positively related to AMU in European broiler farms for several antimicrobial classes.

Proof of concept: used malaria rapid diagnostic tests applied for parallel sequencing for surveillance of molecular markers of anti-malarial resistance in Bissau, Guinea-Bissau during 2014-2017.

BACKGROUND: Large-scale surveillance of molecular markers of anti-malarial drug resistance is an attractive method of resistance monitoring, to complement therapeutic efficacy studies in settings where the latter are logistically challenging. METHODS: Between 2014 and 2017, this study sampled malaria rapid diagnostic tests (RDTs), used in routine clinical care, from two health centres in Bissau, Guinea-Bissau. In order to obtain epidemiological insights, RDTs were collected together with patient data on age and sex. A subset of positive RDTs from one of the two sites (n = 2184) were tested for Plasmodium DNA content. Those testing positive for Plasmodium DNA by PCR (n = 1390) were used for library preparation, custom designed dual indexing and next generation Miseq targeted sequencing of Plasmodium falciparum genes pfcrt, pfmdr1, pfdhfr, pfdhps and pfk13. RESULTS: The study found a high frequency of the pfmdr1 codon 86N at 88-97%, a significant decrease of the pfcrt wildtype CVMNK haplotype and elevated levels of the pfdhfr/pfdhps quadruple mutant ranging from 33 to 51% between 2014 and 2017. No polymorphisms indicating artemisinin tolerance were discovered. The demographic data indicate a large proportion of young adults (66%, interquartile range 11-28 years) presenting with P. falciparum infections. While a total of 5532 gene fragments were successfully analysed on a single Illumina Miseq flow cell, PCR-positivity from the library preparation varied considerably from 13 to 87% for different amplicons. Furthermore, pre-screening of samples for Plasmodium DNA content proved necessary prior to library preparation. CONCLUSIONS: This study serves as a proof of concept for using leftover clinical material (used RDTs) for large-scale molecular surveillance, encompassing the inherent complications regarding to methodology and analysis when doing so. Factors such as RDT storage prior to DNA extraction and parasitaemia of the infection are likely to have an effect on whether or not parasite DNA can be successfully analysed, and are considered part of the reason the data yield is suboptimal. However, given the necessity of molecular surveillance of anti-malarial resistance in settings where poor infrastructure, poor economy, lack of educated staff and even surges of political instability remain major obstacles to performing clinical studies, obtaining the necessary data from used RDTs, despite suboptimal output, becomes a feasible, affordable and hence a justifiable method.