The complete mitochondrial genome of the yeast Kluyveromyces thermotolerans

Discordant evolution of mitochondrial and nuclear yeast genomes at population level

Background

Mitochondria are essential organelles partially regulated by their own genomes. The mitochondrial genome maintenance and inheritance differ from the nuclear genome, potentially uncoupling their evolutionary trajectories. Here, we analysed mitochondrial sequences obtained from the 1011 Saccharomyces cerevisiae strain collection and identified pronounced differences with their nuclear genome counterparts.

Results

In contrast with pre-whole genome duplication fungal species, S. cerevisiae mitochondrial genomes show higher genetic diversity compared to the nuclear genomes. Strikingly, mitochondrial genomes appear to be highly admixed, resulting in a complex interconnected phylogeny with a weak grouping of isolates, whereas interspecies introgressions are very rare. Complete genome assemblies revealed that structural rearrangements are nearly absent with rare inversions detected. We tracked intron variation in COX1 and COB to infer gain and loss events throughout the species evolutionary history. Mitochondrial genome copy number is connected with the nuclear genome and linearly scale up with ploidy. We observed rare cases of naturally occurring mitochondrial DNA loss, petite, with a subset of them that do not suffer the expected growth defect in fermentable rich media.

Conclusions

Overall, our results illustrate how differences in the biology of two genomes coexisting in the same cells can lead to discordant evolutionary histories.

Molecular Characterization and Enological Potential of A High Lactic Acid-Producing Lachancea thermotolerans Vineyard Strain

Lactic acid production is an important feature of the yeast Lachancea thermotolerans that has gained increasing interest in winemaking. In particular, in light of climate change, the biological acidification and ethanol reduction by the use of selected yeast strains may counteract the effect of global warming in wines. Here, the enological potential of a high lactate-producing L. thermotolerans strain (P-HO1) in mixed fermentations with S. cerevisiae was examined. Among the different inoculation schemes evaluated, the most successful implantation of L. thermotolerans was accomplished by sequential inoculation of S. cerevisiae, i.e., at 1% vol. ethanol. P-HO1produced the highest levels of lactic acid ever recorded in mixed fermentations (10.4 g/L), increasing thereby the acidity and reducing ethanol by 1.6% vol. L. thermotolerans was also associated with increases in ethyl isobutyrate (strawberry aroma), free SO₂, organoleptically perceived citric nuances and aftertaste. To start uncovering the molecular mechanisms of lactate biosynthesis in L. thermotolerans, the relative expressions of the three lactate dehydrogenase (LDH) paralogous genes, which encode the key enzyme for lactate biosynthesis, along with the alcohol dehydrogenase paralogs (ADHs) were determined. Present results point to the possible implication of LDH2, but not of other LDH or ADH genes, in the high production of lactic acid in certain strains at the expense of ethanol. Taken together, the important enological features of P-HO1 highlighted here, and potentially of other L. thermotolerans strains, indicate its great importance in modern winemaking, particularly in the light of the upcoming climate change and its consequences in the grape/wine system.

The evolutionary history of Saccharomyces species inferred from completed mitochondrial genomes and revision in the 'yeast mitochondrial genetic code'

The yeast Saccharomyces are widely used to test ecological and evolutionary hypotheses. A large number of nuclear genomic DNA sequences are available, but mitochondrial genomic data are insufficient. We completed mitochondrial DNA (mtDNA) sequencing from Illumina MiSeq reads for all Saccharomyces species. All are circularly mapped molecules decreasing in size with phylogenetic distance from Saccharomyces cerevisiae but with similar gene content including regulatory and selfish elements like origins of replication, introns, free-standing open reading frames or GC clusters. Their most profound feature is species-specific alteration in gene order. The genetic code slightly differs from well-established yeast mitochondrial code as GUG is used rarely as the translation start and CGA and CGC code for arginine. The multilocus phylogeny, inferred from mtDNA, does not correlate with the trees derived from nuclear genes. mtDNA data demonstrate that Saccharomyces cariocanus should be assigned as a separate species and Saccharomyces bayanus CBS 380T should not be considered as a distinct species due to mtDNA nearly identical to Saccharomyces uvarum mtDNA. Apparently, comparison of mtDNAs should not be neglected in genomic studies as it is an important tool to understand the origin and evolutionary history of some yeast species.

