Sequencing of random Euplotes crassus macronuclear genes supports a high frequency of +1 translational frameshifting

Eukaryotic release factor 1 from Euplotes promotes frameshifting at premature stop codons in human cells

Human physiology is highly susceptible to frameshift mutations within coding regions, and many hereditary diseases and cancers are caused by such indels. Presently, therapeutic options to counteract them are limited and, in the case of direct genome editing, risky. Here, we show that release factor 1 (eRF1) from Euplotes, an aquatic protist known for frequent +1 frameshifts in its coding regions, can enhance +1 ribosomal frameshifting at slippery heptameric sequences in human cells without an apparent requirement for an mRNA secondary structure. We further show an increase in frameshifting rate at the premature termination sequence found in the HEXA gene of Tay-Sachs disease patients, or a breast cancer cell line that harbors a tumor-driving frameshift mutation in GATA3. Although the overall increase in frameshifting would need further improvement for clinical applications, our results underscore the potential of exogenous factors, such as Eu eRF1, to increase frameshifting in human cells.

eIF5A promotes +1 programmed ribosomal frameshifting in Euplotes octocarinatus

Programmed ribosomal frameshifting (PRF) exists in all branches of life that regulate gene expression at the translational level. The single-celled eukaryote Euplotes exhibit high frequency of PRF. However, the molecular mechanism of modulating Euplotes PRF remains largely unknown. Here, we identified two novel eIF5A genes, eIF5A1 and eIF5A2, in Euplotes octocarinatus and found that the Eo-eIF5A2 gene requires a -1 PRF to produce complete protein product. Although both Eo-eIF5As showed significant structural similarity with yeast eIF5A, neither of them could functionally replace yeast eIF5A. Eo-eIF5A knockdown inhibited +1 PRF of the η-tubulin gene. Using an in vitro reconstituted translation system, we found that hypusinated Eo-eIF5A (Eo-eIF5AH) can promote +1 PRF at the canonical AAA_UAA frameshifting site of Euplotes. The results showed eIF5A is a novel trans-regulator of PRF in Euplotes and has an evolutionary conserved role in regulating +1 PRF in eukaryotes.

Nontriplet feature of genetic code in Euplotes ciliates is a result of neutral evolution

The triplet nature of the genetic code is considered a universal feature of known organisms. However, frequent stop codons at internal mRNA positions in Euplotes ciliates ultimately specify ribosomal frameshifting by one or two nucleotides depending on the context, thus posing a nontriplet feature of the genetic code of these organisms. Here, we sequenced transcriptomes of eight Euplotes species and assessed evolutionary patterns arising at frameshift sites. We show that frameshift sites are currently accumulating more rapidly by genetic drift than they are removed by weak selection. The time needed to reach the mutational equilibrium is several times longer than the age of Euplotes and is expected to occur after a several-fold increase in the frequency of frameshift sites. This suggests that Euplotes are at an early stage of the spread of frameshifting in expression of their genome. In addition, we find the net fitness burden of frameshift sites to be noncritical for the survival of Euplotes. Our results suggest that fundamental genome-wide changes such as a violation of the triplet character of genetic code can be introduced and maintained solely by neutral evolution.

The macronuclear genome of the Antarctic psychrophilic marine ciliate Euplotes focardii reveals new insights on molecular cold adaptation

The macronuclear (MAC) genomes of ciliates belonging to the genus Euplotes species are comprised of numerous small DNA molecules, nanochromosomes, each typically encoding a single gene. These genomes are responsible for all gene expression during vegetative cell growth. Here, we report the analysis of the MAC genome from the Antarctic psychrophile Euplotes focardii. Nanochromosomes containing bacterial sequences were not found, suggesting that phenomena of horizontal gene transfer did not occur recently, even though this ciliate species has a substantial associated bacterial consortium. As in other euplotid species, E. focardii MAC genes are characterized by a high frequency of translational frameshifting. Furthermore, in order to characterize differences that may be consequent to cold adaptation and defense to oxidative stress, the main constraints of the Antarctic marine microorganisms, we compared E. focardii MAC genome with those available from mesophilic Euplotes species. We focussed mainly on the comparison of tubulin, antioxidant enzymes and heat shock protein (HSP) 70 families, molecules which possess peculiar characteristic correlated with cold adaptation in E. focardii. We found that α-tubulin genes and those encoding SODs and CATs antioxidant enzymes are more numerous than in the mesophilic Euplotes species. Furthermore, the phylogenetic trees showed that these molecules are divergent in the Antarctic species. In contrast, there are fewer hsp70 genes in E. focardii compared to mesophilic Euplotes and these genes do not respond to thermal stress but only to oxidative stress. Our results suggest that molecular adaptation to cold and oxidative stress in the Antarctic environment may not only be due to particular amino acid substitutions but also due to duplication and divergence of paralogous genes.

