A mechanism for a single nucleotide intron shift

An ascomycete H4 variant with an unknown function

Histone variants leading to altered nucleosome structure, dynamics and DNA accessibility occur frequently, albeit rarely for H4. We carried out a comprehensive in silico scrutiny of fungal genomes, which revealed the presence of a novel H4 variant (H4E) in the ascomycetes, throughout the Pezizomycotina, in basal species of the Taphrinomycotina and also in the Glomeromycota. The coding cognate genes show a specific intron/exon organization, different from H4 canonical genes. H4Es diverge from canonical H4s mainly in the N- and C-terminal extensions, showing marked differences in the distribution and number of Lys and Arg residues, which may result in novel post-translational modifications. In Aspergillus nidulans (Pezizomycotina, Eurotiomycetes) the H4E variant protein level is low in mycelia. However, the encoding gene is well expressed at 37°C under nitrogen starvation. H4E localizes to the nucleus and interacts with H3, but its absence or overexpression does not result in any detectable phenotype. Deletion of only one of the of the two canonical H4 genes results in a strikingly impaired growth phenotype, which indicates that H4E cannot replace this canonical histone. Thus, an H4 variant is present throughout a whole subphylum of the ascomycetes, but with hitherto no experimentally detectable function.

Bi-alignments with affine gaps costs

Background

Commonly, sequence and structure elements are assumed to evolve congruently, such that homologous sequence positions correspond to homologous structural features. Assuming congruent evolution, alignments based on sequence and structure similarity can therefore optimize both similarities at the same time in a single alignment. To model incongruent evolution, where sequence and structural features diverge positionally, we recently introduced bi-alignments. This generalization of sequence and structure-based alignments is best understood as alignments of two distinct pairwise alignments of the same entities: one modeling sequence similarity, the other structural similarity.

Results

Optimal bi-alignments with affine gap costs (or affine shift cost) for two constituent alignments can be computed exactly in quartic space and time. Even bi-alignments with affine shift and gap cost, as well as bi-alignment with sub-additive gap cost are optimized efficiently. Affine gap-cost bi-alignment of large proteins (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 930$$\end{document}∼930 aa) can be computed.

Conclusion

Affine cost bi-alignments are of practical interest to study shifts of protein sequences and protein structures relative to each other.

Availability

The affine cost bi-alignment algorithm has been implemented in Python 3 and Cython. It is available as free software from https://github.com/s-will/BiAlign/releases/tag/v0.3 and as bioconda package bialign.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13015-022-00219-7.

Unique and Repeated Stwintrons (Spliceosomal Twin Introns) in the Hypoxylaceae

Introns are usually non-coding sequences interrupting open reading frames in pre-mRNAs [D1,2]. Stwintrons are nested spliceosomal introns, where an internal intron splits a second donor sequence into two consecutive splicing reactions leading to mature mRNA. In Hypoxylon sp. CO27-5, 36 highly sequence-similar [D1,2] stwintrons are extant (sister stwintrons). An additional 81 [D1,2] sequence-unrelated stwintrons are described here. Most of them are located at conserved gene positions rooted deep in the Hypoxylaceae. Absence of exonic sequence bias at the exon–stwintron junctions and a very similar phase distribution were noted for both groups. The presence of an underlying sequence symmetry in all 117 stwintrons was striking. This symmetry, more pronounced near the termini of most of the full-length sister stwintrons, may lead to a secondary structure that brings into close proximity the most distal splice sites, the donor of the internal and the acceptor of the external intron. The Hypoxylon stwintrons were overwhelmingly excised by consecutive splicing reactions precisely removing the whole intervening sequence, whereas one excision involving the distal splice sites led to a frameshift. Alternative (mis)splicing took place for both sister and uniquely occurring stwintrons. The extraordinary symmetry of the sister stwintrons thus seems dispensable for the infrequent, direct utilisation of the distal splice sites.

Internally Symmetrical Stwintrons and Related Canonical Introns in Hypoxylaceae Species