The mitochondrial genome of the plant-pathogenic fungus Stemphylium lycopersici uncovers a dynamic structure due to repetitive and mobile elements

Stemphylium lycopersici (Pleosporales) is a plant-pathogenic fungus that has been associated with a broad range of plant-hosts worldwide. It is one of the causative agents of gray leaf spot disease in tomato and pepper. The aim of this work was to characterize the mitochondrial genome of S. lycopersici CIDEFI-216, to use it to trace taxonomic relationships with other fungal taxa and to get insights into the evolutionary history of this phytopathogen. The complete mitochondrial genome was assembled into a circular double-stranded DNA molecule of 75,911 bp that harbors a set of 37 protein-coding genes, 2 rRNA genes (rns and rnl) and 28 tRNA genes, which are transcribed from both sense and antisense strands. Remarkably, its gene repertoire lacks both atp8 and atp9, contains a free-standing gene for the ribosomal protein S3 (rps3) and includes 13 genes with homing endonuclease domains that are mostly located within its 15 group I introns. Strikingly, subunits 1 and 2 of cytochrome oxidase are encoded by a single continuous open reading frame (ORF). A comparative mitogenomic analysis revealed the large extent of structural rearrangements among representatives of Pleosporales, showing the plasticity of their mitochondrial genomes. Finally, an exhaustive phylogenetic analysis of the subphylum Pezizomycotina based on mitochondrial data reconstructed their relationships in concordance with several studies based on nuclear data. This is the first report of a mitochondrial genome belonging to a representative of the family Pleosporaceae.

Development of microsatellite markers for Lachancea thermotolerans typing and population structure of wine-associated isolates

Lachancea (Kluyveromyces) thermotolerans is an important member of the grape/wine yeast community with great technological potential for the wine industry. Although several molecular marker techniques have been developed for typing different yeast species, no one has been designed so far for L. thermotolerans. Here we present a simple and efficient method based on a multilocus SSR analysis for molecular typing and genetic diversity assessment of L. thermotolerans isolates. Following whole genome screening, five polymorphic microsatellite markers were selected and tested on a panel of grape isolates from different vineyards of two geographically separated viticultural zones, Nemea and Peza, in Greece. The SSR method proved quite discriminatory as compared to tandem repeat-tRNA-PCR, a fingerprinting method for typing non-Saccharomyces yeasts. Genetic analysis based on SSR data revealed a clear structure between the populations of the two zones. Furthermore, significant differences were also detected in a number of phenotypic characters of enological interest. A positive correlation was observed between phenotypic and genotypic diversity. Taking together, present results support the microbial terroir concept in the case of L. thermotolerans in Greece, which is an important prerequisite for the exploitation of selected genotypes as fermentation starters with region-specific characters.

The mitochondrial genome of the ethanol-metabolizing, wine cellar mold Zasmidium cellare is the smallest for a filamentous ascomycete

Fungi in the class Dothideomycetes often live in extreme environments or have unusual physiology. One of these, the wine cellar mold Zasmidium cellare, produces thick curtains of mycelia in cellars with high humidity, and its ability to metabolize volatile organic compounds is thought to improve air quality. Whether these abilities have affected its mitochondrial genome is not known. To fill this gap, the circular-mapping mitochondrial genome of Z. cellare was sequenced and, at only 23 743 bp, is the smallest reported for a filamentous fungus. Genes were encoded on both strands with a single change of direction, different from most other fungi but consistent with the Dothideomycetes. Other than its small size, the only unusual feature of the Z. cellare mitochondrial genome was two copies of a 110-bp sequence that were duplicated, inverted and separated by approximately 1 kb. This inverted-repeat sequence confused the assembly program but appears to have no functional significance. The small size of the Z. cellare mitochondrial genome was due to slightly smaller genes, lack of introns and non-essential genes, reduced intergenic spacers and very few ORFs relative to other fungi rather than a loss of essential genes. Whether this reduction facilitates its unusual biology remains unknown.

The Genetics of Non-conventional Wine Yeasts: Current Knowledge and Future Challenges

