plasmidSPAdes: assembling plasmids from whole genome sequencing data

Whole-Genome Sequence of Pseudomonas putida Strain 1312, a Potential Biostimulant Developed for Agriculture

We report the draft genome sequence of strain 1312 of Pseudomonas putida, which could be interesting to develop as a biostimulant for agriculture and soil depollution treatments.

We report the draft genome sequence of strain 1312 of Pseudomonas putida, which could be interesting to develop as a biostimulant for agriculture and soil depollution treatments.

Genomic Alterations of Staphylococcus aureus ATCC 25923 after Prolonged Passage in the Laboratory

Staphylococcus aureus reference strain ATCC 25923 has been maintained for more than a decade in our laboratory. Genomic study revealed that the resulting strain AFIPCBER_B_8.4 has lost a 37-kb genomic fragment of the ATCC 25923 parental strain.

Staphylococcus aureus reference strain ATCC 25923 has been maintained for more than a decade in our laboratory. Genomic study revealed that the resulting strain AFIPCBER_B_8.4 has lost a 37-kb genomic fragment of the ATCC 25923 parental strain. The missing fragment showed sequence similarity to genes of bacteriophage proteins.

Characterization of a colistin-resistant Avian Pathogenic Escherichia coli ST69 isolate recovered from a broiler chicken in Germany

In recent years, several plasmids harbouring genes encoding phosphoethanolamine transferases conferring colistin resistance have been described in multiple Enterobacteriaceae species. Avian Pathogenic E. coli (APEC) causes colibacillosis and is responsible for a considerable proportion of the disease burden in commercial poultry flocks, and may be linked to zoonotic infections in humans. Here, we describe the genotypic and phenotypic characteristics of a multidrug-resistant APEC ST69 isolate (APECA2), recovered in 2016 from a diseased broiler at post-mortem examination in Germany. The isolate was resistant to several antibiotics of human and veterinary importance, including colistin. The mcr-1 gene was detected on a mobile genetic element located on an IncHI2/ST4 plasmid, which was characterized using long-read Nanopore and short-read Illumina sequencing of purified plasmid. Isolate APECA2 displayed resistance to chicken serum and harbours numerous virulence genes. This study highlights the public health importance of enhanced antimicrobial resistance surveillance and strict antimicrobial stewardship in human and veterinary healthcare.

Conjugative Transfer of a Novel Staphylococcal Plasmid Encoding the Biocide Resistance Gene, qacA

Staphylococcus aureus is the leading cause of skin and soft tissue infections (SSTI). Some S. aureus strains harbor plasmids that carry genes that affect resistance to biocides. Among these genes, qacA encodes the QacA Multidrug Efflux Pump that imparts decreased susceptibility to chlorhexidine, a biocide used ubiquitously in healthcare facilities. Furthermore, chlorhexidine has been considered as a S. aureus decolonization strategy in community settings. We previously conducted a chlorhexidine-based SSTI prevention trial among Ft. Benning Army trainees. Analysis of a clinical isolate (C02) from that trial identified a novel qacA-positive plasmid, pC02. Prior characterization of qacA-containing plasmids is limited and conjugative transfer of those plasmids has not been demonstrated. Given the implications of increased biocide resistance, herein we characterized pC02. In silico analysis identified genes typically associated with conjugative plasmids. Moreover, pC02 was efficiently transferred to numerous S. aureus strains and to Staphylococcus epidermidis. We screened additional qacA-positive S. aureus clinical isolates and pC02 was present in 27% of those strains; other unique qacA-harboring plasmids were also identified. Ten strains were subjected to whole genome sequencing. Sequence analysis combined with plasmid screening studies suggest that qacA-containing strains are transmitted among military personnel at Ft. Benning and that strains carrying qacA are associated with SSTIs within this population. The identification of a novel mechanism of qacA conjugative transfer among Staphylococcal strains suggests a possible future increase in the prevalence of antiseptic tolerant bacterial strains, and an increase in the rate of infections in settings where these agents are commonly used.

Exploring the genome of Arctic Psychrobacter sp. DAB_AL32B and construction of novel Psychrobacter-specific cloning vectors of an increased carrying capacity

Cold-active bacteria are currently of great interest in biotechnology, and their genomic and physiological features have been extensively studied. One of the model psychrotolerant bacteria are Psychrobacter spp. Analysis of Arctic psychrophilic Psychrobacter sp. DAB_AL32B genome content provided an insight into its overall stress response, and genes conferring protection against various life-limiting factors (i.e., low temperature, increased ultraviolet radiation, oxidative stress and osmotic pressure) were recognized and described. Moreover, it was revealed that the strain carries a large plasmid pP32BP2. Its replication system was used for the construction of two novel shuttle vectors (pPS-NR-Psychrobacter-Escherichia coli-specific plasmid and pPS-BR-Psychrobacter-various Proteobacteria-specific plasmid) of an increased carrying capacity, which may be used for genetic engineering of Psychrobacter spp.

mlplasmids: a user-friendly tool to predict plasmid- and chromosome-derived sequences for single species

Assembly of bacterial short-read whole-genome sequencing data frequently results in hundreds of contigs for which the origin, plasmid or chromosome, is unclear. Complete genomes resolved by long-read sequencing can be used to generate and label short-read contigs. These were used to train several popular machine learning methods to classify the origin of contigs from Enterococcus faecium, Klebsiella pneumoniae and Escherichia coli using pentamer frequencies. We selected support-vector machine (SVM) models as the best classifier for all three bacterial species (F1-score E. faecium=0.92, F1-score K. pneumoniae=0.90, F1-score E. coli=0.76), which outperformed other existing plasmid prediction tools using a benchmarking set of isolates. We demonstrated the scalability of our models by accurately predicting the plasmidome of a large collection of 1644 E. faecium isolates and illustrate its applicability by predicting the location of antibiotic-resistance genes in all three species. The SVM classifiers are publicly available as an R package and graphical-user interface called 'mlplasmids'. We anticipate that this tool may significantly facilitate research on the dissemination of plasmids encoding antibiotic resistance and/or contributing to host adaptation.

First Report of an OXA-48- and CTX-M-213-Producing Kluyvera Species Clone Recovered from Patients Admitted in a University Hospital in Madrid, Spain

Enterobacterales species other than Klebsiella pneumoniae also contribute to OXA-48 carbapenemase endemicity. We studied the emergence of an OXA-48-producing Kluyvera species clone, which expresses the novel CTX-M-213 enzyme, colonizing patients in our hospital. Rectal swabs from patients admitted in four wards (March 2014 to March 2016; R-GNOSIS project) were seeded onto Chromo ID-ESBL) and Chrom-CARB/OXA-48 chromogenic agar plates. Carbapenemases and extended-spectrum β-lactamases (ESBLs) were characterized (PCR, sequencing, cloning, and site-directed mutagenesis), and antibiotic susceptibility was determined. Clonal relatedness was established (XbaI pulsed-field gel electrophoresis [XbaI-PFGE]), and plasmid content was studied (transformation, S1 nuclease digestion-PFGE, SB-hybridization, restriction fragment length polymorphism [RFLP] analysis [DraI and HpaI], and PCR [incompatibility group and repA, traU, and parA genes]). Whole-genome sequencing (WGS) (Illumina HiSeq-2500) and further bioinformatics analysis of plasmids (PLACNET and plasmidSPAdes) were performed. Patients' charts were reviewed. Six unrelated patients (median age, 75 years [range, 59 to 81 years]; 4/6 male patients) colonized with OXA-48-producing Kluyvera species isolates (>95% similarity of the PFGE pattern) were identified. Nosocomial acquisition was demonstrated. In two patients, OXA-48-producing Kluyvera species isolates coexisted with OXA-48-producing Raoultella ornithinolytica, K. pneumoniae, and Escherichia coli The bla _OXA-48 gene was located on an ∼60-kb IncL plasmid related to IncL/M-pOXA-48a and the novel bla _CTX-M-213 gene in a conserved chromosomal region of Kluyvera species isolates. CTX-M-213, different from CTX-M-13 (K56E) but conferring a similar β-lactam resistance profile, was identified. Genomic analysis also revealed a 177-kb IncF plasmid (class I integron harboring sul1 and aadA2) and an 8-kb IncQ plasmid (IS4-bla _FOX-8). We describe the first bla _OXA-48 plasmid in Kluyvera spp. and the novel chromosomal CTX-M-213 enzyme and highlight further nosocomial dissemination of bla _OXA-48 through clonal lineages or plasmids related to IncL/M-pOXA-48a.

