Mitochondrial group I and group II introns in the sponge orders Agelasida and Axinellida

Confirmation of the first intronic sequence in the bivalvian mitochondrial genome of Macoma balthica (Linnaeus, 1758)

In 2020, the first male-type mitochondrial genome from the clam Macoma balthica was published. Apart from the unusual doubly uniparental inheritance of mtDNA, scientists observed a unique (over 4k bp long) extension in the middle of the cox2 gene. We have attempted to replicate these data by NGS DNA sequencing and explore further the expression of the long cox2 gene. In our study, we report an even longer cox2 gene (over 5.5 kbp) with no stop codon separating conserved cox2 domains, as well as, based on the rtPCR, a lower relative gene expression pattern of the middle part of the gene (5' = 1; mid = 0.46; 3' = 0.89). Lastly, we sequenced the cox2 gene transcript proving the excision of the intronic sequence.

Unprecedented frequency of mitochondrial introns in colonial bilaterians

Animal mitogenomes are typically devoid of introns. Here, we report the largest number of mitochondrial introns ever recorded from bilaterian animals. Mitochondrial introns were identified for the first time from the phylum Bryozoa. They were found in four species from three families (Order Cheilostomatida). A total of eight introns were found in the complete mitogenome of Exechonella vieirai, and five, 17 and 18 introns were found in the partial mitogenomes of Parantropora penelope, Discoporella cookae and Cupuladria biporosa, respectively. Intron-encoded protein domains reverse transcriptase and intron maturase (RVT-IM) were identified in all species. Introns in E. vieirai and P. penelope had conserved Group II intron ribozyme domains V and VI. Conserved domains were lacking from introns in D. cookae and C. biporosa, preventing their further categorization. Putative origins of metazoan introns were explored in a phylogenetic context, using an up-to-date alignment of mitochondrial RVT-IM domains. Results confirmed previous findings of multiple origins of annelid, placozoan and sponge RVT-IM domains and provided evidence for common intron donor sources across metazoan phyla. Our results corroborate growing evidence that some metazoans with regenerative abilities (i.e. placozoans, sponges, annelids and bryozoans) are susceptible to intron integration, most likely via horizontal gene transfer.

The Pet127 protein is a mitochondrial 5'-to-3' exoribonuclease from the PD-(D/E)XK superfamily involved in RNA maturation and intron degradation in yeasts

Pet127 is a mitochondrial protein found in multiple eukaryotic lineages, but absent from several taxa, including plants and animals. Distant homology suggests that it belongs to the divergent PD-(D/E)XK superfamily which includes various nucleases and related proteins. Earlier yeast genetics experiments suggest that it plays a nonessential role in RNA degradation and 5' end processing. Our phylogenetic analysis suggests that it is a primordial eukaryotic invention that was retained in diverse groups, and independently lost several times in the evolution of other organisms. We demonstrate for the first time that the fungal Pet127 protein in vitro is a processive 5'-to-3' exoribonuclease capable of digesting various substrates in a sequence nonspecific manner. Mutations in conserved residues essential in the PD-(D/E)XK superfamily active site abolish the activity of Pet127. Deletion of the PET127 gene in the pathogenic yeast Candida albicans results in a moderate increase in the steady-state levels of several transcripts and in accumulation of unspliced precursors and intronic sequences of three introns. Mutations in the active site residues result in a phenotype identical to that of the deletant, confirming that the exoribonuclease activity is related to the physiological role of the Pet127 protein. Pet127 activity is, however, not essential for maintaining the mitochondrial respiratory activity in C. albicans.

Mitogenome of a stink worm (Annelida: Travisiidae) includes degenerate group II intron that is also found in five congeneric species

Mitogenomes are useful for inferring phylogenetic relationships between organisms. Although the mitogenomes of Annelida, one of the most morphologically and ecologically diverse metazoan groups have been well sequenced, those of several families remain unexamined. This study determined the first mitogenome from the family Travisiidae (Travisia sanrikuensis), analyzed its mitogenomic features, and reconstructed a phylogeny of Sedentaria. The monophyly of the Terebellida + Arenicolida + Travisiidae clade is supported by molecular phylogenetic analysis. The placement of Travisiidae is unclear because of the lack of mitogenomes from closely related lineages. An unexpected intron appeared within the cox1 gene of T. sanrikuensis and in the same positions of five undescribed Travisia spp. Although the introns are shorter (790–1386 bp) than other group II introns, they can be considered degenerate group II introns due to type II intron maturase open reading frames, found in two of the examined species, and motifs characteristic of group II introns. This is likely the first known case in metazoans where mitochondrial group II introns obtained by a common ancestor are conserved in several descendants. Insufficient evolutionary time for intron loss in Travisiidae, or undetermined mechanisms may have helped maintain the degenerate introns.

Organellar Introns in Fungi, Algae, and Plants

