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.

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.

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.

Prokaryotic reverse transcriptases: from retroelements to specialized defense systems

Reverse transcriptases (RTs) catalyze the polymerization of DNA from an RNA template. These enzymes were first discovered in RNA tumor viruses in 1970, but it was not until 1989 that they were found in prokaryotes as a key component of retrons. Apart from RTs encoded by the 'selfish' mobile retroelements known as group II introns, prokaryotic RTs are extraordinarily diverse, but their function has remained elusive. However, recent studies have revealed that different lineages of prokaryotic RTs, including retrons, those associated with CRISPR-Cas systems, Abi-like RTs and other yet uncharacterized RTs, are key components of different lines of defense against phages and other mobile genetic elements. Prokaryotic RTs participate in various antiviral strategies, including abortive infection (Abi), in which the infected cell is induced to commit suicide to protect the host population, adaptive immunity, in which a memory of previous infection is used to build an efficient defense, and other as yet unidentified mechanisms. These prokaryotic enzymes are attracting considerable attention, both for use in cutting-edge technologies, such as genome editing, and as an emerging research topic. In this review, we discuss what is known about prokaryotic RTs, and the exciting evidence for their domestication from retroelements to create specialized defense systems.

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.

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.

Recruitment of Reverse Transcriptase-Cas1 Fusion Proteins by Type VI-A CRISPR-Cas Systems

Type VI CRISPR-Cas systems contain a single effector nuclease (Cas13) that exclusively targets single-stranded RNA. It remains unknown how these systems acquire spacers. It has been suggested that type VI systems with adaptation modules can acquire spacers from RNA bacteriophages, but sequence similarities suggest that spacers may provide immunity to DNA phages. We searched databases for Cas13 proteins with linked RTs. We identified two different type VI-A systems with adaptation modules including an RT-Cas1 fusion and Cas2 proteins. Phylogenetic reconstruction analyses revealed that these adaptation modules were recruited by different effector Cas13a proteins, possibly from RT-associated type III-D systems within the bacterial classes Alphaproteobacteria and Clostridia. These type VI-A systems are predicted to acquire spacers from RNA molecules, paving the way for future studies investigating their role in bacterial adaptive immunity and biotechnological applications.

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.

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.

Metabarcoding reveals that rhizospheric microbiota of Quercus pyrenaica is composed by a relatively small number of bacterial taxa highly abundant

Melojo oak (Quercus pyrenaica Willd.) is a key tree species of Mediterranean forests; however, these forests show an advanced stage of deterioration in the Iberian Peninsula. Plant-associated microorganisms play an essential role improving their host’s fitness, hence, a better understanding of oak rhizospheric microbiome, especially of those active members, could be the first step towards microbiome-based approaches for oak-forest improvement. Here we reported, for the first time, the diversity of total (DNA-based) and potentially active (RNA-based) bacterial communities of different melojo-oak forest formations through pyrosequencing of 16S rRNA gene amplicons. We found that potentially active bacterial communities were as rich and diverse as total bacterial communities, but different in terms of relative abundance patterns in some of the studied areas. Both core microbiomes were dominated by a relatively small percentage of OTUs, most of which showed positive correlation between both libraries. However, the uncoupling between abundance (rDNA) and potential activity (rRNA) for some taxa suggests that the most abundant taxa are not always the most active, and that low-abundance OTUs may have a strong influence on oak’s rhizospheric ecology. Thus, measurement of rRNA:rDNA ratio could be helpful in identifying major players for the development of bacterial bioinoculants.

On the Origin and Evolutionary Relationships of the Reverse Transcriptases Associated With Type III CRISPR-Cas Systems

Reverse transcriptases (RTs) closely related to those encoded by group II introns but lacking the intron RNA structure have been found associated with type III clustered regularly interspaced short palindromic repeats (CRISPR)-Cas systems, a prokaryotic immune system against invading viruses and foreign genetic elements. Two models have been proposed to explain the origin and evolutionary relationships of these RTs: (i) the “single point of origin” model, according to which these RTs originated from a single acquisition event in bacterial, with the various protein domains (RT, RT-Cas1, and Cas6-RT-Cas1 fusions) corresponding to single points in evolution; and (ii) the “various origins” model, according to which, independent acquisition events in different evolutionary episodes led to these fusions. We tested these alternative hypotheses, by analyzing and integrating published datasets of RT sequences associated with CRISPR-Cas systems and inferring phylogenetic trees by maximum likelihood (ML) methods. The RTs studied could be grouped into 13 clades, mostly in bacteria, in which they probably evolved. The various clades appear to form three independent lineages in bacteria and a recent lineage in archaea. Our data show that the Cas6 domain was acquired twice, independently, through RT-Cas1 fusion, in the bacterial lineages. Taken together, there more evidence to support the “various origins” hypothesis.

