Exploring bacterial insertion sequences with ISfinder: objectives, uses, and future developments

Nature-based bioreactors: Tackling antibiotic resistance in urban wastewater treatment

The overuse and misuse of antibiotics have accelerated the selection of antibiotic-resistant bacteria, significantly impacting human, animal, and environmental health. As aquatic environments are vulnerable to antibiotic resistance, suitable management practices should be adopted to tackle this phenomenon. Here we show an effective, nature-based solution for reducing antibiotic resistance from actual wastewater. We utilize a bioreactor that relies on benthic (biofilms) and planktonic microbial communities to treat secondary effluent from a small urban wastewater treatment plant (<10,000 population equivalent). This treated effluent is eventually released into the local aquatic ecosystem. We observe high removal efficiency for genes that provide resistance to commonly used antibiotic families, as well as for mobile genetic elements that could potentially aid in their spread. Importantly, we notice a buildup of sulfonamide (sul1 and sul2) and tetracycline (tet(C), tet(G), and tetR) resistance genes specifically in biofilms. This advancement marks the initial step in considering this bioreactor as a nature-based, cost-effective tertiary treatment option for small UWWTPs facing antibiotic resistance challenges.

Exploring the Potential Probiotic Properties of Bifidobacterium breve DSM 32583-A Novel Strain Isolated from Human Milk

Human milk is the best nutrition for infants, providing optimal support for the developing immune system and gut microbiota. Hence, it has been used as source for probiotic strain isolation, including members of the genus Bifidobacterium, in an effort to provide beneficial effects to infants who cannot be exclusively breastfed. However, not all supplemented bifidobacteria can effectively colonise the infant gut, nor confer health benefits to the individual infant host; therefore, new isolates are needed to develop a range of dietary products for this specific age group. Here, we investigated the beneficial potential of Bifidobacterium breve DSM 32583 isolated from human milk. We show that in vitro B. breve DSM 32583 exhibited several characteristics considered fundamental for beneficial bacteria, including survival in conditions simulating those present in the digestive tract, adherence to human epithelial cell lines, and inhibition of growth of potentially pathogenic microorganisms. Its antibiotic resistance patterns were comparable to those of known beneficial bifidobacterial strains, and its genome did not contain plasmids nor virulence-associated genes. These results suggest that B. breve DSM 32583 is a potential probiotic candidate.

Genomic analysis of antibiotic resistance genes and mobile genetic elements in eight strains of nontyphoid Salmonella

