Plasmid - Chromosome interplay in natural and non-natural hosts: global transcription study of three Bacillus cereus group strains carrying pCER270 plasmid.

The Bacillus cereus group comprises genetically related Gram-positive spore-forming bacteria that colonize a wide range of ecological niches and hosts. Despite their high degree of genome conservation, extrachromosomal genetic material diverges between these species. The discriminating properties of the B. cereus group strains are mainly due to plasmid-borne toxins, reflecting the importance of horizontal gene transfers in bacterial evolution and species definition. To investigate how a newly acquired megaplasmid can impact the transcriptome of its host, we transferred the pCER270 from the emetic B. cereus strains to phylogenetically distant B. cereus group strains. RNA-sequencing experiments allowed us to determine the transcriptional influence of the plasmid on host gene expression and the impact of the host genomic background on the pCER270 gene expression. Our results show a transcriptional cross-regulation between the megaplasmid and the host genome. pCER270 impacted carbohydrate metabolism and sporulation genes expression, with a higher effect in the natural host of the plasmid, suggesting a role of the plasmid in the adaptation of the carrying strain to its environment. In addition, the host genomes also modulated the expression of pCER270 genes. Altogether, these results provide an example of the involvement of megaplasmids in the emergence of new pathogenic strains.

ORFograph: search for novel insecticidal protein genes in genomic and metagenomic assembly graphs.

BACKGROUND: Since the prolonged use of insecticidal proteins has led to toxin resistance, it is important to search for novel insecticidal protein genes (IPGs) that are effective in controlling resistant insect populations. IPGs are usually encoded in the genomes of entomopathogenic bacteria, especially in large plasmids in strains of the ubiquitous soil bacteria, Bacillus thuringiensis (Bt). Since there are often multiple similar IPGs encoded by such plasmids, their assemblies are typically fragmented and many IPGs are scattered through multiple contigs. As a result, existing gene prediction tools (that analyze individual contigs) typically predict partial rather than complete IPGs, making it difficult to conduct downstream IPG engineering efforts in agricultural genomics. METHODS: Although it is difficult to assemble IPGs in a single contig, the structure of the genome assembly graph often provides clues on how to combine multiple contigs into segments encoding a single IPG. RESULTS: We describe ORFograph, a pipeline for predicting IPGs in assembly graphs, benchmark it on (meta)genomic datasets, and discover nearly a hundred novel IPGs. This work shows that graph-aware gene prediction tools enable the discovery of greater diversity of IPGs from (meta)genomes. CONCLUSIONS: We demonstrated that analysis of the assembly graphs reveals novel candidate IPGs. ORFograph identified both already known genes "hidden" in assembly graphs and potential novel IPGs that evaded existing tools for IPG identification. As ORFograph is fast, one could imagine a pipeline that processes many (meta)genomic assembly graphs to identify even more novel IPGs for phenotypic testing than would previously be inaccessible by traditional gene-finding methods. While here we demonstrated the results of ORFograph only for IPGs, the proposed approach can be generalized to any class of genes. Video abstract.

Long inverted repeats around the chromosome replication terminus in the model strain Bacillus thuringiensis serovar israelensis BGSC 4Q7.

Bacillus thuringiensis serovar israelensis is the most widely used natural biopesticide against mosquito larvae worldwide. Its lineage has been actively studied and a plasmid-free strain, B. thuringiensis serovar israelensis BGSC 4Q7 (4Q7), has been produced. Previous sequencing of the genome of this strain has revealed the persistent presence of a 235 kb extrachromosomal element, pBtic235, which has been shown to be an inducible prophage, although three putative chromosomal prophages have been lost. Moreover, a 492 kb region, potentially including the standard replication terminus, has also been deleted in the 4Q7 strain, indicating an absence of essential genes in this area. We reanalysed the genome coverage distribution of reads for the previously sequenced variant strain, and sequenced two independently maintained samples of the 4Q7 strain. A 553 kb area, close to the 492 kb deletion, was found to be duplicated. This duplication presumably restored the equal sizes of the replichores, and a balanced functioning of replication termination. An analysis of genome assembly graphs revealed a transient association of the host chromosome with the pBtic235 element. This association may play a functional role in the replication of the bacterial chromosome, and the termination of this process in particular. The genome-restructuring events detected may modify the genetic status of cytotoxic or haemolytic toxins, potentially influencing strain virulence. Twelve of the single-nucleotide variants identified in 4Q7 were probably due to the procedure used for strain construction or were present in the precursor of this strain. No sequence variants were found in pBtic235, but the distribution of the corresponding 4Q7 reads indicates a significant difference from counterparts in natural B. thuringiensis serovar israelensis strains, suggesting a duplication or over-replication in 4Q7. Thus, the 4Q7 strain is not a pure plasmid-less offshoot, but a highly genetically modified derivative of its natural ancestor. In addition to potentially influencing virulence, genome-restructuring events can modify the replication termination machinery. These findings have potential implications for the conclusions of virulence studies on 4Q7 as a model, but they also raise interesting fundamental questions about the functioning of the Bacillus genome.

