Evolution and dynamics of megaplasmids with genome sizes larger than 100 kb in the Bacillus cereus group

Genome Sequence Analysis of Native Xenorhabdus Strains Isolated from Entomopathogenic Nematodes in Argentina

Entomopathogenic nematodes from the genus Steinernema (Nematoda: Steinernematidae) are capable of causing the rapid killing of insect hosts, facilitated by their association with symbiotic Gram-negative bacteria in the genus Xenorhabdus (Enterobacterales: Morganellaceae), positioning them as interesting candidate tools for the control of insect pests. In spite of this, only a limited number of species from this bacterial genus have been identified from their nematode hosts and their insecticidal properties documented. This study aimed to perform the genome sequence analysis of fourteen Xenorhabdus strains that were isolated from Steinernema nematodes in Argentina. All of the strains were found to be able of killing 7th instar larvae of Galleria mellonella (L.) (Lepidoptera: Pyralidae). Their sequenced genomes harbour 110 putative insecticidal proteins including Tc, Txp, Mcf, Pra/Prb and App homologs, plus other virulence factors such as putative nematocidal proteins, chitinases and secondary metabolite gene clusters for the synthesis of different bioactive compounds. Maximum-likelihood phylogenetic analysis plus average nucleotide identity calculations strongly suggested that three strains should be considered novel species. The species name for strains PSL and Reich (same species according to % ANI) is proposed as Xenorhabdus littoralis sp. nov., whereas strain 12 is proposed as Xenorhabdus santafensis sp. nov. In this work, we present a dual insight into the biocidal potential and diversity of the Xenorhabdus genus, demonstrated by different numbers of putative insecticidal genes and biosynthetic gene clusters, along with a fresh exploration of the species within this genus.

Get to Know Your Neighbors: Characterization of Close Bacillus anthracis Isolates and Toxin Profile Diversity in the Bacillus cereus Group

Unexpected atypical isolates of Bacillus cereus s.l. occasionally challenge conventional microbiology and even the most advanced techniques for anthrax detection. For anticipating and gaining trust, 65 isolates of Bacillus cereus s.l. of diverse origin were sequenced and characterized. The BTyper3 tool was used for assignation to genomospecies B. mosaicus (34), B. cereus s.s (29) and B. toyonensis (2), as well as virulence factors and toxin profiling. None of them carried any capsule or anthrax-toxin genes. All harbored the non-hemolytic toxin nheABC and sphygomyelinase spH genes, whereas 41 (63%), 30 (46%), 11 (17%) and 6 (9%) isolates harbored cytK-2, hblABCD, cesABCD and at least one insecticidal toxin gene, respectively. Matrix-assisted laser desorption ionization-time of flight mass spectrometry confirmed the production of cereulide (ces genes). Phylogeny inferred from single-nucleotide polymorphisms positioned isolates relative to the B. anthracis lineage. One isolate (BC38B) was of particular interest as it appeared to be the closest B. anthracis neighbor described so far. It harbored a large plasmid similar to other previously described B. cereus s.l. megaplasmids and at a lower extent to pXO1. Whereas bacterial collection is enriched, these high-quality public genetic data offer additional knowledge for better risk assessment using future NGS-based technologies of detection.

DNA Nanomachine (DNM) Biplex Assay for Differentiating Bacillus cereus Species

Conventional methods for the detection and differentiation of Bacillus cereus group species have drawbacks mostly due to the complexity of genetic discrimination between the Bacillus cereus species. Here, we describe a simple and straightforward assay based on the detected unamplified bacterial 16S rRNA by DNA nanomachine (DNM). The assay uses a universal fluorescent reporter and four all-DNA binding fragments, three of which are responsible for "opening up" the folded rRNA while the fourth stand is responsible for detecting single nucleotide variation (SNV) with high selectivity. Binding of the DNM to 16S rRNA results in the formation of the 10-23 deoxyribozyme catalytic core that cleaves the fluorescent reporter and produces a signal, which is amplified over time due to catalytic turnover. This developed biplex assay enables the detection of B. thuringiensis 16S rRNA at fluorescein and B. mycoides at Cy5 channels with a limit of detection of 30 × 10³ and 35 × 10³ CFU/mL, respectively, after 1.5 h with a hands-on time of ~10 min. The new assay may simplify the analysis of biological RNA samples and might be useful for environmental monitoring as a simple and inexpensive alternative to amplification-based nucleic acid analysis. The DNM proposed here may become an advantageous tool for detecting SNV in clinically significant DNA or RNA samples and can easily differentiate SNV under broadly variable experimental conditions and without prior amplification.

Comparative genomics of Bacillus cereus sensu lato spp. biocontrol strains in correlation to in-vitro phenotypes and plant pathogen antagonistic capacity

Bacillus cereus sensu lato (Bcsl) strains are widely explored due to their capacity to antagonize a broad range of plant pathogens. These include B. cereus sp. UW85, whose antagonistic capacity is attributed to the secondary metabolite Zwittermicin A (ZwA). We recently isolated four soil and root-associated Bcsl strains (MO2, S-10, S-25, LSTW-24) that displayed different growth profiles and in-vitro antagonistic effects against three soilborne plant pathogens models: Pythium aphanidermatum (oomycete) Rhizoctonia solani (basidiomycete), and Fusarium oxysporum (ascomycete). To identify genetic mechanisms potentially responsible for the differences in growth and antagonistic phenotypes of these Bcsl strains, we sequenced and compared their genomes, and that of strain UW85 using a hybrid sequencing pipeline. Despite similarities, specific Bcsl strains had unique secondary metabolite and chitinase-encoding genes that could potentially explain observed differences in in-vitro chitinolytic potential and anti-fungal activity. Strains UW85, S-10 and S-25 contained a (~500 Kbp) mega-plasmid that harbored the ZwA biosynthetic gene cluster. The UW85 mega-plasmid contained more ABC transporters than the other two strains, whereas the S-25 mega-plasmid carried a unique cluster containing cellulose and chitin degrading genes. Collectively, comparative genomics revealed several mechanisms that can potentially explain differences in in-vitro antagonism of Bcsl strains toward fungal plant pathogens.