Firefly genomes illuminate parallel origins of bioluminescence in beetles

Fireflies and their luminous courtships have inspired centuries of scientific study. Today firefly luciferase is widely used in biotechnology, but the evolutionary origin of bioluminescence within beetles remains unclear. To shed light on this long-standing question, we sequenced the genomes of two firefly species that diverged over 100 million-years-ago: the North American Photinus pyralis and Japanese Aquatica lateralis. To compare bioluminescent origins, we also sequenced the genome of a related click beetle, the Caribbean Ignelater luminosus, with bioluminescent biochemistry near-identical to fireflies, but anatomically unique light organs, suggesting the intriguing hypothesis of parallel gains of bioluminescence. Our analyses support independent gains of bioluminescence in fireflies and click beetles, and provide new insights into the genes, chemical defenses, and symbionts that evolved alongside their luminous lifestyle.

Genome sequencing and assessment of plant growth-promoting properties of a Serratia marcescens strain isolated from vermicompost

Background

Plant-bacteria associations have been extensively studied for their potential in increasing crop productivity in a sustainable manner. Serratia marcescens is a species of Enterobacteriaceae found in a wide range of environments, including soil.

Results

Here we describe the genome sequencing and assessment of plant growth-promoting abilities of S. marcescens UENF-22GI, a strain isolated from mature cattle manure vermicompost. In vitro, S. marcescens UENF-22GI is able to solubilize P and Zn, to produce indole compounds (likely IAA), to colonize hyphae and counter the growth of two phytopathogenic fungi. Inoculation of maize with this strain remarkably increased seedling growth and biomass under greenhouse conditions. The S. marcescens UENF-22GI genome has 5 Mb, assembled in 17 scaffolds comprising 4662 genes (4528 are protein-coding). No plasmids were identified. S. marcescens UENF-22GI is phylogenetically placed within a clade comprised almost exclusively of non-clinical strains. We identified genes and operons that are likely responsible for the interesting plant-growth promoting features that were experimentally described. The S. marcescens UENF-22GI genome harbors a horizontally-transferred genomic island involved in antibiotic production, antibiotic resistance, and anti-phage defense via a novel ADP-ribosyltransferase-like protein and possible modification of DNA by a deazapurine base, which likely contributes to its competitiveness against other bacteria.

Conclusions

Collectively, our results suggest that S. marcescens UENF-22GI is a strong candidate to be used in the enrichment of substrates for plant growth promotion or as part of bioinoculants for agriculture.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5130-y) contains supplementary material, which is available to authorized users.

Draft Genome Sequence of Xylella fastidiosa subsp. fastidiosa Strain IVIA5235, Isolated from Prunus avium in Mallorca Island, Spain

We report the complete annotated genome sequence of the plant-pathogenic bacterium Xylella fastidiosa subsp. fastidiosa strain IVIA5235.

We report the complete annotated genome sequence of the plant-pathogenic bacterium Xylella fastidiosa subsp. fastidiosa strain IVIA5235. This strain was recovered from a cherry tree in Mallorca, Spain.

Comparative Genomics of the First and Complete Genome of "Actinobacillus porcitonsillarum" Supports the Novel Species Hypothesis

“Actinobacillus porcitonsillarum” is considered a nonpathogenic member of the Pasteurellaceae family, which phenotypically resembles the pathogen Actinobacillus pleuropneumoniae. Previous studies suggested that “A. porcitonsillarum” may represent a new species closely related to Actinobacillus minor, yet no full genome has been sequenced so far. We implemented the Oxford Nanopore and Illumina sequencing technologies to obtain the highly accurate and complete genome sequence of the “A. porcitonsillarum” strain 9953L55. After validating our de novo assembly strategy by comparing the A. pleuropneumoniae S4074^T genome sequence obtained by Oxford Nanopore Technology combined with Illumina reads with a PacBio-sequenced S4074^T genome from the NCBI database, we performed comparative analyses of the 9953L55 genome with the A. minor type strain NM305^T, A. minor strain 202, and A. pleuropneumoniae S4074^T. The 2,263,191 bp circular genome of 9953L55 consisted of 2168 and 2033 predicted genes and proteins, respectively. The lipopolysaccharide cluster resembled the genetic organization of A. pleuropneumoniae serotypes 1, 9, and 11, possibly explaining the positive reactions observed previously in serotyping tests. In contrast to NM305^T, we confirmed the presence of a complete apxIICABD operon in 9953L55 and 202 accounting for their hemolytic phenotype and Christie-Atkins-Munch-Petersen (CAMP) reaction positivity. Orthologous gene cluster analysis provided insight into the differential ability of strains of the A. minor/“porcitonsillarum” complex and A. pleuropneumoniae to ferment lactose, raffinose, trehalose, and mannitol. The four strains showed distinct and shared transposable elements, CRISPR/Cas systems, and integrated prophages. Genome comparisons based on average nucleotide identity and in silico DNA-DNA hybridization confirmed the close relationship among strains belonging to the A. minor/“porcitonsillarum” complex compared to other Actinobacillus spp., but also suggested that 9953L55 and 202 belong to the same novel species closely related to A. minor, namely, “A. porcitonsillarum.” Recognition of the taxon as a separate species would improve diagnostics and control strategies of pig pleuropneumonia.

Whole-Genome Sequence of Pseudomonas aeruginosa Strain 4014, Isolated from Soil in France

We report here the draft genome sequence of strain 4014 of Pseudomonas aeruginosa, a common human pathogen, isolated from soil in France. This sequence predicts resistance to multiple antibiotics, including vancomycin.

We report here the draft genome sequence of strain 4014 of Pseudomonas aeruginosa, a common human pathogen, isolated from soil in France. This sequence predicts resistance to multiple antibiotics, including vancomycin.

Construction and Characterization of Synthetic Bacterial Community for Experimental Ecology and Evolution

Experimental microbial ecology and evolution have yielded foundational insights into ecological and evolutionary processes using simple microcosm setups and phenotypic assays with one- or two-species model systems. The fields are now increasingly incorporating more complex systems and exploration of the molecular basis of observations. For this purpose, simplified, manageable and well-defined multispecies model systems are required that can be easily investigated using culturing and high-throughput sequencing approaches, bridging the gap between simpler and more complex synthetic or natural systems. Here we address this need by constructing a completely synthetic 33 bacterial strain community that can be cultured in simple laboratory conditions. We provide whole-genome data for all the strains as well as metadata about genomic features and phenotypic traits that allow resolving individual strains by amplicon sequencing and facilitate a variety of envisioned mechanistic studies. We further show that a large proportion of the strains exhibit coexistence in co-culture over serial transfer for 48 days in the absence of any experimental manipulation to maintain diversity. The constructed bacterial community can be a valuable resource in future experimental work.

MOB-suite: software tools for clustering, reconstruction and typing of plasmids from draft assemblies

