Members of the Mesorhizobium genus are soil bacteria that often form nitrogen-fixing symbioses with legumes. Most characterised Mesorhizobium spp. genomes are ~8 Mb in size and harbour extensive pangenomes including large integrative and conjugative elements (ICEs) carrying genes required for symbiosis (ICESyms). Here, we document and compare the conjugative mobilome of 41 complete Mesorhizobium genomes. We delineated 56 ICEs and 24 integrative and mobilizable elements (IMEs) collectively occupying 16 distinct integration sites, along with 24 plasmids. We also demonstrated horizontal transfer of the largest (853,775 bp) documented ICE, the tripartite ICEMspSymAA22. The conjugation systems of all identified ICEs and several plasmids were related to those of the paradigm ICESym ICEMlSymR7A, with each carrying conserved genes for conjugative pilus formation (trb), excision (rdfS), DNA transfer (rlxS) and regulation (fseA). ICESyms have likely evolved from a common ancestor, despite occupying a variety of distinct integration sites and specifying symbiosis with diverse legumes. We found extensive evidence for recombination between ICEs and particularly ICESyms, which all uniquely lack the conjugation entry-exclusion factor gene trbK. Frequent duplication, replacement and pseudogenization of genes for quorum-sensing-mediated activation and antiactivation of ICE transfer suggests ICE transfer regulation is constantly evolving. Pangenome-wide association analysis of the ICE identified genes potentially involved in symbiosis, rhizosphere colonisation and/or adaptation to distinct legume hosts. In summary, the Mesorhizobium genus has accumulated a large and dynamic pangenome that evolves through ongoing horizontal gene transfer of large conjugative elements related to ICEMlSymR7A.

Many public health laboratories across the world have implemented whole-genome sequencing (WGS) for the surveillance and outbreak detection of foodborne pathogens. PulseNet-affiliated laboratories have determined that most single-strain foodborne outbreaks are contained within 0–10 multi-locus sequence typing (MLST)-based allele differences and/or core genome single-nucleotide variants (SNVs). In addition to being a food- and travel-associated outbreak pathogen, most Shigella spp. cases occur through continuous person-to-person transmission, predominantly involving men who have sex with men (MSM), leading to long-term and recurrent outbreaks. Continuous transmission patterns coupled to genetic evolution under antibiotic treatment pressure require an assessment of existing WGS-based subtyping methods and interpretation criteria for cluster inclusion/exclusion. An evaluation of 4 WGS-based subtyping methods [SNVPhyl, coreMLST, core genome MLST (cgMLST) and whole-genome MLST (wgMLST)] was performed on 9 foodborne-, travel- and MSM-related retrospective outbreaks from a collection of 91 Shigella flexneri and 232  Shigella sonnei isolates to determine the methods’ epidemiological concordance, discriminatory power, robustness and ability to generate stable interpretation criteria. The discriminatory powers were ranked as follows: coreMLST

Macrolides are broad-spectrum antibiotics used to treat a range of infections. Resistance to macrolides is often conferred by mobile resistance genes encoding Erm methyltransferases or Mph phosphotransferases. New erm and mph genes keep being discovered in clinical settings but their origins remain unknown, as is the type of macrolide resistance genes that will appear in the future. In this study, we used optimized hidden Markov models to characterize the macrolide resistome. Over 16 terabases of genomic and metagenomic data, representing a large taxonomic diversity (11 030 species) and diverse environments (1944 metagenomic samples), were searched for the presence of erm and mph genes. From this data, we predicted 28 340 macrolide resistance genes encoding 2892 unique protein sequences, which were clustered into 663 gene families (<70 % amino acid identity), of which 619 (94 %) were previously uncharacterized. This included six new resistance gene families, which were located on mobile genetic elements in pathogens. The function of ten predicted new resistance genes were experimentally validated in Escherichia coli using a growth assay. Among the ten tested genes, seven conferred increased resistance to erythromycin, with five genes additionally conferring increased resistance to azithromycin, showing that our models can be used to predict new functional resistance genes. Our analysis also showed that macrolide resistance genes have diverse origins and have transferred horizontally over large phylogenetic distances into human pathogens. This study expands the known macrolide resistome more than ten-fold, provides insights into its evolution, and demonstrates how computational screening can identify new resistance genes before they become a significant clinical problem.

