Deviation from the universal code shown by the gene for surface protein 51A in Paramecium

Cracking the Code: Reprogramming the Genetic Script in Prokaryotes and Eukaryotes to Harness the Power of Noncanonical Amino Acids

Over 500 natural and synthetic amino acids have been genetically encoded in the last two decades. Incorporating these noncanonical amino acids into proteins enables many powerful applications, ranging from basic research to biotechnology, materials science, and medicine. However, major challenges remain to unleash the full potential of genetic code expansion across disciplines. Here, we provide an overview of diverse genetic code expansion methodologies and systems and their final applications in prokaryotes and eukaryotes, represented by Escherichia coli and mammalian cells as the main workhorse model systems. We highlight the power of how new technologies can be first established in simple and then transferred to more complex systems. For example, whole-genome engineering provides an excellent platform in bacteria for enabling transcript-specific genetic code expansion without off-targets in the transcriptome. In contrast, the complexity of a eukaryotic cell poses challenges that require entirely new approaches, such as striving toward establishing novel base pairs or generating orthogonally translating organelles within living cells. We connect the milestones in expanding the genetic code of living cells for encoding novel chemical functionalities to the most recent scientific discoveries, from optimizing the physicochemical properties of noncanonical amino acids to the technological advancements for their in vivo incorporation. This journey offers a glimpse into the promising developments in the years to come.

Noncanonical usage of stop codons in ciliates expands proteins with structurally flexible Q-rich motifs

Serine(S)/threonine(T)-glutamine(Q) cluster domains (SCDs), polyglutamine (polyQ) tracts and polyglutamine/asparagine (polyQ/N) tracts are Q-rich motifs found in many proteins. SCDs often are intrinsically disordered regions that mediate protein phosphorylation and protein-protein interactions. PolyQ and polyQ/N tracts are structurally flexible sequences that trigger protein aggregation. We report that due to their high percentages of STQ or STQN amino acid content, four SCDs and three prion-causing Q/N-rich motifs of yeast proteins possess autonomous protein expression-enhancing activities. Since these Q-rich motifs can endow proteins with structural and functional plasticity, we suggest that they represent useful toolkits for evolutionary novelty. Comparative Gene Ontology (GO) analyses of the near-complete proteomes of 26 representative model eukaryotes reveal that Q-rich motifs prevail in proteins involved in specialized biological processes, including Saccharomyces cerevisiae RNA-mediated transposition and pseudohyphal growth, Candida albicans filamentous growth, ciliate peptidyl-glutamic acid modification and microtubule-based movement, Tetrahymena thermophila xylan catabolism and meiosis, Dictyostelium discoideum development and sexual cycles, Plasmodium falciparum infection, and the nervous systems of Drosophila melanogaster, Mus musculus and Homo sapiens. We also show that Q-rich-motif proteins are expanded massively in 10 ciliates with reassigned TAAQ and TAGQ codons. Notably, the usage frequency of CAGQ is much lower in ciliates with reassigned TAAQ and TAGQ codons than in organisms with expanded and unstable Q runs (e.g. D. melanogaster and H. sapiens), indicating that the use of noncanonical stop codons in ciliates may have coevolved with codon usage biases to avoid triplet repeat disorders mediated by CAG/GTC replication slippage.

Flatworm Transcriptomes Reveal Widespread Parasitism by Histophagous Ciliates

Unicellular ciliates like Tetrahymena are best known as free-living bacteriovores, but many species are facultative or obligate parasites. These "histophages" feed on the tissues of hosts ranging from planarian flatworms to commercially important fish and the larvae of imperiled freshwater mussels. Here, we developed a novel bioinformatics pipeline incorporating the nonstandard ciliate genetic code and used it to search for Ciliophora sequences in 34 publicly available Platyhelminthes EST libraries. From 2,615,036 screened ESTs, we identified nearly 6,000 high-confidence ciliate transcripts, supporting parasitism of seven additional flatworm species. We also cultured and identified Tetrahymena from nine terrestrial and freshwater planarians, including invasive earthworm predators from the genus Bipalium and the widely studied regeneration models Dugesia japonica and Schmidtea mediterranea. A co-phylogenetic reconstruction provides strong evidence for the coevolution of histophagous Ciliophora with their Platyhelminthes hosts. We further report the antiprotozoal aminoglycoside paromomycin expels Tetrahymena from S. mediterranea, providing new opportunities to investigate the effects of this relationship on planarian biology. Together, our findings raise the possibility that invasive flatworms constitute a novel dispersal mechanism for Tetrahymena parasites and position the Platyhelminthes as an ideal model phylum for studying the ecology and evolution of histophagous ciliates.

FLATWORM TRANSCRIPTOMES REVEAL WIDESPREAD PARASITISM BY HISTOPHAGOUS CILIATES

Unicellular ciliates like Tetrahymena are best known as free-living bacteriovores, but many species are facultative or obligate parasites. These 'histophages' feed on the tissues of hosts ranging from planarian flatworms to commercially important fish and the larvae of imperiled freshwater mussels. Here, we developed a novel bioinformatics pipeline incorporating the nonstandard ciliate genetic code and used it to search for Ciliophora sequences in 34 publicly available Platyhelminthes EST libraries. From 2,615,036 screened ESTs, we identified nearly 6,000 high-confidence ciliate transcripts, supporting parasitism of seven additional flatworm species. We also cultured and identified Tetrahymena from nine terrestrial and freshwater planarians, including invasive earthworm predators from the genus Bipalium and the widely studied regeneration models Dugesia japonica and Schmidtea mediterranea. A cophylogenetic reconstruction provides strong evidence for coevolution of histophagous Ciliophora with their Platyhelminthes hosts. We further report the antiprotozoal aminoglycoside paromomycin expels Tetrahymena from S. mediterranea, providing new opportunities to investigate the effects of this relationship on planarian biology. Together, our findings raise the possibility that invasive flatworms constitute a novel dispersal mechanism for Tetrahymena parasites and position the Platyhelminthes as an ideal model phylum for studying the ecology and evolution of histophagous ciliates.

RNA interference reveals the escape response mechanism of Paramecium to mechanical stimulation