Genomic Sequence of the Prophage-Free Lactococcus lactis Strain IL6288.

We report here the genome sequence of IL6288, a prophage-free derivative of Lactococcus lactis subsp. lactis strain IL1403, and confirm precise deletion of all prophages. Several single-nucleotide variations and an extra copy of the IS981 element, apparently having a minor influence on cell physiology, were also detected in the IL6288 genome.

Role of plasmid plasticity and mobile genetic elements in the entomopathogen Bacillus thuringiensis serovar israelensis.

Bacillus thuringiensis is a well-known biopesticide that has been used for more than 80 years. This spore-forming bacterium belongs to the group of Bacillus cereus that also includes, among others, emetic and diarrheic pathotypes of B. cereus, the animal pathogen Bacillus anthracis and the psychrotolerant Bacillus weihenstephanensis. Bacillus thuringiensis is rather unique since it has adapted its lifestyle as an efficient pathogen of specific insect larvae. One of the peculiarities of B. thuringiensis strains is the extent of their extrachromosomal pool, with strains harbouring more than 10 distinct plasmid molecules. Among the numerous serovars of B. thuringiensis, 'israelensis' is certainly emblematic since its host spectrum is apparently restricted to dipteran insects like mosquitoes and black flies, vectors of human and animal diseases such as malaria, yellow fever, or river blindness. In this review, the putative role of the mobile gene pool of B. thuringiensis serovar israelensis in its pathogenicity and dedicated lifestyle is reviewed, with specific emphasis on the nature, diversity, and potential mobility of its constituents. Variations among the few related strains of B. thuringiensis serovar israelensis will also be reported and discussed in the scope of this specialised insect pathogen, whose lifestyle in the environment remains largely unknown.

Molecular Mechanisms That Contribute to Horizontal Transfer of Plasmids by the Bacteriophage SPP1.

Natural transformation and viral-mediated transduction are the main avenues of horizontal gene transfer in Firmicutes. Bacillus subtilis SPP1 is a generalized transducing bacteriophage. Using this lytic phage as a model, we have analyzed how viral replication and recombination systems contribute to the transfer of plasmid-borne antibiotic resistances. Phage SPP1 DNA replication relies on essential phage-encoded replisome organizer (G38P), helicase loader (G39P), hexameric replicative helicase (G40P), recombinase (G35P) and in less extent on the partially dispensable 5'→3' exonuclease (G34.1P), the single-stranded DNA binding protein (G36P) and the Holliday junction resolvase (G44P). Correspondingly, the accumulation of linear concatemeric plasmid DNA, and the formation of transducing particles were blocked in the absence of G35P, G38P, G39P, and G40P, greatly reduced in the G34.1P, G36P mutants, and slightly reduced in G44P mutants. In contrast, establishment of injected linear plasmid DNA in the recipient host was independent of viral-encoded functions. DNA homology between SPP1 and the plasmid, rather than a viral packaging signal, enhanced the accumulation of packagable plasmid DNA. The transfer efficiency was also dependent on plasmid copy number, and rolling-circle plasmids were encapsidated at higher frequencies than theta-type replicating plasmids.

Comparative genomics of extrachromosomal elements in Bacillus thuringiensis subsp. israelensis.

Bacillus thuringiensis subsp. israelensis is one of the most important microorganisms used against mosquitoes. It was intensively studied following its discovery and became a model bacterium of the B. thuringiensis species. Those studies focused on toxin genes, aggregation-associated conjugation, linear genome phages, etc. Recent announcements of genomic sequences of different strains have not been explicitly related to the biological properties studied. We report data on plasmid content analysis of four strains using ultra-high-throughput sequencing. The strains were commercial product isolates, with their putative ancestor and type B. thuringiensis subsp. israelensis strain sequenced earlier. The assembled contigs corresponding to published and novel data were assigned to plasmids described earlier in B. thuringiensis subsp. israelensis and other B. thuringiensis strains. A new 360 kb plasmid was identified, encoding multiple transporters, also found in most of the earlier sequenced strains. Our genomic data show the presence of two toxin-coding plasmids of 128 and 100 kb instead of the reported 225 kb plasmid, a co-integrate of the former two. In two of the sequenced strains, only a 100 kb plasmid was present. Some heterogeneity exists in the small plasmid content and structure between strains. These data support the perception of active plasmid exchange among B. thuringiensis subsp. israelensis strains in nature.

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.

