Analysis of expressed sequence tags (ESTs) in the ciliated protozoan Tetrahymena thermophila

The transcriptome of the rumen ciliate Entodinium caudatum reveals some of its metabolic features

Background

Rumen ciliates play important roles in rumen function by digesting and fermenting feed and shaping the rumen microbiome. However, they remain poorly understood due to the lack of definitive direct evidence without influence by prokaryotes (including symbionts) in co-cultures or the rumen. In this study, we used RNA-Seq to characterize the transcriptome of Entodinium caudatum, the most predominant and representative rumen ciliate species.

Results

Of a large number of transcripts, > 12,000 were annotated to the curated genes in the NR, UniProt, and GO databases. Numerous CAZymes (including lysozyme and chitinase) and peptidases were represented in the transcriptome. This study revealed the ability of E. caudatum to depolymerize starch, hemicellulose, pectin, and the polysaccharides of the bacterial and fungal cell wall, and to degrade proteins. Many signaling pathways, including the ones that have been shown to function in E. caudatum, were represented by many transcripts. The transcriptome also revealed the expression of the genes involved in symbiosis, detoxification of reactive oxygen species, and the electron-transport chain. Overall, the transcriptomic evidence is consistent with some of the previous premises about E. caudatum. However, the identification of specific genes, such as those encoding lysozyme, peptidases, and other enzymes unique to rumen ciliates might be targeted to develop specific and effective inhibitors to improve nitrogen utilization efficiency by controlling the activity and growth of rumen ciliates. The transcriptomic data will also help the assembly and annotation in future genomic sequencing of E. caudatum.

Conclusion

As the first transcriptome of a single species of rumen ciliates ever sequenced, it provides direct evidence for the substrate spectrum, fermentation pathways, ability to respond to various biotic and abiotic stimuli, and other physiological and ecological features of E. caudatum. The presence and expression of the genes involved in the lysis and degradation of microbial cells highlight the dependence of E. caudatum on engulfment of other rumen microbes for its survival and growth. These genes may be explored in future research to develop targeted control of Entodinium species in the rumen. The transcriptome can also facilitate future genomic studies of E. caudatum and other related rumen ciliates.

Pharmacological characterization of N-methyl-d-aspartic acid (NMDA)-like receptors in the single-celled organism Paramecium primaurelia

Paramecium primaurelia is a unicellular eukaryote that moves in freshwater by ciliary beating and responds to environmental stimuli by altering motile behaviour. The movements of the cilia are controlled by the electrical changes of the cell membrane: when the intraciliary Ca2+ concentration associated with plasma membrane depolarization increases, the ciliary beating reverses its direction, and consequently the swimming direction changes. The ciliary reversal duration is correlated with the amount of Ca2+ influx. Here we evaluated the effects due to the activation or blockade of NMDA receptors on swimming behaviour in Paramecium. Paramecia normally swim forward drawing almost linear tracks. We observed that the simultaneous administration of NMDA and glycine induced a partial ciliary reversal (PaCR) leading to a continuous spiral-like swim. Furthermore, the duration of continuous ciliary reversal (CCR), triggered by high external KCl concentrations, was longer in NMDA/glycine treated cells. NMDA action required the presence of Ca2+, as the normal forward swimming was restored when the ion was omitted from the extracellular milieu. The PaCR and the enhancement of CCR duration significantly decreased when the antagonists of the glutamate site D-AP5 or CGS19755, the NMDA channel blocker MK-801, or the glycine site antagonist DCKA were added. The action of NMDA/glycine was also abolished by Zn2+ or ifenprodil, the GluN2A and the GluN2B NMDA-containing subunit blockers, respectively. Searches of the Paramecium genome database currently available indicate that the NMDA-like receptor with ligand binding characteristics of an NMDA receptor-like complex, purified from rat brain synaptic membranes and found in some metazoan genome, is also present in Paramecium. These results provide evidence that functional NMDA receptors similar to those typical of mammalian neuronal cells are present in the single-celled organism Paramecium and thus suggest that the glutamatergic NMDA system is a phylogenetically old behaviour-controlling mechanism.

