Whole genome studies of Tetrahymena

An Optimized and Versatile Counter-Flow Centrifugal Elutriation Workflow to Obtain Synchronized Eukaryotic Cells

Cell synchronization is a powerful tool to understand cell cycle events and its regulatory mechanisms. Counter-flow centrifugal elutriation (CCE) is a more generally desirable method to synchronize cells because it does not significantly alter cell behavior and/or cell cycle progression, however, adjusting specific parameters in a cell type/equipment-dependent manner can be challenging. In this paper, we used the unicellular eukaryotic model organism, Tetrahymena thermophila as a testing system for optimizing CCE workflow. Firstly, flow cytometry conditions were identified that reduced nuclei adhesion and improved the assessment of cell cycle stage. We then systematically examined how to achieve the optimal conditions for three critical factors affecting the outcome of CCE, including loading flow rate, collection flow rate and collection volume. Using our optimized workflow, we obtained a large population of highly synchronous G1-phase Tetrahymena as measured by 5-ethynyl-2'-deoxyuridine (EdU) incorporation into nascent DNA strands, bulk DNA content changes by flow cytometry, and cell cycle progression by light microscopy. This detailed protocol can be easily adapted to synchronize other eukaryotic cells.

A strategy for complete telomere-to-telomere assembly of ciliate macronuclear genome using ultra-high coverage Nanopore data

Ciliates contain two kinds of nuclei: the germline micronucleus (MIC) and the somatic macronucleus (MAC) in a single cell. The MAC usually have fragmented chromosomes. These fragmented chromosomes, capped with telomeres at both ends, could be gene size to several megabases in length among different ciliate species. So far, no telomere-to-telomere assembly of entire MAC genome in ciliate species has been finished. Development of the third generation sequencing technologies allows to generate sequencing reads up to megabases in length that could possibly span an entire MAC chromosome. Taking advantage of the ultra-long Nanopore reads, we established a simple strategy for the complete assembly of ciliate MAC genomes. Using this strategy, we assembled the complete MAC genomes of two ciliate species Tetrahymena thermophila and Tetrahymena shanghaiensis, composed of 181 and 214 chromosomes telomere-to-telomere respectively. The established strategy as well as the high-quality genome data will provide a useful approach for ciliate genome assembly, and a valuable community resource for further biological, evolutionary and population genomic studies.

The macronuclear genome of anaerobic ciliate Entodinium caudatum reveals its biological features adapted to the distinct rumen environment

Entodinium caudatum is an anaerobic binucleated ciliate representing the most dominant protozoal species in the rumen. However, its biological features are largely unknown due to the inability to establish an axenic culture. In this study, we primally sequenced its macronucleus (MAC) genome to aid the understanding of its metabolism, physiology, ecology. We isolated the MAC of E. caudatum strain MZG-1 and sequenced the MAC genome using Illumina MiSeq, MinION, and PacBio RSII systems. De novo assembly of the MiSeq sequence reads followed with subsequent scaffolding with MinION and PacBio reads resulted in a draft MAC genome about 117 Mbp. A large number of carbohydrate-active enzymes were likely acquired through horizontal gene transfer. About 8.74% of the E. caudatum predicted proteome was predicted as proteases. The MAC genome of E. caudatum will help better understand its important roles in rumen carbohydrate metabolism, and interaction with other members of the rumen microbiome.

The Compact Macronuclear Genome of the Ciliate Halteria grandinella: A Transcriptome-Like Genome with 23,000 Nanochromosomes

How to achieve protein diversity by genome and transcriptome processing is essential for organismal complexity and adaptation. The present work identifies that the macronuclear genome of Halteria grandinella, a cosmopolitan unicellular eukaryote, is composed almost entirely of gene-sized nanochromosomes with extremely short nongenic regions.

