gplas: a comprehensive tool for plasmid analysis using short-read graphs.

SUMMARY: Plasmids can horizontally transmit genetic traits, enabling rapid bacterial adaptation to new environments and hosts. Short-read whole-genome sequencing data are often applied to large-scale bacterial comparative genomics projects but the reconstruction of plasmids from these data is facing severe limitations, such as the inability to distinguish plasmids from each other in a bacterial genome. We developed gplas, a new approach to reliably separate plasmid contigs into discrete components using sequence composition, coverage, assembly graph information and network partitioning based on a pruned network of plasmid unitigs. Gplas facilitates the analysis of large numbers of bacterial isolates and allows a detailed analysis of plasmid epidemiology based solely on short-read sequence data. AVAILABILITY AND IMPLEMENTATION: Gplas is written in R, Bash and uses a Snakemake pipeline as a workflow management system. Gplas is available under the GNU General Public License v3.0 at https://gitlab.com/sirarredondo/gplas.git. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Sequence-based epidemiology of an OXA-48 plasmid during a hospital outbreak.

Objectives. A large OXA-48 outbreak in the Netherlands involved the spread of OXA-48producing Enterobacteriaceae among at least 118 patients, suggesting horizontal transfer of this resistance gene through one or more plasmids. Elucidating transmission dynamics of resistance plasmids is hampered by the low resolution of classic typing methods. This study aimed to investigate the molecular epidemiology of plasmids carrying OXA-48 carbapenemase using a next-generation sequencing approach.Methods. A total of 68 OXA-48-producing Enterobacteriaceae isolated from the hospital outbreak, as well as 22 non-outbreak related OXA-48-producing Enterobacteriaceae from the Netherlands, Libya and Turkey were selected. Plasmids were sequenced using the Illumina Miseq platform, and read sets were assembled and analysed.Results. In all plasmids bla OXA-48 was embedded in transposon Tn1999.2 and located on a ca. 62 kb IncL/M conjugative plasmid in 14 different species. There were a maximum of 2 SNPs (single nucleotide polymorphisms) between the core sequence alignment of all plasmids. Closely related sequence variants of this plasmid were detected in non-outbreak isolates from the Netherlands and other countries. Thirty-one of 89 OXA-48-producing isolates also harboured bla CTX-M-15, which was not located on the bla OXA-48-carrying plasmid. Sequencing of four plasmids harbouring bla CTX-M15 revealed extensive plasmid heterogeneity.Conclusions. A ca 62 kb plasmid was responsible for the OXA-48 outbreak in a Dutch hospital. Our findings provide strong evidence for both within-host inter-species and between host dissemination of plasmid-based OXA-48 during a nosocomial outbreak. These findings exemplify the complex epidemiology of carbapenemase producing Enterobacteriaceae (CPE).

Prevalence, risk factors and genetic characterisation of extended-spectrum beta-lactamase and carbapenemase-producing Enterobacteriaceae (ESBL-E and CPE): a community-based cross-sectional study, the Netherlands, 2014 to 2016.