Identification and characterization of polyubiquitin gene from cDNA library of aspergillus fumigatus

Aspergillus fumigatus (Afu) causes allergic and invasive forms of diseases in humans. In order to identify genes relevant for pathogenesis, a total of 235 cDNA clones were randomly selected and sequenced from cDNA library of Afu. One of the partially sequenced cDNA clones was homologous to polyubiquitin. Sequencing of the complete cDNA clone showed an open reading frame of 912 bases. Comparison with genomic sequence of Afu using BlastN program, revealed that polyubiquitin gene comprises of 992 bases and contains one intron of 80 bases. The recombinant expression of fusion protein showed an approximately molecular weight of 43-kDa on SDS-PAGE. The translation product of the cDNA sequence showed four tandem repeats of 76 amino acid residues in a single polyubiquitin protein and showed 100% identity with polyubiquitin protein sequences of S. cerevisiae, N. crassa, C. albicans, S. pombe, and M. grisae. Polyubiquitin gene is known to play important role in a variety of cellular processes and recently have been implicated in fungal pathogenesis. Identification of polyubiquitin gene of Afu has opened up scope to study its role in understanding Aspergillus biology and pathogenesis.

Toxicogenomic investigation of Tetrahymena thermophila exposed to dichlorodiphenyltrichloroethane (DDT), tributyltin (TBT), and 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD)

Dichlorodiphenyltrichloroethane (DDT), tributyltin (TBT), and 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) are persistent in the environment and cause continuous toxic effects in humans and aquatic life. Tetrahymena thermophila has the potential for use as a model for research regarding toxicants. In this study, this organism was used to analyze a genome-wide microarray generated from cells exposed to DDT, TBT and TCDD. To accomplish this, genes differentially expressed when treated with each toxicant were identified, after which their functions were categorized using GO enrichment analysis. The results suggested that the responses of T. thermophila were similar to those of multicellular organisms. Additionally, the context likelihood of relatedness method (CLR) was applied to construct a TCDD-relevant network. The T-shaped network obtained could be functionally divided into two subnetworks. The general functions of both subnetworks were related to the epigenetic mechanism of TCDD. Based on analysis of the networks, a model of the TCDD effect on T. thermophila was inferred. Thus, Tetrahymena has the potential to be a good unicellular eukaryotic model for toxic mechanism research at the genome level.

Gene network landscape of the ciliate Tetrahymena thermophila

Background

Genome-wide expression data of gene microarrays can be used to infer gene networks. At a cellular level, a gene network provides a picture of the modules in which genes are densely connected, and of the hub genes, which are highly connected with other genes. A gene network is useful to identify the genes involved in the same pathway, in a protein complex or that are co-regulated. In this study, we used different methods to find gene networks in the ciliate Tetrahymena thermophila, and describe some important properties of this network, such as modules and hubs.

Methodology/Principal Findings

Using 67 single channel microarrays, we constructed the Tetrahymena gene network (TGN) using three methods: the Pearson correlation coefficient (PCC), the Spearman correlation coefficient (SCC) and the context likelihood of relatedness (CLR) algorithm. The accuracy and coverage of the three networks were evaluated using four conserved protein complexes in yeast. The CLR network with a Z-score threshold 3.49 was determined to be the most robust. The TGN was partitioned, and 55 modules were found. In addition, analysis of the arbitrarily determined 1200 hubs showed that these hubs could be sorted into six groups according to their expression profiles. We also investigated human disease orthologs in Tetrahymena that are missing in yeast and provide evidence indicating that some of these are involved in the same process in Tetrahymena as in human.

Conclusions/Significance

This study constructed a Tetrahymena gene network, provided new insights to the properties of this biological network, and presents an important resource to study Tetrahymena genes at the pathway level.