Introns are ubiquitous in eukaryotic genomes and have long been considered as ‘junk RNA’ but the huge energy expenditure in their transcription, removal, and degradation indicate that they may have functional significance and can offer evolutionary advantages. In fungi, plants and algae introns make a significant contribution to the size of the organellar genomes. Organellar introns are classified as catalytic self-splicing introns that can be categorized as either Group I or Group II introns. There are some biases, with Group I introns being more frequently encountered in fungal mitochondrial genomes, whereas among plants Group II introns dominate within the mitochondrial and chloroplast genomes. Organellar introns can encode a variety of proteins, such as maturases, homing endonucleases, reverse transcriptases, and, in some cases, ribosomal proteins, along with other novel open reading frames. Although organellar introns are viewed to be ribozymes, they do interact with various intron- or nuclear genome-encoded protein factors that assist in the intron RNA to fold into competent splicing structures, or facilitate the turn-over of intron RNAs to prevent reverse splicing. Organellar introns are also known to be involved in non-canonical splicing, such as backsplicing and trans-splicing which can result in novel splicing products or, in some instances, compensate for the fragmentation of genes by recombination events. In organellar genomes, Group I and II introns may exist in nested intronic arrangements, such as introns within introns, referred to as twintrons, where splicing of the external intron may be dependent on splicing of the internal intron. These nested or complex introns, with two or three-component intron modules, are being explored as platforms for alternative splicing and their possible function as molecular switches for modulating gene expression which could be potentially applied towards heterologous gene expression. This review explores recent findings on organellar Group I and II introns, focusing on splicing and mobility mechanisms aided by associated intron/nuclear encoded proteins and their potential roles in organellar gene expression and cross talk between nuclear and organellar genomes. Potential application for these types of elements in biotechnology are also discussed.

Recent and Ongoing Horizontal Transfer of Mitochondrial Introns Between Two Fungal Tree Pathogens

Two recently introduced fungal plant pathogens (Ceratocystis lukuohia and Ceratocystis huliohia) are responsible for Rapid ‘ōhi‘a Death (ROD) in Hawai‘i. Despite being sexually incompatible, the two pathogens often co-occur in diseased ‘ōhi‘a sapwood, where genetic interaction is possible. We sequenced and annotated 33 mitochondrial genomes of the two pathogens and related species, and investigated 35 total Ceratocystis mitogenomes. Ten mtDNA regions [one group I intron, seven group II introns, and two autonomous homing endonuclease (HE) genes] were heterogeneously present in C. lukuohia mitogenomes, which were otherwise identical. Molecular surveys with specific primers showed that the 10 regions had uneven geographic distribution amongst populations of C. lukuohia. Conversely, identical orthologs of each region were present in every studied isolate of C. huliohia regardless of geographical origin. Close relatives of C. lukuohia lacked or, rarely, had few and dissimilar orthologs of the 10 regions, whereas most relatives of C. huliohia had identical or nearly identical orthologs. Each region included or worked in tandem with HE genes or reverse transcriptase/maturases that could facilitate interspecific horizontal transfers from intron-minus to intron-plus alleles. These results suggest that the 10 regions originated in C. huliohia and are actively moving to populations of C. lukuohia, perhaps through transient cytoplasmic contact of hyphal tips (anastomosis) in the wound surface of ‘ōhi‘a trees. Such contact would allow for the transfer of mitochondria followed by mitochondrial fusion or cytoplasmic exchange of intron intermediaries, which suggests that further genomic interaction may also exist between the two pathogens.

Mitogenomics reveals phylogenetic relationships of Arcoida (Mollusca, Bivalvia) and multiple independent expansions and contractions in mitochondrial genome size

Arcoida, comprising about 570 species of blood cockles, is an ecologically and economically important lineage of bivalve molluscs. Current classification of arcoids is largely based on morphology, which shows widespread homoplasy. Despite two recent studies employing multi-locus analyses with broad sampling of Arcoida, evolutionary relationships among major lineages remain controversial. Interestingly, mitochondrial genomes of several ark shell species are 2-3 times larger than those found in most bilaterians, and are among the largest bilaterian mitochondrial genomes reported to date. These results highlight the need of detailed phylogenetic study to explore evolutionary relationships within Arcoida so that the evolution of mitochondrial genome size can be understood. To this end, we sequenced 17 mitochondrial genomes and compared them with publicly available data, including those from other lineages of Arcoida with emphasis on the subclade Arcoidea species. Our phylogenetic analyses indicate that Noetiidae, Cucullaeidae and Glycymerididae are nested within a polyphyletic Arcidae. Moreover, we find multiple independent expansions and potential contractions of mitochondrial genome size, suggesting that the large mitochondrial genome is not a shared ancestral feature in Arcoida. We also examined tandem repeats and inverted repeats in non-coding regions and investigated the presence of such repeats with relation to genome size variation. Our results suggest that tandem repeats might facilitate intraspecific mitochondrial genome size variation, and that inverted repeats, which could be derived from transposons, might be responsible for mitochondrial genome expansions and contractions. We show that mitochondrial genome size in Arcoida is more dynamic than previously understood and provide insights into evolution of mitochondrial genome size variation in metazoans.