BackgroundThe epidemiology of carriage of extended-spectrum beta-lactamase-producing (ESBL-E) and carbapenemase-producing Enterobacteriaceae (CPE) in the general population is unknown.AimIn this observational study, the prevalence and risk factors for intestinal ESBL-E and CPE carriage in the Dutch general population were determined. ESBL-E were characterised.MethodsFrom 2014 to 2016, ca 2,000 residents were invited monthly to complete a questionnaire and provide a faecal sample, which was tested for ESBL-E. The first 1,758 samples were also tested for CPE. Risk factors for ESBL-E carriage were identified by multivariable logistic regression analysis. ESBL-E isolates underwent whole genome sequencing.ResultsOf 47,957 individuals invited, 4,177 (8.7%) completed the questionnaire and provided a faecal sample. ESBL-E were detected in 186 (4.5%) individuals, resulting in an adjusted prevalence of 5.0% (95% confidence interval (CI):3.4-6.6%). Risk factors were: born outside the Netherlands (odds ratio (OR): 1.99; 95% CI: 1.16-4.54), eating in restaurants > 20 times/year (OR: 1.70; 95% CI: 1.04-2.76), antibiotic use < 6 months ago (OR: 2.05; 95% CI: 1.05-4.03), swimming in sea/ocean < 12 months ago (OR: 1.63; 95% CI: 1.11-2.39), travelling to Africa (OR: 3.03; 95% CI: 1.23-7.46) or Asia (OR: 2.00; 95% CI: 1.02-3.90) < 12 months ago, and not changing kitchen towels daily (OR: 2.19; 95% CI: 1.24-3.87). The last had the largest population attributable risk (PAR) (47.5%). Eighty-four of 189 (44.4%) ESBL-E isolates carried bla CTX-M-15. Escherichia coli isolates belonged to 70 different sequence types (ST)s, of which ST131 (42/178 isolates; 23.6%) was most prevalent. Associations were observed between IncFIA plasmids and ST131 and bla CTX-M-27, and between IncI1 and ST88 and bla CTX-M-1. No CPE were detected.ConclusionsThe prevalence of ESBL-E carriage in the Netherlands' community-dwelling population is 5.0%. Identified risk factors were mostly travelling (particularly to Asia and Africa) and kitchen hygiene. CPE were not detected.

CRISPR/Cas9-Mediated Genome Editing of Herpesviruses Limits Productive and Latent Infections.

Herpesviruses infect the majority of the human population and can cause significant morbidity and mortality. Herpes simplex virus (HSV) type 1 causes cold sores and herpes simplex keratitis, whereas HSV-2 is responsible for genital herpes. Human cytomegalovirus (HCMV) is the most common viral cause of congenital defects and is responsible for serious disease in immuno-compromised individuals. Epstein-Barr virus (EBV) is associated with infectious mononucleosis and a broad range of malignancies, including Burkitt's lymphoma, nasopharyngeal carcinoma, Hodgkin's disease, and post-transplant lymphomas. Herpesviruses persist in their host for life by establishing a latent infection that is interrupted by periodic reactivation events during which replication occurs. Current antiviral drug treatments target the clinical manifestations of this productive stage, but they are ineffective at eliminating these viruses from the infected host. Here, we set out to combat both productive and latent herpesvirus infections by exploiting the CRISPR/Cas9 system to target viral genetic elements important for virus fitness. We show effective abrogation of HCMV and HSV-1 replication by targeting gRNAs to essential viral genes. Simultaneous targeting of HSV-1 with multiple gRNAs completely abolished the production of infectious particles from human cells. Using the same approach, EBV can be almost completely cleared from latently infected EBV-transformed human tumor cells. Our studies indicate that the CRISPR/Cas9 system can be effectively targeted to herpesvirus genomes as a potent prophylactic and therapeutic anti-viral strategy that may be used to impair viral replication and clear latent virus infection.

New viruses in idiopathic human diarrhea cases, the Netherlands.

Emerging viral infections can be identified by using a viral metagenomics approach for clinical human material. Diarrhea samples of patients with unexplained gastroenteritis from the Netherlands were analyzed by using viral metagenomics. Novel circular DNA viruses, bufaviruses, and genogroup III picobirnaviruses were identified. These data expand our knowledge of the human virome.

Novel divergent nidovirus in a python with pneumonia.

The order Nidovirales contains large, enveloped viruses with a non-segmented positive-stranded RNA genome. Nidoviruses have been detected in man and various animal species, but, to date, there have been no reports of nidovirus in reptiles. In the present study, we describe the detection, characterization, phylogenetic analyses and disease association of a novel divergent nidovirus in the lung of an Indian python (Python molurus) with necrotizing pneumonia. Characterization of the partial genome (>33 000 nt) of this virus revealed several genetic features that are distinct from other nidoviruses, including a very large polyprotein 1a, a putative ribosomal frameshift signal that was identical to the frameshift signal of astroviruses and retroviruses and an accessory ORF that showed some similarity with the haemagglutinin-neuraminidase of paramyxoviruses. Analysis of genome organization and phylogenetic analysis of polyprotein 1ab suggests that this virus belongs to the subfamily Torovirinae. Results of this study provide novel insights into the genetic diversity within the order Nidovirales.