In Paramecium, a mechanical stimulus applied to the posterior portion of the cell causes a transient increase in membrane permeability to potassium ions, transiently rendering the membrane in a hyperpolarized state. Hyperpolarization causes a transient increase in Cyclic adenosine monophosphate (cAMP) concentration in the cilia, resulting in a transient fast-forward swimming of the cell. Schultz and coworkers (1992) reported that a unique adenylate cyclase (AC)-coupled potassium channel is involved in the reaction underlying this response, which is known as the "escape response." However, the AC responsible for this reaction remains to be identified. Moreover, the molecular linkage between mechanoreception and AC activation has not been elucidated adequately. Currently, we can perform an efficient and simple gene-knockdown technique in Paramecium using RNA interference (RNAi). Paramecium is one of the several model organisms for which whole-genome sequences have been elucidated. The RNAi technique can be applied to whole genome sequences derived from the Paramecium database (ParameciumDB) to investigate the types of proteins that elicit specific biological responses and compare them with those of other model organisms. In this review, we describe the applications of the RNAi technique in elucidating the molecular mechanism underlying the escape response and identifying the AC involved in this reaction. The findings of this study highlight the advantages of the RNAi technique and ParameciumDB.

Species-Specific Duplication of Surface Antigen Genes in Paramecium

Paramecium is a free-living ciliate that undergoes antigenic variation and still the functions of these variable surface antigen coats in this non-pathogenic ciliate remain elusive. Only a few surface antigen genes have been described, mainly in the two model species P. tetraurelia strain 51 and P. primaurelia strain 156. Given the lack of suitable sequence data to allow for phylogenetics and deeper sequence comparisons, we screened the genomes of six different Paramecium species for serotype genes and isolated 548 candidates. Our approach identified the subfamilies of the isogenes of individual serotypes that were mostly represented by intrachromosomal gene duplicates. These showed different duplication levels, and chromosome synteny suggested rather young duplication events after the emergence of the P. aurelia species complex, indicating a rapid evolution of surface antigen genes. We were able to identify the different subfamilies of the surface antigen genes with internal tandem repeats, which showed consensus motifs across species. The individual isogene families showed additional consensus motifs, indicating that the selection pressure holds individual amino acids constant in these repeats. This may be a hint of the receptor function of these antigens rather than a presentation of random epitopes, generating the variability of these surface molecules.

Improving the Efficiency and Orthogonality of Genetic Code Expansion

The site-specific incorporation of the noncanonical amino acid (ncAA) into proteins via genetic code expansion (GCE) has enabled the development of new and powerful ways to learn, regulate, and evolve biological functions in vivo. However, cellular biosynthesis of ncAA-containing proteins with high efficiency and fidelity is a formidable challenge. In this review, we summarize up-to-date progress towards improving the efficiency and orthogonality of GCE and enhancing intracellular compatibility of introduced translation machinery in the living cells by creation and optimization of orthogonal translation components, constructing genomically recoded organism (GRO), utilization of unnatural base pairs (UBP) and quadruplet codons (four-base codons), and spatial separation of orthogonal translation.

Widespread stop-codon recoding in bacteriophages may regulate translation of lytic genes

Bacteriophages (phages) are obligate parasites that use host bacterial translation machinery to produce viral proteins. However, some phages have alternative genetic codes with reassigned stop codons that are predicted to be incompatible with bacterial translation systems. We analysed 9,422 phage genomes and found that stop-codon recoding has evolved in diverse clades of phages that infect bacteria present in both human and animal gut microbiota. Recoded stop codons are particularly over-represented in phage structural and lysis genes. We propose that recoded stop codons might function to prevent premature production of late-stage proteins. Stop-codon recoding has evolved several times in closely related lineages, which suggests that adaptive recoding can occur over very short evolutionary timescales.

Construction and Evaluation of Recombinant Attenuated Edwardsiella piscicida Vaccine (RAEV) Vector System Encoding Ichthyophthirius multifiliis (Ich) Antigen IAG52B

We have successfully designed and constructed a RAEV vector system with regulated-delayed attenuation in vivo attributes that synthesizes Ichthyophthirius multifiliis (Ich) protective antigen IAG52B to enable vaccination of fish susceptible to edwardsiellosis and white spot disease. The first feature of this vaccine delivery system is an Edwardsiella piscicida strain carrying genomic deletions of asdA. AsdA is an enzyme necessary for the synthesis of diaminopimelic acid (DAP), which is an essential component of the peptidoglycan layer of the cell wall of Gram-negative bacteria. asdA mutant strains have obligate growth requirements for DAP in the medium or a plasmid vector with the wild-type asdA gene enabling synthesis of DAP. This balanced-lethal plasmid vector-host system in E. piscicida enables as a second feature the synthesis of recombinant antigens to induce protective immunity against fish pathogens. Recombinant protective antigen IAG52B from the fish pathogen I. multifiliis was synthesized by RAEV strains harboring the AsdA⁺ plasmid pG8R8029. The third feature of this vaccine strain is a regulated-delayed attenuation in vivo phenotype that is based on the replacement of an arabinose-regulated araC P_araBAD cassette for the promoters of the fur and crp genes of E. piscicida such that the expression of these genes is dependent on arabinose provided during growth. Thus, following colonization, the Fur and Crp proteins stop being synthesized due to the lack of arabinose and attenuation is progressively achieved in vivo to prevent generation of diseases symptoms. Our vaccine strain χ16022 with the genotype ΔasdA10 ΔP_fur170::TT araC P_araBAD fur ΔP_crp68::TT araC P_araBAD crp contains the AsdA⁺ plasmid, pG8R8029, which encodes the IAG52B antigen. Vaccine strain χ16022(pG8R8029) is attenuated and induces systemic and mucosal IgM titer against E. piscicida and Ich in zebrafish. In addition, transcript levels of tnf-α, il-1β, il-6 and il-8 were significantly increased in different tissues of vaccinated zebrafish compared to unimmunized fish. Zebrafish vaccinated with χ16022(pG8R8029) showed 60% survival upon intracoelomic (i.c.) challenge with a lethal dose of virulent E. piscicida strain J118. Our RAEV system could be used as a generalized vaccine-vector system to protect teleost fish against multiple bacterial, viral and parasitic infectious diseases.

Programmed genome rearrangements in ciliates