Evidence of structural rearrangements in ESBL-positive pESI(like) megaplasmids of S.Infantis

The increasing prevalence of pESI(like)-positive, multidrug-resistant (MDR) S. Infantis in Europe is a cause of major concern. As previously demonstrated, the pESI(like) megaplasmid is not only a carrier of antimicrobial resistant (AMR) genes (at least tet, dfr, and sul genes), but also harbours several virulence and fitness genes, and toxin/antitoxin systems that enhance its persistence in the S. Infantis host. In this study, five prototype pESI(like) plasmids, of either CTX-M-1 or CTX-M-65 ESBL-producing strains, were long-read sequenced using Oxford Nanopore Technology (ONT), and their complete sequences were resolved. Comparison of the structure and gene content of the five sequenced plasmids, and further comparison with previously published pESI(like) sequences, indicated that although the sequence of such pESI(like) 'mosaic' plasmids remains almost identical, their structures appear different and composed of regions inserted or transposed after different events. The results obtained in this study are essential to better understand the plasticity and the evolution of the pESI(like) megaplasmid, and therefore to better address risk management options and policy decisions to fight against AMR and MDR in Salmonella and other food-borne pathogens. Graphical representation of the pESI-like plasmid complete sequence (ID 12037823/11). Block colours indicate the function of the genes: red: repB gene; pink: class I integrons (IntI); yellow; mobile elements; blue: resistance genes; green: toxin/anti-toxin systems; grey: mer operon; light green: genes involve in conjugation.

The hazard of carbapenemase (OXA-181)-producing Escherichia coli spreading in pig and veal calf holdings in Italy in the genomics era: Risk of spill over and spill back between humans and animals

Carbapenemase-producing Enterobacterales (CPE) are considered a major public health issue. In the frame of the EU Harmonized AMR Monitoring program conducted in Italy in 2021, 21 epidemiological units of fattening pigs (6.98%; 95% CI 4.37-10.47%; 21/301) and four epidemiological units of bovines <12 months (1.29%; 95% CI 0.35-3.27%, 4/310) resulted positive to OXA-48-like-producing E. coli (n = 24 OXA-181, n = 1 OXA-48). Whole Genome Sequencing (WGS) for in-depth characterization, genomics and cluster analysis of OXA-181-(and one OXA-48) producing E. coli isolated, was performed. Tracing-back activities at: (a) the fattening holding of origin of one positive slaughter batch, (b) the breeding holding, and (c) one epidemiologically related dairy cattle holding, allowed detection of OXA-48-like-producing E. coli in different units and comparison of further human isolates from fecal samples of farm workers. The OXA-181-producing isolates were multidrug resistant (MDR), belonged to different Sequence Types (STs), harbored the IncX and IncF plasmid replicons and multiple virulence genes. Bioinformatics analysis of combined Oxford Nanopore Technologies (ONT) long reads and Illumina short reads identified bla OXA-181 as part of a transposon in IncX1, IncX3, and IncFII fully resolved plasmids from 16 selected E. coli, mostly belonging to ST5229, isolated during the survey at slaughter and tracing-back activities. Although human source could be the most likely cause for the introduction of the bla OXA-181-carrying IncX1 plasmid in the breeding holding, concerns arise from carbapenemase OXA-48-like-producing E. coli spreading in 2021 in Italian fattening pigs and, to a lesser extent, in veal calf holdings.

Conjugation-mediated transfer of pXO16, a large plasmid from Bacillus thuringiensis sv. israelensis, across the Bacillus cereus group and its impact on host phenotype

pXO16, the 350 kb-conjugative plasmid from Bacillus thuringiensis sv. israelensis promotes its own transfer at high efficiency, triggers the transfer of mobilizable and non-mobilizable plasmids, as well as the transfer of host chromosomal loci. Naturally found in B. thuringiensis sv. israelensis, pXO16 transfers to various strains of Bacillus cereus sensu lato (s.l.) at a wide range of frequencies. Despite this host diversity, a paradox remains between the relatively large host spectrum and the natural occurrence of pXO16, so far restricted to B. thuringiensis sv. israelensis. Proposing first insights exploring this paradox, we investigated the behaviour of pXO16 amongst different members of the B. cereus group. We first looked at the transfer of pXO16 to two new host clusters of B. cereus s.l., Bacillus mycoides and Bacillus anthracis clusters. This examination brought to light the impairment of the characteristic rhizoidal phenotype of B. mycoides in presence of pXO16. We also explored the stability of pXO16 at different temperatures as some B. cereus group members are well-known for their psychro- or thermo-tolerance. This shed light on the thermo-sensitivity of the plasmid. The influence of pXO16 on its host cell growth and on swimming capacity also revealed no or limited impact on its natural host B. thuringiensis sv. israelensis. On the contrary, pXO16 affected more strongly both the growth and swimming capacity of other B. cereus s.l. hosts. This reinforced the running hypothesis of a co-evolution between pXO16 and B. thuringiensis sv. israelensis, enabling the plasmid maintenance without impairing the host strain development.