Spliceosomal introns are pervasive in eukaryotes. Intron gains and losses have occurred throughout evolution, but the origin of new introns is unclear. Stwintrons are complex intervening sequences where one of the sequence elements (5′-donor, lariat branch point element or 3′-acceptor) necessary for excision of a U2 intron (external intron) is itself interrupted by a second (internal) U2 intron. In Hypoxylaceae, a family of endophytic fungi, we uncovered scores of donor-disrupted stwintrons with striking sequence similarity among themselves and also with canonical introns. Intron–exon structure comparisons suggest that these stwintrons have proliferated within diverging taxa but also give rise to proliferating canonical introns in some genomes. The proliferated (stw)introns have integrated seamlessly at novel gene positions. The recently proliferated (stw)introns appear to originate from a conserved ancestral stwintron characterised by terminal inverted repeats (45–55 nucleotides), a highly symmetrical structure that may allow the formation of a double-stranded intron RNA molecule. No short tandem duplications flank the putatively inserted intervening sequences, which excludes a DNA transposition-based mechanism of proliferation. It is tempting to suggest that this highly symmetrical structure may have a role in intron proliferation by (an)other mechanism(s).

Spliceosomal Introns: Features, Functions, and Evolution

Spliceosomal introns, which have been found in most eukaryotic genes, are non-coding sequences excised from pre-mRNAs by a special complex called spliceosome during mRNA splicing. Introns occur in both protein- and RNA-coding genes and can be found in coding and untranslated gene regions. Because intron sequences vary greatly due to a high rate of polymorphism, the functions of intron had been for a long time associated only with alternative splicing, while intron evolution had been viewed not as an evolution of an individual genomic element, but rather considered within a framework of the evolution of the gene intron-exon structure. Here, we review the theories of intron origin, evolutionary events in the exon-intron structure, such as intron gain, loss, and sliding, intron functions known to date, and mechanisms by which changes in the intron features (length and phase) can affect the regulation of gene-mediated processes.

Is there any intron sliding in mammals?

BACKGROUND

Eukaryotic protein-coding genes consist of exons and introns. Exon-intron borders are conserved between species and thus their changes might be observed only on quite long evolutionary distances. One of the rarest types of change, in which intron relocates over a short distance, is called "intron sliding", but the reality of this event has been debated for a long time. The main idea of a search for intron sliding is to use the most accurate genome annotation and genome sequence, as well as high-quality transcriptome data. We applied them in a search for sliding introns in mammals in order to widen knowledge about the presence or absence of such phenomena in this group.

RESULTS

We didn't find any significant evidence of intron sliding in the primate group (human, chimpanzee, rhesus macaque, crab-eating macaque, green monkey, marmoset). Only one possible intron sliding event supported by a set of high quality transcriptomes was observed between EIF1AX human and sheep gene orthologs. Also, we checked a list of previously observed intron sliding events in mammals and showed that most likely they are artifacts of genome annotations and are not shown in subsequent annotation versions as well as are not supported by transcriptomic data.

CONCLUSIONS

We assume that intron sliding is indeed a very rare evolutionary event if it exists at all. Every case of intron sliding needs a lot of supportive data for detection and confirmation.

Phylogenetic Analyses of Glycosyl Hydrolase Family 6 Genes in Tunicates: Possible Horizontal Transfer

Horizontal gene transfer (HGT) is the movement of genetic material between different species. Although HGT is less frequent in eukaryotes than in bacteria, several instances of HGT have apparently shaped animal evolution. One well-known example is the tunicate cellulose synthase gene, CesA, in which a gene, probably transferred from bacteria, greatly impacted tunicate evolution. A Glycosyl Hydrolase Family 6 (GH6) hydrolase-like domain exists at the C-terminus of tunicate CesA, but not in cellulose synthases of other organisms. The recent discovery of another GH6 hydrolase-like gene (GH6-1) in tunicate genomes further raises the question of how tunicates acquired GH6. To examine the probable origin of these genes, we analyzed the phylogenetic relationship of GH6 proteins in tunicates and other organisms. Our analyses show that tunicate GH6s, the GH6-1 gene, and the GH6 part of the CesA gene, form two independent, monophyletic gene groups. We also compared their sequence signatures and exon splice sites. All tunicate species examined have shared splice sites in GH6-containing genes, implying ancient intron acquisitions. It is likely that the tunicate CesA and GH6-1 genes existed in the common ancestor of all extant tunicates.

A functional twintron, 'zombie' twintrons and a hypermobile group II intron invading itself in plant mitochondria

