Structure and expression of Euglena gracilis nuclear rbcS genes encoding the small subunits of the ribulose 1,5-bisphoshate carboxylase/oxygenase: A novel splicing process for unusual intervening sequences?

Unique features of conventional and nonconventional introns in Euglena gracilis

BACKGROUND

Nuclear introns in Euglenida have been understudied. This study aimed to investigate nuclear introns in Euglenida by identifying a large number of introns in Euglena gracilis (E. gracilis), including cis-spliced conventional and nonconventional introns, as well as trans-spliced outrons. We also examined the sequence characteristics of these introns.

RESULTS

A total of 28,337 introns and 11,921 outrons were identified. Conventional and nonconventional introns have distinct splice site features; the former harbour canonical GT/C-AG splice sites, whereas the latter are capable of forming structured motifs with their terminal sequences. We observed that short introns had a preference for canonical GT-AG introns. Notably, conventional introns and outrons in E. gracilis exhibited a distinct cytidine-rich polypyrimidine tract, in contrast to the thymidine-rich tracts observed in other organisms. Furthermore, the SL-RNAs in E. gracilis, as well as in other trans-splicing species, can form a recently discovered motif called the extended U6/5' ss duplex with the respective U6s. We also describe a novel type of alternative splicing pattern in E. gracilis. The tandem repeat sequences of introns in this protist were determined, and their contents were comparable to those in humans.

CONCLUSIONS

Our findings highlight the unique features of E. gracilis introns and provide insights into the splicing mechanism of these introns, as well as the genomics and evolution of Euglenida.

Insights into the Structure of Rubisco from Dinoflagellates-In Silico Studies

Ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) is one of the best studied enzymes. It is crucial for photosynthesis, and thus for all of biosphere’s productivity. There are four isoforms of this enzyme, differing by amino acid sequence composition and quaternary structure. However, there is still a group of organisms, dinoflagellates, single-cell eukaryotes, that are confirmed to possess Rubisco, but no successful purification of the enzyme of such origin, and hence a generation of a crystal structure was reported to date. Here, we are using in silico tools to generate the possible structure of Rubisco from a dinoflagellate representative, Symbiodinium sp. We selected two templates: Rubisco from Rhodospirillum rubrum and Rhodopseudomonas palustris. Both enzymes are the so-called form II Rubiscos, but the first is exclusively a homodimer, while the second one forms homo-hexamers. Obtained models show no differences in amino acids crucial for Rubisco activity. The variation was found at two closely located inserts in the C-terminal domain, of which one extends a helix and the other forms a loop. These inserts most probably do not play a direct role in the enzyme’s activity, but may be responsible for interaction with an unknown protein partner, possibly a regulator or a chaperone. Analysis of the possible oligomerization interface indicated that Symbiodinium sp. Rubisco most likely forms a trimer of homodimers, not just a homodimer. This hypothesis was empowered by calculation of binding energies. Additionally, we found that the protein of study is significantly richer in cysteine residues, which may be the cause for its activity loss shortly after cell lysis. Furthermore, we evaluated the influence of the loop insert, identified exclusively in the Symbiodinium sp. protein, on the functionality of the recombinantly expressed R. rubrum Rubisco. All these findings shed new light onto dinoflagellate Rubisco and may help in future obtainment of a native, active enzyme.

A New Type of Circular RNA derived from Nonconventional Introns in Nuclear Genes of Euglenids

Nuclear protein-coding genes of euglenids (Discoba, Euglenozoa, Euglenida) contain conventional (spliceosomal) and nonconventional introns. The latter have been found only in euglenozoans. A unique feature of nonconventional introns is the ability to form a stable and slightly conserved RNA secondary structure bringing together intron ends and placing adjacent exons in proximity. To date, little is known about the mechanism of their excision (e.g. whether it involves the spliceosome or not). The tubA gene of Euglena gracilis harbors three conventional and three nonconventional introns. While the conventional introns are excised as lariats, nonconventional introns are present in the cell solely as circular RNAs with full-length ends. Based on this discovery as well as on previous observations indicating that nonconventional introns are observed frequently at unique positions of genes, we suggest that this new type of intronic circRNA might play a role in intron mobility.

Comparative molecular cell biology of phototrophic euglenids and parasitic trypanosomatids sheds light on the ancestor of Euglenozoa