How to achieve protein diversity by genome and transcriptome processing is essential for organismal complexity and adaptation. The present work identifies that the macronuclear genome of Halteria grandinella, a cosmopolitan unicellular eukaryote, is composed almost entirely of gene-sized nanochromosomes with extremely short nongenic regions. This challenges our usual understanding of chromosomal structure and suggests the possibility of novel mechvanisms in transcriptional regulation. Comprehensive analysis of multiple data sets reveals that Halteria transcription dynamics are influenced by: (i) nonuniform nanochromosome copy numbers correlated with gene-expression level; (ii) dynamic alterations at both the DNA and RNA levels, including alternative internal eliminated sequence (IES) deletions during macronucleus formation and large-scale alternative splicing in transcript maturation; and (iii) extremely short 5′ and 3′ untranslated regions (UTRs) and universal TATA box-like motifs in the compact 5′ subtelomeric regions of most chromosomes. This study broadens the view of ciliate biology and the evolution of unicellular eukaryotes, and identifies Halteria as one of the most compact known eukaryotic genomes, indicating that complex cell structure does not require complex gene architecture.

Programmed genome rearrangements in ciliates

Ciliates are a highly divergent group of unicellular eukaryotes with separate somatic and germline genomes found in distinct dimorphic nuclei. This characteristic feature is tightly linked to extremely laborious developmentally regulated genome rearrangements in the development of a new somatic genome/nuclei following sex. The transformation from germline to soma genome involves massive DNA elimination mediated by non-coding RNAs, chromosome fragmentation, as well as DNA amplification. In this review, we discuss the similarities and differences in the genome reorganization processes of the model ciliates Paramecium and Tetrahymena (class Oligohymenophorea), and the distantly related Euplotes, Stylonychia, and Oxytricha (class Spirotrichea).

The completed macronuclear genome of a model ciliate Tetrahymena thermophila and its application in genome scrambling and copy number analyses

The ciliate Tetrahymena thermophila has been a powerful model system for molecular and cellular biology. However, some investigations have been limited due to the incomplete closure and sequencing of the macronuclear genome assembly, which for many years has been stalled at 1,158 scaffolds, with large sections of unknown sequences (available in Tetrahymena Genome Database, TGD, http://ciliate.org/ ). Here we completed the first chromosome-level Tetrahymena macronuclear genome assembly, with approximately 300× long Single Molecule, Real-Time reads of the wild-type SB210 cells-the reference strain for the initial macronuclear genome sequencing project. All 181 chromosomes were capped with two telomeres and gaps were entirely closed. The completed genome shows significant improvements over the current assembly (TGD 2014) in both chromosome structure and sequence integrity. The majority of previously identified gene models shown in TGD were retained, with the addition of 36 new genes and 883 genes with modified gene models. The new genome and annotation were incorporated into TGD. This new genome allows for pursuit in some underexplored areas that were far more challenging previously; two of them, genome scrambling and chromosomal copy number, were investigated in this study. We expect that the completed macronuclear genome will facilitate many studies in Tetrahymena biology, as well as multiple lines of research in other eukaryotes.

Evolution under pH stress and high population densities leads to increased density-dependent fitness in the protist Tetrahymena thermophila

Abiotic stress is a major force of selection that organisms are constantly facing. While the evolutionary effects of various stressors have been broadly studied, it is only more recently that the relevance of interactions between evolution and underlying ecological conditions, that is, eco-evolutionary feedbacks, have been highlighted. Here, we experimentally investigated how populations adapt to pH-stress under high population densities. Using the protist species Tetrahymena thermophila, we studied how four different genotypes evolved in response to stressfully low pH conditions and high population densities. We found that genotypes underwent evolutionary changes, some shifting up and others shifting down their intrinsic rates of increase (r₀ ). Overall, evolution at low pH led to the convergence of r₀ and intraspecific competitive ability (α) across the four genotypes. Given the strong correlation between r₀ and α, we argue that this convergence was a consequence of selection for increased density-dependent fitness at low pH under the experienced high density conditions. Increased density-dependent fitness was either attained through increase in r₀ , or decrease of α, depending on the genetic background. In conclusion, we show that demography can influence the direction of evolution under abiotic stress.

Insights into an Extensively Fragmented Eukaryotic Genome: De Novo Genome Sequencing of the Multinuclear Ciliate Uroleptopsis citrina

Ciliated protists are a large group of single-celled eukaryotes with separate germline and somatic nuclei in each cell. The somatic genome is developed from the zygotic nucleus through a series of chromosomal rearrangements, including fragmentation, DNA elimination, de novo telomere addition, and DNA amplification. This unique feature makes them perfect models for research in genome biology and evolution. However, genomic research of ciliates has been limited to a few species, owing to problems with DNA contamination and obstacles in cultivation. Here, we introduce a method combining telomere-primer PCR amplification and high-throughput sequencing, which can reduce DNA contamination and obtain genomic data efficiently. Based on this method, we report a draft somatic genome of a multimacronuclear ciliate, Uroleptopsis citrina. 1) The telomeric sequence in U. citrina is confirmed to be C4A4C4A4C4 by directly blunt-end cloning. 2) Genomic analysis of the resulting chromosomes shows a "one-gene one-chromosome" pattern, with a small number of multiple-gene chromosomes. 3) Amino acid usage is analyzed, and reassignment of stop codons is confirmed. 4) Chromosomal analysis shows an obvious asymmetrical GC skew and high bias between A and T in the subtelomeric regions of the sense-strand, with the detection of an 11-bp high AT motif region in the 3' subtelomeric region. 5) The subtelomeric sequence also has an obvious 40 nt strand oscillation of nucleotide ratio. 6) In the 5' subtelomeric region of the coding strand, the distribution of potential TATA-box regions is illustrated, which accumulate between 30 and 50 nt. This work provides a valuable reference for genomic research and furthers our understanding of the dynamic nature of unicellular eukaryotic genomes.