Mitochondrial genome and polymorphic microsatellite markers from the abyssal sponge Plenaster craigi Lim & Wiklund, 2017: tools for understanding the impact of deep-sea mining

The abyssal demosponge Plenaster craigi is endemic to the Clarion - Clipperton Zone (CCZ) in the NE Pacific, a region with abundant seafloor polymetallic nodules and of potential interest for mining. Plenaster craigi encrusts on these nodules and is an abundant component of the ecosystem. To assess the impact of mining operations, it is crucial to understand the genetics of this species, because its genetic diversity and connectivity across the area may be representative of other nodule-encrusting invertebrate epifauna. Here we describe and characterize 14 polymorphic microsatellite markers from this keystone species using Illumina MiSeq, tested for 75 individuals from three different areas across the CCZ, including an Area of Particular Environmental Interest (APEI-6) and two areas within the adjacent UK1 mining exploration area. The number of alleles per locus ranged from 3 to 30 (13.33 average alleles for all loci across areas). Observed and expected heterozygosity ranged from 0.909-0.048 and from 0.954-0.255, respectively. Several loci displayed significant deviation from the Hardy-Weinberg equilibrium, which appears to be common in other sponge studies. The microsatellite loci described here will be used to assess the genetic structure and connectivity on populations of the sponge across the CCZ, which will be invaluable for monitoring the impact of mining operations on its habitat. Also, we provide the annotated mitochondrial genome of P. craigi, compare its arrangement with other closely related species, and discuss the phylogenetic framework for the sponge after Maximum Likelihood and Bayesian Inference analyses using nucleotide and amino acid sequences data sets separately.

Multiple introns in a deep-sea Annelid (Decemunciger: Ampharetidae) mitochondrial genome

Wood falls provide episodic fluxes of energy to the sea floor that are degraded by a species-rich benthic fauna. Part of this rich diversity includes annelid polychaetes but unfortunately, our understanding of such fauna is limited and their genetic variability and evolutionary origins remain poorly known. In this study, we sequenced complete mitochondrial genomes from three congeneric Decemunciger (Ampharetidae) individuals that had colonized multiple wood falls in the deep (~1600 m) NE Pacific Ocean. Mitochondrial gene order within Decemunciger was similar to the three other available Terebellomorpha genomes, consistent with the relatively conserved nature of mitochondrial genomes within annelids. Unexpectedly, we found introns within the cox1, nad1 and nad4 genes of all three genomes assembled. This is the greatest number of introns observed in annelid mtDNA genomes, and possibly in bilaterians. Interestingly, the introns were of variable sizes suggesting possible evolutionary differences in the age and origins of introns. The sequence of the introns within cox1 is similar to Group II introns previously identified, suggesting that introns in the mitochondrial genome of annelids may be more widespread then realized. Phylogenetically, Decemunciger appears to be a sister clade among current vent and seep deep-sea Ampharetinae.

Mitochondrial RNA processing in absence of tRNA punctuations in octocorals

BACKGROUND

Mitogenome diversity is staggering among early branching animals with respect to size, gene density, content and order, and number of tRNA genes, especially in cnidarians. This last point is of special interest as tRNA cleavage drives the maturation of mitochondrial mRNAs and is a primary mechanism for mt-RNA processing in animals. Mitochondrial RNA processing in non-bilaterian metazoans, some of which possess a single tRNA gene in their mitogenomes, is essentially unstudied despite its importance in understanding the evolution of mitochondrial transcription in animals.

RESULTS

We characterized the mature mitochondrial mRNA transcripts in a species of the octocoral genus Sinularia (Alcyoniidae: Octocorallia), and defined precise boundaries of transcription units using different molecular methods. Most mt-mRNAs were polycistronic units containing two or three genes and 5' and/or 3' untranslated regions of varied length. The octocoral specific, mtDNA-encoded mismatch repair gene, the mtMutS, was found to undergo alternative polyadenylation, and exhibited differential expression of alternate transcripts suggesting a unique regulatory mechanism for this gene. In addition, a long noncoding RNA complementary to the ATP6 gene (lncATP6) potentially involved in antisense regulation was detected.

CONCLUSIONS

Mt-mRNA processing in octocorals possessing a single mt-tRNA is complex. Considering the variety of mitogenome arrangements known in cnidarians, and in general among non-bilaterian metazoans, our findings provide a first glimpse into the complex mtDNA transcription, mt-mRNA processing, and regulation among early branching animals and represent a first step towards understanding its functional and evolutionary implications.

