Root Primordium Defective 1 Encodes an Essential PORR Protein Required for the Splicing of Mitochondria-Encoded Group II Introns and for Respiratory Complex I Biogenesis

Cellular respiration involves complex organellar metabolic activities that are pivotal for plant growth and development. Mitochondria contain their own genetic system (mitogenome, mtDNA), which encodes key elements of the respiratory machinery. Plant mtDNAs are notably larger than their counterparts in Animalia, with complex genome organization and gene expression characteristics. The maturation of the plant mitochondrial transcripts involves extensive RNA editing, trimming and splicing events. These essential processing steps rely on the activities of numerous nuclear-encoded cofactors, which may also play key regulatory roles in mitochondrial biogenesis and function and hence in plant physiology. Proteins that harbor the plant organelle RNA recognition (PORR) domain are represented in a small gene family in plants. Several PORR members, including WTF1, WTF9 and LEFKOTHEA, are known to act in the splicing of organellar group II introns in angiosperms. The AT4G33495 gene locus encodes an essential PORR protein in Arabidopsis, termed ROOT PRIMORDIUM DEFECTIVE 1 (RPD1). A null mutation of At.RPD1 causes arrest in early embryogenesis, while the missense mutant lines, rpd1.1 and rpd1.2, exhibit a strong impairment in root development and retarded growth phenotypes, especially under high-temperature conditions. Here, we further show that RPD1 functions in the splicing of introns that reside in the coding regions of various complex I (CI) subunits (i.e. nad2, nad4, nad5 and nad7), as well as in the maturation of the ribosomal rps3 pre-RNA in Arabidopsis mitochondria. The altered growth and developmental phenotypes and modified respiration activities are tightly correlated with respiratory chain CI defects in rpd1 mutants.

Categorizing 161 plant (streptophyte) mitochondrial group II introns into 29 families of related paralogues finds only limited links between intron mobility and intron-borne maturases

Group II introns are common in the two endosymbiotic organelle genomes of the plant lineage. Chloroplasts harbor 22 positionally conserved group II introns whereas their occurrence in land plant (embryophyte) mitogenomes is highly variable and specific for the seven major clades: liverworts, mosses, hornworts, lycophytes, ferns, gymnosperms and flowering plants. Each plant group features "signature selections" of ca. 20-30 paralogues from a superset of altogether 105 group II introns meantime identified in embryophyte mtDNAs, suggesting massive intron gains and losses along the backbone of plant phylogeny. We report on systematically categorizing plant mitochondrial group II introns into "families", comprising evidently related paralogues at different insertion sites, which may even be more similar than their respective orthologues in phylogenetically distant taxa. Including streptophyte (charophyte) algae extends our sampling to 161 and we sort 104 streptophyte mitochondrial group II introns into 25 core families of related paralogues evidently arising from retrotransposition events. Adding to discoveries of only recently created intron paralogues, hypermobile introns and twintrons, our survey led to further discoveries including previously overlooked "fossil" introns in spacer regions or e.g., in the rps8 pseudogene of lycophytes. Initially excluding intron-borne maturase sequences for family categorization, we added an independent analysis of maturase phylogenies and find a surprising incongruence between intron mobility and the presence of intron-borne maturases. Intriguingly, however, we find that several examples of nuclear splicing factors meantime characterized simultaneously facilitate splicing of independent paralogues now placed into the same intron families. Altogether this suggests that plant group II intron mobility, in contrast to their bacterial counterparts, is not intimately linked to intron-encoded maturases.

Group II intron and repeat-rich red algal mitochondrial genomes demonstrate the dynamic recent history of autocatalytic RNAs

Background

Group II introns are mobile genetic elements that can insert at specific target sequences, however, their origins are often challenging to reconstruct because of rapid sequence decay following invasion and spread into different sites. To advance understanding of group II intron spread, we studied the intron-rich mitochondrial genome (mitogenome) in the unicellular red alga, Porphyridium.

Results

Analysis of mitogenomes in three closely related species in this genus revealed they were 3–6-fold larger in size (56–132 kbp) than in other red algae, that have genomes of size 21–43 kbp. This discrepancy is explained by two factors, group II intron invasion and expansion of repeated sequences in large intergenic regions. Phylogenetic analysis demonstrates that many mitogenome group II intron families are specific to Porphyridium, whereas others are closely related to sequences in fungi and in the red alga-derived plastids of stramenopiles. Network analysis of intron-encoded proteins (IEPs) shows a clear link between plastid and mitochondrial IEPs in distantly related species, with both groups associated with prokaryotic sequences.