Saccharomyces cerevisiae is by far the most widely used yeast in oenology. However, during the last decade, several other yeasts species has been purposed for winemaking as they could positively impact wine quality. Some of these non-conventional yeasts (Torulaspora delbrueckii, Metschnikowia pulcherrima, Pichia kluyveri, Lachancea thermotolerans, etc.) are now proposed as starters culture for winemakers in mixed fermentation with S. cerevisiae, and several others are the subject of various studies (Hanseniaspora uvarum, Starmerella bacillaris, etc.). Along with their biotechnological use, the knowledge of these non-conventional yeasts greatly increased these last 10 years. The aim of this review is to describe the last updates and the current state-of-art of the genetics of non-conventional yeasts (including S. uvarum, T. delbrueckii, S. bacillaris, etc.). We describe how genomics and genetics tools provide new data into the population structure and biodiversity of non-conventional yeasts in winemaking environments. Future challenges will lie on the development of selection programs and/or genetic improvement of these non-conventional species. We discuss how genetics, genomics and the advances in next-generation sequencing will help the wine industry to develop the biotechnological use of non-conventional yeasts to improve the quality and differentiation of wines.

Mitochondrial genome evolution in yeasts: an all-encompassing view

Mitochondria are important organelles that harbor their own genomes encoding a key set of proteins that ensure respiration and provide the eukaryotic cell with energy. Recent advances in high-throughput sequencing technologies present a unique opportunity to explore mitochondrial (mt) genome evolution. The Saccharomycotina yeasts have proven to be the leading organisms for mt comparative and population genomics. In fact, the explosion of complete yeast mt genome sequences has allowed for a broader view of the mt diversity across this incredibly diverse subphylum, both within and between closely related species. Here, we summarize the present state of yeast mitogenomics, including the currently available data and what it reveals concerning the diversity of content, organization, structure and evolution of mt genomes.

Population genomic analysis reveals highly conserved mitochondrial genomes in the yeast species Lachancea thermotolerans

The increasing availability of mitochondrial (mt) sequence data from various yeasts provides a tool to study genomic evolution within and between different species. While the genomes from a range of lineages are available, there is a lack of information concerning intraspecific mtDNA diversity. Here, we analyzed the mt genomes of 50 strains from Lachancea thermotolerans, a protoploid yeast species that has been isolated from several locations (Europe, Asia, Australia, South Africa, and North / South America) and ecological sources (fruit, tree exudate, plant material, and grape and agave fermentations). Protein-coding genes from the mtDNA were used to construct a phylogeny, which reflected a similar, yet less resolved topology than the phylogenetic tree of 50 nuclear genes. In comparison to its sister species Lachancea kluyveri, L. thermotolerans has a smaller mt genome. This is due to shorter intergenic regions and fewer introns, of which the latter are only found in COX1. We revealed that L. kluyveri and L. thermotolerans share similar levels of intraspecific divergence concerning the nuclear genomes. However, L. thermotolerans has a more highly conserved mt genome with the coding regions characterized by low rates of nonsynonymous substitution. Thus, in the mt genomes of L. thermotolerans, stronger purifying selection and lower mutation rates potentially shape genome diversity in contract to what was found for L. kluyveri, demonstrating that the factors driving mt genome evolution are different even between closely related species.

Comparative mitochondrial genomics within and among yeast species of the Lachancea genus

Yeasts are leading model organisms for mitochondrial genome studies. The explosion of complete sequence of yeast mitochondrial (mt) genomes revealed a wide diversity of organization and structure between species. Recently, genome-wide polymorphism survey on the mt genome of isolates of a single species, Lachancea kluyveri, was also performed. To compare the mitochondrial genome evolution at two hierarchical levels: within and among closely related species, we focused on five species of the Lachancea genus, which are close relatives of L. kluyveri. Hence, we sequenced the complete mt genome of L. dasiensis, L. nothofagi, L. mirantina, L. fantastica and L. meyersii. The phylogeny of the Lachancea genus was explored using these data. Analysis of intra- and interspecific variability across the whole Lachancea genus led to the same conclusions regarding the mitochondrial genome evolution. These genomes exhibit a similar architecture and are completely syntenic. Nevertheless, genome sizes vary considerably because of the variations of the intergenic regions and the intron content, contributing to mitochondrial genome plasticity. The high variability of the intergenic regions stands in contrast to the high level of similarity of protein sequences. Quantification of the selective constraints clearly revealed that most of the mitochondrial genes are under purifying selection in the whole genus.

Evolutionary patterns of non-coding RNAs