Contribution of Mobile Group II Introns to Sinorhizobium meliloti Genome Evolution

Mobile group II introns are ribozymes and retroelements that probably originate from bacteria. Sinorhizobium meliloti, the nitrogen-fixing endosymbiont of legumes of genus Medicago, harbors a large number of these retroelements. One of these elements, RmInt1, has been particularly successful at colonizing this multipartite genome. Many studies have improved our understanding of RmInt1 and phylogenetically related group II introns, their mobility mechanisms, spread and dynamics within S. meliloti and closely related species. Although RmInt1 conserves the ancient retroelement behavior, its evolutionary history suggests that this group II intron has played a role in the short- and long-term evolution of the S. meliloti genome. We will discuss its proposed role in genome evolution by controlling the spread and coexistence of potentially harmful mobile genetic elements, by ectopic transposition to different genetic loci as a source of early genomic variation and by generating sequence variation after a very slow degradation process, through intron remnants that may have continued to evolve, contributing to bacterial speciation.

Functionality of In vitro Reconstituted Group II Intron RmInt1-Derived Ribonucleoprotein Particles

The functional unit of mobile group II introns is a ribonucleoprotein particle (RNP) consisting of the intron-encoded protein (IEP) and the excised intron RNA. The IEP has reverse transcriptase activity but also promotes RNA splicing, and the RNA-protein complex triggers site-specific DNA insertion by reverse splicing, in a process called retrohoming. In vitro reconstituted ribonucleoprotein complexes from the Lactococcus lactis group II intron Ll.LtrB, which produce a double strand break, have recently been studied as a means of developing group II intron-based gene targeting methods for higher organisms. The Sinorhizobium meliloti group II intron RmInt1 is an efficient mobile retroelement, the dispersal of which appears to be linked to transient single-stranded DNA during replication. The RmInt1IEP lacks the endonuclease domain (En) and cannot cut the bottom strand to generate the 3' end to initiate reverse transcription. We used an Escherichia coli expression system to produce soluble and active RmInt1 IEP and reconstituted RNPs with purified components in vitro. The RNPs generated were functional and reverse-spliced into a single-stranded DNA target. This work constitutes the starting point for the use of group II introns lacking DNA endonuclease domain-derived RNPs for highly specific gene targeting methods.

Host Factors Influencing the Retrohoming Pathway of Group II Intron RmInt1, Which Has an Intron-Encoded Protein Naturally Devoid of Endonuclease Activity

Bacterial group II introns are self-splicing catalytic RNAs and mobile retroelements that have an open reading frame encoding an intron-encoded protein (IEP) with reverse transcriptase (RT) and RNA splicing or maturase activity. Some IEPs carry a DNA endonuclease (En) domain, which is required to cleave the bottom strand downstream from the intron-insertion site for target DNA-primed reverse transcription (TPRT) of the inserted intron RNA. Host factors complete the insertion of the intron. By contrast, the major retrohoming pathway of introns with IEPs naturally lacking endonuclease activity, like the Sinorhizobium meliloti intron RmInt1, is thought to involve insertion of the intron RNA into the template for lagging strand DNA synthesis ahead of the replication fork, with possible use of the nascent strand to prime reverse transcription of the intron RNA. The host factors influencing the retrohoming pathway of such introns have not yet been described. Here, we identify key candidates likely to be involved in early and late steps of RmInt1 retrohoming. Some of these host factors are common to En⁺ group II intron retrohoming, but some have different functions. Our results also suggest that the retrohoming process of RmInt1 may be less dependent on the intracellular free Mg²⁺ concentration than those of other group II introns.

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.

Bacterial Group II Introns: Identification and Mobility Assay