The Sub-Chromosomic Macronuclear Pheromone Genes of the Ciliate Euplotes raikovi: Comparative Structural Analysis and Insights into the Mechanism of Expression

In Euplotes raikovi, we have determined the full-length sequences of a family of macronuclear genes that are the transcriptionally active versions of codominant alleles inherited at the mating-type (mat) locus of the micronuclear genome, and encode cell type-distinctive signaling pheromones. These genes include a 225-231-bp coding region flanked by a conserved 544-bp 5'-leader region and a more variable 3'-trailer region. Two transcription initiation start sites and two polyadenylation sites associated with nonconventional signals cooperate with a splicing phenomenon of a 326-bp intron residing in the 5'-leader region in the generation of multiple transcripts from the same gene. In two of them, the synthesis of functional products depends on the reassignment to a sense codon, or readthrough of a strictly conserved leaky UAG stop codon. That this reassignment may take place is suggested by the position this codon occupies in the transcripts, close to the transcript extremity and far from the poly(A) tail. In such a case, one product is a 69-amino acid protein in search of function and the second product is a 126-amino acid protein that represents a membrane-bound pheromone isoform candidate to function as a cell type-specific binding site (receptor) of the soluble pheromones.

Expression of a MORN repeat protein from Euplotes octocarinatus requires a +1 programmed ribosomal frameshifting

Analysis of transcriptome revealed that a membrane occupation and recognition nexus (MORN) repeat protein-encoding gene of Euplotes octocarinatus (Eo-morn-9-31) was a candidate for programmed +1 ribosomal frameshifting (+1 PRF). In this study, a dual-luciferase assay was performed to detect its expression. The result showed that the MORN repeat protein (Eo-MORN-9-31) could be produced by the +1 PRF event during the process of translation in yeast and the frameshifting efficiency was about 4-5%. We further confirmed its reality by western blot and mass spectrometry. This study provided experimental evidence indicating that the expression of the Eo-MORN-9-31 of E. octocarinatus required the +1 PRF.

Large-scale mass spectrometry-based analysis of Euplotes octocarinatus supports the high frequency of +1 programmed ribosomal frameshift

Programmed ribosomal frameshifting (PRF) is commonly used to express many viral and some cellular genes. We conducted a genome-wide investigation of +1 PRF in ciliate Euplotes octocarinatus through genome and transcriptome sequencing and our results demonstrated that approximately 11.4% of genes require +1 PRF to produce complete gene products. While nucleic acid-based evidence for candidate genes with +1 PRF is strong, only very limited information is available at protein levels to date. In this study, E. octocarinatus was subjected to large-scale mass spectrometry-based analysis to verify the high frequency of +1 PRF and 226 +1 PRF gene products were identified. Based on the amino acid sequences of the peptides spanning the frameshift sites, typical frameshift motif AAA-UAR for +1 PRF in Euplotes was identified. Our data in this study provide very useful insight into the understanding of the molecular mechanism of +1 PRF.

Ribosomal frameshifting and transcriptional slippage: From genetic steganography and cryptography to adventitious use

Genetic decoding is not ‘frozen’ as was earlier thought, but dynamic. One facet of this is frameshifting that often results in synthesis of a C-terminal region encoded by a new frame. Ribosomal frameshifting is utilized for the synthesis of additional products, for regulatory purposes and for translational ‘correction’ of problem or ‘savior’ indels. Utilization for synthesis of additional products occurs prominently in the decoding of mobile chromosomal element and viral genomes. One class of regulatory frameshifting of stable chromosomal genes governs cellular polyamine levels from yeasts to humans. In many cases of productively utilized frameshifting, the proportion of ribosomes that frameshift at a shift-prone site is enhanced by specific nascent peptide or mRNA context features. Such mRNA signals, which can be 5′ or 3′ of the shift site or both, can act by pairing with ribosomal RNA or as stem loops or pseudoknots even with one component being 4 kb 3′ from the shift site. Transcriptional realignment at slippage-prone sequences also generates productively utilized products encoded trans-frame with respect to the genomic sequence. This too can be enhanced by nucleic acid structure. Together with dynamic codon redefinition, frameshifting is one of the forms of recoding that enriches gene expression.