Genome-wide identification and evolution of ATP-binding cassette transporters in the ciliate Tetrahymena thermophila: A case of functional divergence in a multigene family

BACKGROUND

In eukaryotes, ABC transporters that utilize the energy of ATP hydrolysis to expel cellular substrates into the environment are responsible for most of the efflux from cells. Many members of the superfamily of ABC transporters have been linked with resistance to multiple drugs or toxins. Owing to their medical and toxicological importance, members of the ABC superfamily have been studied in several model organisms and warrant examination in newly sequenced genomes.

RESULTS

A total of 165 ABC transporter genes, constituting a highly expanded superfamily relative to its size in other eukaryotes, were identified in the macronuclear genome of the ciliate Tetrahymena thermophila. Based on ortholog comparisons, phylogenetic topologies and intron characterizations, each highly expanded ABC transporter family of T. thermophila was classified into several distinct groups, and hypotheses about their evolutionary relationships are presented. A comprehensive microarray analysis revealed divergent expression patterns among the members of the ABC transporter superfamily during different states of physiology and development. Many of the relatively recently formed duplicate pairs within individual ABC transporter families exhibit significantly different expression patterns. Further analysis showed that multiple mechanisms have led to functional divergence that is responsible for the preservation of duplicated genes.

CONCLUSION

Gene duplications have resulted in an extensive expansion of the superfamily of ABC transporters in the Tetrahymena genome, making it the largest example of its kind reported in any organism to date. Multiple independent duplications and subsequent divergence contributed to the formation of different families of ABC transporter genes. Many of the members within a gene family exhibit different expression patterns. The combination of gene duplication followed by both sequence divergence and acquisition of new patterns of expression likely plays a role in the adaptation of Tetrahymen a to its environment.

Genome-wide identification and characterization of cytochrome P450 monooxygenase genes in the ciliate Tetrahymena thermophila

Background

Cytochrome P450 monooxygenases play key roles in the metabolism of a wide variety of substrates and they are closely associated with endocellular physiological processes or detoxification metabolism under environmental exposure. To date, however, none has been systematically characterized in the phylum Ciliophora. T. thermophila possess many advantages as a eukaryotic model organism and it exhibits rapid and sensitive responses to xenobiotics, making it an ideal model system to study the evolutionary and functional diversity of the P450 monooxygenase gene family.

Results

A total of 44 putative functional cytochrome P450 genes were identified and could be classified into 13 families and 21 sub-families according to standard nomenclature. The characteristics of both the conserved intron-exon organization and scaffold localization of tandem repeats within each P450 family clade suggested that the enlargement of T. thermophila P450 families probably resulted from recent separate small duplication events. Gene expression patterns of all T. thermophila P450s during three important cell physiological stages (vegetative growth, starvation and conjugation) were analyzed based on EST and microarray data, and three main categories of expression patterns were postulated. Evolutionary analysis including codon usage preference, site-specific selection and gene-expression evolution patterns were investigated and the results indicated remarkable divergences among the T. thermophila P450 genes.

Conclusion

The characterization, expression and evolutionary analysis of T. thermophila P450 monooxygenase genes in the current study provides useful information for understanding the characteristics and diversities of the P450 genes in the Ciliophora, and provides the baseline for functional analyses of individual P450 isoforms in this model ciliate species.

Refined annotation and assembly of the Tetrahymena thermophila genome sequence through EST analysis, comparative genomic hybridization, and targeted gap closure

BACKGROUND

Tetrahymena thermophila, a widely studied model for cellular and molecular biology, is a binucleated single-celled organism with a germline micronucleus (MIC) and somatic macronucleus (MAC). The recent draft MAC genome assembly revealed low sequence repetitiveness, a result of the epigenetic removal of invasive DNA elements found only in the MIC genome. Such low repetitiveness makes complete closure of the MAC genome a feasible goal, which to achieve would require standard closure methods as well as removal of minor MIC contamination of the MAC genome assembly. Highly accurate preliminary annotation of Tetrahymena's coding potential was hindered by the lack of both comparative genomic sequence information from close relatives and significant amounts of cDNA evidence, thus limiting the value of the genomic information and also leaving unanswered certain questions, such as the frequency of alternative splicing.