Evolution of group I introns in Porifera: new evidence for intron mobility and implications for DNA barcoding

BACKGROUND

Mitochondrial introns intermit coding regions of genes and feature characteristic secondary structures and splicing mechanisms. In metazoans, mitochondrial introns have only been detected in sponges, cnidarians, placozoans and one annelid species. Within demosponges, group I and group II introns are present in six families. Based on different insertion sites within the cox1 gene and secondary structures, four types of group I and two types of group II introns are known, which can harbor up to three encoding homing endonuclease genes (HEG) of the LAGLIDADG family (group I) and/or reverse transcriptase (group II). However, only little is known about sponge intron mobility, transmission, and origin due to the lack of a comprehensive dataset. We analyzed the largest dataset on sponge mitochondrial group I introns to date: 95 specimens, from 11 different sponge genera which provided novel insights into the evolution of group I introns.

RESULTS

For the first time group I introns were detected in four genera of the sponge family Scleritodermidae (Scleritoderma, Microscleroderma, Aciculites, Setidium). We demonstrated that group I introns in sponges aggregate in the most conserved regions of cox1. We showed that co-occurrence of two introns in cox1 is unique among metazoans, but not uncommon in sponges. However, this combination always associates an active intron with a degenerating one. Earlier hypotheses of HGT were confirmed and for the first time VGT and secondary losses of introns conclusively demonstrated.

CONCLUSION

This study validates the subclass Spirophorina (Tetractinellida) as an intron hotspot in sponges. Our analyses confirm that most sponge group I introns probably originated from fungi. DNA barcoding is discussed and the application of alternative primers suggested.

Evolutionary and biogeographical implications of degraded LAGLIDADG endonuclease functionality and group I intron occurrence in stony corals (Scleractinia) and mushroom corals (Corallimorpharia)

Group I introns and homing endonuclease genes (HEGs) are mobile genetic elements, capable of invading target sequences in intron-less genomes. LAGLIDADG HEGs are the largest family of endonucleases, playing a key role in the mobility of group I introns in a process known as ‘homing’. Group I introns and HEGs are rare in metazoans, and can be mainly found inserted in the COXI gene of some sponges and cnidarians, including stony corals (Scleractinia) and mushroom corals (Corallimorpharia). Vertical and horizontal intron transfer mechanisms have been proposed as explanations for intron occurrence in cnidarians. However, the central role of LAGLIDADG motifs in intron mobility mechanisms remains poorly understood. To resolve questions regarding the evolutionary origin and distribution of group I introns and HEGs in Scleractinia and Corallimorpharia, we examined intron/HEGs sequences within a comprehensive phylogenetic framework. Analyses of LAGLIDADG motif conservation showed a high degree of degradation in complex Scleractinia and Corallimorpharia. Moreover, the two motifs lack the respective acidic residues necessary for metal-ion binding and catalysis, potentially impairing horizontal intron mobility. In contrast, both motifs are highly conserved within robust Scleractinia, indicating a fully functional endonuclease capable of promoting horizontal intron transference. A higher rate of non-synonymous substitutions (K_a) detected in the HEGs of complex Scleractinia and Corallimorpharia suggests degradation of the HEG, whereas lower K_a rates in robust Scleractinia are consistent with a scenario of purifying selection. Molecular-clock analyses and ancestral inference of intron type indicated an earlier intron insertion in complex Scleractinia and Corallimorpharia in comparison to robust Scleractinia. These findings suggest that the lack of horizontal intron transfers in the former two groups is related to an age-dependent degradation of the endonuclease activity. Moreover, they also explain the peculiar geographical patterns of introns in stony and mushroom corals.

Animal Mitochondrial DNA as We Do Not Know It: mt-Genome Organization and Evolution in Nonbilaterian Lineages

Animal mitochondrial DNA (mtDNA) is commonly described as a small, circular molecule that is conserved in size, gene content, and organization. Data collected in the last decade have challenged this view by revealing considerable diversity in animal mitochondrial genome organization. Much of this diversity has been found in nonbilaterian animals (phyla Cnidaria, Ctenophora, Placozoa, and Porifera), which, from a phylogenetic perspective, form the main branches of the animal tree along with Bilateria. Within these groups, mt-genomes are characterized by varying numbers of both linear and circular chromosomes, extra genes (e.g. atp9, polB, tatC), large variation in the number of encoded mitochondrial transfer RNAs (tRNAs) (0-25), at least seven different genetic codes, presence/absence of introns, tRNA and mRNA editing, fragmented ribosomal RNA genes, translational frameshifting, highly variable substitution rates, and a large range of genome sizes. This newly discovered diversity allows a better understanding of the evolutionary plasticity and conservation of animal mtDNA and provides insights into the molecular and evolutionary mechanisms shaping mitochondrial genomes.