Role of arginine kinase in Paramecium tetraurelia (Ciliophora, Peniculida): Subcellular localization of AK3 and phosphoarginine shuttle system in cilia.

The ciliate Paramecium tetraurelia has four arginine kinase genes (AK1, AK2, AK3, and AK4). Of these genes, only AK3 has a signal sequence for farnesylation, a post-translational modification that enables anchoring of the modified enzyme to the ciliary membrane. To confirm this modification, AK3 was synthesized using a cell-free protein synthesis system and the peptide masses were analyzed using peptide mass fingerprinting (PMF). The PMF analysis indicated that the C-terminal peptide of AK3 is farnesylated. Thus, AK3 can be farnesylated under physiologically appropriate conditions. To determine the subcellular localization of P. tetraurelia AK3, Western blot analysis was performed using an AK3 polyclonal antibody for the proteins extracted from intact cells and ciliary fractions. When extraction was performed using Triton X-100, AK3 was detected the ciliary fraction. This result suggested that the ciliary fraction contains AK3. In addition, we investigated the role of P. tetraurelia AKs in ciliary movement using the feeding RNA interference method. The swimming velocity of AK1- and AK3-silenced cells was significantly reduced to half the value of that control cells. In summary, P. tetraurelia AK3 is likely to be located in the ciliary membrane and influences swimming velocity, presumably through the phosphoarginine shuttle system present in cilia.

Characterization of four arginine kinases in the ciliate Paramecium tetraurelia: Investigation on the substrate inhibition mechanism.

The ciliate Paramecium tetraurelia contains four arginine kinase genes (AK1-4). We detected cDNA for only three of the AKs (AK1-3) via PCR. Recombinant AK1-4 were expressed in Escherichia coli and their kinetics parameters determined. AK3 showed typical substrate inhibition toward arginine, and enzymatic activity markedly decreased when arginine concentration increased. This is the first example of substrate inhibition in wild-type phosphagen kinases. To explore the substrate inhibition mechanism, site-directed mutations were generated, targeting the amino acid sequence D-D-S-Q-V at positions 77-81 in P. tetraurelia AK3. Among the mutants, substrate inhibition was lost remarkably in the S79A mutant. In spite of high amino acid sequence identity (91%) between P. tetraurelia AK3 and AK4, the enzymatic activity of AK4 was less by 3% than that of AK3. We noticed that the conservative G298 was unusually replaced by R in P. tetraurelia AK4, and we constructed two mutants, R298G/AK4 and G298R/AK3. Enzymatic activity of the former mutant was comparable with that of the wild-type AK3, whereas that of the latter mutant was dramatically reduced. Thus, we concluded that the significantly low activity of P. tetraurelia AK4 is due to the residue R298.

Helical arrays of U-shaped ATP synthase dimers form tubular cristae in ciliate mitochondria.

F1Fo-ATP synthases are universal energy-converting membrane protein complexes that synthesize ATP from ADP and inorganic phosphate. In mitochondria of yeast and mammals, the ATP synthase forms V-shaped dimers, which assemble into rows along the highly curved ridges of lamellar cristae. Using electron cryotomography and subtomogram averaging, we have determined the in situ structure and organization of the mitochondrial ATP synthase dimer of the ciliate Paramecium tetraurelia. The ATP synthase forms U-shaped dimers with parallel monomers. Each complex has a prominent intracrista domain, which links the c-ring of one monomer to the peripheral stalk of the other. Close interaction of intracrista domains in adjacent dimers results in the formation of helical ATP synthase dimer arrays, which differ from the loose dimer rows in all other organisms observed so far. The parameters of the helical arrays match those of the cristae tubes, suggesting the unique features of the P. tetraurelia ATP synthase are directly responsible for generating the helical tubular cristae. We conclude that despite major structural differences between ATP synthase dimers of ciliates and other eukaryotes, the formation of ATP synthase dimer rows is a universal feature of mitochondria and a fundamental determinant of cristae morphology.

Functional diversification of Dicer-like proteins and small RNAs required for genome sculpting.