High frequency of +1 programmed ribosomal frameshifting in Euplotes octocarinatus

Programmed -1 ribosomal frameshifting (-1 PRF) has been identified as a mechanism to regulate the expression of many viral genes and some cellular genes. The slippery site of -1 PRF has been well characterized, whereas the +1 PRF signal and the mechanism involved in +1 PRF remain poorly understood. Previous study confirmed that +1 PRF is required for the synthesis of protein products in several genes of ciliates from the genus Euplotes. To accurately assess the frequency of genes requiring frameshift in Euplotes, the macronuclear genome and transcriptome of Euplotes octocarinatus were analyzed in this study. A total of 3,700 +1 PRF candidate genes were identified from 32,353 transcripts, and the gene products of these putative +1 PRFs were mainly identified as protein kinases. Furthermore, we reported a putative suppressor tRNA of UAA which may provide new insights into the mechanism of +1 PRF in euplotids. For the first time, our transcriptome-wide survey of +1 PRF in E. octocarinatus provided a dataset which serves as a valuable resource for the future understanding of the mechanism underlying +1 PRF.

Characterization and expression of the gene encoding En-MAPK1, an intestinal cell kinase (ICK)-like kinase activated by the autocrine pheromone-signaling loop in the Polar Ciliate, Euplotes nobilii

In the protozoan ciliate Euplotes, a transduction pathway resulting in a mitogenic cell growth response is activated by autocrine receptor binding of cell type-specific, water-borne signaling protein pheromones. In Euplotes raikovi, a marine species of temperate waters, this transduction pathway was previously shown to involve the phosphorylation of a nuclear protein kinase structurally similar to the intestinal-cell and male germ cell-associated kinases described in mammals. In E. nobilii, which is phylogenetically closely related to E. raikovi but inhabits Antarctic and Arctic waters, we have now characterized a gene encoding a structurally homologous kinase. The expression of this gene requires +1 translational frameshifting and a process of intron splicing for the production of the active protein, designated En-MAPK1, which contains amino acid substitutions of potential significance for cold-adaptation.

Genomes on the edge: programmed genome instability in ciliates

Ciliates are an ancient and diverse group of microbial eukaryotes that have emerged as powerful models for RNA-mediated epigenetic inheritance. They possess extensive sets of both tiny and long noncoding RNAs that, together with a suite of proteins that includes transposases, orchestrate a broad cascade of genome rearrangements during somatic nuclear development. This Review emphasizes three important themes: the remarkable role of RNA in shaping genome structure, recent discoveries that unify many deeply diverged ciliate genetic systems, and a surprising evolutionary "sign change" in the role of small RNAs between major species groups.

Methionine sulfoxide reduction in ciliates: characterization of the ready-to-use methionine sulfoxide-R-reductase genes in Euplotes

Genes encoding the enzyme methionine sulfoxide reductase type B, specific to the reduction of the oxidized methionine-R form, were characterized from the expressed (macronuclear) genome of two ecologically separate marine species of Euplotes, i.e. temperate water E. raikovi and polar water E. nobilii. Both species were found to contain a single msrB gene with a very simple structural organization encoding a protein of 127 (E. raikovi) or 126 (E. nobilii) amino acid residues that belongs to the group of zinc-containing enzymes. Both msrB genes are constitutively expressed, suggesting that the MsrB enzyme plays an essential role in repairing oxidative damages that appear to be primarily caused by physiological cell aging in E. raikovi and by interactions with an O(2) saturated environment in E. nobilii.

Extensive frameshift at all AGG and CCC codons in the mitochondrial cytochrome c oxidase subunit 1 gene of Perkinsus marinus (Alveolata; Dinoflagellata)