Group II introns are large catalytic RNAs and mobile retroelements that encode a reverse transcriptase. Here, we provide methods for their identification in bacterial genomes and further analysis of their splicing and mobility capacities.

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.

Genomic characterization of Sinorhizobium meliloti AK21, a wild isolate from the Aral Sea Region

The symbiotic, nitrogen-fixing bacterium Sinorhizobium meliloti has been widely studied due to its ability to improve crop yields through direct interactions with leguminous plants. S. meliloti AK21 is a wild type strain that forms nodules on Medicago plants in saline and drought conditions in the Aral Sea Region. The aim of this work was to establish the genetic similarities and differences between S. meliloti AK21 and the reference strain S. meliloti 1021. Comparative genome hybridization with the model reference strain S. meliloti 1021 yielded 365 variable genes, grouped into 11 regions in the three main replicons in S. meliloti AK21. The most extensive regions of variability were found in the symbiotic plasmid pSymA, which also contained the largest number of orthologous and polymorphic sequences identified by suppression subtractive hybridization. This procedure identified a large number of divergent sequences and others without homology in the databases, the further investigation of which could provide new insight into the alternative metabolic pathways present in S. meliloti AK21. We identified a plasmid replication module from the repABC replicon family, together with plasmid mobilization-related genes (traG and a VirB9-like protein), which suggest that this indigenous isolate harbors an accessory plasmid. Furthermore, the transcriptomic profiles reflected differences in gene content and regulation between S. meliloti AK21 and S. meliloti 1021 (ExpR and PhoB regulons), but provided evidence for an as yet unknown, alternative mechanism involving activation of the cbb3 terminal oxidase. Finally, phenotypic microarrays characterization revealed a greater versatility of substrate use and chemical degradation than for S. meliloti 1021.

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.

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.

Spread of the group II intron RmInt1 and its insertion sequence target sites in the plant endosymbiont Sinorhizobium meliloti

RmInt1 is a mobile group II intron from Sinorhizobium meliloti that is exceptionally abundant in this bacterial species. We compared the presence of RmInt1 and two of its insertion sequence homing sites (ISRm2011-2 and ISRm10-2) in two phylogenetic clusters (I and II) identified by AFLP analysis in a collection of S. meliloti field isolates from Italy. Both clusters contained several copies of the ISRm2011-2 element, which is present at high copy number in almost all S. meliloti isolates. By contrast, isolates from cluster I harbored no copies of ISRm10-2 and only a truncated copy of RmInt1, despite the absence of constraints on intron mobility in this genetic background, whereas cluster II strains harbored several copies of this intron. The absence of ISRm10-2 from one of the strains of this cluster suggests that this element was acquired more recently than the other two elements. Furthermore, studies of insertional polymorphisms in cluster II strains revealed the acquisition of ISRm10-2 and subsequent retrohoming of RmInt1 to this homing site. These results highlight the role of intron homing sites (ISs) in facilitating intron dispersal and the dynamic and ongoing nature of the spread of the group II intron RmInt1 in S. meliloti.

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.

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.

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.

Localization of a bacterial group II intron-encoded protein in eukaryotic nuclear splicing-related cell compartments

Some bacterial group II introns are widely used for genetic engineering in bacteria, because they can be reprogrammed to insert into the desired DNA target sites. There is considerable interest in developing this group II intron gene targeting technology for use in eukaryotes, but nuclear genomes present several obstacles to the use of this approach. The nuclear genomes of eukaryotes do not contain group II introns, but these introns are thought to have been the progenitors of nuclear spliceosomal introns. We investigated the expression and subcellular localization of the bacterial RmInt1 group II intron-encoded protein (IEP) in Arabidopsis thaliana protoplasts. Following the expression of translational fusions of the wild-type protein and several mutant variants with EGFP, the full-length IEP was found exclusively in the nucleolus, whereas the maturase domain alone targeted EGFP to nuclear speckles. The distribution of the bacterial RmInt1 IEP in plant cell protoplasts suggests that the compartmentalization of eukaryotic cells into nucleus and cytoplasm does not prevent group II introns from invading the host genome. Furthermore, the trafficking of the IEP between the nucleolus and the speckles upon maturase inactivation is consistent with the hypothesis that the spliceosomal machinery evolved from group II introns.

Independent activity of the homologous small regulatory RNAs AbcR1 and AbcR2 in the legume symbiont Sinorhizobium meliloti