The roseobacter group of marine bacteria is characterized by a mosaic distribution of ecologically important phenotypes. These are often encoded on mobile extrachromosomal replicons. So far, conjugation had only been experimentally proven between the two model organisms Phaeobacter inhibens and Dinoroseobacter shibae. Here, we show that two large natural RepABC-type plasmids from D. shibae can be transferred into representatives of all known major Rhodobacterales lineages. Complete genome sequencing of the newly established Phaeobacter inhibens transconjugants confirmed their genomic integrity. The conjugated plasmids were stably maintained as single copy number replicons in the genuine as well as the new host. Co-cultivation of Phaeobacter inhibens and the transconjugants with the dinoflagellate Prorocentrum minimum demonstrated that Phaeobacter inhibens is a probiotic strain that improves the yield and stability of the dinoflagellate culture. The transconjugant carrying the 191 kb plasmid, but not the 126 kb sister plasmid, killed the dinoflagellate in co-culture.

Microbial genomes are highly adaptable, with mobile genetic elements (MGEs) such as integrative conjugative elements (ICEs) mediating the dissemination of new genetic information throughout bacterial populations. This is countered by defence mechanisms such as CRISPR-Cas systems, which limit invading MGEs by sequence-specific targeting. Here we report the distribution of the pVir, pTet and PCC42 plasmids and a new 70–129 kb ICE (CampyICE1) in the foodborne bacterial pathogens Campylobacter jejuni and Campylobacter coli . CampyICE1 contains a degenerated Type II-C CRISPR system consisting of a sole Cas9 protein, which is distinct from the previously described Cas9 proteins from C. jejuni and C. coli . CampyICE1 is conserved in structure and gene order, containing blocks of genes predicted to be involved in recombination, regulation and conjugation. CampyICE1 was detected in 134/5829 (2.3 %) C . jejuni genomes and 92/1347 (6.8 %) C . coli genomes. Similar ICEs were detected in a number of non-jejuni/coli Campylobacter species, although these lacked a CRISPR-Cas system. CampyICE1 carries three separate short CRISPR spacer arrays containing a combination of 108 unique spacers and 16 spacer-variant families. A total of 69 spacers and 10 spacer-variant families (63.7 %) were predicted to target Campylobacter plasmids. The presence of a functional CampyICE1 Cas9 protein and matching anti-plasmid spacers was associated with the absence of the pVir, pTet and pCC42 plasmids (188/214 genomes, 87.9 %), suggesting that the CampyICE1-encoded CRISPR-Cas has contributed to the exclusion of competing plasmids. In conclusion, the characteristics of the CRISPR-Cas9 system on CampyICE1 suggests a history of plasmid warfare in Campylobacter .

Streptomyces clavuligerus is an industrially important actinomycete whose genetic manipulation is limited by low transformation and conjugation efficiencies, low levels of recombination of introduced DNA, and difficulty in obtaining consistent sporulation. We describe the construction and application of versatile vectors for Cas9-mediated genome editing of this strain. To design spacer sequences with confidence, we derived a highly accurate genome assembly for an isolate of the type strain (ATCC 27064). This yielded a chromosome assembly (6.75 Mb) plus assemblies for pSCL4 (1795 kb) and pSCL2 (149 kb). The strain also carries pSCL1 (12 kb), but its small size resulted in only partial sequence coverage. The previously described pSCL3 (444 kb) is not present in this isolate. Using our Cas9 vectors, we cured pSCL4 with high efficiency by targeting the plasmid's parB gene. Five of the resulting pSCL4-cured isolates were characterized and all showed impaired sporulation. Shotgun genome sequencing of each of these derivatives revealed large deletions at the ends of the chromosomes in all of them, and for two clones sufficient sequence data was obtained to show that the chromosome had circularized. Taken together, these data indicate that pSCL4 is essential for the structural stability of the linear chromosome.