Diverse mitochondrial (mt) genetic systems have evolved independently of the more uniform nuclear system and often employ modified genetic codes. The organization and genetic system of dinoflagellate mt genomes are particularly unusual and remain an evolutionary enigma. We determined the sequence of full-length cytochrome c oxidase subunit 1 (cox1) mRNA of the earliest diverging dinoflagellate Perkinsus and show that this gene resides in the mt genome. Apparently, this mRNA is not translated in a single reading frame with standard codon usage. Our examination of the nucleotide sequence and three-frame translation of the mRNA suggest that the reading frame must be shifted 10 times, at every AGG and CCC codon, to yield a consensus COX1 protein. We suggest two possible mechanisms for these translational frameshifts: a ribosomal frameshift in which stalled ribosomes skip the first bases of these codons or specialized tRNAs recognizing non-triplet codons, AGGY and CCCCU. Regardless of the mechanism, active and efficient machinery would be required to tolerate the frameshifts predicted in Perkinsus mitochondria. To our knowledge, this is the first evidence of translational frameshifts in protist mitochondria and, by far, is the most extensive case in mitochondria.

Sequences Promoting Recoding Are Singular Genomic Elements

The distribution of sequences which induce non-standard decoding, especially of shift-prone sequences, is very unusual. On one hand, since they can disrupt standard genetic readout, they are avoided within the coding regions of most genes. On the other hand, they play important regulatory roles for the expression of those genes where they do occur. As a result, they are preserved among homologs and exhibit deep phylogenetic conservation. The combination of these two constraints results in a characteristic distribution of recoding sequences across genomes: they are highly conserved at specific locations while they are very rare in other locations. We term such sequences singular genomic elements to signify their rare occurrence and biological importance.

Codon size reduction as the origin of the triplet genetic code

The genetic code appears to be optimized in its robustness to missense errors and frameshift errors. In addition, the genetic code is near-optimal in terms of its ability to carry information in addition to the sequences of encoded proteins. As evolution has no foresight, optimality of the modern genetic code suggests that it evolved from less optimal code variants. The length of codons in the genetic code is also optimal, as three is the minimal nucleotide combination that can encode the twenty standard amino acids. The apparent impossibility of transitions between codon sizes in a discontinuous manner during evolution has resulted in an unbending view that the genetic code was always triplet. Yet, recent experimental evidence on quadruplet decoding, as well as the discovery of organisms with ambiguous and dual decoding, suggest that the possibility of the evolution of triplet decoding from living systems with non-triplet decoding merits reconsideration and further exploration. To explore this possibility we designed a mathematical model of the evolution of primitive digital coding systems which can decode nucleotide sequences into protein sequences. These coding systems can evolve their nucleotide sequences via genetic events of Darwinian evolution, such as point-mutations. The replication rates of such coding systems depend on the accuracy of the generated protein sequences. Computer simulations based on our model show that decoding systems with codons of length greater than three spontaneously evolve into predominantly triplet decoding systems. Our findings suggest a plausible scenario for the evolution of the triplet genetic code in a continuous manner. This scenario suggests an explanation of how protein synthesis could be accomplished by means of long RNA-RNA interactions prior to the emergence of the complex decoding machinery, such as the ribosome, that is required for stabilization and discrimination of otherwise weak triplet codon-anticodon interactions.

Connection between stop codon reassignment and frequent use of shifty stop frameshifting

Ciliated protozoa of the genus Euplotes have undergone genetic code reassignment, redefining the termination codon UGA to encode cysteine. In addition, Euplotes spp. genes very frequently employ shifty stop frameshifting. Both of these phenomena involve noncanonical events at a termination codon, suggesting they might have a common cause. We recently demonstrated that Euplotes octocarinatus peptide release factor eRF1 ignores UGA termination codons while continuing to recognize UAA and UAG. Here we show that both the Tetrahymena thermophila and E. octocarinatus eRF1 factors allow efficient frameshifting at all three termination codons, suggesting that UGA redefinition also impaired UAA/UAG recognition. Mutations of the Euplotes factor restoring a phylogenetically conserved motif in eRF1 (TASNIKS) reduced programmed frameshifting at all three termination codons. Mutation of another conserved residue, Cys124, strongly reduces frameshifting at UGA while actually increasing frameshifting at UAA/UAG. We will discuss these results in light of recent biochemical characterization of these mutations.