Selenocysteine tRNA identification in the model organisms Dictyostelium discoideum and Tetrahymena thermophila

Molecular Characterization, Protein-Protein Interaction Network, and Evolution of Four Glutathione Peroxidases from Tetrahymena thermophila

Glutathione peroxidases (GPxs) form a broad family of antioxidant proteins essential for maintaining redox homeostasis in eukaryotic cells. In this study, we used an integrative approach that combines bioinformatics, molecular biology, and biochemistry to investigate the role of GPxs in reactive oxygen species detoxification in the unicellular eukaryotic model organism Tetrahymena thermophila. Both phylogenetic and mechanistic empirical model analyses provided indications about the evolutionary relationships among the GPXs of Tetrahymena and the orthologous enzymes of phylogenetically related species. In-silico gene characterization and text mining were used to predict the functional relationships between GPxs and other physiologically-relevant processes. The GPx genes contain conserved transcriptional regulatory elements in the promoter region, which suggest that transcription is under tight control of specialized signaling pathways. The bioinformatic findings were next experimentally validated by studying the time course of gene transcription and enzymatic activity after copper (Cu) exposure. Results emphasize the role of GPxs in the detoxification pathways that, by complex regulation of GPx gene expression, enable Tethraymena to survive in high Cu concentrations and the associated redox environment.

Tetrahymena Glutathione Peroxidase Family: A Comparative Analysis of These Antioxidant Enzymes and Differential Gene Expression to Metals and Oxidizing Agents

In the present work, an extensive analysis of the putative glutathione peroxidases (GPx) of the eukaryotic microorganism model Tetrahymena thermophila is carried out. A comparative analysis with GPx present in other Tetrahymena species and other very taxonomically diverse ciliates is also performed. A majority of ciliate GPx have replaced the selenocysteine (Sec) by Cys in its catalytic center, so they can be considered as phospholipid hydroperoxide glutathione peroxidases (PHGPx). Selenocysteine insertion sequence (SECIS) elements have been detected in several ciliate GPx that do not incorporate Sec in their amino acid sequences, and conversely, in other ciliate GPx with Sec, no SECIS elements are detected. These anomalies are analyzed and discussed. From the phylogenetic analysis using the ciliate GPx amino acid sequences, the existence of extensive intra- and interspecific gene duplications that produced multiple GPx isoforms in each species is inferred. The ancestral character of the selenoproteins is also corroborated. The analysis by qRT-PCR of six selected T. thermophila GPx genes has shown a quantitative differential expression between them, depending on the stressor (oxidizing agents, apoptotic inducer or metals) and the time of exposure.

Mutual detoxification of mercury and selenium in unicellular Tetrahymena

Selenium (Se) is commonly recognized as a protective element with an antagonistic effect against mercury (Hg) toxicity. However, the mechanisms of this Hg-Se antagonism are complex and remain controversial. To gain insight into the Hg-Se antagonism, a type of unicellular eukaryotic protozoa (Tetrahymena malaccensis, T. malaccensis) was selected and individually or jointly exposed to two Hg and three Se species. We found that Se species showed different toxic effects on the proliferation of T. malaccensis with the toxicity following the order: selenite (Se(IV))>selenomethionine (SeMeth)>selenate (Se(VI)). The Hg-Se antagonism in Tetrahymena was observed because the joint toxicity significantly decreased under co-exposure to highly toxic dosages of Hg and Se versus individual toxicity. Unlike Se(IV) and Se(VI), non-toxic dosage of SeMeth significantly decreased the Hg toxicity, revealing the influence of the Se species and dosages on the Hg-Se antagonism. Unexpectedly, inorganic divalent Hg (Hg²⁺) and monomethylmercury (MeHg) also displayed detoxification towards extremely highly toxic dosages of Se, although their detoxifying efficiency was discrepant. These results suggested mutual Hg-Se detoxification in T. malaccensis, which was highly dependent on the dosages and species of both elements. As compared to other species, SeMeth and MeHg promoted the Hg-Se joint effects to a higher degree. Additionally, the Hg contents decreased for all the Hg-Se co-exposed groups, revealing a sequestering effect of Se towards Hg in T. malaccensis.

Selenocysteine tRNA[Ser]Sec, the Central Component of Selenoprotein Biosynthesis: Isolation, Identification, Modification, and Sequencing

