Complete Mitochondrial Genome Sequence of Three Tetrahymena Species Reveals Mutation Hot Spots and Accelerated Nonsynonymous Substitutions in Ymf Genes

The ciliate Tetrahymena, a model organism, contains divergent mitochondrial (Mt) genome with unusual properties, where half of its 44 genes still remain without a definitive function. These genes could be categorized into two major groups of KPC (known protein coding) and Ymf (genes without an identified function). To gain insights into the mechanisms underlying gene divergence and molecular evolution of Tetrahymena (T.) Mt genomes, we sequenced three Mt genomes of T.paravorax, T.pigmentosa, and T.malaccensis. These genomes were aligned and the analyses were carried out using several programs that calculate distance, nucleotide substitution (dn/ds), and their rate ratios (ω) on individual codon sites and via a sliding window approach. Comparative genomic analysis indicated a conserved putative transcription control sequence, a GC box, in a region where presumably transcription and replication initiate. We also found distinct features in Mt genome of T.paravorax despite similar genome organization among these ∼47 kb long linear genomes. Another significant finding was the presence of at least one or more highly variable regions in Ymf genes where majority of substitutions were concentrated. These regions were mutation hotspots where elevated distances and the dn/ds ratios were primarily due to an increase in the number of nonsynonymous substitutions, suggesting relaxed selective constraint. However, in a few Ymf genes, accelerated rates of nonsynonymous substitutions may be due to positive selection. Similarly, on protein level the majority of amino acid replacements occurred in these regions. Ymf genes comprise half of the genes in Tetrahymena Mt genomes, so understanding why they have not been assigned definitive functions is an important aspect of molecular evolution. Importantly, nucleotide substitution types and rates suggest possible reasons for not being able to find homologues for Ymf genes. Additionally, comparative genomic analysis of complete Mt genomes is essential in identifying biologically significant motifs such as control regions

Tetrahymena Metallothioneins Fall into Two Discrete Subfamilies

BACKGROUND: Metallothioneins are ubiquitous small, cysteine-rich, multifunctional proteins which can bind heavy metals. METHODOLOGY/PRINCIPAL FINDINGS: We report the results of phylogenetic and gene expression analyses that include two new Tetrahymena thermophila metallothionein genes (MTT3 and MTT5). Sequence alignments of all known Tetrahymena metallothioneins have allowed us to rationalize the structure of these proteins. We now formally subdivide the known metallothioneins from the ciliate genus Tetrahymena into two well defined subfamilies, 7a and 7b, based on phylogenetic analysis, on the pattern of clustering of Cys residues, and on the pattern of inducibility by the heavy metals Cd and Cu. Sequence alignment also reveals a remarkably regular, conserved and hierarchical modular structure of all five subfamily 7a MTs, which include MTT3 and MTT5. The former has three modules, while the latter has only two. Induction levels of the three T. thermophila genes were determined using quantitative real time RT-PCR. Various stressors (including heavy metals) brought about dramatically different fold-inductions for each gene; MTT5 showed the highest fold-induction. Conserved DNA motifs with potential regulatory significance were identified, in an unbiased way, upstream of the start codons of subfamily 7a MTs. EST evidence for alternative splicing in the 3′ UTR of the MTT5 mRNA with potential regulatory activity is reported. CONCLUSION/SIGNIFICANCE: The small number and remarkably regular structure of Tetrahymena MTs, coupled with the experimental tractability of this model organism for studies of in vivo function, make it an attractive system for the experimental dissection of the roles, structure/function relationships, regulation of gene expression, and adaptive evolution of these proteins, as well as for the development of biotechnological applications for the environmental monitoring of toxic substances

Distinct Functional Roles of β-Tubulin Isotypes in Microtubule Arrays of Tetrahymena thermophila, a Model Single-Celled Organism

<div><h3>Background</h3><p>The multi-tubulin hypothesis proposes that each tubulin isotype performs a unique role, or subset of roles, in the universe of microtubule function(s). To test this hypothesis, we are investigating the functions of the recently discovered, noncanonical β-like tubulins (BLTs) of the ciliate, <em>Tetrahymena thermophila</em>. <em>Tetrahymena</em> forms 17 distinct microtubular structures whose assembly had been thought to be based on single α- and β-isotypes. However, completion of the macronuclear genome sequence of <em>Tetrahymena</em> demonstrated that this ciliate possessed a β-tubulin multigene family: two synonymous genes (<em>BTU1</em> and <em>BTU2</em>) encode the canonical β-tubulin, BTU2, and six genes (<em>BLT1-6</em>) yield five divergent β-tubulin isotypes. In this report, we examine the structural features and functions of two of the BLTs (BLT1 and BLT4) and compare them to those of BTU2.</p> <h3>Methodology/Principal Findings</h3><p>With respect to BTU2, BLT1 and BLT4 had multiple sequence substitutions in their GTP-binding sites, in their interaction surfaces, and in their microtubule-targeting motifs, which together suggest that they have specialized functions. To assess the roles of these tubulins <em>in vivo</em>, we transformed <em>Tetrahymena</em> with expression vectors that direct the synthesis of GFP-tagged versions of the isotypes. We show that GFP-BLT1 and GFP-BLT4 were not detectable in somatic cilia and basal bodies, whereas GFP-BTU2 strongly labeled these structures. During cell division, GFP-BLT1 and GFP-BLT4, but not GFP-BTU2, were incorporated into the microtubule arrays of the macronucleus and into the mitotic apparatus of the micronucleus. GFP-BLT1 also participated in formation of the microtubules of the meiotic apparatus of the micronucleus during conjugation. Partitioning of the isotypes between nuclear and ciliary microtubules was confirmed biochemically.</p> <h3>Conclusion/Significance</h3><p>We conclude that <em>Tetrahymena</em> uses a family of distinct β-tubulin isotypes to construct subsets of functionally different microtubules, a result that provides strong support for the multi-tubulin hypothesis.</p> </div

Comprehensive Analysis Reveals Dynamic and Evolutionary Plasticity of Rab GTPases and Membrane Traffic in Tetrahymena thermophila

Cellular sophistication is not exclusive to multicellular organisms, and unicellular eukaryotes can resemble differentiated animal cells in their complex network of membrane-bound structures. These comparisons can be illuminated by genome-wide surveys of key gene families. We report a systematic analysis of Rabs in a complex unicellular Ciliate, including gene prediction and phylogenetic clustering, expression profiling based on public data, and Green Fluorescent Protein (GFP) tagging. Rabs are monomeric GTPases that regulate membrane traffic. Because Rabs act as compartment-specific determinants, the number of Rabs in an organism reflects intracellular complexity. The Tetrahymena Rab family is similar in size to that in humans and includes both expansions in conserved Rab clades as well as many divergent Rabs. Importantly, more than 90% of Rabs are expressed concurrently in growing cells, while only a small subset appears specialized for other conditions. By localizing most Rabs in living cells, we could assign the majority to specific compartments. These results validated most phylogenetic assignments, but also indicated that some sequence-conserved Rabs were co-opted for novel functions. Our survey uncovered a rare example of a nuclear Rab and substantiated the existence of a previously unrecognized core Rab clade in eukaryotes. Strikingly, several functionally conserved pathways or structures were found to be associated entirely with divergent Rabs. These pathways may have permitted rapid evolution of the associated Rabs or may have arisen independently in diverse lineages and then converged. Thus, characterizing entire gene families can provide insight into the evolutionary flexibility of fundamental cellular pathways