Two for the price of one: a bifunctional intron-encoded DNA endonuclease-RNA maturase

Interconnected roles of fungal nuclear- and intron-encoded maturases: at the crossroads of mitochondrial intron splicing

Group I and II introns are large catalytic RNAs (ribozymes) that are frequently encountered in fungal mitochondrial genomes. The discovery of respiratory mutants linked to intron splicing defects demonstrated that for the efficient removal of organellar introns there appears to be a requirement of protein splicing factors. These splicing factors can be intron-encoded proteins with maturase activities that usually promote the splicing of the introns that encode them (cis-acting) and/or nuclear-encoded factors that can promote the splicing of a range of different introns (trans-acting). Compared to plants organellar introns, fungal mitochondrial intron splicing is still poorly explored, especially in terms of the synergy of nuclear factors with intron-encoded maturases that has direct impact on splicing through their association with intron RNA. In addition, nuclear-encoded accessory factors might drive the splicing impetus through translational activation, mitoribosome assembly, and phosphorylation-mediated RNA turnover. This review explores protein-assisted splicing of introns by nuclear and mitochondrial-encoded maturases as a means of mitonuclear interplay that could respond to environmental and developmental factors promoting phenotypic adaptation and potentially speciation. It also highlights key evolutionary events that have led to changes in structure and ATP-dependence to accommodate the dual functionality of nuclear and organellar splicing factors.

Comprehensive analysis of insertion sequences within rRNA genes of CPR bacteria and biochemical characterization of a homing endonuclease encoded by these sequences

The Candidate Phyla Radiation (CPR) represents an extensive bacterial clade comprising primarily uncultured lineages and is distinguished from other bacteria by a significant prevalence of insertion sequences (ISs) within their rRNA genes. However, our understanding of the taxonomic distribution and characteristics of these ISs remains limited. In this study, we used a comprehensive approach to systematically determine the nature of the rRNA ISs in CPR bacteria. The analysis of hundreds of rRNA gene sequences across 65 CPR phyla revealed that ISs are present in 48% of 16S rRNA genes and 82% of 23S rRNA genes, indicating a broad distribution across the CPR clade, with exceptions in the 16S and 23S rRNA genes of Candidatus (Ca.) Saccharibacteria and the 16S rRNA genes of Ca. Peregrinibacteria. Over half the ISs display a group-I-intron-like structure, whereas specific 16S rRNA gene ISs display features reminiscent of group II introns. The ISs frequently encode proteins with homing endonuclease (HE) domains, centered around the LAGLIDADG motif. The LAGLIDADG HE (LHE) proteins encoded by the rRNA ISs of CPR bacteria predominantly have a single-domain structure, deviating from the usual single- or double-domain configuration observed in typical prokaryotic LHEs. Experimental analysis of one LHE protein, I-ShaI from Ca. Shapirobacteria, confirmed that its endonuclease activity targets the DNA sequence of its insertion site, and chemical cross-linking experiments demonstrated its capacity to form homodimers. These results provide robust evidence supporting the hypothesis that the explosive proliferation of rRNA ISs in CPR bacteria was facilitated by mechanisms involving LHEs.

IMPORTANCE

Insertion sequences (ISs) in rRNA genes are relatively limited and infrequent in most bacterial phyla. With a comprehensive bioinformatic analysis, we show that in CPR bacteria, these ISs occur in 48% of 16S rRNA genes and 82% of 23S rRNA genes. We also report the systematic and biochemical characterization of the LAGLIDADG homing endonucleases (LHEs) encoded by these ISs in the first such analysis of the CPR bacteria. This study significantly extends our understanding of the phylogenetic positions of rRNA ISs within CPR bacteria and the biochemical features of their LHEs.

Efficacy of maturase K and rpL20 protein extracted from C. procera leaves on Anophelesstephensi

The present work is carried out to protein isolation, purification, and characterization from leaves, stem, and seed of C. procera and to evaluate the larvicidal potential on Anopheles stephensi. The whole protein was isolated using protein extraction buffer and precipitated by ammonium sulphate and larvicidal active protein was purified by the column chromatography. The homogeneity of larvicidal protein was confirmed by the SDS-PAGE. The identification of protein was done by the HPLC and LC-MS/ESI-MS. The crude protein from leaves showed 100% mortality of 3rd instar larvae of An. stephensi at the concentration of 5.5 mg/ml after 24 h of exposure. The crude protein from stem showed 25% mortality and no mortality observed was observed in seed protein. The leaves crude protein was further purified by ion exchange chromatography and eluted fractions were tested for larvicidal potential. The purified single protein fractions L2 and L3 from C. procera leaves showed 100% mortality at concentration of 0.06 mg/ml. The homogeneity of purified protein was confirmed by SDS-PAGE and two bands of 26 kDa and 15 kDa protein were observed. The peptide sequence "R.SQMLENSFLIENVMKR.L" was identified in the trypsin digested homogenous protein fraction L2 and "R.DRGSQKR.N" peptide sequence in L3 fraction by LC-MS/ESI-MS. The CprL2 peptide showed the sequence similarity with the protein maturase K and CprL3 peptide showed the sequence similarity with ribosomal protein L20 of C. procera. The conserved functional domain was also identified in both the CprL2 and CprL3 peptide. The identified proteins showed strong larvicidal efficacy at very low concentration. The identified proteins are novel and natural larvicidal agents against An. stephensi and hence can be used to control the malaria.

Group I introns: Structure, splicing and their applications in medical mycology

Group I introns are small RNAs (250-500 nt) capable of catalyzing their own splicing from the precursor RNA. They are widely distributed across the tree of life and have intricate relationships with their host genomes. In this work, we review its basic structure, self-splicing and its mechanisms of gene mobility. As they are widely found in unicellular eukaryotes, especially fungi, we gathered information regarding their possible impact on the physiology of fungal cells and the possible application of these introns in medical mycology.

RNA processing and degradation mechanisms shaping the mitochondrial transcriptome of budding yeasts