Ciliates are a highly divergent group of unicellular eukaryotes with separate somatic and germline genomes found in distinct dimorphic nuclei. This characteristic feature is tightly linked to extremely laborious developmentally regulated genome rearrangements in the development of a new somatic genome/nuclei following sex. The transformation from germline to soma genome involves massive DNA elimination mediated by non-coding RNAs, chromosome fragmentation, as well as DNA amplification. In this review, we discuss the similarities and differences in the genome reorganization processes of the model ciliates Paramecium and Tetrahymena (class Oligohymenophorea), and the distantly related Euplotes, Stylonychia, and Oxytricha (class Spirotrichea).

Rewriting the Genetic Code

The genetic code-the language used by cells to translate their genomes into proteins that perform many cellular functions-is highly conserved throughout natural life. Rewriting the genetic code could lead to new biological functions such as expanding protein chemistries with noncanonical amino acids (ncAAs) and genetically isolating synthetic organisms from natural organisms and viruses. It has long been possible to transiently produce proteins bearing ncAAs, but stabilizing an expanded genetic code for sustained function in vivo requires an integrated approach: creating recoded genomes and introducing new translation machinery that function together without compromising viability or clashing with endogenous pathways. In this review, we discuss design considerations and technologies for expanding the genetic code. The knowledge obtained by rewriting the genetic code will deepen our understanding of how genomes are designed and how the canonical genetic code evolved.

Nuclear codon reassignments in the genomics era and mechanisms behind their evolution

The canonical genetic code ubiquitously translates nucleotide into peptide sequence with several alterations known in viruses, bacteria, mitochondria, plastids, and single-celled eukaryotes. A new hypothesis to explain genetic code changes, termed tRNA loss driven codon reassignment, has been proposed recently when the polyphyly of the yeast codon reassignment events has been uncovered. According to this hypothesis, the driving force for genetic code changes are tRNA or translation termination factor loss-of-function mutations or loss-of-gene events. The free codon can subsequently be captured by all tRNAs that have an appropriately mutated anticodon and are efficiently charged. Thus, codon capture most likely happens by near-cognate tRNAs and tRNAs whose anticodons are not part of the recognition sites of the respective aminoacyl-tRNA-synthetases. This hypothesis comprehensively explains the CTG codon translation as alanine in Pachysolen yeast together with the long known translation of the same codon as serine in Candida albicans and related species, and can also be applied to most other known reassignments.

Nuclear genetic codes with a different meaning of the UAG and the UAA codon

BACKGROUND

Departures from the standard genetic code in eukaryotic nuclear genomes are known for only a handful of lineages and only a few genetic code variants seem to exist outside the ciliates, the most creative group in this regard. Most frequent code modifications entail reassignment of the UAG and UAA codons, with evidence for at least 13 independent cases of a coordinated change in the meaning of both codons. However, no change affecting each of the two codons separately has been documented, suggesting the existence of underlying evolutionary or mechanistic constraints.

RESULTS

Here, we present the discovery of two new variants of the nuclear genetic code, in which UAG is translated as an amino acid while UAA is kept as a termination codon (along with UGA). The first variant occurs in an organism noticed in a (meta)transcriptome from the heteropteran Lygus hesperus and demonstrated to be a novel insect-dwelling member of Rhizaria (specifically Sainouroidea). This first documented case of a rhizarian with a non-canonical genetic code employs UAG to encode leucine and represents an unprecedented change among nuclear codon reassignments. The second code variant was found in the recently described anaerobic flagellate Iotanema spirale (Metamonada: Fornicata). Analyses of transcriptomic data revealed that I. spirale uses UAG to encode glutamine, similarly to the most common variant of a non-canonical code known from several unrelated eukaryotic groups, including hexamitin diplomonads (also a lineage of fornicates). However, in these organisms, UAA also encodes glutamine, whereas it is the primary termination codon in I. spirale. Along with phylogenetic evidence for distant relationship of I. spirale and hexamitins, this indicates two independent genetic code changes in fornicates.

CONCLUSIONS

Our study documents, for the first time, that evolutionary changes of the meaning of UAG and UAA codons in nuclear genomes can be decoupled and that the interpretation of the two codons by the cytoplasmic translation apparatus is mechanistically separable. The latter conclusion has interesting implications for possibilities of genetic code engineering in eukaryotes. We also present a newly developed generally applicable phylogeny-informed method for inferring the meaning of reassigned codons.

Genetic Codes with No Dedicated Stop Codon: Context-Dependent Translation Termination

The prevailing view of the nuclear genetic code is that it is largely frozen and unambiguous. Flexibility in the nuclear genetic code has been demonstrated in ciliates that reassign standard stop codons to amino acids, resulting in seven variant genetic codes, including three previously undescribed ones reported here. Surprisingly, in two of these species, we find efficient translation of all 64 codons as standard amino acids and recognition of either one or all three stop codons. How, therefore, does the translation machinery interpret a “stop” codon? We provide evidence, based on ribosomal profiling and “stop” codon depletion shortly before coding sequence ends, that mRNA 3′ ends may contribute to distinguishing stop from sense in a context-dependent manner. We further propose that such context-dependent termination/readthrough suppression near transcript ends enables genetic code evolution.

•

Alternative nuclear genetic codes continue to be discovered in ciliates

•

Genetic codes with stops and all their codons encoding standard amino acids exist

•

Transcript ends may distinguish stop codons as such in ambiguous genetic codes

•

The ability to resolve genetic code ambiguity may enable genetic code evolution

Programmed Rearrangement in Ciliates: Paramecium

Programmed genome rearrangements in the ciliate Paramecium provide a nice illustration of the impact of transposons on genome evolution and plasticity. During the sexual cycle, development of the somatic macronucleus involves elimination of ∼30% of the germline genome, including repeated DNA (e.g., transposons) and ∼45,000 single-copy internal eliminated sequences (IES). IES excision is a precise cut-and-close process, in which double-stranded DNA cleavage at IES ends depends on PiggyMac, a domesticated piggyBac transposase. Genome-wide analysis has revealed that at least a fraction of IESs originate from Tc/mariner transposons unrelated to piggyBac. Moreover, genomic sequences with no transposon origin, such as gene promoters, can be excised reproducibly as IESs, indicating that genome rearrangements contribute to the control of gene expression. How the system has evolved to allow elimination of DNA sequences with no recognizable conserved motif has been the subject of extensive research during the past two decades. Increasing evidence has accumulated for the participation of noncoding RNAs in epigenetic control of elimination for a subset of IESs, and in trans-generational inheritance of alternative rearrangement patterns. This chapter summarizes our current knowledge of the structure of the germline and somatic genomes for the model species Paramecium tetraurelia, and describes the DNA cleavage and repair factors that constitute the IES excision machinery. We present an overview of the role of specialized RNA interference machineries and their associated noncoding RNAs in the control of DNA elimination. Finally, we discuss how RNA-dependent modification and/or remodeling of chromatin may guide PiggyMac to its cognate cleavage sites.

