Ciliate Genes Research Articles

Parallel complexity of signed graphs for gene assembly in ciliates

We consider a graph-based model for the study of parallelism in ciliate gene assembly, where a signed graph is associated to each micronuclear gene and the gene assembly is modeled as a graph rewriting process. A natural measure of complexity for gene assembly counts the minimal number of parallel steps needed to reduce the associated signed graph. We investigate the complexity of several classes of the graphs, so far found graphs of parallel complexity up to six. The general problem of whether there exists a finite upper bound for the graph parallel complexity still remains open.

Cloning, site-directed mutagenesis and expression of cathepsin L-like cysteine protease from Uronema marinum (Ciliophora: Scuticociliatida)

A cysteine protease gene (ScCtL) homologous to the cathepsin L genes was isolated from a cDNA library of the scuticociliate parasite (Uronema marinum). To express the ScCtL recombinant protein in heterologous system, 17 codons were redesigned to conform to the standard eukaryotic genetic code using PCR-based site-directed mutagenesis. The synthetic U. marinum procathepsin L (proScCtL) was expressed at high levels in E. coli BL21 (DE3) with pGEX-4T-1 vector, and successfully refolded and purified into a functional and enzymatically active form. The optimal pH for protease activity was found to be 4.5. Like any typical cysteine protease, the enzyme was inhibited by E-64 and leupeptin. A dot-blot immunoassay was conducted in an attempt to determine the reaction abilities and sensitivity of the anti-proScCtL polyclonal antibody to the cytosol and to the membrane fraction from the scuticociliate. Our results suggest that the biochemical characteristics of the recombinant ciliate proScCtL protein are similar to that of the cathepsin L-like cysteine protease, and that the PCR-based site-direct mutated ciliate gene was successfully expressed in a biochemically active form.

Horizontal gene transfer from Bacteria to rumen Ciliates indicates adaptation to their anaerobic, carbohydrates-rich environment

BackgroundThe horizontal transfer of expressed genes from Bacteria into Ciliates which live in close contact with each other in the rumen (the foregut of ruminants) was studied using ciliate Expressed Sequence Tags (ESTs). More than 4000 ESTs were sequenced from representatives of the two major groups of rumen Cilates: the order Entodiniomorphida (Entodinium simplex, Entodinium caudatum, Eudiplodinium maggii, Metadinium medium, Diploplastron affine, Polyplastron multivesiculatum and Epidinium ecaudatum) and the order Vestibuliferida, previously called Holotricha (Isotricha prostoma, Isotricha intestinalis and Dasytricha ruminantium).ResultsA comparison of the sequences with the completely sequenced genomes of Eukaryotes and Prokaryotes, followed by large-scale construction and analysis of phylogenies, identified 148 ciliate genes that specifically cluster with genes from the Bacteria and Archaea. The phylogenetic clustering with bacterial genes, coupled with the absence of close relatives of these genes in the Ciliate Tetrahymena thermophila, indicates that they have been acquired via Horizontal Gene Transfer (HGT) after the colonization of the gut by the rumen Ciliates.ConclusionAmong the HGT candidates, we found an over-representation (>75%) of genes involved in metabolism, specifically in the catabolism of complex carbohydrates, a rich food source in the rumen. We propose that the acquisition of these genes has greatly facilitated the Ciliates' colonization of the rumen providing evidence for the role of HGT in the adaptation to new niches.

Developmentally and environmentally regulated expression of gamone 1: the trigger molecule for sexual reproduction inBlepharisma japonicum

Sexual reproduction (conjugation) in protozoan ciliates is induced by specific cell-cell interactions between cells of complementary mating types. The ancestral ciliate Blepharisma japonicum has two mating types, I and II. The substances that act as signaling molecules in this extracellular interaction for conjugation are called gamones. The glycoprotein gamone 1, produced by mating type I cells, is a key factor that triggers this interaction. We have previously isolated gamone 1 and determined its complete amino acid sequence. To elucidate the mechanism of initiation of conjugation in ciliates, we investigated the transcription of the gamone 1 gene and found that it is controlled by various internal and external factors. The gamone 1 gene transcript appeared specifically when sexually mature mating type I cells were starved. It was not detected in immature cells, mating type II cells or proliferating cells. The level of transcription was markedly increased in type I cells when they were stimulated with gamone 2, which is secreted by type II cells. This is the first report that the transcription of gamone genes in ciliates is strictly regulated by developmental and environmental factors. This study suggests that the onset of transcription of gamone 1 is linked to the switching mechanism that converts mitotically proliferating cells to differentiated preconjugants, the mechanism of differentiation from immature to mature cells in clonal development, and the mechanism that ensures mating type-specific gene silencing.

