Birth of new spliceosomal introns in fungi by multiplication of introner-like elements

Spliceosomal introns are noncoding sequences that separate exons in eukaryotic genes and are removed from pre-messenger RNAs by the splicing machinery. Their origin has remained a mystery in biology since their discovery because intron gains seem to be infrequent in many eukaryotic lineages. Although a few recent intron gains have been reported, none of the proposed gain mechanisms can convincingly explain the high number of introns in present-day eukaryotic genomes. Here we report on particular spliceosomal introns that share high sequence similarity and are reminiscent of introner elements. These elements multiplied in unrelated genes of six fungal genomes and account for the vast majority of intron gains in these fungal species. Such introner-like elements (ILEs) contain all typical characteristics of regular spliceosomal introns (RSIs) but are longer and predicted to harbor more stable secondary structures. However, dating of multiplication events showed that they degenerate in sequence and length within 100,000 years to eventually become indistinguishable from RSIs. We suggest that ILEs not only account for intron gains in six fungi but also in ancestral eukaryotes to give rise to most RSIs by a yet unknown multiplication mechanism.     

Frequent Gain and Loss of Introns in Fungal Cytochrome b Genes

In this study, all available cytochrome b (Cyt b) genes from the GOBASE database were compiled and the evolutionary dynamics of the Cyt b gene introns was assessed. Cyt b gene introns were frequently present in the fungal kingdom and some lower plants, but generally absent or rare in Chromista, Protozoa, and Animalia. Fungal Cyt b introns were found at 35 positions in Cyt b genes and the number of introns varied at individual positions from a single representative to 32 different introns at position 131, showing a wide and patchy distribution. Many homologous introns were present at the same position in distantly related species but absent in closely related species, suggesting that introns of the Cyt b genes were frequently lost. On the other hand, highly similar intron sequences were observed in some distantly related species rather than in closely related species, suggesting that these introns were gained independently, likely through lateral transfers. The intron loss-and-gain events could be mediated by transpositions that might have occurred between nuclear and mitochondria. Southern hybridization analysis confirmed that some introns contained repetitive sequences and might be transposable elements. An intron gain in Botryotinia fuckeliana prevented the development of QoI fungicide resistance, suggesting that intron loss-and-gain events were not necessarily beneficial to their host organisms.     

Analysis of Fungal Genomes Reveals Commonalities of Intron Gain or Loss and Functions in Intron-Poor Species

Previous evolutionary reconstructions have concluded that early eukaryotic ancestors including both the last common ancestor of eukaryotes and of all fungi had intron-rich genomes. By contrast, some extant eukaryotes have few introns, underscoring the complex histories of intron–exon structures, and raising the question as to why these few introns are retained. Here, we have used recently available fungal genomes to address a variety of questions related to intron evolution. Evolutionary reconstruction of intron presence and absence using 263 diverse fungal species supports the idea that massive intron reduction through intron loss has occurred in multiple clades. The intron densities estimated in various fungal ancestors differ from zero to 7.6 introns per 1 kb of protein-coding sequence. Massive intron loss has occurred not only in microsporidian parasites and saccharomycetous yeasts, but also in diverse smuts and allies. To investigate the roles of the remaining introns in highly-reduced species, we have searched for their special characteristics in eight intron-poor fungi. Notably, the introns of ribosome-associated genes RPL7 and NOG2 have conserved positions; both intron-containing genes encoding snoRNAs. Furthermore, both the proteins and snoRNAs are involved in ribosome biogenesis, suggesting that the expression of the protein-coding genes and noncoding snoRNAs may be functionally coordinated. Indeed, these introns are also conserved in three-quarters of fungi species. Our study shows that fungal introns have a complex evolutionary history and underappreciated roles in gene expression.     

Conservation versus parallel gains in intron evolution