Genomic attributes of thermophilic and hyperthermophilic bacteria and archaea

Thermophiles and hyperthermophiles are immensely useful in understanding the evolution of life, besides their utility in environmental and industrial biotechnology. Advancements in sequencing technologies have revolutionized the field of microbial genomics. The massive generation of data enhances the sequencing coverage multi-fold and allows to analyse the entire genomic features of microbes efficiently and accurately. The mandate of a pure isolate can also be bypassed where whole metagenome-assembled genomes and single cell-based sequencing have fulfilled the majority of the criteria to decode various attributes of microbial genomes. A boom has, therefore, been seen in analysing the extremophilic bacteria and archaea using sequence-based approaches. Due to extensive sequence analysis, it becomes easier to understand the gene flow and their evolution among the members of bacteria and archaea. For instance, sequencing unveiled that Thermotoga maritima shares around 24% of genes of archaeal origin. Comparative and functional genomics provide an analytical view to understanding the microbial diversity of thermophilic bacteria and archaea, their interactions with other microbes, their adaptations, gene flow, and evolution over time. In this review, the genomic features of thermophilic bacteria and archaea are dealt with comprehensively.

Identification of a New Pathogenicity Island Within the Large pAH187_270 Plasmid Involved in Bacillus cereus Virulence

Objectives

Bacillus cereus is responsible for food poisoning and rare but severe clinical infections. The pathogenicity of B. cereus strains varies from harmless to lethal strains. The objective of this study was to characterize three B. cereus isolates isolated from the same patient and identify their virulence potentials.

Methods

Three isolates of B. cereus were isolated from various blood samples from a patient who developed sepsis following a central venous catheter infection. The three isolates were compared by WGS, genotyping and SNP analysis. Furthermore, the isolates were compared by phenotypical analysis including bacterial growth, morphology, germination efficacy, toxin production, antibiotic susceptibility and virulence in an insect model of infection.

Results

According to WGS and genotyping, the 3 isolates were shown to be identical strains. However, the last recovered strain had lost the mega pAH187_270 plasmid. This last strain showed different phenotypes compared to the first isolated strain, such as germination delay, different antibiotic susceptibility and a decreased virulence capacity towards insects. A 50- kbp region of pAH187_270 plasmid was involved in the virulence potential and could thus be defined as a new pathogenicity island of B. cereus.

Conclusions

These new findings help in the understanding of B. cereus pathogenic potential and complexity and provide further hints into the role of large plasmids in the virulence of B. cereus strains. This may provide tools for a better assessment of the risks associated with B. cereus hospital contamination to improve hygiene procedure and patient health.

What makes a megaplasmid?

Naturally occurring plasmids come in different sizes. The smallest are less than a kilobase of DNA, while the largest can be over three orders of magnitude larger. Historically, research has tended to focus on smaller plasmids that are usually easier to isolate, manipulate and sequence, but with improved genome assemblies made possible by long-read sequencing, there is increased appreciation that very large plasmids—known as megaplasmids—are widespread, diverse, complex, and often encode key traits in the biology of their host microorganisms. Why are megaplasmids so big? What other features come with large plasmid size that could affect bacterial ecology and evolution? Are megaplasmids 'just' big plasmids, or do they have distinct characteristics? In this perspective, we reflect on the distribution, diversity, biology, and gene content of megaplasmids, providing an overview to these large, yet often overlooked, mobile genetic elements.

This article is part of the theme issue ‘The secret lives of microbial mobile genetic elements’.

Positive Selection Inhibits Plasmid Coexistence in Bacterial Genomes

Plasmids play an important role in bacterial evolution by transferring niche-adaptive functional genes between lineages, thus driving genomic diversification. Bacterial genomes commonly contain multiple, coexisting plasmid replicons, which could fuel adaptation by increasing the range of gene functions available to selection and allowing their recombination. However, plasmid coexistence is difficult to explain because the acquisition of plasmids typically incurs high fitness costs for the host cell. Here, we show that plasmid coexistence was stably maintained without positive selection for plasmid-borne gene functions and was associated with compensatory evolution to reduce fitness costs. In contrast, with positive selection, plasmid coexistence was unstable despite compensatory evolution. Positive selection discriminated between differential fitness benefits of functionally redundant plasmid replicons, retaining only the more beneficial plasmid. These data suggest that while the efficiency of negative selection against plasmid fitness costs declines over time due to compensatory evolution, positive selection to maximize plasmid-derived fitness benefits remains efficient. Our findings help to explain the forces structuring bacterial genomes: coexistence of multiple plasmids in a genome is likely to require either rare positive selection in nature or nonredundancy of accessory gene functions among the coexisting plasmids.

Genomic and transcriptomic insights into the habitat adaptation of the diazotrophic paddy-field cyanobacterium Nostoc sphaeroides