The selenocysteine (Sec) tRNA^[Ser]Sec population consists of two isoforms that differ from each other by a single 2'-O-methylribosyl moiety at position 34 (Um34). These two isoforms, which are encoded in a single gene, Trsp, and modified posttranscriptionally, are involved individually in the synthesis of two subclasses of selenoproteins, designated housekeeping and stress-related selenoproteins. Techniques used in obtaining these isoforms for their characterization include extraction of RNA from mammalian cells and tissues, purifying the tRNA^[Ser]Sec population by one or more procedures, and finally resolving the two isoforms from each other. Since some of the older techniques for isolating tRNA^[Ser]Sec and resolving the isoforms are used in only a few laboratories, these procedures will be discussed briefly and references provided for more detailed information, while the more recently developed procedures are discussed in detail. In addition, a novel technique that was developed in sequencing tRNA^[Ser]Sec for identifying their occurrence in other organisms is also presented.

Genome-wide evolutionary analysis of the noncoding RNA genes and noncoding DNA of Paramecium tetraurelia

The compact genome of the unicellular eukaryote Paramecium tetraurelia contains noncoding DNA (ncDNA) distributed into >39,000 intergenic sequences and >90,000 introns of 390 base pairs (bp) and 25 bp on average, respectively. Here we analyzed the molecular features of the ncRNA genes, introns, and intergenic sequences of this genome. We mainly used computational programs and comparative genomics possible because the P. tetraurelia genome had formed throughout whole-genome duplications (WGDs). We characterized 417 5S rRNA, snRNA, snoRNA, SRP RNA, and tRNA putative genes, 415 of which map within intergenic sequences, and two, within introns. The evolution of these ncRNA genes appears to have mainly involved purifying selection and gene deletion. We then compared the introns that interrupt the protein-coding gene duplicates arisen from the recent WGD and identified a population of a few thousands of introns having evolved under most stringent constraints (>95% of identity). We also showed that low nucleotide substitution levels characterize the 50 and 80-115 base pairs flanking, respectively, the stop and start codons of the protein-coding genes. Lower substitution levels mark the base pairs flanking the highly transcribed genes, or the start codons of the genes of the sets with a high number of WGD-related sequences. Finally, adjacent to protein-coding genes, we characterized 32 DNA motifs able to encode stable and evolutionary conserved RNA secondary structures and defining putative expression controlling elements. Fourteen DNA motifs with similar properties map distant from protein-coding genes and may encode regulatory ncRNAs.

Microarray analyses of gene expression during the Tetrahymena thermophila life cycle

BACKGROUND

The model eukaryote, Tetrahymena thermophila, is the first ciliated protozoan whose genome has been sequenced, enabling genome-wide analysis of gene expression.

METHODOLOGY/PRINCIPAL FINDINGS

A genome-wide microarray platform containing the predicted coding sequences (putative genes) for T. thermophila is described, validated and used to study gene expression during the three major stages of the organism's life cycle: growth, starvation and conjugation.

CONCLUSIONS/SIGNIFICANCE

Of the approximately 27,000 predicted open reading frames, transcripts homologous to only approximately 5900 are not detectable in any of these life cycle stages, indicating that this single-celled organism does indeed contain a large number of functional genes. Transcripts from over 5000 predicted genes are expressed at levels >5x corrected background and 95 genes are expressed at >250x corrected background in all stages. Transcripts homologous to 91 predicted genes are specifically expressed and 155 more are highly up-regulated in growing cells, while 90 are specifically expressed and 616 are up-regulated during starvation. Strikingly, transcripts homologous to 1068 predicted genes are specifically expressed and 1753 are significantly up-regulated during conjugation. The patterns of gene expression during conjugation correlate well with the developmental stages of meiosis, nuclear differentiation and DNA elimination. The relationship between gene expression and chromosome fragmentation is analyzed. Genes encoding proteins known to interact or to function in complexes show similar expression patterns, indicating that co-ordinate expression with putative genes of known function can identify genes with related functions. New candidate genes associated with the RNAi-like process of DNA elimination and with meiosis are identified and the late stages of conjugation are shown to be characterized by specific expression of an unexpectedly large and diverse number of genes not involved in nuclear functions.

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.

Evolutionary dynamics of eukaryotic selenoproteomes: large selenoproteomes may associate with aquatic life and small with terrestrial life

In silico and metabolic labeling studies of the selenoproteomes of several eukaryotes revealed distinct selenoprotein patterns as well as an ancient origin of selenoproteins and massive, independent losses in land plants, fungi, nematodes, insects and some protists, suggesting that the environment plays an important role in selenoproteome evolution.

Background