Template-guided DNA recombination

The family of stichotrichous ciliates have received a great deal of study due to the presence of scrambled genes in their genomes. The mechanism by which these genes are descrambled is of interest both as a biological process and as a model of natural computation. Several formal models of this process have been proposed, the most recent of which involves the recombination of DNA strands based on template guides. We generalize this template-guided DNA recombination model proposed by Prescott, Ehrenfeucht and Rozenberg to an operation on strings and languages. We then proceed to investigate the properties of this operation with the intention of viewing ciliate gene descrambling as a computational process.

Gene Unscrambler for detangling scrambled genes in ciliates.

Scrambled genes are surprisingly common in some species of ciliates. Until now there was no software available to analyze automatically these genes. We present here a program that can automatically align the macronuclear and micronuclear forms of a gene, outputting the location of the macronuclear destined segments and pointer sequences. A web version of the program is available free of charge and can be accessed at http://oxytricha.princeton.edu/GeneUnscrambler.htm

Myb Genes in Ciliates: A Common Origin with the myb Protooncogene?

In animals, the protooncogene myb family is characterized by a DNA-binding domain (so-called MYB domain), which consists of 3 imperfect tandem repeats of a helix-turn-helix motif. Homologous genes have been characterized in plants and also in Dictyostelium discoideum. However, in plants, the myb family is more diverse and displays 2 types of MYB domains: the animal-like 3 repeats (MYB-3R) and the 2 repeats (MYB-2R) domains. The question is therefore raised as to the putative existence of genes with MYB-3R and/or MYB-2R domains in their last common unicellular ancestor. Here, we present evidence that in ciliates like in plants, both types of domain exist. A gene having a MYB-3R domain has been identified in the oxytrichid Sterkiella histriomuscorum and a gene having a MYB-2R domain has been identified in the euplotid Euplotes aediculatus. Both genes are expressed during the vegetative growth of the cells. A conserved intron exists in the gene of Sterkiella and phylogenetical analyses show that the 2 ciliate genes belong to the myb protooncogene family as deeply split lineages. This is the first report of a myb homolog in a ciliated protist, thus, confirming its origin in strict unicellular eukaryotes.

Template-guided recombination for IES elimination and unscrambling of genes in stichotrichous ciliates

The micronuclear versions of genes in stichotrichous ciliates are interrupted by multiple, short, non-coding DNA segments called internal eliminated segments, or IESs. IESs divide a gene into macronuclear destined segments, or MDSs. In some micronuclear genes MDSs are in a scrambled disorder. During development of a micronucleus into a macronucleus after cell mating the IESs are excised from micronuclear genes and the MDSs are spliced in the sequentially correct order. Pairs of short repeat sequences in the ends of MDSs undergo homologous recombination to excise IESs and splice MDSs. However, the repeat sequences are too short to guide unambiguously their own alignment in preparation for recombination. Based on experiments by others on the distantly related ciliate, Paramecium, we propose a molecular model of template-guided recombination to explain the excision of the 100 000–150 000 IESs and splicing of MDSs, including unscrambling, in the genome of stichotrichous ciliates. The model solves the problem of correct pairing of pointers, precisely identifies MDS–IES junctions, and provides for irreversible recombination.

The use of synthetic genes for the expression of ciliate proteins in heterologous systems

The common fish parasite, Ichthyophthirius multifiliis, expresses abundant glycosylated phosphatidylinositol (GPI)-anchored membrane proteins known as immobilization antigens, or i-antigens. These proteins are targets of the host immune response, and have been identified as potential candidates for recombinant subunit vaccine development. Nevertheless, because Ichthyophthirius utilizes a non-standard genetic code, expression of the corresponding gene products, either as subunit antigens in conventional protein expression systems, or as vector-encoded antigens in the case of DNA vaccines, is far from straightforward. To overcome this problem, we utilized ‘assembly polymerase chain reaction’ to manufacture synthetic versions of two genes (designated IAG52A[G5/CC] and IAG52B[G5/CC]) encoding ∼52/55 kDa i-antigens from parasite strain G5. This approach made it possible to eliminate unwanted stop codons and substitute the preferred codon usage of channel catfish for the native sequences of the genes. To determine whether the synthetic alleles could be expressed in cells that use the standard genetic code, we introduced IAG52A[G5/CC] into a variety of heterologous cell types and tested for expression either by immunofluorescence light microscopy or Western blotting. When cloned downstream of appropriate promoters, IAG52A[G5/CC] was expressed in Escherichia coli, mammalian COS-7 cells, and channel catfish where it elicited antigen-specific immune responses. Interestingly, the localization pattern of the corresponding gene product in COS-7 cells indicated that while the protein was correctly folded, it was not present on the cell membrane, suggesting that the signal peptides required for GPI-anchor addition differ in ciliate and mammalian systems. Construction of synthetic alleles should have practical utility in the development of vaccines against Ichthyophthirius, and at the same time, provide a general method for the expression of ciliate genes in heterologous systems.