The legume symbiont Sinorhizobium meliloti expresses a plethora of small noncoding RNAs (sRNAs) whose function is mostly unknown. Here, we have functionally characterized two tandemly encoded S. meliloti Rm1021 sRNAs that are similar in sequence and structure. Homologous sRNAs (designated AbcR1 and AbcR2) have been shown to regulate several ABC transporters in the related α-proteobacteria Agrobacterium tumefaciens and Brucella abortus. In Rm1021, AbcR1 and AbcR2 exhibit divergent unlinked regulation and are stabilized by the RNA chaperone Hfq. AbcR1 is transcribed in actively dividing bacteria, either in culture, rhizosphere or within the invasion zone of mature alfalfa nodules. Conversely, AbcR2 expression is induced upon entry into stationary phase and under abiotic stress. Only deletion of AbcR1 resulted into a discrete growth delay in rich medium, but both are dispensable for symbiosis. Periplasmic proteome profiling revealed down-regulation of the branched-chain amino acid binding protein LivK by AbcR1, but not by AbcR2. A double-plasmid reporter assay confirmed the predicted specific targeting of the 5′-untranslated region of the livK mRNA by AbcR1 in vivo. Our findings provide evidences of independent regulatory functions of these sRNAs, probably to fine-tune nutrient uptake in free-living and undifferentiated symbiotic rhizobia.

Comprehensive phylogenetic analysis of bacterial group II intron-encoded ORFs lacking the DNA endonuclease domain reveals new varieties

Group II introns are self-splicing RNAs that act as mobile retroelements in the organelles of plants, fungi and protists. They are also widely distributed in bacteria, and are generally assumed to be the ancestors of nuclear spliceosomal introns. Most bacterial group II introns have a multifunctional intron-encoded protein (IEP) ORF within the ribozyme domain IV (DIV). This ORF encodes an N-terminal reverse transcriptase (RT) domain, followed by a putative RNA-binding domain with RNA splicing or maturase activity and, in some cases, a C-terminal DNA-binding (D) region followed by a DNA endonuclease (En) domain. In this study, we focused on bacterial group II intron ORF phylogenetic classes containing only reverse transcriptase/maturase open reading frames, with no recognizable D/En region (classes A, C, D, E, F and unclassified introns). On the basis of phylogenetic analyses of the maturase domain and its C-terminal extension, which appears to be a signature characteristic of ORF phylogenetic class, with support from the phylogeny inferred from the RT domain, we have revised the proposed new class F, defining new intron ORF varieties. Our results increase knowledge of the lineage of group II introns encoding proteins lacking the En-domain.

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.

Search for a new gauge boson in electron-nucleus fixed-target scattering by the APEX experiment

We present a search at the Jefferson Laboratory for new forces mediated by sub-GeV vector bosons with weak coupling α' to electrons. Such a particle A' can be produced in electron-nucleus fixed-target scattering and then decay to an e + e- pair, producing a narrow resonance in the QED trident spectrum. Using APEX test run data, we searched in the mass range 175-250 MeV, found no evidence for an A'→ e+ e- reaction, and set an upper limit of α'/α ~/= 10(-6). Our findings demonstrate that fixed-target searches can explore a new, wide, and important range of masses and couplings for sub-GeV forces.

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.

Characterization of novel antibiotic resistance genes identified by functional metagenomics on soil samples

The soil microbial community is highly complex and contains a high density of antibiotic-producing bacteria, making it a likely source of diverse antibiotic resistance determinants. We used functional metagenomics to search for antibiotic resistance genes in libraries generated from three different soil samples, containing 3.6 Gb of DNA in total. We identified 11 new antibiotic resistance genes: 3 conferring resistance to ampicillin, 2 to gentamicin, 2 to chloramphenicol and 4 to trimethoprim. One of the clones identified was a new trimethoprim resistance gene encoding a 26.8 kDa protein closely resembling unassigned reductases of the dihydrofolate reductase group. This protein, Tm8-3, conferred trimethoprim resistance in Escherichia coli and Sinorhizobium meliloti (γ- and α-proteobacteria respectively). We demonstrated that this gene encoded an enzyme with dihydrofolate reductase activity, with kinetic constants similar to other type I and II dihydrofolate reductases (K(m) of 8.9 µM for NADPH and 3.7 µM for dihydrofolate and IC(50) of 20 µM for trimethoprim). This is the first description of a new type of reductase conferring resistance to trimethoprim. Our results indicate that soil bacteria display a high level of genetic diversity and are a reservoir of antibiotic resistance genes, supporting the use of this approach for the discovery of novel enzymes with unexpected activities unpredictable from their amino acid sequences.

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.