Large-scale bacterial population genetics studies are now routine due to cost-effective Illumina short-read sequencing. However, analysing plasmid content remains difficult due to incomplete assembly of plasmids. Bacterial isolates can contain any number of plasmids and assembly remains complicated due to the presence of repetitive elements. Numerous tools have been developed to analyse plasmids but the performance and functionality of the tools are variable. The MOB-suite was developed as a set of modular tools for reconstruction and typing of plasmids from draft assembly data to facilitate characterization of plasmids. Using a set of closed genomes with publicly available Illumina data, the MOB-suite identified contigs of plasmid origin with both high sensitivity and specificity (95 and 88 %, respectively). In comparison, plasmidfinder demonstrated high specificity (99 %) but limited sensitivity (50 %). Using the same dataset of 377 known plasmids, MOB-recon accurately reconstructed 207 plasmids so that they were assigned to a single grouping without other plasmid or chromosomal sequences, whereas plasmidSPAdes was only able to accurately reconstruct 102 plasmids. In general, plasmidSPAdes has a tendency to merge different plasmids together, with 208 plasmids undergoing merge events. The MOB-suite reduces the number of errors but produces more hybrid plasmids, with 84 plasmids undergoing both splits and merges. The MOB-suite also provides replicon typing similar to plasmidfinder but with the inclusion of relaxase typing and prediction of conjugation potential. The MOB-suite is written in Python 3 and is available from https://github.com/phac-nml/mob-suite.

Genome Variation in the Model Halophilic Bacterium Salinibacter ruber

The halophilic bacterium Salinibacter ruber is an abundant and ecologically important member of halophilic communities worldwide. Given its broad distribution and high intraspecific genetic diversity, S. ruber is considered one of the main models for ecological and evolutionary studies of bacterial adaptation to hypersaline environments. However, current insights on the genomic diversity of this species is limited to the comparison of the genomes of two co-isolated strains. Here, we present a comparative genomic analysis of eight S. ruber strains isolated at two different time points in each of two different Mediterranean solar salterns. Our results show an open pangenome with contrasting evolutionary patterns in the core and accessory genomes. We found that the core genome is shaped by extensive homologous recombination (HR), which results in limited sequence variation within population clusters. In contrast, the accessory genome is modulated by horizontal gene transfer (HGT), with genomic islands and plasmids acting as gateways to the rest of the genome. In addition, both types of genetic exchange are modulated by restriction and modification (RM) or CRISPR-Cas systems. Finally, genes differentially impacted by such processes reveal functional processes potentially relevant for environmental interactions and adaptation to extremophilic conditions. Altogether, our results support scenarios that conciliate “Neutral” and “Constant Diversity” models of bacterial evolution.

Staphylococcus aureus viewed from the perspective of 40,000+ genomes

Low-cost Illumina sequencing of clinically-important bacterial pathogens has generated thousands of publicly available genomic datasets. Analyzing these genomes and extracting relevant information for each pathogen and the associated clinical phenotypes requires not only resources and bioinformatic skills but organism-specific knowledge. In light of these issues, we created Staphopia, an analysis pipeline, database and application programming interface, focused on Staphylococcus aureus, a common colonizer of humans and a major antibiotic-resistant pathogen responsible for a wide spectrum of hospital and community-associated infections. Written in Python, Staphopia's analysis pipeline consists of submodules running open-source tools. It accepts raw FASTQ reads as an input, which undergo quality control filtration, error correction and reduction to a maximum of approximately 100× chromosome coverage. This reduction significantly reduces total runtime without detrimentally affecting the results. The pipeline performs de novo assembly-based and mapping-based analysis. Automated gene calling and annotation is performed on the assembled contigs. Read-mapping is used to call variants (single nucleotide polymorphisms and insertion/deletions) against a reference S. aureus chromosome (N315, ST5). We ran the analysis pipeline on more than 43,000 S. aureus shotgun Illumina genome projects in the public European Nucleotide Archive database in November 2017. We found that only a quarter of known multi-locus sequence types (STs) were represented but the top 10 STs made up 70% of all genomes. methicillin-resistant S. aureus (MRSA) were 64% of all genomes. Using the Staphopia database we selected 380 high quality genomes deposited with good metadata, each from a different multi-locus ST, as a non-redundant diversity set for studying S. aureus evolution. In addition to answering basic science questions, Staphopia could serve as a potential platform for rapid clinical diagnostics of S. aureus isolates in the future. The system could also be adapted as a template for other organism-specific databases.

Heterogeneous and Flexible Transmission of mcr-1 in Hospital-Associated Escherichia coli

The recent emergence of a transferable colistin resistance mechanism, MCR-1, has gained global attention because of its threat to clinical treatment of infections caused by multidrug-resistant Gram-negative bacteria. However, the possible transmission route of mcr-1 among Enterobacteriaceae species in clinical settings is largely unknown. Here, we present a comprehensive genomic analysis of Escherichia coli isolates collected in a hospital in Hangzhou, China. We found that mcr-1-carrying isolates from clinical infections and feces of inpatients and healthy volunteers were genetically diverse and were not closely related phylogenetically, suggesting that clonal expansion is not involved in the spread of mcr-1. The mcr-1 gene was found on either chromosomes or plasmids, but in most of the E. coli isolates, mcr-1 was carried on plasmids. The genetic context of the plasmids showed considerable diversity as evidenced by the different functional insertion sequence (IS) elements, toxin-antitoxin (TA) systems, heavy metal resistance determinants, and Rep proteins of broad-host-range plasmids. Additionally, the genomic analysis revealed nosocomial transmission of mcr-1 and the coexistence of mcr-1 with other genes encoding β-lactamases and fluoroquinolone resistance in the E. coli isolates. These findings indicate that mcr-1 is heterogeneously disseminated in both commensal and pathogenic strains of E. coli, suggest the high flexibility of this gene in its association with diverse genetic backgrounds of the hosts, and provide new insights into the genome epidemiology of mcr-1 among hospital-associated E. coli strains.

Colistin represents one of the very few available drugs for treating infections caused by extensively multidrug-resistant Gram-negative bacteria. The recently emergent mcr-1 colistin resistance gene threatens the clinical utility of colistin and has gained global attention. How mcr-1 spreads in hospital settings remains unknown and was investigated by whole-genome sequencing of mcr-1-carrying Escherichia coli in this study. The findings revealed extraordinary flexibility of mcr-1 in its spread among genetically diverse E. coli hosts and plasmids, nosocomial transmission of mcr-1-carrying E. coli, and the continuous emergence of novel Inc types of plasmids carrying mcr-1 and new mcr-1 variants. Additionally, mcr-1 was found to be frequently associated with other genes encoding β-lactams and fluoroquinolone resistance. These findings provide important information on the transmission and epidemiology of mcr-1 and are of significant public health importance as the information is expected to facilitate the control of this significant antibiotic resistance threat.

Comparative genomics of Campylobacter concisus: Analysis of clinical strains reveals genome diversity and pathogenic potential

In recent years, an increasing number of Campylobacter species have been associated with human gastrointestinal (GI) diseases including gastroenteritis, inflammatory bowel disease, and colorectal cancer. Campylobacter concisus, an oral commensal historically linked to gingivitis and periodontitis, has been increasingly detected in the lower GI tract. In the present study, we generated robust genome sequence data from C. concisus strains and undertook a comprehensive pangenome assessment to identify C. concisus virulence properties and to explain potential adaptations acquired while residing in specific ecological niche(s) of the GI tract. Genomes of 53 new C. concisus strains were sequenced, assembled, and annotated including 36 strains from gastroenteritis patients, 13 strains from Crohn’s disease patients and four strains from colitis patients (three collagenous colitis and one lymphocytic colitis). When compared with previous published sequences, strains clustered into two main groups/genomospecies (GS) with phylogenetic clustering explained neither by disease phenotype nor sample location. Paired oral/faecal isolates, from the same patient, indicated that there are few genetic differences between oral and gut isolates which suggests that gut isolates most likely reflect oral strain relocation. Type IV and VI secretion systems genes, genes known to be important for pathogenicity in the Campylobacter genus, were present in the genomes assemblies, with 82% containing Type VI secretion system genes. Our findings indicate that C. concisus strains are genetically diverse, and the variability in bacterial secretion system content may play an important role in their virulence potential.

Genome analysis of methicillin resistance in Macrococcus caseolyticus from dairy cattle in England and Wales