Genome Sequence of the Atypical Symbiotic Frankia R43 Strain, a Nitrogen-Fixing and Hydrogen-Producing Actinobacterium.

Frankia strain R43 is a nitrogen-fixing and hydrogen-producing symbiotic actinobacterium that was isolated from nodules of Casuarina cunninghamiana but infects only Elaeagnaceae. This communication reports the genome of the strain R43 and provides insights into the microbe genomics and physiological potentials.

Genome Sequence of "Candidatus Arthromitus" sp. Strain SFB-Mouse-NL, a Commensal Bacterium with a Key Role in Postnatal Maturation of Gut Immune Functions.

"Candidatus Arthromitus" sp. strain SFB-mouse-NL (SFB, segmented filamentous bacteria) is a commensal bacterium necessary for inducing the postnatal maturation of homeostatic innate and adaptive immune responses in the mouse gut. Here, we report the genome sequence of this bacterium, which sets it apart from earlier sequenced mouse SFB isolates.

Bacillus cytotoxicus sp. nov. is a novel thermotolerant species of the Bacillus cereus Group occasionally associated with food poisoning.

An aerobic endospore-forming bacillus (NVH 391-98(T)) was isolated during a severe food poisoning outbreak in France in 1998, and four other similar strains have since been isolated, also mostly from food poisoning cases. Based on 16S rRNA gene sequence similarity, these strains were shown to belong to the Bacillus cereus Group (over 97% similarity with the current Group species) and phylogenetic distance from other validly described species of the genus Bacillus was less than 95%. Based on 16S rRNA gene sequence similarity and MLST data, these novel strains were shown to form a robust and well-separated cluster in the B. cereus Group, and constituted the most distant cluster from species of this Group. Major fatty acids (iso-C(15:0), C(16:0), iso-C(17:0), anteiso-C(15 : 0), iso-C(16:0), iso-C(13:0)) supported the affiliation of these strains to the genus Bacillus, and more specifically to the B. cereus Group. NVH 391-98(T) taxon was more specifically characterized by an abundance of iso-C(15:0) and low amounts of iso-C(13:0) compared with other members of the B. cereus Group. Genome similarity together with DNA-DNA hybridization values and physiological and biochemical tests made it possible to genotypically and phenotypically differentiate NVH 391-98(T) taxon from the six current B. cereus Group species. NVH 391-98(T) therefore represents a novel species, for which the name Bacillus cytotoxicus sp. nov. is proposed, with the type strain NVH 391-98(T) (= DSM 22905(T) = CIP 110041(T)).

Sequence Analysis of Inducible Prophage phIS3501 Integrated into the Haemolysin II Gene of Bacillus thuringiensis var israelensis ATCC35646.

Diarrheic food poisoning by bacteria of the Bacillus cereus group is mostly due to several toxins encoded in the genomes. One of them, cytotoxin K, was recently identified as responsible for severe necrotic syndromes. Cytotoxin K is similar to a class of proteins encoded by genes usually annotated as haemolysin II (hlyII) in the majority of genomes of the B. cereus group. The partially sequenced genome of Bacillus thuringiensis var israelensis ATCC35646 contains several potentially induced prophages, one of them integrated into the hlyII gene. We determined the complete sequence and established the genomic organization of this prophage-designated phIS3501. During induction of excision of this prophage with mitomycin C, intact hlyII gene is formed, thus providing to cells a genetic ability to synthesize the active toxin. Therefore, this prophage, upon its excision, can be implicated in the regulation of synthesis of the active toxin and thus in the virulence of bacterial host. A generality of selection for such systems in bacterial pathogens is indicated by the similarity of this genetic arrangement to that of Staphylococcus aureus ß-haemolysin.

Complete genome sequence of the highly hemolytic strain Bacillus cereus F837/76.

Highly hemolytic strain Bacillus cereus F837/76 was isolated in 1976 from a contaminated prostate wound. The complete nucleotide sequence of this strain reported here counts nearly 36,500 single-nucleotide differences from the closest sequenced strain, Bacillus thuringiensis Al Hakam. F827/76 also contains a 10-kb plasmid that was not detected in the Al Hakam strain.