A plethora of new functions of non-coding RNAs (ncRNAs) have been discovered in past few years. In fact, RNA is emerging as the central player in cellular regulation, taking on active roles in multiple regulatory layers from transcription, RNA maturation, and RNA modification to translational regulation. Nevertheless, very little is known about the evolution of this "Modern RNA World" and its components. In this contribution, we attempt to provide at least a cursory overview of the diversity of ncRNAs and functional RNA motifs in non-translated regions of regular messenger RNAs (mRNAs) with an emphasis on evolutionary questions. This survey is complemented by an in-depth analysis of examples from different classes of RNAs focusing mostly on their evolution in the vertebrate lineage. We present a survey of Y RNA genes in vertebrates and study the molecular evolution of the U7 snRNA, the snoRNAs E1/U17, E2, and E3, the Y RNA family, the let-7 microRNA (miRNA) family, and the mRNA-like evf-1 gene. We furthermore discuss the statistical distribution of miRNAs in metazoans, which suggests an explosive increase in the miRNA repertoire in vertebrates. The analysis of the transcription of ncRNAs suggests that small RNAs in general are genetically mobile in the sense that their association with a hostgene (e.g. when transcribed from introns of a mRNA) can change on evolutionary time scales. The let-7 family demonstrates, that even the mode of transcription (as intron or as exon) can change among paralogous ncRNA.

Mitochondrial genome evolution in a single protoploid yeast species

Mitochondria are organelles, which play a key role in some essential functions, including respiration, metabolite biosynthesis, ion homeostasis, and apoptosis. The vast numbers of mitochondrial DNA (mtDNA) sequences of various yeast species, which have recently been published, have also helped to elucidate the structural diversity of these genomes. Although a large corpus of data are now available on the diversity of yeast species, little is known so far about the mtDNA diversity in single yeast species. To study the genetic variations occurring in the mtDNA of wild yeast isolates, we performed a genome-wide polymorphism survey on the mtDNA of 18 Lachancea kluyveri (formerly Saccharomyces kluyveri) strains. We determined the complete mt genome sequences of strains isolated from various geographical locations (in North America, Asia, and Europe) and ecological niches (Drosophila, tree exudates, soil). The mt genome of the NCYC 543 reference strain is 51,525 bp long. It contains the same core of genes as Lachancea thermotolerans, the nearest relative to L. kluyveri. To explore the mt genome variations in a single yeast species, we compared the mtDNAs of the 18 isolates. The phylogeny and population structure of L. kluyveri provide clear-cut evidence for the existence of well-defined geographically isolated lineages. Although these genomes are completely syntenic, their size and the intron content were found to vary among the isolates studied. These genomes are highly polymorphic, showing an average diversity of 28.5 SNPs/kb and 6.6 indels/kb. Analysis of the SNP and indel patterns showed the existence of a particularly high overall level of polymorphism in the intergenic regions. The dN/dS ratios obtained are consistent with purifying selection in all these genes, with the noteworthy exception of the VAR1 gene, which gave a very high ratio. These data suggest that the intergenic regions have evolved very fast in yeast mitochondrial genomes.

Mitochondrial inheritance in budding yeasts: towards an integrated understanding

Recent advances in yeast mitogenomics have significantly contributed to our understanding of the diversity of organization, structure and topology in the mitochondrial genome of budding yeasts. In parallel, new insights on mitochondrial DNA (mtDNA) inheritance in the model organism Saccharomyces cerevisiae highlighted an integrated scenario where recombination, replication and segregation of mtDNA are intricately linked to mitochondrial nucleoid (mt-nucleoid) structure and organelle sorting. In addition to this, recent discoveries of bifunctional roles of some mitochondrial proteins have interesting implications on mito-nuclear genome interactions and the relationship between mtDNA inheritance, yeast fitness and speciation. This review summarizes the current knowledge on yeast mitogenomics, mtDNA inheritance with regard to mt-nucleoid structure and organelle dynamics, and mito-nuclear genome interactions.

The subunit composition of mitochondrial NADH:ubiquinone oxidoreductase (complex I) from Pichia pastoris