Analysis of Rhizobium meliloti Sym Mutants Obtained by Heat Treatment

Deletions in the pSym megaplasmid of Rhizobium meliloti were produced at a high frequency, and their lengths varied according to incubation temperature. Morphological differentiation into large and small colonies occurred after heat treatment. Small colonies elicited pseudonodules on alfalfa roots.

The Sinorhizobium meliloti RNA chaperone Hfq influences central carbon metabolism and the symbiotic interaction with alfalfa

Background

The bacterial Hfq protein is able to interact with diverse RNA molecules, including regulatory small non-coding RNAs (sRNAs), and thus it is recognized as a global post-transcriptional regulator of gene expression. Loss of Hfq has an extensive impact in bacterial physiology which in several animal pathogens influences virulence. Sinorhizobium meliloti is a model soil bacterium known for its ability to establish a beneficial nitrogen-fixing intracellular symbiosis with alfalfa. Despite the predicted general involvement of Hfq in the establishment of successful bacteria-eukaryote interactions, its function in S. meliloti has remained unexplored.

Results

Two independent S. meliloti mutants, 2011-3.4 and 1021Δhfq, were obtained by disruption and deletion of the hfq gene in the wild-type strains 2011 and 1021, respectively, both exhibiting similar growth defects as free-living bacteria. Transcriptomic profiling of 1021Δhfq revealed a general down-regulation of genes of sugar transporters and some enzymes of the central carbon metabolism, whereas transcripts specifying the uptake and metabolism of nitrogen sources (mainly amino acids) were more abundant than in the wild-type strain. Proteomic analysis of the 2011-3.4 mutant independently confirmed these observations. Symbiotic tests showed that lack of Hfq led to a delayed nodulation, severely compromised bacterial competitiveness on alfalfa roots and impaired normal plant growth. Furthermore, a large proportion of nodules (55%-64%) elicited by the 1021Δhfq mutant were non-fixing, with scarce content in bacteroids and signs of premature senescence of endosymbiotic bacteria. RT-PCR experiments on RNA from bacteria grown under aerobic and microoxic conditions revealed that Hfq contributes to regulation of nifA and fixK1/K2, the genes controlling nitrogen fixation, although the Hfq-mediated regulation of fixK is only aerobiosis dependent. Finally, we found that some of the recently identified S. meliloti sRNAs co-inmunoprecipitate with a FLAG-epitope tagged Hfq protein.

Conclusions

Our results support that the S. meliloti RNA chaperone Hfq contributes to the control of central metabolic pathways in free-living bacteria and influences rhizospheric competence, survival of the microsymbiont within the nodule cells and nitrogen fixation during the symbiotic interaction with its legume host alfalfa. The identified S. meliloti Hfq-binding sRNAs are predicted to participate in the Hfq regulatory network.

Structural features in the C-terminal region of the Sinorhizobium meliloti RmInt1 group II intron-encoded protein contribute to its maturase and intron DNA-insertion function

Group II introns are both catalytic RNAs and mobile retroelements that move through a process catalyzed by a RNP complex consisting of an intron-encoded protein and the spliced intron lariat RNA. Group II intron-encoded proteins are multifunctional and contain an N-terminal reverse transcriptase domain, followed by a putative RNA-binding domain (domain X) associated with RNA splicing or maturase activity and a C-terminal DNA binding/DNA endonuclease region. The intron-encoded protein encoded by the mobile group II intron RmInt1, which lacks the DNA binding/DNA endonuclease region, has only a short C-terminal extension (C-tail) after a typical domain X, apparently unrelated to the C-terminal regions of other group II intron-encoded proteins. Multiple sequence alignments identified features of the C-terminal portion of the RmInt1 intron-encoded protein that are conserved throughout evolution in the bacterial ORF class D, suggesting a group-specific functionally important protein region. The functional importance of these features was demonstrated by analyses of deletions and mutations affecting conserved amino acid residues. We found that the C-tail of the RmInt1 intron-encoded protein contributes to the maturase function of this reverse transcriptase protein. Furthermore, within the C-terminal region, we identified, in a predicted alpha-helical region and downstream, conserved residues that are specifically required for the insertion of the intron into DNA targets in the orientation that would make it possible to use the nascent leading strand as a primer. These findings suggest that these group II intron intron-encoded proteins may have adapted to function in mobility by different mechanisms to make use of either leading or lagging-oriented targets in the absence of an endonuclease domain.

