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.

Results of the 2020 Genomic Proficiency Test for the network of European Union Reference Laboratory for Antimicrobial Resistance assessing whole-genome-sequencing capacities.

The global surveillance and outbreak investigation of antimicrobial resistance (AMR) is amidst a paradigm shift from traditional biology to bioinformatics. This is due to developments in whole-genome-sequencing (WGS) technologies, bioinformatics tools, and reduced costs. The increased use of WGS is accompanied by challenges such as standardization, quality control (QC), and data sharing. Thus, there is global need for inter-laboratory WGS proficiency test (PT) schemes to evaluate laboratories' capacity to produce reliable genomic data. Here, we present the results of the first iteration of the Genomic PT (GPT) organized by the Global Capacity Building Group at the Technical University of Denmark in 2020. Participating laboratories sequenced two isolates and corresponding DNA of Salmonella enterica, Escherichia coli and Campylobacter coli, using WGS methodologies routinely employed at their laboratories. The participants' ability to obtain consistently good-quality WGS data was assessed based on several QC WGS metrics. A total of 21 laboratories from 21 European countries submitted WGS and meta-data. Most delivered high-quality sequence data with only two laboratories identified as overall underperforming. The QC metrics, N50 and number of contigs, were identified as good indicators for high-sequencing quality. We propose QC thresholds for N50 greater than 20 000 and 25 000 for Campylobacter coli and Escherichia coli, respectively, and number of contigs >200 bp greater than 225, 265 and 100 for Salmonella enterica, Escherichia coli and Campylobacter coli, respectively. The GPT2020 results confirm the importance of systematic QC procedures, ensuring the submission of reliable WGS data for surveillance and outbreak investigation to meet the requirements of the paradigm shift in methodology.

Infant Formula With a Specific Blend of Five Human Milk Oligosaccharides Drives the Gut Microbiota Development and Improves Gut Maturation Markers: A Randomized Controlled Trial.

Background: Human milk oligosaccharides (HMOs) have important biological functions for a healthy development in early life. Objective: This study aimed to investigate gut maturation effects of an infant formula containing five HMOs (2'-fucosyllactose, 2',3-di-fucosyllactose, lacto-N-tetraose, 3'-sialyllactose, and 6'-sialyllactose). Methods: In a multicenter study, healthy infants (7-21 days old) were randomly assigned to a standard cow's milk-based infant formula (control group, CG); the same formula with 1.5 g/L HMOs (test group 1, TG1); or with 2.5 g/L HMOs (test group 2, TG2). A human milk-fed group (HMG) was enrolled as a reference. Fecal samples collected at baseline (n∼150/formula group; HMG n = 60), age 3 (n∼140/formula group; HMG n = 65) and 6 (n∼115/formula group; HMG n = 60) months were analyzed for microbiome (shotgun metagenomics), metabolism, and biomarkers. Results: At both post-baseline visits, weighted UniFrac analysis indicated different microbiota compositions in the two test groups (TGs) compared to CG (P < 0.01) with coordinates closer to that of HMG. The relative abundance of Bifidobacterium longum subsp. infantis (B. infantis) was higher in TGs vs. CG (P < 0.05; except at 6 months: TG2 vs. CG P = 0.083). Bifidobacterium abundance was higher by ∼45% in TGs vs. CG at 6-month approaching HMG. At both post-baseline visits, toxigenic Clostridioides difficile abundance was 75-85% lower in TGs vs. CG (P < 0.05) and comparable with HMG. Fecal pH was significantly lower in TGs vs. CG, and the overall organic acid profile was different in TGs vs. CG, approaching HMG. At 3 months, TGs (vs. CG) had higher secretory immunoglobulin A (sIgA) and lower alpha-1-antitrypsin (P < 0.05). At 6 months, sIgA in TG2 vs. CG remained higher (P < 0.05), and calprotectin was lower in TG1 (P < 0.05) vs. CG. Conclusion: Infant formula with a specific blend of five HMOs supports the development of the intestinal immune system and gut barrier function and shifts the gut microbiome closer to that of breastfed infants with higher bifidobacteria, particularly B. infantis, and lower toxigenic Clostridioides difficile. Clinical Trial Registration: [https://clinicaltrials.gov/ct2/show/], identifier [NCT03722550].

Microbiome Compositions and Resistome Levels after Antibiotic Treatment of Critically Ill Patients: An Observational Cohort Study.