Complete genome sequence of Lactococcus lactis subsp. cremoris A76.

We report the complete genome sequence of Lactococcus lactis subsp. cremoris A76, a dairy strain isolated from a cheese production outfit. Genome analysis detected two contiguous islands fitting to the L. lactis subsp. lactis rather than to the L. lactis subsp. cremoris lineage. This indicates the existence of genetic exchange between the diverse subspecies, presumably related to the technological process.

Functional and morphological adaptation to peptidoglycan precursor alteration in Lactococcus lactis.

Cell wall peptidoglycan assembly is a tightly regulated process requiring the combined action of multienzyme complexes. In this study we provide direct evidence showing that substrate transformations occurring at the different stages of this process play a crucial role in the spatial and temporal coordination of the cell wall synthesis machinery. Peptidoglycan substrate alteration was investigated in the Gram-positive bacterium Lactococcus lactis by substituting the peptidoglycan precursor biosynthesis genes of this bacterium for those of the vancomycin-resistant bacterium Lactobacillus plantarum. A set of L. lactis mutant strains in which the normal d-Ala-ended precursors were partially or totally replaced by d-Lac-ended precursors was generated. Incorporation of the altered precursor into the cell wall induced morphological changes arising from a defect in cell elongation and cell separation. Structural analysis of the muropeptides confirmed that the activity of multiple enzymes involved in peptidoglycan synthesis was altered. Optimization of this altered pathway was necessary to increase the level of vancomycin resistance conferred by the utilization of d-Lac-ended peptidoglycan precursors in the mutant strains. The implications of these findings on the control of bacterial cell morphogenesis and the mechanisms of vancomycin resistance are discussed.

The genetically remote pathogenic strain NVH391-98 of the Bacillus cereus group is representative of a cluster of thermophilic strains.

Bacteria of the Bacillus cereus group are known to cause food poisoning. A rare phylogenetically remote strain, NVH391-98, was recently characterized to encode a particularly efficient cytotoxin K presumably responsible for food poisoning. This pathogenic strain and its close relatives can be phenotypically distinguished from other strains of the B. cereus group by the inability to grow at temperatures below 17 degrees C and by the ability to grow at temperatures from 48 to 53 degrees C. A temperate phage, phBC391A2, residing in the genome of NVH391-98 allows us to distinguish the three known members of this thermophilic strain cluster.

Ecological diversification in the Bacillus cereus Group.

The Bacillus cereus Group comprises organisms that are widely distributed in the environment and are of health and economic interest. We demonstrate an 'ecotypic' structure of populations in the B. cereus Group using (i) molecular data from Fluorescent Amplified Fragment Length Polymorphism patterns, ribosomal gene sequences, partial panC gene sequences, 'psychrotolerant' DNA sequence signatures and (ii) phenotypic and descriptive data from range of growth temperature, psychrotolerance and thermal niches. Seven major phylogenetic groups (I to VII) were thus identified, with ecological differences that provide evidence for a multiemergence of psychrotolerance in the B. cereus Group. A moderate thermotolerant group (VII) was basal to the mesophilic group I, from which in turn distinct thermal lineages have emerged, comprising two mesophilic groups (III, IV), an intermediate group (V) and two psychrotolerant groups (VI, II). This stepwise evolutionary transition toward psychrotolerance was particularly well illustrated by the relative abundance of the 'psychrotolerant' rrs signature (as defined by Pruss et al.) copies accumulated in strains that varied according to the phylogenetic group. The 'psychrotolerant' cspA signature (as defined by Francis et al.) was specific to group VI and provided a useful way to differentiate it from the psychrotolerant group II. This study illustrates how adaptation to novel environments by the modification of temperature tolerance limits has shaped historical patterns of global ecological diversification in the B. cereus Group. The implications for the taxonomy of this Group and for the human health risk are discussed.

Extending the Bacillus cereus group genomics to putative food-borne pathogens of different toxicity.