Parasitic trypanosomatids and phototrophic euglenids are among the most extensively studied euglenozoans. The phototrophic euglenid lineage arose relatively recently through secondary endosymbiosis between a phagotrophic euglenid and a prasinophyte green alga that evolved into the euglenid secondary chloroplast. The parasitic trypanosomatids (i.e. Trypanosoma spp. and Leishmania spp.) and the freshwater phototrophic euglenids (i.e. Euglena gracilis) are the most evolutionary distant lineages in the Euglenozoa phylogenetic tree. The molecular and cell biological traits they share can thus be considered as ancestral traits originating in the common euglenozoan ancestor. These euglenozoan ancestral traits include common mitochondrial presequence motifs, respiratory chain complexes containing various unique subunits, a unique ATP synthase structure, the absence of mitochondria-encoded transfer RNAs (tRNAs), a nucleus with a centrally positioned nucleolus, closed mitosis without dissolution of the nuclear membrane and nucleoli, a nuclear genome containing the unusual 'J' base (β-D-glucosyl-hydroxymethyluracil), processing of nucleus-encoded precursor messenger RNAs (pre-mRNAs) via spliced-leader RNA (SL-RNA) trans-splicing, post-transcriptional gene silencing by the RNA interference (RNAi) pathway and the absence of transcriptional regulation of nuclear gene expression. Mitochondrial uridine insertion/deletion RNA editing directed by guide RNAs (gRNAs) evolved in the ancestor of the kinetoplastid lineage. The evolutionary origin of other molecular features known to be present only in either kinetoplastids (i.e. polycistronic transcripts, compaction of nuclear genomes) or euglenids (i.e. monocistronic transcripts, huge genomes, many nuclear cis-spliced introns, polyproteins) is unclear.

Order of removal of conventional and nonconventional introns from nuclear transcripts of Euglena gracilis

Nuclear genes of euglenids and marine diplonemids harbor atypical, nonconventional introns which are not observed in the genomes of other eukaryotes. Nonconventional introns do not have the conserved borders characteristic for spliceosomal introns or the sequence complementary to U1 snRNA at the 5' end. They form a stable secondary structure bringing together both exon/intron junctions, nevertheless, this conformation does not resemble the form of self-splicing or tRNA introns. In the genes studied so far, frequent nonconventional introns insertions at new positions have been observed, whereas conventional introns have been either found at the conserved positions, or simply lost. In this work, we examined the order of intron removal from Euglena gracilis transcripts of the tubA and gapC genes, which contain two types of introns: nonconventional and spliceosomal. The relative order of intron excision was compared for pairs of introns belonging to different types. Furthermore, intermediate products of splicing were analyzed using the PacBio Next Generation Sequencing system. The analysis led to the main conclusion that nonconventional introns are removed in a rapid way but later than spliceosomal introns. Moreover, the observed accumulation of transcripts with conventional introns removed and nonconventional present may suggest the existence of a time gap between the two types of splicing.

The existence of conventional spliceosomal introns in genes of eukaryotic organisms is a well-known theorem. However, genes of the unicellular algae group, euglenids, contain also another type of introns, so-called nonconventional ones. They lack canonical borders, a feature most characteristic for conventional introns and form a stable secondary structure bringing together their ends. Along with the increasing popularity of whole genome studies, nonconventional introns were also disclosed in the genes of other protists, diplonemids. In this study we were particularly interested which introns–conventional or nonconventional–are removed earlier from euglenids’ pre-mRNA. To track this process we analyzed transcripts of the two Euglena gracilis genes. The relative order of intron excision was compared for pairs of introns belonging to different types. We also surveyed thousands of intermediate products of splicing using the Next-Generation Sequencing system. Summarizing the results of both experiments, we proved that spliceosomal introns are removed at an earlier stage of pre-mRNA maturation than nonconventional ones.

Intermediate introns in nuclear genes of euglenids - are they a distinct type?

Background

Nuclear genes of euglenids contain two major types of introns: conventional spliceosomal and nonconventional introns. The latter are characterized by variable non-canonical borders, RNA secondary structure that brings intron ends together, and an unknown mechanism of removal. Some researchers also distinguish intermediate introns, which combine features of both types. They form a stable RNA secondary structure and are classified into two subtypes depending on whether they contain one (intermediate/nonconventional subtype) or both (conventional/intermediate subtype) canonical spliceosomal borders. However, it has been also postulated that most introns classified as intermediate could simply be special cases of conventional or nonconventional introns.

Results

Sequences of tubB, hsp90 and gapC genes from six strains of Euglena agilis were obtained. They contain four, six, and two or three introns, respectively (the third intron in the gapC gene is unique for just one strain). Conventional introns were present at three positions: two in the tubB gene (at one position conventional/intermediate introns were also found) and one in the gapC gene. Nonconventional introns are present at ten positions: two in the tubB gene (at one position intermediate/nonconventional introns were also found), six in hsp90 (at four positions intermediate/nonconventional introns were also found), and two in the gapC gene.