Augmented genetic decoding: global, local and temporal alterations of decoding processes and codon meaning

The non-universality of the genetic code is now widely appreciated. Codes differ between organisms, and certain genes are known to alter the decoding rules in a site-specific manner. Recently discovered examples of decoding plasticity are particularly spectacular. These examples include organisms and organelles with disruptions of triplet continuity during the translation of many genes, viruses that alter the entire genetic code of their hosts and organisms that adjust their genetic code in response to changing environments. In this Review, we outline various modes of alternative genetic decoding and expand existing terminology to accommodate recently discovered manifestations of this seemingly sophisticated phenomenon.

Genetic code evolution reveals the neutral emergence of mutational robustness, and information as an evolutionary constraint

The standard genetic code (SGC) is central to molecular biology and its origin and evolution is a fundamental problem in evolutionary biology, the elucidation of which promises to reveal much about the origins of life. In addition, we propose that study of its origin can also reveal some fundamental and generalizable insights into mechanisms of molecular evolution, utilizing concepts from complexity theory. The first is that beneficial traits may arise by non-adaptive processes, via a process of "neutral emergence". The structure of the SGC is optimized for the property of error minimization, which reduces the deleterious impact of point mutations. Via simulation, it can be shown that genetic codes with error minimization superior to the SGC can emerge in a neutral fashion simply by a process of genetic code expansion via tRNA and aminoacyl-tRNA synthetase duplication, whereby similar amino acids are added to codons related to that of the parent amino acid. This process of neutral emergence has implications beyond that of the genetic code, as it suggests that not all beneficial traits have arisen by the direct action of natural selection; we term these "pseudaptations", and discuss a range of potential examples. Secondly, consideration of genetic code deviations (codon reassignments) reveals that these are mostly associated with a reduction in proteome size. This code malleability implies the existence of a proteomic constraint on the genetic code, proportional to the size of the proteome (P), and that its reduction in size leads to an "unfreezing" of the codon - amino acid mapping that defines the genetic code, consistent with Crick's Frozen Accident theory. The concept of a proteomic constraint may be extended to propose a general informational constraint on genetic fidelity, which may be used to explain variously, differences in mutation rates in genomes with differing proteome sizes, differences in DNA repair capacity and genome GC content between organisms, a selective pressure in the evolution of sexual reproduction, and differences in translational fidelity. Lastly, the utility of the concept of an informational constraint to other diverse fields of research is explored.

Genomic characterization of variable surface antigens reveals a telomere position effect as a prerequisite for RNA interference-mediated silencing in Paramecium tetraurelia

Antigenic or phenotypic variation is a widespread phenomenon of expression of variable surface protein coats on eukaryotic microbes. To clarify the mechanism behind mutually exclusive gene expression, we characterized the genetic properties of the surface antigen multigene family in the ciliate Paramecium tetraurelia and the epigenetic factors controlling expression and silencing. Genome analysis indicated that the multigene family consists of intrachromosomal and subtelomeric genes; both classes apparently derive from different gene duplication events: whole-genome and intrachromosomal duplication. Expression analysis provides evidence for telomere position effects, because only subtelomeric genes follow mutually exclusive transcription. Microarray analysis of cultures deficient in Rdr3, an RNA-dependent RNA polymerase, in comparison to serotype-pure wild-type cultures, shows cotranscription of a subset of subtelomeric genes, indicating that the telomere position effect is due to a selective occurrence of Rdr3-mediated silencing in subtelomeric regions. We present a model of surface antigen evolution by intrachromosomal gene duplication involving the maintenance of positive selection of structurally relevant regions. Further analysis of chromosome heterogeneity shows that alternative telomere addition regions clearly affect transcription of closely related genes. Consequently, chromosome fragmentation appears to be of crucial importance for surface antigen expression and evolution. Our data suggest that RNAi-mediated control of this genetic network by trans-acting RNAs allows rapid epigenetic adaptation by phenotypic variation in combination with long-term genetic adaptation by Darwinian evolution of antigen genes.

IMPORTANCE

Alternating surface protein structures have been described for almost all eukaryotic microbes, and a broad variety of functions have been described, such as virulence factors, adhesion molecules, and molecular camouflage. Mechanisms controlling gene expression of variable surface proteins therefore represent a powerful tool for rapid phenotypic variation across kingdoms in pathogenic as well as free-living eukaryotic microbes. However, the epigenetic mechanisms controlling synchronous expression and silencing of individual genes are hardly understood. Using the ciliate Paramecium tetraurelia as a (epi)genetic model, we showed that a subtelomeric gene position effect is associated with the selective occurrence of RNAi-mediated silencing of silent surface protein genes, suggesting small interfering RNA (siRNA)-mediated epigenetic cross talks between silent and active surface antigen genes. Our integrated genomic and molecular approach discloses the correlation between gene position effects and siRNA-mediated trans-silencing, thus providing two new parameters for regulation of mutually exclusive gene expression and the genomic organization of variant gene families.

Identification of amino acids responsible for stop codon recognition for polypeptide chain release factor

One factor involved in eukaryotic translation termination is class 1 release factor in eukaryotes (eRF1), which functions to decode stop codons. Variant code species, such as ciliates, frequently exhibit altered stop codon recognition. Studies revealed that some class-specific residues in the eRF1 N-terminal domain are responsible for stop codon reassignment in ciliates. Here, we investigated the effects on stop codon recognition of chimeric eRF1s containing the N-terminal domain of Euplotes octocarinatus and Blepharisma japonicum eRF1 fused to Saccharomyces cerevisiae M and C domains using dual luciferase read-through assays. Mutation of class-specific residues in different eRF1 classes was also studied to identify key residues and motifs involved in stop codon decoding. As expected, our results demonstrate that 3 pockets within the eRF1 N-terminal domain were involved in decoding stop codon nucleotides. However, allocation of residues to each pocket was revalued. Our data suggest that hydrophobic and class-specific surface residues participate in different functions: modulation of pocket conformation and interaction with stop codon nucleotides, respectively. Residues conserved across all eRF1s determine the relative orientation of the 3 pockets according to stop codon nucleotides. However, quantitative analysis of variant ciliate and yeast eRF1 point mutants did not reveal any correlation between evolutionary conservation of class-specific residues and termination-related functional specificity and was limited in elucidating a detailed mechanism for ciliate stop codon reassignment. Thus, based on isolation of suppressor tRNAs from Euplotes and Tetrahymena, we propose that stop codon reassignment in ciliates may be controlled by cooperation between eRF1 and suppressor tRNAs.