Respiratory complex I (NADH:quinone oxidoreductase) is an entry point to the electron transport chain in the mitochondria of many eukaryotes. It is a large, multisubunit enzyme with a hydrophilic domain in the matrix and a hydrophobic domain in the mitochondrial inner membrane. Here we present a comprehensive analysis of the protein composition and post-translational modifications of complex I from Pichia pastoris, using a combination of proteomic and bioinformatic approaches. Forty-one subunits were identified in P. pastoris complex I, comprising the 14 core (conserved) subunits and 27 supernumerary subunits; seven of the core subunits are mitochondrial encoded. Three of the supernumerary subunits (named NUSM, NUTM, and NUUM) have not been observed previously in any species of complex I. However, homologues to all three of them are present in either Yarrowia lipolytica or Pichia angusta complex I. P. pastoris complex I has 39 subunits in common with Y. lipolytica complex I, 37 in common with N. crassa complex I, and 35 in common with the bovine enzyme. The mitochondrial encoded subunits (translated by the mold mitochondrial genetic code) retain their N-α-formyl methionine residues. At least eight subunits are N-α-acetylated, but the N-terminal modifications of the nuclear encoded subunits are not well-conserved. A combination of two methods of protein separation (SDS-PAGE and HPLC) and three different mass spectrometry techniques (peptide mass fingerprinting, tandem MS and molecular mass measurements) were required to define the protein complement of P. pastoris complex I. This requirement highlights the need for inclusive and comprehensive strategies for the characterization of challenging membrane-bound protein complexes containing both hydrophilic and hydrophobic components.

Mitochondrial genome from the facultative anaerobe and petite-positive yeast Dekkera bruxellensis contains the NADH dehydrogenase subunit genes

The progenitor of the Dekkera/Brettanomyces clade separated from the Saccharomyces/Kluyveromyces clade over 200 million years ago. However, within both clades, several lineages developed similar physiological traits. Both Saccharomyces cerevisiae and Dekkera bruxellensis are facultative anaerobes; in the presence of excess oxygen and sugars, they accumulate ethanol (Crabtree effect) and they both spontaneously generate respiratory-deficient mutants (petites). In order to understand the role of respiratory metabolism, the mitochondrial DNA (mtDNA) molecules of two Dekkera/Brettanomyces species were analysed. Dekkera bruxellensis mtDNA shares several properties with S. cerevisiae, such as the large genome size (76 453 bp), and the organization of the intergenic sequences consisting of spacious AT-rich regions containing a number of hairpin GC-rich cluster-like elements. In addition to a basic set of the mitochondrial genes coding for the components of cytochrome oxidase, cytochrome b, subunits of ATPase, two rRNA subunits and 25 tRNAs, D. bruxellensis also carries genes for the NADH dehydrogenase complex. Apparently, in yeast, the loss of this complex is not a precondition to develop a petite-positive, Crabtree-positive and anaerobic nature. On the other hand, mtDNA from a petite-negative Brettanomyces custersianus is much smaller (30 058 bp); it contains a similar gene set and has only short intergenic sequences.

Mitochondrial DNA heteroplasmy in Candida glabrata after mitochondrial transformation

Genetic manipulation of mitochondrial DNA (mtDNA) is the most direct method for investigating mtDNA, but until now, this has been achieved only in the diploid yeast Saccharomyces cerevisiae. In this study, the ATP6 gene on mtDNA of the haploid yeast Candida glabrata (Torulopsis glabrata) was deleted by biolistic transformation of DNA fragments with a recoded ARG8(m) mitochondrial genetic marker, flanked by homologous arms to the ATP6 gene. Transformants were identified by arginine prototrophy. However, in the transformants, the original mtDNA was not lost spontaneously, even under arginine selective pressure. Moreover, the mtDNA transformants selectively lost the transformed mtDNA under aerobic conditions. The mtDNA heteroplasmy in the transformants was characterized by PCR, quantitative PCR, and Southern blotting, showing that the heteroplasmy was relatively stable in the absence of arginine. Aerobic conditions facilitated the loss of the original mtDNA, and anaerobic conditions favored loss of the transformed mtDNA. Moreover, detailed investigations showed that increases in reactive oxygen species in mitochondria lacking ATP6, along with their equal cell division, played important roles in determining the dynamics of heteroplasmy. Based on our analysis of mtDNA heteroplasmy in C. glabrata, we were able to generate homoplasmic Deltaatp6 mtDNA strains.

Comparative genomics of protoploid Saccharomycetaceae