Conclusions

Sequence and RNA secondary structure analyses of nonconventional introns confirmed that their most strongly conserved elements are base pairing nucleotides at positions +4, +5 and +6/ -8, −7 and −6 (in most introns CAG/CTG nucleotides were observed). It was also confirmed that the presence of the 5' GT/C end in intermediate/nonconventional introns is not the result of kinship with conventional introns, but is due to evolutionary pressure to preserve the purine at the 5' end. However, an example of a nonconventional intron with GC-AG ends was shown, suggesting the possibility of intron type conversion between nonconventional and conventional. Furthermore, an analysis of conventional introns revealed that the ability to form a stable RNA secondary structure by some introns is probably not a result of their relationship with nonconventional introns. It was also shown that acquisition of new nonconventional introns is an ongoing process and can be observed at the level of a single species. In the recently acquired intron in the gapC gene an extended direct repeats at the intron-exon junctions are present, suggesting that double-strand break repair process could be the source of new nonconventional introns.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0620-5) contains supplementary material, which is available to authorized users.

Distribution of conventional and nonconventional introns in tubulin (α and β) genes of euglenids

The nuclear genomes of euglenids contain three types of introns: conventional spliceosomal introns, nonconventional introns for which a splicing mechanism is unknown (variable noncanonical borders, RNA secondary structure bringing together intron ends), and so-called intermediate introns, which combine features of conventional and nonconventional introns. Analysis of two genes, tubA and tubB, from 20 species of euglenids reveals contrasting distribution patterns of conventional and nonconventional introns--positions of conventional introns are conserved, whereas those of the nonconventional ones are unique to individual species or small groups of closely related taxa. Moreover, in the group of phototrophic euglenids, 11 events of conventional intron loss versus 15 events of nonconventional intron gain were identified. A comparison of all nonconventional intron sequences highlighted the most conserved elements in their sequence and secondary structure. Our results led us to put forward two hypotheses. 1) The first one posits that mutational changes in intron sequence could lead to a change in their excision mechanism--intermediate introns would then be a transitional form between the conventional and nonconventional introns. 2) The second hypothesis concerns the origin of nonconventional introns--because of the presence of inverted repeats near their ends, insertion of MITE-like transposon elements is proposed as a possible source of new introns.

A possible role for short introns in the acquisition of stroma-targeting peptides in the flagellate Euglena gracilis

The chloroplasts of Euglena gracilis bounded by three membranes arose via secondary endosymbiosis of a green alga in a heterotrophic euglenozoan host. Many genes were transferred from symbiont to the host nucleus. A subset of Euglena nuclear genes of predominately symbiont, but also host, or other origin have obtained complex presequences required for chloroplast targeting. This study has revealed the presence of short introns (41–93 bp) either in the second half of presequence-encoding regions or shortly downstream of them in nine nucleus-encoded E. gracilis genes for chloroplast proteins (Eno29, GapA, PetA, PetF, PetJ, PsaF, PsbM, PsbO, and PsbW). In addition, the E. gracilis Pbgd gene contains two introns in the second half of presequence-encoding region and one at the border of presequence-mature peptide-encoding region. Ten of 12 introns present within presequence-encoding regions or shortly downstream of them identified in this study have typical eukaryotic GT/AG borders, are T-rich, 45–50 bp long, and pairwise sequence identities range from 27 to 61%. Thus single recombination events might have been mediated via these cis-spliced introns. A double crossing over between these cis-spliced introns and trans-spliced introns present in 5′-UTRs of Euglena nuclear genes is also likely to have occurred. Thus introns and exon-shuffling could have had an important role in the acquisition of chloroplast targeting signals in E. gracilis. The results are consistent with a late origin of photosynthetic euglenids.

Euglena light-harvesting complexes are encoded by multifarious polyprotein mRNAs that evolve in concert

Light-harvesting complexes (LHCs) are a superfamily of chlorophyll- and carotenoid-binding proteins that are responsible for the capture of light energy and its transfer to the photosynthetic reaction centers. Unlike those of most eukaryotes, the LHCs of Euglena gracilis are translated from large mRNAs, producing polyprotein precursors consisting of multiple concatenated LHC subunits that are separated by conserved decapeptide linkers. These precursors are posttranslationally targeted to the chloroplast and cleaved into individual proteins. We analyzed expressed sequence tags from Euglena to further characterize the structural features of the LHC polyprotein-coding genes and to examine the evolution of this multigene family. Of the 19 different LHC transcriptional units we detected, 17 encoded polyproteins composed of both tandem and nontandem repeats of LHC subunits; organizations that likely occurred through unequal crossing-over. Of the 2 nonpolyprotein-encoding LHC transcripts detected, 1 evolved from the truncation of a polyprotein-coding gene. Duplication of LHC polyprotein-coding genes was particularly important in the LHCI gene family where multiple paralogous sequences were detected. Intriguingly, several of the individual LHC-coding subunits both within and between transcriptional units appeared to be evolving in concert, suggesting that gene conversion has been a significant mechanism for LHC evolution in Euglena.