Conclusion

Our analysis of group II introns in Porphyridium mitogenomes demonstrates the dynamic nature of group II intron evolution, strongly supports the lateral movement of group II introns among diverse eukaryotes, and reveals their ability to proliferate, once integrated in mitochondrial DNA.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01200-3.

Easy identification of insertion sequence mobilization events in related bacterial strains with ISCompare

Bacterial genomes are composed of core and accessory genomes. The first is composed of housekeeping and essential genes, while the second is highly enriched in mobile genetic elements, including transposable elements (TEs). Insertion sequences (ISs), the smallest TEs, have an important role in genome evolution, and contribute to bacterial genome plasticity and adaptability. ISs can spread in a genome, presenting different locations in nearly related strains, and producing phenotypic variations. Few tools are available which can identify differentially located ISs (DLISs) on assembled genomes. Here, we introduce ISCompare, a new program to profile IS mobilization events in related bacterial strains using complete or draft genome assemblies. ISCompare was validated using artificial genomes with simulated random IS insertions and real sequences, achieving the same or better results than other available tools, with the advantage that ISCompare can analyze multiple ISs at the same time and outputs a list of candidate DLISs. ISCompare provides an easy and straightforward approach to look for differentially located ISs on bacterial genomes.

The association of group IIB intron with integrons in hypersaline environments

Background

Group II introns are mobile genetic elements used as efficient gene targeting tools. They function as both ribozymes and retroelements. Group IIC introns are the only class reported so far to be associated with integrons. In order to identify group II introns linked with integrons and CALINS (cluster of attC sites lacking a neighboring integron integrase) within halophiles, we mined for integrons in 28 assembled metagenomes from hypersaline environments and publically available 104 halophilic genomes using Integron Finder followed by blast search for group II intron reverse transcriptases (RT)s.

Results

We report the presence of different group II introns associated with integrons and integron-related sequences denoted by UHB.F1, UHB.I2, H.ha.F1 and H.ha.F2. The first two were identified within putative integrons in the metagenome of Tanatar-5 hypersaline soda lake, belonging to IIC and IIB intron classes, respectively at which the first was a truncated intron. Other truncated introns H.ha.F1 and H.ha.F2 were also detected in a CALIN within the extreme halophile Halorhodospira halochloris, both belonging to group IIB introns. The intron-encoded proteins (IEP) s identified within group IIB introns belonged to different classes: CL1 class in UHB.I2 and bacterial class E in H.ha.Fa1 and H.ha.F2. A newly identified insertion sequence (ISHahl1) of IS200/605 superfamily was also identified adjacent to H. halochloris CALIN. Finally, an abundance of toxin-antitoxin (TA) systems was observed within the identified integrons.

Conclusion

So far, this is the first investigation of group II introns within integrons in halophilic genomes and metagenomes from hypersaline environments. We report the presence of group IIB introns associated with integrons or CALINs. This study provides the basis for understanding the role of group IIB introns in the evolution of halophiles and their potential biotechnological role.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-021-00234-2.

PPR-SMR1 is required for the splicing of multiple mitochondrial introns, interacts with Zm-mCSF1, and is essential for seed development in maize

Group II introns are ribozymes that can excise themselves from precursor-RNA transcripts, but plant organellar group II introns have structural deviations that inhibit ribozyme activity. Therefore, splicing of these introns requires the assistance of nuclear- and/or organellar-encoded splicing factors; however, how these splicing factors function remains unclear. In this study, we report the functions and interactions of two splicing factors, PPR-SMR1 and Zm-mCSF1, in intron splicing in maize mitochondria. PPR-SMR1 is a SMR domain-containing pentatricopeptide repeat (PPR) protein and Zm-mCSF1 is a CRM domain-containing protein, and both are targeted to mitochondria. Loss-of-function mutations in each of them severely arrests embryogenesis and endosperm development in maize. Functional analyses indicate that PPR-SMR1 and Zm-mCSF1 are required for the splicing of most mitochondrial group II introns. Among them, nad2-intron 2 and 3, and nad5-intron 1 are PPR-SMR1/Zm-mCSF1-dependent introns. Protein interaction assays suggest that PPR-SMR1 can interact with Zm-mCSF1 through its N-terminus, and that Zm-mCSF1 is self-interacting. Our findings suggest that PPR-SMR1, a novel splicing factor, acts in the splicing of multiple group II introns in maize mitochondria, and the protein-protein interaction between it and Zm-mCSF1 might allow the formation of large macromolecular splicing complexes.

Diversity of Bacillus cereus sensu lato mobilome