Nitrogen-fixing cyanobacteria are common in paddy fields, one of the most productive wetland ecosystems. Here, we present the complete genome of Nostoc sphaeroides, a paddy-field diazotroph used for food and medicine for more than 1700 years and deciphered the transcriptional regulation during the developmental transition from hormogonia to vegetative filaments with heterocysts. The genome of N. sphaeroides consists of one circular chromosome (6.48 Mb), one of the largest ever reported megaplasmids (2.34 Mb), and seven plasmids. Multiple gene families involved in the adaption to high solar radiation and water fluctuation conditions were found expanded, while genes involved in anoxic adaptation and phosphonate utilization are located on the megaplasmid, suggesting its indispensable role in environmental adaptation. Distinct gene expression patterns were observed during the light-intensity-regulated transition from hormogonia to vegetative filaments, specifically, genes encoding proteins involved in photosynthetic light reaction, carbon fixation, nitrogen metabolism and heterocyst differentiation were significantly upregulated, whereas genes related to cell motility were down-regulated. Our results provide genomic and transcriptomic insights into the adaptation of a filamentous nitrogen-fixing cyanobacterium to the highly dynamic paddy-field habitat, suggesting N. sphaeroides as an excellent system to understand the transition from aquatic to terrestrial habitats and to support sustainable rice production.

Coevolving Plasmids Drive Gene Flow and Genome Plasticity in Host-Associated Intracellular Bacteria

Plasmids are important in microbial evolution and adaptation to new environments. Yet, carrying a plasmid can be costly, and long-term association of plasmids with their hosts is poorly understood. Here, we provide evidence that the Chlamydiae, a phylum of strictly host-associated intracellular bacteria, have coevolved with their plasmids since their last common ancestor. Current chlamydial plasmids are amalgamations of at least one ancestral plasmid and a bacteriophage. We show that the majority of plasmid genes are also found on chromosomes of extant chlamydiae. The most conserved plasmid gene families are predominantly vertically inherited, while accessory plasmid gene families show significantly increased mobility. We reconstructed the evolutionary history of plasmid gene content of an entire bacterial phylum over a period of around one billion years. Frequent horizontal gene transfer and chromosomal integration events illustrate the pronounced impact of coevolution with these extrachromosomal elements on bacterial genome dynamics in host-dependent microbes.

Coexistence of cry9 with the vip3A Gene in an Identical Plasmid of Bacillus thuringiensis Indicates Their Synergistic Insecticidal Toxicity

Bacillus thuringiensis (Bt) strains may express several insecticidal proteins with synergistic features, achieving high insecticidal toxicity and delaying development of resistance in insect pests. Previous work showed that Cry9Aa and Vip3Aa proteins present synergistic activity against Chilo suppressalis. In this study, genome-wide analysis of 489 Bt genomes revealed that cry9A was associated with the vip3A gene in seven Bt strains. Among all Bt genomes analyzed, not a single strain was found to have the cry9A gene alone without the presence of the vip3A gene. The complete genome sequencing of two Bt strains, 4AP1 and 4AO1, revealed that cry9A and vip3A genes were located in the same plasmid in both strains. The genome context analysis suggested a recombination mechanism responsible for the insertion of the cry9A gene into the plasmid containing vip3A. The coexistence of Cry9A with Vip3A proteins in strain 4AP1 was confirmed by liquid chromatography-tandem mass spectrometry and western blot analyses. Furthermore, another Cry9 protein codified by the gene in the identical plasmid also showed synergistic activity with the Vip3A protein. Overall, our results support that cry9 genes coexisted with vip3A and that complete genome sequencing combined with protein expression analysis may be used to identify associations of insecticidal proteins with potential synergistic toxicity.

The CRISPR-Cas systems were selectively inactivated during evolution of Bacillus cereus group for adaptation to diverse environments

CRISPR-Cas systems are considered as barriers to horizontal gene transfer (HGT). However, the influence of such systems on HGT within species is unclear. Also, little is known about the impact of CRISPR-Cas systems on bacterial evolution at the population level. Here, using Bacillus cereus sensu lato as model, we investigate the interplay between CRISPR-Cas systems and HGT at the population scale. We found that only a small fraction of the strains have CRISPR-Cas systems (13.9% of 1871), and most of such systems are defective based on their gene content analysis. Comparative genomic analysis revealed that the CRISPR-Cas systems are barriers to HGT within this group, since strains harboring active systems contain less mobile genetic elements (MGEs), have lower fraction of unique genes and also display limited environmental distributions than strains without active CRISPR-Cas systems. The introduction of a functional CRISPR-Cas system into a strain lacking the system resulted in reduced adaptability to various stresses and decreased pathogenicity of the transformant strain, indicating that B. cereus group strains could benefit from inactivating such systems. Our work provides a large-scale case to support that the CRISPR-Cas systems are barriers to HGT within species, and that in the B. cereus group the inactivation of CRISPR-Cas systems correlated with acquisition of MGEs that could result in better adaptation to diverse environments.

Genetic characterisation of variants of the virulence plasmid, pSLT, in Salmonella enterica serovar Typhimurium provides evidence of a variety of evolutionary directions consistent with vertical rather than horizontal transmission