Exploring the potential of next-generation sequencing in detection of respiratory viruses.

Efficient detection of human respiratory viral pathogens is crucial in the management of patients with acute respiratory tract infection. Sequence-independent amplification of nucleic acids combined with next-generation sequencing technology and bioinformatics analyses is a promising strategy for identifying pathogens in clinical and public health settings. It allows the characterization of hundreds of different known pathogens simultaneously and of novel pathogens that elude conventional testing. However, major hurdles for its routine use exist, including cost, turnaround time, and especially sensitivity of the assay, as the detection limit is dependent on viral load, host genetic material, and sequencing depth. To obtain insights into these aspects, we analyzed nasopharyngeal aspirates from a cohort of 81 Thai children with respiratory disease for the presence of respiratory viruses using a sequence-independent next-generation sequencing approach and routinely used diagnostic real-time reverse transcriptase PCR (real-time RT-PCR) assays. With respect to the detection of rhinovirus and human metapneumovirus, the next-generation sequencing approach was at least as sensitive as diagnostic real-time RT-PCR in this small cohort, whereas for bocavirus and enterovirus, next-generation sequencing was less sensitive than real-time RT-PCR. The advantage of the sequencing approach over real-time RT-PCR was the immediate availability of virus-typing information. Considering the development of platforms capable of generating more output data at declining costs, next-generation sequencing remains of interest for future virus diagnosis in clinical and public health settings and certainly as an additional tool when screening results from real-time RT-PCR are negative.

Optimization of standard in-house 24-locus variable-number tandem-repeat typing for Mycobacterium tuberculosis and its direct application to clinical material.

Variable-number tandem-repeat (VNTR) typing with a panel of 24 loci is the current gold standard in the molecular typing of Mycobacterium tuberculosis complex isolates. However, because of technical problems, a part of the loci often cannot be amplified by multiplex PCRs. Therefore, a considerable number of single-locus PCRs have to be performed for the loci with missing results, which impairs the laboratory work flow. Therefore, the original in-house method described by Supply et al. in 2006 was reevaluated. We modified seven primers and the PCR master mixture and obtained a strongly optimized in-house 24-locus VNTR typing method. The percentage of instantly complete 24-locus VNTR patterns detected in the routine flow of typing activities increased to 84.7% from the 72.3% obtained with the typing conducted with the commercially available Genoscreen MIRU-VNTR typing kit. The analytical sensitivity of the optimized in-house method was assessed by serial dilutions of M. tuberculosis in bronchoalveolar lavage fluid. A 1:10 dilution of the different strains tested was the lowest dilution for the detection of a complete 24-locus VNTR pattern. The optimized in-house 24-locus VNTR typing method will reduce the turnaround time of typing significantly and also the financial burden of these activities.

Metagenomic assembly is the main bottleneck in the identification of mobile genetic elements.

Antimicrobial resistance genes (ARG) are commonly found on acquired mobile genetic elements (MGEs) such as plasmids or transposons. Understanding the spread of resistance genes associated with mobile elements (mARGs) across different hosts and environments requires linking ARGs to the existing mobile reservoir within bacterial communities. However, reconstructing mARGs in metagenomic data from diverse ecosystems poses computational challenges, including genome fragment reconstruction (assembly), high-throughput annotation of MGEs, and identification of their association with ARGs. Recently, several bioinformatics tools have been developed to identify assembled fragments of plasmids, phages, and insertion sequence (IS) elements in metagenomic data. These methods can help in understanding the dissemination of mARGs. To streamline the process of identifying mARGs in multiple samples, we combined these tools in an automated high-throughput open-source pipeline, MetaMobilePicker, that identifies ARGs associated with plasmids, IS elements and phages, starting from short metagenomic sequencing reads. This pipeline was used to identify these three elements on a simplified simulated metagenome dataset, comprising whole genome sequences from seven clinically relevant bacterial species containing 55 ARGs, nine plasmids and five phages. The results demonstrated moderate precision for the identification of plasmids (0.57) and phages (0.71), and moderate sensitivity of identification of IS elements (0.58) and ARGs (0.70). In this study, we aim to assess the main causes of this moderate performance of the MGE prediction tools in a comprehensive manner. We conducted a systematic benchmark, considering metagenomic read coverage, contig length cutoffs and investigating the performance of the classification algorithms. Our analysis revealed that the metagenomic assembly process is the primary bottleneck when linking ARGs to identified MGEs in short-read metagenomics sequencing experiments rather than ARGs and MGEs identification by the different tools.

