Erratum: Evolution of the mating type gene pair and multiple sexes in Tetrahymena.

[This corrects the article DOI: 10.1016/j.isci.2020.101950.].

Evolution of the mating type gene pair and multiple sexes in Tetrahymena.

The multiple mating type system of the Ciliate Tetrahymena thermophila is a self/non-self recognition system, whose specificity resides in a head-to-head, functionally distinct pair of genes, MTA and MTB. We have now sequenced and analyzed these mating type genes in nine additional Tetrahymena species. We conclude that MTA and MTB are derived from a common ancestral gene and have co-evolved for at least ∼150 Myr. We show that T. shanghaiensis, a perpetual selfer (unisexual) species, has a single mating type gene pair, whose MTA and MTB genes likely have different mating type specificity. We document the recent replacement of a complete different set of mating type specificities for another, illustrating how quickly this can happen. We discuss how varying conditions of reproductive stress could result in evolutionary co-adaptations of MTA and MTB genes and changes in mating type determination mechanisms.

A germline-limited piggyBac transposase gene is required for precise excision in Tetrahymena genome rearrangement.

Developmentally programmed genome rearrangement accompanies differentiation of the silent germline micronucleus into the transcriptionally active somatic macronucleus in the ciliated protozoan Tetrahymena thermophila. Internal eliminated sequences (IES) are excised, followed by rejoining of MAC-destined sequences, while fragmentation occurs at conserved chromosome breakage sequences, generating macronuclear chromosomes. Some macronuclear chromosomes, referred to as non-maintained chromosomes (NMC), are lost soon after differentiation. Large NMC contain genes implicated in development-specific roles. One such gene encodes the domesticated piggyBac transposase TPB6, required for heterochromatin-dependent precise excision of IES residing within exons of functionally important genes. These conserved exonic IES determine alternative transcription products in the developing macronucleus; some even contain free-standing genes. Examples of precise loss of some exonic IES in the micronucleus and retention of others in the macronucleus of related species suggest an evolutionary analogy to introns. Our results reveal that germline-limited sequences can encode genes with specific expression patterns and development-related functions, which may be a recurring theme in eukaryotic organisms experiencing programmed genome rearrangement during germline to soma differentiation.

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 Tetrahymena's germline/soma differentiation to that of other characterized ciliates, illustrating the wide diversity of adaptations that have occurred within this phylum.

Selecting one of several mating types through gene segment joining and deletion in Tetrahymena thermophila.

The unicellular eukaryote Tetrahymena thermophila has seven mating types. Cells can mate only when they recognize cells of a different mating type as non-self. As a ciliate, Tetrahymena separates its germline and soma into two nuclei. During growth the somatic nucleus is responsible for all gene transcription while the germline nucleus remains silent. During mating, a new somatic nucleus is differentiated from a germline nucleus and mating type is decided by a stochastic process. We report here that the somatic mating type locus contains a pair of genes arranged head-to-head. Each gene encodes a mating type-specific segment and a transmembrane domain that is shared by all mating types. Somatic gene knockouts showed both genes are required for efficient non-self recognition and successful mating, as assessed by pair formation and progeny production. The germline mating type locus consists of a tandem array of incomplete gene pairs representing each potential mating type. During mating, a complete new gene pair is assembled at the somatic mating type locus; the incomplete genes of one gene pair are completed by joining to gene segments at each end of germline array. All other germline gene pairs are deleted in the process. These programmed DNA rearrangements make this a fascinating system of mating type determination.

Tetrahymena thermophila, a unicellular eukaryote with separate germline and somatic genomes.

Tetrahymena thermophila is a ciliate--a unicellular eukaryote. Remarkably, every cell maintains differentiated germline and somatic genomes: one silent, the other expressed. Moreover, the two genomes undergo diverse processes, some as extreme as life and death, simultaneously in the same cytoplasm. Conserved eukaryotic mechanisms have been modified in ciliates to selectively deal with the two genomes. We describe research in several areas of Tetrahymena biology, including meiosis, amitosis, genetic assortment, selective nuclear pore transport, somatic RNAi-guided heterochromatin formation, DNA excision and programmed nuclear death by autophagy, which has enriched and broadened knowledge of those mechanisms.

Genome-Scale Analysis of Programmed DNA Elimination Sites in Tetrahymena thermophila.