Yeast mitochondrial genes are expressed as polycistronic transcription units that contain RNAs from different classes and show great evolutionary variability. The promoters are simple, and transcriptional control is rudimentary. Posttranscriptional mechanisms involving RNA maturation, stability, and degradation are thus the main force shaping the transcriptome and determining the expression levels of individual genes. Primary transcripts are fragmented by tRNA excision by RNase P and tRNase Z, additional processing events occur at the dodecamer site at the 3' end of protein-coding sequences. groups I and II introns are excised in a self-splicing reaction that is supported by protein splicing factors encoded by the nuclear genes, or by the introns themselves. The 3'-to-5' exoribonucleolytic complex called mtEXO is the main RNA degradation activity involved in RNA turnover and processing, supported by an auxiliary 5'-to-3' exoribonuclease Pet127p. tRNAs and, to a lesser extent, rRNAs undergo several different base modifications. This complex gene expression system relies on the coordinated action of mitochondrial and nuclear genes and undergoes rapid evolution, contributing to speciation events. Moving beyond the classical model yeast Saccharomyces cerevisiae to other budding yeasts should provide important insights into the coevolution of both genomes that constitute the eukaryotic genetic system.

Distribution and Polymorphisms of Group I Introns in Mitochondrial Genes from Cryptococcus neoformans and Cryptococcus gattii

The species complexes Cryptococcus neoformans and Cryptococcus gattii are the causative agents of cryptococcosis. Virulence and susceptibility to antifungals may vary within each species according to the fungal genotype. Therefore, specific and easily accessible molecular markers are required to distinguish cryptic species and/or genotypes. Group I introns are potential markers for this purpose because they are polymorphic concerning their presence and sequence. Therefore, in this study, we evaluated the presence of group I introns in the mitochondrial genes cob and cox1 in different Cryptococcus isolates. Additionally, the origin, distribution, and evolution of these introns were investigated by phylogenetic analyses, including previously sequenced introns for the mtLSU gene. Approximately 80.5% of the 36 sequenced introns presented homing endonucleases, and phylogenetic analyses revealed that introns occupying the same insertion site form monophyletic clades. This suggests that they likely share a common ancestor that invaded the site prior to species divergence. There was only one case of heterologous invasion, probably through horizontal transfer to C. decagattii (VGIV genotype) from another fungal species. Our results showed that the C. neoformans complex has fewer introns compared to C. gattii. Additionally, there is significant polymorphism in the presence and size of these elements, both among and within genotypes. As a result, it is impossible to differentiate the cryptic species using a single intron. However, it was possible to differentiate among genotypes within each species complex, by combining PCRs of mtLSU and cox1 introns, for C. neoformans species, and mtLSU and cob introns for C. gattii species.

An overview of RNA splicing and functioning of splicing factors in land plant chloroplasts

RNA splicing refers to a process by which introns of a pre-mRNA are excised and the exons at both ends are joined together. Chloroplast introns are inherently self-splicing ribozymes, but over time, they have lost self-splicing ability due to the degeneration of intronic elements. Thus, the splicing of chloroplast introns relies heavily on nuclear-encoded splicing factors, which belong to diverse protein families. Different splicing factors and their shared intron targets are supposed to form ribonucleoprotein particles (RNPs) to facilitate intron splicing. As characterized in a previous review, around 14 chloroplast intron splicing factors were identified until 2010. However, only a few genetic and biochemical evidence has shown that these splicing factors are required for the splicing of one or several introns. The roles of splicing factors are generally believed to facilitate intron folding; however, the precise role of each protein in RNA splicing remains ambiguous. This may be because the precise binding site of most of these splicing factors remains unexplored. In the last decade, several new splicing factors have been identified. Also, several splicing factors were found to bind to specific sequences within introns, which enhanced the understanding of splicing factors. Here, we summarize recent progress on the splicing factors in land plant chloroplasts and discuss their possible roles in chloroplast RNA splicing based on previous studies.

Global Characterization of Fungal Mitogenomes: New Insights on Genomic Diversity and Dynamism of Coding Genes and Accessory Elements

Fungi comprise a great diversity of species with distinct ecological functions and lifestyles. Similar to other eukaryotes, fungi rely on interactions with prokaryotes and one of the most important symbiotic events was the acquisition of mitochondria. Mitochondria are organelles found in eukaryotic cells whose main function is to generate energy through aerobic respiration. Mitogenomes (mtDNAs) are double-stranded circular or linear DNA from mitochondria that may contain core genes and accessory elements that can be replicated, transcribed, and independently translated from the nuclear genome. Despite their importance, investigative studies on the diversity of fungal mitogenomes are scarce. Herein, we have evaluated 788 curated fungal mitogenomes available at NCBI database to assess discrepancies and similarities among them and to better understand the mechanisms involved in fungal mtDNAs variability. From a total of 12 fungal phyla, four do not have any representative with available mitogenomes, which highlights the underrepresentation of some groups in the current available data. We selected representative and non-redundant mitogenomes based on the threshold of 90% similarity, eliminating 81 mtDNAs. Comparative analyses revealed considerable size variability of mtDNAs with a difference of up to 260 kb in length. Furthermore, variation in mitogenome length and genomic composition are generally related to the number and length of accessory elements (introns, HEGs, and uORFs). We identified an overall average of 8.0 (0–39) introns, 8.0 (0–100) HEGs, and 8.2 (0–102) uORFs per genome, with high variation among phyla. Even though the length of the core protein-coding genes is considerably conserved, approximately 36.3% of the mitogenomes evaluated have at least one of the 14 core coding genes absent. Also, our results revealed that there is not even a single gene shared among all mitogenomes. Other unusual genes in mitogenomes were also detected in many mitogenomes, such as dpo and rpo, and displayed diverse evolutionary histories. Altogether, the results presented in this study suggest that fungal mitogenomes are diverse, contain accessory elements and are absent of a conserved gene that can be used for the taxonomic classification of the Kingdom Fungi.

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.

The mitochondrial genome of the grape powdery mildew pathogen Erysiphe necator is intron rich and exhibits a distinct gene organization

Powdery mildews are notorious fungal plant pathogens but only limited information exists on their genomes. Here we present the mitochondrial genome of the grape powdery mildew fungus Erysiphe necator and a high-quality mitochondrial gene annotation generated through cloning and Sanger sequencing of full-length cDNA clones. The E. necator mitochondrial genome consists of a circular DNA sequence of 188,577 bp that harbors a core set of 14 protein-coding genes that are typically present in fungal mitochondrial genomes, along with genes encoding the small and large ribosomal subunits, a ribosomal protein S3, and 25 mitochondrial-encoded transfer RNAs (mt-tRNAs). Interestingly, it also exhibits a distinct gene organization with atypical bicistronic-like expression of the nad4L/nad5 and atp6/nad3 gene pairs, and contains a large number of 70 introns, making it one of the richest in introns mitochondrial genomes among fungi. Sixty-four intronic ORFs were also found, most of which encoded homing endonucleases of the LAGLIDADG or GIY-YIG families. Further comparative analysis of five E. necator isolates revealed 203 polymorphic sites, but only five were located within exons of the core mitochondrial genes. These results provide insights into the organization of mitochondrial genomes of powdery mildews and represent valuable resources for population genetic and evolutionary studies.