Our knowledge of yeast genomes remains largely dominated by the extensive studies on Saccharomyces cerevisiae and the consequences of its ancestral duplication, leaving the evolution of the entire class of hemiascomycetes only partly explored. We concentrate here on five species of Saccharomycetaceae, a large subdivision of hemiascomycetes, that we call "protoploid" because they diverged from the S. cerevisiae lineage prior to its genome duplication. We determined the complete genome sequences of three of these species: Kluyveromyces (Lachancea) thermotolerans and Saccharomyces (Lachancea) kluyveri (two members of the newly described Lachancea clade), and Zygosaccharomyces rouxii. We included in our comparisons the previously available sequences of Kluyveromyces lactis and Ashbya (Eremothecium) gossypii. Despite their broad evolutionary range and significant individual variations in each lineage, the five protoploid Saccharomycetaceae share a core repertoire of approximately 3300 protein families and a high degree of conserved synteny. Synteny blocks were used to define gene orthology and to infer ancestors. Far from representing minimal genomes without redundancy, the five protoploid yeasts contain numerous copies of paralogous genes, either dispersed or in tandem arrays, that, altogether, constitute a third of each genome. Ancient, conserved paralogs as well as novel, lineage-specific paralogs were identified.

GenomeVx: simple web-based creation of editable circular chromosome maps

We describe GenomeVx, a web-based tool for making editable, publication-quality, maps of mitochondrial and chloroplast genomes and of large plasmids. These maps show the location of genes and chromosomal features as well as a position scale. The program takes as input either raw feature positions or GenBank records. In the latter case, features are automatically extracted and colored, an example of which is given. Output is in the Adobe Portable Document Format (PDF) and can be edited by programs such as Adobe Illustrator.

AVAILABILITY

GenomeVx is available at http://wolfe.gen.tcd.ie/GenomeVx

The mechanisms of codon reassignments in mitochondrial genetic codes

Many cases of nonstandard genetic codes are known in mitochondrial genomes. We carry out analysis of phylogeny and codon usage of organisms for which the complete mitochondrial genome is available, and we determine the most likely mechanism for codon reassignment in each case. Reassignment events can be classified according to the gain-loss framework. The “gain” represents the appearance of a new tRNA for the reassigned codon or the change of an existing tRNA such that it gains the ability to pair with the codon. The “loss” represents the deletion of a tRNA or the change in a tRNA so that it no longer translates the codon. One possible mechanism is codon disappearance (CD), where the codon disappears from the genome prior to the gain and loss events. In the alternative mechanisms the codon does not disappear. In the unassigned codon mechanism, the loss occurs first, whereas in the ambiguous intermediate mechanism, the gain occurs first. Codon usage analysis gives clear evidence of cases where the codon disappeared at the point of the reassignment and also cases where it did not disappear. CD is the probable explanation for stop to sense reassignments and a small number of reassignments of sense codons. However, the majority of sense-to-sense reassignments cannot be explained by CD. In the latter cases, by analysis of the presence or absence of tRNAs in the genome and of the changes in tRNA sequences, it is sometimes possible to distinguish between the unassigned codon and the ambiguous intermediate mechanisms. We emphasize that not all reassignments follow the same scenario and that it is necessary to consider the details of each case carefully.

Hybrid E. coli--Mitochondrial ribonuclease P RNAs are catalytically active

RNase P is a ribonucleoprotein that cleaves tRNA precursors at their 5'-end. Mitochondrion-encoded RNA subunits of mitochondrial RNase P (mtP-RNA) have been identified in jakobid flagellates such as Reclinomonas americana, in the prasinophyte alga Nephroselmis olivacea, and in several ascomycete and zygomycete fungi. While the structures of ascomycete mtP-RNAs are highly reduced, those of jakobids, prasinophytes, and zygomycetes retain most conserved features of their bacterial counterparts. Therefore, these mtP-RNAs might be active in vitro in the absence of a protein subunit, as are bacterial P-RNAs. Here we present a comparative structural analysis including seven newly characterized jakobid mtP-RNAs. We investigate ribozyme activities of mtP-RNAs and find that even the most bacteria-like molecules of jakobids are inactive in vitro. However, when certain domains of jakobid and N. olivacea mtP-RNAs are replaced with those from Escherichia coli, these hybrid RNAs show catalytic activity. In vitro mutagenesis of these hybrid mtP-RNAs shows that various structural elements play a critical role in ribozyme catalysis and provide further support for the presence of these elements in mtP-RNAs. These include GNRA tetraloops in helix P14 and P18 of Jakoba libera, and a remnant P3 pairing in Seculamonas ecuadoriensis. Finally, we will discuss reasons for the failure of mtP-RNAs to show catalytic activity in the absence of P-proteins based on our mutagenesis analysis.