The virulence plasmid pSLT as exemplified by the 94 Kb plasmid in Salmonella Typhimurium strain LT2 is only found in isolates of serovar Typhimurium. While it occurs commonly among such isolates recent genotyping methods have shown that it is mostly confined to certain genotypes. Although pSLT plasmids are capable of self-transmissibility under experimental conditions their confinement to certain host genotypes suggests that in practice they are maintained by vertical rather than by horizontal transmission. This would imply that evolution of the pSLT plasmid proceeds in parallel with evolution of its host. The development of a phylogenetic evolutionary framework for genotypes of S. Typhimurium based on single-nucleotide-polymorphism (SNPs) typing provided an opportunity to test whether the pSLT plasmid coevolves with its host genotype. Accordingly SNPs analysis was applied to the pSLT plasmids from 71 strains S. Typhimurium of Australian and international origins representing most of the genotypes which commonly have a pSLT. The phylogenetic tree showed that pSLT sequences clustered into almost the same groups as the host chromosomes so that each pSLT genotype was associated with a single host genotype. A search for tandem repeats in pSLT plasmids showed that a 9 bp VNTR in the traD gene occurred in the pSLT from all isolates belonging to Clade II but not from isolates belonging to Clade I. Another 9 bp repeat occurred only in three Clade I genotypes with a recent common ancestor. The evidence relating to both of these VNTRs supports the proposition that the pSLT plasmid is only transmitted vertically. Some isolates belonging to one S. Typhimurium genotype were found to have pSLTs which have lost a large block of genes when a resistance gene cassette has been acquired. Examples were found of pSLT plasmids which have recombined with other plasmids to form fusion plasmids sometimes with loss of some pSLT genes. In all cases the underlying genotype of the modified pSLT was the same as the genotype of regular pSLTs with the same host genotype implying that these changes have occurred within the host cell of the pSLT plasmid.

Complete Sequences of Two Plasmids Found in a Brazilian Bacillus thuringiensis Serovar israelensis Strain

Plasmids play a crucial role in the evolution of bacterial genomes by mediating horizontal gene transfer. In this work, we sequenced two plasmids found in a Brazilian Bacillus thuringiensis serovar israelensis strain which showed 100% nucleotide identities with Bacillus thuringiensis serovar kurstaki plasmids.

Plasmids play a crucial role in the evolution of bacterial genomes by mediating horizontal gene transfer. In this work, we sequenced two plasmids found in a Brazilian Bacillus thuringiensis serovar israelensis strain which showed 100% nucleotide identities with Bacillus thuringiensis serovar kurstaki plasmids.

Antimicrobial Resistance Profile of mcr-1 Positive Clinical Isolates of Escherichia coli in China From 2013 to 2016

Multidrug-resistant (MDR) Escherichia coli poses a great challenge for public health in recent decades. Polymyxins have been reconsidered as a valuable therapeutic option for the treatment of infections caused by MDR E. coli. A plasmid-encoded colistin resistance gene mcr-1 encoding phosphoethanolamine transferase has been recently described in Enterobacteriaceae. In this study, a total of 123 E. coli isolates obtained from patients with diarrheal diseases in China were used for the genetic analysis of colistin resistance in clinical isolates. Antimicrobial resistance profile of polymyxin B (PB) and 11 commonly used antimicrobial agents were determined. Among the 123 E. coli isolates, 9 isolates (7.3%) were resistant to PB and PCR screening showed that seven (5.7%) isolates carried the mcr-1 gene. A hybrid sequencing analysis using single-molecule, real-time (SMRT) sequencing and Illumina sequencing was then performed to resolve the genomes of the seven mcr-1 positive isolates. These seven isolates harbored multiple plasmids and are MDR, with six isolates carrying one mcr-1 positive plasmid and one isolate (14EC033) carrying two mcr-1 positive plasmids. These eight mcr-1 positive plasmids belonged to the IncX4, IncI2, and IncP1 types. In addition, the mcr-1 gene was the solo antibiotic resistance gene identified in the mcr-1 positive plasmids, while the rest of the antibiotic resistance genes were mostly clustered into one or two plasmids. Interestingly, one mcr-1 positive isolate (14EC047) was susceptible to PB, and we showed that the activity of MCR-1-mediated colistin resistance was not phenotypically expressed in 14EC047 host strain. Furthermore, three isolates exhibited resistance to PB but did not carry previously reported mcr-related genes. Multilocus sequence typing (MLST) showed that these mcr-1 positive E. coli isolates belonged to five different STs, and three isolates belonged to ST301 which carried multiple virulence factors related to diarrhea. Additionally, the mcr-1 positive isolates were all susceptible to imipenem (IMP), suggesting that IMP could be used to treat infection caused by mcr-1 positive E. coli isolates. Collectively, this study showed a high occurrence of mcr-1 positive plasmids in patients with diarrheal diseases of Guangzhou in China and the abolishment of the MCR-1 mediated colistin resistance in one E. coli isolate.

Bacillus wiedmannii biovar thuringiensis: A Specialized Mosquitocidal Pathogen with Plasmids from Diverse Origins

Bacillus cereus sensu lato also known as B. cereus group is composed of an ecologically diverse bacterial group with an increasing number of related species, some of which are medically or agriculturally important. Numerous efforts have been undertaken to allow presumptive differentiation of B. cereus group species from one another. FCC41 is a Bacillus sp. strain toxic against mosquito species like Aedes aegypti, Aedes (Ochlerotatus) albifasciatus, Culex pipiens, Culex quinquefasciatus, and Culex apicinus, some of them responsible for the transmission of vector-borne diseases. Here, we report the complete genome sequence of FCC41 strain, which consists of one circular chromosome and eight circular plasmids ranging in size from 8 to 490 kb. This strain harbors six crystal protein genes, including cry24Ca, two cry4-like and two cry52-like, a cry41-like parasporin gene and multiple virulence factors. The phylogenetic analysis of the whole-genome sequence of this strain with molecular approaches places this strain into the Bacillus wiedmannii cluster. However, according with phenotypical characteristics such as the mosquitocidal activity due to the presence of Cry proteins found in the parasporal body and cry genes encoded in plasmids of different sizes, indicate that this strain could be renamed as B. wiedmannii biovar thuringiensis strain FCC41.

