Elucidation of host diversity of the VanD-carrying genomic islands in enterococci and anaerobes

Genomic investigation of the emergence of vanD vancomycin-resistant Enterococcus faecium

Vancomycin-resistant Enterococcus (VRE) is an increasingly identified cause of human disease, with most infections resulting from the vanA and vanB genotypes; less is known about other clinically relevant genotypes. Here we report a genomic exploration of a vanD VRE faecium (VREfm), which arose de novo during a single infectious episode. The genomes of the vancomycin-susceptible E. faecium (VSEfm) recipient and resulting VREfm were subjected to long-read sequencing and closed, with whole-genome alignments, cross-mapping and orthologue clustering used to identify genomic variation. Three key differences were identified. (i) The VREfm chromosome gained a 142.6 kb integrative conjugative element (ICE) harbouring the vanD locus. (ii) The native ligase (ddl) was disrupted by an ISEfm1 insertion. (iii) A large 1.74 Mb chromosomal inversion of unknown consequence occurred. Alignment and phylogenetic-based comparisons of the VREfm with a global collection of vanD-harbouring genomes identified strong similarities in the 120-160 kb genomic region surrounding vanD, suggestive of a common mobile element and integration site, irrespective of the diverse taxonomic, geographical and host origins of the isolates. This isolate diversity revealed that this putative ICE (and its source) is globally disseminated and is capable of being acquired by different genera. Although the incidence of vanD VREfm is low, understanding its emergence and potential for spread is crucial for the ongoing efforts to reduce antimicrobial resistance.

Noise reduction strategies in metagenomic chromosome confirmation capture to link antibiotic resistance genes to microbial hosts

The gut microbiota is a reservoir for antimicrobial resistance genes (ARGs). With current sequencing methods, it is difficult to assign ARGs to their microbial hosts, particularly if these ARGs are located on plasmids. Metagenomic chromosome conformation capture approaches (meta3C and Hi-C) have recently been developed to link bacterial genes to phylogenetic markers, thus potentially allowing the assignment of ARGs to their hosts on a microbiome-wide scale. Here, we generated a meta3C dataset of a human stool sample and used previously published meta3C and Hi-C datasets to investigate bacterial hosts of ARGs in the human gut microbiome. Sequence reads mapping to repetitive elements were found to cause problematic noise in, and may importantly skew interpretation of, meta3C and Hi-C data. We provide a strategy to improve the signal-to-noise ratio by discarding reads that map to insertion sequence elements and to the end of contigs. We also show the importance of using spike-in controls to quantify whether the cross-linking step in meta3C and Hi-C protocols has been successful. After filtering to remove artefactual links, 87 ARGs were assigned to their bacterial hosts across all datasets, including 27 ARGs in the meta3C dataset we generated. We show that commensal gut bacteria are an important reservoir for ARGs, with genes coding for aminoglycoside and tetracycline resistance being widespread in anaerobic commensals of the human gut.