Hospitalization and treatment with antibiotics increase the risk of acquiring multidrug-resistant bacteria due to antibiotic-mediated changes in patient microbiota. This study aimed to investigate how broad- and narrow-spectrum antibiotics affect the gut microbiome and the resistome in antibiotic naïve patients during neurointensive care. Patients admitted to the neurointensive care unit were treated with broad-spectrum (meropenem or piperacillin/tazobactam) or narrow-spectrum antibiotic treatment (including ciprofloxacin, cefuroxime, vancomycin and dicloxacillin) according to clinical indications. A rectal swab was collected from each patient before and after 5-7 days of antibiotic therapy (N = 34), respectively. Shotgun metagenomic sequencing was performed and the composition of metagenomic species (MGS) was determined. The resistome was characterized with CARD RGI software and the CARD database. As a measure for selection pressure in the patient, we used the sum of the number of days with each antibiotic (antibiotic days). We observed a significant increase in richness and a tendency for an increase in the Shannon index after narrow-spectrum treatment. For broad-spectrum treatment the effect was more diverse, with some patients increasing and some decreasing in richness and Shannon index. This was studied further by comparison of patients who had gained or lost >10 MGS, respectively. Selection pressure was significantly higher in patients with decreased richness and a decreased Shannon index who received the broad treatment. A decrease in MGS richness was significantly correlated to the number of drugs administered and the selection pressure in the patient. Bray-Curtis dissimilarities were significant between the pre- and post-treatment of samples in the narrow group, indicating that the longer the narrow-spectrum treatment, the higher the differences between the pre- and the post-treatment microbial composition. We did not find significant differences between pre- and post-treatment for both antibiotic spectrum treatments; however, we observed that most of the antibiotic class resistance genes were higher in abundance in post-treatment after broad-spectrum treatment.

Sialylated human milk oligosaccharides program cognitive development through a non-genomic transmission mode.