A homologue of CROC-1 in a ciliated protist (Sterkiella histriomuscorum) testifies to the ancient origin of the ubiquitin-conjugating enzyme variant family.

Resting cysts of Sterkiella histriomuscorum (Ciliophora, Oxytrichidae) have been shown to contain messenger RNA, one of which codes for a protein significantly similar to CROC-1. CROC-1 is a human regulatory protein capable of transactivating the promoter of c-fos and belongs to a newly characterized family of ubiquitin-conjugating enzyme (E2) variants (UEV). We have determined the corresponding macronuclear gene sequence, which is the first protistan UEV sequence available. The phylogenetic analysis indicates the deep separation and solid clustering of all the UEV sequences within the E2 tree showing the ancient origin of these regulatory genes and their high structural conservation during evolution. Furthermore, overexpression of the ciliate UEV is able to rescue the Saccharomyces cerevisiae mms2 null mutant from killing by DNA damaging agents, implying that the UEV family proteins are functionally conserved. In S. histriomuscorum, expression of UEV is correlated with the growth of the cells as transcripts are present in excysting and vegetative cells but are rapidly down-regulated during starvation. These data support the high conservation of the UEV family in eukaryotes, and a regulatory role of the gene is discussed in relation to known functions of UEVs. This analysis may promote the search for homologues of other regulatory genes (metazoan regulators of differentiation) in ciliates.

Evolution of IESs and scrambling in the actin I gene in hypotrichous ciliates.

Germ-line (micronuclear) genes in hypotrichous ciliates are interrupted by numerous, short, noncoding, AT-rich segments called internal eliminated segments, or IESs. IESs divide a gene into macronuclear destined segments, or MDSs. IESs are excised from micronuclear genes, and the MDSs are spliced when a micronuclear genome is processed into a macronuclear genome after cell mating. In the micronuclear version of the actin I gene intramolecular recombination between IESs during evolution has put MDSs into a scrambled disorder in some but not all hypotrichs. Studies using rDNA sequences to define phylogenetic relationships among eight hypotrichs suggests that evolution of the micronuclear actin I gene proceeds by successive addition of IESs in earlier diverging species, without MDS scrambling. Continued addition of IESs and recombination among IESs in later diverging species produced actin I genes with scrambled MDSs. Subsequent to MDS scrambling, additional IESs were inserted into the more recently evolved species. Thus, IES insertions and gene scrambling occur in a progressive manner during species evolution to produce micronuclear actin I genes of increasing structural complexity.

Messenger RNA in Dormant Cells of Sterkiella histriomuscorum (Oxytrichidae): Identification of Putative Regulatory Gene Transcripts

In the absence of food, the oxytrichid Sterkiella histriomuscorum, like many ciliates, enters into dormancy and transforms into a round and walled encysted cell. When transferred back into a feeding medium, the cyst re-transforms into a vegetative cell in a few hours. This encystment-excystment pathway, which is common to many free-living and parasitic protists, is still poorly understood at the molecular level. In order to identify potential dormant transcripts in the cysts of Sterkiella, we have constructed cDNA libraries from mature cysts. Transcripts have been isolated confirming the presence of a mRNA pool in the dormant cells. The sequence analysis of two cDNA indicates open reading frames which show significant similarities to known proteins involved in mechanisms of regulation: 1) nifR3, an element of the nitrogen regulatory system in bacteria and 2) CROC-1, a newly identified human transcription factor. The two corresponding macronuclear genes represent the first putative regulatory genes isolated in ciliates. From a differential screening of the cDNA library against vegetative cDNA, one cyst-specific (and very abundant) transcript has been isolated but the product has not yet been identified. The possible involvment of these new ciliate genes in the excystment process is discussed.

On the possibility of constructive neutral evolution.

The neutral theory often is presented as a theory of "noise" or silent changes at an isolated "molecular level," relevant to marking the steady pace of divergence, but not to the origin of biological structure, function, or complexity. Nevertheless, precisely these issues can be addressed in neutral models, such as those elaborated here with regard to scrambled ciliate genes, gRNA-mediated RNA editing, the transition from self-splicing to spliceosomal splicing, and the retention of duplicate genes. All of these are instances of a more general scheme of "constructive neutral evolution" that invokes biased variation, epistatic interactions, and excess capacities to account for a complex series of steps giving rise to novel structures or operations. The directional and constructive outcomes of these models are due not to neutral allele fixations per se, but to these other factors. Neutral models of this type may help to clarify the poorly understood role of nonselective factors in evolutionary innovation and directionality.