Species of the genus Macrococcus are widespread commensals of animals but are becoming increasingly recognised as veterinary pathogens. They can encode methicillin resistance and are implicated in its spread to the closely-related, but more pathogenic, staphylococci. In this study we have identified 33 isolates of methicillin-resistant Macrococcus caseolyticus from bovine bulk tank milk from England and Wales. These isolates were characterised to provide insight into the molecular epidemiology of M. caseolyticus and to discern the genetic basis for their methicillin resistance. Antimicrobial susceptibility testing was performed by Vitek2 and disc diffusion. Isolates were whole-genome sequenced to evaluate phylogenetic relationships and the presence of methicillin resistance determinants, mecA–D. All 33 isolates were phenotypically methicillin-resistant according to cefoxitin disc diffusion, cefoxitin Etest and oxacillin resistance assessed by Vitek2. In contrast only a single isolate was resistant in the Vitek2 cefoxitin screen. Twenty-seven isolates were positive for mecD and six were positive for mecB. mecA and mecC were not detected. The results of phylogenetic analysis indicated that these methicillin-resistant isolates represented a heterogeneous population with both mecB and mecD found in diverse isolates. Isolates had a widespread distribution across the sampled region. Taken together with the role of M. caseolyticus in veterinary infections, including bovine mastitis, and in the potential spread of methicillin resistance to more pathogenic staphylococci, this work highlights the need to better understand their epidemiology and for increased awareness among veterinary microbiology laboratories.

First Complete Genome Sequence of Yersinia massiliensis

Using a combination of Illumina paired-end sequencing, Pacific Biosciences RS II sequencing, and OpGen Argus whole-genome optical mapping, we report here the first complete genome sequence of Yersinia massiliensis. The completed genome consists of a 4.99-Mb chromosome, a 121-kb megaplasmid, and a 57-kb plasmid.

The Use of a Combined Bioinformatics Approach to Locate Antibiotic Resistance Genes on Plasmids From Whole Genome Sequences of Salmonella enterica Serovars From Humans in Ghana

In the current study, we identified plasmids carrying antimicrobial resistance genes in draft whole genome sequences of 16 selected Salmonella enterica isolates representing six different serovars from humans in Ghana. The plasmids and the location of resistance genes in the genomes were predicted using a combination of PlasmidFinder, ResFinder, plasmidSPAdes and BLAST genomic analysis tools. Subsequently, S1-PFGE was employed for analysis of plasmid profiles. Whole genome sequencing confirmed the presence of antimicrobial resistance genes in Salmonella isolates showing multidrug resistance phenotypically. ESBL, either blaTEM52-B or blaCTX-M15 were present in two cephalosporin resistant isolates of S. Virchow and S. Poona, respectively. The systematic genome analysis revealed the presence of different plasmids in different serovars, with or without insertion of antimicrobial resistance genes. In S. Enteritidis, resistance genes were carried predominantly on plasmids of IncN type, in S. Typhimurium on plasmids of IncFII(S)/IncFIB(S)/IncQ1 type. In S. Virchow and in S. Poona, resistance genes were detected on plasmids of IncX1 and TrfA/IncHI2/IncHI2A type, respectively. The latter two plasmids were described for the first time in these serovars. The combination of genomic analytical tools allowed nearly full mapping of the resistance plasmids in all Salmonella strains analyzed. The results suggest that the improved analytical approach used in the current study may be used to identify plasmids that are specifically associated with resistance phenotypes in whole genome sequences. Such knowledge would allow the development of rapid multidrug resistance tracking tools in Salmonella populations using WGS.

Short genome report of cellulose-producing commensal Escherichia coli 1094

Bacterial surface colonization and biofilm formation often rely on the production of an extracellular polymeric matrix that mediates cell-cell and cell-surface contacts. In Escherichia coli and many Betaproteobacteria and Gammaproteobacteria cellulose is often the main component of the extracellular matrix. Here we report the complete genome sequence of the cellulose producing strain E. coli 1094 and compare it with five other closely related genomes within E. coli phylogenetic group A. We present a comparative analysis of the regions encoding genes responsible for cellulose biosynthesis and discuss the changes that could have led to the loss of this important adaptive advantage in several E. coli strains. Data deposition: The annotated genome sequence has been deposited at the European Nucleotide Archive under the accession number PRJEB21000.

Acquisition and dissemination of cephalosporin-resistant E. coli in migratory birds sampled at an Alaska landfill as inferred through genomic analysis

Antimicrobial resistance (AMR) in bacterial pathogens threatens global health, though the spread of AMR bacteria and AMR genes between humans, animals, and the environment is still largely unknown. Here, we investigated the role of wild birds in the epidemiology of AMR Escherichia coli. Using next-generation sequencing, we characterized cephalosporin-resistant E. coli cultured from sympatric gulls and bald eagles inhabiting a landfill habitat in Alaska to identify genetic determinants conferring AMR, explore potential transmission pathways of AMR bacteria and genes at this site, and investigate how their genetic diversity compares to isolates reported in other taxa. We found genetically diverse E. coli isolates with sequence types previously associated with human infections and resistance genes of clinical importance, including bla_CTX-M and bla_CMY. Identical resistance profiles were observed in genetically unrelated E. coli isolates from both gulls and bald eagles. Conversely, isolates with indistinguishable core-genomes were found to have different resistance profiles. Our findings support complex epidemiological interactions including bacterial strain sharing between gulls and bald eagles and horizontal gene transfer among E. coli harboured by birds. Results suggest that landfills may serve as a source for AMR acquisition and/or maintenance, including bacterial sequence types and AMR genes relevant to human health.

Enterococcus faecalis AHG0090 is a Genetically Tractable Bacterium and Produces a Secreted Peptidic Bioactive that Suppresses Nuclear Factor Kappa B Activation in Human Gut Epithelial Cells

Enterococcus faecalis is an early coloniser of the human infant gut and contributes to the development of intestinal immunity. To better understand the functional capacity of E. faecalis, we constructed a broad host range RP4 mobilizable vector, pEHR513112, that confers chloramphenicol resistance and used a metaparental mating approach to isolate E. faecalis AHG0090 from a fecal sample collected from a healthy human infant. We demonstrated that E. faecalis AHG0090 is genetically tractable and could be manipulated using traditional molecular microbiology approaches. E. faecalis AHG0090 was comparable to the gold-standard anti-inflammatory bacterium Faecalibacterium prausnitzii A2-165 in its ability to suppress cytokine-mediated nuclear factor kappa B (NF-κB) activation in human gut-derived LS174T goblet cell like and Caco-2 enterocyte-like cell lines. E. faecalis AHG0090 and F. prausnitzii A2-165 produced secreted low molecular weight NF-κB suppressive peptidic bioactives. Both bioactives were sensitive to heat and proteinase K treatments although the E. faecalis AHG0090 bioactive was more resilient to both forms of treatment. As expected, E. faecalis AHG0090 suppressed IL-1β-induced NF-κB-p65 subunit nuclear translocation and expression of the NF-κB regulated genes IL-6, IL-8 and CXCL-10. Finally, we determined that E. faecalis AHG0090 is distantly related to other commensal strains and likely encodes niche factors that support effective colonization of the infant gut.

Draft Genome Sequences of Four Clinical Legionella pneumophila Isolates from Ontario, Canada

Legionella pneumophila outbreak investigations require the development of reliable typing methods to better understand the genetic relationships of the isolates involved. Here, we report the draft genome sequences of four clinical Legionella pneumophila isolates obtained between 2000 and 2012 in Ontario, Canada.

Genomic analysis of oral Campylobacter concisus strains identified a potential bacterial molecular marker associated with active Crohn's disease