PlasmidEC and gplas2: an optimized short-read approach to predict and reconstruct antibiotic resistance plasmids in Escherichia coli.

Accurate reconstruction of Escherichia coli antibiotic resistance gene (ARG) plasmids from Illumina sequencing data has proven to be a challenge with current bioinformatic tools. In this work, we present an improved method to reconstruct E. coli plasmids using short reads. We developed plasmidEC, an ensemble classifier that identifies plasmid-derived contigs by combining the output of three different binary classification tools. We showed that plasmidEC is especially suited to classify contigs derived from ARG plasmids with a high recall of 0.941. Additionally, we optimized gplas, a graph-based tool that bins plasmid-predicted contigs into distinct plasmid predictions. Gplas2 is more effective at recovering plasmids with large sequencing coverage variations and can be combined with the output of any binary classifier. The combination of plasmidEC with gplas2 showed a high completeness (median=0.818) and F1-Score (median=0.812) when reconstructing ARG plasmids and exceeded the binning capacity of the reference-based method MOB-suite. In the absence of long-read data, our method offers an excellent alternative to reconstruct ARG plasmids in E. coli.

Novel cyclovirus in human cerebrospinal fluid, Malawi, 2010-2011.

To identify unknown human viruses, we analyzed serum and cerebrospinal fluid samples from patients with unexplained paraplegia from Malawi by using viral metagenomics. A novel cyclovirus species was identified and subsequently found in 15% and 10% of serum and cerebrospinal fluid samples, respectively. These data expand our knowledge of cyclovirus diversity and tropism.

Inferring patient to patient transmission of Mycobacterium tuberculosis from whole genome sequencing data.

BACKGROUND: Mycobacterium tuberculosis is characterised by limited genomic diversity, which makes the application of whole genome sequencing particularly attractive for clinical and epidemiological investigation. However, in order to confidently infer transmission events, an accurate knowledge of the rate of change in the genome over relevant timescales is required. METHODS: We attempted to estimate a molecular clock by sequencing 199 isolates from epidemiologically linked tuberculosis cases, collected in the Netherlands spanning almost 16 years. RESULTS: Multiple analyses support an average mutation rate of ~0.3 SNPs per genome per year. However, all analyses revealed a very high degree of variation around this mean, making the confirmation of links proposed by epidemiology, and inference of novel links, difficult. Despite this, in some cases, the phylogenetic context of other strains provided evidence supporting the confident exclusion of previously inferred epidemiological links. CONCLUSIONS: This in-depth analysis of the molecular clock revealed that it is slow and variable over short time scales, which limits its usefulness in transmission studies. However, the superior resolution of whole genome sequencing can provide the phylogenetic context to allow the confident exclusion of possible transmission events previously inferred via traditional DNA fingerprinting techniques and epidemiological cluster investigation. Despite the slow generation of variation even at the whole genome level we conclude that the investigation of tuberculosis transmission will benefit greatly from routine whole genome sequencing.

Proteogenomic analysis of polymorphisms and gene annotation divergences in prokaryotes using a clustered mass spectrometry-friendly database.