In eukaryotes, small RNAs (sRNAs) have key roles in development, gene expression regulation, and genome integrity maintenance. In ciliates, such as Paramecium, sRNAs form the heart of an epigenetic system that has evolved from core eukaryotic gene silencing components to selectively target DNA for deletion. In Paramecium, somatic genome development from the germline genome accurately eliminates the bulk of typically gene-interrupting, noncoding DNA. We have discovered an sRNA class (internal eliminated sequence [IES] sRNAs [iesRNAs]), arising later during Paramecium development, which originates from and precisely delineates germline DNA (IESs) and complements the initial sRNAs ("scan" RNAs [scnRNAs]) in targeting DNA for elimination. We show that whole-genome duplications have facilitated successive differentiations of Paramecium Dicer-like proteins, leading to cooperation between Dcl2 and Dcl3 to produce scnRNAs and to the production of iesRNAs by Dcl5. These innovations highlight the ability of sRNA systems to acquire capabilities, including those in genome development and integrity.

Genome-wide analysis of the phosphoinositide kinome from two ciliates reveals novel evolutionary links for phosphoinositide kinases in eukaryotic cells.

BACKGROUND: The complexity of phosphoinositide signaling in higher eukaryotes is partly due to expansion of specific families and types of phosphoinositide kinases (PIKs) that can generate all phosphoinositides via multiple routes. This is particularly evident in the PI3Ks and PIPKs, and it is considered an evolutionary trait associated with metazoan diversification. Yet, there are limited comprehensive studies on the PIK repertoire of free living unicellular organisms. METHODOLOGY/PRINCIPAL FINDINGS: We undertook a genome-wide analysis of putative PIK genes in two free living ciliated cells, Tetrahymena and Paramecium. The Tetrahymena thermophila and Paramecium tetraurelia genomes were probed with representative kinases from all families and types. Putative homologs were verified by EST, microarray and deep RNA sequencing database searches and further characterized for domain structure, catalytic efficiency, expression patterns and phylogenetic relationships. In total, we identified and characterized 22 genes in the Tetrahymena thermophila genome and 62 highly homologues genes in Paramecium tetraurelia suggesting a tight evolutionary conservation in the ciliate lineage. Comparison to the kinome of fungi reveals a significant expansion of PIK genes in ciliates. CONCLUSIONS/SIGNIFICANCE: Our study highlights four important aspects concerning ciliate and other unicellular PIKs. First, ciliate-specific expansion of PI4KIII-like genes. Second, presence of class I PI3Ks which, at least in Tetrahymena, are associated with a metazoan-type machinery for PIP3 signaling. Third, expansion of divergent PIPK enzymes such as the recently described type IV transmembrane PIPKs. Fourth, presence of possible type II PIPKs and presumably inactive PIKs (hence, pseudo-PIKs) not previously described. Taken together, our results provide a solid framework for future investigation of the roles of PIKs in ciliates and indicate that novel functions and novel regulatory pathways of phosphoinositides may be more widespread than previously thought in unicellular organisms.

Extraordinary genome stability in the ciliate Paramecium tetraurelia.

Mutation plays a central role in all evolutionary processes and is also the basis of genetic disorders. Established base-substitution mutation rates in eukaryotes range between ∼5 × 10(-10) and 5 × 10(-8) per site per generation, but here we report a genome-wide estimate for Paramecium tetraurelia that is more than an order of magnitude lower than any previous eukaryotic estimate. Nevertheless, when the mutation rate per cell division is extrapolated to the length of the sexual cycle for this protist, the measure obtained is comparable to that for multicellular species with similar genome sizes. Because Paramecium has a transcriptionally silent germ-line nucleus, these results are consistent with the hypothesis that natural selection operates on the cumulative germ-line replication fidelity per episode of somatic gene expression, with the germ-line mutation rate per cell division evolving downward to the lower barrier imposed by random genetic drift. We observe ciliate-specific modifications of widely conserved amino acid sites in DNA polymerases as one potential explanation for unusually high levels of replication fidelity.