Plasmids are extrachromosomal genetic elements that replicate independently of the chromosome and play a vital role in the environmental adaptation of bacteria. Due to potential mobilization or conjugation capabilities, plasmids are important genetic vehicles for antimicrobial resistance genes and virulence factors with huge and increasing clinical implications. They are therefore subject to large genomic studies within the scientific community worldwide. As a result of rapidly improving next-generation sequencing methods, the quantity of sequenced bacterial genomes is constantly increasing, in turn raising the need for specialized tools to (i) extract plasmid sequences from draft assemblies, (ii) derive their origin and distribution, and (iii) further investigate their genetic repertoire. Recently, several bioinformatic methods and tools have emerged to tackle this issue; however, a combination of high sensitivity and specificity in plasmid sequence identification is rarely achieved in a taxon-independent manner. In addition, many software tools are not appropriate for large high-throughput analyses or cannot be included in existing software pipelines due to their technical design or software implementation. In this study, we investigated differences in the replicon distributions of protein-coding genes on a large scale as a new approach to distinguish plasmid-borne from chromosome-borne contigs. We defined and computed statistical discrimination thresholds for a new metric: the replicon distribution score (RDS), which achieved an accuracy of 96.6 %. The final performance was further improved by the combination of the RDS metric with heuristics exploiting several plasmid-specific higher-level contig characterizations. We implemented this workflow in a new high-throughput taxon-independent bioinformatics software tool called Platon for the recruitment and characterization of plasmid-borne contigs from short-read draft assemblies. Compared to PlasFlow, Platon achieved a higher accuracy (97.5 %) and more balanced predictions (F1=82.6 %) tested on a broad range of bacterial taxa and better or equal performance against the targeted tools PlasmidFinder and PlaScope on sequenced Escherichia coli isolates. Platon is available at: http://platon.computational.bio/.

Abstract The emergence of plasmid-mediated tigecycline-resistant strains is posing a serious threat to food safety and human health, which has attracted worldwide attention. The tigecycline resistance gene tet(X4) has been found in diverse sources, but the distribution of tet(X4) and its genetic background in the animal farming environment is not fully understood. Thirty-two tet(X)-positive Escherichia coli strains isolated from 159 samples collected from swine farms showed resistance to tigecycline. The tet(X)-positive strains were characterized by antimicrobial susceptibility testing, conjugation assay, PCR, Illumina and long-read Nanopore sequencing, and bioinformatics analysis. A total of 11 different sequence types (STs) were identified and most of them belonged to phylogroup A, except ST641. In total, 196 possible prophage sequences were identified and some of the prophage regions were found to carry resistance genes, including tet(X4). Furthermore, our results showed possible correlations between CRISPR spacer sequences and serotypes or STs. The co-existence of tigecycline-resistant tet(A) variants and tet(X4) complicates the evolution of vital resistance genes in farming environments. Further, four reorganization plasmids carrying tet(X4) were observed, and the formation mechanism mainly involved homologous recombination. These findings contribute significantly to a better understanding of the diversity and complexity of tet(X4)-bearing plasmids, an emerging novel public health concern.

Enterohaemorrhagic Escherichia coli (EHEC) is an important human pathogen worldwide. Although serotype O157 is currently the most dominant and important EHEC strain, serotypes O26, O111, O91, O103 and O121 are also recognized as serious pathogens that affect public health. EHEC outbreaks often occur in nurseries and elderly care facilities. In 2012, a nursery outbreak of EHEC O121 occurred during which the bacterium acquired a plasmid-borne extended-spectrum β-lactamase (ESBL) gene. ESBL-producing E. coli O86 was concurrently isolated from one of the EHEC patients. Therefore, we investigated the isolates by whole-genome sequence (WGS) analysis to elucidate the transmission dynamics of the EHEC strains and the ESBL plasmid. According to WGS-based phylogeny, all 17 EHEC O121 isolates were clonal, while E. coli O86 was genetically distant from the EHEC O121 isolates. The complete sequence of an ESBL plasmid encoding the CTX-M-55 β-lactamase was determined using S1-PFGE bands, and subsequent mapping of the WGS reads confirmed that the plasmid sequences from EHEC O121 and E. coli O86 were identical. Furthermore, conjugation experiments showed that the plasmid was capable of conjugative transfer. These results support the hypothesis that EHEC O121 acquired an ESBL-producing plasmid from E. coli O86 during the outbreak. This report demonstrates the importance of implementing preventive measures during EHEC outbreaks to control both secondary infection and the spread of antimicrobial resistance factors.