Whats, hows and whys of programmed DNA elimination in Tetrahymena

Programmed genome rearrangements in ciliates provide fascinating examples of flexible epigenetic genome regulations and important insights into the interaction between transposable elements (TEs) and host genomes. DNA elimination in Tetrahymena thermophila removes approximately 12 000 internal eliminated sequences (IESs), which correspond to one-third of the genome, when the somatic macronucleus (MAC) differentiates from the germline micronucleus (MIC). More than half of the IESs, many of which show high similarity to TEs, are targeted for elimination in cis by the small RNA-mediated genome comparison of the MIC to the MAC. Other IESs are targeted for elimination in trans by the same small RNAs through repetitive sequences. Furthermore, the small RNA–heterochromatin feedback loop ensures robust DNA elimination. Here, we review an updated picture of the DNA elimination mechanism, discuss the physiological and evolutionary roles of DNA elimination, and outline the key questions that remain unanswered.

Programmed Genome Rearrangements in the Ciliate Oxytricha

The ciliate Oxytricha is a microbial eukaryote with two genomes, one of which experiences extensive genome remodeling during development. Each round of conjugation initiates a cascade of events that construct a transcriptionally active somatic genome from a scrambled germline genome, with considerable help from both long and small noncoding RNAs. This process of genome remodeling entails massive DNA deletion and reshuffling of remaining DNA segments to form functional genes from their interrupted and scrambled germline precursors. The use of Oxytricha as a model system provides an opportunity to study an exaggerated form of programmed genome rearrangement. Furthermore, studying the mechanisms that maintain nuclear dimorphism and mediate genome rearrangement has demonstrated a surprising plasticity and diversity of noncoding RNA pathways, with new roles that go beyond conventional gene silencing. Another aspect of ciliate genetics is their unorthodox patterns of RNA-mediated, epigenetic inheritance that rival Mendelian inheritance. This review takes the reader through the key experiments in a model eukaryote that led to fundamental discoveries in RNA biology and pushes the biological limits of DNA processing.

Phosphorylation of an HP1-like Protein Regulates Heterochromatin Body Assembly for DNA Elimination

Heterochromatic loci are often assembled into higher-order heterochromatin bodies in diverse eukaryotes. However, the formation and biological roles of heterochromatin bodies are poorly understood. In the ciliated protozoan Tetrahymena, de novo heterochromatin body formation is accompanied by programmed DNA elimination. Here, we show that the heterochromatin body component Jub1p promotes heterochromatin body formation and dephosphorylation of the Heterochromatin Protein 1-like protein Pdd1p. Through the mutagenesis of the phosphorylated residues of Pdd1p, we demonstrate that Pdd1p dephosphorylation promotes the electrostatic interaction between Pdd1p and RNA in vitro and heterochromatin body formation in vivo. We therefore propose that heterochromatin body is assembled by the Pdd1p-RNA interaction. Pdd1p dephosphorylation and Jub1p are required for heterochromatin body formation and DNA elimination but not for local heterochromatin assembly, indicating that heterochromatin body plays an essential role in DNA elimination.

<mds_ies_db>: a database of ciliate genome rearrangements