Metabolic engineering of Saccharomyces cerevisiae for production of Eicosapentaenoic Acid, using a novel {Delta}5-Desaturase from Paramecium tetraurelia.

Very-long-chain polyunsaturated fatty acids, such as arachidonic acid (ARA), eicosapentaenoic acid (EPA), and docosahexaenoic acid (DHA), have well-documented importance in human health and nutrition. Sustainable production in robust host organisms that do not synthesize them naturally requires the coordinated expression of several heterologous desaturases and elongases. In the present study we show production of EPA in Saccharomyces cerevisiae using glucose as the sole carbon source through expression of five heterologous fatty acid desaturases and an elongase. Novel Δ5-desaturases from the ciliate protozoan Paramecium tetraurelia and from the microalgae Ostreococcus tauri and Ostreococcus lucimarinus were identified via a BLAST search, and their substrate preferences and desaturation efficiencies were assayed in a yeast strain producing the ω6 and ω3 fatty acid substrates for Δ5-desaturation. The Δ5-desaturase from P. tetraurelia was up-to-2-fold more efficient than the microalgal desaturases and was also more efficient than Δ5-desaturases from Mortierella alpina and Leishmania major. In vivo investigation of acyl carrier substrate specificities showed that the Δ5-desaturases from P. tetraurelia, O. lucimarinus, O. tauri, and M. alpina are promiscuous toward the acyl carrier substrate but prefer phospholipid-bound substrates. In contrast, the Δ5-desaturase from L. major showed no activity on phospholipid-bound substrate and thus appears to be an exclusively acyl coenzyme A-dependent desaturase.

Marked amplification and diversification of products of ras genes from rat brain, Rab GTPases, in the ciliates Tetrahymena thermophila and Paramecium tetraurelia.

Small GTPase Rab (products of ras genes from rat brain) is a widely conserved molecular switch among eukaryotes and regulates membrane trafficking pathways. It is generally considered that the number of Rab encoded in the genome correlates with multicellularity; however, we found that unicellular ciliates Tetrahymena thermophila (Tt) and Paramecium tetraurelia (Pt) possess many more Rab genes in their genome than the 64 HsRab genes in the human genome. We succeeded in isolating 86 cDNA clones of 88 TtRab genes in the Tetrahymena genome. By comparing the amino acid sequence of Rab in humans and the budding yeast Saccharomyces cerevisiae, 42 TtRab belonged to subfamilies functionally characterized and designated as conventional Rab, while the remaining 44 TtRab were considered to be species-specific. To examine the diversity of Rab in ciliates, we searched for Rab genes in the genome database of P. tetraurelia. Overall, 229 PtRab genes were found and categorized as 157 conventional and 72 species-specific PtRab, respectively. Among them, nine PtRab genes showed high homology to seven TtRab, suggesting the conservation of ciliate-specific Rab. These data suggested that the range of Rab is markedly amplified and diversified in ciliates, which may support the elaborate cellular structures and vigorous phagocytosis of those organisms.

Protein phosphatase 2B (PP2B, calcineurin) in Paramecium: partial characterization reveals that two members of the unusually large catalytic subunit family have distinct roles in calcium-dependent processes.

We characterized the calcineurin (CaN) gene family, including the subunits CaNA and CaNB, based upon sequence information obtained from the Paramecium genome project. Paramecium tetraurelia has seven subfamilies of the catalytic CaNA subunit and one subfamily of the regulatory CaNB subunit, with each subfamily having two members of considerable identity on the amino acid level (>or=55% between subfamilies, >or=94% within CaNA subfamilies, and full identity in the CaNB subfamily). Within CaNA subfamily members, the catalytic domain and the CaNB binding region are highly conserved and molecular modeling revealed a three-dimensional structure almost identical to a human ortholog. At 14 members, the size of the CaNA family is unprecedented, and we hypothesized that the different CaNA subfamily members were not strictly redundant and that at least some fulfill different roles in the cell. This was tested by selecting two phylogenetically distinct members of this large family for posttranscriptional silencing by RNA interference. The two targets resulted in differing effects in exocytosis, calcium dynamics, and backward swimming behavior that supported our hypothesis that the large, highly conserved CaNA family members are not strictly redundant and that at least two members have evolved diverse but overlapping functions. In sum, the occurrence of CaN in Paramecium spp., although disputed in the past, has been established on a molecular level. Its role in exocytosis and ciliary beat regulation in a protozoan, as well as in more complex organisms, suggests that these roles for CaN were acquired early in the evolution of this protein family.