Identification of differentially expressed small non-coding RNAs in the legume endosymbiont Sinorhizobium meliloti by comparative genomics

Bacterial small non-coding RNAs (sRNAs) are being recognized as novel widespread regulators of gene expression in response to environmental signals. Here, we present the first search for sRNA-encoding genes in the nitrogen-fixing endosymbiont Sinorhizobium meliloti, performed by a genome-wide computational analysis of its intergenic regions. Comparative sequence data from eight related α-proteobacteria were obtained, and the interspecies pairwise alignments were scored with the programs eQRNA and RNAz as complementary predictive tools to identify conserved and stable secondary structures corresponding to putative non-coding RNAs. Northern experiments confirmed that eight of the predicted loci, selected among the original 32 candidates as most probable sRNA genes, expressed small transcripts. This result supports the combined use of eQRNA and RNAz as a robust strategy to identify novel sRNAs in bacteria. Furthermore, seven of the transcripts accumulated differentially in free-living and symbiotic conditions. Experimental mapping of the 5′-ends of the detected transcripts revealed that their encoding genes are organized in autonomous transcription units with recognizable promoter and, in most cases, termination signatures. These findings suggest novel regulatory functions for sRNAs related to the interactions of α-proteobacteria with their eukaryotic hosts.

Bacterial group II introns: not just splicing

Group II introns are both catalytic RNAs (ribozymes) and mobile retroelements that were discovered almost 14 years ago. It has been suggested that eukaryotic mRNA introns might have originated from the group II introns present in the alphaproteobacterial progenitor of the mitochondria. Bacterial group II introns are of considerable interest not only because of their evolutionary significance, but also because they could potentially be used as tools for genetic manipulation in biotechnology and for gene therapy. This review summarizes what is known about the splicing mechanisms and mobility of bacterial group II introns, and describes the recent development of group II intron-based gene-targetting methods. Bacterial group II intron diversity, evolutionary relationships, and behaviour in bacteria are also discussed.

Identification and characterization of a nodH ortholog from the alfalfa-nodulating Or191-like rhizobia

Nodulation of Medicago sativa (alfalfa) is known to be restricted to Sinorhizobium meliloti and a few other rhizobia that include the poorly characterized isolates related to Rhizobium sp. strain Or191. Distinctive features of the symbiosis between alfalfa and S. meliloti are the marked specificity from the plant to the bacteria and the strict requirement for the presence of sulfated lipochitooligosaccharides (Nod factors [NFs]) at its reducing end. Here, we present evidence of the presence of a functional nodH-encoded NF sulfotransferase in the Or191-like rhizobia. The nodH gene, present in single copy, maps to a high molecular weight megaplasmid. As in S. meliloti, a nodF homolog was identified immediately upstream of nodH that was transcribed in the opposite direction (local synteny). This novel nodH ortholog was cloned and shown to restore both NF sulfation and the Nif+Fix+ phenotypes when introduced into an S. meliloti nodH mutant. Unexpectedly, however, nodH disruption in the Or191-like bacteria did not abolish their ability to nodulate alfalfa, resulting instead in a severely delayed nodulation. In agreement with evidence from other authors, the nodH sequence analysis strongly supports the idea that the Or191-like rhizobia most likely represent a genetic mosaic resulting from the horizontal transfer of symbiotic genes from a sinorhizobial megaplasmid to a not yet clearly identified ancestor.

Rhizosphere-bacterial community in Eperua falcata (Caesalpiniaceae) a putative nitrogen-fixing tree from French Guiana rainforest