Carbapenem-hydrolysing enzymes belonging to the OXA-48-like group are encoded by blaOXA-48-like alleles and are abundant among Enterobacterales in the Netherlands. Therefore, the objective here was to investigate the characteristics, gene content and diversity of the blaOXA-48-like carrying plasmids and chromosomes of Escherichia coli and Klebsiella pneumoniae collected in the Dutch national surveillance from 2014 to 2019 in comparison with genome sequences from 29 countries. A combination of short-read genome sequencing with long-read sequencing enabled the reconstruction of 47 and 132 complete blaOXA-48-like plasmids for E. coli and K. pneumoniae, respectively. Seven distinct plasmid groups designated as pOXA-48-1 to pOXA-48-5, pOXA-181 and pOXA-232 were identified in the Netherlands which were similar to internationally reported plasmids obtained from countries from North and South America, Europe, Asia and Oceania. The seven plasmid groups varied in size, G+C content, presence of antibiotic resistance genes, replicon family and gene content. The pOXA-48-1 to pOXA-48-5 plasmids were variable, and the pOXA-181 and pOXA-232 plasmids were conserved. The pOXA-48-1, pOXA-48-2, pOXA-48-3 and pOXA-48-5 groups contained a putative conjugation system, but this was absent in the pOXA-48-4, pOXA-181 and pOXA-232 plasmid groups. pOXA-48 plasmids contained the PemI antitoxin, while the pOXA-181 and pOXA-232 plasmids did not. Furthermore, the pOXA-181 plasmids carried a virB2-virB3-virB9-virB10-virB11 type IV secretion system, while the pOXA-48 plasmids and pOXA-232 lacked this system. A group of non-related pOXA-48 plasmids from the Netherlands contained different resistance genes, non-IncL-type replicons or no replicons. Whole genome multilocus sequence typing revealed that the blaOXA-48-like plasmids were found in a wide variety of genetic backgrounds in contrast to chromosomally encoded blaOXA-48-like alleles. Chromosomally localized blaOXA-48 and blaOXA-244 alleles were located on genetic elements of variable sizes and comprised regions of pOXA-48 plasmids. The blaOXA-48-like genetic element was flanked by a direct repeat upstream of IS1R, and was found at multiple locations in the chromosomes of E. coli. Lastly, K. pneumoniae isolates carrying blaOXA-48 or blaOXA-232 were mostly resistant for meropenem, whereas E. coli blaOXA-48, blaOXA-181 and chromosomal blaOXA-48 or blaOXA-244 isolates were mostly sensitive. In conclusion, the overall blaOXA-48-like plasmid population in the Netherlands is conserved and similar to that reported for other countries, confirming global dissemination of blaOXA-48-like plasmids. Variations in size, presence of antibiotic resistance genes and gene content impacted pOXA-48, pOXA-181 and pOXA-232 plasmid architecture.