Distinct RNA-dependent RNA polymerases are required for RNAi triggered by double-stranded RNA versus truncated transgenes in Paramecium tetraurelia.

In many eukaryotes, RNA-dependent RNA polymerases (RdRPs) play key roles in the RNAi pathway. They have been implicated in the recognition and processing of aberrant transcripts triggering the process, and in amplification of the silencing response. We have tested the functions of RdRP genes from the ciliate Paramecium tetraurelia in experimentally induced and endogenous mechanisms of gene silencing. In this organism, RNAi can be triggered either by high-copy, truncated transgenes or by directly feeding cells with double-stranded RNA (dsRNA). Surprisingly, dsRNA-induced silencing depends on the putatively functional RDR1 and RDR2 genes, which are required for the accumulation of both primary siRNAs and a distinct class of small RNAs suggestive of secondary siRNAs. In contrast, a third gene with a highly divergent catalytic domain, RDR3, is required for siRNA accumulation when RNAi is triggered by truncated transgenes. Our data further implicate RDR3 in the accumulation of previously described endogenous siRNAs and in the regulation of the surface antigen gene family. While only one of these genes is normally expressed in any clonal cell line, the knockdown of RDR3 leads to co-expression of multiple antigens. These results provide evidence for a functional specialization of Paramecium RdRP genes in distinct RNAi pathways operating during vegetative growth.

Subtraction by addition: domesticated transposases in programmed DNA elimination.

The ciliate Paramecium tetraurelia must eliminate approximately 60,000 short sequences from its genome to generate uninterrupted coding sequences in its somatic macronucleus. In this issue of Genes & Development, Baudry and colleagues (pp. 2478-2483) identify the protein that excises these noncoding sequences: a domesticated piggyBac transposase that has been adapted to remove what are likely the remnants of transposon insertions. This new study reveals how addition of a transposase to small RNA-directed silencing machinery can guide major genome reorganization.

PiggyMac, a domesticated piggyBac transposase involved in programmed genome rearrangements in the ciliate Paramecium tetraurelia.

Programmed genome rearrangements drive functional gene assembly in ciliates during the development of the somatic macronucleus. The elimination of germline sequences is directed by noncoding RNAs and is initiated by DNA double-strand breaks, but the enzymes responsible for DNA cleavage have not been identified. We show here that PiggyMac (Pgm), a domesticated piggyBac transposase, is required for these rearrangements in Paramecium tetraurelia. A GFP-Pgm fusion localizes in developing macronuclei, where rearrangements take place, and RNAi-mediated silencing of PGM abolishes DNA cleavage. This is the first in vivo evidence suggesting an essential endonucleolytic function of a domesticated piggyBac transposase.

A kinome of 2600 in the ciliate Paramecium tetraurelia.

Protein kinases play a crucial role in the regulation of cellular processes. Most eukaryotes reserve about 2.5% of their genes for protein kinases. We analysed the genome of the single-celled ciliate Paramecium tetraurelia and identified 2606 kinases, about 6.6% of its genes, representing the largest kinome to date. A gene tree combined with human kinases revealed a massive expansion of the calcium calmodulin regulated subfamily, underlining the importance of calcium in the physiology of P. tetraurelia. The kinases are embedded in only 40 domain architectures, contrasting 134 in human. This might indicate different mechanisms to achieve target specificity.

Aziridide-based inhibitors of cathepsin L: synthesis, inhibition activity, and docking studies.