RNA-Targeting Splicing Modifiers: Drug Development and Screening Assays

RNA splicing is an essential step in producing mature messenger RNA (mRNA) and other RNA species. Harnessing RNA splicing modifiers as a new pharmacological modality is promising for the treatment of diseases caused by aberrant splicing. This drug modality can be used for infectious diseases by disrupting the splicing of essential pathogenic genes. Several antisense oligonucleotide splicing modifiers were approved by the U.S. Food and Drug Administration (FDA) for the treatment of spinal muscular atrophy (SMA) and Duchenne muscular dystrophy (DMD). Recently, a small-molecule splicing modifier, risdiplam, was also approved for the treatment of SMA, highlighting small molecules as important warheads in the arsenal for regulating RNA splicing. The cellular targets of these approved drugs are all mRNA precursors (pre-mRNAs) in human cells. The development of novel RNA-targeting splicing modifiers can not only expand the scope of drug targets to include many previously considered “undruggable” genes but also enrich the chemical-genetic toolbox for basic biomedical research. In this review, we summarized known splicing modifiers, screening methods for novel splicing modifiers, and the chemical space occupied by the small-molecule splicing modifiers.

The mitochondrial genome of Ophiostoma himal-ulmi and comparison with other fungi causing Dutch elm disease

The mitochondrial genome of Ophiostoma himal-ulmi, a species endemic to the Western Himalayas and one of the fungi that cause Dutch elm disease, has been sequenced and characterized. The mitochondrial genome was compared with other available genomes for members of the Ophiostomatales, including other agents of Dutch elm disease (Ophiostoma ulmi, Ophiostoma novo-ulmi subspecies novo-ulmi, and Ophiostoma novo-ulmi subspecies americana), and it was observed that gene synteny is highly conserved, and variability among members of the fungi that cause Dutch-elm disease is primarily due to the number of intron insertions. Among the fungi that cause Dutch elm disease that we examined, O. himal-ulmi has the largest mitochondrial genomes (ranging from 94 934 to 111 712 bp), owing to the expansion of the number of introns.

The Mitogenomes of Ophiostoma minus and Ophiostoma piliferum and Comparisons With Other Members of the Ophiostomatales

Fungi assigned to the Ophiostomatales are of economic concern as many are blue-stain fungi and some are plant pathogens. The mitogenomes of two blue-stain fungi, Ophiostoma minus and Ophiostoma piliferum, were sequenced and compared with currently available mitogenomes for other members of the Ophiostomatales. Species representing various genera within the Ophiostomatales have been examined for gene content, gene order, phylogenetic relationships, and the distribution of mobile elements. Gene synteny is conserved among the Ophiostomatales but some members were missing the atp9 gene. A genome wide intron landscape has been prepared to demonstrate the distribution of the mobile genetic elements (group I and II introns and homing endonucleases) and to provide insight into the evolutionary dynamics of introns among members of this group of fungi. Examples of complex introns or nested introns composed of two or three intron modules have been observed in some species. The size variation among the mitogenomes (from 23.7 kb to about 150 kb) is mostly due to the presence and absence of introns. Members of the genus Sporothrix sensu stricto appear to have the smallest mitogenomes due to loss of introns. The taxonomy of the Ophiostomatales has recently undergone considerable revisions; however, some lineages remain unresolved. The data showed that genera such as Raffaelea appear to be polyphyletic and the separation of Sporothrix sensu stricto from Ophiostoma is justified.

Exploring the Relationship Among Divergence Time and Coding and Non-coding Elements in the Shaping of Fungal Mitochondrial Genomes

The order Hypocreales (Ascomycota) is composed of ubiquitous and ecologically diverse fungi such as saprobes, biotrophs, and pathogens. Despite their phylogenetic relationship, these species exhibit high variability in biomolecules production, lifestyle, and fitness. The mitochondria play an important role in the fungal biology, providing energy to the cells and regulating diverse processes, such as immune response. In spite of its importance, the mechanisms that shape fungal mitogenomes are still poorly understood. Herein, we investigated the variability and evolution of mitogenomes and its relationship with the divergence time using the order Hypocreales as a study model. We sequenced and annotated for the first time Trichoderma harzianum mitochondrial genome (mtDNA), which was compared to other 34 mtDNAs species that were publicly available. Comparative analysis revealed a substantial structural and size variation on non-coding mtDNA regions, despite the conservation of copy number, length, and structure of protein-coding elements. Interestingly, we observed a highly significant correlation between mitogenome length, and the number and size of non-coding sequences in mitochondrial genome. Among the non-coding elements, group I and II introns and homing endonucleases genes (HEGs) were the main contributors to discrepancies in mitogenomes structure and length. Several intronic sequences displayed sequence similarity among species, and some of them are conserved even at gene position, and were present in the majority of mitogenomes, indicating its origin in a common ancestor. On the other hand, we also identified species-specific introns that advocate for the origin by different mechanisms. Investigation of mitochondrial gene transfer to the nuclear genome revealed that nuclear copies of the nad5 are the most frequent while atp8, atp9, and cox3 could not be identified in any of the nuclear genomes analyzed. Moreover, we also estimated the divergence time of each species and investigated its relationship with coding and non-coding elements as well as with the length of mitogenomes. Altogether, our results demonstrated that introns and HEGs are key elements on mitogenome shaping and its presence on fast-evolving mtDNAs could be mostly explained by its divergence time, although the intron sharing profile suggests the involvement of other mechanisms on the mitochondrial genome evolution, such as horizontal transference.

The fungal mitochondrial Nad5 pan-genic intron landscape

An intron landscape was prepared for the fungal mitochondrial nad5 gene. A hundred and eighty-eight fungal species were examined and a total of 265 introns were noted to be located in 29 intron insertion sites within the examined nad5 genes. Two hundred and sixty-three introns could be classified as group I types and two group II introns were noted. One additional group II intron module was identified nested within a composite group I intron. Based on features related to RNA secondary structures, introns can be classified into different subtypes and it was observed that intron insertion-sites are biased towards phase 0 and they appear to be specific to an intron type. Intron landscapes could be used as a guide map to predict the location of fungal mtDNA mobile introns, which are composite elements that include a ribozyme component and in some instances open reading frames encoding homing endonucleases or reverse transcriptases and all of these have applications in biotechnology.