Orthologous genes from distant eukaryotic species, e.g. animals and plants, share up to 25–30% intron positions. However, the relative contributions of evolutionary conservation and parallel gain of new introns into this pattern remain unknown. Here, the extent of independent insertion of introns in the same sites (parallel gain) in orthologous genes from phylogenetically distant eukaryotes is assessed within the framework of the protosplice site model. It is shown that protosplice sites are no more conserved during evolution of eukaryotic gene sequences than random sites. Simulation of intron insertion into protosplice sites with the observed protosplice site frequencies and intron densities shows that parallel gain can account but for a small fraction (5–10%) of shared intron positions in distantly related species. Thus, the presence of numerous introns in the same positions in orthologous genes from distant eukaryotes, such as animals, fungi and plants, appears to reflect mostly bona fide evolutionary conservation.     

GROUP I INTRON VERSUS ITS SEQUENCES IN PHYLOGENY OF CETRARIOID LICHENS

Abstract: Phylogenetic trees based on group I intron sequences and on internal transcribed spacer (ITS) sequences of mycobiont ribosomal genes were calculated and compared. Eight cetrarioid and four non-cetrarioid species of the Parmeliaceae were compared. The phylogeny based on group I intron sequences is partly congruent with the ITS sequence phylogeny. Group I intron sequences are presumably less informative for infragenic studies. The introns have a length of 214–233 nucleotides, and differ at up to 33% of the bases between species. All introns analysed are located between the positions 1516 and 1517 of the fungal 18S ribosomal RNA gene. Cetrarioid lichens form a non-homogeneous group within theParmeliaceae according to both group I intron and ITS sequences.     

Molecular evolution of eukaryotic genomes: hemiascomycetous yeast spliceosomal introns

As part of the exploratory sequencing program Génolevures, visual scrutinisation and bioinformatic tools were used to detect spliceosomal introns in seven hemiascomycetous yeast species. A total of 153 putative novel introns were identified. Introns are rare in yeast nuclear genes (<5% have an intron), mainly located at the 5′ end of ORFs, and not highly conserved in sequence. They all share a clear non-random vocabulary: conserved splice sites and conserved nucleotide contexts around splice sites. Homologues of metazoan snRNAs and putative homologues of SR splicing factors were identified, confirming that the spliceosomal machinery is highly conserved in eukaryotes. Several introns’ features were tested as possible markers for phylogenetic analysis. We found that intron sizes vary widely within each genome, and according to the phylogenetic position of the yeast species. The evolutionary origin of spliceosomal introns was examined by analysing the degree of conservation of intron positions in homologous yeast genes. Most introns appeared to exist in the last common ancestor of present day yeast species, and then to have been differentially lost during speciation. However, in some cases, it is difficult to exclude a possible sliding event affecting a pre-existing intron or a gain of a novel intron. Taken together, our results indicate that the origin of spliceosomal introns is complex within a given genome, and that present day introns may have resulted from a dynamic flux between intron conservation, intron loss and intron gain during the evolution of hemiascomycetous yeasts.     

Evidence for extensive recent intron transposition in closely related fungi

Though spliceosomal introns are a major structural component of most eukaryotic genes and intron density varies by more than three orders of magnitude among eukaryotes [1-3], the origins of introns are poorly understood, and only a few cases of unambiguous intron gain are known [4-8]. We utilized population genomic comparisons of three closely related fungi to identify crucial transitory phases of intron gain and loss. We found 74 intron positions showing intraspecific presence-absence polymorphisms (PAPs) for the entire intron. Population genetic analyses identified intron PAPs at different stages of fixation and showed that intron gain or loss was very recent. We found direct support for extensive intron transposition among unrelated genes. A substantial proportion of highly similar introns in the genome either were recently gained or showed a transient phase of intron PAP. We also identified an intron transfer among paralogous genes that created a new intron. Intron loss was due mainly to homologous recombination involving reverse-transcribed mRNA. The large number of intron positions in transient phases of either intron gain or loss shows that intron evolution is much faster than previously thought and provides an excellent model to study molecular mechanisms of intron gain.     