Genetic Environment of cry1 Genes Indicates Their Common Origin

Although in Bacillus thuringiensis the cry genes coding for the insecticidal crystal proteins are plasmid-borne and are usually associated with mobile genetic elements, several aspects related to their genomic organization, diversification, and transmission remain to be elucidated. Plasmids of B. thuringiensis and other members of the Bacillus cereus group (n = 364) deposited in GenBank were screened for the presence of cry1 genes, and their genetic environment was analyzed using a comparative bioinformatic approach. The cry1 genes were identified in 27 B. thuringiensis plasmids ranging from 64 to 761 kb, and were predominantly associated with the ori44, ori60, or double orf156/orf157 and pXO1-16/pXO1-14 replication systems. In general, the cry1 genes occur individually or as a part of an insecticidal pathogenicity island (PAI), and are preceded by genes coding for an N-acetylmuramoyl-l-alanine amidase and a putative K⁺(Na⁺)/H⁺ antiporter. However, except in the case of the PAI, the latter gene is disrupted by the insertion of IS231B. Similarly, numerous mobile elements were recognized in the region downstream of cry1, except for cry1I that follows cry1A in the PAI. Therefore, the cassette involving cry1 and these two genes, flanked by transposable elements, named as the cry1 cassette, was the smallest cry1-carrying genetic unit recognized in the plasmids. Conservation of the genomic environment of the cry1 genes carried by various plasmids strongly suggests a common origin, possibly from an insecticidal PAI carried by B. thuringiensis megaplasmids.

In Vivo Transfer and Microevolution of Avian Native IncA/C2blaNDM-1-Carrying Plasmid pRH-1238 during a Broiler Chicken Infection Study

The emergence and spread of carbapenemase-producing Enterobacteriaceae (CPE) in wildlife and livestock animals pose an important safety concern for public health. With our in vivo broiler chicken infection study, we investigated the transfer and experimental microevolution of the blaNDM-1-carrying IncA/C2 plasmid (pRH-1238) introduced by avian native Salmonella enterica subsp. enterica serovar Corvallis without inducing antibiotic selection pressure. We evaluated the dependency of the time point of inoculation on donor (S Corvallis [12-SA01738]) and plasmid-free Salmonella recipient [d-tartrate-fermenting (d-Ta+) S Paratyphi B (13-SA01617), referred to here as S Paratyphi B (d-Ta+)] excretion by quantifying their excretion dynamics. Using plasmid profiling by S1 nuclease-restricted pulsed-field gel electrophoresis, we gained insight into the variability of the native plasmid content among S Corvallis reisolates as well as plasmid acquisition in S Paratyphi B (d-Ta+) and the enterobacterial gut microflora. Whole-genome sequencing enabled us to gain an in-depth insight into the microevolution of plasmid pRH-1238 in S Corvallis and enterobacterial recipient isolates. Our study revealed that the fecal excretion of avian native carbapenemase-producing S Corvallis is significantly higher than that of S Paratyphi (d-Ta+) and is not hampered by S Paratyphi (d-Ta+). Acquisition of pRH-1238 in other Enterobacteriaceae and several events of plasmid pRH-1238 transfer to different Escherichia coli sequence types and Klebsiella pneumoniae demonstrated an interspecies broad host range. Regardless of the microevolutionary structural deletions in pRH-1238, the single carbapenem resistance marker blaNDM-1 was maintained on pRH-1238 throughout the trial. Furthermore, we showed the importance of the gut E. coli population as a vector of pRH-1238. In a potential scenario of the introduction of NDM-1-producing S Corvallis into a broiler flock, the pRH-1238 plasmid could persist and spread to a broad host range even in the absence of antibiotic pressure.

Lineage-specific plasmid acquisition and the evolution of specialized pathogens in Bacillus thuringiensis and the Bacillus cereus group

Bacterial plasmids can vary from small selfish genetic elements to large autonomous replicons that constitute a significant proportion of total cellular DNA. By conferring novel function to the cell, plasmids may facilitate evolution but their mobility may be opposed by co-evolutionary relationships with chromosomes or encouraged via the infectious sharing of genes encoding public goods. Here, we explore these hypotheses through large-scale examination of the association between plasmids and chromosomal DNA in the phenotypically diverse Bacillus cereus group. This complex group is rich in plasmids, many of which encode essential virulence factors (Cry toxins) that are known public goods. We characterized population genomic structure, gene content and plasmid distribution to investigate the role of mobile elements in diversification. We analysed coding sequence within the core and accessory genome of 190 B. cereus group isolates, including 23 novel sequences and genes from 410 reference plasmid genomes. While cry genes were widely distributed, those with invertebrate toxicity were predominantly associated with one sequence cluster (clade 2) and phenotypically defined Bacillus thuringiensis. Cry toxin plasmids in clade 2 showed evidence of recent horizontal transfer and variable gene content, a pattern of plasmid segregation consistent with transfer during infectious cooperation. Nevertheless, comparison between clades suggests that co-evolutionary interactions may drive association between plasmids and chromosomes and limit wider transfer of key virulence traits. Proliferation of successful plasmid and chromosome combinations is a feature of specialized pathogens with characteristic niches (Bacillus anthracis, B. thuringiensis) and has occurred multiple times in the B. cereus group.

A minireplicon of plasmid pBMB26 represents a new typical replicon in the megaplasmids of Bacillus cereus group