Genome streamlining via complete loss of introns has occurred multiple times in lichenized fungal mitochondria

Reductions in genome size and complexity are a hallmark of obligate symbioses. The mitochondrial genome displays clear examples of these reductions, with the ancestral alpha-proteobacterial genome size and gene number having been reduced by orders of magnitude in most descendent modern mitochondrial genomes. Here, we examine patterns of mitochondrial evolution specifically looking at intron size, number, and position across 58 species from 21 genera of lichenized Ascomycete fungi, representing a broad range of fungal diversity and niches. Our results show that the cox1gene always contained the highest number of introns out of all the mitochondrial protein-coding genes, that high intron sequence similarity (>90%) can be maintained between different genera, and that lichens have undergone at least two instances of complete, genome-wide intron loss consistent with evidence for genome streamlining via loss of parasitic, noncoding DNA, in Phlyctis boliviensisand Graphis lineola. Notably, however, lichenized fungi have not only undergone intron loss but in some instances have expanded considerably in size due to intron proliferation (e.g., Alectoria fallacina and Parmotrema neotropicum), even between closely related sister species (e.g., Cladonia). These results shed light on the highly dynamic mitochondrial evolution that is occurring in lichens and suggest that these obligate symbiotic organisms are in some cases undergoing recent, broad-scale genome streamlining via loss of protein-coding genes as well as noncoding, parasitic DNA elements.

Acetonic Fraction of Bidens pilosa Enriched for Maturase K Is Able to Control Cerebral Parasite Burden in Mice Experimentally Infected With Toxoplasma gondii

Toxoplasma gondii infection can cause abortions or congenital infection for a vast number of domestic animals and humans, leading to economic loss in veterinary sciences, as well as severe consequences for immunocompromised patients. Bidens pilosa Linné has been used in ethnopharmacology for treatment of diseases, as malaria, diabetes and hepatitis, in addition to its use as antioxidant, antiallergic, anti-inflammatory, and antiviral. The components of this plant have never been studied before for treatment of toxoplasmosis, and the conventional drugs currently used to treat this disease have high degree of toxicity. Thus, the aim of this study was to evaluate the effect of B. pilosa against T. gondii, by analyzing a total extract of this plant in parallel with a fraction obtained by precipitation in acetone. Also, it was assessed if the acetonic fraction could present lectinic activity, followed by its identification by mass spectrometry. It was observed with the experimental models designed that both total extract and acetonic fraction of B. pilosa were able to control T. gondii infection by in vitro and in vivo experiments, in addition to their low toxicity to host cells. Both total extract and acetonic fraction of this plant display capacity to impair replication of T. gondii tachyzoites. Interesting, the B. pilosa acetonic fraction treatment for 10 days after infection decreases significantly the number of T. gondii brain cyst in comparison with controls. The protein isolated from B. pilosa acetonic fraction was characterized as a novel lectin identified as maturase K. Taken together, these findings open new perspectives to treat patients infected by T. gondii. Future studies will be necessary to investigate the precise mechanism underlying the control of T. gondii infection to impair the replication of this parasite in the host cells after treatment with B. pilosa maturase K.

The complete mitochondrial genome of the Chan-hua fungus Isaria cicadae: a tale of intron evolution in Cordycipitaceae

Isaria cicadae is an entomogenous fungus of great medicinal value. Its nuclear genome has been reported, while its mitogenome remains unknown. Herein, we first described its mitogenome and then inferred intron evolution from both intraspecific and interspecific perspectives. The fungus represented the largest mitogenome (56.6 kb in strain CCAD02) known in Cordycipitaceae due to the presence of 25 introns interrupting nine genes. Comparison of three I. cicadae strains revealed intron presence/absence dynamics at six intron loci plus a few indels and single nucleotide polymorphisms. Phylogenetic analyses confirmed the placement of I. cicadae in Cordycipitaceae. Comparison of 10 Cordycipitaceae species revealed a high degree of synteny and conserved genetic content. They, however, varied in intron numbers (1-25 per species) with overall 34 intron loci identified, which resulted in more than twofold variations in mitogenome sizes (24.5-56.6 kb). An rnl intron encoding ribosomal protein S3 was present in all species, suggesting its early invasion in Cordycipitaceae, while further divergence occurred for this intron. The other introns identified in this study were present in some, but not all of the species and have undergone multiple gains and losses in Cordycipitaceae. This study greatly enhanced our understanding of intron evolution in Cordycipitaceae.

The mitochondrial genome of Endoconidiophora resinifera is intron rich

Endoconidiophora resinifera (=Ceratocystis resinifera) is a blue-stain fungus that occurs on conifers. The data showed that the Endoconidiophora resinifera mitochondrial genome is one of the largest mitochondrial genomes (>220 kb) so far reported among members of the Ascomycota. An exceptional large number of introns (81) were noted and differences among the four strains were restricted to minor variations in intron numbers and a few indels and single nucleotide polymorphisms. The major differences among the four strains examined are due to size polymorphisms generated by the absence or presence of mitochondrial introns. Also, these mitochondrial genomes encode the largest cytochrome oxidase subunit 1 gene (47.5 kb) reported so far among the fungi. The large size for this gene again can be attributed to the large number of intron insertions. This study reports the first mitochondrial genome for the genus Endoconidiophora, previously members of this genus were assigned to Ceratocystis. The latter genus has recently undergone extensive taxonomic revisions and the mitochondrial genome might provide loci that could be applied as molecular markers assisting in the identification of taxa within this group of economically important fungi. The large mitochondrial genome also may provide some insight on mechanisms that can lead to mitochondrial genome expansion.

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.

Population structure of mitochondrial genomes in Saccharomyces cerevisiae

Background

Rigorous study of mitochondrial functions and cell biology in the budding yeast, Saccharomyces cerevisiae has advanced our understanding of mitochondrial genetics. This yeast is now a powerful model for population genetics, owing to large genetic diversity and highly structured populations among wild isolates. Comparative mitochondrial genomic analyses between yeast species have revealed broad evolutionary changes in genome organization and architecture. A fine-scale view of recent evolutionary changes within S. cerevisiae has not been possible due to low numbers of complete mitochondrial sequences.

Results

To address challenges of sequencing AT-rich and repetitive mitochondrial DNAs (mtDNAs), we sequenced two divergent S. cerevisiae mtDNAs using a single-molecule sequencing platform (PacBio RS). Using de novo assemblies, we generated highly accurate complete mtDNA sequences. These mtDNA sequences were compared with 98 additional mtDNA sequences gathered from various published collections. Phylogenies based on mitochondrial coding sequences and intron profiles revealed that intraspecific diversity in mitochondrial genomes generally recapitulated the population structure of nuclear genomes. Analysis of intergenic sequence indicated a recent expansion of mobile elements in certain populations. Additionally, our analyses revealed that certain populations lacked introns previously believed conserved throughout the species, as well as the presence of introns never before reported in S. cerevisiae.