RESULTS

We addressed the problem of MIC contamination using comparative genomic hybridization with purified MIC and MAC DNA probes against a whole genome oligonucleotide microarray, allowing the identification of 763 genome scaffolds likely to contain MIC-limited DNA sequences. We also employed standard genome closure methods to essentially finish over 60% of the MAC genome. For the improvement of annotation, we have sequenced and analyzed over 60,000 verified EST reads from a variety of cellular growth and development conditions. Using this EST evidence, a combination of automated and manual reannotation efforts led to updates that affect 16% of the current protein-coding gene models. By comparing EST abundance, many genes showing apparent differential expression between these conditions were identified. Rare instances of alternative splicing and uses of the non-standard amino acid selenocysteine were also identified.

CONCLUSION

We report here significant progress in genome closure and reannotation of Tetrahymena thermophila. Our experience to date suggests that complete closure of the MAC genome is attainable. Using the new EST evidence, automated and manual curation has resulted in substantial improvements to the over 24,000 gene models, which will be valuable to researchers studying this model organism as well as for comparative genomics purposes.

Studies on the toxic effects of La(3+) to Tetrahymena thermophila by microcalorimetry

The toxic effects of La(3+) on Tetrahymena thermophila have been studied by microcalorimetry at 28 degrees C. The metabolic rate constant (r) and peak time were linked to the concentration of La(3+). The changes of metabolic rate constant indicated that low-concentration La(3+) (0-75 mg/L) had no significant effects on the metabolism of Tetrahymena cells but high-concentration La(3+) (100-175 mg/L) could inhibit their metabolism. From the results obtained by cell counting and fluorescence depolarization measurements, the inhibition of metabolism resulted from the decrease in cell number and the reduction in cell membrane fluidity. According to the results, it is clear that the metabolic mechanism of Tetrahymena cells has been changed with the addition of high-concentration La(3+). In addition, microcalorimetry of Tetrahymena could be a sensible, easy-to-use, and convenient method for monitoring the potential effects of rare earth elements on cells and the freshwater ecosystem.

Differentially expressed genes of Tetrahymena thermophila in response to tributyltin (TBT) identified by suppression subtractive hybridization and real time quantitative PCR

Tributyltin (TBT) is widely used as antifouling paints, agriculture biocides, and plastic stabilizers around the world, resulting in great pollution problem in aquatic environments. However, it has been short of the biomonitor to detect TBT in freshwater. We constructed the suppression subtractive hybridization library of Tetrahymena thermophila exposed to TBT, and screened out 101 Expressed Sequence Tags whose expressions were significantly up- or down-regulated with TBT treatment. From this, a series of genes related to the TBT toxicity were discovered, such as glutathione-S-transferase gene (down-regulated), plasma membrane Ca2+ ATPase isoforms 3 gene (up-regulated) and NgoA (up-regulated). Furthermore, their expressions under different concentrations of TBT treatment (0.5-40 ppb) were detected by real time fluorescent quantitative PCR. The differentially expressed genes of T. thermophila in response to TBT were identified, which provide the basic to make Tetrahymena as a sensitive, rapid and convenient TBT biomonitor in freshwater based on rDNA inducible expression system.

Molecular genetic analysis of an SNF2/brahma-related gene in Tetrahymena thermophila suggests roles in growth and nuclear development