Precise annotation of genes or open reading frames is still a difficult task that results in divergence even for data generated from the same genomic sequence. This has an impact in further proteomic studies, and also compromises the characterization of clinical isolates with many specific genetic variations that may not be represented in the selected database. We recently developed software called multistrain mass spectrometry prokaryotic database builder (MSMSpdbb) that can merge protein databases from several sources and be applied on any prokaryotic organism, in a proteomic-friendly approach. We generated a database for the Mycobacterium tuberculosis complex (using three strains of Mycobacterium bovis and five of M. tuberculosis), and analyzed data collected from two laboratory strains and two clinical isolates of M. tuberculosis. We identified 2561 proteins, of which 24 were present in M. tuberculosis H37Rv samples, but not annotated in the M. tuberculosis H37Rv genome. We were also able to identify 280 nonsynonymous single amino acid polymorphisms and confirm 367 translational start sites. As a proof of concept we applied the database to whole-genome DNA sequencing data of one of the clinical isolates, which allowed the validation of 116 predicted single amino acid polymorphisms and the annotation of 131 N-terminal start sites. Moreover we identified regions not present in the original M. tuberculosis H37Rv sequence, indicating strain divergence or errors in the reference sequence. In conclusion, we demonstrated the potential of using a merged database to better characterize laboratory or clinical bacterial strains.

Mge-cluster: a reference-free approach for typing bacterial plasmids.

Extrachromosomal elements of bacterial cells such as plasmids are notorious for their importance in evolution and adaptation to changing ecology. However, high-resolution population-wide analysis of plasmids has only become accessible recently with the advent of scalable long-read sequencing technology. Current typing methods for the classification of plasmids remain limited in their scope which motivated us to develop a computationally efficient approach to simultaneously recognize novel types and classify plasmids into previously identified groups. Here, we introduce mge-cluster that can easily handle thousands of input sequences which are compressed using a unitig representation in a de Bruijn graph. Our approach offers a faster runtime than existing algorithms, with moderate memory usage, and enables an intuitive visualization, classification and clustering scheme that users can explore interactively within a single framework. Mge-cluster platform for plasmid analysis can be easily distributed and replicated, enabling a consistent labelling of plasmids across past, present, and future sequence collections. We underscore the advantages of our approach by analysing a population-wide plasmid data set obtained from the opportunistic pathogen Escherichia coli, studying the prevalence of the colistin resistance gene mcr-1.1 within the plasmid population, and describing an instance of resistance plasmid transmission within a hospital environment.

Characterization of blaKPC-2 and blaNDM-1 Plasmids of a K. pneumoniae ST11 Outbreak Clone.

The most common resistance mechanism to carbapenems is the production of carbapenemases. In 2021, the Pan American Health Organization warned of the emergence and increase in new carbapenemase combinations in Enterobacterales in Latin America. In this study, we characterized four Klebsiella pneumoniae isolates harboring blaKPC and blaNDM from an outbreak during the COVID-19 pandemic in a Brazilian hospital. We assessed their plasmids' transference ability, fitness effects, and relative copy number in different hosts. The K. pneumoniae BHKPC93 and BHKPC104 strains were selected for whole genome sequencing (WGS) based on their pulsed-field gel electrophoresis profile. The WGS revealed that both isolates belong to ST11, and 20 resistance genes were identified in each isolate, including blaKPC-2 and blaNDM-1. The blaKPC gene was present on a ~56 Kbp IncN plasmid and the blaNDM-1 gene on a ~102 Kbp IncC plasmid, along with five other resistance genes. Although the blaNDM plasmid contained genes for conjugational transfer, only the blaKPC plasmid conjugated to E. coli J53, without apparent fitness effects. The minimum inhibitory concentrations (MICs) of meropenem/imipenem against BHKPC93 and BHKPC104 were 128/64 and 256/128 mg/L, respectively. Although the meropenem and imipenem MICs against E. coli J53 transconjugants carrying the blaKPC gene were 2 mg/L, this was a substantial increment in the MIC relative to the original J53 strain. The blaKPC plasmid copy number was higher in K. pneumoniae BHKPC93 and BHKPC104 than in E. coli and higher than that of the blaNDM plasmids. In conclusion, two ST11 K. pneumoniae isolates that were part of a hospital outbreak co-harbored blaKPC-2 and blaNDM-1. The blaKPC-harboring IncN plasmid has been circulating in this hospital since at least 2015, and its high copy number might have contributed to the conjugative transfer of this particular plasmid to an E. coli host. The observation that the blaKPC-containing plasmid had a lower copy number in this E. coli strain may explain why this plasmid did not confer phenotypic resistance against meropenem and imipenem.