Conclusions

Our results revealed that the extensive variation in S. cerevisiae mtDNAs is often population specific, thus offering a window into the recent evolutionary processes shaping these genomes. In addition, we offer an effective strategy for sequencing these challenging AT-rich mitochondrial genomes for small scale projects.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1664-4) contains supplementary material, which is available to authorized users.

A survey of transposable element classification systems--a call for a fundamental update to meet the challenge of their diversity and complexity

The increase of publicly available sequencing data has allowed for rapid progress in our understanding of genome composition. As new information becomes available we should constantly be updating and reanalyzing existing and newly acquired data. In this report we focus on transposable elements (TEs) which make up a significant portion of nearly all sequenced genomes. Our ability to accurately identify and classify these sequences is critical to understanding their impact on host genomes. At the same time, as we demonstrate in this report, problems with existing classification schemes have led to significant misunderstandings of the evolution of both TE sequences and their host genomes. In a pioneering publication Finnegan (1989) proposed classifying all TE sequences into two classes based on transposition mechanisms and structural features: the retrotransposons (class I) and the DNA transposons (class II). We have retraced how ideas regarding TE classification and annotation in both prokaryotic and eukaryotic scientific communities have changed over time. This has led us to observe that: (1) a number of TEs have convergent structural features and/or transposition mechanisms that have led to misleading conclusions regarding their classification, (2) the evolution of TEs is similar to that of viruses by having several unrelated origins, (3) there might be at least 8 classes and 12 orders of TEs including 10 novel orders. In an effort to address these classification issues we propose: (1) the outline of a universal TE classification, (2) a set of methods and classification rules that could be used by all scientific communities involved in the study of TEs, and (3) a 5-year schedule for the establishment of an International Committee for Taxonomy of Transposable Elements (ICTTE).

Intronic and plasmid-derived regions contribute to the large mitochondrial genome sizes of Agaricomycetes

Sizes of mitochondrial genomes vary extensively between fungal species although they typically contain a conserved set of core genes. We have characterised the mitochondrial genome of the conifer root rot pathogen Heterobasidion irregulare and compared the size, gene content and structure of 20 Basidiomycete mitochondrial genomes. The mitochondrial genome of H. irregulare was 114, 193 bp and contained a core set of 15 protein coding genes, two rRNA genes and 26 tRNA genes. In addition, we found six non-conserved open reading frames (ORFs) and four putative plasmid genes clustered in three separate regions together with 24 introns and 14 intronic homing endonuclease genes, unequally spread across seven of the core genes. The size differences among the 20 Basidiomycetes can largely be explained by length variation of intergenic regions and introns. The Agaricomycetes contained the nine largest mitochondrial genomes in the Basidiomycete group and Agaricomycete genomes are significantly (p < 0.001) larger than the other Basidiomycetes. A feature of the Agaricomycete mitochondrial genomes in this study was the simultaneous occurrence of putative plasmid genes and non-conserved ORFs, with Cantharellus cibarius as only exception, where no non-conserved ORF was identified. This indicates a mitochondrial plasmid origin of the non-conserved ORFs or increased mitochondrial genome dynamics of species harbouring mitochondrial plasmids. We hypothesise that two independent factors are the driving forces for large mitochondrial genomes: the homing endonuclease genes in introns and integration of plasmid DNA.

Bacterial group I introns: mobile RNA catalysts

Group I introns are intervening sequences that have invaded tRNA, rRNA and protein coding genes in bacteria and their phages. The ability of group I introns to self-splice from their host transcripts, by acting as ribozymes, potentially renders their insertion into genes phenotypically neutral. Some group I introns are mobile genetic elements due to encoded homing endonuclease genes that function in DNA-based mobility pathways to promote spread to intronless alleles. Group I introns have a limited distribution among bacteria and the current assumption is that they are benign selfish elements, although some introns and homing endonucleases are a source of genetic novelty as they have been co-opted by host genomes to provide regulatory functions. Questions regarding the origin and maintenance of group I introns among the bacteria and phages are also addressed.

Homing endonucleases: from genetic anomalies to programmable genomic clippers

Homing endonucleases are strong drivers of genetic exchange and horizontal transfer of both their own genes and their local genetic environment. The mechanisms that govern the function and evolution of these genetic oddities have been well documented over the past few decades at the genetic, biochemical, and structural levels. This wealth of information has led to the manipulation and reprogramming of the endonucleases and to their exploitation in genome editing for use as therapeutic agents, for insect vector control and in agriculture. In this chapter we summarize the molecular properties of homing endonucleases and discuss their strengths and weaknesses in genome editing as compared to other site-specific nucleases such as zinc finger endonucleases, TALEN, and CRISPR-derived endonucleases.

Two independent activities define Ccm1p as a moonlighting protein in Saccharomyces cerevisiae

Ccm1p is a nuclear-encoded PPR (pentatricopeptide repeat) protein that localizes into mitochondria of Saccharomyces cerevisiae. It was first defined as an essential factor to remove the bI4 [COB (cytochrome b) fourth intron)] and aI4 [COX1 (cytochrome c oxidase subunit 1) fourth intron] of pre-mRNAs, along with bI4 maturase, a protein encoded by part of bI4 and preceding exons that removes the intronic RNA sequence that codes for it. Later on, Ccm1p was described as key to maintain the steady-state levels of the mitoribosome small subunit RNA (15S rRNA). bI4 maturase is produced inside the mitochondria and therefore its activity depends on the functionality of mitochondrial translation. This report addresses the dilemma of whether Ccm1p supports bI4 maturase activity by keeping steady-state levels of 15S rRNA or separately and directly supports bI4 maturase activity per se. Experiments involving loss of Ccm1p, SMDC (sudden mitochondrial deprivation of Ccm1p) and mutations in one of the PPR (pentatricopeptide repeat) motifs revealed that the failure of bI4 maturase activity in CCM1 deletion mutants was not due to a malfunction of the translational machinery. Both functions were found to be independent, defining Ccm1p as a moonlighting protein. bI4 maturase activity was significantly more dependent on Ccm1p levels than the maintenance of 15S rRNA. The novel strategy of SMDC described here allowed the study of immediate short-term effects, before the mutant phenotype was definitively established. This approach can be also applied for further studies on 15S rRNA stability and mitoribosome assembly.