Campylobacter concisus is an oral bacterium that is associated with inflammatory bowel disease (IBD) including Crohn's disease (CD) and ulcerative colitis (UC). C. concisus consists of two genomospecies (GS) and diverse strains. This study aimed to identify molecular markers to differentiate commensal and IBD-associated C. concisus strains. The genomes of 63 oral C. concisus strains isolated from patients with IBD and healthy controls were examined, of which 38 genomes were sequenced in this study. We identified a novel secreted enterotoxin B homologue, Csep1. The csep1 gene was found in 56% of GS2 C. concisus strains, presented in the plasmid pICON or the chromosome. A six-nucleotide insertion at the position 654-659 bp in csep1 (csep1-6bpi) was found. The presence of csep1-6bpi in oral C. concisus strains isolated from patients with active CD (47%, 7/15) was significantly higher than that in strains from healthy controls (0/29, P = 0.0002), and the prevalence of csep1-6bpi positive C. concisus strains was significantly higher in patients with active CD (67%, 4/6) as compared to healthy controls (0/23, P = 0.0006). Proteomics analysis detected the Csep1 protein. A csep1 gene hot spot in the chromosome of different C. concisus strains was found. The pICON plasmid was only found in GS2 strains isolated from the two relapsed CD patients with small bowel complications. This study reports a C. concisus molecular marker (csep1-6bpi) that is associated with active CD.

PlasmidSeeker: identification of known plasmids from bacterial whole genome sequencing reads

Background

Plasmids play an important role in the dissemination of antibiotic resistance, making their detection an important task. Using whole genome sequencing (WGS), it is possible to capture both bacterial and plasmid sequence data, but short read lengths make plasmid detection a complex problem.

Results

We developed a tool named PlasmidSeeker that enables the detection of plasmids from bacterial WGS data without read assembly. The PlasmidSeeker algorithm is based on k-mers and uses k-mer abundance to distinguish between plasmid and bacterial sequences. We tested the performance of PlasmidSeeker on a set of simulated and real bacterial WGS samples, resulting in 100% sensitivity and 99.98% specificity.

Conclusion

PlasmidSeeker enables quick detection of known plasmids and complements existing tools that assemble plasmids de novo. The PlasmidSeeker source code is stored on GitHub: https://github.com/bioinfo-ut/PlasmidSeeker.

The global distribution and spread of the mobilized colistin resistance gene mcr-1

Colistin represents one of the few available drugs for treating infections caused by carbapenem-resistant Enterobacteriaceae. As such, the recent plasmid-mediated spread of the colistin resistance gene mcr-1 poses a significant public health threat, requiring global monitoring and surveillance. Here, we characterize the global distribution of mcr-1 using a data set of 457 mcr-1-positive sequenced isolates. We find mcr-1 in various plasmid types but identify an immediate background common to all mcr-1 sequences. Our analyses establish that all mcr-1 elements in circulation descend from the same initial mobilization of mcr-1 by an ISApl1 transposon in the mid 2000s (2002–2008; 95% highest posterior density), followed by a marked demographic expansion, which led to its current global distribution. Our results provide the first systematic phylogenetic analysis of the origin and spread of mcr-1, and emphasize the importance of understanding the movement of antibiotic resistance genes across multiple levels of genomic organization.

The recent plasmid-mediated spread of the mobilized colistin resistance gene mcr-1 poses a significant public health threat, requiring worldwide monitoring and surveillance. Here, Wang et al. compile and analyze a data set of 457 mcr-1-positive sequenced isolates to investigate the origin and global distribution of mcr-1.

Evolution of a Dominant Natural Isolate of Escherichia coli in the Human Gut over the Course of a Year Suggests a Neutral Evolution with Reduced Effective Population Size

In vitro and in vivo evolution experiments on Escherichia coli revealed several principles of bacterial adaptation. However, few data are available in the literature describing the behavior of E. coli in its natural environment. We attempted here to study the evolution in the human gut of a commensal dominant E. coli clone, ED1a belonging to the B2 phylogroup, through a longitudinal genomic study. We sequenced 24 isolates sampled at three different time points within a healthy individual over almost a year. We computed a mutation rate of 6.90 × 10^-7 mutations per base per year of the chromosome for E. coli ED1a in healthy human gut. We observed very limited genomic diversity and could not detect any evidence of selection, in contrast to what is observed in experimental evolution over a similar length of time. We therefore suggest that ED1a, being well adapted to the healthy human gut, evolves mostly neutrally with a low effective population size (N_e of ≈500 to 1,700).IMPORTANCE In this study, we follow the genomic fate of a dominant clone of Escherichia coli in the human gut of a healthy individual over about a year. We could compute a low annual mutation rate that supports low diversity, and we could not retrieve any clear signature of selection. These observations support a neutral evolution of E. coli in the human gut, compatible with a very limited effective population size that deviates drastically with the observations made previously in experimental evolution.

PlasmidTron: assembling the cause of phenotypes and genotypes from NGS data

Increasingly rich metadata are now being linked to samples that have been whole-genome sequenced. However, much of this information is ignored. This is because linking this metadata to genes, or regions of the genome, usually relies on knowing the gene sequence(s) responsible for the particular trait being measured and looking for its presence or absence in that genome. Examples of this would be the spread of antimicrobial resistance genes carried on mobile genetic elements (MGEs). However, although it is possible to routinely identify the resistance gene, identifying the unknown MGE upon which it is carried can be much more difficult if the starting point is short-read whole-genome sequence data. The reason for this is that MGEs are often full of repeats and so assemble poorly, leading to fragmented consensus sequences. Since mobile DNA, which can carry many clinically and ecologically important genes, has a different evolutionary history from the host, its distribution across the host population will, by definition, be independent of the host phylogeny. It is possible to use this phenomenon in a genome-wide association study to identify both the genes associated with the specific trait and also the DNA linked to that gene, for example the flanking sequence of the plasmid vector on which it is encoded, which follows the same patterns of distribution as the marker gene/sequence itself. We present PlasmidTron, which utilizes the phenotypic data normally available in bacterial population studies, such as antibiograms, virulence factors, or geographical information, to identify traits that are likely to be present on DNA that can randomly reassort across defined bacterial populations. It is also possible to use this methodology to associate unknown genes/sequences (e.g. plasmid backbones) with a specific molecular signature or marker (e.g. resistance gene presence or absence) using PlasmidTron. PlasmidTron uses a k-mer-based approach to identify reads associated with a phylogenetically unlinked phenotype. These reads are then assembled de novo to produce contigs in a fast and scalable-to-large manner. PlasmidTron is written in Python 3 and is available under the open source licence GNU GPL3 from https://github.com/sanger-pathogens/plasmidtron.

Transfer of a bla CTX-M-1-carrying plasmid between different Escherichia coli strains within the human gut explored by whole genome sequencing analyses

Horizontal transfer of antibiotic resistance determinants contributes to dissemination of antibiotic resistance. Such transfer of resistance genes within the human gut has been documented in some in vivo studies. The present study investigated seven bla _CTX-M-1-carrying Escherichia coli isolates from three consecutive faecal samples collected from one cystic fibrosis patient in a nine-months period, by analysing whole genome sequencing data. The analyses showed that the seven E. coli isolates represented three genetically diverse strains. All isolates contained bla _CTX-M-1-carrying Incl1 plasmids that shared a common 101 kb backbone differing by only four SNPs. The plasmids harboured by the three different E. coli strains varied within limited regions suggestive of recombination events, according to the phylogenetic topology of the genomes of the isolates harbouring them. The findings strongly suggest that horizontal transfer of a bla _CTX-M-1-carrying plasmid had occurred within the patient´s gut. The study illustrates the within-host diversity of faecally carried resistant E. coli isolates and highlights the value of collecting multiple bacterial colonies from longitudinally collected samples to assess faecal carriage of resistant enterobacteria. The clustering of the plasmids with the corresponding E. coli strains carrying them indicates that the plasmids appear to have adapted to their respective E. coli hosts.

Within-host evolution of Enterococcus faecium during longitudinal carriage and transition to bloodstream infection in immunocompromised patients

BACKGROUND

Enterococcus faecium is a leading cause of hospital-acquired infection, particularly in the immunocompromised. Here, we use whole genome sequencing of E. faecium to study within-host evolution and the transition from gut carriage to invasive disease.

METHODS

We isolated and sequenced 180 E. faecium from four immunocompromised patients who developed bloodstream infection during longitudinal surveillance of E. faecium in stool and their immediate environment.

RESULTS