Separate Introns Gained within Short and Long Soluble Peridinin-Chlorophyll a-Protein Genes during Radiation of Symbiodinium (Dinophyceae) Clade A and B Lineages

Here we document introns in two Symbiodinium clades that were most likely gained following divergence of this genus from other peridinin-containing dinoflagellate lineages. Soluble peridinin-chlorophyll a-proteins (sPCP) occur in short and long forms in different species. Duplication and fusion of short sPCP genes produced long sPCP genes. All short and long sPCP genes characterized to date, including those from free living species and Symbiodinium sp. 203 (clade C/type C2) are intronless. However, we observed that long sPCP genes from two Caribbean Symbiodinium clade B isolates each contained two introns. To test the hypothesis that introns were gained during radiation of clade B, we compared sPCP genomic and cDNA sequences from 13 additional distinct Caribbean and Pacific Symbiodinium clade A, B, and F isolates. Long sPCP genes from all clade B/B1 and B/B19 descendants contain orthologs of both introns. Short sPCP genes from S. pilosum (A/A2) and S. muscatinei (B/B4) plus long sPCP genes from S. microadriaticum (A/A1) and S. kawagutii (F/F1) are intronless. Short sPCP genes of S. microadriaticum have a third unique intron. Symbiodinium clade B long sPCP sequences are useful for assessing divergence among B1 and B19 descendants. Phylogenetic analyses of coding sequences from four dinoflagellate orders indicate that introns were gained independently during radiation of Symbiodinium clades A and B. Long sPCP introns were present in the most recent common ancestor of Symbiodinium clade B core types B1 and B19, which apparently diverged sometime during the Miocene. The clade A short sPCP intron was either gained by S. microadriaticum or possibly by the ancestor of Symbiodinium types A/A1, A3, A4 and A5. The timing of short sPCP intron gain in Symbiodinium clade A is less certain. But, all sPCP introns were gained after fusion of ancestral short sPCP genes, which we confirm as occurring once in dinoflagellate evolution.     

125 Changing Environmental Conditions and Changes in the Abundances of Rocky Intertidal Seaweeds on Southern California Shores

In a recent study of the North American biogeography of the red algae genus Hildenbrandia , the presence of group I introns were noted in the nuclear SSU rRNA gene of the marine species H. rubra (Hildenbrandiales). Group I introns in the nuclear encoded rRNAs have been previously reported in the Hildenbrandiales as well as the Bangiales. All reported introns within the red algae have been identified as belonging to the IC1 subclass and occur at two insertion sites in the nuclear small subunit rRNA (516 and 1506). However, an unclassified intron was discovered at position 989 in the nuclear SSU rRNA gene of a collection of H. rubra from British Columbia, Canada. We have determined that the intron is a member of the IE subclass and this is the first report of an IE intron and an intron in position 989 in the red algae. Phylogenetic analyses of the intron sequences reveal a close relationship between this group IE intron and similar ascomycete and basidiomycete fungal IE introns in the nuclear SSU rRNA genes at positions 989 and 1199. In addition, a common unique helix (structural signature) in the P13 domain of the Hildenbrandia intron and those of the fungi at the 989 and 1199 IE positions in the nuclear SSU rRNA gene also indicates a close relationship. Hence, this study provides evidence for a possible lateral transfer of the IE intron in position 989 between fungal and red algal nuclear SSU rRNA genes.     

124 Evidence for Lateral Transfer of a IE Intron between Fungal and Red Algal Small Subunit Ribosomal RNA Genes