The brown algae Pl.LSU/2 group II intron-encoded protein has functional reverse transcriptase and maturase activities

Group II introns are self-splicing mobile elements found in prokaryotes and eukaryotic organelles. These introns propagate by homing into precise genomic locations, following assembly of a ribonucleoprotein complex containing the intron-encoded protein (IEP) and the spliced intron RNA. Engineered group II introns are now commonly used tools for targeted genomic modifications in prokaryotes but not in eukaryotes. We speculate that the catalytic activation of currently known group II introns is limited in eukaryotic cells. The brown algae Pylaiella littoralis Pl.LSU/2 group II intron is uniquely capable of in vitro ribozyme activity at physiological level of magnesium but this intron remains poorly characterized. We purified and characterized recombinant Pl.LSU/2 IEP. Unlike most IEPs, Pl.LSU/2 IEP displayed a reverse transcriptase activity without intronic RNA. The Pl.LSU/2 intron could be engineered to splice accurately in Saccharomyces cerevisiae and splicing efficiency was increased by the maturase activity of the IEP. However, spliced transcripts were not expressed. Furthermore, intron splicing was not detected in human cells. While further tool development is needed, these data provide the first functional characterization of the PI.LSU/2 IEP and the first evidence that the Pl.LSU/2 group II intron splicing occurs in vivo in eukaryotes in an IEP-dependent manner.

The mitochondrial LSU rRNA group II intron of Ustilago maydis encodes an active homing endonuclease likely involved in intron mobility

Background

The a2 mating type locus gene lga2 is critical for uniparental mitochondrial DNA inheritance during sexual development of Ustilago maydis. Specifically, the absence of lga2 results in biparental inheritance, along with efficient transfer of intronic regions in the large subunit rRNA gene between parental molecules. However, the underlying role of the predicted LAGLIDADG homing endonuclease gene I-UmaI located within the group II intron LRII1 has remained unresolved.

Methodology/Principal Findings

We have investigated the enzymatic activity of I-UmaI in vitro based on expression of a tagged full-length and a naturally occurring mutant derivative, which harbors only the N-terminal LAGLIDADG domain. This confirmed Mg²⁺-dependent endonuclease activity and cleavage at the LRII1 insertion site to generate four base pair extensions with 3′ overhangs. Specifically, I-UmaI recognizes an asymmetric DNA sequence with a minimum length of 14 base pairs (5′-GACGGGAAGACCCT-3′) and tolerates subtle base pair substitutions within the homing site. Enzymatic analysis of the mutant variant indicated a correlation between the activity in vitro and intron homing. Bioinformatic analyses revealed that putatively functional or former functional I-UmaI homologs are confined to a few members within the Ustilaginales and Agaricales, including the phylogenetically distant species Lentinula edodes, and are linked to group II introns inserted into homologous positions in the LSU rDNA.

Conclusions/Significance

The present data provide strong evidence that intron homing efficiently operates under conditions of biparental inheritance in U. maydis. Conversely, uniparental inheritance may be critical to restrict the transmission of mobile introns. Bioinformatic analyses suggest that I-UmaI-associated introns have been acquired independently in distant taxa and are more widespread than anticipated from available genomic data.

A complete sequence of Saccharomyces paradoxus mitochondrial genome that restores the respiration in S. cerevisiae

We determined the complete sequence of 71 355-bp-long mitochondrial genome from Saccharomyces paradoxus entirely by direct sequencing of purified mitochondrial DNA (mtDNA). This mtDNA possesses the same features as its close relative Saccharomyces cerevisiae - A + T content 85.9%, set of genes coding for the three components of cytochrome oxidase, cytochrome b, three subunits of ATPase, both ribosomal subunits, gene for ribosomal protein, rnpB gene, tRNA package (24) and yeast genetic code. Genes are interrupted by nine group I and group II introns, two of which are in positions unknown in S. cerevisiae, but recognized in Saccharomyces pastorianus. The gene products are related to S. cerevisiae, and the identity of amino acid residues varies from 100% for cox2 to 83% for rps3. The remarkable differences from S. cerevisiae are (1) different gene order (translocation of trnF-trnT1-trnV-cox3-trnfM-rnpb-trnP and transposition of trnW-rns), (2) occurrence of two unusual GI introns, (3) eight active ori elements, and (4) reduced number of GC clusters and divergent intergenic spacers. Despite these facts, the sequenced S. paradoxus mtDNA introduced to S. cerevisiae was able to support the respiratory function to the same extent as the original mtDNAs.

Homing endonucleases: DNA scissors on a mission

Buried within the genomes of many microorganisms are genetic elements that encode rare-cutting homing endonucleases that assist in the mobility of the elements that encode them, such as the self-splicing group I and II introns and in some cases inteins. There are several different families of homing endonucleases and their ability to initiate and target specific sequences for lateral transfers makes them attractive reagents for gene targeting. Homing endonucleases have been applied in promoting DNA modification or genome editing such as gene repair or "gene knockouts". This review examines the categories of homing endonucleases that have been described so far and their possible applications to biotechnology. Strategies to engineer homing endonucleases to alter target site specificities will also be addressed. Alternatives to homing endonucleases such as zinc finger nucleases, transcription activator-like effector nucleases, triplex forming oligonucleotide nucleases, and targetrons are also briefly discussed.

Phosphatase activity in barley proteins tightly bound to DNA and its development-dependent changes

The tightly bound proteins (TBPs), a protein group that remains attached to DNA either covalently or noncovalently after deproteinization, have been found in numerous eukaryotic species. Some TBPs isolated from mammalian and yeast cells possess phosphatase or kinase activity. The aim of this study was to characterize further TBPs in barley (Hordeum vulgare) cells. The spectra of TBPs varied in different organs of barley shoots (first leaves, coleoptile, and roots) and at different developmental stages of the plant. Some barley TBPs manifested phosphatase, probably Ser/Thr or dual Ser/Thr/Tyr activity. MALDI-TOF mass spectrometry of barley TBPs identified several proteins involved in chromatin rearrangement and regulation processes, including transcription factors, serpins, protein phosphatases and protein kinases, RNA helicases, and DNA topoisomerase II.

Learning to live together: mutualism between self-splicing introns and their hosts