The evolutionary scrambling and developmental unscrambling of germline genes in hypotrichous ciliates.

Genes in the germline (micronuclear) genome of hypotrichous ciliates are interrupted by multiple, short, non-coding, AT-rich sequences called internal eliminated segments, or IESs. During conversion of a micronucleus to a somatic nucleus (macronucleus) after cell mating, all IESs are excised from the germline genes and the gene segments, called macronuclear-destined segments, or MDSs, are spliced. Excision of the approximately 150 000 IESs from a haploid germline genome in Oxytricha nova requires approximately 150 000 recombinant events. In three of 10 genes the MDSs are scrambled. During macronuclear development the MDSs are unscrambled, possibly by folding of the DNA to allow MDSs to ligate in the correct order. The nine MDSs in the actin I gene of O.nova are scrambled in the random order, 3-4-6-5-7-9-2-1-8, and MDS 2 is inverted. The 14 MDSs in the alphaTP gene of O.nova and Stylonychia mytilus are scrambled in the non-random order, 1-3-5-7-9-11-2-4-6-8-10-12-13-14. The 45 MDSs in the DNA pol alpha gene are non-randomly scrambled into an odd/even series, with an inversion of one-third of the gene. Additional IESs have been inserted into these three genes during evolution of Oxytricha trifallax, slightly modifying scrambling patterns. The non-random scrambled patterns in the alphaTP and DNA pol alpha genes are explained by multiple, simultaneous IES insertions. The randomly scrambled pattern in the actin I gene may arise from an initially non-randomly scrambled pattern by recombination among multiple IESs. Alternatively, IESs inserted sporadically (individually) in a non-scrambled configuration might subsequently recombine, converting a non-scrambled gene into a randomly scrambled one. IESs shift along a DNA molecule, most likely as a result of mutations at MDS/IES junctions. Shifting of IESs has the effect of 'transferring' nucleotides from one MDS to another, but does not change the overall sequence of nucleotides in the combined MDSs. In addition to shifting in position, IESs accumulate mutations at a high rate and increase and decrease in length within a species and during speciation. The phenomena of IESs and of MDS scrambling represent remarkable flexibility of the hypotrich genome, possibly reflecting a process of MDS shuffling that facilitates the evolution of genes.

Internal eliminated sequences interrupting the Oxytricha 81 locus: allelic divergence, conservation, conversions, and possible transposon origins.

Internal eliminated sequences (IESs) often interrupt ciliate genes in the silent germline nucleus but are exactly excised and eliminated from the developing somatic nucleus from which genes are then expressed. Some long IESs are transposons, supporting the hypothesis that short IESs are ancient transposon relics. In light of that hypothesis and to explore the evolutionary history of a collection of IESs, we have compared various alleles of a particular locus (the 81 locus) of the ciliated protozoa Oxytricha trifallax and O. fallax. Three short IESs that interrupt two genes of the locus are found in alleles from both species, and thus must be relatively ancient, consistent with the hypothesis that short IESs are transposon relics. In contrast, TBE1 transposon interruptions of the locus are allele-specific and probably the results of recent transpositions. These IESs (and the TBE1s) are precisely excised from the DNA of the developing somatic macronucleus. Each IES interrupts a highly conserved sequence. A few nucleotides at the ends of each IES are also conserved, suggesting that they interact critically with IES excision machinery. However, most IES nucleotide positions have evolved at high rates, showing little or no selective constraint for function. Nonetheless, the length of each IES has been maintained (+/- 3 bp). While one IES is approximately 33 bp long, three other IESs have very similar sizes, approximately 70 bp long. Two IESs are surrounded by direct repeats of the sequence TTCTT. No other sequence similarities were found between any of the four IESs. However, the ends of one IES do match the inverted terminal repeat consensus sequence of the "TA" IESs of Paramecium. Three O. trifallax alleles appear to have been recipients in recent conversion events that could have been provoked by double-strand breaks associated with IES ends subsequent to IES transposition. Our findings support the hypothesis that short IESs evolved from ancient transposons that have lost most of their sequences, except those necessary for precise excision during macronuclear development.

Sequence, codon usage and cysteine periodicity of the SerH1 gene and in the encoded surface protein of Tetrahymena thermophila