Ciliated protists exhibit nuclear dimorphism through the presence of somatic macronuclei (MAC) and germline micronuclei (MIC). In some ciliates, DNA from precursor segments in the MIC genome rearranges to form transcriptionally active genes in the mature MAC genome, making these ciliates model organisms to study the process of somatic genome rearrangement. Similar broad scale, somatic rearrangement events occur in many eukaryotic cells and tumors. The <mds_ies_db> (http://oxytricha.princeton.edu/mds_ies_db) is a database of genome recombination and rearrangement annotations, and it provides tools for visualization and comparative analysis of precursor and product genomes. The database currently contains annotations for two completely sequenced ciliate genomes: Oxytricha trifallax and Tetrahymena thermophila.

DNA-guided establishment of nucleosome patterns within coding regions of a eukaryotic genome

A conserved hallmark of eukaryotic chromatin architecture is the distinctive array of well-positioned nucleosomes downstream from transcription start sites (TSS). Recent studies indicate that trans-acting factors establish this stereotypical array. Here, we present the first genome-wide in vitro and in vivo nucleosome maps for the ciliate Tetrahymena thermophila. In contrast with previous studies in yeast, we find that the stereotypical nucleosome array is preserved in the in vitro reconstituted map, which is governed only by the DNA sequence preferences of nucleosomes. Remarkably, this average in vitro pattern arises from the presence of subsets of nucleosomes, rather than the whole array, in individual Tetrahymena genes. Variation in GC content contributes to the positioning of these sequence-directed nucleosomes and affects codon usage and amino acid composition in genes. Given that the AT-rich Tetrahymena genome is intrinsically unfavorable for nucleosome formation, we propose that these “seed” nucleosomes—together with trans-acting factors—may facilitate the establishment of nucleosome arrays within genes in vivo, while minimizing changes to the underlying coding sequences.

Small-RNA-Mediated Genome-wide trans-Recognition Network in Tetrahymena DNA Elimination

Small RNAs are used to silence transposable elements (TEs) in many eukaryotes, which use diverse evolutionary solutions to identify TEs. In ciliated protozoans, small-RNA-mediated comparison of the germline and somatic genomes underlies identification of TE-related sequences, which are then eliminated from the soma. Here, we describe an additional mechanism of small-RNA-mediated identification of TE-related sequences in the ciliate Tetrahymena. We show that a limited set of internal eliminated sequences (IESs) containing potentially active TEs produces a class of small RNAs that recognize not only the IESs from which they are derived, but also other IESs in trans. This trans recognition triggers the expression of yet another class of small RNAs that identify other IESs. Therefore, TE-related sequences in Tetrahymena are robustly targeted for elimination by a genome-wide trans-recognition network accompanied by a chain reaction of small RNA production.

•

Two types of siRNAs (scnRNAs) are expressed in Tetrahymena sexual reproduction

•

Early-scnRNAs are produced from Type-A IESs containing potentially active transposons

•

Early-scnRNAs trans-recognize Type-B IESs and trigger late-scnRNA production in cis

•

Late-scnRNA production forms trans-recognition network for robust IES elimination

A domesticated PiggyBac transposase interacts with heterochromatin and catalyzes reproducible DNA elimination in Tetrahymena

The somatic genome of the ciliated protist Tetrahymena undergoes DNA elimination of defined sequences called internal eliminated sequences (IESs), which account for ~30% of the germline genome. During DNA elimination, IES regions are heterochromatinized and assembled into heterochromatin bodies in the developing somatic nucleus. The domesticated piggyBac transposase Tpb2p is essential for the formation of heterochromatin bodies and DNA elimination. In this study, we demonstrate that the activities of Tpb2p involved in forming heterochromatin bodies and executing DNA elimination are genetically separable. The cysteine-rich domain of Tpb2p, which interacts with the heterochromatin-specific histone modifications, is necessary for both heterochromatin body formation and DNA elimination, whereas the endonuclease activity of Tpb2p is only necessary for DNA elimination. Furthermore, we demonstrate that the endonuclease activity of Tpb2p in vitro and the endonuclease activity that executes DNA elimination in vivo have similar substrate sequence preferences. These results strongly indicate that Tpb2p is the endonuclease that directly catalyzes the excision of IESs and that the boundaries of IESs are at least partially determined by the combination of Tpb2p-heterochromatin interaction and relaxed sequence preference of the endonuclease activity of Tpb2p.

Role of class III phosphatidylinositol 3-kinase during programmed nuclear death of Tetrahymena thermophila

Programmed nuclear death (PND) in the ciliate protozoan Tetrahymena thermophila is a novel type of autophagy that occurs during conjugation, in which only the parental somatic macronucleus is destined to die and is then eliminated from the progeny cytoplasm. Other coexisting nuclei, however, such as new micro- and macronuclei are unaffected. PND starts with condensation in the nucleus followed by apoptotic DNA fragmentation, lysosomal acidification, and final resorption. Because of the peculiarity in the process and the absence of some ATG genes in this organism, the mechanism of PND has remained unclear. In this study, we focus on the role of class III phosphatidylinositol 3-kinase (PtdIns3K, corresponding to yeast Vps34) in order to identify central regulators of PND. We identified the sole Tetrahymena thermophila ortholog (TtVPS34) to yeast Vps34 and human PIK3C3 (the catalytic subunit of PtdIns3K), through phylogenetic analysis, and generated the gene knockdown mutant for functional analysis. Loss of TtVPS34 activity prevents autophagosome formation on the parental macronucleus, and this nucleus escapes from the lysosomal pathway. In turn, DNA fragmentation and final resorption of the nucleus are drastically impaired. These phenotypes are similar to the situation in the ATG8Δ mutants of Tetrahymena, implying an inextricable link between TtVPS34 and TtATG8s in controlling PND as well as general macroautophagy. On the other hand, TtVPS34 does not appear responsible for the nuclear condensation and does not affect the progeny nuclear development. These results demonstrate that TtVPS34 is critically involved in the nuclear degradation events of PND in autophagosome formation rather than with an involvement in commitment to the death program.

Transposon domestication versus mutualism in ciliate genome rearrangements

Ciliated protists rearrange their genomes dramatically during nuclear development via chromosome fragmentation and DNA deletion to produce a trimmer and highly reorganized somatic genome. The deleted portion of the genome includes potentially active transposons or transposon-like sequences that reside in the germline. Three independent studies recently showed that transposase proteins of the DDE/DDD superfamily are indispensible for DNA processing in three distantly related ciliates. In the spirotrich Oxytricha trifallax, high copy-number germline-limited transposons mediate their own excision from the somatic genome but also contribute to programmed genome rearrangement through a remarkable transposon mutualism with the host. By contrast, the genomes of two oligohymenophorean ciliates, Tetrahymena thermophila and Paramecium tetraurelia, encode homologous PiggyBac-like transposases as single-copy genes in both their germline and somatic genomes. These domesticated transposases are essential for deletion of thousands of different internal sequences in these species. This review contrasts the events underlying somatic genome reduction in three different ciliates and considers their evolutionary origins and the relationships among their distinct mechanisms for genome remodeling.

Genomes on the edge: programmed genome instability in ciliates

Ciliates are an ancient and diverse group of microbial eukaryotes that have emerged as powerful models for RNA-mediated epigenetic inheritance. They possess extensive sets of both tiny and long noncoding RNAs that, together with a suite of proteins that includes transposases, orchestrate a broad cascade of genome rearrangements during somatic nuclear development. This Review emphasizes three important themes: the remarkable role of RNA in shaping genome structure, recent discoveries that unify many deeply diverged ciliate genetic systems, and a surprising evolutionary "sign change" in the role of small RNAs between major species groups.

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

In Tetrahymena, a multi-sexed single-celled organism, the sex of the progeny is randomly determined by site-specific recombination events that assemble one complete gene pair and delete all others.

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 is a single-celled organism with seven sexes. After two cells of different sexes mate, the progeny cells can be of any one of the seven sexes. In this article we show how this sex decision is made. Every cell has two genomes, each contained within a separate nucleus. The germline genome is analogous to that in our ovaries or testes, containing all the genetic information for the sexual progeny; the somatic or working genome controls the operation of the cell (including its sex). We show that the germline genome contains a tandem array of similarly organized but incomplete gene pairs, one for each sex. Sex is chosen after fertilization when a new somatic genome is generated by rearrangement of a copy of the germline genome. One complete sex gene pair is assembled when the cell joins DNA segments at opposite ends of the array to each end of one incomplete gene pair; this gene pair is thus completed and becomes fully functional, while the remaining sex gene pairs are excised and lost. The process involves programmed, site-specific genome rearrangements, and the physically independent rearrangements that occur at opposite ends of the selected gene pair happen with high reliability and precision.

Piwi-interacting RNAs protect DNA against loss during Oxytricha genome rearrangement

Genome duality in ciliated protozoa offers a unique system to showcase their epigenome as a model of inheritance. In Oxytricha, the somatic genome is responsible for vegetative growth, whereas the germline contributes DNA to the next sexual generation. Somatic nuclear development removes all transposons and other so-called "junk" DNA, which comprise ~95% of the germline. We demonstrate that Piwi-interacting small RNAs (piRNAs) from the maternal nucleus can specify genomic regions for retention in this process. Oxytricha piRNAs map primarily to the somatic genome, representing the ~5% of the germline that is retained. Furthermore, injection of synthetic piRNAs corresponding to normally deleted regions leads to their retention in later generations. Our findings highlight small RNAs as powerful transgenerational carriers of epigenetic information for genome programming.

The Oxytricha trifallax macronuclear genome: a complex eukaryotic genome with 16,000 tiny chromosomes

With more chromosomes than any other sequenced genome, the macronuclear genome of Oxytricha trifallax has a unique and complex architecture, including alternative fragmentation and predominantly single-gene chromosomes.

The macronuclear genome of the ciliate Oxytricha trifallax displays an extreme and unique eukaryotic genome architecture with extensive genomic variation. During sexual genome development, the expressed, somatic macronuclear genome is whittled down to the genic portion of a small fraction (∼5%) of its precursor “silent” germline micronuclear genome by a process of “unscrambling” and fragmentation. The tiny macronuclear “nanochromosomes” typically encode single, protein-coding genes (a small portion, 10%, encode 2–8 genes), have minimal noncoding regions, and are differentially amplified to an average of ∼2,000 copies. We report the high-quality genome assembly of ∼16,000 complete nanochromosomes (∼50 Mb haploid genome size) that vary from 469 bp to 66 kb long (mean ∼3.2 kb) and encode ∼18,500 genes. Alternative DNA fragmentation processes ∼10% of the nanochromosomes into multiple isoforms that usually encode complete genes. Nucleotide diversity in the macronucleus is very high (SNP heterozygosity is ∼4.0%), suggesting that Oxytricha trifallax may have one of the largest known effective population sizes of eukaryotes. Comparison to other ciliates with nonscrambled genomes and long macronuclear chromosomes (on the order of 100 kb) suggests several candidate proteins that could be involved in genome rearrangement, including domesticated MULE and IS1595-like DDE transposases. The assembly of the highly fragmented Oxytricha macronuclear genome is the first completed genome with such an unusual architecture. This genome sequence provides tantalizing glimpses into novel molecular biology and evolution. For example, Oxytricha maintains tens of millions of telomeres per cell and has also evolved an intriguing expansion of telomere end-binding proteins. In conjunction with the micronuclear genome in progress, the O. trifallax macronuclear genome will provide an invaluable resource for investigating programmed genome rearrangements, complementing studies of rearrangements arising during evolution and disease.

The macronuclear genome of the ciliate Oxytricha trifallax, contained in its somatic nucleus, has a unique genome architecture. Unlike its diploid germline genome, which is transcriptionally inactive during normal cellular growth, the macronuclear genome is fragmented into at least 16,000 tiny (∼3.2 kb mean length) chromosomes, most of which encode single actively transcribed genes and are differentially amplified to a few thousand copies each. The smallest chromosome is just 469 bp, while the largest is 66 kb and encodes a single enormous protein. We found considerable variation in the genome, including frequent alternative fragmentation patterns, generating chromosome isoforms with shared sequence. We also found limited variation in chromosome amplification levels, though insufficient to explain mRNA transcript level variation. Another remarkable feature of Oxytricha's macronuclear genome is its inordinate fondness for telomeres. In conjunction with its possession of tens of millions of chromosome-ending telomeres per macronucleus, we show that Oxytricha has evolved multiple putative telomere-binding proteins. In addition, we identified two new domesticated transposase-like protein classes that we propose may participate in the process of genome rearrangement. The macronuclear genome now provides a crucial resource for ongoing studies of genome rearrangement processes that use Oxytricha as an experimental or comparative model.