The Bacillus cereus group represents sporulating soil bacteria containing pathogenic strains which may cause diarrheic or emetic food poisoning outbreaks. Multiple locus sequence typing revealed a presence in natural samples of these bacteria of about 30 clonal complexes. Application of genomic methods to this group was however biased due to the major interest for representatives closely related to Bacillus anthracis. Albeit the most important food-borne pathogens were not yet defined, existing data indicate that they are scattered all over the phylogenetic tree. The preliminary analysis of the sequences of three genomes discussed in this paper narrows down the gaps in our knowledge of the B. cereus group. The strain NVH391-98 is a rare but particularly severe food-borne pathogen. Sequencing revealed that the strain should be a representative of a novel bacterial species, for which the name Bacillus cytotoxis or Bacillus cytotoxicus is proposed. This strain has a reduced genome size compared to other B. cereus group strains. Genome analysis revealed absence of sigma B factor and the presence of genes encoding diarrheic Nhe toxin, not detected earlier. The strain B. cereus F837/76 represents a clonal complex close to that of B. anthracis. Including F837/76, three such B. cereus strains had been sequenced. Alignment of genomes suggests that B. anthracis is their common ancestor. Since such strains often emerge from clinical cases, they merit a special attention. The third strain, KBAB4, is a typical facultative psychrophile generally found in soil. Phylogenic studies show that in nature it is the most active group in terms of gene exchange. Genomic sequence revealed high presence of extra-chromosomal genetic material (about 530kb) that may account for this phenomenon. Genes coding Nhe-like toxin were found on a big plasmid in this strain. This may indicate a potential mechanism of toxicity spread from the psychrophile strain community. The results of this genomic work and ecological compartments of different strains incite to consider a necessity of creating prophylactic vaccines against bacteria closely related to NVH391-98 and F837/76. Presumably developing of such vaccines can be based on the properties of non-pathogenic strains such as KBAB4 or ATCC14579 reported here or earlier. By comparing the protein coding genes of strains being sequenced in this project to others we estimate the shared proteome, or core genome, in the B. cereus group to be 3000+/-200 genes and the total proteome, or pan-genome, to be 20-25,000 genes.

Multiple-locus sequence typing analysis of Bacillus cereus and Bacillus thuringiensis reveals separate clustering and a distinct population structure of psychrotrophic strains.

We used multilocus sequence typing (MLST) to characterize phylogenetic relationships for a collection of Bacillus cereus group strains isolated from forest soil in the Paris area during a mild winter. This collection contains multiple strains isolated from the same soil sample and strains isolated from samples from different sites. We characterized 115 strains of this collection and 19 other strains based on the sequences of the clpC, dinB, gdpD, panC, purF, and yhfL loci. The number of alleles ranged from 36 to 53, and a total of 93 allelic profiles or sequence types were distinguished. We identified three major strain clusters-C, T, and W-based on the comparison of individual gene sequences or concatenated sequences. Some less representative clusters and subclusters were also distinguished. Analysis of the MLST data using the concept of clonal complexes led to the identification of two, five, and three such groups in clusters C, T, and W, respectively. Some of the forest isolates were closely related to independently isolated psychrotrophic strains. Systematic testing of the strains of this collection showed that almost all the strains that were able to grow at a low temperature (6 degrees C) belonged to cluster W. Most of these strains, including three independently isolated strains, belong to two clonal complexes and are therefore very closely related genetically. These clonal complexes represent strains corresponding to the previously identified species Bacillus weihenstephanensis. Most of the other strains of our collection, including some from the W cluster, are not psychrotrophic. B. weihenstephanensis (cluster W) strains appear to comprise an effectively sexual population, whereas Bacillus thuringiensis (cluster T) and B. cereus (cluster C) have clonal population structures.

Clustered regularly interspaced short palindrome repeats (CRISPRs) have spacers of extrachromosomal origin.

Numerous prokaryote genomes contain structures known as clustered regularly interspaced short palindromic repeats (CRISPRs), composed of 25-50 bp repeats separated by unique sequence spacers of similar length. CRISPR structures are found in the vicinity of four genes named cas1 to cas4. In silico analysis revealed another cluster of three genes associated with CRISPR structures in many bacterial species, named here as cas1B, cas5 and cas6, and also revealed a certain number of spacers that have homology with extant genes, most frequently derived from phages, but also derived from other extrachromosomal elements. Sequence analysis of CRISPR structures from 24 strains of Streptococcus thermophilus and Streptococcus vestibularis confirmed the homology of spacers with extrachromosomal elements. Phage sensitivity of S. thermophilus strains appears to be correlated with the number of spacers in the CRISPR locus the strain carries. The authors suggest that the spacer elements are the traces of past invasions by extrachromosomal elements, and hypothesize that they provide the cell immunity against phage infection, and more generally foreign DNA expression, by coding an anti-sense RNA. The presence of gene fragments in CRISPR structures and the nuclease motifs in cas genes of both cluster types suggests that CRISPR formation involves a DNA degradation step.