Selenocysteine (Sec) is a selenium-containing amino acid that is co-translationally inserted into nascent polypeptides by recoding UGA codons. Selenoproteins occur in both eukaryotes and prokaryotes, but the selenoprotein content of organisms (selenoproteome) is highly variable and some organisms do not utilize Sec at all.

Results

We analyzed the selenoproteomes of several model eukaryotes and detected 26 and 29 selenoprotein genes in the green algae Ostreococcus tauri and Ostreococcus lucimarinus, respectively, five in the social amoebae Dictyostelium discoideum, three in the fly Drosophila pseudoobscura, and 16 in the diatom Thalassiosira pseudonana, including several new selenoproteins. Distinct selenoprotein patterns were verified by metabolic labeling of O. tauri and D. discoideum with ⁷⁵Se. More than half of the selenoprotein families were shared by unicellular eukaryotes and mammals, consistent with their ancient origin. Further analyses identified massive, independent selenoprotein losses in land plants, fungi, nematodes, insects and some protists. Comparative analyses of selenoprotein-rich and -deficient organisms revealed that aquatic organisms generally have large selenoproteomes, whereas several groups of terrestrial organisms reduced their selenoproteomes through loss of selenoprotein genes and replacement of Sec with cysteine.

Conclusion

Our data suggest many selenoproteins originated at the base of the eukaryotic domain and show that the environment plays an important role in selenoproteome evolution. In particular, aquatic organisms apparently retained and sometimes expanded their selenoproteomes, whereas the selenoproteomes of some terrestrial organisms were reduced or completely lost. These findings suggest a hypothesis that, with the exception of vertebrates, aquatic life supports selenium utilization, whereas terrestrial habitats lead to reduced use of this trace element due to an unknown environmental factor.

De novo search for non-coding RNA genes in the AT-rich genome of Dictyostelium discoideum: performance of Markov-dependent genome feature scoring

Genome data are increasingly important in the computational identification of novel regulatory non-coding RNAs (ncRNAs). However, most ncRNA gene-finders are either specialized to well-characterized ncRNA gene families or require comparisons of closely related genomes. We developed a method for de novo screening for ncRNA genes with a nucleotide composition that stands out against the background genome based on a partial sum process. We compared the performance when assuming independent and first-order Markov-dependent nucleotides, respectively, and used Karlin-Altschul and Karlin-Dembo statistics to evaluate the significance of hits. We hypothesized that a first-order Markov-dependent process might have better power to detect ncRNA genes since nearest-neighbor models have been shown to be successful in predicting RNA structures. A model based on a first-order partial sum process (analyzing overlapping dinucleotides) had better sensitivity and specificity than a zeroth-order model when applied to the AT-rich genome of the amoeba Dictyostelium discoideum. In this genome, we detected 94% of previously known ncRNA genes (at this sensitivity, the false positive rate was estimated to be 25% in a simulated background). The predictions were further refined by clustering candidate genes according to sequence similarity and/or searching for an ncRNA-associated upstream element. We experimentally verified six out of 10 tested ncRNA gene predictions. We conclude that higher-order models, in combination with other information, are useful for identification of novel ncRNA gene families in single-genome analysis of D. discoideum. Our generalizable approach extends the range of genomic data that can be searched for novel ncRNA genes using well-grounded statistical methods.

A genomic survey of the fish parasite Spironucleus salmonicida indicates genomic plasticity among diplomonads and significant lateral gene transfer in eukaryote genome evolution

BACKGROUND

Comparative genomic studies of the mitochondrion-lacking protist group Diplomonadida (diplomonads) has been lacking, although Giardia lamblia has been intensively studied. We have performed a sequence survey project resulting in 2341 expressed sequence tags (EST) corresponding to 853 unique clones, 5275 genome survey sequences (GSS), and eleven finished contigs from the diplomonad fish parasite Spironucleus salmonicida (previously described as S. barkhanus).

RESULTS

The analyses revealed a compact genome with few, if any, introns and very short 3' untranslated regions. Strikingly different patterns of codon usage were observed in genes corresponding to frequently sampled ESTs versus genes poorly sampled, indicating that translational selection is influencing the codon usage of highly expressed genes. Rigorous phylogenomic analyses identified 84 genes--mostly encoding metabolic proteins--that have been acquired by diplomonads or their relatively close ancestors via lateral gene transfer (LGT). Although most acquisitions were from prokaryotes, more than a dozen represent likely transfers of genes between eukaryotic lineages. Many genes that provide novel insights into the genetic basis of the biology and pathogenicity of this parasitic protist were identified including 149 that putatively encode variant-surface cysteine-rich proteins which are candidate virulence factors. A number of genomic properties that distinguish S. salmonicida from its human parasitic relative G. lamblia were identified such as nineteen putative lineage-specific gene acquisitions, distinct mutational biases and codon usage and distinct polyadenylation signals.