We used a reverse genetic approach to identify three members of the SNF2 superfamily of chromatin remodeling genes in the ciliated protozoan Tetrahymena thermophila in order to investigate possible functions of ATP-dependent chromatin remodeling factors in growth and nuclear development. Comparative sequence analysis of the gene product of the Tetrahymena brahma-related gene (TtBRG1) indicates it is a member of the SNF2/BRM subgroup of the SNF2 superfamily. Northern analysis suggests that TtBRG1 has roles in growth and nuclear development in Tetrahymena. Indirect immunofluorescence analysis during nuclear development indicates that TtBrg1p localizes to both the parental and developing macronucleus of Tetrahymena during the time period corresponding to genome rearrangements. We generated germ line knockout heterokaryons for TtBRG1 and demonstrated that expression of the gene is required to complete nuclear development of Tetrahymena. In addition, the formation of distinct Pdd1p-containing structures is disturbed during the late stages of conjugation in TtBRG1 germ line knockout heterokaryons. We discuss these results in light of possible roles of SNF2-related proteins in growth and nuclear development of Tetrahymena.

Environmental cDNA analysis of the genes involved in lignocellulose digestion in the symbiotic protist community of Reticulitermes speratus

To clarify the lignocellulolytic process of the lower termite symbiotic protistan system, we constructed a cDNA library from an as yet uncultivated symbiotic protist community of the lower termite Reticulitermes speratus. The library was constructed by the biotinylated CAP trapper method and analyzed by one-pass sequencing. Phylogenetic analysis of actin orthologs confirmed that the resulting library reflected the intact organismal and mRNA composition of the symbiotic system. The contents of the library included abundant numbers of lignocellulolytic genes of the glycosyl hydrolase family orthologs (families 3, 5, 7, 8, 10, 11, 26, 43, 45 and 62). Our results clearly indicated that a multiple family of glycosyl hydrolase enzymes was involved in the protistan cellulose degradation system. The data also suggested that the most extensively expressed enzyme was glycosyl hydrolase family 7, a cellobiohydrolase ortholog. This family of enzymes enables the degradation of crystalline cellulose, the principal component of wood biomass.

Endoplasmic reticulum retention signal-dependent glycylation of the Hsp70/Grp170-related Pgp1p in Tetrahymena

Glycylation is an uncommon posttranslational modification. It has been found that tubulin glycylation is essential for cell survival in Tetrahymena. Here we describe PGP1, a Tetrahymena gene encoding an Hsp70 homologue that is a novel glycylated protein. Pgp1p is a conserved glycoprotein that localizes within the lumen of the endoplasmic reticulum (ER). We demonstrate that PGP1 is essential for viability and present evidence that both glycosylation and ER retention are necessary but not sufficient for glycylation.

Use of HAPPY mapping for the higher order assembly of the Tetrahymena genome

Tetrahymena thermophila is the best studied of the ciliates, a diversified and successful lineage of eukaryotic protists. Mirroring the way in which many metazoans partition their germ line and soma into distinct cell types, ciliates separate germ line and soma into two distinct nuclei in a single cell. The diploid, transcriptionally silent micronucleus undergoes meiosis and fertilization during sexual reproduction and determines the genotype of the progeny; in contrast, the expressed macronucleus contains many copies of hundreds of small chromosomes, determines the cell's phenotype, and is inherited only through vegetative reproduction. Here we demonstrate the power of HAPPY physical mapping to aid the complete assembly of T. thermophila macronuclear chromosomes from shotgun sequence scaffolds. The finished genome, one of only two ciliate genomes shotgun sequenced, will shed valuable additional light upon the biology of this extraordinary, diverse, and, from a genomics standpoint, as yet largely unexplored evolutionary branch of eukaryotes.

Identification of differentially expressed genes in Tetrahymena thermophila in response to dichlorodiphenyltrichloroethane (DDT) by suppression subtractive hybridization