BACKGROUND

Bacillus cereus sensu lato s.l.) is a group of bacteria displaying close phylogenetic relationships but a high ecological diversity. The three most studied species are Bacillus anthracis, Bacillus cereus sensu stricto and Bacillus thuringiensis. While some species are pathogenic to mammals or associated with food poisoning, Bacillus thuringiensis is a well-known entomopathogenic bacterium used as biopesticide worldwide. B. cereus s.l. also contains a large variety of mobile genetic elements (MGEs).

RESULTS

In this study, we detail the occurrence and plasmid vs. chromosome distribution of several MGEs in 102 complete and annotated genomes of B. cereus s.l. These MGEs include 16 Insertion Sequence (IS) families, the Tn3 family, 18 different Bacillus cereus repeats (BCRs) and 30 known group II introns.

CONCLUSIONS

Our analysis not only shows the diversity of these MGEs among strains of the same species and between different species within the B. cereus s.l. group, but also highlights the potential impact of these elements on the plasticity of the plasmid pool, and the TEs (Transposable Elements) - species relationship within B. cereus s.l.

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.

Insertion of Group II Intron-Based Ribozyme Switches into Homing Endonuclease Genes

Fungal mitochondrial genomes act as "reservoirs" for homing endonucleases. These enzymes with their DNA site-specific cleavage activities are attractive tools for genome editing, targeted mutagenesis and gene therapy applications. Herein, we present strategies where homing endonuclease open reading frames (HEases ORFs) are interrupted with group II intron sequences. The ultimate goal is to achieve in vivo expression of HEases that can be regulated by manipulating the splicing efficiency of the HEase ORF-embedded group II introns. That addition of exogenous magnesium chloride (MgCl₂) appears to stimulate splicing of nonnative group II introns in Escherichia coli and the addition of cobalt chloride (CoCl₂) to the growth medium antagonizes the expression of HEase activity (i.e., splicing). Group II introns are potentially autocatalytic self-splicing elements and thus can be used as molecular switches that allow for temporal regulated HEase expression. This should be useful in precision genome engineering, mutagenesis, and minimizing off-target activities.

Proliferation of group II introns in the chloroplast genome of the green alga Oedocladium carolinianum (Chlorophyceae)

Background

The chloroplast genome sustained extensive changes in architecture during the evolution of the Chlorophyceae, a morphologically and ecologically diverse class of green algae belonging to the Chlorophyta; however, the forces driving these changes are poorly understood. The five orders recognized in the Chlorophyceae form two major clades: the CS clade consisting of the Chlamydomonadales and Sphaeropleales, and the OCC clade consisting of the Oedogoniales, Chaetophorales, and Chaetopeltidales. In the OCC clade, considerable variations in chloroplast DNA (cpDNA) structure, size, gene order, and intron content have been observed. The large inverted repeat (IR), an ancestral feature characteristic of most green plants, is present in Oedogonium cardiacum (Oedogoniales) but is lacking in the examined members of the Chaetophorales and Chaetopeltidales. Remarkably, the Oedogonium 35.5-kb IR houses genes that were putatively acquired through horizontal DNA transfer. To better understand the dynamics of chloroplast genome evolution in the Oedogoniales, we analyzed the cpDNA of a second representative of this order, Oedocladium carolinianum.

Methods

The Oedocladium cpDNA was sequenced and annotated. The evolutionary distances separating Oedocladium and Oedogonium cpDNAs and two other pairs of chlorophycean cpDNAs were estimated using a 61-gene data set. Phylogenetic analysis of an alignment of group IIA introns from members of the OCC clade was performed. Secondary structures and insertion sites of oedogonialean group IIA introns were analyzed.

Results

The 204,438-bp Oedocladium genome is 7.9 kb larger than the Oedogonium genome, but its repertoire of conserved genes is remarkably similar and gene order differs by only one reversal. Although the 23.7-kb IR is missing the putative foreign genes found in Oedogonium, it contains sequences coding for a putative phage or bacterial DNA primase and a hypothetical protein. Intergenic sequences are 1.5-fold longer and dispersed repeats are more abundant, but a smaller fraction of the Oedocladium genome is occupied by introns. Six additional group II introns are present, five of which lack ORFs and carry highly similar sequences to that of the ORF-less IIA intron shared with Oedogonium. Secondary structure analysis of the group IIA introns disclosed marked differences in the exon-binding sites; however, each intron showed perfect or nearly perfect base pairing interactions with its target site.

Discussion