The rainforest of French Guiana is still largely unaffected by human activity. Various pristine sites like the Paracou Research Station are devoted to study this tropical ecosystem. We used culture-independent techniques, like polymerase chain reaction-temperature gradient gel electrophoresis, and construction of clone libraries of partial 16S rRNA and nifH genes, to analyze the composition of the bacterial community in the rhizosphere of mature trees of Eperua falcata and Dicorynia guianensis, both species within the Caesalpiniaceae family. E. falcata is one of the more abundant pioneer tree species in this ecosystem and so far, no root nodules have ever been found. However, its nitrogen-fixing status is regarded as "uncertain", whereas D. guianensis is clearly considered a non-nitrogen-fixing plant. The rhizospheres of these mature trees contain specific bacterial communities, including several currently found uncultured microorganisms. In these communities, there are putative nitrogen-fixing bacteria specifically associated to each tree: D. guianensis harbors several Rhizobium spp. and E. falcata members of the genera Burkholderia and Bradyrhizobium. In addition, nifH sequences in the rhizosphere of the latter tree were very diverse. Retrieved sequences were related to bacteria belonging to the alpha-, beta-, and gamma-Proteobacteria in the E. falcata rhizoplane, whereas only two sequences related to gamma-Proteobacteria were found in D. guianensis. Differences in the bacterial communities and the abundance and diversity of nifH sequences in E. falcata rhizosphere suggest that this tree could obtain nitrogen through a nonnodulating bacterial interaction.

Dispersion of the RmInt1 group II intron in the Sinorhizobium meliloti genome upon acquisition by conjugative transfer

RmInt1 is a self-splicing and mobile group II intron initially identified in the bacterium Sinorhizobium meliloti, which encodes a reverse transcriptase–maturase (Intron Encoded Protein, IEP) lacking the C-terminal DNA binding (D) and DNA endonuclease domains (En). RmInt1 invades cognate intronless homing sites (ISRm2011-2) by a mechanism known as retrohoming. This work describes how the RmInt1 intron spreads in the S.meliloti genome upon acquisition by conjugation. This process was revealed by using the wild-type intron RmInt1 and engineered intron-donor constructs based on ribozyme coding sequence (ΔORF)-derivatives with higher homing efficiency than the wild-type intron. The data demonstrate that RmInt1 propagates into the S.meliloti genome primarily by retrohoming with a strand bias related to replication of the chromosome and symbiotic megaplasmids. Moreover, we show that when expressed in trans from a separate plasmid, the IEP is able to mobilize genomic ΔORF ribozymes that afterward displayed wild-type levels of retrohoming. Our results contribute to get further understanding of how group II introns spread into bacterial genomes in nature.

Excision of the Sinorhizobium meliloti group II intron RmInt1 as circles in vivo

Excision of group II introns as circles has been described only for a few eukaryotic introns and little is known about the mechanisms involved, the relevance or consequences of the process. We report that splicing of the bacterial group II intron RmInt1 in vivo leads to the formation of both intron lariat and intron RNA circles. We determined that besides being required for the intron splicing reaction, the maturase domain of the intron-encoded protein also controls the balance between lariat and RNA intron circle production. Furthermore, comparison with in vitro self-splicing products indicates that in vivo, the intron-encoded protein appears to promote the use of a correct EBS1/IBS1 intron-exon interaction as well as cleavage at, or next to, the expected 3' splice site. These findings provide new insights on the mechanism of excision of group II introns as circles.

Potential for alternative intron-exon pairings in group II intron RmInt1 from Sinorhizobium meliloti and its relatives

Ribozyme constructs derived from group II intron RmInt1 of Sinorhizobium meliloti self-splice in vitro when incubated under permissive conditions, but exon ligation is unusually inefficient when the 5' exon is truncated close to the IBS2 intron-binding site. One plausible explanation for this observation is the presence of an alternative intron-exon pairing between an intron segment that overlaps with the EBS2 exon-binding site and a 5' exon site located just distal of IBS2 relative to the splice junction. Strikingly, the existence of this pairing is supported by comparative sequence analysis of introns related to RmInt1.

An alternative intron-exon pairing scheme implied by unexpected in vitro activities of group II intron RmInt1 from Sinorhizobium meliloti