A phylogenetic tree based on single nucleotide polymorphisms (SNPs) in the core genome of the 180 isolates demonstrated several distinct clones. This was highly concordant with the population structure inferred by Bayesian methods, which contained four main BAPS (Bayesian Analysis of Population Structure) groups. The majority of isolates from each patient resided in a single group, but all four patients also carried minority populations in stool from multiple phylogenetic groups. Bloodstream isolates from each case belonged to a single BAPS group, which differed in all four patients. Analysis of 87 isolates (56 from blood) belonging to a single BAPS group that were cultured from the same patient over 54 days identified 30 SNPs in the core genome (nine intergenic, 13 non-synonymous, eight synonymous), and 250 accessory genes that were variably present. Comparison of these genetic variants in blood isolates versus those from stool or environment did not identify any variants associated with bloodstream infection. The substitution rate for these isolates was estimated to be 128 (95% confidence interval 79.82 181.77) mutations per genome per year, more than ten times higher than previous estimates for E. faecium. Within-patient variation in vancomycin resistance associated with vanA was common and could be explained by plasmid loss, or less often by transposon loss.

CONCLUSIONS

These findings demonstrate the diversity of E. faecium carriage by individual patients and significant within-host diversity of E. faecium, but do not provide evidence for adaptive genetic variation associated with invasion.

Occurrence and Genomic Characterization of ESBL-Producing, MCR-1-Harboring Escherichia coli in Farming Soil

The emergence and spread of the mobile colistin resistance gene (mcr-1) has become a major global public health concern. So far, this gene has been widely detected in food animals, pets, food, and humans. However, there is little information on the contamination of mcr-1-containing bacteria in farming soils. In August 2016, a survey of ESBL-producing Escherichia coli isolated from farming soils was conducted in Shandong Province, China. We observed colistin resistance in 12 of 53 (22.6%) ESBL-producing Enterobacteriaceae isolates from farming soil. Six mcr-1-positive E. coli strains originating from a livestock-intensive area were found. The isolates belonged to four different STs (ST2060, ST3014, ST6756, and ST1560) and harbored extensive additional resistance genes. An E. coli with bla_NDM-1 was also detected in a soil sample from the same area. Comparative whole genome sequencing and S1-PFGE analysis indicated that mcr-1 was chromosomally encoded in four isolates and located on IncHI2 plasmids in two isolates. To our knowledge, we report the first isolation of mcr-1 in ESBL-producing E. coli from farming soils. This work highlights the importance of active surveillance of colistin-resistant organisms in soil. Moreover, investigations addressing the influence of animal manure application on the transmission of mcr-1-producing bacteria are also warranted.

Recycler: an algorithm for detecting plasmids from de novo assembly graphs

Motivation

Plasmids and other mobile elements are central contributors to microbial evolution and genome innovation. Recently, they have been found to have important roles in antibiotic resistance and in affecting production of metabolites used in industrial and agricultural applications. However, their characterization through deep sequencing remains challenging, in spite of rapid drops in cost and throughput increases for sequencing. Here, we attempt to ameliorate this situation by introducing a new circular element assembly algorithm, leveraging assembly graphs provided by a conventional de novo assembler and alignments of paired-end reads to assemble cyclic sequences likely to be plasmids, phages and other circular elements.

Results

We introduce Recycler, the first tool that can extract complete circular contigs from sequence data of isolate microbial genomes, plasmidome and metagenome sequence data. We show that Recycler greatly increases the number of true plasmids recovered relative to other approaches while remaining highly accurate. We demonstrate this trend via simulations of plasmidomes, comparisons of predictions with reference data for isolate samples, and assessments of annotation accuracy on metagenome data. In addition, we provide validation by DNA amplification of 77 plasmids predicted by Recycler from the different sequenced samples in which Recycler showed mean accuracy of 89% across all data types—isolate, microbiome and plasmidome.

Availability and Implementation

Recycler is available at http://github.com/Shamir-Lab/Recycler

Supplementary information

Supplementary data are available at Bioinformatics online.

Genomic characterisation of clinical and environmental Pseudomonas putida group strains and determination of their role in the transfer of antimicrobial resistance genes to Pseudomonas aeruginosa

Background

Pseudomonas putida is a Gram-negative, non-fermenting bacterium frequently encountered in various environmental niches. P. putida rarely causes disease in humans, though serious infections and outbreaks have been reported from time to time. Some have suggested that P. putida functions as an exchange platform for antibiotic resistance genes (ARG), and thus represents a serious concern in the spread of ARGs to more pathogenic organisms within a hospital. Though poorly understood, the frequency of ARG exchange between P. putida and the more virulent Pseudomonas aeruginosa and its clinical relevance are particularly important for designing efficient infection control strategies, such as deciding whether high-risk patients colonized with a multidrug resistant but typically low pathogenic P. putida strain should be contact isolated or not.

Results

In this study, 21,373 screening samples (stool, rectal and throat swab) were examined to determine the presence of P. putida in a high-risk group of haemato-oncology patients during a 28-month period. A total of 89 P. putida group strains were isolated from 85 patients, with 41 of 89 (46.1%) strains harbouring the metallo-beta-lactamase gene bla_VIM. These 41 clinical isolates, plus 18 bla_VIM positive environmental P. putida isolates, and 17 bla_VIM positive P. aeruginosa isolates, were characterized by whole genome sequencing (WGS).

We constructed a maximum-likelihood tree to separate the 59 bla_VIM positive P. putida group strains into eight distinct phylogenetic clusters. Bla_VIM-1 was present in 6 clusters while bla_VIM-2 was detected in 4 clusters. Five P. putida group strains contained both, bla_VIM-1 and bla_VIM-2 genes.

In contrast, all P. aeruginosa strains belonged to a single genetic cluster and contained the same ARGs. Apart from bla_VIM-2 and sul genes, no other ARGs were shared between P. aeruginosa and P. putida. Furthermore, the bla_VIM-2 gene in P. aeruginosa was predicted to be only chromosomally located.

Conclusion

These data provide evidence that no exchange of comprehensive ARG harbouring mobile genetic elements had occurred between P. aeruginosa and P. putida group strains during the study period, thus eliminating the need to implement enhanced infection control measures for high-risk patients colonized with a bla_VIM positiv P. putida group strains in our clinical setting.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4216-2) contains supplementary material, which is available to authorized users.

Complete Genome Sequence of Exiguobacterium sp. Strain N4-1P, a Psychrophilic Bioemulsifier Producer Isolated from a Cold Marine Environment in North Atlantic Canada

Here, we present the complete genome sequence of Exiguobacterium sp. strain N4-1P, a psychrophilic bacterium that produces bioemulsifier, isolated for the first time from petroleum hydrocarbon-contaminated sediment samples from shoreline Newfoundland, Canada. Many strains of the genus Exiguobacterium are extremophiles and have properties of biotechnological interest.

Evolutionary Thrift: Mycobacteria Repurpose Plasmid Diversity during Adaptation of Type VII Secretion Systems

Mycobacteria have a distinct secretion system, termed type VII (T7SS), which is encoded by paralogous chromosomal loci (ESX) and associated with pathogenesis, conjugation, and metal homeostasis. Evolution of paralogous gene families is of interest because duplication is an important mechanism by which novel genes evolve, but there are potential conflicts between adaptive forces that stabilize duplications and those that enable evolution of new functions. Our objective was to delineate the adaptive forces underlying diversification of T7SS. Plasmid-borne ESX were described recently, and we found evidence that the initial duplication and divergence of ESX systems occurred on plasmids and was driven by selection for advantageous mutations. Plasmid conjugation has been linked to T7SS and type IV secretion systems (T4SS) in mycobacteria, and we discovered that T7SS and T4SS genes evolved in concert on the plasmids. We hypothesize that differentiation of plasmid ESX helps to prevent conjugation among cells harboring incompatible plasmids. Plasmid ESX appear to have been repurposed following migration to the chromosome, and there is evidence of positive selection driving further differentiation of chromosomal ESX. We hypothesize that ESX loci were initially stabilized on the chromosome by mediating their own transfer. These results emphasize the diverse adaptive paths underlying evolution of novelty, which in this case involved plasmid duplications, selection for advantageous mutations in the mobile and core genomes, migration of the loci between plasmids and chromosomes, and lateral transfer among chromosomes. We discuss further implications for the choice of model organism to study ESX functions in Mycobacterium tuberculosis.

