Microarray Analyses of Gene Expression during the Tetrahymena thermophila Life Cycle

The model eukaryote, Tetrahymena thermophila, is the first ciliated protozoan whose genome has been sequenced, enabling genome-wide analysis of gene expression.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.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

<p>Abstract</p> <p>Background</p> <p><it>Tetrahymena thermophila</it>, 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 <it>Tetrahymena</it>'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.</p> <p>Results</p> <p>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.</p> <p>Conclusion</p> <p>We report here significant progress in genome closure and reannotation of <it>Tetrahymena thermophila</it>. 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.</p

The environmental impact of fertilizer embodied in a wheat-to-bread supply chain

Food production and consumption cause approximately one-third of total greenhouse gas emissions, and therefore delivering food security challenges not only the capacity of our agricultural system, but also its environmental sustainability. Knowing where and at what level environmental impacts occur within particular food supply chains is necessary if farmers, agri-food industries and consumers are to share responsibility to mitigate these impacts. Here we present an analysis of a complete supply chain for a staple of the global diet, a loaf of bread. We obtained primary data for all the processes involved in the farming, production and transport systems that lead to the manufacture of a particular brand of 800 g loaf. The data were analysed using an advanced life cycle assessment (LCA) tool, yielding metrics of environmental impact, including greenhouse gas emissions. We show that more than half of the environmental impact of producing the loaf of bread arises directly from wheat cultivation, with the use of ammonium nitrate fertilizer alone accounting for around 40%. These findings reveal the dependency of bread production on the unsustainable use of fertilizer and illustrate the detail needed if the actors in the supply chain are to assume shared responsibility for achieving sustainable food production

Kin5 Knockdown in Tetrahymena thermophila Using RNAi Blocks Cargo Transport of Gef1

A critical process that builds and maintains the eukaryotic cilium is intraflagellar transport (IFT). This process utilizes members of the kinesin-2 superfamily to transport cargo into the cilium (anterograde transport) and a dynein motor for the retrograde traffic. Using a novel RNAi knockdown method, we have analyzed the function of the homodimeric IFT kinesin-2, Kin5, in Tetrahymena ciliary transport. In RNAi transformants, Kin5 was severely downregulated and disappeared from the cilia, but cilia did not resorb, although tip structure was affected. After deciliation of the knockdown cell, cilia regrew and cells swam, which suggested that Kin5 is not responsible for the trafficking of axonemal precursors to build the cilium, but could be transporting molecules that act in ciliary signal transduction, such as guanine nucleotide exchange proteins (GEFs). Gef1 is a Tetrahymena ciliary protein, and current coimmunoprecipitation and immunofluorescence studies showed that it is absent in regrowing cilia of the knockdown cells lacking ciliary Kin5. We suggest that one important cargo of Kin5 is Gef1 and knockdown of Kin5 results in cell lethality

Transcriptome Analysis of the Model Protozoan, Tetrahymena thermophila, Using Deep RNA Sequencing

Background: The ciliated protozoan Tetrahymena thermophila is a well-studied single-celled eukaryote model organism for cellular and molecular biology. However, the lack of extensive T. thermophila cDNA libraries or a large expressed sequence tag (EST) database limited the quality of the original genome annotation. Methodology/Principal Findings: This RNA-seq study describes the first deep sequencing analysis of the T. thermophila transcriptome during the three major stages of the life cycle: growth, starvation and conjugation. Uniquely mapped reads covered more than 96 % of the 24,725 predicted gene models in the somatic genome. More than 1,000 new transcribed regions were identified. The great dynamic range of RNA-seq allowed detection of a nearly six order-of-magnitude range of measurable gene expression orchestrated by this cell. RNA-seq also allowed the first prediction of transcript untranslated regions (UTRs) and an updated (larger) size estimate of the T. thermophila transcriptome: 57 Mb, or about 55 % of the somatic genome. Our study identified nearly 1,500 alternative splicing (AS) events distributed over 5.2 % of T. thermophila genes. This percentage represents a two order-of-magnitude increase over previous EST-based estimates in Tetrahymena. Evidence of stage-specific regulation of alternative splicing was also obtained. Finally, our study allowed us to completely confirm about 26.8 % of the genes originally predicted by the gene finder, to correct coding sequence boundaries an

Structure of the germline genome of Tetrahymena thermophila and relationship to the massively rearranged somatic genome

The germline genome of the binucleated ciliate Tetrahymena thermophila undergoes programmed chromosome breakage and massive DNA elimination to generate the somatic genome. Here, we present a complete sequence assembly of the germline genome and analyze multiple features of its structure and its relationship to the somatic genome, shedding light on the mechanisms of genome rearrangement as well as the evolutionary history of this remarkable germline/soma differentiation. Our results strengthen the notion that a complex, dynamic, and ongoing interplay between mobile DNA elements and the host genome have shaped Tetrahymena chromosome structure, locally and globally. Non-standard outcomes of rearrangement events, including the generation of short-lived somatic chromosomes and excision of DNA interrupting protein-coding regions, may represent novel forms of developmental gene regulation. We also compare Tetrahymenas germline/soma differentiation to that of other characterized ciliates, illustrating the wide diversity of adaptations that have occurred within this phylum.</p

The Recombinases Rad51 and Dmc1 Play Distinct Roles in DNA Break Repair and Recombination Partner Choice in the Meiosis of Tetrahymena

Repair of programmed DNA double-strand breaks (DSBs) by meiotic recombination relies on the generation of flanking 3′ single-stranded DNA overhangs and their interaction with a homologous double-stranded DNA template. In various common model organisms, the ubiquitous strand exchange protein Rad51 and its meiosis-specific homologue Dmc1 have been implicated in the joint promotion of DNA–strand exchange at meiotic recombination sites. However, the division of labor between these two recombinases is still a puzzle. Using RNAi and gene-disruption experiments, we have studied their roles in meiotic recombination and chromosome pairing in the ciliated protist Tetrahymena as an evolutionarily distant meiotic model. Cytological and electrophoresis-based assays for DSBs revealed that, without Rad51p, DSBs were not repaired. However, in the absence of Dmc1p, efficient Rad51p-dependent repair took place, but crossing over was suppressed. Immunostaining and protein tagging demonstrated that only Dmc1p formed strong DSB–dependent foci on meiotic chromatin, whereas the distribution of Rad51p was diffuse within nuclei. This suggests that meiotic nucleoprotein filaments consist primarily of Dmc1p. Moreover, a proximity ligation assay confirmed that little if any Rad51p forms mixed nucleoprotein filaments with Dmc1p. Dmc1p focus formation was independent of the presence of Rad51p. The absence of Dmc1p did not result in compensatory assembly of Rad51p repair foci, and even artificial DNA damage by UV failed to induce Rad51p foci in meiotic nuclei, while it did so in somatic nuclei within one and the same cell. The observed interhomologue repair deficit in dmc1Δ meiosis is consistent with a requirement for Dmc1p in promoting the homologue as the preferred recombination partner. We propose that relatively short and/or transient Rad51p nucleoprotein filaments are sufficient for intrachromosomal recombination, whereas long nucleoprotein filaments consisting primarily of Dmc1p are required for interhomolog recombination

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

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

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