Breastmilk contains bioactive molecules essential for brain and cognitive development. While sialylated human milk oligosaccharides (HMOs) have been implicated in phenotypic programming, their selective role and underlying mechanisms remained elusive. Here, we investigated the long-term consequences of a selective lactational deprivation of a specific sialylated HMO in mice. We capitalized on a knock-out (KO) mouse model (B6.129-St6gal1tm2Jxm/J) lacking the gene responsible for the synthesis of sialyl(alpha2,6)lactose (6'SL), one of the two sources of sialic acid (Neu5Ac) to the lactating offspring. Neu5Ac is involved in the formation of brain structures sustaining cognition. To deprive lactating offspring of 6'SL, we cross-fostered newborn wild-type (WT) pups to KO dams, which provide 6'SL-deficient milk. To test whether lactational 6'SL deprivation affects cognitive capabilities in adulthood, we assessed attention, perseveration, and memory. To detail the associated endophenotypes, we investigated hippocampal electrophysiology, plasma metabolomics, and gut microbiota composition. To investigate the underlying molecular mechanisms, we assessed gene expression (at eye-opening and in adulthood) in two brain regions mediating executive functions and memory (hippocampus and prefrontal cortex, PFC). Compared to control mice, WT offspring deprived of 6'SL during lactation exhibited consistent alterations in all cognitive functions addressed, hippocampal electrophysiology, and in pathways regulating the serotonergic system (identified through gut microbiota and plasma metabolomics). These were associated with a site- (PFC) and time-specific (eye-opening) reduced expression of genes involved in central nervous system development. Our data suggest that 6'SL in maternal milk adjusts cognitive development through a short-term upregulation of genes modulating neuronal patterning in the PFC.

Addressing Learning Needs on the Use of Metagenomics in Antimicrobial Resistance Surveillance.

One Health surveillance of antimicrobial resistance (AMR) depends on a harmonized method for detection of AMR. Metagenomics-based surveillance offers the possibility to compare resistomes within and between different target populations. Its potential to be embedded into policy in the future calls for a timely and integrated knowledge dissemination strategy. We developed a blended training (e-learning and a workshop) on the use of metagenomics in surveillance of pathogens and AMR. The objectives were to highlight the potential of metagenomics in the context of integrated surveillance, to demonstrate its applicability through hands-on training and to raise awareness to bias factors. The target participants included staff of competent authorities responsible for AMR monitoring and academic staff. The training was organized in modules covering the workflow, requirements, benefits and challenges of surveillance by metagenomics. The training had 41 participants. The face-to-face workshop was essential to understand the expectations of the participants about the transition to metagenomics-based surveillance. After revision of the e-learning, we released it as a Massive Open Online Course (MOOC), now available at https://www.coursera.org/learn/metagenomics. This course has run in more than 20 sessions, with more than 3,000 learners enrolled, from more than 120 countries. Blended learning and MOOCs are useful tools to deliver knowledge globally and across disciplines. The released MOOC can be a reference knowledge source for international players in the application of metagenomics in surveillance.

Similar genomic patterns of clinical infective endocarditis and oral isolates of Streptococcus sanguinis and Streptococcus gordonii.

Streptococcus gordonii and Streptococcus sanguinis belong to the Mitis group streptococci, which mostly are commensals in the human oral cavity. Though they are oral commensals, they can escape their niche and cause infective endocarditis, a severe infection with high mortality. Several virulence factors important for the development of infective endocarditis have been described in these two species. However, the background for how the commensal bacteria, in some cases, become pathogenic is still not known. To gain a greater understanding of the mechanisms of the pathogenic potential, we performed a comparative analysis of 38 blood culture strains, S. sanguinis (n = 20) and S. gordonii (n = 18) from patients with verified infective endocarditis, along with 21 publicly available oral isolates from healthy individuals, S. sanguinis (n = 12) and S. gordonii (n = 9). Using whole genome sequencing data of the 59 streptococci genomes, functional profiles were constructed, using protein domain predictions based on the translated genes. These functional profiles were used for clustering, phylogenetics and machine learning. A clear separation could be made between the two species. No clear differences between oral isolates and clinical infective endocarditis isolates were found in any of the 675 translated core-genes. Additionally, random forest-based machine learning and clustering of the pan-genome data as well as amino acid variations in the core-genome could not separate the clinical and oral isolates. A total of 151 different virulence genes was identified in the 59 genomes. Among these homologs of genes important for adhesion and evasion of the immune system were found in all of the strains. Based on the functional profiles and virulence gene content of the genomes, we believe that all analysed strains had the ability to become pathogenic.

Pathogen surveillance in the informal settlement, Kibera, Kenya, using a metagenomics approach.

BACKGROUND: Worldwide, the number of emerging and re-emerging infectious diseases is increasing, highlighting the importance of global disease pathogen surveillance. Traditional population-based methods may fail to capture important events, particularly in settings with limited access to health care, such as urban informal settlements. In such environments, a mixture of surface water runoff and human feces containing pathogenic microorganisms could be used as a surveillance surrogate. METHOD: We conducted a temporal metagenomic analysis of urban sewage from Kibera, an urban informal settlement in Nairobi, Kenya, to detect and quantify bacterial and associated antimicrobial resistance (AMR) determinants, viral and parasitic pathogens. Data were examined in conjunction with data from ongoing clinical infectious disease surveillance. RESULTS: A large variation of read abundances related to bacteria, viruses, and parasites of medical importance, as well as bacterial associated antimicrobial resistance genes over time were detected. Significant increased abundances were observed for a number of bacterial pathogens coinciding with higher abundances of AMR genes. Vibrio cholerae as well as rotavirus A, among other virus peaked in several weeks during the study period whereas Cryptosporidium spp. and Giardia spp, varied more over time. CONCLUSION: The metagenomic surveillance approach for monitoring circulating pathogens in sewage was able to detect putative pathogen and resistance loads in an urban informal settlement. Thus, valuable if generated in real time to serve as a comprehensive infectious disease agent surveillance system with the potential to guide disease prevention and treatment. The approach may lead to a paradigm shift in conducting real-time global genomics-based surveillance in settings with limited access to health care.

Bifidobacteriumbreve Bif195 Protects Against Small-Intestinal Damage Caused by Acetylsalicylic Acid in Healthy Volunteers.

BACKGROUND & AIMS: Enteropathy and small-intestinal ulcers are common adverse effects of nonsteroidal anti-inflammatory drugs such as acetylsalicylic acid (ASA). Safe, cytoprotective strategies are needed to reduce this risk. Specific bifidobacteria might have cytoprotective activities, but little is known about these effects in humans. We used serial video capsule endoscopy (VCE) to assess the efficacy of a specific Bifidobacterium strain in healthy volunteers exposed to ASA. METHODS: We performed a single-site, double-blind, parallel-group, proof-of-concept analysis of 75 heathy volunteers given ASA (300 mg) daily for 6 weeks, from July 31 through October 24, 2017. The participants were randomly assigned (1:1) to groups given oral capsules of Bifidobacterium breve (Bif195) (≥5 × 1010 colony-forming units) or placebo daily for 8 weeks. Small-intestinal damage was analyzed by serial VCE at 6 visits. The area under the curve (AUC) for intestinal damage (Lewis score) and the AUC value for ulcers were the primary and first-ranked secondary end points of the trial, respectively. RESULTS: Efficacy data were obtained from 35 participants given Bif195 and 31 given placebo. The AUC for Lewis score was significantly lower in the Bif195 group (3040 ± 1340 arbitrary units) than the placebo group (4351 ± 3195) (P = .0376). The AUC for ulcer number was significantly lower in the Bif195 group (50.4 ± 53.1 arbitrary units) than in the placebo group (75.2 ± 85.3 arbitrary units) (P = .0258). Twelve adverse events were reported from the Bif195 group and 20 from the placebo group. None of the events was determined to be related to Bif195 intake. CONCLUSIONS: In a randomized, double-blind trial of healthy volunteers, we found oral Bif195 to safely reduce the risk of small-intestinal enteropathy caused by ASA. ClinicalTrials.gov no: NCT03228589.

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.

Global monitoring of antimicrobial resistance based on metagenomics analyses of urban sewage.

Antimicrobial resistance (AMR) is a serious threat to global public health, but obtaining representative data on AMR for healthy human populations is difficult. Here, we use metagenomic analysis of untreated sewage to characterize the bacterial resistome from 79 sites in 60 countries. We find systematic differences in abundance and diversity of AMR genes between Europe/North-America/Oceania and Africa/Asia/South-America. Antimicrobial use data and bacterial taxonomy only explains a minor part of the AMR variation that we observe. We find no evidence for cross-selection between antimicrobial classes, or for effect of air travel between sites. However, AMR gene abundance strongly correlates with socio-economic, health and environmental factors, which we use to predict AMR gene abundances in all countries in the world. Our findings suggest that global AMR gene diversity and abundance vary by region, and that improving sanitation and health could potentially limit the global burden of AMR. We propose metagenomic analysis of sewage as an ethically acceptable and economically feasible approach for continuous global surveillance and prediction of AMR.

Metagenomic analysis of viruses in toilet waste from long distance flights-A new procedure for global infectious disease surveillance.

Human viral pathogens are a major public health threat. Reliable information that accurately describes and characterizes the global occurrence and transmission of human viruses is essential to support national and global priority setting, public health actions, and treatment decisions. However, large areas of the globe are currently without surveillance due to limited health care infrastructure and lack of international cooperation. We propose a novel surveillance strategy, using metagenomic analysis of toilet material from international air flights as a method for worldwide viral disease surveillance. The aim of this study was to design, implement, and evaluate a method for viral analysis of airplane toilet waste enabling simultaneous detection and quantification of a wide range of human viral pathogens. Toilet waste from 19 international airplanes was analyzed for viral content, using viral capture probes followed by high-throughput sequencing. Numerous human pathogens were detected including enteric and respiratory viruses. Several geographic trends were observed with samples originating from South Asia having significantly higher viral species richness as well as higher abundances of salivirus A, aichivirus A and enterovirus B, compared to samples originating from North Asia and North America. In addition, certain city specific trends were observed, including high numbers of rotaviruses in airplanes departing from Islamabad. Based on this study we believe that central sampling and analysis at international airports could be a useful supplement for global viral surveillance, valuable for outbreak detection and for guiding public health resources.

Author Correction: Abundance and diversity of the faecal resistome in slaughter pigs and broilers in nine European countries.

In the version of this Article originally published, the surname of author Oksana Lukjancenko was spelt incorrectly as 'Lukjacenko'. This has now been corrected.

Abundance and diversity of the faecal resistome in slaughter pigs and broilers in nine European countries.

Antimicrobial resistance (AMR) in bacteria and associated human morbidity and mortality is increasing. The use of antimicrobials in livestock selects for AMR that can subsequently be transferred to humans. This flow of AMR between reservoirs demands surveillance in livestock and in humans. We quantified and characterized the acquired resistance gene pools (resistomes) of 181 pig and 178 poultry farms from nine European countries, sequencing more than 5,000 Gb of DNA using shotgun metagenomics. We quantified acquired AMR using the ResFinder database and a second database constructed for this study, consisting of AMR genes identified through screening environmental DNA. The pig and poultry resistomes were very different in abundance and composition. There was a significant country effect on the resistomes, more so in pigs than in poultry. We found higher AMR loads in pigs, whereas poultry resistomes were more diverse. We detected several recently described, critical AMR genes, including mcr-1 and optrA, the abundance of which differed both between host species and between countries. We found that the total acquired AMR level was associated with the overall country-specific antimicrobial usage in livestock and that countries with comparable usage patterns had similar resistomes. However, functionally determined AMR genes were not associated with total drug use.

Hospital Epidemiology of Methicillin-Resistant Staphylococcus aureus in a Tertiary Care Hospital in Moshi, Tanzania, as Determined by Whole Genome Sequencing.

OBJECTIVE: To determine molecular epidemiology of methicillin-resistant S. aureus in Tanzania using whole genome sequencing. METHODS: DNA from 33 Staphylococcus species was recovered from subcultured archived Staphylococcus isolates. Whole genome sequencing was performed on Illumina Miseq using paired-end 2 × 250 bp protocol. Raw sequence data were analyzed using online tools. RESULTS: Full susceptibility to vancomycin and chloramphenicol was observed. Thirteen isolates (43.3%) resisted cefoxitin and other antimicrobials tested. Multilocus sequence typing revealed 13 different sequence types among the 30 S. aureus isolates, with ST-8 (n = seven, 23%) being the most common. Gene detection in S. aureus stains were as follows: mecA, 10 (33.3%); pvl, 5 (16.7%); tst, 2 (6.7%). The SNP difference among the six Tanzanian ST-8 MRSA isolates ranged from 24 to 196 SNPs and from 16 to 446 SNPs when using the USA300_FPR3757 or the USA500_2395 as a reference, respectively. The mutation rate was 1.38 × 10-11 SNPs/site/year or 1.4 × 10-6 SNPs/site/year as estimated by USA300_FPR3757 or the USA500_2395, respectively. CONCLUSION: S. aureus isolates causing infections in hospitalized patients in Moshi are highly diverse and epidemiologically unrelated. Temporal phylogenetic analysis provided better resolution on transmission and introduction of MRSA and it may be important to include this in future routines.

In silico assessment of virulence factors in strains of Streptococcus oralis and Streptococcus mitis isolated from patients with Infective Endocarditis.

Purpose. Streptococcus oralis and Streptococcus mitis belong to the Mitis group, which are mostly commensals in the human oral cavity. Even though S. oralis and S. mitis are oral commensals, they can be opportunistic pathogens causing infective endocarditis. A recent taxonomic re-evaluation of the Mitis group has embedded the species Streptococcus tigurinus and Streptococcus dentisani into the species S. oralis as subspecies. In this study, the distribution of virulence factors that contribute to bacterial immune evasion, colonization and adhesion was assessed in clinical strains of S. oralis (subsp. oralis, subsp. tigurinus and subsp. dentisani) and S. mitis. Methodology. Forty clinical S. oralis (subsp. oralis, subsp. dentisani and subsp. tigurinus) and S. mitis genomes were annotated with the pipeline PanFunPro and aligned against the VFDB database for assessment of virulence factors.Results/Key findings. Three homologues of pavA, psaA and lmb, encoding adhesion proteins, were present in all strains. Seven homologues of nanA, nanB, ply, lytA, lytB, lytC and iga, of importance regarding survival in blood and modulation of the human immune system, were variously present in the genomes. Few S. oralis subspecies specific differences were observed. iga homologues were identified in S. oralis subsp. oralis, whereas lytA homologues were identified in S. oralis subsp. oralis and subsp. tigurinus. Conclusion. Differences in the presence of virulence factors among the three S. oralis subspecies were observed. The virulence gene profiles of the 40 S. mitis and S. oralis (subsp. oralis, subsp. dentisani and subsp. tigurinus) contribute with important new knowledge regarding these species and new subspecies.

Correction: MGmapper: Reference based mapping and taxonomy annotation of metagenomics sequence reads.

[This corrects the article DOI: 10.1371/journal.pone.0176469.].

MGmapper: Reference based mapping and taxonomy annotation of metagenomics sequence reads.

An increasing amount of species and gene identification studies rely on the use of next generation sequence analysis of either single isolate or metagenomics samples. Several methods are available to perform taxonomic annotations and a previous metagenomics benchmark study has shown that a vast number of false positive species annotations are a problem unless thresholds or post-processing are applied to differentiate between correct and false annotations. MGmapper is a package to process raw next generation sequence data and perform reference based sequence assignment, followed by a post-processing analysis to produce reliable taxonomy annotation at species and strain level resolution. An in-vitro bacterial mock community sample comprised of 8 genuses, 11 species and 12 strains was previously used to benchmark metagenomics classification methods. After applying a post-processing filter, we obtained 100% correct taxonomy assignments at species and genus level. A sensitivity and precision at 75% was obtained for strain level annotations. A comparison between MGmapper and Kraken at species level, shows MGmapper assigns taxonomy at species level using 84.8% of the sequence reads, compared to 70.5% for Kraken and both methods identified all species with no false positives. Extensive read count statistics are provided in plain text and excel sheets for both rejected and accepted taxonomy annotations. The use of custom databases is possible for the command-line version of MGmapper, and the complete pipeline is freely available as a bitbucked package (https://bitbucket.org/genomicepidemiology/mgmapper). A web-version (https://cge.cbs.dtu.dk/services/MGmapper) provides the basic functionality for analysis of small fastq datasets.

A sampling and metagenomic sequencing-based methodology for monitoring antimicrobial resistance in swine herds.

OBJECTIVES: Reliable methods for monitoring antimicrobial resistance (AMR) in livestock and other reservoirs are essential to understand the trends, transmission and importance of agricultural resistance. Quantification of AMR is mostly done using culture-based techniques, but metagenomic read mapping shows promise for quantitative resistance monitoring. METHODS: We evaluated the ability of: (i) MIC determination for Escherichia coli; (ii) cfu counting of E. coli; (iii) cfu counting of aerobic bacteria; and (iv) metagenomic shotgun sequencing to predict expected tetracycline resistance based on known antimicrobial consumption in 10 Danish integrated slaughter pig herds. In addition, we evaluated whether fresh or manure floor samples constitute suitable proxies for intestinal sampling, using cfu counting, qPCR and metagenomic shotgun sequencing. RESULTS: Metagenomic read-mapping outperformed cultivation-based techniques in terms of predicting expected tetracycline resistance based on antimicrobial consumption. Our metagenomic approach had sufficient resolution to detect antimicrobial-induced changes to individual resistance gene abundances. Pen floor manure samples were found to represent rectal samples well when analysed using metagenomics, as they contain the same DNA with the exception of a few contaminating taxa that proliferate in the extraintestinal environment. CONCLUSIONS: We present a workflow, from sampling to interpretation, showing how resistance monitoring can be carried out in swine herds using a metagenomic approach. We propose metagenomic sequencing should be part of routine livestock resistance monitoring programmes and potentially of integrated One Health monitoring in all reservoirs.

Evaluation of Methods for the Concentration and Extraction of Viruses from Sewage in the Context of Metagenomic Sequencing.

Viral sewage metagenomics is a novel field of study used for surveillance, epidemiological studies, and evaluation of waste water treatment efficiency. In raw sewage human waste is mixed with household, industrial and drainage water, and virus particles are, therefore, only found in low concentrations. This necessitates a step of sample concentration to allow for sensitive virus detection. Additionally, viruses harbor a large diversity of both surface and genome structures, which makes universal viral genomic extraction difficult. Current studies have tackled these challenges in many different ways employing a wide range of viral concentration and extraction procedures. However, there is limited knowledge of the efficacy and inherent biases associated with these methods in respect to viral sewage metagenomics, hampering the development of this field. By the use of next generation sequencing this study aimed to evaluate the efficiency of four commonly applied viral concentrations techniques (precipitation with polyethylene glycol, organic flocculation with skim milk, monolithic adsorption filtration and glass wool filtration) and extraction methods (Nucleospin RNA XS, QIAamp Viral RNA Mini Kit, NucliSENS® miniMAG®, or PowerViral® Environmental RNA/DNA Isolation Kit) to determine the viriome in a sewage sample. We found a significant influence of concentration and extraction protocols on the detected viriome. The viral richness was largest in samples extracted with QIAamp Viral RNA Mini Kit or PowerViral® Environmental RNA/DNA Isolation Kit. Highest viral specificity were found in samples concentrated by precipitation with polyethylene glycol or extracted with Nucleospin RNA XS. Detection of viral pathogens depended on the method used. These results contribute to the understanding of method associated biases, within the field of viral sewage metagenomics, making evaluation of the current literature easier and helping with the design of future studies.