A comprehensive screening of N-acylated aziridine (aziridide) based cysteine protease inhibitors containing either Boc-Leu-Caa (Caa=cyclic amino acid), Boc-Gly-Caa, or Boc-Phe-Ala attached to the aziridine nitrogen atom revealed Boc-(S)-Leu-(S)-Azy-(S,S)-Azi(OBn)(2) (18 a) as a highly potent cathepsin L (CL) inhibitor (K(i)=13 nM) (Azy=aziridine-2-carboxylate, Azi=aziridine-2,3-dicarboxylate). Docking studies, which also accounted for the unusual bonding situations (the flexibility and hybridization of the aziridides) predict that the inhibitor adopts a Y shape and spans across the entire active site cleft, binding into both the nonprimed and primed sites of CL.

Seventeen a-subunit isoforms of paramecium V-ATPase provide high specialization in localization and function.

In the Paramecium tetraurelia genome, 17 genes encoding the 100-kDa-subunit (a-subunit) of the vacuolar-proton-ATPase were identified, representing by far the largest number of a-subunit genes encountered in any organism investigated so far. They group into nine clusters, eight pairs with >82% amino acid identity and one single gene. Green fluorescent protein-tagging of representatives of the nine clusters revealed highly specific targeting to at least seven different compartments, among them dense core secretory vesicles (trichocysts), the contractile vacuole complex, and phagosomes. RNA interference for two pairs confirmed their functional specialization in their target compartments: silencing of the trichocyst-specific form affected this secretory pathway, whereas silencing of the contractile vacuole complex-specific form altered organelle structure and functioning. The construction of chimeras between selected a-subunits surprisingly revealed the targeting signal to be located in the C terminus of the protein, in contrast with the N-terminal targeting signal of the a-subunit in yeast. Interestingly, some chimeras provoked deleterious effects, locally in their target compartment, or remotely, in the compartment whose specific a-subunit N terminus was used in the chimera.

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.

Ca(2+)-binding proteins of cilia and infraciliary lattice of Paramecium tetraurelia: their phosphorylation by purified endogenous Ca(2+)-dependent protein kinases.

We purified two small, acidic calcium-binding proteins (Paramecium Ca(2+)-binding proteins, PCBP-25alpha and PCBP-25beta) from Paramecium tetraurelia by Ca(2+)-dependent chromatography on phenyl-Sepharose and by anion-exchange chromatography. The proteins were immunologically distinct. Monoclonal antibodies against PCBP-25beta did not react with PCBP-25alpha, and antibodies against centrin from Chlamydomonas reacted with PCBP-25alpha but not with PCBP-25beta. Like the centrins described previously, both PCBPs were associated with the infraciliary lattice (ICL), a fibrillar cytoskeletal element in Paramecium. Both were also present in isolated cilia, from which they could be released (with dynein) by a high-salt wash, and both PCBPs cosedimented with dynein in a sucrose gradient. PCBP-25beta was especially prominent in cilia and in the deciliation supernatant, a soluble fraction released during the process of deciliation. The results of immunoreactivity and localization experiments suggest that PCBP-25alpha is a Paramecium centrin and that PCBP-25beta is a distinct Ca(2+)-binding protein that confers Ca(2+) sensitivity on some component of the cilium, ciliary basal body or ICL. We characterized these proteins and Paramecium calmodulin as substrates for two Ca(2+)-dependent protein kinases purified from Paramecium. PCBP-25alpha and calmodulin were in vitro substrates for one of the two Ca(2+)-dependent protein kinases (CaPK-2), but only PCBP-25alpha was phosphorylated by CaPK-1. These results raise the possibility that the biological activities of PCBP-25alpha and calmodulin are regulated by phosphorylation.

A novel member of the cyclin-dependent kinase family in Paramecium tetraurelia.