Group I and II introns can be considered as molecular parasites that interrupt protein-coding and structural RNA genes in all domains of life. They function as self-splicing ribozymes and thereby limit the phenotypic costs associated with disruption of a host gene while they act as mobile DNA elements to promote their spread within and between genomes. Once considered purely selfish DNA elements, they now seem, in the light of recent work on the molecular mechanisms regulating bacterial and phage group I and II intron dynamics, to show evidence of co-evolution with their hosts. These previously underappreciated relationships serve the co-evolving entities particularly well in times of environmental stress.

Evolutionary dynamics of the mS952 intron: a novel mitochondrial group II intron encoding a LAGLIDADG homing endonuclease gene

Examination of the mitochondrial small subunit ribosomal RNA (rns) gene of five species of the fungal genus Leptographium revealed that the gene has been invaded at least once at position 952 by a group II intron encoding a LAGLIDADG homing endonuclease gene. Phylogenetic analyses of the intron and homing endonuclease sequences indicated that each element in Leptographium species forms a single clade and is closely related to the group II intron/homing endonuclease gene composite element previously reported at position 952 of the mitochondrial rns gene of Cordyceps species and of Cryphonectria parasitica. The results of an intron survey of the mt rns gene of Leptographium species superimposed onto the phylogenetic analysis of the host organisms suggest that the composite element was transmitted vertically in Leptographium lundbergii. However, its stochastic distribution among strains of L. wingfieldii, L. terebrantis, and L. truncatum suggests that it has been horizontally transmitted by lateral gene transfer among these species, although the random presence of the intron may reflect multiple random loss events. A model is proposed describing the initial invasion of the group II intron in the rns gene of L. lundbergii by a LAGLIDADG homing endonuclease gene and subsequent evolution of this gene to recognize a novel DNA target site, which may now promote the mobility of the intron and homing endonuclease gene as a composite element.

A maturase that specifically stabilizes and activates its cognate group I intron at high temperatures

Folding of large structured RNAs into their functional tertiary structures at high temperatures is challenging. Here we show that I-TnaI protein, a small LAGLIDADG homing endonuclease encoded by a group I intron from a hyperthermophilic bacterium, acts as a maturase that is essential for the catalytic activity of this intron at high temperatures and physiological cationic conditions. I-TnaI specifically binds to and induces tertiary packing of the P4-P6 domain of the intron; this RNA-protein complex might serve as a thermostable platform for active folding of the entire intron. Interestingly, the binding affinity of I-TnaI to its cognate intron RNA largely increases with temperature; over 30-fold stronger binding at higher temperatures relative to 37 °C correlates with a switch from an entropy-driven (37 °C) to an enthalpy-driven (55-60 °C) interaction mode. This binding mode may represent a novel strategy how an RNA binding protein can promote the function of its target RNA specifically at high temperatures.

Diversity of sponge mitochondrial introns revealed by cox 1 sequences of Tetillidae

Background

Animal mitochondrial introns are rare. In sponges and cnidarians they have been found in the cox 1 gene of some spirophorid and homosclerophorid sponges, as well as in the cox 1 and nad 5 genes of some Hexacorallia. Their sporadic distribution has raised a debate as to whether these mobile elements have been vertically or horizontally transmitted among their hosts. The first sponge found to possess a mitochondrial intron was a spirophorid sponge from the Tetillidae family. To better understand the mode of transmission of mitochondrial introns in sponges, we studied cox 1 intron distribution among representatives of this family.

Results

Seventeen tetillid cox 1 sequences were examined. Among these sequences only six were found to possess group I introns. Remarkably, three different forms of introns were found, named introns 714, 723 and 870 based on their different positions in the cox 1 alignment. These introns had distinct secondary structures and encoded LAGLIDADG ORFs belonging to three different lineages. Interestingly, sponges harboring the same intron form did not always form monophyletic groups, suggesting that their introns might have been transferred horizontally. To evaluate whether the introns were vertically or horizontally transmitted in sponges and cnidarians we used a host parasite approach. We tested for co-speciation between introns 723 (the introns with the highest number of sponge representatives) and their nesting cox 1 sequences. Reciprocal AU tests indicated that the intron and cox 1 tree are significantly different, while a likelihood ratio test was not significant. A global test of co-phylogeny had significant results; however, when cnidarian sequences were analyzed separately the results were not significant.

Conclusions

The co-speciation analyses thus suggest that a vertical transmission of introns in the ancestor of sponges and cnidarians, followed by numerous independent losses, cannot solely explain the current distribution of metazoan group I introns. An alternative scenario that includes horizontal gene transfer events appears to be more suitable to explain the incongruence between the intron 723 and the cox 1 topologies. In addition, our results suggest that three different intron forms independently colonized the cox 1 gene of tetillids. Among sponges, the Tetillidae family seems to be experiencing an unusual number of intron insertions.

A group II intron encodes a functional LAGLIDADG homing endonuclease and self-splices under moderate temperature and ionic conditions

A group II intron encoding a protein belonging to the LAGLIDADG family of homing endonucleases was identified in the mitochondrial rns gene of the filamentous fungus Leptographium truncatum, and the catalytic activities of both the intron and its encoded protein were characterized. A model of the RNA secondary structure indicates that the intron is a member of the IIB1 subclass and the open reading frame is inserted in ribozyme domain III. In vitro assays carried out with two versions of the intron, one in which the open reading frame was removed and the other in which it was present, demonstrate that both versions of the intron readily self-splice at 37 degrees C and at a concentration of MgCl(2) as low as 6 mM. The open reading frame encodes a functional LAGLIDADG homing endonuclease that cleaves 2 (top strand) and 6 (bottom strand) nucleotides (nt) upstream of the intron insertion site, generating 4 nt 3' OH overhangs. In vitro splicing assays carried out in the absence and presence of the intron-encoded protein indicate that the protein does not enhance intron splicing, and RNA-binding assays show that the protein does not appear to bind to the intron RNA precursor transcript. These findings raise intriguing questions concerning the functional and evolutionary relationships of the two components of this unique composite element.

The Cryphonectria parasitica mitochondrial rns gene: plasmid-like elements, introns and homing endonucleases

The mt-rns gene of Cryphonectria parasitica is 9872bp long and includes two group I and two group II introns. An analysis of intronic protein-encoding sequences revealed that LAGLIDADG ORFs, which usually are associated with group I introns, were transferred at least twice into group II introns. A plasmid-like mitochondrial element (plME) that appears in high amounts in previously mutagen-induced mit1 and mit2 hypovirulent mutants of the Ep155 standard virulent strain of C. parasitica was found to be derived from a short region of the mt-rns gene, including the exon 1 and most of the first intron. The plME is a 4.2-kb circular, multimeric DNA and an autonomously-replicating mtDNA fragment. Although sexual transmission experiments indicate that the plME does not directly cause hypovirulence, its emergence is one manifestation of the many complex molecular and genetic events that appear to underlie this phenotype.