The occurrence of group II introns in plant mitochondrial genomes is strikingly different between the six major land plant clades, contrasting their highly conserved counterparts in chloroplast DNA. Their present distribution likely reflects numerous ancient intron gains and losses during early plant evolution before the emergence of seed plants. As a novelty for plant organelles, we here report on five cases of twintrons, introns-within-introns, in the mitogenomes of lycophytes and hornworts. An internal group II intron interrupts an intron-borne maturase of an atp9 intron in Lycopodiaceae, whose splicing precedes splicing of the external intron. An invasive, hypermobile group II intron in cox1, has conquered nine further locations including a previously overlooked sdh3 intron and, most surprisingly, also itself. In those cases, splicing of the external introns does not depend on splicing of the internal introns. Similar cases are identified in the mtDNAs of hornworts. Although disrupting a group I intron-encoded protein in one case, we could not detect splicing of the internal group II intron in this ‘mixed’ group I/II twintron. We suggest the name ‘zombie’ twintrons (half-dead, half-alive) for such cases where splicing of external introns does not depend any more on prior splicing of fossilized internal introns.

Complex intron generation in the yeast genus Lipomyces

In primary transcripts of eukaryotic nuclear genes, coding sequences are often interrupted by U2-type introns. Such intervening sequences can constitute complex introns excised by consecutive splicing reactions. The origin of spliceosomal introns is a vexing problem. Sequence variation existent across fungal taxa provides means to study their structure and evolution. In one class of complex introns called [D] stwintrons, an (internal) U2 intron is nested within the 5'-donor element of another (external) U2 intron. In the gene for a reticulon-like protein in species of the ascomycete yeast genus Lipomyces, the most 5' terminal intron position is occupied by one of three complex intervening sequences consistent of differently nested U2 intron units, as demonstrated in L. lipofer, L. suomiensis, and L. starkeyi. In L. starkeyi, the donor elements of the constituent introns are abutting and the complex intervening sequence can be excised alternatively either with one standard splicing reaction or, as a [D] stwintron, by two consecutive reactions. Our work suggests how [D] stwintrons could emerge by the appearance of new functional splice sites within an extant intron. The stepwise stwintronisation mechanism may involve duplication of the functional intron donor element of the ancestor intron.

Comprehensive genomic analyses with 115 plastomes from algae to seed plants: structure, gene contents, GC contents, and introns

BACKGROUND

Chloroplasts are a common character in plants. The chloroplasts in each plant lineage have shaped their own genomes, plastomes, by structural changes and transferring many genes to nuclear genomes during plant evolution. Some plastid genes have introns that are mostly group II introns.

OBJECTIVE

This study aimed to get genomic and evolutionary insights on the plastomes from green algae to flowering plants.

METHODS

Plastomes of 115 species from green algae, bryophytes, pteridophytes (spore bearing vascular plants), gymnosperms, and angiosperms were mined from NCBI organelle genome database. Plastome structure, gene contents and GC contents were analyzed by the in-house developed Phyton code. Intronic features including presence/absence, length, intron phases were analyzed by manually in the annotated information in NCBI.

RESULTS

The canonical quadripartite structures were retained in most plastomes except of a few plastomes that had lost an invert repeat (IR). Expansion or reduction or deletion of IRs resulted in the length variation of the plastomes. The number of protein coding genes ranged from 40 to 92 with an average 79.43 ± 5.84 per plastome and gene losses were apparent in specific lineages. The number of trn genes ranged from 13 to 33 with an average 21.19 ± 2.42 per plastome. Ribosomal RNA genes, rrn, were located in the IRs so that they were present in a duplicate except of the species that had lost one of the IR. GC contents were variable from 24.9 to 51.0% with an average 38.21 ± 3.27%, indicating bias to high AT contents. Plastid introns were present in 18 protein coding genes, six trn genes, and one rrn gene. Intron losses occurred among the orthologous genes in different plant lineages. The plastid introns were long compared with the nuclear introns, which might be related with the spliceosome nuclear introns and self-splicing group II plastid introns. The trnK-UUU intron contained the maturase encoding matK gene except in the chlorophyte algae and monilophyte ferns in which the trnK-UUU was lost, but matK retained. There were many annotation artefacts in the intron positions in the NCBI database. In the analysis of intron phases, phase 0 introns were more frequent than those of phase 2 and 3 introns. Phase polymorphism was observed in the introns of clpP which was derived from nucleotide insertion. Plastid trn introns were long compared to the archaeal or eukaryotic nuclear tRNA introns. Of the six plastid trn introns, one was at the D loop and other five were at the anticodon loop. The insertion sites were conserved among the trn genes in archaea, eukaryotic nuclear and plastid tRNA genes.