Preferential deletion events in the direct repeat locus of Mycobacterium tuberculosis.

The "Harlingen" IS6110 restriction fragment length polymorphism (RFLP) cluster has linked over 100 tuberculosis cases in The Netherlands since 1993. Four Mycobacterium tuberculosis isolates that were epidemiologically linked to this cluster had different spoligotype patterns, as well as slightly divergent IS6110 profiles, compared to the majority of the isolates. Sequencing of the direct repeat (DR) locus revealed sequence polymorphisms at the putative deletion sites. These deletion footprints provided evidence for independent deletions of the central region of the DR locus in three isolates, while the different genotype of the fourth isolate was explained by transmission. Our finding suggests that convergent deletions in the DR locus occur frequently. However, deletion footprints are not suitable to detect convergent deletions in the DR because they seem to be exceptional. Deletion footprints in the DR were not described previously, and we did not observe them in any public M. tuberculosis complex sequences. We conclude that preferential deletions in the DR loci of closely related strains are usually an unnoted event that interferes with clustering of closely related strains.

Mode and dynamics of vanA-type vancomycin resistance dissemination in Dutch hospitals.

BACKGROUND: Enterococcus faecium is a commensal of the gastrointestinal tract of animals and humans but also a causative agent of hospital-acquired infections. Resistance against glycopeptides and to vancomycin has motivated the inclusion of E. faecium in the WHO global priority list. Vancomycin resistance can be conferred by the vanA gene cluster on the transposon Tn1546, which is frequently present in plasmids. The vanA gene cluster can be disseminated clonally but also horizontally either by plasmid dissemination or by Tn1546 transposition between different genomic locations. METHODS: We performed a retrospective study of the genomic epidemiology of 309 vancomycin-resistant E. faecium (VRE) isolates across 32 Dutch hospitals (2012-2015). Genomic information regarding clonality and Tn1546 characterization was extracted using hierBAPS sequence clusters (SC) and TETyper, respectively. Plasmids were predicted using gplas in combination with a network approach based on shared k-mer content. Next, we conducted a pairwise comparison between isolates sharing a potential epidemiological link to elucidate whether clonal, plasmid, or Tn1546 spread accounted for vanA-type resistance dissemination. RESULTS: On average, we estimated that 59% of VRE cases with a potential epidemiological link were unrelated which was defined as VRE pairs with a distinct Tn1546 variant. Clonal dissemination accounted for 32% cases in which the same SC and Tn1546 variants were identified. Horizontal plasmid dissemination accounted for 7% of VRE cases, in which we observed VRE pairs belonging to a distinct SC but carrying an identical plasmid and Tn1546 variant. In 2% of cases, we observed the same Tn1546 variant in distinct SC and plasmid types which could be explained by mixed and consecutive events of clonal and plasmid dissemination. CONCLUSIONS: In related VRE cases, the dissemination of the vanA gene cluster in Dutch hospitals between 2012 and 2015 was dominated by clonal spread. However, we also identified outbreak settings with high frequencies of plasmid dissemination in which the spread of resistance was mainly driven by horizontal gene transfer (HGT). This study demonstrates the feasibility of distinguishing between modes of dissemination with short-read data and provides a novel assessment to estimate the relative contribution of nested genomic elements in the dissemination of vanA-type resistance.

A high-throughput multiplexing and selection strategy to complete bacterial genomes.