The temperature-regulated SerH1 gene coding for an immunodominant surface glycoprotein (i-Ag H1) of Tetrahymena thermophila has been sequenced. The gene is reproducibly rearranged during macronuclear development and steady state mRNA levels are present at < 36° C. The deduced i-Ag H1 amino acid (aa) sequence is rich in Ser, Thr and Cys, and contains three periods each consisting of 85 aa punctuated by eight Cys with the general formula, CX 6CX 17CX 2CX 18CX 2CX 11CX 2CX 19 (where X = any aa). Such Cys periodicity is common to ciliate i-Ag. Codon usage in Tt, Paramecium primaurelia and P. tetraurelia i-Ag encoding genes is similar, with approx. 80% A + T in the 3′ position which is in marked contrast to the approx. 54% 3′ A + T in other ciliate genes.

Expression of a ciliate gene in Escherichia coli using a suppressor tRNA to read the UAA and UAG glutamine codons

Most ciliates use a particular genetic code where the standard stop codons UAA and UAG encode glutamine. Ciliate genes cannot therefore be expressed in heterologous systems such as Escherichia coli. To overcome this problem, we worked out a system of inducible suppression to permit efficient readthrough of UAAs and UAGs: a strong UAA tRNA suppressor that inserts glutamic acid was cloned downstream from a tac promoter whose efficiency was reduced by a transcription terminator. This system proved to be operational (1) to suppress UAG mutations by wobble pairing in an E. coli lacI-lacZ gene fusion and (2) to read through at least eight UAA glutamine codons in a Paramecium alpha-tubulin gene, as detected by Western blotting and colony hybridization. This work opens the way for cloning Ciliate genes from expression libraries and for expressing particular sequences without extended in vitro mutagenesis. A similar approach can be envisaged for expression of genes from Mycoplasma, mitochondria or other genomes that use non-standard genetic codes.

Characterizaton of the promoter region of Tetrahymena genes

The regions between adjacent histone H3 and H4 genes, as well as portions of the genes, from 22 species of Tetrahymena have been amplified using the polymerase chain reaction and sequenced. Both histone genes are transcribed divergently with initiation occurring within the intergenic region, thus 2 sets of 22 homologous Tetrahymena promoters can be compared. A sequence comparison of these regions reveals a single putative promoter element, with a consensus sequence TATCCAATTCARA, present in front of each gene. This sequence contains a 'CCAAT' box, which also occurs at 8 locations preceding other ciliate genes. No other putative promoter sequences are found in front of these sets of histone genes. Sequences searched for include 'TATA' boxes, 'GC' boxes and other sequences suggested as putative promoter elements for ciliate genes. The coding strand immediately preceding ciliate genes is very high in A content and the consensus sequence at the site of protein synthesis is AAAATGG.

Codon Usage in Tetrahymena and Other Ciliates

Codon usage in ciliates was examined by analyzing the coding regions of 22 ciliate genes corresponding to a total of 26,142 nucleotides (8,714 codons). It was found that Tetrahymena, Paramecium and the hypotrichs (Oxytricha and Stylonychia) differed in which synonymous codons were used most frequently by their genes. In fact, the codon choices in highly expressed Tetrahymena genes were more similar to those of yeast genes than those of Paramecium genes. The ciliates do not appear to have unusually strong biases in codon usage frequency when compared to other protists such as yeast. The analysis of the Tetrahymena genes indicated that genes which are highly expressed during normal cell growth have a stronger bias towards using the "preferred" codons than those expressed at lower levels during growth or for brief periods during processes such as conjugation. This conforms to what is found in other protists.

Multiple introns in a conjugation-specific gene from Tetrahymena thermophila.

Multiple introns have been found in a gene from a ciliated protozoan. This Tetrahymena thermophila gene (cnjB) is large (7.5 kb mRNA) and active only during conjugation, the organism's sexual cycle. Six introns ranging in size from 62 bp to 676 bp were found when we sequenced a 3.1 kb segment of the cnjB gene together with its corresponding cDNA. We estimate, by extrapolation of our current data, a total of approximately 30 introns in this gene with a total gene size (introns plus exons) of 15 kb or more. The number of introns is surprising given the scarcity of introns in ciliate genes examined to date. Our findings constitute the first example of multiple introns in a ciliate gene. Having the sequence of several introns has allowed us to construct consensus sequences for T. thermophila mRNA introns. The 5' and 3' intron junctions resemble those of general nuclear mRNA (GT/AG rule is followed) but differences are seen. In particular, stretches of 10 or more adenines and thymines are found adjacent to the conserved GT and AGs at the junctions. Unusual aspects of the coding region of this gene are discussed.