Survey of group I and group II introns in 29 sequenced genomes of the Bacillus cereus group: insights into their spread and evolution

Group I and group II introns are different catalytic self-splicing and mobile RNA elements that contribute to genome dynamics. In this study, we have analyzed their distribution and evolution in 29 sequenced genomes from the Bacillus cereus group of bacteria. Introns were of different structural classes and evolutionary origins, and a large number of nearly identical elements are shared between multiple strains of different sources, suggesting recent lateral transfers and/or that introns are under a strong selection pressure. Altogether, 73 group I introns were identified, inserted in essential genes from the chromosome or newly described prophages, including the first elements found within phages in bacterial plasmids. Notably, bacteriophages are an important source for spreading group I introns between strains. Furthermore, 77 group II introns were found within a diverse set of chromosomal and plasmidic genes. Unusual findings include elements located within conserved DNA metabolism and repair genes and one intron inserted within a novel retroelement. Group II introns are mainly disseminated via plasmids and can subsequently invade the host genome, in particular by coupling mobility with host cell replication. This study reveals a very high diversity and variability of mobile introns in B. cereus group strains.     

Selection-Driven Extinction Dynamics for Group II Introns in Enterobacteriales

Transposable elements (TEs) are one of the major driving forces of genome evolution, raising the question of the long-term dynamics underlying their evolutionary success. Some TEs were proposed to evolve under a pattern of periodic extinctions-recolonizations, in which elements recurrently invade and quickly proliferate within their host genomes, then start to disappear until total extinction. Depending on the model, TE extinction is assumed to be driven by purifying selection against colonized host genomes (Sel-DE model) or by saturation of host genomes (Sat-DE model). Bacterial group II introns are suspected to follow an extinction-recolonization model of evolution, but whether they follow Sel-DE or Sat-DE dynamics is not known. Our analysis of almost 200 group II intron copies from 90 sequenced Enterobacteriales genomes confirms their extinction-recolonization dynamics: patchy element distributions among genera and even among strains within genera, acquisition of new group II introns through plasmids or other mobile genetic elements, and evidence for recent proliferations in some genomes. Distributions of recent and past proliferations and of their respective homing sites further provide strong support for the Sel-DE model, suggesting that group II introns are deleterious to their hosts. Overall, our observations emphasize the critical impact of host properties on TE dynamics.     

IS150: distribution, nucleotide sequence and phylogenetic relationships of a new E. coli insertion element.

Recently we identified the new insertion (IS) sequence IS150 in various strains of Escherichia coli K-12. We have screened other strains of E. coli and Salmonella typhimurium for the presence of homologous sequences. The strains of E. coli K-12 and W tested contain one or more copies of homology to IS150. We have also determined the complete nucleotide sequence of a copy of IS150 inserted into IS1. Comparison of nucleotide and deduced amino acid sequences of IS150, IS2, IS3, IS51, IS600 and IS629 reveals significant homologies suggesting that these elements are members of a family of phylogenetically related insertion sequences.     

The molecular evolution and structural organization of self-splicing group I introns at position 516 in nuclear SSU rDNA of myxomycetes

Group I introns are relatively common within nuclear ribosomal DNA of eukaryotic microorganisms, especially in myxomycetes. Introns at position S516 in the small subunit ribosomal RNA gene are particularly common, but have a sporadic occurrence in myxomycetes. Fuligo septica, Badhamia gracilis, and Physarum flavicomum, all members of the family Physaraceae, contain related group IC1 introns at this site. The F. septica intron was studied at the molecular level and found to self-splice as naked RNA and to generate full-length intron RNA circles during incubation. Group I introns at position S516 appear to have a particularly widespread distribution among protists and fungi. Secondary structural analysis of more than 140 S516 group I introns available in the database revealed five different types of organization, including IC1 introns with and without His-Cys homing endonuclease genes, complex twin-ribozyme introns, IE introns, and degenerate group I-like introns. Both intron structural and phylogenetic analyses indicate a multiple origin of the S516 introns during evolution. The myxomycete introns are related to S516 introns in the more distantly related brown algae and Acanthamoeba species. Possible mechanisms of intron transfer both at the RNA- and DNA-levels are discussed in order to explain the observed widespread, but scattered, phylogenetic distribution.     

Method for selection of transposable DNA and characterization of a new insertion sequence, IS493, from Streptomyces lividans.

A method to select for transposable elements from Streptomyces spp. by using insertional inactivation of a repressor gene that functions in Escherichia coli was developed. Plasmid pCZA126, which can replicate in Streptomyces spp. or E. coli, contains a gene coding for the lambda cI857 repressor and a gene, under repressor control, coding for apramycin resistance. E. coli cells containing the plasmid are apramycin sensitive but become apramycin resistant if the cI857 repressor gene is disrupted. Plasmids propagated in Streptomyces spp. can be screened for transposable elements that have disrupted the cI857 gene by transforming E. coli cells to apramycin resistance. This method was used to isolate a new 1.6-kilobase insertion sequence, IS493, from Streptomyces lividans CT2. IS493 duplicated host DNA at the target site, had inverted repeats at its ends, and contained two tandem open reading frames on each strand. IS493 was present in three copies in the same genomic locations in several S. lividans strains. Two of the copies appeared to be present in regions of similar DNA context that extended at least 11.5 kilobases. Several other Streptomyces spp. did not appear to contain copies of IS493.     