BACKGROUND: Bacterial whole-genome sequencing based on short-read technologies often results in a draft assembly formed by contiguous sequences. The introduction of long-read sequencing technologies permits those contiguous sequences to be unambiguously bridged into complete genomes. However, the elevated costs associated with long-read sequencing frequently limit the number of bacterial isolates that can be long-read sequenced. Here we evaluated the recently released 96 barcoding kit from Oxford Nanopore Technologies (ONT) to generate complete genomes on a high-throughput basis. In addition, we propose an isolate selection strategy that optimizes a representative selection of isolates for long-read sequencing considering as input large-scale bacterial collections. RESULTS: Despite an uneven distribution of long reads per barcode, near-complete chromosomal sequences (assembly contiguity = 0.89) were generated for 96 Escherichia coli isolates with associated short-read sequencing data. The assembly contiguity of the plasmid replicons was even higher (0.98), which indicated the suitability of the multiplexing strategy for studies focused on resolving plasmid sequences. We benchmarked hybrid and ONT-only assemblies and showed that the combination of ONT sequencing data with short-read sequencing data is still highly desirable (i) to perform an unbiased selection of isolates for long-read sequencing, (ii) to achieve an optimal genome accuracy and completeness, and (iii) to include small plasmids underrepresented in the ONT library. CONCLUSIONS: The proposed long-read isolate selection ensures the completion of bacterial genomes that span the genome diversity inherent in large collections of bacterial isolates. We show the potential of using this multiplexing approach to close bacterial genomes on a high-throughput basis.

Recovering Escherichia coli Plasmids in the Absence of Long-Read Sequencing Data.

The incidence of infections caused by multidrug-resistant E. coli strains has risen in the past years. Antibiotic resistance in E. coli is often mediated by acquisition and maintenance of plasmids. The study of E. coli plasmid epidemiology and genomics often requires long-read sequencing information, but recently a number of tools that allow plasmid prediction from short-read data have been developed. Here, we reviewed 25 available plasmid prediction tools and categorized them into binary plasmid/chromosome classification tools and plasmid reconstruction tools. We benchmarked six tools (MOB-suite, plasmidSPAdes, gplas, FishingForPlasmids, HyAsP and SCAPP) that aim to reliably reconstruct distinct plasmids, with a special focus on plasmids carrying antibiotic resistance genes (ARGs) such as extended-spectrum beta-lactamase genes. We found that two thirds (n = 425, 66.3%) of all plasmids were correctly reconstructed by at least one of the six tools, with a range of 92 (14.58%) to 317 (50.23%) correctly predicted plasmids. However, the majority of plasmids that carried antibiotic resistance genes (n = 85, 57.8%) could not be completely recovered as distinct plasmids by any of the tools. MOB-suite was the only tool that was able to correctly reconstruct the majority of plasmids (n = 317, 50.23%), and performed best at reconstructing large plasmids (n = 166, 46.37%) and ARG-plasmids (n = 41, 27.9%), but predictions frequently contained chromosome contamination (40%). In contrast, plasmidSPAdes reconstructed the highest fraction of plasmids smaller than 18 kbp (n = 168, 61.54%). Large ARG-plasmids, however, were frequently merged with sequences derived from distinct replicons. Available bioinformatic tools can provide valuable insight into E. coli plasmids, but also have important limitations. This work will serve as a guideline for selecting the most appropriate plasmid reconstruction tool for studies focusing on E. coli plasmids in the absence of long-read sequencing data.

Apparent nosocomial adaptation of Enterococcus faecalis predates the modern hospital era.

Enterococcus faecalis is a commensal and nosocomial pathogen, which is also ubiquitous in animals and insects, representing a classical generalist microorganism. Here, we study E. faecalis isolates ranging from the pre-antibiotic era in 1936 up to 2018, covering a large set of host species including wild birds, mammals, healthy humans, and hospitalised patients. We sequence the bacterial genomes using short- and long-read techniques, and identify multiple extant hospital-associated lineages, with last common ancestors dating back as far as the 19th century. We find a population cohesively connected through homologous recombination, a metabolic flexibility despite a small genome size, and a stable large core genome. Our findings indicate that the apparent hospital adaptations found in hospital-associated E. faecalis lineages likely predate the "modern hospital" era, suggesting selection in another niche, and underlining the generalist nature of this nosocomial pathogen.