In a recent study of the North American biogeography of the red algae genus Hildenbrandia, the presence of group I introns were noted in the nuclear SSU rRNA gene of the marine species H. rubra (Hildenbrandiales). Group I introns in the nuclear encoded rRNAs have been previously reported in the Hildenbrandiales as well as the Bangiales. All reported introns within the red algae have been identified as belonging to the IC1 subclass and occur at two insertion sites in the nuclear small subunit rRNA (516 and 1506). However, an unclassified intron was discovered at position 989 in the nuclear SSU rRNA gene of a collection of H. rubra from British Columbia, Canada. We have determined that the intron is a member of the IE subclass and this is the first report of an IE intron and an intron in position 989 in the red algae. Phylogenetic analyses of the intron sequences reveal a close relationship between this group IE intron and similar ascomycete and basidiomycete fungal IE introns in the nuclear SSU rRNA genes at positions 989 and 1199. In addition, a common unique helix (structural signature) in the P13 domain of the Hildenbrandia intron and those of the fungi at the 989 and 1199 IE positions in the nuclear SSU rRNA gene also indicates a close relationship. Hence, this study provides evidence for a possible lateral transfer of the IE intron in position 989 between fungal and red algal nuclear SSU rRNA genes.     

Evolutionary relics dominate the small number of secondary metabolism genes in the hemibiotrophic fungus Dothistroma septosporum

Fungal secondary metabolites have important functions for the fungi that produce them, such as roles in virulence and competition. The hemibiotrophic pine needle pathogen Dothistroma septosporum has one of the lowest complements of secondary metabolite (SM) backbone genes of plant pathogenic fungi, indicating that this fungus produces a limited range of SMs. Amongst these SMs is dothistromin, a well-characterised polyketide toxin and virulence factor that is required for expansion of disease lesions in Dothistroma needle blight disease. Dothistromin genes are dispersed across six loci on one chromosome, rather than being clustered as for most SM genes. We explored other D. septosporum SM genes to determine if they are associated with gene clusters, and to predict what their likely products and functions might be. Of nine functional SM backbone genes in the D. septosporum genome, only four were expressed under a range of in planta and in culture conditions, one of which was the dothistromin PKS backbone gene. Of the other three expressed genes, gene knockout studies suggested that DsPks1 and DsPks2 are not required for virulence and attempts to determine a functional squalestatin-like SM product for DsPks2 were not successful. However preliminary evidence suggested that DsNps3, the only SM backbone gene to be most highly expressed in the early stage of disease, appears to be a virulence factor. Thus, despite the small number of SM backbone genes in D. septosporum, most of them appear to be poorly expressed or dispensable for virulence in planta. This work contributes to a growing body of evidence that many fungal secondary metabolite gene clusters might be non-functional and may be evolutionary relics.     

Gene invasion in distant eukaryotic lineages: discovery of mutually exclusive genetic elements reveals marine biodiversity

Inteins are rare, translated genetic parasites mainly found in bacteria and archaea, while spliceosomal introns are distinctly eukaryotic features abundant in most nuclear genomes. Using targeted metagenomics, we discovered an intein in an Atlantic population of the photosynthetic eukaryote, Bathycoccus, harbored by the essential spliceosomal protein PRP8 (processing factor 8 protein). Although previously thought exclusive to fungi, we also identified PRP8 inteins in parasitic (Capsaspora) and predatory (Salpingoeca) protists. Most new PRP8 inteins were at novel insertion sites that, surprisingly, were not in the most conserved regions of the gene. Evolutionarily, Dikarya fungal inteins at PRP8 insertion site a appeared more related to the Bathycoccus intein at a unique insertion site, than to other fungal and opisthokont inteins. Strikingly, independent analyses of Pacific and Atlantic samples revealed an intron at the same codon as the Bathycoccus PRP8 intein. The two elements are mutually exclusive and neither was found in cultured Bathycoccus or other picoprasinophyte genomes. Thus, wild Bathycoccus contain one of few non-fungal eukaryotic inteins known and a rare polymorphic intron. Our data indicate at least two Bathycoccus ecotypes exist, associated respectively with oceanic or mesotrophic environments. We hypothesize that intein propagation is facilitated by marine viruses; and, while intron gain is still poorly understood, presence of a spliceosomal intron where a locus lacks an intein raises the possibility of new, intein-primed mechanisms for intron gain. The discovery of nucleus-encoded inteins and associated sequence polymorphisms in uncultivated marine eukaryotes highlights their diversity and reveals potential sexual boundaries between populations indistinguishable by common marker genes.     