A new minireplicon (rep26 minireplicon) from pBMB26, the 188 kb indigenous plasmid related to spore-crystal association (SCA) phenotype in Bacillus thuringiensis strain YBT-020, was characterized. A 12 kb EcoRI fragment, which encoded 10 putative open reading frames (ORFs), was capable of supporting replication when cloned in a replication probe vector. Deletion and frame shift mutation analysis showed that a 4.1 kb region encompassing two putative ORFs (orf21 and orf22) was essential for the plasmid replication in B. thuringiensis. Gene orf21 encoding a 49.8 kDa protein (named Rep26) with a helix-turn-helix motif showed no homology with known replication proteins and gene orf22 encoding a protein of 82.6 kDa showed homology to bacterial PcrA helicase. The replication origin of rep26 minireplicon was proved to be located in the coding region of orf21. Plasmid stability experiments indicated that the recombinant plasmid containing rep26 minireplicon has excellent segregational stability. BLASTP analysis revealed that amino acid sequences of ORF21 and ORF22 were well conserved among Bacillus cereus group strains. The rep26 minireplicon was widely distributed and could be defined as a new typical replicon in the megaplasmids of B. cereus group.

Comparative Genomics of Bacillus thuringiensis Reveals a Path to Specialized Exploitation of Multiple Invertebrate Hosts

Understanding the genetic basis of host shifts is a key genomic question for pathogen and parasite biology. The Bacillus cereus group, which encompasses Bacillus thuringiensis and Bacillus anthracis, contains pathogens that can infect insects, nematodes, and vertebrates. Since the target range of the essential virulence factors (Cry toxins) and many isolates is well known, this group presents a powerful system for investigating how pathogens can diversify and adapt to phylogenetically distant hosts. Specialization to exploit insects occurs at the level of the major clade and is associated with substantial changes in the core genome, and host switching between insect orders has occurred repeatedly within subclades. The transfer of plasmids with linked cry genes may account for much of the adaptation to particular insect orders, and network analysis implies that host specialization has produced strong associations between key toxin genes with similar targets. Analysis of the distribution of plasmid minireplicons shows that plasmids with orf156 and orf157, which carry genes encoding toxins against Lepidoptera or Diptera, were contained only by B. thuringiensis in the specialized insect clade (clade 2), indicating that tight genome/plasmid associations have been important in adaptation to invertebrate hosts. Moreover, the accumulation of multiple virulence factors on transposable elements suggests that cotransfer of diverse virulence factors is advantageous in terms of expanding the insecticidal spectrum, overcoming insect resistance, or through gains in pathogenicity via synergistic interactions between toxins.IMPORTANCE Population genomics have provided many new insights into the formation, evolution, and dynamics of bacterial pathogens of humans and other higher animals, but these pathogens usually have very narrow host ranges. As a pathogen of insects and nematodes, Bacillus thuringiensis, which produces toxins showing toxicity to many orders of insects and other invertebrates, can be used as a model to study the evolution of pathogens with wide host ranges. Phylogenomic analysis revealed that host specialization and switching occur at the level of the major clade and subclade, respectively. A toxin gene co-occurrence network indicates that multiple toxins with similar targets were accumulated by the same cell in the whole species. This accumulation may be one of the strategies that B. thuringiensis has used to fight against host resistance. This kind of formation and evolution of pathogens represents a different path used against multiple invertebrate hosts from that used against higher animals.

Detection of the cryptic prophage-like molecule pBtic235 in Bacillus thuringiensis subsp. israelensis

Bacillus thuringiensis has long been recognized to carry numerous extrachromosomal molecules. Of particular interest are the strains belonging to the B. thuringiensis subsp. israelensis lineage, as they can harbor at least seven extrachromosomal molecules. One of these elements seems to be a cryptic molecule that may have been disregarded in strains considered plasmid-less. Therefore, this work focused on this cryptic molecule, named pBtic235. Using different approaches that included transposition-tagging, large plasmid gel electrophoresis and Southern blotting, conjugation and phage-induction experiments, in combination with bioinformatics analyses, it was found that pBtic235 is a hybrid molecule of 235,425 bp whose genome displays potential plasmid- and phage-like modules. The sequence of pBtic235 has been identified in all sequenced genomes of B. thuringiensis subsp. israelensis strains. Here, the pBtic235 sequence was considered identical to that of plasmid pBTHD789-2 from strain HD-789. Despite the fact that the pBtic235 genome possesses 240 putative CDSs, many of them have no homologs in the databases. However, CDSs coding for potential proteins involved in replication, genome packaging and virion structure, cell lysis, regulation of lytic-lysogenic cycles, metabolite transporters, stress and metal resistance, were identified. The candidate plasmidial prophage pBtic235 exemplifies the notable diversity of the extrachromosomal realm found in B. thuringiensis.

Origin and Evolution of Rickettsial Plasmids

Background

Rickettsia species are strictly intracellular bacteria that have undergone a reductive genomic evolution. Despite their allopatric lifestyle, almost half of the 26 currently validated Rickettsia species have plasmids. In order to study the origin, evolutionary history and putative roles of rickettsial plasmids, we investigated the evolutionary processes that have shaped 20 plasmids belonging to 11 species, using comparative genomics and phylogenetic analysis between rickettsial, microbial and non-microbial genomes.

Results