CONCLUSION

Our results highlight the power of comparative genomic studies to yield insights into the biology of parasitic protists and the evolution of their genomes, and suggest that genetic exchange between distantly-related protist lineages may be occurring at an appreciable rate in eukaryote genome evolution.

Macronuclear genome sequence of the ciliate Tetrahymena thermophila, a model eukaryote

The ciliate Tetrahymena thermophila is a model organism for molecular and cellular biology. Like other ciliates, this species has separate germline and soma functions that are embodied by distinct nuclei within a single cell. The germline-like micronucleus (MIC) has its genome held in reserve for sexual reproduction. The soma-like macronucleus (MAC), which possesses a genome processed from that of the MIC, is the center of gene expression and does not directly contribute DNA to sexual progeny. We report here the shotgun sequencing, assembly, and analysis of the MAC genome of T. thermophila, which is approximately 104 Mb in length and composed of approximately 225 chromosomes. Overall, the gene set is robust, with more than 27,000 predicted protein-coding genes, 15,000 of which have strong matches to genes in other organisms. The functional diversity encoded by these genes is substantial and reflects the complexity of processes required for a free-living, predatory, single-celled organism. This is highlighted by the abundance of lineage-specific duplications of genes with predicted roles in sensing and responding to environmental conditions (e.g., kinases), using diverse resources (e.g., proteases and transporters), and generating structural complexity (e.g., kinesins and dyneins). In contrast to the other lineages of alveolates (apicomplexans and dinoflagellates), no compelling evidence could be found for plastid-derived genes in the genome. UGA, the only T. thermophila stop codon, is used in some genes to encode selenocysteine, thus making this organism the first known with the potential to translate all 64 codons in nuclear genes into amino acids. We present genomic evidence supporting the hypothesis that the excision of DNA from the MIC to generate the MAC specifically targets foreign DNA as a form of genome self-defense. The combination of the genome sequence, the functional diversity encoded therein, and the presence of some pathways missing from other model organisms makes T. thermophila an ideal model for functional genomic studies to address biological, biomedical, and biotechnological questions of fundamental importance.

Selenoprotein synthesis: UGA does not end the story

It is well established that the beneficial effects of the trace element selenium are mediated by its major biological product, the amino acid selenocysteine, present in the active site of selenoproteins. These fulfill different functions, as varied as oxidation-reduction of metabolites in bacteria, reduction of reactive oxygen species, control of the redox status of the cell or thyroid hormone maturation. This review will focus on the singularities of the selenocysteine biosynthesis pathway and its unique incorporation mechanism into eukaryal selenoproteins. Selenocysteine biosynthesis from serine is achieved on tRNA(Sec) and requires four proteins. As this amino acid is encoded by an in-frame UGA codon, otherwise signaling termination of translation, ribosomes must be told not to stop at this position in the mRNA. Several molecular partners acting in cis or in trans have been identified, but their knowledge has not enabled yet to firmly establish the molecular events underlying this mechanism. Data suggest that other, so far uncharacterized factors might exist. In this survey, we attempted to compile all the data available in the literature and to describe the latest developments in the field.

The Plasmodium selenoproteome

The use of selenocysteine (Sec) as the 21st amino acid in the genetic code has been described in all three major domains of life. However, within eukaryotes, selenoproteins are only known in animals and algae. In this study, we characterized selenoproteomes and Sec insertion systems in protozoan Apicomplexa parasites. We found that among these organisms, Plasmodium and Toxoplasma utilized Sec, whereas Cryptosporidium did not. However, Plasmodium had no homologs of known selenoproteins. By searching computationally for evolutionarily conserved selenocysteine insertion sequence (SECIS) elements, which are RNA structures involved in Sec insertion, we identified four unique Plasmodium falciparum selenoprotein genes. These selenoproteins were incorrectly annotated in PlasmoDB, were conserved in other Plasmodia and had no detectable homologs in other species. We provide evidence that two Plasmodium SECIS elements supported Sec insertion into parasite and endogenous selenoproteins when they were expressed in mammalian cells, demonstrating that the Plasmodium SECIS elements are functional and indicating conservation of Sec insertion between Apicomplexa and animals. Dependence of the plasmodial parasites on selenium suggests possible strategies for antimalarial drug development.