Isolation and sequence analysis of a β-tubulin gene from arbuscular mycorrhizal fungi

A full-length β-tubulin gene has been cloned and sequenced from Gigaspora gigantea and Glomus clarum, two arbuscular mycorrhizal fungi (AMF) species in the phylum Glomeromyota. The gene in both species is organized into five exons and four introns. Both genes are 94.9% similar and encode a 447 amino acid protein. In comparison with other fungal groups, the amino acid sequence is most similar to that of fungi in the Chytridiomycota. The codon usage of the gene in both AMF species is broad and biased in favor of an A or a T in the third position. The four introns varied in length from 87 to 168 bp for G. gigantea and from 90 to 136 bp for G. clarum. Of all fungi in which full-length sequences have been published, only AMF do not have an intron before codon 174. The introns positioned at codons 174 and 257 in AMF match the position of different introns in β-tubulin genes of some Zygomycete, Basidiomycete, and Ascomycete fungi. The 5′ and 3′ splice site consensus sequences are similar to those found in introns of most fungi. Sequence analysis from single-strand conformation polymorphism analysis confirmed the presence of two β-tubulin gene copies in G. clarum, but only one copy was evident in G. gigantea based on Southern hybridization analysis.     

The evolutionary gain of spliceosomal introns: sequence and phase preferences

Theories regarding the evolution of spliceosomal introns differ in the extent to which the distribution of introns reflects either a formative role in the evolution of protein-coding genes or the adventitious gain of genetic elements. Here, systematic methods are used to assess the causes of the present-day distribution of introns in 10 families of eukaryotic protein-coding genes comprising 1,868 introns in 488 distinct alignment positions. The history of intron evolution inferred using a probabilistic model that allows ancestral inheritance of introns, gain of introns, and loss of introns reveals that the vast majority of introns in these eukaryotic gene families were not inherited from the most recent common ancestral genes, but were gained subsequently. Furthermore, among inferred events of intron gain that meet strict criteria of reliability, the distribution of sites of gain with respect to reading-frame phase shows a 5:3:2 ratio of phases 0, 1 and 2, respectively, and exhibits a nucleotide preference for MAG GT (positions -3 to +2 relative to the site of gain). The nucleotide preferences of intron gain may prove to be the ultimate cause for the phase bias. The phase bias of intron gain is sufficient to account quantitatively for the well-known 5:3:2 bias in phase frequencies among extant introns, a conclusion that holds even when taxonomic heterogeneity in phase patterns is considered. Thus, intron gain accounts for the vast majority of extant introns and for the bias toward phase 0 introns that previously was interpreted as evidence for ancient formative introns.     

Intron gain and loss in the evolution of the conserved eukaryotic recombination machinery

Intron conservation, intron gain or loss and putative intron sliding events were determined for a set of three genes (SPO11, MRE11 and DMC1) involved in basic aspects of recombination in eukaryotes. These are ancient genes and present in nearly all of the major kingdoms. MRE11 is of bacterial origin and can be found in all kingdoms. DMC1 is a specialized homolog of the bacterial RecA protein, whereas the SPO11 gene is of archaebacterial origin. Only unique homologs of SPO11 are found in animals and fungi whereas three distantly related SPO11 copies are present in plant genomes. A comparison of the respective intron positions and phases of all genes was performed, demonstrating that a quarter of the intron positions were perfectly conserved over more than 1000000000 years. Regarding the remaining three quarters of the introns we found insertions to be about three times more frequent than deletions. Aligning the introns of the three different SPO11 homologs of Arabidopsis thaliana we propose a conclusive model of their evolution. We postulate that at least one duplication event occurred shortly after the divergence of plants from animals and fungi and that a respective homolog has been retained in a protist group, the apicomplexa.     

Fungal origin by horizontal transfer of a plant mitochondrial group I intron in the chimeric coxI gene of Peperomia