Plasmids were differentially present among Rickettsia species. The 11 species had 1 to 4 plasmid (s) with a size ranging from 12 kb to 83 kb. We reconstructed pRICO, the last common ancestor of the current rickettsial plasmids. pRICO was vertically inherited mainly from Rickettsia/Orientia chromosomes and diverged vertically into a single or multiple plasmid(s) in each species. These plasmids also underwent a reductive evolution by progressive gene loss, similar to that observed in rickettsial chromosomes, possibly leading to cryptic plasmids or complete plasmid loss. Moreover, rickettsial plasmids exhibited ORFans, recent gene duplications and evidence of horizontal gene transfer events with rickettsial and non-rickettsial genomes mainly from the α/γ-proteobacteria lineages. Genes related to maintenance and plasticity of plasmids, and to adaptation and resistance to stress mostly evolved under vertical and/or horizontal processes. Those involved in nucleotide/carbohydrate transport and metabolism were under the influence of vertical evolution only, whereas genes involved in cell wall/membrane/envelope biogenesis, cycle control, amino acid/lipid/coenzyme and secondary metabolites biosynthesis, transport and metabolism underwent mainly horizontal transfer events.

Conclusion

Rickettsial plasmids had a complex evolution, starting with a vertical inheritance followed by a reductive evolution associated with increased complexity via horizontal gene transfer as well as gene duplication and genesis. The plasmids are plastic and mosaic structures that may play biological roles similar to or distinct from their co-residing chromosomes in an obligate intracellular lifestyle.

Plasmids are vectors for redundant chromosomal genes in the Bacillus cereus group

Background

Prokaryotic plasmids have played significant roles in the evolution of bacterial genomes and have a great impact on the metabolic functions of the host cell. Many bacterial strains contain multiple plasmids, but the relationships between bacterial plasmids and chromosomes are unclear. We focused on plasmids from the Bacillus cereus group because most strains contain several plasmids.

Results

We collected the genome sequences of 104 plasmids and 20 chromosomes from B. cereus group strains, and we studied the relationships between plasmids and chromosomes by focusing on the pan-genomes of these plasmids and chromosomes. In terms of basic features (base composition and codon usage), the genes on plasmids were more similar to the chromosomal variable genes (distributed genes and unique genes) than to the chromosomal core genes. Although all the functional categories of the chromosomal genes were exhibited by the plasmid genes, the proportions of each category differed between these two gene sets. The 598 gene families shared between chromosomes and plasmids displayed a uniform distribution between the two groups. A phylogenetic analysis of the shared genes, including the chromosomal core gene set, indicated that gene exchange events between plasmids and chromosomes occurred frequently during the evolutionary histories of the strains and species in this group. Moreover, the shared genes between plasmids and chromosomes usually had different promoter and terminator sequences, suggesting that they are regulated by different elements at the transcriptional level.

Conclusions

We speculate that for the entire B. cereus group, adaptive genes are preserved on both plasmids and chromosomes; however, in a single cell, homologous genes on plasmids and the chromosome are controlled by different regulators to reduce the burden of maintaining redundant genes.

Electronic supplementary material

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

Plasmid flux in Escherichia coli ST131 sublineages, analyzed by plasmid constellation network (PLACNET), a new method for plasmid reconstruction from whole genome sequences

Bacterial whole genome sequence (WGS) methods are rapidly overtaking classical sequence analysis. Many bacterial sequencing projects focus on mobilome changes, since macroevolutionary events, such as the acquisition or loss of mobile genetic elements, mainly plasmids, play essential roles in adaptive evolution. Existing WGS analysis protocols do not assort contigs between plasmids and the main chromosome, thus hampering full analysis of plasmid sequences. We developed a method (called plasmid constellation networks or PLACNET) that identifies, visualizes and analyzes plasmids in WGS projects by creating a network of contig interactions, thus allowing comprehensive plasmid analysis within WGS datasets. The workflow of the method is based on three types of data: assembly information (including scaffold links and coverage), comparison to reference sequences and plasmid-diagnostic sequence features. The resulting network is pruned by expert analysis, to eliminate confounding data, and implemented in a Cytoscape-based graphic representation. To demonstrate PLACNET sensitivity and efficacy, the plasmidome of the Escherichia coli lineage ST131 was analyzed. ST131 is a globally spread clonal group of extraintestinal pathogenic E. coli (ExPEC), comprising different sublineages with ability to acquire and spread antibiotic resistance and virulence genes via plasmids. Results show that plasmids flux in the evolution of this lineage, which is wide open for plasmid exchange. MOB_F12/IncF plasmids were pervasive, adding just by themselves more than 350 protein families to the ST131 pangenome. Nearly 50% of the most frequent γ–proteobacterial plasmid groups were found to be present in our limited sample of ten analyzed ST131 genomes, which represent the main ST131 sublineages.

Plasmids are difficult to analyze in WGS datasets, due to the fragmented nature of the obtained sequences. We developed a method, called PLACNET, which greatly facilitates this analysis. As an example, we analyzed the plasmidome of E. coli ST131, an ExPEC clonal group involved in human urinary tract infections and septicemia. Relevant variation within this clone (e.g., antibiotic resistance and virulence) is frequently caused by the acquisition and loss of plasmids and other mobile genetic elements. Nevertheless, our knowledge of the ST131 plasmidome is limited to a few antibiotic resistance plasmids and to identification of replicons from known plasmid groups. PLACNET analysis extends the number of sequenced plasmids in ST131, which can be used for comparative genomics, from 11 to 50. The ST131 plasmidome is seemingly huge, encompassing roughly 50% of the main plasmid groups of γ–proteobacteria. MOB_F12/IncF plasmids are apparently the most active players in the dissemination of relevant genetic information.