PADLOC: a web server for the identification of antiviral defence systems in microbial genomes.

Most bacteria and archaea possess multiple antiviral defence systems that protect against infection by phages, archaeal viruses and mobile genetic elements. Our understanding of the diversity of defence systems has increased greatly in the last few years, and many more systems likely await discovery. To identify defence-related genes, we recently developed the Prokaryotic Antiviral Defence LOCator (PADLOC) bioinformatics tool. To increase the accessibility of PADLOC, we describe here the PADLOC web server (freely available at https://padloc.otago.ac.nz), allowing users to analyse whole genomes, metagenomic contigs, plasmids, phages and archaeal viruses. The web server includes a more than 5-fold increase in defence system types detected (since the first release) and expanded functionality enabling detection of CRISPR arrays and retron ncRNAs. Here, we provide user information such as input options, description of the multiple outputs, limitations and considerations for interpretation of the results, and guidance for subsequent analyses. The PADLOC web server also houses a precomputed database of the defence systems in > 230,000 RefSeq genomes. These data reveal two taxa, Campylobacterota and Spriochaetota, with unusual defence system diversity and abundance. Overall, the PADLOC web server provides a convenient and accessible resource for the detection of antiviral defence systems.

UG/Abi: a highly diverse family of prokaryotic reverse transcriptases associated with defense functions.

Reverse transcriptases (RTs) are enzymes capable of synthesizing DNA using RNA as a template. Within the last few years, a burst of research has led to the discovery of novel prokaryotic RTs with diverse antiviral properties, such as DRTs (Defense-associated RTs), which belong to the so-called group of unknown RTs (UG) and are closely related to the Abortive Infection system (Abi) RTs. In this work, we performed a systematic analysis of UG and Abi RTs, increasing the number of UG/Abi members up to 42 highly diverse groups, most of which are predicted to be functionally associated with other gene(s) or domain(s). Based on this information, we classified these systems into three major classes. In addition, we reveal that most of these groups are associated with defense functions and/or mobile genetic elements, and demonstrate the antiphage role of four novel groups. Besides, we highlight the presence of one of these systems in novel families of human gut viruses infecting members of the Bacteroidetes and Firmicutes phyla. This work lays the foundation for a comprehensive and unified understanding of these highly diverse RTs with enormous biotechnological potential.

Systematic prediction of genes functionally associated with bacterial retrons and classification of the encoded tripartite systems.

Bacterial retrons consist of a reverse transcriptase (RT) and a contiguous non-coding RNA (ncRNA) gene. One third of annotated retrons carry additional open reading frames (ORFs), the contribution and significance of which in retron biology remains to be determined. In this study we developed a computational pipeline for the systematic prediction of genes specifically associated with retron RTs based on a previously reported large dataset representative of the diversity of prokaryotic RTs. We found that retrons generally comprise a tripartite system composed of the ncRNA, the RT and an additional protein or RT-fused domain with diverse enzymatic functions. These retron systems are highly modular, and their components have coevolved to different extents. Based on the additional module, we classified retrons into 13 types, some of which include additional variants. Our findings provide a basis for future studies on the biological function of retrons and for expanding their biotechnological applications.

Spacer acquisition from RNA mediated by a natural reverse transcriptase-Cas1 fusion protein associated with a type III-D CRISPR-Cas system in Vibrio vulnificus.

The association of reverse transcriptases (RTs) with CRISPR-Cas system has recently attracted interest because the RT activity appears to facilitate the RT-dependent acquisition of spacers from RNA molecules. However, our understanding of this spacer acquisition process remains limited. We characterized the in vivo acquisition of spacers mediated by an RT-Cas1 fusion protein linked to a type III-D system from Vibrio vulnificus strain YJ016, and showed that the adaptation module, consisting of the RT-Cas1 fusion, two different Cas2 proteins (A and B) and one of the two CRISPR arrays, was completely functional in a heterologous host. We found that mutations of the active site of the RT domain significantly decreased the acquisition of new spacers and showed that this RT-Cas1-associated adaptation module was able to incorporate spacers from RNA molecules into the CRISPR array. We demonstrated that the two Cas2 proteins of the adaptation module were required for spacer acquisition. Furthermore, we found that several sequence-specific features were required for the acquisition and integration of spacers derived from any region of the genome, with no bias along the 5'and 3'ends of coding sequences. This study provides new insight into the RT-Cas1 fusion protein-mediated acquisition of spacers from RNA molecules.

A group II intron-encoded protein interacts with the cellular replicative machinery through the ß-sliding clamp.

Group II introns are self-splicing mobile genetic retroelements. The spliced intron RNA and the intron-encoded protein (IEP) form ribonucleoprotein particles (RNPs) that recognize and invade specific DNA target sites. The IEP is a reverse transcriptase/maturase that may bear a C-terminal endonuclease domain enabling the RNP to cleave the target DNA strand to prime reverse transcription. However, some mobile introns, such as RmInt1, lack the En domain but nevertheless retrohome efficiently to transient single-stranded DNA target sites at a DNA replication fork. Their mobility is associated with host DNA replication, and they use the nascent lagging strand as a primer for reverse transcription. We searched for proteins that interact with RmInt1 RNPs and direct these RNPs to the DNA replication fork. Co-immunoprecipitation assays suggested that DnaN (the ß-sliding clamp), a component of DNA polymerase III, interacts with the protein component of the RmInt1 RNP. Pulldown assays, far-western blots and biolayer interferometry supported this interaction. Peptide binding assays also identified a putative DnaN-interacting motif in the RmInt1 IEP structurally conserved in group II intron IEPs. Our results suggest that intron RNP interacts with the ß-sliding clamp of the DNA replication machinery, favouring reverse splicing into the transient ssDNA at DNA replication forks.

DNA cleavage and reverse splicing of ribonucleoprotein particles reconstituted in vitro with linear RmInt1 RNA.

The RmInt1 group II intron is an efficient self-splicing mobile retroelement that catalyzes its own excision as lariat, linear and circular molecules. In vivo, the RmInt1 lariat and the reverse transcriptase (IEP) it encodes form a ribonucleoprotein particle (RNP) that recognizes the DNA target for site-specific full intron insertion via a two-step reverse splicing reaction. RNPs containing linear group II intron RNA are generally thought to be unable to complete the reverse splicing reaction. Here, we show that reconstituted in vitro RNPs containing linear RmInt1 ΔORF RNA can mediate the cleavage of single-stranded DNA substrates in a very precise manner with the attachment of the intron RNA to the 3´exon as the first step of a reverse splicing reaction. Notably, we also observe molecules in which the 5´exon is linked to the RmInt1 RNA, suggesting the completion of the reverse splicing reaction, albeit rather low and inefficiently. That process depends on DNA target recognition and can be successful completed by RmInt1 RNPs with linear RNA displaying 5´ modifications.

Multiple origins of reverse transcriptases linked to CRISPR-Cas systems.

Prokaryotic genomes harbour a plethora of uncharacterized reverse transcriptases (RTs). RTs phylogenetically related to those encoded by group-II introns have been found associated with type III CRISPR-Cas systems, adjacent or fused at the C-terminus to Cas1. It is thought that these RTs may have a relevant function in the CRISPR immune response mediating spacer acquisition from RNA molecules. The origin and relationships of these RTs and the ways in which the various protein domains evolved remain matters of debate. We carried out a large survey of annotated RTs in databases (198,760 sequences) and constructed a large dataset of unique representative sequences (9,141). The combined phylogenetic reconstruction and identification of the RTs and their various protein domains in the vicinity of CRISPR adaptation and effector modules revealed three different origins for these RTs, consistent with their emergence on multiple occasions: a larger group that have evolved from group-II intron RTs, and two minor lineages that may have arisen more recently from Retron/retron-like sequences and Abi-P2 RTs, the latter associated with type I-C systems. We also identified a particular group of RTs associated with CRISPR-cas loci in clade 12, fused C-terminally to an archaeo-eukaryotic primase (AEP), a protein domain (AE-Prim_S_like) forming a particular family within the AEP proper clade. Together, these data provide new insight into the evolution of CRISPR-Cas/RT systems.

The early events underlying genome evolution in a localized Sinorhizobium meliloti population.

BACKGROUND: Population genetic analyses based on genome-wide sequencing data have been carried out for Sinorhizobium medicae and S. meliloti, two closely related bacterial species forming nitrogen-fixing symbioses with plants of the genus Medicago. However, genome coverage was low or the isolates had a broad geographic distribution, making it difficult to interpret the estimated diversity and to unravel the early events underlying population genetic variations and ecological differentiation. RESULTS: Here, to gain insight into the early genome level variation and diversification within S. meliloti populations, we first used Illumina paired-end reads technology to sequence a new clone of S. meliloti strain GR4, a highly competitive strain for alfalfa nodulation. The Illumina data and the GR4 genome sequence previously obtained with 454 technology were used to generate a high-quality reference genome sequence. We then used Illumina technology to sequence the genomes of 13 S. meliloti isolates representative of the genomic variation within the GR4-type population, obtained from a single field site with a high degree of coverage. The genome sequences obtained were analyzed to determine nucleotide diversity, divergence times, polymorphism and genomic variation. Similar low levels of nucleotide diversity were observed for the chromosome, pSymB and pSymA replicons. The isolates displayed other types of variation, such as indels, recombination events, genomic island excision and the transposition of mobile elements. CONCLUSIONS: Our results suggest that the GR4-type population has experienced a process of demographic expansion and behaves as a stable genotypic cluster of genome-wide similarity, with most of the genome following a clonal pattern of evolution. Although some of genetic variation detected within the GR4-type population is probably due to genetic drift, others might be important in diversification and environmental adaptation.

In vitro characterization of the splicing efficiency and fidelity of the RmInt1 group II intron as a means of controlling the dispersion of its host mobile element.

Group II introns are catalytic RNAs that are excised from their precursors in a protein-dependent manner in vivo. Certain group II introns can also react in a protein-independent manner under nonphysiological conditions in vitro. The efficiency and fidelity of the splicing reaction is crucial, to guarantee the correct formation and expression of the protein-coding mRNA. RmInt1 is an efficient mobile intron found within the ISRm2011-2 insertion sequence in the symbiotic bacterium Sinorhizobium meliloti. The RmInt1 intron self-splices in vitro, but this reaction generates side products due to a predicted cryptic IBS1* sequence within the 3' exon. We engineered an RmInt1 intron lacking the cryptic IBS1* sequence, which improved the fidelity of the splicing reaction. However, atypical circular forms of similar electrophoretic mobility to the lariat intron were nevertheless observed. We analyzed a run of four cytidine residues at the 3' splice site potentially responsible for a lack of fidelity at this site leading to the formation of circular intron forms. We showed that mutations of residues base-pairing in the tertiary EBS3-IBS3 interaction increased the efficiency and fidelity of the splicing reaction. Our results indicate that RmInt1 has developed strategies for decreasing its splicing efficiency and fidelity. RmInt1 makes use of unproductive splicing reactions to limit the transposition of the insertion sequence into which it inserts itself in its natural context, thereby preventing potentially harmful dispersion of ISRm2011-2 throughout the genome of its host.

Bacterial Diversity in the Soda Saline Crater Lake from Isabel Island, Mexico.

Isabel Lake is a moderate saline soda crater lake located in Isabel Island in the eastern tropical Pacific coast of Mexico. Lake is mainly formed by rainfall and is strongly affected by evaporation and high input of nutrients derived from excretions of a large bird community inhabiting the island. So far, only the island macrobiota has been studied. The knowledge of the prokaryotic biota inhabiting the upper layers of this meromictic lake can give clues for the maintenance of this ecosystem. We assessed the diversity and composition of prokaryotic community in sediments and water of the lake by DGGE profiling, 16S rRNA gene amplicon pyrosequencing, and cultivation techniques. The bacterial community is largely dominated by halophilic and halotolerant microorganisms. Alpha diversity estimations reveal higher value in sediments than in water (P > 0.005). The lake water is dominated by γ-Proteobacteria belonging to four main families where Halomonadaceae presents the highest abundance. Aerobic, phototrophic, and halotolerant prokaryotes such as Cyanobacteria GPIIa, Halomonas, Alcanivorax, Idiomarina, and Cyclobacterium genera are commonly found. However, in sediment samples, Formosa, Muricauda, and Salegentibacter genera corresponding to Flavobacteriaceae family accounted for 15-20 % of the diversity. Heterotrophs like those involved in sulfur cycle, Desulfotignum, Desulfuromonas, Desulfofustis, and Desulfopila, appear to play an important role in sediments. Finally, a collection of aerobic halophilic bacterial isolates was created from these samples; members of the genus Halomonas were predominantly isolated from lake water. This study contributes to state the bacterial diversity present in this particular soda saline crater lake.

Metagenomic assessment of the potential microbial nitrogen pathways in the rhizosphere of a mediterranean forest after a wildfire.

Wildfires are frequent in the forests of the Mediterranean Basin and have greatly influenced this ecosystem. Changes to the physical and chemical properties of the soil, due to fire and post-fire conditions, result in alterations of both the bacterial communities and the nitrogen cycle. We explored the effects of a holm oak forest wildfire on the rhizospheric bacterial communities involved in the nitrogen cycle. Metagenomic data of the genes involved in the nitrogen cycle showed that both the undisturbed and burned rhizospheres had a conservative nitrogen cycle with a larger number of sequences related to the nitrogen incorporation pathways and a lower number for nitrogen output. However, the burned rhizosphere showed a statistically significant increase in the number of sequences for nitrogen incorporation (allantoin utilization and nitrogen fixation) and a significantly lower number of sequences for denitrification and dissimilatory nitrite reductase subsystems, possibly in order to compensate for nitrogen loss from the soil after burning. The genetic potential for nitrogen incorporation into the ecosystem was assessed through the diversity of the nitrogenase reductase enzyme, which is encoded by the nifH gene. We found that nifH gene diversity and richness were lower in burned than in undisturbed rhizospheric soils. The structure of the bacterial communities involved in the nitrogen cycle showed a statistically significant increase of Actinobacteria and Firmicutes phyla after the wildfire. Both approaches showed the important role of gram-positive bacteria in the ecosystem after a wildfire.

Use of the computer-retargeted group II intron RmInt1 of Sinorhizobium meliloti for gene targeting.

Gene-targeting vectors derived from mobile group II introns capable of forming a ribonucleoprotein (RNP) complex containing excised intron lariat RNA and an intron-encoded protein (IEP) with reverse transcriptase (RT), maturase, and endonuclease (En) activities have been described. RmInt1 is an efficient mobile group II intron with an IEP lacking the En domain. We performed a comprehensive study of the rules governing RmInt1 target site recognition based on selection experiments with donor and recipient plasmid libraries, with randomization of the elements of the intron RNA involved in target recognition and the wild-type target site. The data obtained were used to develop a computer algorithm for identifying potential RmInt1 targets in any DNA sequence. Using this algorithm, we modified RmInt1 for the efficient recognition of DNA target sites at different locations in the Sinorhizobium meliloti chromosome. The retargeted RmInt1 integrated efficiently into the chromosome, regardless of the location of the target gene. Our results suggest that RmInt1 could be efficiently adapted for gene targeting.

Insights into the strategies used by related group II introns to adapt successfully for the colonisation of a bacterial genome.

Group II introns are self-splicing RNAs and site-specific mobile retroelements found in bacterial and organellar genomes. The group II intron RmInt1 is present at high copy number in Sinorhizobium meliloti species, and has a multifunctional intron-encoded protein (IEP) with reverse transcriptase/maturase activities, but lacking the DNA-binding and endonuclease domains. We characterized two RmInt1-related group II introns RmInt2 from S. meliloti strain GR4 and Sr.md.I1 from S. medicae strain WSM419 in terms of splicing and mobility activities. We used both wild-type and engineered intron-donor constructs based on ribozyme ΔORF-coding sequence derivatives, and we determined the DNA target requirements for RmInt2, the element most distantly related to RmInt1. The excision and mobility patterns of intron-donor constructs expressing different combinations of IEP and intron RNA provided experimental evidence for the co-operation of IEPs and intron RNAs from related elements in intron splicing and, in some cases, in intron homing. We were also able to identify the DNA target regions recognized by these IEPs lacking the DNA endonuclease domain. Our results provide new insight into the versatility of related group II introns and the possible co-operation between these elements to facilitate the colonization of bacterial genomes.

Genome-wide profiling of Hfq-binding RNAs uncovers extensive post-transcriptional rewiring of major stress response and symbiotic regulons in Sinorhizobium meliloti.

The RNA chaperone Hfq is a global post-transcriptional regulator in bacteria. Here, we used RNAseq to analyze RNA populations from the legume symbiont Sinorhizobium meliloti that were co-immunoprecipitated (CoIP-RNA) with a FLAG-tagged Hfq in five growth/stress conditions. Hfq-bound transcripts (1315) were largely identified in stressed bacteria and derived from small RNAs (sRNAs), both trans-encoded (6.4%) and antisense (asRNAs; 6.3%), and mRNAs (86%). Pull-down with Hfq recovered a small proportion of annotated S. meliloti sRNAs (14% of trans-sRNAs and 2% of asRNAs) suggesting a discrete impact of this protein in sRNA pathways. Nonetheless, Hfq selectively stabilized CoIP-enriched sRNAs, anticipating that these interactions are functionally significant. Transcription of 26 Hfq-bound sRNAs was predicted to occur from promoters recognized by the major stress σ factors σ(E2) or σ(H1/2). Recovery rates of sRNAs in each of the CoIP-RNA libraries suggest a large impact of Hfq-assisted riboregulation in S. meliloti osmoadaptation. Hfq directly targeted 18% of the predicted S. meliloti mRNAs, which encode functionally diverse proteins involved in transport and metabolism, σ(E2)-dependent stress responses, quorum sensing, flagella biosynthesis, ribosome, and membrane assembly or symbiotic nitrogen fixation. Canonical targeting of the 5' regions of two of the ABC transporter mRNAs by the homologous Hfq-binding AbcR1 and AbcR2 sRNAs leading to inhibition of protein synthesis was confirmed in vivo. We therefore provide a comprehensive resource for the systems-level deciphering of hitherto unexplored S. meliloti stress and symbiotic post-transcriptional regulons and the identification of Hfq-dependent sRNA-mRNA regulatory pairs.

Splicing of the Sinorhizobium meliloti RmInt1 group II intron provides evidence of retroelement behavior.

Group II introns act as both large catalytic RNAs and mobile retroelements. They are found in organelle and bacterial genomes and are spliced via a lariat intermediate, in a mechanism similar to that of spliceosomal introns. However, their distribution and insertion patterns, particularly for bacterial group II introns, suggest that they function and behave more like retroelements than organelle introns. RmInt1 is an efficient mobile intron found within the ISRm2011-2 insertion sequence in the symbiotic bacterium Sinorhizobium meliloti. This group II intron is excised, in vivo and in vitro, as intron lariats. However, the complete splicing reaction in vivo remains to be elucidated. A lacZ reporter gene system, northern blotting and real-time reverse transcription were carried out to investigate RmInt1 splicing activity. Splicing efficiency of 0.07 ± 0.02% was recorded. These findings suggest that bacterial group II introns function more like retroelements than spliceosomal introns. Their location is consistent with a role for these introns in preventing the spread of other potentially harmful mobile elements in bacteria.

Relevance of the branch point adenosine, coordination loop, and 3' exon binding site for in vivo excision of the Sinorhizobium meliloti group II intron RmInt1.

Excision of the bacterial group II intron RmInt1 has been demonstrated in vivo, resulting in the formation of both intron lariat and putative intron RNA circles. We show here that the bulged adenosine in domain VI of RmInt1 is required for splicing via the branching pathway, but branch site mutants produce small numbers of RNA molecules in which the first G residue of the intron is linked to the last C residue. Mutations in the coordination loop in domain I reduced splicing efficiency, but branched templates clearly predominated among splicing products. We also found that a single substitution at the EBS3 position (G329C), preventing EBS3-IBS3 pairing, resulted in the production of 50 to 100 times more RNA molecules in which the 5' and 3' extremities were joined. We provide evidence that these intron molecules may correspond to both, intron circles linked by a 2'-5' phosphodiester bond, and tandem, head-to-tail intron copies.

A survey of sRNA families in α-proteobacteria.

We have performed a computational comparative analysis of six small non-coding RNA (sRNA) families in α-proteobacteria. Members of these families were first identified in the intergenic regions of the nitrogen-fixing endosymbiont S. meliloti by a combined bioinformatics screen followed by experimental verification. Consensus secondary structures inferred from covariance models for each sRNA family evidenced in some cases conserved motifs putatively relevant to the function of trans-encoded base-pairing sRNAs i.e., Hfq-binding signatures and exposed anti Shine-Dalgarno sequences. Two particular family models, namely αr15 and αr35, shared own sub-structural modules with the Rfam model suhB (RF00519) and the uncharacterized sRNA family αr35b, respectively. A third sRNA family, termed αr45, has homology to the cis-acting regulatory element speF (RF00518). However, new experimental data further confirmed that the S. meliloti αr45 representative is an Hfq-binding sRNA processed from or expressed independently of speF, thus refining the Rfam speF model annotation. All the six families have members in phylogenetically related plant-interacting bacteria and animal pathogens of the order of the Rhizobiales, some occurring with high levels of paralogy in individual genomes. In silico and experimental evidences predict differential regulation of paralogous sRNAs in S. meliloti 1021. The distribution patterns of these sRNA families suggest major contributions of vertical inheritance and extensive ancestral duplication events to the evolution of sRNAs in plant-interacting bacteria.

Bacterial community structure in the rhizosphere of three cactus species from semi-arid highlands in central Mexico.

The nature reserve of Tehuacan-Cuicatlan in central Mexico is known for its diversity and endemism mainly in cactus plants. Although the xerophytic flora is reasonably documented, the bacterial communities associated with these species have been largely neglected. We assessed the diversity and composition of bacterial communities in bulk (non-rhizospheric) soil and the rhizosphere of three cactus plant species: Mammillaria carnea, Opuntia pilifera and Stenocereus stellatus, approached using cultivation and molecular techniques, considering the possible effect of dry and rainy seasons. Cultivation-dependent methods were focused on putative N(2)-fixers and heterotrophic aerobic bacteria, in the two media tested the values obtained for dry season samples grouped together regardless of the sample type (rhizospheric or non-rhizospheric), these groups also included the non-rhizospheric sample for rainy season, on each medium. These CFU values were smaller and significantly different from those obtained on rhizospheric samples from rainy season. Genera composition among isolates of the rhizospheric samples was very similar for each season, the most abundant taxa being α-Proteobacteria, Actinobacteria and Firmicutes. Interestingly, the genus Ochrobactrum was highly represented among rhizospheric samples, when cultured in N-free medium. The structure of the bacterial communities was approached with molecular techniques targeting partial 16S rRNA sequences such as denaturing gradient gel electrophoresis and serial analysis of ribosomal sequence tags. Under these approaches, the most represented bacterial phyla were Actinobacteria, Proteobacteria and Acidobacteria. The first two were also highly represented when using isolation techniques.

Identification of Group II Intron RmInt1 Binding Sites in a Bacterial Genome.

RmInt1 is a group II intron encoding a reverse transcriptase protein (IEP) lacking the C-terminal endonuclease domain. RmInt1 is an efficient mobile retroelement that predominantly reverse splices into the transient single-stranded DNA at the template for lagging strand DNA synthesis during host replication, a process facilitated by the interaction of the RmInt1 IEP with DnaN at the replication fork. It has been suggested that group II intron ribonucleoprotein particles bind DNA nonspecifically, and then scan for their correct target site. In this study, we investigated RmInt1 binding sites throughout the Sinorhizobium meliloti genome, by chromatin-immunoprecipitation coupled with next-generation sequencing. We found that RmInt1 binding sites cluster around the bidirectional replication origin of each of the three replicons comprising the S. meliloti genome. Our results provide new evidence linking group II intron mobility to host DNA replication.

Exon sequence requirements for excision in vivo of the bacterial group II intron RmInt1.

BACKGROUND: Group II intron splicing proceeds through two sequential transesterification reactions in which the 5' and 3'-exons are joined together and the lariat intron is released. The intron-encoded protein (IEP) assists the splicing of the intron in vivo and remains bound to the excised intron lariat RNA in a ribonucleoprotein particle (RNP) that promotes intron mobility. Exon recognition occurs through base-pairing interactions between two guide sequences on the ribozyme domain dI known as EBS1 and EBS2 and two stretches of sequence known as IBS1 and IBS2 on the 5' exon, whereas the 3' exon is recognized through interaction with the sequence immediately upstream from EBS1 [(δ-δ' interaction (subgroup IIA)] or with a nucleotide [(EBS3-IBS3 interaction (subgroup IIB and IIC))] located in the coordination-loop of dI. The δ nucleotide is involved in base pairing with another intron residue (δ') in subgroup IIB introns and this interaction facilitates base pairing between the 5' exon and the intron. RESULTS: In this study, we investigated nucleotide requirements in the distal 5'- and 3' exon regions, EBS-IBS interactions and δ-δ' pairing for excision of the group IIB intron RmInt1 in vivo. We found that the EBS1-IBS1 interaction was required and sufficient for RmInt1 excision. In addition, we provide evidence for the occurrence of canonical δ-δ' pairing and its importance for the intron excision in vivo. CONCLUSIONS: The excision in vivo of the RmInt1 intron is a favored process, with very few constraints for sequence recognition in both the 5' and 3'-exons. Our results contribute to understand how group II introns spread in nature, and might facilitate the use of RmInt1 in gene targeting.