In vivo selection of a multidrug-resistant Aeromonas salmonicida during medicinal leech therapy

We report the selection in a 15-year-old boy of a multidrug-resistant, extended-spectrum β-lactamase (ESBL)-producing Aeromonas salmonicida after medicinal leech therapy that required an antibiotic prophylaxis based on piperacillin/tazobactam and cotrimoxazole. Whole genome sequencing of the strain indeed revealed 13 antibiotic resistance genes, including the ESBL CTX-M-3 and the unusual β-lactamase SCO-1.

Genome and Plasmid Sequences of Escherichia coli KV7, an Extended-Spectrum β-Lactamase Isolate Derived from Feces of a Healthy Pig

We present single-contig assemblies for Escherichia coli strain KV7 (serotype O27, phylogenetic group D) and its six plasmids, isolated from a healthy pig, as determined by PacBio RS II and Illumina MiSeq sequencing. The chromosome of 4,997,475 bp and G+C content of 50.75% harbored 4,540 protein-encoding genes.

On the (im)possibility of reconstructing plasmids from whole-genome short-read sequencing data

To benchmark algorithms for automated plasmid sequence reconstruction from short-read sequencing data, we selected 42 publicly available complete bacterial genome sequences spanning 12 genera, containing 148 plasmids. We predicted plasmids from short-read data with four programs (PlasmidSPAdes, Recycler, cBar and PlasmidFinder) and compared the outcome to the reference sequences. PlasmidSPAdes reconstructs plasmids based on coverage differences in the assembly graph. It reconstructed most of the reference plasmids (recall=0.82), but approximately a quarter of the predicted plasmid contigs were false positives (precision=0.75). PlasmidSPAdes merged 84 % of the predictions from genomes with multiple plasmids into a single bin. Recycler searches the assembly graph for sub-graphs corresponding to circular sequences and correctly predicted small plasmids, but failed with long plasmids (recall=0.12, precision=0.30). cBar, which applies pentamer frequency analysis to detect plasmid-derived contigs, showed a recall and precision of 0.76 and 0.62, respectively. However, cBar categorizes contigs as plasmid-derived and does not bin the different plasmids. PlasmidFinder, which searches for replicons, had the highest precision (1.0), but was restricted by the contents of its database and the contig length obtained from de novo assembly (recall=0.36). PlasmidSPAdes and Recycler detected putative small plasmids (<10 kbp), which were also predicted as plasmids by cBar, but were absent in the original assembly. This study shows that it is possible to automatically predict small plasmids. Prediction of large plasmids (>50 kbp) containing repeated sequences remains challenging and limits the high-throughput analysis of plasmids from short-read whole-genome sequencing data.

De Novo Assembly of Complete Chloroplast Genomes from Non-model Species Based on a K-mer Frequency-Based Selection of Chloroplast Reads from Total DNA Sequences

Whole Genome Shotgun (WGS) sequences of plant species often contain an abundance of reads that are derived from the chloroplast genome. Up to now these reads have generally been identified and assembled into chloroplast genomes based on homology to chloroplasts from related species. This re-sequencing approach may select against structural differences between the genomes especially in non-model species for which no close relatives have been sequenced before. The alternative approach is to de novo assemble the chloroplast genome from total genomic DNA sequences. In this study, we used k-mer frequency tables to identify and extract the chloroplast reads from the WGS reads and assemble these using a highly integrated and automated custom pipeline. Our strategy includes steps aimed at optimizing assemblies and filling gaps which are left due to coverage variation in the WGS dataset. We have successfully de novo assembled three complete chloroplast genomes from plant species with a range of nuclear genome sizes to demonstrate the universality of our approach: Solanum lycopersicum (0.9 Gb), Aegilops tauschii (4 Gb) and Paphiopedilum henryanum (25 Gb). We also highlight the need to optimize the choice of k and the amount of data used. This new and cost-effective method for de novo short read assembly will facilitate the study of complete chloroplast genomes with more accurate analyses and inferences, especially in non-model plant genomes.

Whole-Genome Sequences of 14 Strains of Bradyrhizobium canariense and 1 Strain of Bradyrhizobium japonicum Isolated from Lupinus spp. in Algeria

We report here the whole-genome sequences of 14 strains of Bradyrhizobium canariense, isolated from root nodules of Lupinus microanthus and Lupinus angustifolius, and 1 strain of Bradyrhizobium japonicum isolated from root nodules from Lupinus angustifolius in Algeria. These sequences add to the known diversity of this agronomically important genus.

A Case Study into Microbial Genome Assembly Gap Sequences and Finishing Strategies

This study characterized regions of DNA which remained unassembled by either PacBio and Illumina sequencing technologies for seven bacterial genomes. Two genomes were manually finished using bioinformatics and PCR/Sanger sequencing approaches and regions not assembled by automated software were analyzed. Gaps present within Illumina assemblies mostly correspond to repetitive DNA regions such as multiple rRNA operon sequences. PacBio gap sequences were evaluated for several properties such as GC content, read coverage, gap length, ability to form strong secondary structures, and corresponding annotations. Our hypothesis that strong secondary DNA structures blocked DNA polymerases and contributed to gap sequences was not accepted. PacBio assemblies had few limitations overall and gaps were explained as cumulative effect of lower than average sequence coverage and repetitive sequences at contig termini. An important aspect of the present study is the compilation of biological features that interfered with assembly and included active transposons, multiple plasmid sequences, phage DNA integration, and large sequence duplication. Our targeted genome finishing approach and systematic evaluation of the unassembled DNA will be useful for others looking to close, finish, and polish microbial genome sequences.

Biofilm Formation Potential of Heat-Resistant Escherichia coli Dairy Isolates and the Complete Genome of Multidrug-Resistant, Heat-Resistant Strain FAM21845

We tested the biofilm formation potential of 30 heat-resistant and 6 heat-sensitive Escherichia coli dairy isolates. Production of curli and cellulose, static biofilm formation on polystyrene (PS) and stainless steel surfaces, biofilm formation under dynamic conditions (Bioflux), and initial adhesion rates (IAR) were evaluated. Biofilm formation varied greatly between strains, media, and assays. Our results highlight the importance of the experimental setup in determining biofilm formation under conditions of interest, as correlation between different assays was often not a given. The heat-resistant, multidrug-resistant (MDR) strain FAM21845 showed the strongest biofilm formation on PS and the highest IAR and was the only strain that formed significant biofilms on stainless steel under conditions relevant to the dairy industry, and it was therefore fully sequenced. Its chromosome is 4.9 Mb long, and it harbors a total of five plasmids (147.2, 54.2, 5.8, 2.5, and 1.9 kb). The strain carries a broad range of genes relevant to antimicrobial resistance and biofilm formation, including some on its two large conjugative plasmids, as demonstrated in plate mating assays.IMPORTANCE In biofilms, cells are embedded in an extracellular matrix that protects them from stresses, such as UV radiation, osmotic shock, desiccation, antibiotics, and predation. Biofilm formation is a major bacterial persistence factor of great concern in the clinic and the food industry. Many tested strains formed strong biofilms, and especially strains such as the heat-resistant, MDR strain FAM21845 may pose a serious issue for food production. Strong biofilm formation combined with diverse resistances (some encoded on conjugative plasmids) may allow for increased persistence, coselection, and possible transfer of these resistance factors. Horizontal gene transfer may conceivably occur in the food production setting or the gastrointestinal tract after consumption.

Resolving plasmid structures in Enterobacteriaceae using the MinION nanopore sequencer: assessment of MinION and MinION/Illumina hybrid data assembly approaches