The evolution of homing endonuclease genes and group I introns in nuclear rDNA

Group I introns are autonomous genetic elements that can catalyze their own excision from pre-RNA. Understanding how group I introns move in nuclear ribosomal (r)DNA remains an important question in evolutionary biology. Two models are invoked to explain group I intron movement. The first is termed homing and results from the action of an intron-encoded homing endonuclease that recognizes and cleaves an intronless allele at or near the intron insertion site. Alternatively, introns can be inserted into RNA through reverse splicing. Here, we present the sequences of two large group I introns from fungal nuclear rDNA, which both encode putative full-length homing endonuclease genes (HEGs). Five remnant HEGs in different fungal species are also reported. This brings the total number of known nuclear HEGs from 15 to 22. We determined the phylogeny of all known nuclear HEGs and their associated introns. We found evidence for intron-independent HEG invasion into both homologous and heterologous introns in often distantly related lineages, as well as the "switching" of HEGs between different intron peripheral loops and between sense and antisense strands of intron DNA. These results suggest that nuclear HEGs are frequently mobilized. HEG invasion appears, however, to be limited to existing introns in the same or neighboring sites. To study the intron-HEG relationship in more detail, the S943 group I intron in fungal small-subunit rDNA was used as a model system. The S943 HEG is shown to be widely distributed as functional, inactivated, or remnant ORFs in S943 introns.     

Characterization of the lactococcal group II intron target site in its native host

The Lactococcus lactis group II intron (Ll.ltrB) retrohomes into the ltrB gene at high efficiency. To date, the critical DNA bases recognized in vivo by the Ll.ltrB ribonucleoprotein (RNP) have been exclusively elucidated in Escherichia coli. However, recent evidence indicates host-dependant differences in Ll.ltrB mobility, raising the possibility of limitations of the current model for RNP-homing site recognition in the native L. lactis host. In this work, intron retargeting experiments in L. lactis have demonstrated that adherence to specific target site critical bases is not sufficient to predict success or failure of chromosomal invasion, as in E. coli. Accordingly, a quantitative real-time PCR (QPCR) assay was developed to test target site nucleotides previously demonstrated as critical for homing in E. coli, for relevance in its native host. This two-plasmid QPCR homing assay is highly sensitive and, unlike previous E. coli-based assays, resolves differential homing efficiencies in the absence of selection. As in E. coli, deviation from wild type at target site positions -23, -21, -20, -19, and +5 resulted in lower homing efficiencies in L. lactis. Furthermore, the same trends are observed when assaying select variants in Enterococcus faecalis. Our results suggest that these target site positions are critical in both E. coli and L. lactis.     

Nonrandom location of IS1 elements in the genomes of natural isolates of Escherichia coli

We have studied the spatial distribution of IS1 elements in the genomes of natural isolates comprising the ECOR reference collection of Escherichia coli. We find evidence for nonrandomness at three levels. Many pairs of IS1 elements are in much closer proximity (< 10 kb) than can be accounted for by chance. IS1 elements in close proximity were identified by long-range PCR amplification of the genomic sequence between them. Each amplified region was sequenced and its map location determined by database screening of DNA hybridization. Among the ECOR strains with at least two IS1 elements, 54% had one or more pairs of elements separated by < 10 kb. We propose that this type of clustering is a result of "local hopping," in which we assume that a significant proportion of tranposition events leads to the insertion of a daughter IS element in the vicinity of the parental element. A second level of nonrandomness is found in strains with a modest number of IS1 elements that are mapped through the use of inverse PCR to amplify flanking genomic sequences: in these strains, the insertion sites tend to be clustered over a smaller region of chromosome than would be expected by chance. A third level of nonrandomness is observed in the composite distribution of IS elements across strains: among 20 mapped IS1 elements, none were found in the region of 48-77 minutes, a significant gap. One region of the E. coli chromosome, at 98 min, had a cluster of IS1 elements in seven ECOR strains of diverse phylogenetic origin. We deduce from sequence analysis that this pattern of distribution is a result of initial insertion in the most recent common ancestor of these strains and therefore not a hot spot of insertion. Analysis using long- range PCR with primers for IS2 and IS3 also yielded pairs of elements in close proximity, suggesting that these elements may also occasionally transpose by local hopping.      

Diversity of group II introns in the genome of <Emphasis Type="Italic"> Sinorhizobium meliloti </Emphasis>strain 1021: splicing and mobility of RmInt1

The number and diversity of known group II introns in eubacteria are continually increasing with the addition of new data from sequencing projects, but the significance of these introns in the evolution of bacterial genomes is unknown. We analyzed the main features of the group II introns present in the genome of the soil microorganism Sinorhizobium meliloti (strain 1021), the nitrogen-fixing symbiont of alfalfa, the DNA sequence of which was recently determined. Strain 1021 harbors three different classes of group II introns: RmInt1, of bacterial class D; SMb2147/SMb21167, which cluster within bacterial class C; and SMa1875, the phylogenetic class of which is uncertain. The group II introns SMb2147/SMb21167 and SMa1875 are widely distributed in S. meliloti, but are present in lower copy numbers than RmInt1. Strain 1021 harbors three copies of RmInt1, which is pSym-specific. Although RmInt1 is spliced in strain 1021, mobility assays suggested that, in contrast to other S. meliloti strains, the genetic background of strain 1021 does not support intron homing events.