Nontyphoidal Salmonella (NTS) is the main etiological agent of human nontyphoidal salmonellosis. The aim of this study was to analyze the epidemiological characteristics and horizontal transfer mechanisms of antimicrobial resistance (AMR) genes from eight strains of NTS detected in Zhenjiang City, Jiangsu Province, China. Fecal samples from outpatients with food-borne diarrhea were collected in 2022. The NTS isolates were identified, and their susceptibility was tested with the Vitek 2 Compact system. The genomes of the NTS isolates were sequenced with the Illumina NovaSeq platform and Oxford Nanopore Technologies platform. The AMR genes and mobile genetic elements (MGEs) were predicted with the relevant open access resources. Eight strains of NTS were isolated from 153 specimens, and Salmonella Typhimurium ST19 was the most prevalent serotype. The AMR gene with the highest detection rate was AAC(6')-Iaa (10.5%) followed by TEM-1 (7.9%), sul2 (6.6%), and tet(A) (5.3%). Eleven MGEs carrying 34 AMR genes were identified on the chromosomes of 3 of the 8 NTS, including 3 resistance islands, 6 composite transposons (Tns), and 2 integrons. Eighteen plasmids carrying 40 AMR genes were detected in the 8 NTS strains, including 6 mobilizable plasmids, 3 conjugative plasmids, and 9 nontransferable plasmids, 7 of which carried 10 composite Tns and 3 integrons. This study provided a theoretical basis, from a genetic perspective, for the prevention and control of NTS resistance in Zhenjiang City.

IMPORTANCE

Human nontyphoidal salmonellosis is one of the common causes of bacterial food-borne illnesses, with significant social and economic impacts, especially those caused by invasive multidrug-resistant nontyphoidal Salmonella, which entails high morbidity and mortality. Antimicrobial resistance is mainly mediated by drug resistance genes, and mobile genetic elements play key roles in the capture, accumulation, and dissemination of antimicrobial resistance genes. Therefore, it is necessary to study the epidemiological characteristics and horizontal transfer mechanisms of antimicrobial resistance genes of nontyphoidal Salmonella to prevent the spread of multidrug-resistant nontyphoidal Salmonella.

Capture of armA by a novel ISCR element, ISCR28

ISCR28 is a fully functional and active member of the IS91-like family of insertion sequences. ISCR28 is 1,708-bp long and contains a 1,293-bp long putative open reading frame that codes a transposase. Sixty ISCR28-containing sequences from GenBank generated 27 non-repeat genetic contexts, all of which represented naturally occurring biological events that had occurred in a wide range of gram-negative organisms. Insertion of ISCR28 into target DNA preferred the presence of a 5'-GXXT-3' sequence at its terIS (replication terminator) end. Loss of the first 4 bp of its oriIS (origin of replication) likely caused ISCR28 to be trapped in ISApl1-based transposons or similar structures. Loss of terIS and fusion with a mobile element upstream likely promoted co-transfer of ISCR28 and the downstream resistance genes. ArmA and its downstream intact ISCR28 can be excised from recombinant pKD46 plasmids forming circular intermediates, further elucidating its activity as a transposase.

Annotation and Comparative Genomics of Prokaryotic Transposable Elements

The data generated in nearly 30 years of bacterial genome sequencing has revealed the abundance of transposable elements (TE) and their importance in genome and transcript remodeling through the mediation of DNA insertions and deletions, structural rearrangements, and regulation of gene expression. Furthermore, what we have learned from studying transposition mechanisms and their regulation in bacterial TE is fundamental to our current understanding of TE in other organisms because much of what has been observed in bacteria is conserved in all domains of life. However, unlike eukaryotic TE, prokaryotic TE sequester and transmit important classes of genes that impact host fitness, such as resistance to antibiotics and heavy metals and virulence factors affecting animals and plants, among other acquired traits. This provides dynamism and plasticity to bacteria, which would otherwise be propagated clonally. The insertion sequences (IS), the simplest form of prokaryotic TE, are autonomous and compact mobile genetic elements. These can be organized into compound transposons, in which two similar IS can flank any DNA segment and render it transposable. Other more complex structures, called unit transposons, can be grouped into four major families (Tn3, Tn7, Tn402, Tn554) with specific genetic characteristics. This chapter will revisit the prominent structural features of these elements, focusing on a genomic annotation framework and comparative analysis. Relevant aspects of TE will also be presented, stressing their key position in genome impact and evolution, especially in the emergence of antimicrobial resistance and other adaptive traits.

Mobilizable plasmids drive the spread of antimicrobial resistance genes and virulence genes in Klebsiella pneumoniae

BACKGROUND

Klebsiella pneumoniae is a notorious clinical pathogen and frequently carries various plasmids, which are the main carriers of antimicrobial resistance and virulence genes. In comparison to self-transmissible conjugative plasmids, mobilizable plasmids have received much less attention due to their defects in conjugative elements. However, the contribution of mobilizable plasmids to the horizontal transfer of antimicrobial resistance genes and virulence genes of K. pneumoniae remains unclear. In this study, the transfer, stability, and cargo genes of the mobilizable plasmids of K. pneumoniae were examined via genetic experiments and genomic analysis.

METHODS

Carbapenem-resistant (CR) plasmid pHSKP2 and multidrug-resistant (MDR) plasmid pHSKP3 of K. pneumoniae HS11286, virulence plasmid pRJF293 of K. pneumoniae RJF293 were employed in conjugation assays to assess the transfer ability of mobilizable plasmids. Mimic mobilizable plasmids and genetically modified plasmids were constructed to confirm the cotransfer models. The plasmid morphology was evaluated through XbaI and S1 nuclease pulsed-field gel electrophoresis and/or complete genome sequencing. Mobilizable plasmid stability in transconjugants was analyzed via serial passage culture. In addition, in silico genome analysis of 3923 plasmids of 1194 completely sequenced K. pneumoniae was performed to investigate the distribution of the conjugative elements, the cargo genes, and the targets of the CRISPR-Cas system. The mobilizable MDR plasmid and virulence plasmid of K. pneumoniae were investigated, which carry oriT but lack other conjugative elements.

RESULTS

Our results showed that mobilizable MDR and virulence plasmids carrying oriT but lacking the relaxase gene were able to cotransfer with a helper conjugative CR plasmid across various Klebsiella and Escherichia coli strains. The transfer and stability of mobilizable plasmids rather than conjugative plasmids were not interfered with by the CRISPR-Cas system of recipient strains. According to the in silico analysis, the mobilizable plasmids carry about twenty percent of acquired antimicrobial resistance genes and more than seventy-five percent of virulence genes in K. pneumoniae.

CONCLUSIONS

Our work observed that a mobilizable MDR or virulence plasmid that carries oriT but lacks the relaxase genes transferred with the helper CR conjugative plasmid and mobilizable plasmids escaped from CRISPR-Cas defence and remained stable in recipients. These results highlight the threats of mobilizable plasmids as vital vehicles in the dissemination of antibiotic resistance and virulence genes in K. pneumoniae.

Whole genome sequencing of increased number of azithromycin-resistant Shigella flexneri 1b isolates in Ontario

Azithromycin (AZM) resistance among Shigella is a major public health concern. Here, we investigated the epidemiology of Shigella flexneri serotype 1b recovered during 2016-2018 in Ontario, to describe the prevalence and spread of AZM resistance. We found that 72.3% (47/65) of cases were AZM-resistant (AZMR), of which 95.7% (45/47) were among males (P < 0.001). Whole-genome based phylogenetic analysis showed three major clusters, and 56.9% of isolates grouped within a single closely-related cluster (0-10 ∆SNP). A single AZMR clonal population was persistent over 3 years and involved 67.9% (36/53) of all male cases, and none reported international travel. In 2018, a different AZMR cluster appeared among adult males not reporting travel. A proportion of isolates (10.7%) with reduced susceptibility to ciprofloxacin (CIP) due to S83L mutation in gyrA were AZM susceptible, and 71.4% reported international travel. Resistance to AZM was due to the acquisition of mph gene-bearing incFII plasmids having > 95% nucleotide similarity to pKSR100. Plasmid-borne resistance limiting treatment options to AZM, ceftriaxone (CRO) and CIP was noted in a single isolate. We characterized AZMR isolates circulating locally among males and found that genomic analysis can support targeted prevention and mitigation strategies against antimicrobial-resistance.

Glabridin inhibited the spread of polymyxin-resistant Enterobacterium carrying ICEMmoMP63

Introduction

The role of integrative and conjugative elements (ICEs) in antibiotic resistance in Morganella morganii is unknown. This study aimed to determine whether an ICE identified in the M. morganii genome contributed to the polymyxin resistance.

Methods

Whole-genome sequencing was performed followed by bioinformatics analyses to identify ICEs and antibiotic resistance genes. Conjugation assays were performed to analyze the transferability of a discovered ICE. A drug transporter encoded on the ICE was heterogeneously expressed in Escherichia coli, minimum inhibitory concentrations of antibiotics were determined, and a traditional Chinese medicine library was screened for potential efflux pump inhibitors.

Results

An antibiotic resistance-conferring ICE, named ICEMmoMP63, was identified. ICEMmoMP63 was verified to be horizontally transferred among Enterobacteriaceae bacteria. G3577_03020 in ICEMmoMP63 was found to mediate multiple antibiotic resistances, especially polymyxin resistance. However, natural compound glabridin was demonstrated to inhibit polymyxin resistance.

Discussion

Our findings support the need for monitoring dissemination of ICEMmoMP63 in Enterobacteriaceae bacteria. Combined glabridin and polymyxin may have therapeutic potential for treating infections from multi-drug resistant bacteria carrying ICEMmoMP63.

Phylogenetic Relationships among TnpB-Containing Mobile Elements in Six Bacterial Species

Some families of mobile elements in bacterial genomes encode not only a transposase but also an accessory TnpB gene. This gene has been shown to encode an RNA-guided DNA endonuclease, co-evolving with Y1 transposase and serine recombinase in mobile elements IS605 and IS607. In this paper, we reveal the evolutionary relationships among TnpB-containing mobile elements (TCMEs) in well-assembled genomes of six bacterial species: Bacillus cereus, Clostridioides difficile, Deinococcus radiodurans, Escherichia coli, Helicobacter pylori and Salmonella enterica. In total, 9996 TCMEs were identified in 4594 genomes. They belonged to 39 different insertion sequences (ISs). Based on their genetic structures and sequence identities, the 39 TCMEs were classified into three main groups and six subgroups. According to our phylogenetic analysis, TnpBs include two main branches (TnpB-A and TnpB-B) and two minor branches (TnpB-C and TnpB-D). The key TnpB motifs and the associated Y1 and serine recombinases were highly conserved across species, even though their overall sequence identities were low. Substantial variation was observed for the rate of invasion across bacterial species and strains. Over 80% of the genomes of B. cereus, C. difficile, D. radiodurans and E. coli contained TCMEs; however, only 64% of the genomes of H. pylori and 44% of S. enterica genomes contained TCMEs. IS605 showed the largest rate of invasion in these species, while IS607 and IS1341 had a relatively narrow distribution. Co-invasions of IS605, IS607 and IS1341 elements were observed in various genomes. The largest average copy number was observed for IS605b elements in C. difficile. The average copy numbers of most other TCMEs were smaller than four. Our findings have important implications for understanding the co-evolution of TnpB-containing mobile elements and their biological roles in host genome evolution.

Genomic Epidemiological Analysis of Antimicrobial-Resistant Bacteria with Nanopore Sequencing

Antimicrobial-resistant (AMR) bacterial infections caused by clinically important bacteria, including ESKAPE pathogens (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter species) and mycobacteria (Mycobacterium tuberculosis and nontuberculous mycobacteria), have become a global public health threat. Their epidemic and pandemic clones often accumulate useful accessory genes in their genomes, such as AMR genes (ARGs) and virulence factor genes (VFGs). This process is facilitated by horizontal gene transfer among microbial communities via mobile genetic elements (MGEs), such as plasmids and phages. Nanopore long-read sequencing allows easy and inexpensive analysis of complex bacterial genome structures, although some aspects of sequencing data calculation and genome analysis methods are not systematically understood. Here we describe the latest and most recommended experimental and bioinformatics methods available for the construction of complete bacterial genomes from nanopore sequencing data and the detection and classification of genotypes of bacterial chromosomes, ARGs, VFGs, plasmids, and other MGEs based on their genomic sequences for genomic epidemiological analysis of AMR bacteria.

Identification and genetic characterization of two conjugative plasmids that confer azithromycin resistance in Salmonella

With the development of multidrug resistance in Salmonella spp. in recent years, ciprofloxacin, ceftriaxone, and azithromycin have become the principal antimicrobial agents used for the treatment of Salmonella infections. The underlying mechanisms of plasmid-mediated ciprofloxacin and ceftriaxone resistance have attracted extensive research interest, but not much is focused on azithromycin resistance in Salmonella. In this study, we investigated the genetic features of two conjugative plasmids and a non-transferable virulence plasmid that encode azithromycin resistance in food-borne Salmonella strains. We showed that the azithromycin resistance phenotype of these strains was conferred by erm(B) gene and/or the complete genetic structure IS26-mph(A)-mrx-mphR-IS6100. Comparative genetic analysis showed that these conjugative plasmids might originate from Escherichia coli and play a role in the rapid dissemination of azithromycin resistance in Salmonella. These conjugative plasmids may also serve as a reservoir of antimicrobial resistance (AMR) genes in Salmonella in which these AMR genes may be acquired by the virulence plasmids of Salmonella via genetic transposition events. Importantly, the formation of a novel macrolide-resistance and virulence-encoding plasmid, namely pS1380-118 kb, was observed in this study. This plasmid was found to exhibit transmission potential and pose a serious health threat as the extensive transmission of azithromycin resistant and virulent Salmonella strains would further compromise the effectiveness of treatment for salmonellosis. Further surveillance and research on the dissemination and evolution routes of pS1380-118kb-like plasmids in potential human pathogens of the family of Enterobacteriaceae are necessary.

Methicillin Resistant Staphylococci Isolated from Goats and Their Farm Environments in Saudi Arabia Genotypically Linked to Known Human Clinical Isolates: a Pilot Study

We conducted a pilot whole genome sequencing (WGS) study to characterize the genotypes of nine methicillin resistant staphylococci (MRS) isolates recovered from goats and their farm environments in Eastern Province, Saudi Arabia, between November 2019 to August 2020. Seven out of nine isolates were methicillin resistant Staphylococcus aureus (MRSA), and two were methicillin resistant Staphylococcus epidermidis (MRSE). All MRSA isolates possessed genotypes previously identified to infect humans, including isolates harboring ST6-SCCmec IV-t304 (n = 4), ST5-SCCmec VI- t688 (n = 2) and ST5-SCCmec V-t311 (n = 1). 2 MRSA isolates possessed plasmids that were genetically similar to those identified in S. aureus isolates recovered from humans and poultry. In contrast, plasmids found in three MRSA isolates and one MRSE isolate were genetically similar to those recovered from humans. All MRSA isolates harbored the host innate modulate genes sak and scn previously associated with human infections. The genotypes of MRSE isolates were determined as ST35, a well-known zoonotic sequence type and ST153, which has been associated with humans. However, the MRSE isolates were untypeable due to extra ccr complexes identified in their SCCmec elements. Moreover, we identified in ST153 isolate SCCmec element also harbored the Arginine Catabolic Mobile Element (ACME) IV. All MRS isolates were phenotypically resistant to trimethoprim-sulfamethoxazole, an antibiotic for the decolonization of MRS. Three isolates carried antibiotic resistance genes in their SCCmec elements that were not previously described, including those encoding fusidic acid resistance (fusC) and trimethoprim resistance (dfrC) incorporated in the MRSA SCCmec VI. IMPORTANCE Our findings demonstrate a possible cross-transmission of methicillin resistant staphylococci between goats and their local environments and between goats and humans. Due to ever increasing resistance to multiple antibiotics, the burden of MRS has a significant impact on livestock farming, public health, and the economy worldwide. This study highlights that implementing a holistic approach to whole genome sequencing surveillance in livestock and farm environments would aid our understanding of the transmission of methicillin resistant staphylococci and, most importantly, allow us to implement appropriate infection control and hygiene practices.

ExplorePipolin: reconstruction and annotation of piPolB-encoding bacterial mobile elements from draft genomes

Motivation

Detailed and accurate analysis of mobile genetic elements (MGEs) in bacteria is essential to deal with the current threat of multiresistant microbes. The overwhelming use of draft, contig-based genomes hinder the delineation of the genetic structure of these plastic and variable genomic stretches, as in the case of pipolins, a superfamily of MGEs that spans diverse integrative and plasmidic elements, characterized by the presence of a primer-independent DNA polymerase.

Results

ExplorePipolin is a Python-based pipeline that screens for the presence of the element and performs its reconstruction and annotation. The pipeline can be used on virtually any genome from diverse organisms and of diverse quality, obtaining the highest-scored possible structure and reconstructed out of different contigs if necessary. Then, predicted pipolin boundaries and pipolin encoded genes are subsequently annotated using a custom database, returning the standard file formats suitable for comparative genomics of this mobile element.

Availability and implementation

All code is available and can be accessed here: github.com/pipolinlab/ExplorePipolin.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

High Prevalence of bla CTXM-1/IncI1-Iγ/ST3 Plasmids in Extended-Spectrum β-Lactamase-Producing Escherichia coli Isolates Collected From Domestic Animals in Guadeloupe (French West Indies)

Extended-spectrum β-lactamase-producing Enterobacteriaceae (ESBL-E) have been classified in the group of resistant bacteria of highest priority. We determined the prevalence of ESBL-E collected in feces from household and shelter pets in Guadeloupe (French West Indies). A single rectal swab was taken from 125 dogs and 60 cats between June and September 2019. The prevalence of fecal carriage of ESBL-E was 7.6% (14/185, 95% CI: 4.2-12.4), within the range observed worldwide. The only risk factor associated with a higher prevalence of ESBL-E rectal carriage was a stay in a shelter, suggesting that refuges could be hotspots for their acquisition. All but one (Klebsiella pneumoniae from a cat) were Escherichia coli. We noted the presence of a bla _CTX-M-1/IncI1-Iγ/sequence type (ST3) plasmid in 11 ESBL-producing E. coli isolates belonging to ST328 (n = 6), ST155 (n = 4) and ST953 (n = 1). A bla _CTX-M-15 gene was identified in the three remaining ESBL-E isolates. The bla _CTX-M-1 and most of the antimicrobial resistance genes were present in a well-conserved large conjugative IncI1-Iγ/ST3 plasmid characterized by two accessory regions containing antibiotic resistance genes. The plasmid has been detected worldwide in E. coli isolates from humans and several animal species, such as food-producing animals, wild birds and pets, and from the environment. This study shows the potential role of pets as a reservoir of antimicrobial-resistant bacteria or genes for humans and underlines the importance of basic hygiene measures by owners of companion animals.

Identification of Promoter Region Markers Associated With Altered Expression of Resistance-Nodulation-Division Antibiotic Efflux Pumps in Acinetobacter baumannii

Genetic alterations leading to the constitutive upregulation of specific efflux pumps contribute to antibacterial resistance in multidrug resistant bacteria. The identification of such resistance markers remains one of the most challenging tasks of genome-level resistance predictors. In this study, 487 non-redundant genetic events were identified in upstream zones of three operons coding for resistance-nodulation-division (RND) efflux pumps of 4,130 Acinetobacter baumannii isolates. These events included insertion sequences, small indels, and single nucleotide polymorphisms. In some cases, alterations explicitly modified the expression motifs described for these operons, such as the promoter boxes, operators, and Shine-Dalgarno sequences. In addition, changes in DNA curvature and mRNA secondary structures, which are structural elements that regulate expression, were also calculated. According to their influence on RND upregulation, the catalog of upstream modifications were associated with “experimentally verified,” “presumed,” and “probably irrelevant” degrees of certainty. For experimental verification, DNA of upstream sequences independently carrying selected markers, three for each RND operon, were fused to a luciferase reporter plasmid system. Five out of the nine selected markers tested showed significant increases in expression with respect to the wild-type sequence control. In particular, a 25-fold expression increase was observed with the ISAba1 insertion sequence upstream the adeABC pump. Next, overexpression of each of the three multi-specific RND pumps was linked to their respective antibacterial substrates by a deep A. baumannii literature screen. Consequently, a data flow framework was then developed to link genomic upregulatory RND determinants to potential antibiotic resistance. Assignment of potential increases in minimal inhibitory concentrations at the “experimentally verified” level was permitted for 42 isolates to 7–8 unrelated antibacterial agents including tigecycline, which is overlooked by conventional resistome predictors. Thus, our protocol may represent a time-saving filter step prior to laborious confirmation experiments for efflux-driven resistance. Altogether, a computational-experimental pipeline containing all components required for identifying the upstream regulatory resistome is proposed. This schema may provide the foundational stone for the elaboration of tools approaching antibiotic efflux that complement routine resistome predictors for preventing antimicrobial therapy failure against difficult-to-threat bacteria.

The Integrative and Conjugative Element ICECspPOL2 Contributes to the Outbreak of Multi-Antibiotic-Resistant Bacteria for Chryseobacterium Spp. and Elizabethkingia Spp

Antibiotic resistance genes (ARGs) and horizontal transfer of ARGs among bacterial species in the environment can have serious clinical implications as such transfers can lead to disease outbreaks from multidrug-resistant (MDR) bacteria. Infections due to antibiotic-resistant Chryseobacterium and Elizabethkingia in intensive care units have been increasing in recent years. In this study, the multi-antibiotic-resistant strain Chryseobacterium sp. POL2 was isolated from the wastewater of a livestock farm. Whole-genome sequencing and annotation revealed that the POL2 genome encodes dozens of ARGs. The integrative and conjugative element (ICE) ICECspPOL2, which encodes ARGs associated with four types of antibiotics, including carbapenem, was identified in the POL2 genome, and phylogenetic affiliation analysis suggested that ICECspPOL2 evolved from related ICEEas of Elizabethkingia spp. Conjugation assays verified that ICECspPOL2 can horizontally transfer to Elizabethkingia species, suggesting that ICECspPOL2 contributes to the dissemination of multiple ARGs among Chryseobacterium spp. and Elizabethkingia spp. Because Elizabethkingia spp. is associated with clinically significant infections and high mortality, there would be challenges to clinical treatment if these bacteria acquire ICECspPOL2 with its multiple ARGs, especially the carbapenem resistance gene. Therefore, the results of this study support the need for monitoring the dissemination of this type of ICE in Chryseobacterium and Elizabethkingia strains to prevent further outbreaks of MDR bacteria. IMPORTANCE Infections with multiple antibiotic-resistant Chryseobacterium and Elizabethkingia in intensive care units have been increasing in recent years. In this study, the mobile integrative and conjugative element ICECspPOL2, which was associated with the transmission of a carbapenem resistance gene, was identified in the genome of the multi-antibiotic-resistant strain Chryseobacterium sp. POL2. ICECspPOL2 is closely related to the ICEEas from Elizabethkingia species, and ICECspPOL2 can horizontally transfer to Elizabethkingia species with the tRNA-Glu-TTC gene as the insertion site. Because Elizabethkingia species are associated with clinically significant infections and high mortality, the ability of ICECspPOL2 to transfer carbapenem resistance from environmental strains of Chryseobacterium to Elizabethkingia is of clinical concern.

Characterization of blaKPC-2-Carrying Plasmid pR31-KPC from a Pseudomonas aeruginosa Strain Isolated in China

This work aimed to characterize a 29-kb bla_KPC-2-carrying plasmid, pR31-KPC, from a multidrug resistant strain of Pseudomonas aeruginosa isolated from the sputum of an elderly patient with multiple chronic conditions in China. The backbone of pR31-KPC is closely related to four other bla_KPC-2-carrying plasmids, YLH6_p3, p1011-KPC2, p14057A, and pP23-KPC, none of which have been assigned to any of the known incompatibility groups. Two accessory modules, the IS26-bla_KPC-2-IS26 unit and IS26-ΔTn6376-IS26 region, separated by a 5.9-kb backbone region, were identified in pR31-KPC, which was also shown to carry the unique resistance marker bla_KPC-2. A comparative study of the above five plasmids showed that p1011-KPC2 may be the most complete plasmid of this group to be reported, while pR31-KPC is the smallest plasmid having lost most of its conjugative region. Regions between the iterons and orf207 in the backbone may be hot spots for the acquisition of exogenous resistance entities. The accessory regions of these plasmids have all undergone several biological events when compared with Tn6296. The further transfer of bla_KPC-2 in these plasmids may be initiated by either the Tn3 family or IS26-associated transposition or homologous recombination. The data presented here will contribute to a deeper understanding of bla_KPC-2 carrying plasmids in Pseudomonas.

Development of a multiplex-PCR serotyping assay for characterizing Legionella pneumophila serogroups based on the diversity of LPS biosynthetic loci

Legionella pneumophila, which is the main cause of Legionnaires' disease, comprises at least 15 serogroups (SGs). We show here the diversity of LPS biosynthetic loci among serogroups and describe the development of a PCR serotyping assay for 15 SGs based on the sequences of LPS biosynthetic loci. Using this multiplex-PCR (M-PCR) system, serogroup(s) were detected using primers that specifically amplify the sequences of SG1, SG2, SG5, SG7, SG8, SG9, SG11, SG13, SG3/15, and SG6/12. When PCR products of the expected sizes were not detected, we used primers that identified SG4/10/14. The PCR serotyping system specifically amplified the sequences corresponding SGs of 238 L. pneumophila strains. This method will be very useful for conducting epidemiological studies and investigating outbreak of Legionnaires' disease.

Sequence analysis of tyrosine recombinases allows annotation of mobile genetic elements in prokaryotic genomes

Mobile genetic elements (MGEs) sequester and mobilize antibiotic resistance genes across bacterial genomes. Efficient and reliable identification of such elements is necessary to follow resistance spreading. However, automated tools for MGE identification are missing. Tyrosine recombinase (YR) proteins drive MGE mobilization and could provide markers for MGE detection, but they constitute a diverse family also involved in housekeeping functions. Here, we conducted a comprehensive survey of YRs from bacterial, archaeal, and phage genomes and developed a sequence-based classification system that dissects the characteristics of MGE-borne YRs. We revealed that MGE-related YRs evolved from non-mobile YRs by acquisition of a regulatory arm-binding domain that is essential for their mobility function. Based on these results, we further identified numerous unknown MGEs. This work provides a resource for comparative analysis and functional annotation of YRs and aids the development of computational tools for MGE annotation. Additionally, we reveal how YRs adapted to drive gene transfer across species and provide a tool to better characterize antibiotic resistance dissemination.

A Novel Mobile Element ICERspD18B in Rheinheimera sp. D18 Contributes to Antibiotic and Arsenic Resistance

Antibiotics and organoarsenical compounds are frequently used as feed additives in many countries. However, these compounds can cause serious antibiotic and arsenic (As) pollution in the environment, and the spread of antibiotic and As resistance genes from the environment. In this report, we characterized the 28.5 kb genomic island (GI), named as ICERspD18B, as a novel chromosomal integrative and conjugative element (ICE) in multidrug-resistant Rheinheimera sp. D18. Notably, ICERspD18B contains six antibiotic resistance genes (ARGs) and an arsenic tolerance operon, as well as genes encoding conjugative transfer proteins of a type IV secretion system, relaxase, site-specific integrase, and DNA replication or partitioning proteins. The transconjugant strain 25D18-B4 was generated using Escherichia coli 25DN as the recipient strain. ICERspD18B was inserted into 3'-end of the guaA gene in 25D18-B4. In addition, 25D18-B4 had markedly higher minimum inhibitory concentrations for arsenic compounds and antibiotics when compared to the parental E. coli strain. These findings demonstrated that the integrative and conjugative element ICERspD18B could mediate both antibiotic and arsenic resistance in Rheinheimera sp. D18 and the transconjugant 25D18-B4.

Emergence and Spread of Different ESBL-Producing Salmonella enterica Serovars in Hospitalized Horses Sharing a Highly Transferable IncM2 CTX-M-3-Encoding Plasmid

Salmonella enterica is a major causative pathogen of human and animal gastroenteritis. Antibiotic resistant strains have emerged due to the production of extended-spectrum β-lactamases (ESBLs) posing a major health concern. With the increasing reports on ESBL-producing Enterobacterales that colonize companion animals, we aimed to investigate ESBL dissemination among ESBL-producing Salmonella enterica (ESBL-S) in hospitalized horses. We prospectively collected ESBL-S isolates from hospitalized horses in a Veterinary-Teaching Hospital during Dec 2015-Dec 2017. Selection criteria for ESBL-S were white colonies on CHROMagarESBL plates and an ESBL phenotypic confirmation. Salmonella enterica serovars were determined using the Kaufmann-White-Le-Minor serological scheme. ESBL-encoding plasmids were purified, transformed and compared using restriction fragment length polymorphism (RFLP). Whole genome sequencing (Illumina and MinION platforms) were performed for detailed phylogenetic and plasmid analyses. Twelve ESBL-S were included in this study. Molecular investigation and Sequence Read Archive (SRA) meta-analysis revealed the presence of three unique Salmonella enterica serovars, Cerro, Havana and Liverpool, all reported for the first time in horses. PFGE revealed the clonal spread of S. Cerro between seven horses. All twelve isolates carried bla CTX-M- 3 and showed an identical multidrug resistance profile with co-resistance to trimethoprim/sulfamethoxazole and to aminoglycosides. Plasmid RFLP proved the inter-serovar horizontal spread of a single bla CTX-M- 3-encoding plasmid. Complete sequence of a representative plasmid (S. Havana strain 373.3.1), designated pSEIL-3 was a -86.4 Kb IncM2 plasmid, that encoded nine antibiotic resistance genes. pSEIL-3 was virtually identical to pCTX-M3 from Citrobacter freundii, and showed high identity (>95%) to six other bla CTX-M- 3 or bla NDM- 1 IncM2 broad host range plasmids from various Enterobacterales of human origin. Using a specific six gene-based multiplex PCR, we detected pSEIL-3 in various Enterobacterales species that co-colonized the horses' gut. Together, our findings show the alarming emergence of ESBL-S in hospitalized horses associated with gut shedding and foal morbidity and mortality. We demonstrated the dissemination of CTX-M-3 ESBL among different Salmonella enterica serovars due to transmission of a broad host range plasmid. This report highlights horses as a zoonotic reservoir for ESBL-S, including highly transmissible plasmids that may represent a 'One-Health' hazard. This risk calls for the implementation of infection control measures to monitor and control the spread of ESBL-S in hospitalized horses.

Web-Based Genome Analysis of Bacterial Meningitis Pathogens for Public Health Applications Using the Bacterial Meningitis Genomic Analysis Platform (BMGAP)

Effective laboratory-based surveillance and public health response to bacterial meningitis depends on timely characterization of bacterial meningitis pathogens. Traditionally, characterizing bacterial meningitis pathogens such as Neisseria meningitidis (Nm) and Haemophilus influenzae (Hi) required several biochemical and molecular tests. Whole genome sequencing (WGS) has enabled the development of pipelines capable of characterizing the given pathogen with equivalent results to many of the traditional tests. Here, we present the Bacterial Meningitis Genomic Analysis Platform (BMGAP): a secure, web-accessible informatics platform that facilitates automated analysis of WGS data in public health laboratories. BMGAP is a pipeline comprised of several components, including both widely used, open-source third-party software and customized analysis modules for the specific target pathogens. BMGAP performs de novo draft genome assembly and identifies the bacterial species by whole-genome comparisons against a curated reference collection of 17 focal species including Nm, Hi, and other closely related species. Genomes identified as Nm or Hi undergo multi-locus sequence typing (MLST) and capsule characterization. Further typing information is captured from Nm genomes, such as peptides for the vaccine antigens FHbp, NadA, and NhbA. Assembled genomes are retained in the BMGAP database, serving as a repository for genomic comparisons. BMGAP's species identification and capsule characterization modules were validated using PCR and slide agglutination from 446 bacterial invasive isolates (273 Nm from nine different serogroups, 150 Hi from seven different serotypes, and 23 from nine other species) collected from 2017 to 2019 through surveillance programs. Among the validation isolates, BMGAP correctly identified the species for all 440 isolates (100% sensitivity and specificity) and accurately characterized all Nm serogroups (99% sensitivity and 98% specificity) and Hi serotypes (100% sensitivity and specificity). BMGAP provides an automated, multi-species analysis pipeline that can be extended to include additional analysis modules as needed. This provides easy-to-interpret and validated Nm and Hi genome analysis capacity to public health laboratories and collaborators. As the BMGAP database accumulates more genomic data, it grows as a valuable resource for rapid comparative genomic analyses during outbreak investigations.

Incorporating the plasmidome into antibiotic resistance surveillance in animal agriculture

Mobile genetic elements (MGE) carrying resistance genes represent a unique challenge to risk assessment and surveillance of antimicrobial resistance (AMR). Yet determining the mobility of resistance genes within animal microbiomes is essential to evaluating the potential dissemination from livestock to potential human pathogens, as well as evaluating co-selection mechanisms that may impact persistence of resistance genes with changing antibiotic use patterns. Current surveillance efforts utilize phenotypic testing and sequencing of individual isolates for tracking of AMR in livestock. In this work, we investigated the utility of using long-read sequencing of the plasmids from mixed Enterobacterales enrichments of swine fecal samples as a surveillance strategy for AMR plasmids. Enrichments were performed in either MacConkey broth without selection or with selection by addition of tetracycline or ceftriaxone, and plasmids were extracted and sequenced in order to evaluate the diversity of plasmids enriched by each method. Intact resistance plasmids were successfully assembled, as well as complex resistance transposons carrying multiple repeated elements that would interfere with assembly by short read sequencing technologies. Comparison of the assembled plasmids with representatives from public databases confirmed the quality of the assemblies and also revealed the occurrence of IncI2 plasmids carrying bla_CMY-2 in Ontario swine samples, which have not been found in previous studies.

First report of plasmid-mediated colistin resistance mcr-8.1 gene from a clinical Klebsiella pneumoniae isolate from Lebanon

Colistin is considered as a last resort treatment for infections caused by multidrug-resistant Enterobacteriaceae. Plasmid-mediated mobile colistin resistance (mcr) genes contributed to the global spread of colistin resistance. This is the first report of plasmid-mediated colistin resistance mcr-8 gene from a clinical Klebsiella pneumoniae K9 isolate recovered from Lebanon. The isolate was characterized phenotypically and genotypically through both short and long read whole-genome sequencing, plasmid typing and conjugation assays. k9 belonged to sequence type 15 and harbored 31 antimicrobial resistance genes. The mcr-8.1 variant was carried on a novel ~ 300 kb multireplicon plasmid having IncFIA, IncR and IncHI1B. The plasmid was conjugative and carried a plethora of antimicrobial resistance determinants. The introduction of novel mcr variants in Lebanon poses an alarming health concern. Surveillance and screening for colistin resistant Enterobacteriaceae and mcr in livestock, animal farms, imported meat and poultry is highly recommended along with monitoring antibiotic use.

Evaluating Eucalyptus leaf colonization by Brasilonema octagenarum (Cyanobacteria, Scytonemataceae) using in planta experiments and genomics

Background

Brasilonema is a cyanobacterial genus found on the surface of mineral substrates and plants such as bromeliads, orchids and eucalyptus. B. octagenarum stands out among cyanobacteria due to causing damage to the leaves of its host in an interaction not yet observed in other cyanobacteria. Previous studies revealed that B. octagenaum UFV-E1 is capable of leading eucalyptus leaves to suffer internal tissue damage and necrosis by unknown mechanisms. This work aimed to investigate the effects of B. octagenarum UFV-E1 inoculation on Eucalyptus urograndis and to uncover molecular mechanisms potentially involved in leaf damage by these cyanobacteria using a comparative genomics approach.

Results

Leaves from E. urograndis saplings were exposed for 30 days to B. octagenarum UFV-E1, which was followed by the characterization of its genome and its comparison with the genomes of four other Brasilonema strains isolated from phyllosphere and the surface of mineral substrates. While UFV-E1 inoculation caused an increase in root and stem dry mass of the host plants, the sites colonized by cyanobacteria on leaves presented a significant decrease in pigmentation, showing that the cyanobacterial mats have an effect on leaf cell structure. Genomic analyses revealed that all evaluated Brasilonema genomes harbored genes encoding molecules possibly involved in plant-pathogen interactions, such as hydrolases targeting plant cell walls and proteins similar to known virulence factors from plant pathogens. However, sequences related to the type III secretory system and effectors were not detected, suggesting that, even if any virulence factors could be expressed in contact with their hosts, they would not have the structural means to actively reach plant cytoplasm.

Conclusions

Leaf damage by this species is likely related to the blockage of access to sunlight by the efficient growth of cyanobacterial mats on the phyllosphere, which may hinder the photosynthetic machinery and prevent access to some essential molecules. These results reveal that the presence of cyanobacteria on leaf surfaces is not as universally beneficial as previously thought, since they may not merely provide the products of nitrogen fixation to their hosts in exchange for physical support, but in some cases also hinder regular leaf physiology leading to tissue damage.

New candidates for regulated gene integrity revealed through precise mapping of integrative genetic elements

Integrative genetic elements (IGEs) are mobile multigene DNA units that integrate into and excise from host bacterial genomes. Each IGE usually targets a specific site within a conserved host gene, integrating in a manner that preserves target gene function. However, a small number of bacterial genes are known to be inactivated upon IGE integration and reactivated upon excision, regulating phenotypes of virulence, mutation rate, and terminal differentiation in multicellular bacteria. The list of regulated gene integrity (RGI) cases has been slow-growing because IGEs have been challenging to precisely and comprehensively locate in genomes. We present software (TIGER) that maps IGEs with unprecedented precision and without attB site bias. TIGER uses a comparative genomic, ping-pong BLAST approach, based on the principle that the IGE integration module (i.e. its int-attP region) is cohesive. The resultant IGEs from 2168 genomes, along with integrase phylogenetic analysis and gene inactivation tests, revealed 19 new cases of genes whose integrity is regulated by IGEs (including dut, eccCa1, gntT, hrpB, merA, ompN, prkA, tqsA, traG, yifB, yfaT and ynfE), as well as recovering previously known cases (in sigK, spsM, comK, mlrA and hlb genes). It also recovered known clades of site-promiscuous integrases and identified possible new ones.

Whole-genome comparison between the type strain of Halobacterium salinarum (DSM 3754T ) and the laboratory strains R1 and NRC-1

Halobacterium salinarum is an extremely halophilic archaeon that is widely distributed in hypersaline environments and was originally isolated as a spoilage organism of salted fish and hides. The type strain 91-R6 (DSM 3754T ) has seldom been studied and its genome sequence has only recently been determined by our group. The exact relationship between the type strain and two widely used model strains, NRC-1 and R1, has not been described before. The genome of Hbt. salinarum strain 91-R6 consists of a chromosome (2.17 Mb) and two large plasmids (148 and 102 kb, with 39,230 bp being duplicated). Cytosine residues are methylated (m4 C) within CTAG motifs. The genomes of type and laboratory strains are closely related, their chromosomes sharing average nucleotide identity (ANIb) values of 98% and in silico DNA-DNA hybridization (DDH) values of 95%. The chromosomes are completely colinear, do not show genome rearrangement, and matching segments show <1% sequence difference. Among the strain-specific sequences are three large chromosomal replacement regions (>10 kb). The well-studied AT-rich island (61 kb) of the laboratory strains is replaced by a distinct AT-rich sequence (47 kb) in 91-R6. Another large replacement (91-R6: 78 kb, R1: 44 kb) codes for distinct homologs of proteins involved in motility and N-glycosylation. Most (107 kb) of plasmid pHSAL1 (91-R6) is very closely related to part of plasmid pHS3 (R1) and codes for essential genes (e.g. arginine-tRNA ligase and the pyrimidine biosynthesis enzyme aspartate carbamoyltransferase). Part of pHS3 (42.5 kb total) is closely related to the largest strain-specific sequence (164 kb) in the type strain chromosome. Genome sequencing unraveled the close relationship between the Hbt. salinarum type strain and two well-studied laboratory strains at the DNA and protein levels. Although an independent isolate, the type strain shows a remarkably low evolutionary difference to the laboratory strains.

Meta-analysis of Pandemic Escherichia coli ST131 Plasmidome Proves Restricted Plasmid-clade Associations

Extraintestinal multidrug resistant Escherichia coli sequence type (ST) 131 is a worldwide pandemic pathogen and a major cause of urinary tract and bloodstream infections. The role of this pandemic lineage in multidrug resistance plasmid dissemination is still scarce. We herein performed a meta-analysis on E. coli ST131 whole-genome sequence (WGS) databases to unravel ST131 plasmidome and specifically to decipher CTX-M encoding plasmids-clade associations. We mined 880 ST131 WGS data and proved that CTX-M-27-encoding IncF[F1:A2:B20] (Group1) plasmids are strictly found in clade C1, whereas CTX-M-15-encoding IncF[F2:A1:B-] (Group2) plasmids exist only in clade C2 suggesting strong plasmid-clade adaptations. Specific Col-like replicons (Col156, Col(MG828), and Col8282) were also found to be clade C1-associated. BLAST-based search revealed that Group1 and Group2 plasmids are narrow-host-range and restricted to E.coli. Among a collection of 20 newly sequenced Israeli ST131 CTX-M-encoding plasmids (2003–2016), Group1 and Group2 plasmids were dominant and associated with the expected clades. We found, for the first time in ST131, a CTX-M-15-encoding phage-like plasmid group (Group3) and followed its spread in the WGS data. This study offers a comprehensive way to decipher plasmid-bacterium associations and demonstrates that the CTX-M-encoding ST131 Group1 and Group2 plasmids are clade-restricted and presumably less transmissible, potentially contributing to ST131 clonal superiority.

Mobilome and Resistome Reconstruction from Genomes Belonging to Members of the Bifidobacterium Genus

Specific members of the genus Bifidobacterium are among the first colonizers of the human/animal gut, where they act as important intestinal commensals associated with host health. As part of the gut microbiota, bifidobacteria may be exposed to antibiotics, used in particular for intrapartum prophylaxis, especially to prevent Streptococcus infections, or in the very early stages of life after the birth. In the current study, we reconstructed the in silico resistome of the Bifidobacterium genus, analyzing a database composed of 625 bifidobacterial genomes, including partial assembled strains with less than 100 genomic sequences. Furthermore, we screened bifidobacterial genomes for mobile genetic elements, such as transposases and prophage-like elements, in order to investigate the correlation between the bifido-mobilome and the bifido-resistome, also identifying genetic insertion hotspots that appear to be prone to horizontal gene transfer (HGT) events. These insertion hotspots were shown to be widely distributed among analyzed bifidobacterial genomes, and suggest the acquisition of antibiotic resistance genes through HGT events. These data were further corroborated by growth experiments directed to evaluate bacitracin A resistance in Bifidobacterium spp., a property that was predicted by in silico analyses to be part of the HGT-acquired resistome.

RefSeq: an update on prokaryotic genome annotation and curation

The Reference Sequence (RefSeq) project at the National Center for Biotechnology Information (NCBI) provides annotation for over 95 000 prokaryotic genomes that meet standards for sequence quality, completeness, and freedom from contamination. Genomes are annotated by a single Prokaryotic Genome Annotation Pipeline (PGAP) to provide users with a resource that is as consistent and accurate as possible. Notable recent changes include the development of a hierarchical evidence scheme, a new focus on curating annotation evidence sources, the addition and curation of protein profile hidden Markov models (HMMs), release of an updated pipeline (PGAP-4), and comprehensive re-annotation of RefSeq prokaryotic genomes. Antimicrobial resistance proteins have been reannotated comprehensively, improved structural annotation of insertion sequence transposases and selenoproteins is provided, curated complex domain architectures have given upgraded names to millions of multidomain proteins, and we introduce a new kind of annotation rule—BlastRules. Continual curation of supporting evidence, and propagation of improved names onto RefSeq proteins ensures that the functional annotation of genomes is kept current. An increasing share of our annotation now derives from HMMs and other sets of annotation rules that are portable by nature, and available for download and for reuse by other investigators. RefSeq is found at https://www.ncbi.nlm.nih.gov/refseq/.

Transposition of IS4 Family Insertion Sequences ISTeha3, ISTeha4, and ISTeha5 into the arc Operon Disrupts Arginine Deiminase System in Tetragenococcus halophilus

Tetragenococcus halophilus, a halophilic lactic acid bacterium, is often used as a starter culture in the manufacturing of soy sauce. T. halophilus possesses an arginine deiminase system, which is responsible for the accumulation of citrulline, the main precursor of the potential carcinogen ethyl carbamate. In this study, we generated five derivatives lacking arginine deiminase activity from T. halophilus NBRC 12172 by UV irradiation. Using these derivatives as a fermentation starter prevented arginine deimination in soy sauce. DNA sequence analysis of the derivatives revealed that novel IS4 family insertion sequences, designated ISTeha3, ISTeha4, and ISTeha5, were transposed into the region around the arginine deiminase (arc) operon in the mutants. These insertion sequences contain a single open reading frame encoding a putative transposase and 13- to 15-bp inverted repeats at both termini, which are adjacent to 7- to 9-bp duplications of the target sequence. Investigation of wild strains isolated from soy sauce mash incapable of arginine deimination also indicated that insertion sequences are involved in the disruption of the arginine deiminase system in T. halophilus IMPORTANCE Insertion sequences play important roles in bacterial evolution and are frequently utilized in mutagenesis systems. However, the intrinsic insertion sequences of tetragenococci are not well characterized. Here, we identified three active insertion sequences of T. halophilus by transposition into the region around the arc operon. This report provides an example of insertion sequence-mediated generation and evolution of T. halophilus and primary information about their characteristics.

CRISPR tracking reveals global spreading of antimicrobial resistance genes by Staphylococcus of canine origin

The close contact between pets and their owners is a potential source for microorganisms and genetic material exchange. Staphylococcus species considered as harmless inhabitants of animals' and humans' microbiota can act as reservoirs of antimicrobial resistance genes to more virulent species, thereby increasing their potential to resist drug therapy. This process could be inhibited by the antiplasmid immunity conferred by CRISPR systems. On the other hand, CRISPR spacer sequences can be explored as molecular clocks to track the history of genetic invasion suffered by a bacterial strain. To understand better the role of domestic dogs in human health as an antimicrobial resistance genes source, we analyzed 129 genomes of Staphylococcus strains of canine origin for the presence of CRISPR systems. Only 8% of the strains were positive for CRISPR, which is consistent with Staphylococcus role as gene reservoirs. The plasmidial origin or some spacers confirms the unsuccessful attempt of plasmid exchange in strains carrying CRISPRs. Some of these systems are within a staphylococcal cassette chromosome mec (SCCmec), sharing 98% of identity between their harboring strains. These CRISPRs' spacers reveal that this SCCmec was transferred between canine S. pseudintermedius strains, then to S. schleiferi and to Staphylococcus strains isolated from human beings. Our findings shows genetic evidence for the global spreading of pathogenic bacteria and the antimicrobial resistance genes carried by them and reinforce that, in the age of antimicrobial resistance, it is imperative that drug therapies consider the integrated nature of the relationship between pets and humans.

Genome plasticity favours double chromosomal Tn4401b-blaKPC-2 transposon insertion in the Pseudomonas aeruginosa ST235 clone

BACKGROUND

Pseudomonas aeruginosa Sequence Type 235 is a clone that possesses an extraordinary ability to acquire mobile genetic elements and has been associated with the spread of resistance genes, including genes that encode for carbapenemases. Here, we aim to characterize the genetic platforms involved in resistance dissemination in blaKPC-2-positive P. aeruginosa ST235 in Colombia.

RESULTS

In a prospective surveillance study of infections in adult patients attended in five ICUs in five distant cities in Colombia, 58 isolates of P. aeruginosa were recovered, of which, 27 (46.6%) were resistant to carbapenems. The molecular analysis showed that 6 (22.2%) and 4 (14.8%) isolates harboured the blaVIM and blaKPC-2 genes, respectively. The four blaKPC-2-positive isolates showed a similar PFGE pulsotype and belonged to ST235. Complete genome sequencing of a representative ST235 isolate shows a unique chromosomal contig of 7097.241 bp with eight different resistance genes identified and five transposons: a Tn6162-like with ant(2″)-Ia, two Tn402-like with ant(3″)-Ia and blaOXA-2 and two Tn4401b with blaKPC-2. All transposons were inserted into the genomic islands. Interestingly, the two Tn4401b copies harbouring blaKPC-2 were adjacently inserted into a new genomic island (PAGI-17) with traces of a replicative transposition process. This double insertion was probably driven by several structural changes within the chromosomal region containing PAGI-17 in the ST235 background.

CONCLUSION

This is the first report of a double Tn4401b chromosomal insertion in P. aeruginosa, just within a new genomic island (PAGI-17). This finding indicates once again the great genomic plasticity of this microorganism.

Phenotypic and genomic comparison of Photorhabdus luminescens subsp. laumondii TT01 and a widely used rifampicin-resistant Photorhabdus luminescens laboratory strain

BACKGROUND

Photorhabdus luminescens is an enteric bacterium, which lives in mutualistic association with soil nematodes and is highly pathogenic for a broad spectrum of insects. A complete genome sequence for the type strain P. luminescens subsp. laumondii TT01, which was originally isolated in Trinidad and Tobago, has been described earlier. Subsequently, a rifampicin resistant P. luminescens strain has been generated with superior possibilities for experimental characterization. This strain, which is widely used in research, was described as a spontaneous rifampicin resistant mutant of TT01 and is known as TT01-Rif^R.

RESULTS

Unexpectedly, upon phenotypic comparison between the rifampicin resistant strain and its presumed parent TT01, major differences were found with respect to bioluminescence, pigmentation, biofilm formation, haemolysis as well as growth. Therefore, we renamed the strain TT01-Rif^R to DJC. To unravel the genomic basis of the observed differences, we generated a complete genome sequence for strain DJC using the PacBio long read technology. As strain DJC was supposed to be a spontaneous mutant, only few sequence differences were expected. In order to distinguish these from potential sequencing errors in the published TT01 genome, we re-sequenced a derivative of strain TT01 in parallel, also using the PacBio technology. The two TT01 genomes differed at only 30 positions. In contrast, the genome of strain DJC varied extensively from TT01, showing 13,000 point mutations, 330 frameshifts, and 220 strain-specific regions with a total length of more than 300 kb in each of the compared genomes.

CONCLUSIONS

According to the major phenotypic and genotypic differences, the rifampicin resistant P. luminescens strain, now named strain DJC, has to be considered as an independent isolate rather than a derivative of strain TT01. Strains TT01 and DJC both belong to P. luminescens subsp. laumondii.

Single-strand DNA processing: phylogenomics and sequence diversity of a superfamily of potential prokaryotic HuH endonucleases

Background

Some mobile genetic elements target the lagging strand template during DNA replication. Bacterial examples are insertion sequences IS608 and ISDra2 (IS200/IS605 family members). They use obligatory single-stranded circular DNA intermediates for excision and insertion and encode a transposase, TnpA_IS200, which recognizes subterminal secondary structures at the insertion sequence ends. Similar secondary structures, Repeated Extragenic Palindromes (REP), are present in many bacterial genomes. TnpA_IS200-related proteins, TnpA_REP, have been identified and could be responsible for REP sequence proliferation. These proteins share a conserved HuH/Tyrosine core domain responsible for catalysis and are involved in processes of ssDNA cleavage and ligation. Our goal is to characterize the diversity of these proteins collectively referred as the TnpA_Y1 family.

Results

A genome-wide analysis of sequences similar to TnpA_IS200 and TnpA_REP in prokaryotes revealed a large number of family members with a wide taxonomic distribution. These can be arranged into three distinct classes and 12 subclasses based on sequence similarity. One subclass includes sequences similar to TnpA_IS200. Proteins from other subclasses are not associated with typical insertion sequence features. These are characterized by specific additional domains possibly involved in protein/DNA or protein/protein interactions. Their genes are found in more than 25% of species analyzed. They exhibit a patchy taxonomic distribution consistent with dissemination by horizontal gene transfers followed by loss. The tnpA_REP genes of five subclasses are flanked by typical REP sequences in a REPtron-like arrangement. Four distinct REP types were characterized with a subclass specific distribution. Other subclasses are not associated with REP sequences but have a large conserved domain located in C-terminal end of their sequence. This unexpected diversity suggests that, while most likely involved in processing single-strand DNA, proteins from different subfamilies may play a number of different roles.

Conclusions

We established a detailed classification of TnpA_Y1 proteins, consolidated by the analysis of the conserved core domains and the characterization of additional domains. The data obtained illustrate the unexpected diversity of the TnpA_Y1 family and provide a strong framework for future evolutionary and functional studies. By their potential function in ssDNA editing, they may confer adaptive responses to host cell physiology and metabolism.

Electronic supplementary material

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

Genomic Epidemiology of NDM-1-Encoding Plasmids in Latin American Clinical Isolates Reveals Insights into the Evolution of Multidrug Resistance

Bacteria that produce the broad-spectrum Carbapenem antibiotic New Delhi Metallo-β-lactamase (NDM) place a burden on health care systems worldwide, due to the limited treatment options for infections caused by them and the rapid global spread of this antibiotic resistance mechanism. Although it is believed that the associated resistance gene bla_NDM-1 originated in Acinetobacter spp., the role of Enterobacteriaceae in its dissemination remains unclear. In this study, we used whole genome sequencing to investigate the dissemination dynamics of bla_NDM-1-positive plasmids in a set of 21 clinical NDM-1-positive isolates from Colombia and Mexico (Providencia rettgeri, Klebsiella pneumoniae, and Acinetobacter baumannii) as well as six representative NDM-1-positive Escherichia coli transconjugants. Additionally, the plasmids from three representative P. rettgeri isolates were sequenced by PacBio sequencing and finished. Our results demonstrate the presence of previously reported plasmids from K. pneumoniae and A. baumannii in different genetic backgrounds and geographically distant locations in Colombia. Three new previously unclassified plasmids were also identified in P. rettgeri from Colombia and Mexico, plus an interesting genetic link between NDM-1-positive P. rettgeri from distant geographic locations (Canada, Mexico, Colombia, and Israel) without any reported epidemiological links was discovered. Finally, we detected a relationship between plasmids present in P. rettgeri and plasmids from A. baumannii and K. pneumoniae. Overall, our findings suggest a Russian doll model for the dissemination of bla_NDM-1 in Latin America, with P. rettgeri playing a central role in this process, and reveal new insights into the evolution and dissemination of plasmids carrying such antibiotic resistance genes.

Complete genome sequencing and genomic characterization of two Escherichia coli strains co-producing MCR-1 and NDM-1 from bloodstream infection

We previously described the discovery of two Escherichia coli isolates (EC1002 and EC2474) co-harbouring mcr-1 and bla_NDM-1 genes, which were recovered from bloodstream infection in China. More importantly, these antibiotic resistance genes were located on different plasmids and signaling the potential spread of pandrug-resistant bacteria. Here, the complete genome sequences of both isolates were determined using Pacbio RS II and Illumina HiSeq2000 systems. The genome of EC1002 consists of a 5,177,501 base pair chromosome and four circular plasmids, while the genome of EC2474 consists of a 5,013,813 base pair chromosome and three plasmids. The plasmid replicon type of pEC1002_NDM and pEC2474_NDM were identified as IncA/C2 and IncF, respectively. The genetic environment of bla_NDM-1 in this study was similar to bla_NDM-carrying plasmids detected in China, although the overall nucleotide identity and query coverage were variable. The plasmid replicon type of pEC1002_MCR and pEC2474_MCR were identified as IncI2 and IncHI2, respectively. Two different genetic strategies for mcr-1 gene spread were observed in this study and bla_NDM-1 genes were also found transferred by two different mobile genetic elements in two plasmids. The findings of this study further support that the diversified transfer mechanisms of bla_NDM-1 and mcr-1 present in Enterobacteriaceae.

What matters in chronic Burkholderia cenocepacia infection in cystic fibrosis: Insights from comparative genomics

Burkholderia cenocepacia causes severe pulmonary infections in cystic fibrosis (CF) patients. Since the bacterium is virtually untreatable by antibiotics, chronic infections persist for years and might develop into fatal septic pneumonia (cepacia syndrome, CS). To devise new strategies to combat chronic B. cenocepacia infections, it is essential to obtain comprehensive knowledge about their pathogenesis. We conducted a comparative genomic analysis of 32 Czech isolates of epidemic clone B. cenocepacia ST32 isolated from various stages of chronic infection in 8 CF patients. High numbers of large-scale deletions were found to occur during chronic infection, affecting preferentially genomic islands and nonessential replicons. Recombination between insertion sequences (IS) was inferred as the mechanism behind deletion formation; the most numerous IS group was specific for the ST32 clone and has undergone transposition burst since its divergence. Genes functionally related to transition metal metabolism were identified as hotspots for deletions and IS insertions. This functional category was also represented among genes where nonsynonymous point mutations and indels occurred parallelly among patients. Another category exhibiting parallel mutations was oxidative stress protection; mutations in catalase KatG resulted in impaired detoxification of hydrogen peroxide. Deep sequencing revealed substantial polymorphism in genes of both categories within the sputum B. cenocepacia ST32 populations, indicating extensive adaptive evolution. Neither oxidative stress response nor transition metal metabolism genes were previously reported to undergo parallel evolution during chronic CF infection. Mutations in katG and copper metabolism genes were overrepresented in patients where chronic infection developed into CS. Among professional phagocytes, macrophages use both hydrogen peroxide and copper for their bactericidal activity; our results thus tentatively point to macrophages as suspects in pathogenesis towards the fatal CS.

The large Burkholderia cenocepacia populations which persist in cystic fibrosis lungs during many years of chronic infections have an inherent potential for adaptive evolution. The results provided by comparative genomics are key in understanding the processes involved. Mutational events which have taken place allow us to deductively reconstruct the history of chronic infection and to identify driving forces acting upon the bacteria. Beyond the conventional point mutation analysis of next generation sequencing data, we observed interesting phenomena such as large deletions and transposable element movement which represent another facet of adaptive evolution of B. cenocepacia during chronic infection. We also found, unexpectedly, that adaptive evolution in B. cenocepacia strain ST32 affects a set of genes conspicuously different from related species B. dolosa; these appear to be linked to host immune response. Our study provides clues to the complex puzzle of chronic B. cenocepacia infection establishment, persistence and outcome in cystic fibrosis.

Using bioinformatic and phylogenetic approaches to classify transposable elements and understand their complex evolutionary histories

In recent years, much attention has been paid to comparative genomic studies of transposable elements (TEs) and the ensuing problems of their identification, classification, and annotation. Different approaches and diverse automated pipelines are being used to catalogue and categorize mobile genetic elements in the ever-increasing number of prokaryotic and eukaryotic genomes, with little or no connectivity between different domains of life. Here, an overview of the current picture of TE classification and evolutionary relationships is presented, updating the diversity of TE types uncovered in sequenced genomes. A tripartite TE classification scheme is proposed to account for their replicative, integrative, and structural components, and the need to expand in vitro and in vivo studies of their structural and biological properties is emphasized. Bioinformatic studies have now become front and center of novel TE discovery, and experimental pursuits of these discoveries hold great promise for both basic and applied science.

CRISPR-Cas Systems Features and the Gene-Reservoir Role of Coagulase-Negative Staphylococci

The claimed role of gene reservoir of coagulase-negative staphylococci (CoNS) could be contradicted by estimates that CRISPR/Cas systems are found in the genomes of 40–50% of bacteria, as these systems interfere with plasmid uptake in staphylococci. To further correlate this role with presence of CRISPR, we analyzed, by computational methods, 122 genomes from 15 species of CoNS. Only 15% of them harbored CRISPR/Cas systems, and this proportion was much lower for S. epidermidis and S. haemolyticus, the CoNS most frequently associated with opportunistic infections in humans. These systems are of type II or III, and at least two of them are located within SCCmec, a mobile genetic element of Staphylococcus bacterial species. An analysis of the spacers of these CRISPRs, which come from exogenous origin, allowed us to track the transference of the SCCmec, which was exchanged between different strains, species and hosts. Some of the spacers are derived from plasmids described in Staphylococcus species that are different from those in which the CRISPR are found, evidencing the attempt (and failure) of plasmid transference between them. Based on the polymorphisms of the cas1 gene in CRISPRs of types II and III, we developed a multiplex polymerase chain reaction (PCR) suitable to screen and type CRISPR systems in CoNS. The PCR was tested in 59 S. haemolyticus strains, of which only two contained a type III cas1. This gene was shown to be expressed in the exponential growth, stationary phase and during biofilm formation. The low abundance of CRISPRs in CoNS is in accordance with their role as gene reservoirs, but when present, their spacers sequence evidence and give an insight on the dynamics of horizontal genetic transfer among staphylococci.

Genomic insights into the pathogenicity and environmental adaptability of Enterococcus hirae R17 isolated from pork offered for retail sale

Genetic information about Enterococcus hirae is limited, a feature that has compromised our understanding of these clinically challenging bacteria. In this study, comparative analysis was performed of E. hirae R17, a daptomycin‐resistant strain isolated from pork purchased from a retail market in Beijing, China, and three other enterococcal genomes (Enterococcus faecium DO, Enterococcus faecalis V583, and E. hirae ATCC^™9790). Some 1,412 genes were identified that represented the core genome together with an additional 139 genes that were specific to E. hirae R17. The functions of these R17 strain‐specific coding sequences relate to the COGs categories of carbohydrate transport and metabolism and transcription, a finding that suggests the carbohydrate utilization capacity of E. hirae R17 may be more extensive when compared with the other three bacterial species (spp.). Analysis of genomic islands and virulence genes highlighted the potential that horizontal gene transfer played as a contributor of variations in pathogenicity in this isolate. Drug‐resistance gene prediction and antibiotic susceptibility testing indicated E. hirae R17 was resistant to several antimicrobial compounds, including bacitracin, ciprofloxacin, daptomycin, erythromycin, and tetracycline, thereby limiting chemotherapeutic treatment options. Further, tolerance to biocides and metals may confer a phenotype that facilitates the survival and adaptation of this isolate against food preservatives, disinfectants, and antibacterial coatings. The genomic plasticity, mediated by IS elements, transposases, and tandem repeats, identified in the E. hirae R17 genome may support adaptation to new environmental niches, such as those that are found in hospitalized patients. A predicted transmissible plasmid, pRZ1, was found to carry several antimicrobial determinants, along with some predicted pathogenic genes. These data supported the previously determined phenotype confirming that the foodborne E. hirae R17 is a multidrug‐resistant pathogenic bacterium with evident genome plasticity and environmental adaptability.

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.

Annotation of plasmid genes

Good annotation of plasmid genomes is essential to maximise the value of the rapidly increasing volume of plasmid sequences. This short review highlights some of the current issues and suggests some ways forward. Where a well-studied related plasmid system exists we recommend that new annotation adheres to the convention already established for that system, so long as it is based on sound principles and solid experimental evidence, even if some of the new genes are more similar to homologues in different systems. Where a well-established model does not exist we provide generic gene names that reflect likely biochemical activity rather than overall purpose particularly, for example, where genes clearly belong to a type IV secretion system but it is not known whether they function in conjugative transfer or virulence. We also recommend that annotators use a whole system naming approach to avoid ending up with an illogical mixture of names from other systems based on the highest scoring match from a BLAST search. In addition, where function has not been experimentally established we recommend using just the locus tag, rather than a function-related gene name, while recording possible functions as notes rather than in a provisional name.

Bioinformatics tools and databases for whole genome sequence analysis of Mycobacterium tuberculosis

Tuberculosis (TB) is an infectious disease of global public health importance caused by Mycobacterium tuberculosis complex (MTC) in which M. tuberculosis (Mtb) is the major causative agent. Recent advancements in genomic technologies such as next generation sequencing have enabled high throughput cost-effective generation of whole genome sequence information from Mtb clinical isolates, providing new insights into the evolution, genomic diversity and transmission of the Mtb bacteria, including molecular mechanisms of antibiotic resistance. The large volume of sequencing data generated however necessitated effective and efficient management, storage, analysis and visualization of the data and results through development of novel and customized bioinformatics software tools and databases. In this review, we aim to provide a comprehensive survey of the current freely available bioinformatics software tools and publicly accessible databases for genomic analysis of Mtb for identifying disease transmission in molecular epidemiology and in rapid determination of the antibiotic profiles of clinical isolates for prompt and optimal patient treatment.

A survey of transposable element classification systems--a call for a fundamental update to meet the challenge of their diversity and complexity

The increase of publicly available sequencing data has allowed for rapid progress in our understanding of genome composition. As new information becomes available we should constantly be updating and reanalyzing existing and newly acquired data. In this report we focus on transposable elements (TEs) which make up a significant portion of nearly all sequenced genomes. Our ability to accurately identify and classify these sequences is critical to understanding their impact on host genomes. At the same time, as we demonstrate in this report, problems with existing classification schemes have led to significant misunderstandings of the evolution of both TE sequences and their host genomes. In a pioneering publication Finnegan (1989) proposed classifying all TE sequences into two classes based on transposition mechanisms and structural features: the retrotransposons (class I) and the DNA transposons (class II). We have retraced how ideas regarding TE classification and annotation in both prokaryotic and eukaryotic scientific communities have changed over time. This has led us to observe that: (1) a number of TEs have convergent structural features and/or transposition mechanisms that have led to misleading conclusions regarding their classification, (2) the evolution of TEs is similar to that of viruses by having several unrelated origins, (3) there might be at least 8 classes and 12 orders of TEs including 10 novel orders. In an effort to address these classification issues we propose: (1) the outline of a universal TE classification, (2) a set of methods and classification rules that could be used by all scientific communities involved in the study of TEs, and (3) a 5-year schedule for the establishment of an International Committee for Taxonomy of Transposable Elements (ICTTE).

VisualTE: a graphical interface for transposable element analysis at the genomic scale

Background

Transposable elements are mobile DNA repeat sequences, known to have high impact on genes, genome structure and evolution. This has stimulated broad interest in the detailed biological studies of transposable elements. Hence, we have developed an easy-to-use tool for the comparative analysis of the structural organization and functional relationships of transposable elements, to help understand their functional role in genomes.

Results

We named our new software VisualTE and describe it here. VisualTE is a JAVA stand-alone graphical interface that allows users to visualize and analyze all occurrences of transposable element families in annotated genomes. VisualTE reads and extracts transposable elements and genomic information from annotation and repeat data. Result analyses are displayed in several graphical panels that include location and distribution on the chromosome, the occurrence of transposable elements in the genome, their size distribution, and neighboring genes’ features and ontologies. With these hallmarks, VisualTE provides a convenient tool for studying transposable element copies and their functional relationships with genes, at the whole-genome scale, and in diverse organisms.

Conclusions

VisualTE graphical interface makes possible comparative analyses of transposable elements in any annotated sequence as well as structural organization and functional relationships between transposable elements and other genetic object. This tool is freely available at: http://lcb.cnrs-mrs.fr/spip.php?article867.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1351-5) contains supplementary material, which is available to authorized users.

A TALE of transposition: Tn3-like transposons play a major role in the spread of pathogenicity determinants of Xanthomonas citri and other xanthomonads

Members of the genus Xanthomonas are among the most important phytopathogens. A key feature of Xanthomonas pathogenesis is the translocation of type III secretion system (T3SS) effector proteins (T3SEs) into the plant target cells via a T3SS. Several T3SEs and a murein lytic transglycosylase gene (mlt, required for citrus canker symptoms) are found associated with three transposition-related genes in Xanthomonas citri plasmid pXAC64. These are flanked by short inverted repeats (IRs). The region was identified as a transposon, TnXax1, with typical Tn3 family features, including a transposase and two recombination genes. Two 14-bp palindromic sequences within a 193-bp potential resolution site occur between the recombination genes. Additional derivatives carrying different T3SEs and other passenger genes occur in different Xanthomonas species. The T3SEs include transcription activator-like effectors (TALEs). Certain TALEs are flanked by the same IRs as found in TnXax1 to form mobile insertion cassettes (MICs), suggesting that they may be transmitted horizontally. A significant number of MICs carrying other passenger genes (including a number of TALE genes) were also identified, flanked by the same TnXax1 IRs and delimited by 5-bp target site duplications. We conclude that a large fraction of T3SEs, including individual TALEs and potential pathogenicity determinants, have spread by transposition and that TnXax1, which exhibits all of the essential characteristics of a functional transposon, may be involved in driving MIC transposition. We also propose that TALE genes may diversify by fork slippage during the replicative Tn3 family transposition. These mechanisms may play a crucial role in the emergence of Xanthomonas pathogenicity.

Xanthomonas genomes carry many insertion sequences (IS) and transposons, which play an important role in their evolution and architecture. This study reveals a key relationship between transposons and pathogenicity determinants in Xanthomonas. We propose that several transposition events mediated by a Tn3-like element carrying different sets of passenger genes, such as different type III secretion system effectors (including transcription activation-like effectors [TALEs]), were determinant in the evolution and emergence of Xanthomonas pathogenicity. TALE genes are DNA-binding effectors that modulate plant transcription. We also present a model for generating TALE gene diversity based on fork slippage associated with the replicative transposition mechanism of Tn3-like transposons. This may provide a mechanism for niche adaptation, specialization, host-switching, and other lifestyle changes. These results will also certainly lead to novel insights into the evolution and emergence of the various diseases caused by different Xanthomonas species and pathovars.

Survey of chimeric IStron elements in bacterial genomes: multiple molecular symbioses between group I intron ribozymes and DNA transposons

IStrons are chimeric genetic elements composed of a group I intron associated with an insertion sequence (IS). The group I intron is a catalytic RNA providing the IStron with self-splicing ability, which renders IStron insertions harmless to the host genome. The IS element is a DNA transposon conferring mobility, and thus allowing the IStron to spread in genomes. IStrons are therefore a striking example of a molecular symbiosis between unrelated genetic elements endowed with different functions. In this study, we have conducted the first comprehensive survey of IStrons in sequenced genomes that provides insights into the distribution, diversity, origin and evolution of IStrons. We show that IStrons have a restricted phylogenetic distribution limited to two bacterial phyla, the Firmicutes and the Fusobacteria. Nevertheless, diverse IStrons representing two major groups targeting different insertion site motifs were identified. This taken with the finding that while the intron components of all IStrons belong to the same structural class, they are fused to different IS families, indicates that multiple intron–IS symbioses have occurred during evolution. In addition, introns and IS elements related to those that were at the origin of IStrons were also identified.