Identification of eRF1 residues that play critical and complementary roles in stop codon recognition

The initiation and elongation stages of translation are directed by codon-anticodon interactions. In contrast, a release factor protein mediates stop codon recognition prior to polypeptide chain release. Previous studies have identified specific regions of eukaryotic release factor one (eRF1) that are important for decoding each stop codon. The cavity model for eukaryotic stop codon recognition suggests that three binding pockets/cavities located on the surface of eRF1's domain one are key elements in stop codon recognition. Thus, the model predicts that amino acid changes in or near these cavities should influence termination in a stop codon-dependent manner. Previous studies have suggested that the TASNIKS and YCF motifs within eRF1 domain one play important roles in stop codon recognition. These motifs are highly conserved in standard code organisms that use UAA, UAG, and UGA as stop codons, but are more divergent in variant code organisms that have reassigned a subset of stop codons to sense codons. In the current study, we separately introduced TASNIKS and YCF motifs from six variant code organisms into eRF1 of Saccharomyces cerevisiae to determine their effect on stop codon recognition in vivo. We also examined the consequences of additional changes at residues located between the TASNIKS and YCF motifs. Overall, our results indicate that changes near cavities two and three frequently mediated significant effects on stop codon selectivity. In particular, changes in the YCF motif, rather than the TASNIKS motif, correlated most consistently with variant code stop codon selectivity.

A characterization of Tetrahymena mRNA by in vitro translation: The effects of culture growth on the recruitment of poly (A)(+) and poly (A)(-) RNA

Messenger RNA was extracted from exponentially growing and from stationary cultures of Tetrahymena thermophila. It was separated into polyadenylated and non-polyadenylated fractions which were used as templates in a rabbit reticulocyte protein synthesizing system. The translated proteins were analysed by one and by two dimensional electrophoresis. Our experiments were intended to answer the question to which extent the abundance and the specificity of mRNA facilitates or accompanies the passage of cells through one culture growth cycle. As illustrated by the identification of 113 proteins very few differences between translated messages accompany the transition to the stationary phase, the most obvious feature being a change in the intracellular location of translation activities. These data are discussed with special reference to the prevailing occupation of Tetrahymena which is biomass production.

Isolation and characterization of the actin gene from Tetrahymena thermophila

The macronucleus of Tetrahymena thermophila contains a single actin gene. We have isolated this gene from a partial plasmid library by using the yeast actin gene as a probe. The nucleotide sequence of the gene has been determined and the amino acid sequence of the potential protein deduced. The encoded protein is 375 amino acids long, one amino acid longer than the yeast actin. It is one of the most divergent actins sequenced to date, being only 75% homologous to yeast actin. Unlike the actin genes from most other organisms, it does not contain introns. The coding region contains TAA and TAG codons; the translation termination codon is TGA. Comparison of the amino acid sequence of the Tetrahymena actin with that of actins from other organisms suggests that TAG may code for glutamic acid. The gene is transcribed from multiple initiation sites between 57 and 98 nucleotides upstream of the translation start codon. The 5' flanking region is very A+T-rich and contains numerous "TATA-like" sequences upstream of the transcription start sites.

Identification and characterization of a putative agglutination/immobilization antigen on the surface ofCryptocaryon irritans