RNA processing and decay in bacteriophage T4

Bacteriophage T4 is the archetype of virulent phage. It has evolved very efficient strategies to subvert host functions to its benefit and to impose the expression of its genome. T4 utilizes a combination of host and phage-encoded RNases and factors to degrade its mRNAs in a stage-dependent manner. The host endonuclease RNase E is used throughout the phage development. The sequence-specific, T4-encoded RegB endoribonuclease functions in association with the ribosomal protein S1 to functionally inactivate early transcripts and expedite their degradation. T4 polynucleotide kinase plays a role in this process. Later, the viral factor Dmd protects middle and late mRNAs from degradation by the host RNase LS. T4 codes for a set of eight tRNAs and two small, stable RNA of unknown function that may contribute to phage virulence. Their maturation is assured by host enzymes, but one phage factor, Cef, is required for the biogenesis of some of them. The tRNA gene cluster also codes for a homing DNA endonuclease, SegB, responsible for spreading the tRNA genes to other T4-related phage.

Sequence-based prediction of protein-protein interactions by means of codon usage

A new approach based on similarity in codon usage is used to predict protein-protein interactions.

We introduce a novel approach to predict interaction of two proteins solely by analyzing their coding sequences. We found that similarity in codon usage is a strong predictor of protein-protein interactions and, for high specificity values, is as sensitive as the most powerful current prediction methods. Furthermore, combining codon usage with other predictors results in a 75% increase in sensitivity at a precision of 50%, compared to prediction without considering codon usage.

The sporadic occurrence of a group I intron-like element in the mtDNA rnl gene of Ophiostoma novo-ulmi subsp. americana

The presence of group I intron-like elements within the U7 region of the mtDNA large ribosomal subunit RNA gene (rnl) was investigated in strains of Ophiostoma novo-ulmi subsp. americana from Canada, Europe and Eurasia, and in selected strains of O. ips, O. minus, O. piceae, O. ulmi, and O. himal-ulmi. This insertion is of interest as it has been linked previously to the generation of plasmid-like mtDNA elements in diseased strains of O. novo-ulmi. Among 197 O. novo-ulmi subsp. americana strains tested, 61 contained a 1.6kb insertion within the rnl-U7 region and DNA sequence analysis suggests the presence of a group I intron (IA1 type) that encodes a potential double motif LAGLIDADG homing endonuclease-like gene (HEG). Phylogenetic analysis of rnl-U7 intron encoded HEG-like elements supports the view that double motif HEGs originated from a duplication event of a single-motif HEG followed by a fusion event that combined the two copies into one open reading frame (ORF). The data also show that rnl-U7 intron encoded ORFs belong to a clade that includes ORFs inserted into different types of group I introns, e.g. IB, ID, IC3, IA1, present within a variety of different mtDNA genes, such as the small ribosomal subunit RNA gene (rns), apo-cytochrome b gene (cob), NADH dehydrogenase subunit 5 (nad5), cytochrome oxidase subunit 1 gene (coxI), and ATPase subunit 9 gene (atp9). We also compared the occurrence of the rnl-U7 intron in our collection of 227 strains with the presence of the rnl-U11 group I intron and concluded that the U7 intron appears to be an optional element and the U11 intron is probably essential among the strains tested.

Inteins, introns, and homing endonucleases: recent revelations about the life cycle of parasitic genetic elements

Self splicing introns and inteins that rely on a homing endonuclease for propagation are parasitic genetic elements. Their life-cycle and evolutionary fate has been described through the homing cycle. According to this model the homing endonuclease is selected for function only during the spreading phase of the parasite. This phase ends when the parasitic element is fixed in the population. Upon fixation the homing endonuclease is no longer under selection, and its activity is lost through random processes. Recent analyses of these parasitic elements with functional homing endonucleases suggest that this model in its most simple form is not always applicable. Apparently, functioning homing endonuclease can persist over long evolutionary times in populations and species that are thought to be asexual or nearly asexual. Here we review these recent findings and discuss their implications. Reasons for the long-term persistence of a functional homing endonuclease include: More recombination (sexual and as a result of gene transfer) than previously assumed for these organisms; complex population structures that prevent the element from being fixed; a balance between active spreading of the homing endonuclease and a decrease in fitness caused by the parasite in the host organism; or a function of the homing endonuclease that increases the fitness of the host organism and results in purifying selection for the homing endonuclease activity, even after fixation in a local population. In the future, more detailed studies of the population dynamics of the activity and regulation of homing endonucleases are needed to decide between these possibilities, and to determine their relative contributions to the long term survival of parasitic genes within a population. Two outstanding publications on the amoeba Naegleria group I intron (Wikmark et al. BMC Evol Biol 2006, 6:39) and the PRP8 inteins in ascomycetes (Butler et al.BMC Evol Biol 2006, 6:42) provide important stepping stones towards integrated studies on how these parasitic elements evolve through time together with, or despite, their hosts.

Putative cross-kingdom horizontal gene transfer in sponge (Porifera) mitochondria

BACKGROUND

The mitochondrial genome of Metazoa is usually a compact molecule without introns. Exceptions to this rule have been reported only in corals and sea anemones (Cnidaria), in which group I introns have been discovered in the cox1 and nad5 genes. Here we show several lines of evidence demonstrating that introns can also be found in the mitochondria of sponges (Porifera).

RESULTS

A 2,349 bp fragment of the mitochondrial cox1 gene was sequenced from the sponge Tetilla sp. (Spirophorida). This fragment suggests the presence of a 1143 bp intron. Similar to all the cnidarian mitochondrial introns, the putative intron has group I intron characteristics. The intron is present in the cox1 gene and encodes a putative homing endonuclease. In order to establish the distribution of this intron in sponges, the cox1 gene was sequenced from several representatives of the demosponge diversity. The intron was found only in the sponge order Spirophorida. A phylogenetic analysis of the COI protein sequence and of the intron open reading frame suggests that the intron may have been transmitted horizontally from a fungus donor.

CONCLUSION

Little is known about sponge-associated fungi, although in the last few years the latter have been frequently isolated from sponges. We suggest that the horizontal gene transfer of a mitochondrial intron was facilitated by a symbiotic relationship between fungus and sponge. Ecological relationships are known to have implications at the genomic level. Here, an ecological relationship between sponge and fungus is suggested based on the genomic analysis.