Genetically programmed DNA rearrangements can regulate mRNA expression at an individual locus or, for some organisms, on a genome-wide scale. Ciliates rely on a remarkable process of whole-genome remodeling by DNA elimination to differentiate an expressed macronucleus (MAC) from a copy of the germline micronucleus (MIC) in each cycle of sexual reproduction. Here we describe results from the first high-throughput sequencing effort to investigate ciliate genome restructuring, comparing Sanger long-read sequences from a Tetrahymena thermophila MIC genome library to the MAC genome assembly. With almost 25% coverage of the unique-sequence MAC genome by MIC genome sequence reads, we created a resource for positional analysis of MIC-specific DNA removal that pinpoints MAC genome sites of DNA elimination at nucleotide resolution. The widespread distribution of internal eliminated sequences (IES) in promoter regions and introns suggests that MAC genome restructuring is essential not only for what it removes (for example, active transposons) but also for what it creates (for example, splicing-competent introns). Consistent with the heterogeneous boundaries and epigenetically modulated efficiency of individual IES deletions studied to date, we find that IES sites are dramatically under-represented in the ∼25% of the MAC genome encoding exons. As an exception to this general rule, we discovered a previously unknown class of small (<500 bp) IES with precise elimination boundaries that can contribute the 3' exon of an mRNA expressed during genome restructuring, providing a new mechanism for expanding mRNA complexity in a developmentally regulated manner.

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.

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.

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.

The highly conserved family of Tetrahymena thermophila chromosome breakage elements contains an invariant 10-base-pair core.

As a typical ciliate, Tetrahymena thermophila is a unicellular eukaryote that exhibits nuclear dimorphism: each cell contains a diploid, germ line micronucleus (MICN) and a polyploid, somatic macronucleus (MACN). During conjugation, when a new MACN differentiates from a mitotic descendant of the diploid fertilization nucleus, the five MICN chromosomes are site-specifically fragmented into 250 to 300 MACN chromosomes. The classic chromosome breakage sequence (CBS) is a 15-bp element (TAAACCAACCTCTTT) reported to be necessary and sufficient for chromosome breakage. To determine whether a CBS is present at every site of chromosome fragmentation and to assess the range of sequence variation tolerated, 31 CBSs were isolated without preconception as to the sequence of the chromosome breakage element. Additional CBS-related sequences were identified in the whole-genome sequence by their similarities to the classic CBS. Forty CBS elements behaved as authentic chromosome breakage sites. The CBS nucleotide sequence is more diverse than previously thought: nearly half of the CBS elements identified by unbiased methods have a variant of the classic CBS. Only an internal 10-bp core is completely conserved, but the entire 15-bp chromosome breakage sequence shows significant sequence conservation. Our results suggest that any one member of the CBS family provides a necessary and sufficient cis element for chromosome breakage. No chromosome breakage element totally unrelated to the classic CBS element was found; such elements, if they exist at all, must be rare.

Genome-wide characterization of Tetrahymena thermophila chromosome breakage sites. II. Physical and genetic mapping.

The chromosomes of the macronuclear (expressed) genome of Tetrahymena thermophila are generated by developmental fragmentation of the five micronuclear (germline) chromosomes. This fragmentation is site specific, directed by a conserved chromosome breakage sequence (Cbs element). An accompanying article in this issue reports the development of a successful scheme for the genome-wide cloning and identification of functional chromosome breakage sites. This article reports the physical and genetic characterization of 30 functional chromosome breakage junctions. Unique sequence tags and physical sizes were obtained for the pair of macronuclear chromosomes generated by fragmentation at each Cbs. Cbs-associated polymorphisms were used to genetically map 11 junctions to micronuclear linkage groups and macronuclear coassortment groups. Two pairs of junctions showed statistically significant similarity of the sequences flanking the Cbs, suggestive of relatively recent duplications of entire Cbs junctions during Tetrahymena genome evolution. Two macronuclear chromosomes that lose at least one end in an age-related manner were also identified. The whole-genome shotgun sequencing of the Tetrahymena macronucleus has recently been completed at The Institute for Genome Research (TIGR). By providing unique sequence from natural ends of macronuclear chromosomes, Cbs junctions will provide useful sequence tags for relating macro- and micronuclear genetic, physical, and whole-genome sequence maps.