The insecticide dichlorodiphenyltrichloroethane (DDT) is persistent in the environment, and continues to cause health problems. Tetrahymena has potential as a model organism for assaying low levels of DDT and for analysing the mechanisms of its toxicity. We constructed the suppression subtractive hybridization library of T. thermophila exposed to DDT, and screened out 90 Expressed Sequence Tags whose expressions were significantly up- or downregulated with DDT treatment. From this, a series of important genes related to the DDT metabolism and detoxification were discovered, such as P450 gene, glutathione S-transferase gene and sterol carrier protein 2 gene. Furthermore, their expressions under different concentrations of DDT treatment were detected by real-time fluorescent quantitative PCR. The results show that Tetrahymena is a relevant and useful model organism for detecting DDT in the environment and for discovering biomarkers that can be used to develop specific bio-reporters at the molecular and genomic levels.

Robust method for proteome analysis by MS/MS using an entire translated genome: demonstration on the ciliome of Tetrahymena thermophila

To improve the utility of increasingly large numbers of available unannotated and initially poorly annotated genomic sequences for proteome analysis, we demonstrate that effective protein identification can be made on a large and unannotated genome. The strategy developed is to translate the unannotated genome sequence into amino acid sequence encoding putative proteins in all six reading frames, to identify peptides by tandem mass spectrometry (MS/MS), to localize them on the genome sequence, and to preliminarily annotate the protein via a similarity search by BLAST. These tasks have been optimized and automated. Optimization to obtain multiple peptide matches in effect extends the searchable region and results in more robust protein identification. The viability of this strategy is demonstrated with the identification of 223 cilia proteins in the unicellular eukaryotic model organism Tetrahymena thermophila, whose initial genomic sequence draft was released in November 2003. To the best of our knowledge, this is the first demonstration of large-scale protein identification based on such a large, unannotated genome. Of the 223 cilia proteins, 84 have no similarity to proteins in NCBI's nonredundant (nr) database. This methodology allows identifying the locations of the genes encoding these novel proteins, which is a necessary first step to downstream functional genomic experimentation.

A non-long terminal repeat retrotransposon family is restricted to the germ line micronucleus of the ciliated protozoan Tetrahymena thermophila

The ciliated protozoan Tetrahymena thermophila undergoes extensive programmed DNA rearrangements during the development of a somatic macronucleus from the germ line micronucleus in its sexual cycle. To investigate the relationship between programmed DNA rearrangements and transposable elements, we identified several members of a family of non-long terminal repeat (LTR) retrotransposons (retroposons) in T. thermophila, the first characterized in the ciliated protozoa. This multiple-copy retrotransposon family is restricted to the micronucleus of T. thermophila. The REP (Tetrahymena non-LTR retroposon) elements encode an ORF2 typical of non-LTR elements that contains apurinic/apyrimidinic endonuclease (APE) and reverse transcriptase (RT) domains. Phylogenetic analysis of the RT and APE domains indicates that the element forms a deep-branching clade within the non-LTR retrotransposon family. Northern analysis with a probe to the conserved RT domain indicates that transcripts from the element are small and heterogeneous in length during early macronuclear development. The presence of a repeated transposable element in the genome is consistent with the model that programmed DNA deletion in T. thermophila evolved as a method of eliminating deleterious transposons from the somatic macronucleus.

Hypoxia regulates assembly of cilia in suppressors of Tetrahymena lacking an intraflagellar transport subunit gene

We cloned a Tetrahymena thermophila gene, IFT52, encoding a homolog of the Chlamydomonas intraflagellar transport protein, IFT52. Disruption of IFT52 led to loss of cilia and incomplete cytokinesis, a phenotype indistinguishable from that of mutants lacking kinesin-II, a known ciliary assembly transporter. The cytokinesis failures seem to result from lack of cell movement rather than from direct involvement of ciliary assembly pathway components in cytokinesis. Spontaneous partial suppressors of the IFT52 null mutants occurred, which assembled cilia at high cell density and resorbed cilia at low cell density. The stimulating effect of high cell density on cilia formation is based on the creation of pericellular hypoxia. Thus, at least under certain conditions, ciliary assembly is affected by an extracellular signal and the Ift52p function may be integrated into signaling pathways that regulate ciliogenesis.

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.