The ciliated protozoanCryptocaryon irritans, a parasite of seawater fishes, was found to express an antigen that elicits antibodies in rabbits and tiger puffer (Takifugu ruburipes). Serum from rabbits and fish immunized with theronts had agglutination/immobilization activity against therontsin vitro; fish serum antibody levels (measured by enzyme-linked immunosorbent assays: ELISA) correlated with this activity. Anti-theront antibody levels in fish were significantly higher in the immunized group as compared with control fish at 2 weeks after booster immunization (injection of bovine serum albumin; Student'st-test,P<0·01). Biochemical analyses indicated that a Triton X-114-soluble 32 kDa theront integral membrane protein may be the agglutination/immobilization antigen. Indirect immunofluorescence staining of theronts suggested that this 32 kDa antigen was expressed on the surface of cilia. The full-length 32 kDa antigen cDNA contained 1147 basepairs, encoding a 328-amino acid protein including hydrophobic N- and C-termini. As withTetrahymenaandParameciumspp., TAA and TAG appear to be used as glutamine codons in the 32 kDa antigen gene.

Antigenic Variation in Ciliates: Antigen Structure, Function, Expression

In the past decades, the major focus of antigen variation research has been on parasitic protists. However, antigenic variation occurs also in free-living protists. The antigenic systems of the ciliates Paramecium and Tetrahymena have been studied for more than 100 yr. In spite of different life strategies and distant phylogenetic relationships of free-living ciliates and parasitic protists, their antigenic systems have features in common, such as the presence of repeated protein motifs and multigene families. The function of variable surface antigens in free-living ciliates is still unknown. Up to now no detailed monitoring of antigen expression in free-living ciliates in natural habitats has been performed. Unlike stochastic switching in parasites, antigen expression in ciliates can be directed, e.g. by temperature, which holds great advantages for research on the expression mechanism. Regulated expression of surface antigens occurs in an exclusive way and the responsible mechanism is complex, involving both transcriptional and post-transcriptional features. The involvement of homology-dependent effects has been proposed several times but has not been proved yet.

Variety of Serotypes of Paramecium primaurelia: Single Epitopes are Responsible for Immunological Differentiation

Paramecium primaurelia expresses three major types of surface antigens. We report here the identification of the gene for serotype S, which completes the sequence data of expressed serotypes of P. primaurelia. The complete open reading frame of surface antigen S was identified using a novel technique, based upon the presence of conservative regions in the non-coding areas of the multigene family. We were able to isolate the 7194-bp-long open reading frame from the macronuclear DNA for Serotype 156S. The corresponding mRNA was detected in the two serotype S-expressing stocks, 60 and 156, of P. primaurelia, which clarifies that both stocks are using the same S allele. Comparisons of the nucleic acid and the deduced amino-acid sequence showed high identity to surface antigen 51B of P. tetraurelia, sufficient to cause an immunological cross-reaction in vivo. Immunologically relevant epitopes in vivo were identified in the central regions of the genes, constructed of nearly perfect tandem repeats.

Newly sequenced eRF1s from ciliates: the diversity of stop codon usage and the molecular surfaces that are important for stop codon interactions

The genetic code of nuclear genes in some ciliates was found to differ from that of other organisms in the assignment of UGA, UAG, and UAA codons, which are normally assigned as stop codons. In some ciliate species, the universal stop codons UAA and UAG instead encode glutamine. In some other ciliates, the universal stop codon UGA appears to be translated as cysteine or tryptophan. Eukaryotic release factor 1 (eRF1) is a key protein in stop codon recognition, thus, the protein is believed to play an important role in the stop codon reassignment in ciliates. We have cloned, sequenced, and analyzed the cDNA of eRF1 from four ciliate species of three different classes: Karyorelictea (Loxodes striatus), Heterotrichea (Blepharisma musculus), and Litostomatea (Didinium nasutum, Dileptus margaritifer). Phylogenetic analysis of these eRF1s supports the hypothesis that the genetic code in ciliates has deviated independently several times from the universal genetic code, and that different ciliate eRF1s may have undergone different processes to change the codon specificity. Using computational methods, we have also suggested areas on the surface of eRF1s that are important for stop codon recognition in ciliate eRF1s.

Computational identification and sequence analysis of stop codon readthrough genes in Oryza sativa

Using an approach based on the Readthrough Candidate Extraction System (RCES), we extracted 111 candidates from 9620 gene sequences of rice. The results of homology search and sequence analysis demonstrated that these candidates included actual readthrough genes that would be important for further investigating the mechanism of translation termination regulated by readthrough event, and could also give some useful clues for functional genome annotation. Between the candidates and non-candidates of gene sequences in rice, there exist significant base biases at the positions surrounding the stop codons. These positions, especially both -1 and +4, are referred to as part of an extended stop signal. In candidates, G at position -1, and G or C at position +4 are much more favored than that in non-candidates. Both stop sequence patterns, GUAGC and GUGAG, might drive high readthrough efficiency in rice. Secondary structure analysis revealed that the -1 and +1 amino acids around the first stop codon of candidates have a strong bias toward arginine, particularly the +1 position (20.7%), which indicated that the amino acids at the readthrough region being frequently located in the hydrophilic region of beta-turn might be a determinant for efficient translation termination or not.

Adenylyl cyclases from Plasmodium, Paramecium and Tetrahymena are novel ion channel/enzyme fusion proteins

In Paramecium, cAMP formation is stimulated by a potassium conductance, which is an intrinsic property of the adenylyl cyclase. We cloned a full-length cDNA and several gDNA fragments from Paramecium and Tetrahymena coding for adenylyl cyclases with a novel domain composition. A putative N-terminal ion channel domain contains a canonical S4 voltage-sensor and a canonical potassium pore-loop located C-terminally after the last transmembrane span on the cytoplasmic side. The adenylyl cyclase catalyst is C-terminally located. DNA microinjection of a green fluorescent protein (GFP)-tagged construct into the macronucleus of Paramecium resulted in ciliary localization of the expressed protein. An identical gene coding for an ion-channel adenylyl cyclase was cloned from the malaria parasite Plasmodium falciparum. Expression of the catalytic domain of the latter in Sf9 cells yielded an active homodimeric adenylyl cyclase. The occurrence of this highly unique subtype of adenylyl cyclase appears to be restricted to ciliates and apicomplexa.

Highly divergent actins from karyorelictean, heterotrich, and litostome ciliates

We have cloned, sequenced, and characterized cDNA of actins from five ciliate species of three different classes of the phylum Ciliophora: Karyorelictea (Loxodes striatus), Heterotrichea (Blepharisma japonicum, Blepharisma musculus), and Litostomatea (Didinium nasutum, Dileptus margaritifer). Loxodes striatus uses UGA as the stop codon and has numerous in-frame UAA and UAG, which are translated into glutamine. The other four species use UAA as the stop codon and have no in-frame UAG nor UGA. The putative amino acid sequences of the newly determined actin genes were found to be highly divergent as expected from previous findings of other ciliate actins. These sequences were also highly divergent from other ciliate actins, indicating that actin genes are highly diverse even within the phylum Ciliophora. Phylogenetic analysis showed high evolutionary rate of ciliate actins. Our results suggest that the evolutionary rate was accelerated because of the differences in molecular interactions.

Molecular identification of a calcium-inhibited catalytic subunit of casein kinase type 2 from Paramecium tetraurelia

We have previously described the occurrence in Paramecium of a casein kinase (CK) activity (EC 2.7.1.37) with some unusual properties, including inhibition by Ca(2+) (R. Kissmehl, T. Treptau, K. Hauser, and H. Plattner, FEBS Lett. 402:227-235, 1995). We now have cloned four genes, PtCK2alpha1 to PtCK2alpha4, all of which encode the catalytic alpha subunit of type 2 CK (CK2) with calculated molecular masses ranging from 38.9 to 39.4 kDa and pI values ranging from 8.8 to 9.0. They can be classified into two groups, which differ from each other by 28% on the nucleotide level and by 18% on the derived amino acid level. One of them, PtCK2alpha3, has been expressed in Escherichia coli and characterized in vitro. As we also have observed with the isolated CK, the recombinant protein preferentially phosphorylates casein but also phosphorylates some Paramecium-specific substrates, including the exocytosis-sensitive phosphoprotein pp63/parafusin. Characteristically, Ca(2+) inhibits the phosphorylation at elevated concentrations occurring during stimulation of a cell. Reconstitution with a recombinant form of the regulatory subunit from Xenopus laevis, XlCK2beta, confirms Ca(2+) sensitivity also under conditions of autophosphorylation. This is unusual for CK2 but correlates with the presence of two EF-hand calcium-binding motifs, one of which is located in the N-terminal segment essential for constitutive activity, as well as with an aberrant composition of normally basic domains recognizing acidic substrate domains. Immunogold localization reveals a considerable enrichment in the outermost cell cortex layers, excluding cilia. We discuss a potential role of this Ca(2+)-inhibited PtCK2alpha species in a late step of signal transduction.

Regulation of gene expression by stop codon recoding: selenocysteine

The regulation of gene expression at the translational level not only allows for rapid changes in specific protein levels but also provides an opportunity to alter codon specificity. For the incorporation of selenocysteine (Sec) into protein, the UGA codon is transformed from one that signals translation termination to one specific for Sec. This review provides a look at Sec incorporation from the perspective of the individual steps involved in protein synthesis: initiation, elongation and termination. The roles of the factors known to be required for Sec incorporation are considered in the context of each step in translation including structural modeling of the differences between the standard elongation factor eEF1A and the Sec-specific counterpart, eEFSec.

Stop codon selection in eukaryotic translation termination: comparison of the discriminating potential between human and ciliate eRF1s

During eukaryotic translation termination, eRF1 responds to three stop codons. However, in ciliates with variant genetic codes, only one or two codons function as a stop signal. To localize the region of ciliate eRF1 implicated in stop codon discrimination, we have constructed ciliate-human hybrid eRF1s by swapping regions of human eRF1 for the equivalent region of ciliate Euplotes eRF1. We have examined the formation of a cross-link between recombinant eRF1s and mRNA analogs containing the photoactivable 4-thiouridine (s(4)U) at the first position of stop and control sense codons. With human eRF1, this cross-link can be detected only when either stop or UGG codons are located in the ribosomal A site. Here we show that the cross-link of the Euplotes-human hybrid eRF1 is restricted to mRNAs containing UAG and UAA codons, and that the entire N-terminal domain of Euplotes eRF1 is involved in discriminating against UGA and UGG. On the basis of these results, we discuss the steps of the selection process that determine the accuracy of stop codon recognition in eukaryotes.

Isolation and expression of two genes encoding eukaryotic release factor 1 from Paramecium tetraurelia

Paramecium tetraurelia, like some other ciliate species, uses an alternative nuclear genetic code where UAA and UAG are translated as glutamine and UGA is the only stop codon. It has been postulated that the use of stop codons as sense codons is dependent on the presence of specific tRNAs and on modification of eukaryotic release factor one (eRF1), a factor involved in stop codon recognition during translation termination. We describe here the isolation and characterisation of two genes, eRF1-a and eRF1 b, coding for eRF1 in P. tetraurelia. The two genes are very similar, both in genomic organization and in sequence, and might result from a recent duplication event. The two coding sequences are 1,314 nucleotides long, and encode two putative proteins of 437 amino acids with 98.5% identity. Interestingly, when compared with the eRF1 sequences either of ciliates having the same variant genetic code, or of other eukaryotes, the eRF1 of P. tetraurelia exhibits significant differences in the N-terminal region, which is thought to interact with stop codons. We discuss here the consequences of these changes in the light of recent models proposed to explain the mechanism of stop codon recognition in eukaryotes. Besides, analysis of the expression of the two genes by Northern blotting and primer extension reveals that these genes exhibit a differential expression during vegetative growth and autogamy.

Random sequencing of Paramecium somatic DNA

We report a random survey of 1 to 2% of the somatic genome of the free-living ciliate Paramecium tetraurelia by single-run sequencing of the ends of plasmid inserts. As in all ciliates, the germ line genome of Paramecium (100 to 200 Mb) is reproducibly rearranged at each sexual cycle to produce a somatic genome of expressed or potentially expressed genes, stripped of repeated sequences, transposons, and AT-rich unique sequence elements limited to the germ line. We found the somatic genome to be compact (>68% coding, estimated from the sequence of several complete library inserts) and to feature uniformly small introns (18 to 35 nucleotides). This facilitated gene discovery: 722 open reading frames (ORFs) were identified by similarity with known proteins, and 119 novel ORFs were tentatively identified by internal comparison of the data set. We determined the phylogenetic position of Paramecium with respect to eukaryotes whose genomes have been sequenced by the distance matrix neighbor-joining method by using random combined protein data from the project. The unrooted tree obtained is very robust and in excellent agreement with accepted topology, providing strong support for the quality and consistency of the data set. Our study demonstrates that a random survey of the somatic genome of Paramecium is a good strategy for gene discovery in this organism.

Identification, characterization, and complete amino acid sequence of the conjugation-inducing glycoprotein (blepharmone) in the ciliate Blepharisma japonicum

Conjugation in Blepharisma japonicum is induced by interaction between complementary mating-types I and II, which excrete blepharmone (gamone 1) and blepharismone (gamone 2), respectively. Gamone 1 transforms type II cells such that they can unite, and gamone 2 similarly transforms type I cells. Moreover, each gamone promotes the production of the other gamone. Gamone 2 has been identified as calcium-3-(2'-formylamino-5'-hydroxy-benzoyl) lactate and has been synthesized chemically. Gamone 1 was isolated and characterized as a glycoprotein of 20-30 kDa containing 175 amino acids and 6 sugars. However, the amino acid sequence and arrangement of sugars in this gamone are still unknown. To determine partial amino acid sequences of gamone 1, we established a method of isolation based on the finding that this glycoprotein can be concentrated by a Con A affinity column. Gamone 1 is extremely unstable and loses its biological activity once adsorbed to any of the columns that we tested. By using a Con A affinity column and native PAGE, we detected a 30-kDa protein corresponding to gamone 1 activity and determined the partial amino acid sequences of the four peptides. To isolate gamone 1 cDNA, we isolated mRNA from mating-type I cells stimulated by synthetic gamone 2 and then performed rapid amplification of cDNA ends procedures by using gene-specific primers and cloned cDNA of gamone 1. The cDNA sequence contains an ORF of 305 amino acids and codes a possibly novel protein. We also estimated the arrangement of sugars by comparing the affinity to various lectin columns.

K(+)-channel transgenes reduce K(+) currents in Paramecium, probably by a post-translational mechanism

PAK11 is 1 of more than 15 members in a gene family that encodes K(+)-channel pore-forming subunits in Paramecium tetraurelia. Microinjection of PAK11 DNA into macronuclei of wild-type cells results in clonal transformants that exhibit hyperexcitable swimming behaviors reminiscent of certain loss-of-K(+)-current mutants. PAK2, a distant homolog of PAK11, does not have the same effect. But PAK1, a close homolog of PAK11, induces the same hyperexcitability. Cutting the PAK11 open reading frame (ORF) with restriction enzymes before injection removes this effect entirely. Microinjection of PAK11 ORF flanked by the calmodulin 5' and 3' UTRs also induces the same hyperexcitable phenotype. Direct examination of transformed cells under voltage clamp reveals that two different Ca(2+)-activated K(+)-specific currents are reduced in amplitude. This reduction does not correlate with a deficit of PAK11 message, since RNA is clearly produced from the injected transgenes. Insertion of a single nucleotide at the start of the PAK11 ORF does not affect the RNA level but completely abolishes the phenotypic transformation. Thus, the reduction of K(+) currents by the expression of the K(+)-channel transgenes reported here is likely to be the consequence of a post-translational event. The complexity of behavioral changes, possible mechanisms, and implications in Paramecium biology are discussed.

Biological and molecular characterization of cellular differentiation in Tetrahymena vorax: a potential biocontrol protozoan

Tetrahymena vorax (T. vorax) is an indigenous fresh water protozoan with the natural biological potential to maintain a specific aquatic microbial flora by ingesting and eliminating specific microorganism. To investigate the molecular mechanisms controlling Tetrahymena vorax (T. vorax) cellular differentiation from a small-mouth vegetative cell to a voracious large-mouth carnivore capable of ingesting prey ciliates and bacteria from aquatic environments, we use DNA subtraction and gene discovery techniques to identify and isolate T. vorax differentiation-specific genes. The physiological necessity for one newly discovered gene, SUBII-TG, was determined in vivo using an antisense oligonucleotide directed against the 5' SUBII-TG DNA sequence. The barriers to delivering antisense oligonucleotides to the cytoplasm of T. vorax were circumvented by employing a new but simple procedure of processing the oligonucleotide with the differentiation stimulus, stomatin. In these studies, the antisense oligonucleotide down-regulated SUBII-TG mRNA expression, and blocked differentiation and ingestion of prey ciliates. The ability to down-regulate SUBII-TG expression with the antisense oligonucleotide suggests that the molecular mechanisms controlling the natural biological activities of T. vorax can be manipulated to further study its cellular differentiation and potential as a biocontrol microorganism.

Nuclear-encoded, plastid-targeted genes suggest a single common origin for apicomplexan and dinoflagellate plastids

The phylum Apicomplexa encompasses a large number of intracellular protozoan parasites, including the causative agents of malaria (Plasmodium), toxoplasmosis (Toxoplasma), and many other human and animal diseases. Apicomplexa have recently been found to contain a relic, nonphotosynthetic plastid that has attracted considerable interest as a possible target for therapeutics. This plastid is known to have been acquired by secondary endosymbiosis, but when this occurred and from which type of alga it was acquired remain uncertain. Based on the molecular phylogeny of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) genes, we provide evidence that the apicomplexan plastid is homologous to plastids found in dinoflagellates-close relatives of apicomplexa that contain secondary plastids of red algal origin. Surprisingly, apicomplexan and dinoflagellate plastid-targeted GAPDH sequences were also found to be closely related to the plastid-targeted GAPDH genes of heterokonts and cryptomonads, two other groups that contain secondary plastids of red algal origin. These results address several outstanding issues: (1) apicomplexan and dinoflagellate plastids appear to be the result of a single endosymbiotic event which occurred relatively early in eukaryotic evolution, also giving rise to the plastids of heterokonts and perhaps cryptomonads; (2) apicomplexan plastids are derived from a red algal ancestor; and (3) the ancestral state of apicomplexan parasites was photosynthetic.