We present phylogenetic evidence that a group I intron in an angiosperm mitochondrial gene arose recently by horizontal transfer from a fungal donor species. A 1,716-bp fragment of the mitochondrial coxI gene from the angiosperm Peperomia polybotrya was amplified via the polymerase chain reaction and sequenced. Comparison to other coxI genes revealed a 966-bp group I intron, which, based on homology with the related yeast coxI intron aI4, potentially encodes a 279-amino-acid site-specific DNA endonuclease. This intron, which is believed to function as a ribozyme during its own splicing, is not present in any of 19 coxI genes examined from other diverse vascular plant species. Phylogenetic analysis of intron origin was carried out using three different tree-generating algorithms, and on a variety of nucleotide and amino acid data sets from the intron and its flanking exon sequences. These analyses show that the Peperomia coxI gene intron and exon sequences are of fundamentally different evolutionary origin. The Peperomia intron is more closely related to several fungal mitochondrial introns, two of which are located at identical positions in coxI, than to identically located coxI introns from the land plant Marchantia and the green alga Prototheca. Conversely, the exon sequence of this gene is, as expected, most closely related to other angiosperm coxI genes. These results, together with evidence suggestive of co-conversion of exonic markers immediately flanking the intron insertion site, lead us to conclude that the Peperomia coxI intron probably arose by horizontal transfer from a fungal donor, using the double-strand-break repair pathway. The donor species may have been one of the symbiotic mycorrhizal fungi that live in close obligate association with most plants. Correspondence to: J.C. Vaughn     

Distribution of introns in fungal histone genes

Saccharomycotina and Taphrinomycotina lack intron in their histone genes, except for an intron in one of histone H4 genes of Yarrowia lipolytica. On the other hand, Basidiomycota and Perizomycotina have introns in their histone genes. We compared the distributions of 81, 47, 79, and 98 introns in the fungal histone H2A, H2B, H3, and H4 genes, respectively. Based on the multiple alignments of the amino acid sequences of histones, we identified 19, 13, 31, and 22 intron insertion sites in the histone H2A, H2B, H3, and H4 genes, respectively. Surprisingly only one hot spot of introns in the histone H2A gene is shared between Basidiomycota and Perizomycotina, suggesting that most of introns of Basidiomycota and Perizomycotina were acquired independently. Our findings suggest that the common ancestor of Ascomycota and Basidiomycota maybe had a few introns in the histone genes. In the course of fungal evolution, Saccharomycotina and Taphrinomycotina lost the histone introns; Basidiomycota and Perizomycotina acquired other introns independently. In addition, most of the introns have sequence similarity among introns of phylogenetically close species, strongly suggesting that horizontal intron transfer events between phylogenetically distant species have not occurred recently in the fungal histone genes.     

Tempo and Mode of Spliceosomal Intron Evolution in Actin of Foraminifera

Spliceosomal introns are present in almost all eukaryotic genes, yet little is known about their origin and turnover in the majority of eukaryotic phyla. There is no agreement whether most introns are ancestral and have been lost in some lineage or have been gained recently. We addressed this question by analyzing the spatial and temporal distribution of introns in actins of foraminifera, a group of testate protists whose exceptionally rich fossil record permits the calibration of molecular phylogenies to date intron origins. We identified 24 introns dispersed along the sequence of two foraminiferan actin paralogues and actin deviating proteins, an unconventional type of fast-evolving actin found in some foraminifera. Comparison of intron positions indicates that 20 of 24 introns are specific to foraminifera. Four introns shared between foraminifera and other eukaryotes were interpreted as parallel gains because they have been found only in single species belonging to phylogenetically distinctive lineages. Moreover, additional recent intron gain due to the transfer between the actin paralogues was observed in two cultured species. Based on a relaxed molecular clock timescale, we conclude that intron gains in actin took place throughout the evolution of foraminifera, with the oldest introns inserted between 550 and 500 million years ago and the youngest ones acquired less than 100 million years ago. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Debashish Bhattacharya]