Passage through the cell cycle in eukaryotes requires the successive activation of different cyclin-dependent protein kinases. Here, we describe the identification and characterization of a novel class of cyclin-dependent protein kinase, termed Cdk2, in the ciliate Paramecium tetraurelia. It is 301 amino acids long, 7 amino acids shorter than Cdk1, the CDK that is associated with macronuclear DNA synthesis. All the catalytic domains typical of protein kinases can be located within the sequence and putative regulatory phosphorylation sites equivalent to Thr14, Tyr15, and Thr161 in human CDK1 are also conserved. The 'PSTAIRE' region characteristic of most CDKs is perfectly conserved. Cdk2 shares only 48% homology to Cdk1 at the amino acid level, suggesting that the evolutionary separation of Cdk1 and Cdk2 is ancient, and implying that they have different roles in cell cycle regulation. Like Cdk1, Cdk2 does not bind to yeast p13suc1, even though it has better conservation of p13suc1 binding sites than Cdk1 does. The Cdk2 protein level is relatively constant throughout the vegetative cell cycle. Cdk2 exhibits kinase activity towards bovine histone H1 in vitro with the maximal level late in the cell cycle, suggesting it may be involved in the regulation of cytokinesis. Our results further support the view that an analogue of the cyclin-dependent kinase cell cycle regulatory system like that of yeast and higher eukaryotic cells operates in Paramecium and that a family of cyclin-dependent kinases may control different aspects of the Paramecium cell cycle.

Expression of mutated Paramecium telomerase RNAs in vivo leads to templating errors that resemble those made by retroviral reverse transcriptase.

Telomeric DNA consists of short, tandemly repeated sequences at the ends of chromosomes. Telomeric DNA in the ciliate Paramecium tetraurelia is synthesized by an error-prone telomerase with an RNA template specific for GGGGTT repeats. We have previously shown that misincorporation of TTP residues at the telomerase RNA templating nucleotide C52 accounts for the 30% GGGTTT repeats randomly distributed in wild-type telomeres. To more completely characterize variable repeat synthesis in P. tetraurelia, telomerase RNA genes mutated at C52 (A, U, and G) were expressed in vivo. De novo telomeric repeats from transformants indicate that the predominant TTP misincorporation error seen in the wild-type telomerase is dependent on the presence of a C residue at template position 52. Paradoxically, the effects of various other telomerase RNA template and alignment region mutations on de novo telomeres include significant changes in fidelity, as well as the synthesis of aberrant, 5-nucleotide telomeric repeats. The occurrence of deletion errors and the altered fidelity of mutated P. tetraurelia telomerase, in conjunction with misincorporation by the wild-type enzyme, suggest that the telomerase RNA template domain may be analogous to homopolymeric mutational hot spots that lead to similar errors by the human immunodeficiency virus proofreading-deficient reverse transcriptase.

Glyceraldehyde-3-phosphate dehydrogenase from Tetrahymena pyriformis: enzyme purification and characterization of a gapC gene with primitive eukaryotic features.

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH, EC.1.2.1.12) was purified to electrophoretic homogeneity from an amicronucleated strain of the ciliate Tetrahymena pyriformis using a three-step procedure. The native enzyme is an homotetramer of 145 kDa exhibiting absolute specificity for NAD. In its catalytic properties it is similar to other glycolytic GAPDHs. Chromatofocusing analysis showed the presence of only one basic GAPDH isoform with an isoelectric point of 8.8. Western blots using a monospecific polyclonal antibody raised against the T. pyriformis GAPDH showed a single 36-kDa band corresponding to the enzyme subunit in the cytosolic protein fraction of this strain and the closely related species, both from the class Oligohymenophorea, Paramecium tetraurelia. No bands were immunodetected in the ciliate Colpoda inflata (class Colpodea) and in the diverse eukaryotes and eubacteria tested. A 0.5-kb DNA fragment which corresponds to an internal region of a gapC gene was generated by polymerase chain reaction using cDNA of T. pyriformis as template. This gene codes for a basic GAPDH protein with eukaryotic-diplomonad signatures and exhibits a codon usage biased in the manner typical for T. pyriformis genes. Southern blots performed both under homologous and heterologous conditions using this amplified cDNA fragment as a probe, indicated that it should be the only gapC gene present in the macronuclear genome of this ciliate, its expression being confirmed by Northern blot analysis. These results are discussed in connection with the peculiar genomic organization of ciliates and in the context of protist evolution.