The molecular basis of nuclear genetic code change in ciliates

BACKGROUND

The nuclear genetic code has changed in several lineages of ciliates. These changes, UAR to glutamine and UGA to cysteine, imply that eukaryotic release factor 1 (eRF1), the protein that recognizes stop codons and terminates translation, changes specificity. Here we test whether changes in eRF1 drive genetic code evolution.

RESULTS

Database sequence analysis reveals numerous genetic code alterations in ciliates, including UGA --> tryptophan in Blepharisma americanum and the distantly related Colpoda. We sequenced eRF1 from four ciliates: B. americanum, a heterotrich that independently derived the same eRF1 specificity as Euplotes, and three spirotrichs, Stylonychia lemnae, S. mytilus, and Oxytricha trifallax, that independently derived the same genetic code as Tetrahymena (UAR --> glutamine). Distantly related ciliates with similar codes show characteristic changes in eRF1. We used a sliding window analysis to test associations between changes in specific eRF1 residues and changes in the genetic code. The regions of eRF1 that display convergent substitutions are identical to those identified in a recently reported nonsense suppression mutant screen in yeast.

CONCLUSIONS

Genetic code change by stop codon reassignment is surprisingly frequent in ciliates, with UGA --> tryptophan occurring twice independently. This is the first description of this code, previously found only in bacteria and mitochondria, in a eukaryotic nuclear genome. eRF1 has evolved strikingly convergently in lineages with variant genetic codes. The strong concordance with biochemical data indicates that our methodology may be generally useful for detecting molecular determinants of biochemical changes in evolution.

Molecular mechanism of stop codon recognition by eRF1: a wobble hypothesis for peptide anticodons

We propose that the amino acid residues 57/58 and 60/61 of eukaryotic release factors (eRF1s) (counted from the N-terminal Met of human eRF1) are responsible for stop codon recognition in protein synthesis. The proposal is based on amino acid exchanges in these positions in the eRF1s of two ciliates that reassigned one or two stop codons to sense codons in evolution and on the crystal structure of human eRF1. The proposed mechanism of stop codon recognition assumes that the amino acid residues 57/58 interact with the second and the residues 60/61 with the third position of a stop codon. The fact that conventional eRF1s recognize all three stop codons but not the codon for tryptophan is attributed to the flexibility of the helix containing these residues. We suggest that the helix is able to assume a partly relaxed or tight conformation depending on the stop codon recognized. The restricted codon recognition observed in organisms with unconventional eRF1s is attributed mainly to the loss of flexibility of the helix due to exchanged amino acids.