This study aimed to assess the feasibility of using the Oxford Nanopore Technologies (ONT) MinION long-read sequencer in reconstructing fully closed plasmid sequences from eight Enterobacteriaceae isolates of six different species with plasmid populations of varying complexity. Species represented were Escherichia coli, Klebsiella pneumoniae, Citrobacter freundii, Enterobacter cloacae, Serratia marcescens and Klebsiella oxytoca, with plasmid populations ranging from 1–11 plasmids with sizes of 2–330 kb. Isolates were sequenced using Illumina (short-read) and ONT’s MinION (long-read) platforms, and compared with fully resolved PacBio (long-read) sequence assemblies for the same isolates. We compared the performance of different assembly approaches including SPAdes, plasmidSPAdes, hybridSPAdes, Canu, Canu+Pilon (canuPilon) and npScarf in recovering the plasmid structures of these isolates by comparing with the gold-standard PacBio reference sequences. Overall, canuPilon provided consistently good quality assemblies both in terms of assembly statistics (N50, number of contigs) and assembly accuracy [presence of single nucleotide polymorphisms (SNPs)/indels with respect to the reference sequence]. For plasmid reconstruction, Canu recovered 70 % of the plasmids in complete contigs, and combining three assembly approaches (Canu or canuPilon, hybridSPAdes and plasmidSPAdes) resulted in a total 78 % recovery rate for all the plasmids. The analysis demonstrated the potential of using MinION sequencing technology to resolve important plasmid structures in Enterobacteriaceae species independent of and in conjunction with Illumina sequencing data. A consensus assembly derived from several assembly approaches could present significant benefit in accurately resolving the greatest number of plasmid structures.

Clostridium chauvoei, an Evolutionary Dead-End Pathogen

Full genome sequences of 20 strains of Clostridium chauvoei, the etiological agent of blackleg of cattle and sheep, isolated from four different continents over a period of 64 years (1951–2015) were determined and analyzed. The study reveals that the genome of the species C. chauvoei is highly homogeneous compared to the closely related species C. perfringens, a widespread pathogen that affects human and many animal species. Analysis of the CRISPR locus is sufficient to differentiate most C. chauvoei strains and is the most heterogenous region in the genome, containing in total 187 different spacer elements that are distributed as 30 – 77 copies in the various strains. Some genetic differences are found in the 3 allelic variants of fliC1, fliC2 and fliC3 genes that encode structural flagellin proteins, and certain strains do only contain one or two alleles. However, the major virulence genes including the highly toxic C.chauvoei toxin A, the sialidase and the two hyaluronidases are fully conserved as are the metabolic and structural genes of C. chauvoei. These data indicate that C. chauvoei is a strict ruminant-associated pathogen that has reached a dead end in its evolution.

Genome sequence of Shigella flexneri strain SP1, a diarrheal isolate that encodes an extended-spectrum β-lactamase (ESBL)

BACKGROUND

Shigellosis is the most common cause of gastrointestinal infections in developing countries. In China, the species most frequently responsible for shigellosis is Shigella flexneri. S. flexneri remains largely unexplored from a genomic standpoint and is still described using a vocabulary based on biochemical and serological properties. Moreover, increasing numbers of ESBL-producing Shigella strains have been isolated from clinical samples. Despite this, only a few cases of ESBL-producing Shigella have been described in China. Therefore, a better understanding of ESBL-producing Shigella from a genomic standpoint is required. In this study, a S. flexneri type 1a isolate SP1 harboring bla_CTX-M-14, which was recovered from the patient with diarrhea, was subjected to whole genome sequencing.

RESULTS

The draft genome assembly of S. flexneri strain SP1 consisted of 4,592,345 bp with a G+C content of 50.46%. RAST analysis revealed the genome contained 4798 coding sequences (CDSs) and 100 RNA-encoding genes. We detected one incomplete prophage and six candidate CRISPR loci in the genome. In vitro antimicrobial susceptibility testing demonstrated that strain SP1 is resistant to ampicillin, amoxicillin/clavulanic acid, cefazolin, ceftriaxone and trimethoprim. In silico analysis detected genes mediating resistance to aminoglycosides, β-lactams, phenicol, tetracycline, sulphonamides, and trimethoprim. The bla _CTX-M-14 gene was located on an IncFII2 plasmid. A series of virulence factors were identified in the genome.

CONCLUSIONS

In this study, we report the whole genome sequence of a bla_CTX-M-14-encoding S. flexneri strain SP1. Dozens of resistance determinants were detected in the genome and may be responsible for the multidrug-resistance of this strain, although further confirmation studies are warranted. Numerous virulence factors identified in the strain suggest that isolate SP1 is potential pathogenic. The availability of the genome sequence and comparative analysis with other S. flexneri strains provides the basis to further address the evolution of drug resistance mechanisms and pathogenicity in S. flexneri.

Metagenome-Assembled Genome Sequence of Rhodopseudomonas palustris Strain ELI 1980, Commercialized as a Biostimulant

We report here the draft genome sequence of strain ELI 1980 of Rhodopseudomonas palustris, commercialized as a biostimulant for agriculture. The genome was reconstructed from the metagenome of a commercial product containing this strain as its major component.

Whole genome sequences of three Clade 3 Clostridium difficile strains carrying binary toxin genes in China

Clostridium difficile consists of six clades but studies on Clade 3 are limited. Here, we report genome sequences of three Clade 3 C. difficile strains carrying genes encoding toxin A and B and the binary toxin. Isolates 103 and 133 (both of ST5) and isolate 106 (ST285) were recovered from three ICU patients. Whole genome sequencing using HiSeq 2500 revealed 4.1-Mb genomes with 28–29% GC content. There were ≥1,104 SNP between the isolates, suggesting they were not of a single clone. The toxin A and B gene-carrying pathogenicity locus (PaLoc) of the three isolates were identical and had the insertion of the transposon Tn6218. The genetic components of PaLoc among Clade 3 strains were the same with only a few nucleotide mutations and deletions/insertions, suggesting that the Tn6218 insertion might have occurred before the divergence within Clade 3. The binary toxin-genes carrying CDT locus (CdtLoc) of the three isolates were identical and were highly similar to those of other Clade 3 strains, but were more divergent from those of other clades. In conclusion, Clade 3 has an unusual clade-specific PaLoc characteristic of a Tn6218 insertion which appears to be the main feature to distinguish Clade 3 from other C. difficile.

Plasmid Classification in an Era of Whole-Genome Sequencing: Application in Studies of Antibiotic Resistance Epidemiology

Plasmids are extra-chromosomal genetic elements ubiquitous in bacteria, and commonly transmissible between host cells. Their genomes include variable repertoires of 'accessory genes,' such as antibiotic resistance genes, as well as 'backbone' loci which are largely conserved within plasmid families, and often involved in key plasmid-specific functions (e.g., replication, stable inheritance, mobility). Classifying plasmids into different types according to their phylogenetic relatedness provides insight into the epidemiology of plasmid-mediated antibiotic resistance. Current typing schemes exploit backbone loci associated with replication (replicon typing), or plasmid mobility (MOB typing). Conventional PCR-based methods for plasmid typing remain widely used. With the emergence of whole-genome sequencing (WGS), large datasets can be analyzed using in silico plasmid typing methods. However, short reads from popular high-throughput sequencers can be challenging to assemble, so complete plasmid sequences may not be accurately reconstructed. Therefore, localizing resistance genes to specific plasmids may be difficult, limiting epidemiological insight. Long-read sequencing will become increasingly popular as costs decline, especially when resolving accurate plasmid structures is the primary goal. This review discusses the application of plasmid classification in WGS-based studies of antibiotic resistance epidemiology; novel in silico plasmid analysis tools are highlighted. Due to the diverse and plastic nature of plasmid genomes, current typing schemes do not classify all plasmids, and identifying conserved, phylogenetically concordant genes for subtyping and phylogenetics is challenging. Analyzing plasmids as nodes in a network that represents gene-sharing relationships between plasmids provides a complementary way to assess plasmid diversity, and allows inferences about horizontal gene transfer to be made.