Our results suggest that chloroplast genes rearrange more slowly in the Oedogoniales than in the Chaetophorales and raise questions as to what was the nature of the foreign coding sequences in the IR of the common ancestor of the Oedogoniales. They provide the first evidence for intragenomic proliferation of group IIA introns in the Viridiplantae, revealing that intron spread in the Oedocladium lineage likely occurred by retrohoming after sequence divergence of the exon-binding sites.

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.

Organellar maturases: A window into the evolution of the spliceosome

During the evolution of eukaryotic genomes, many genes have been interrupted by intervening sequences (introns) that must be removed post-transcriptionally from RNA precursors to form mRNAs ready for translation. The origin of nuclear introns is still under debate, but one hypothesis is that the spliceosome and the intron-exon structure of genes have evolved from bacterial-type group II introns that invaded the eukaryotic genomes. The group II introns were most likely introduced into the eukaryotic genome from an α-proteobacterial predecessor of mitochondria early during the endosymbiosis event. These self-splicing and mobile introns spread through the eukaryotic genome and later degenerated. Pieces of introns became part of the general splicing machinery we know today as the spliceosome. In addition, group II introns likely brought intron maturases with them to the nucleus. Maturases are found in most bacterial introns, where they act as highly specific splicing factors for group II introns. In the spliceosome, the core protein Prp8 shows homology to group II intron-encoded maturases. While maturases are entirely intron specific, their descendant of the spliceosomal machinery, the Prp8 protein, is an extremely versatile splicing factor with multiple interacting proteins and RNAs. How could such a general player in spliceosomal splicing evolve from the monospecific bacterial maturases? Analysis of the organellar splicing machinery in plants may give clues on the evolution of nuclear splicing. Plants encode various proteins which are closely related to bacterial maturases. The organellar genomes contain one maturase each, named MatK in chloroplasts and MatR in mitochondria. In addition, several maturase genes have been found in the nucleus as well, which are acting on mitochondrial pre-RNAs. All plant maturases show sequence deviation from their progenitor bacterial maturases, and interestingly are all acting on multiple organellar group II intron targets. Moreover, they seem to function in the splicing of group II introns together with a number of additional nuclear-encoded splicing factors, possibly acting as an organellar proto-spliceosome. Together, this makes them interesting models for the early evolution of nuclear spliceosomal splicing. In this review, we summarize recent advances in our understanding of the role of plant maturases and their accessory factors in plants. This article is part of a Special Issue entitled: Chloroplast Biogenesis.

Recent mobility of plastid encoded group II introns and twintrons in five strains of the unicellular red alga Porphyridium

Group II introns are closely linked to eukaryote evolution because nuclear spliceosomal introns and the small RNAs associated with the spliceosome are thought to trace their ancient origins to these mobile elements. Therefore, elucidating how group II introns move, and how they lose mobility can potentially shed light on fundamental aspects of eukaryote biology. To this end, we studied five strains of the unicellular red alga Porphyridium purpureum that surprisingly contain 42 group II introns in their plastid genomes. We focused on a subset of these introns that encode mobility-conferring intron-encoded proteins (IEPs) and found them to be distributed among the strains in a lineage-specific manner. The reverse transcriptase and maturase domains were present in all lineages but the DNA endonuclease domain was deleted in vertically inherited introns, demonstrating a key step in the loss of mobility. P. purpureum plastid intron RNAs had a classic group IIB secondary structure despite variability in the DIII and DVI domains. We report for the first time the presence of twintrons (introns-within-introns, derived from the same mobile element) in Rhodophyta. The P. purpureum IEPs and their mobile introns provide a valuable model for the study of mobile retroelements in eukaryotes and offer promise for biotechnological applications.

The role of snRNAs in spliceosomal catalysis

The spliceosomes, large ribonucleoprotein (RNP) assemblies that remove the intervening sequences from pre-mRNAs, contain a large number of proteins and five small nuclear RNAs (snRNAs). One snRNA, U6, contains highly conserved sequences that are thought to be the functional counterparts of the RNA elements that form the active site of self-splicing group II intron ribozymes. An in vitro-assembled, protein-free complex of U6 with U2, the base-pairing partner in the spliceosomal catalytic core, can catalyze a two-step splicing reaction in the absence of all other spliceosomal factors, suggesting that the two snRNAs may form all or a large share of the spliceosomal active site. On the other hand, several spliceosomal proteins are thought to help in the formation of functionally required RNA-RNA interactions in the catalytic core. Whether they also contribute functional groups to the spliceosomal active site, and thus whether the spliceosomes are RNA or RNP enzymes remain uncertain.