RmInt1 is a mobile group II intron which interrupts ISRm2011-2, another mobile element from the bacterium Sinorhizobium meliloti. Ribozyme constructs derived from intron RmInt1 self-splice in vitro when incubated under permissive conditions, but the excised intron and ligated exons are largely replaced by unconventional products. These include a slightly shorter, 5'-end truncated 3' exon, truncated variants of the linear and lariat forms of the intron-3' exon reaction intermediate, as well as presumably circular molecules derived from the latter. Two factors explain the abundance of these products: (i) nucleotides 5-11 of the 3' exon (IBS1*) provide a better match to the EBS1 5'-exon-binding site than the authentic IBS1 sequence in the 5' exon; (ii) exon ligation is unusually inefficient, and especially so when the 5' exon is truncated close to the second (IBS2) intron-binding site. We propose that reactions at the IBS1* site play a part in the regulation of the intron ISRm2011-2 host in vivo.

The coexistence of symbiosis and pathogenicity-determining genes in Rhizobium rhizogenes strains enables them to induce nodules and tumors or hairy roots in plants

Bacteria belonging to the family Rhizobiaceae may establish beneficial or harmful relationships with plants. The legume endosymbionts contain nod and nif genes responsible for nodule formation and nitrogen fixation, respectively, whereas the pathogenic strains carry vir genes responsible for the formation of tumors or hairy roots. The symbiotic and pathogenic strains currently belong to different species of the genus Rhizobium and, until now, no strains able to establish symbiosis with legumes and also to induce tumors or hairy roots in plants have been reported. Here, we report for the first time the occurrence of two rhizobial strains (163C and ATCC11325T) belonging to Rhizobium rhizogenes able to induce hairy roots or tumors in plants and also to nodulate Phaseolus vulgaris under natural environmental conditions. Symbiotic plasmids (pSym) containing nod and nif genes and pTi- or pRi-type plasmids containing vir genes were found in these strains. The nodD and nifH genes of the strains from this study are phylogenetically related to those of Sinorhizobium strains nodulating P. vulgaris. The virA and virB4 genes from strain 163C are phylogenetically related to those of R. tumefaciens C58, whereas the same genes from strain ATCC 11325T are related to those of hairy root-inducing strains. These findings may be of high relevance for the better understanding of plant-microbe interactions and knowledge of rhizobial phylogenetic history.

Dispersal and evolution of the Sinorhizobium meliloti group II RmInt1 intron in bacteria that interact with plants

Group II introns are both self-splicing RNAs and mobile retroelements found in bacterial and archaeal genomes and in organelles of eukaryotes. They are thought to be the ancestors of eukaryote spliceosomal introns and non-long terminal repeat retrotransposons. We show here that RmInt1, a bacterial group II intron first described in the nitrogen-fixing symbiont of alfalfa (Medicago sativa) Sinorhizobium meliloti, is also present in other Sinorhizobium and Rhizobium species. The intron-homing sites in these species are IS elements of the ISRm2011-2 group as in S. meliloti, but ectopic insertion is also observed. We present evidence that these related bacteria have acquired RmInt1 by vertical inheritance from a common ancestor and by independent horizontal transfer events. We also show that RmInt1 is mobile in related taxa of bacteria that interact with plants and tends to evolve toward an inactive form by fragmentation, with loss of the 3' terminus including the intron-encoded protein. Our results provide an overview of the evolution and dispersion of a bacterial group II intron.

The RmInt1 group II intron has two different retrohoming pathways for mobility using predominantly the nascent lagging strand at DNA replication forks for priming

Sinorhizobium meliloti RmInt1 is an efficient mobile group II intron that uses an unknown reverse transcriptase priming mechanism as the intron ribonucleoprotein complex can reverse splice into DNA target substrates but cannot carry out site-specific second strand cleavage due to the lack of a C-terminal DNA endonuclease domain. We show here that, like other mobile group II introns, RmInt1 moves around by an efficient RNA-based retrohoming mechanism. We found evidence of two distinct RmInt1 retrohoming pathways for mobility depending on the orientation of the target site relative to the direction of DNA replication. The preferred retrohoming pathway is consistent with reverse splicing of the intron RNA into single-stranded DNA at a replication fork, using a nascent lagging DNA strand as the primer for reverse transcription. This strand bias is the opposite of that reported for mobility of the lactococcal Ll.ltrB intron in the absence of second strand cleavage. The mobility mechanism found here for RmInt1 may be used for dissemination by many bacterial group II introns encoding proteins lacking the DNA endonuclease domain.