Bacterial wilt caused by the Ralstonia solanacearum species complex (RSSC) threatens the cultivation of important crops worldwide. We sequenced 30 RSSC phylotype I (R. pseudosolanacearum) strains isolated from pepper (Capsicum annuum) and tomato (Solanum lycopersicum) across the Republic of Korea. These isolates span the diversity of phylotype I, have extensive effector repertoires and are subject to frequent recombination. Recombination hotspots among South Korean phylotype I isolates include multiple predicted contact-dependent inhibition loci, suggesting that microbial competition plays a significant role in Ralstonia evolution. Rapid diversification of secreted effectors presents challenges for the development of disease-resistant plant varieties. We identified potential targets for disease resistance breeding by testing for allele-specific host recognition of T3Es present among South Korean phyloype I isolates. The integration of pathogen population genomics and molecular plant pathology contributes to the development of location-specific disease control and development of plant cultivars with durable resistance to relevant threats.

Among members of the Bacillales order, there are several species capable of forming a structure called an endospore. Endospores enable bacteria to survive under unfavourable growth conditions and germinate when environmental conditions are favourable again. Spore-coat proteins are found in a multilayered proteinaceous structure encasing the spore core and the cortex. They are involved in coat assembly, cortex synthesis and germination. Here, we aimed to determine the diversity and evolutionary processes that have influenced spore-coat genes in various spore-forming species of Bacillales using an in silico approach. For this, we used sequence similarity searching algorithms to determine the diversity of coat genes across 161 genomes of Bacillales. The results suggest that among Bacillales, there is a well-conserved core genome, composed mainly by morphogenetic coat proteins and spore-coat proteins involved in germination. However, some spore-coat proteins are taxa-specific. The best-conserved genes among different species may promote adaptation to changeable environmental conditions. Because most of the Bacillus species harbour complete or almost complete sets of spore-coat genes, we focused on this genus in greater depth. Phylogenetic reconstruction revealed eight monophyletic groups in the Bacillus genus, of which three are newly discovered. We estimated the selection pressures acting over spore-coat genes in these monophyletic groups using classical and modern approaches and detected horizontal gene transfer (HGT) events, which have been further confirmed by scanning the genomes to find traces of insertion sequences. Although most of the genes are under purifying selection, there are several cases with individual sites evolving under positive selection. Finally, the HGT results confirm that sporulation is an ancestral feature in Bacillus.

Infections with globally disseminated Shiga toxin-producing Escherichia coli (STEC) of the O113:H21 serotype can progress to severe clinical complications, such as hemolytic uremic syndrome (HUS). Two phylogeographically distinct clonal complexes have been established by multi locus sequence typing (MLST). Infections with ST-820 isolates circulating exclusively in Australia have caused severe human disease, such as HUS. Conversely, ST-223 isolates prevalent in the US and outside Australia seem to rarely cause severe human disease but are frequent contaminants. Following a genomic epidemiology approach, we wanted to gain insights into the underlying cause for this disparity. We examined the plasticity in the genome make-up and Shiga toxin production in a collection of 20 ST-820 and ST-223 strains isolated from produce, the bovine reservoir, and clinical cases. STEC are notorious for assembly into fragmented draft sequences when using short-read sequencing technologies due to the extensive and partly homologous phage complement. The application of long-read technology (LRT) sequencing yielded closed reference chromosomes and plasmids for two representative ST-820 and ST-223 strains. The established high-resolution framework, based on whole genome alignments, single nucleotide polymorphism (SNP)-typing and MLST, includes the chromosomes and plasmids of other publicly available O113:H21 sequences and allowed us to refine the phylogeographical boundaries of ST-820 and ST-223 complex isolates and to further identify a historic non-shigatoxigenic strain from Mexico as a quasi-intermediate. Plasmid comparison revealed strong correlations between the strains' featured pO113 plasmid genotypes and chromosomally inferred ST, which suggests coevolution of the chromosome and virulence plasmids. Our pathogenicity assessment revealed statistically significant differences in the Stx2a-production capabilities of ST-820 as compared to ST-223 strains under RecA-induced Stx phage mobilization, a condition that mimics Stx-phage induction. These observations suggest that ST-820 strains may confer an increased pathogenic potential in line with the strain-associated epidemiological metadata. Still, some of the tested ST-223 cultures sourced from contaminated produce or the bovine reservoir also produced Stx at levels comparable to those of ST-820 isolates, which calls for awareness and for continued surveillance of this lineage.