CONCLUSIONS

Current study refurbrished the previous findings of structural variations, gene contents, and GC contents of the chloroplast genomes from green algae to flowering plants. The study also included some noble findings and discussions on the plastome introns including their length variations and phase variation. We also presented and corrected some false annotations on the introns in protein coding and tRNA genes in the genome database, which might be confirmed by the chloroplast transcriptome analysis in the future.

A spliceosomal twin intron (stwintron) participates in both exon skipping and evolutionary exon loss

Spliceosomal twin introns (stwintrons) are introns where any of the three consensus sequences involved in splicing is interrupted by another intron (internal intron). In Aspergillus nidulans, a donor-disrupted stwintron (intron-1) is extant in the transcript encoding a reticulon-like protein. The orthologous transcript of Aspergillus niger can be alternatively spliced; the exon downstream the stwintron could be skipped by excising a sequence that comprises this stwintron, the neighbouring intron-2, and the exon bounded by these. This process involves the use of alternative 3' splice sites for the internal intron, the resulting alternative intervening sequence being a longer 3'-extended stwintron. In 29 species of Onygenales, a multi-step splicing process occurs in the orthologous transcript, in which a complex intervening sequence including the stwintron and neigbouring intron-2, generates by three splicing reactions a "second order intron" which must then be excised with a fourth splicing event. The gene model in two species can be envisaged as one canonical intron (intron-1) evolved from this complex intervening sequence of nested canonical introns found elsewhere in Onygenales. Postulated splicing intermediates were experimentally verified in one or more species. This work illustrates a role of stwintrons in both alternative splicing and the evolution of intron structure.

Emergence and loss of spliceosomal twin introns

BACKGROUND

In the primary transcript of nuclear genes, coding sequences-exons-usually alternate with non-coding sequences-introns. In the evolution of spliceosomal intron-exon structure, extant intron positions can be abandoned and new intron positions can be occupied. Spliceosomal twin introns ("stwintrons") are unconventional intervening sequences where a standard "internal" intron interrupts a canonical splicing motif of a second, "external" intron. The availability of genome sequences of more than a thousand species of fungi provides a unique opportunity to study spliceosomal intron evolution throughout a whole kingdom by means of molecular phylogenetics.

RESULTS

A new stwintron was encountered in Aspergillus nidulans and Aspergillus niger. It is present across three classes of Leotiomyceta in the transcript of a well-conserved gene encoding a putative lipase (lipS). It occupies the same position as a standard intron in the orthologue gene in species of the early divergent classes of the Pezizomycetes and the Orbiliomycetes, suggesting that an internal intron has appeared within a pre-extant intron. On the other hand, the stwintron has been lost from certain taxa in Leotiomycetes and Eurotiomycetes at several occasions, most likely by a mechanism involving reverse transcription and homologous recombination. Another ancient stwintron present across whole Pezizomycotina orders-in the transcript of the bifunctional biotin biosynthesis gene bioDA-occurs at the same position as a standard intron in many species of non-Dikarya. Nevertheless, also the bioDA stwintron has disappeared from certain lineages within the taxa where it occurs, i.e., Sordariomycetes and Botryosphaeriales. Intriguingly, only the internal intron was lost from the Sordariomycetes bioDA stwintron at all but one occasion, leaving a standard intron in the same position, while where the putative lipase stwintron was lost, no intronic sequences remain.

CONCLUSIONS

Molecular phylogeny of the peptide product was used to monitor the existence and fate of a stwintron in the transcripts of two neatly defined fungal genes, encoding well conserved proteins. Both defining events-stwintron emergence and loss-can be explained with extant models for intron insertion and loss. We thus demonstrate that stwintrons can serve as model systems to study spliceosomal intron evolution.