Spliceosomal introns in the deep-branching eukaryote Trichomonas vaginalis

Eukaryotes have evolved elaborate splicing mechanisms to remove introns that would otherwise destroy the protein-coding capacity of genes. Nuclear premRNA splicing requires sequence motifs in the intron and is mediated by a ribonucleoprotein complex, the spliceosome. Here we demonstrate the presence of a splicing apparatus in the protist Trichomonas vaginalis and show that RNA motifs found in yeast and metazoan introns are required for splicing. We also describe the first introns in this deep-branching lineage. The positions of these introns are often conserved in orthologous genes, indicating they were present in a common ancestor of trichomonads, yeast, and metazoa. All examined T. vaginalis introns have a highly conserved 12-nt 3' splice-site motif that encompasses the branch point and is necessary for splicing. This motif is also found in the only described intron in a gene from another deep-branching eukaryote, Giardia intestinalis. These studies demonstrate the conservation of intron splicing signals across large evolutionary distances, reveal unexpected motif conservation in deep-branching lineages that suggest a simplified mechanism of splicing in primitive unicellular eukaryotes, and support the presence of introns in the earliest eukaryote.

Role of the small subunit in ribulose-1,5-bisphosphate carboxylase/oxygenase

Ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) catalyzes the rate-limiting step of CO2 fixation in photosynthesis, but O2 competes with CO2 for substrate ribulose 1,5-bisphosphate, leading to the loss of fixed carbon. Interest in genetically engineering improvements in carboxylation catalytic efficiency and CO2/O2 specificity has focused on the chloroplast-encoded large subunit because it contains the active site. However, there is another type of subunit in the holoenzyme of plants, which, like the large subunit, is present in eight copies. The role of these nuclear-encoded small subunits in Rubisco structure and function is poorly understood. Small subunits may have originated during evolution to concentrate large-subunit active sites, but the extensive divergence of structures among prokaryotes, algae, and land plants seems to indicate that small subunits have more-specialized functions. Furthermore, plants and green algae contain families of differentially expressed small subunits, raising the possibility that these subunits may regulate the structure or function of Rubisco. Studies of interspecific hybrid enzymes have indicated that small subunits are required for maximal catalysis and, in several cases, contribute to CO2/O2 specificity. Although small-subunit genetic engineering remains difficult in land plants, directed mutagenesis of cyanobacterial and green-algal genes has identified specific structural regions that influence catalytic efficiency and CO2/O2 specificity. It is thus apparent that small subunits will need to be taken into account as strategies are developed for creating better Rubisco enzymes.

Membrane heredity and early chloroplast evolution

Membrane heredity was central to the unique symbiogenetic origin from cyanobacteria of chloroplasts in the ancestor of Plantae (green plants, red algae, glaucophytes) and to subsequent lateral transfers of plastids to form even more complex photosynthetic chimeras. Each symbiogenesis integrated disparate genomes and several radically different genetic membranes into a more complex cell. The common ancestor of Plantae evolved transit machinery for plastid protein import. In later secondary symbiogeneses, signal sequences were added to target proteins across host perialgal membranes: independently into green algal plastids (euglenoids, chlorarachneans) and red algal plastids (alveolates, chromists). Conservatism and innovation during early plastid diversification are discussed.

U1 small nuclear RNA and spliceosomal introns in Euglena gracilis

In the flagellated protozoon Euglena gracilis, characterized nuclear genes harbor atypical introns that usually are flanked by short repeats, adopt complex secondary structures in pre-mRNA, and do not obey the GT-AG rule of conventional cis-spliced introns. In the nuclear fibrillarin gene of E. gracilis, we have identified three spliceosomal-type introns that have GT-AG consensus borders. Furthermore, we have isolated a small RNA from E. gracilis and propose, on the basis of primary and secondary structure comparisons, that it is a homolog of U1 small nuclear RNA, an essential component of the cis-spliceosome in higher eukaryotes. Conserved sequences at the 5' splice sites of the fibrillarin introns can potentially base pair with Euglena U1 small nuclear RNA. Our observations demonstrate that spliceosomal GT-AG cis-splicing occurs in Euglena, in addition to the nonconventional cis-splicing and spliced leader trans-splicing previously recognized in this early diverging unicellular eukaryote.