Three different macronuclear DNAs in Oxytricha fallax share a common sequence block

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.

Characterization and taxonomic validity of the ciliate Oxytricha trifallax (Class Spirotrichea) based on multiple gene sequences: limitations in identifying genera solely by morphology

Oxytricha trifallax - an established model organism for studying genome rearrangements, chromosome structure, scrambled genes, RNA-mediated epigenetic inheritance, and other phenomena - has been the subject of a nomenclature controversy for several years. Originally isolated as a sibling species of O. fallax, O. trifallax was reclassified in 1999 as Sterkiella histriomuscorum, a previously identified species, based on morphological similarity. The proper identification of O. trifallax is crucial to resolve in order to prevent confusion in both the comparative genomics and the general scientific communities. We analyzed nine conserved nuclear gene sequences between the two given species and several related ciliates. Phylogenetic analyses suggest that O. trifallax and a bona fide S. histriomuscorum have accumulated significant evolutionary divergence from each other relative to other ciliates such that they should be unequivocally classified as separate species. We also describe the original isolation of O. trifallax, including its comparison to O. fallax, and we provide criteria to identify future isolates of O. trifallax.

New Insights into the Macronuclear Development in Ciliates

During macronuclear differentiation in ciliated protozoa, most amazing "DNA gymnastics" takes place, which includes DNA excision, DNA elimination, DNA reorganization, and DNA-specific amplification. Although the morphological events occurring during macronuclear development are well described, a detailed knowledge of the molecular mechanisms and the regulation of this differentiation process is still missing. However, recently several models have been proposed for the molecular regulation of macronuclear differentiation, but these models have yet to be verified experimentally. The scope of this review is to summarize recent discoveries in different ciliate species and to compare and discuss the different models proposed. Results obtained in these studies are not only relevant for our understanding of nuclear differentiation in ciliates, but also for cellular differentiation in eukaryotic organisms in general as well as for other disciplines such as bioinformatics and computational biology.

Polymorphism, Recombination and Alternative Unscrambling in the DNA Polymerase α Gene of the Ciliate Stylonychia lemnae (Alveolata; class Spirotrichea)

DNA polymerase α is the most highly scrambled gene known in stichotrichous ciliates. In its hereditary micronuclear form, it is broken into >40 pieces on two loci at least 3 kb apart. Scrambled genes must be reassembled through developmental DNA rearrangements to yield functioning macronuclear genes, but the mechanism and accuracy of this process are unknown. We describe the first analysis of DNA polymorphism in the macronuclear version of any scrambled gene. Six functional haplotypes obtained from five Eurasian strains of Stylonychia lemnae were highly polymorphic compared to Drosophila genes. Another incompletely unscrambled haplotype was interrupted by frameshift and nonsense mutations but contained more silent mutations than expected by allelic inactivation. In our sample, nucleotide diversity and recombination signals were unexpectedly high within a region encompassing the boundary of the two micronuclear loci. From this and other evidence we infer that both members of a long repeat at the ends of the loci provide alternative substrates for unscrambling in this region. Incongruent genealogies and recombination patterns were also consistent with separation of the two loci by a large genetic distance. Our results suggest that ciliate developmental DNA rearrangements may be more probabilistic and error prone than previously appreciated and constitute a potential source of macronuclear variation. From this perspective we introduce the nonsense-suppression hypothesis for the evolution of ciliate altered genetic codes. We also introduce methods and software to calculate the likelihood of hemizygosity in ciliate haplotype samples and to correct for multiple comparisons in sliding-window analyses of Tajima's D.

Genome remodeling in ciliated protozoa

The germline genomes of ciliated protozoa are dynamic structures, undergoing massive DNA rearrangement during the formation of a functional macronucleus. Macronuclear development involves chromosome fragmentation coupled with de novo telomere synthesis, numerous DNA splicing events that remove internal segments of DNA, and, in some ciliates, the reordering of scrambled gene segments. Despite the fact that all ciliates share similar forms of DNA rearrangement, there appears to be great diversity in both the nature of the rearranged DNA and the molecular mechanisms involved. Epigenetic effects on rearrangement have also been observed, and recent work suggests that chromatin differentiation plays a role in specifying DNA segments either for rearrangement or for elimination.

Two two-gene macronuclear chromosomes of the hypotrichous ciliates Oxytricha fallax and O. trifallax generated by alternative processing of the 81 locus

We describe the first know macronuclear chromosomes that carry more than one gene in hypotrichous ciliated protozoa. These 4.9- and 2.8-kbp chromosomes each consist almost exclusively of two protein-coding genes, which are conserved and transcribed. The two chromosomes share a common region that consists of a gene that is a member of the family of mitochondrial solute carrier genes (CR-MSC; [Williams and Herrick (1991): Nucleic Acids Res 19:4717-4724]. Each chromosome also carries another gene appended to its common region: The 4.9-kbp chromosome also carries a gene that encodes a protein that is rich in glutamine and charged amino acids and bears regions of heptad repeats characteristic of coiled-coils. Its function is unknown. The second gene of the 2.8 kbp chromosome is a mitochondrial solute carrier gene (LA-MSC); thus, the 2.8-kbp chromosomes consists of two mitochondrial solute carrier paralogs. Phylogenetic analysis indicates that the two genes were duplicated before ciliates diverged from the main eukaryotic lineage and were subsequently juxtaposed. The CR- and LA-MSC genes are each interrupted by three introns. The introns are not in homologous positions, suggesting that they may have originated from multiple group II intron transpositions. These chromosomes and their genes are encoded in the Oxytricha germline by the 81 locus. This locus is alternatively processed to generate a nested set of three macronuclear chromosomes, the 4.9- and 2.8-kbp chromosomes and a third (1.6 kbp) which consists almost exclusively of the shared common gene, CR-MSC. Such alternative processing is common in macronuclear development of O. fallax [Cartinhour and Herrick (1984): Mol Cell Biol 4:931-938]. Possible functions for alternative processing are considered; e.g., it may serve to physically link genes to allow co-regulation or co-replication by a common cis-acting sequence.

A short internal eliminated sequence with central conserved sequences interrupting the LA-MSC gene of the 81 locus in the hypotrichous ciliates Oxytricha fallax and O. trifallax

IES-LA is a short Internal Eliminated Sequence interrupting LA-MSC, a protein-coding gene of the 81 locus of Oxytricha fallax and O. trifallax. IES-LA is precisely excised from the gene during development of the macronucleus. The internal eliminated sequence is bounded by CAAT ... AATG, and thereby resembles a TBE1 transposon internal eliminated sequence insertion that is grossly shortened (4.1 kbp to 52-64 bp), consistent with the hypothesis that short IESs are degenerated ancient transposons. The pattern of sequence conservation between five alleles of IES-LA shows that it differs from previously characterized classes of ciliate short IESs: while many short IESs have conserved ends and diverged centers, IES-LA is more conserved in its center and its ends are diverged. This implies a excision mechanism for IES-LA that is distinct from those for other known Oxytricha IESs.

A high degree of macronuclear chromosome polymorphism is generated by variable DNA rearrangements in Paramecium primaurelia during macronuclear differentiation

DNA rearrangements in Paramecium lead to the formation of macronuclear chromosomes, the sizes of which range from 50 and 800 kb (1 kb is 10(3) base-pairs). This process does not appear to be a simple size reduction of the micronuclear chromosomes by specific and reproducible DNA sequence elimination and chromosomal breakage followed by chromosomal amplification. On the contrary, this process generates a variety of different, but sequence-related, macronuclear chromosomes from a unique set of micronuclear chromosomes. This paper describes an attempt to understand the nature of the diversity of the macronuclear chromosomes and the mechanisms of their production. The structure of three macronuclear chromosomes, 480, 250 and 230 kb in size, have been determined utilizing chromosome-jumping and YAC-cloning techniques. The two smallest chromosomes correspond roughly to the two halves of the longest chromosome. The main contribution to the diversity arises from the chromosomal ends and is due to variable positions of the telomere addition sites and/or to variable rearrangements of DNA sequences. The 480 kb chromosome contains a region of variable length, which is likely to be due to a variable deletion, located at the position of telomerization seen in the two small chromosomes. A model of chromosomal breakage is proposed to rationalize this result where micronuclear DNA is first amplified, broken and degraded to various extent from the newly formed ends, which subsequently are either telomerized or religated. Potential implications of these processes for gene expression is discussed. Known phenotypes that have a macronuclear determinism could be explained by this type of process.

Non-coding DNA in macronuclear chromosomes of hypotrichous ciliates

Massive elimination of sequences occurs in the development of the macronucleus of hypotrichous ciliates. The surviving sequences are presumed to have functions in the macronucleus; what little is known about non-coding macronuclear sequences is reviewed. The 1.7 kbp macronuclear chromosome that carries a histone H4 gene consists primarily of non-coding DNA 5' of the histone gene. This region is shown by sequence comparison to carry several perfectly conserved sequence blocks up to 14 bp long, scattered amongst regions which have evolved greatly since the divergence of Oxytricha nova and Stylonychia lemnae. This result is consistent with the suggestion of Harper and Jahn [Harper, D. S. & Jahn, C. L. 1989. Actin, tubulin and H4 histone genes in three species of hypotrichous ciliated protozoa. Gene, 75:93-107] that this large non-coding 5' region may be involved in the transcriptional regulation of the histone H4 gene carried on the 1.7 kbp chromosome. Very little is known about transcriptional control in hypotrichs; identification of conserved non-coding sequences of orthologous hypotrich genes promises to provide clues to potential cis-acting control signals.

Constancy of adenine methylation in Tetrahymena macronuclear DNA

Macronuclear DNA from Tetrahymena was examined in order to determine whether the pattern of adenine methylation changed with the transcriptional activity of nearby genes. The DNA from growing, starved and conjugating cells was digested with six restriction enzymes which are sensitive to methylation of adenine within their recognition site. Southern blots of the restricted DNAs were probed with seven cDNA clones and one genomic clone which are homologous to polyA+ RNAs, whose transcriptional activity varies with the physiological state of the cell. One of the cDNA clones, BC11, had not been described previously. It hybridized to a 1.3 kb transcript which was present in populations of starved and conjugating, but not in growing cells. On Southern blots of genomic DNA it hybridized to a complex pattern of bands which was highly polymorphic between the DNAs of closely related strains. It was estimated that between 137 and 272 sites were assayed for changes in methylation, including at least 27 sites which were known to be methylated. No differences were seen in the size of restriction fragments from cells in different physiological states. The data suggested that the methylation pattern, which is determined during macronuclear development, does not vary with the physiological state of the cell.

Actin, tubulin and H4 histone genes in three species of hypotrichous ciliated protozoa

In hypotrichous ciliated protozoa, genes are transcribed in the macronucleus where the genome consists of 'gene-sized' linear DNA molecules. We have isolated clones of actin, tubulin and H4 histone macronuclear genes from Oxytricha nova, Stylonychia lemnae and Euplotes crassus in an effort to determine if they possess molecules of similar size for a given coding function, and also to determine the size range of non-coding DNA present on these molecules. Our results indicate that while the length of their non-coding DNA can vary slightly, both between different hypotrichs and within the gene family of a single organism, actin and tubulin macronuclear molecules are similarly sized. The sizes observed for these molecules support the hypothesis that each macronuclear molecule encodes a single gene. However, the H4 histone macronuclear molecules show a much wider size range and generally are much longer than necessary to encode the H4 histone. We therefore sequenced a 1700-bp H4 histone macronuclear molecule from O. nova to determine if it might possibly encode additional gene products. Sequence data reveals the presence of nine open reading frames (ORFs) greater than 100 bp in length; however, Northern hybridization analysis of the products of this DNA molecule reveals only a single transcript.

Organization of the Euplotes crassus micronuclear genome

Euplotes crassus, like other hypotrichous ciliated protozoa, eliminates most of its micronuclear chromosomal DNA in the process of forming the small linear DNA molecules that comprise the macronuclear genome. By characterizing randomly selected lambda phage clones of E. crassus micronuclear DNA, we have determined the distribution of repetitive and unique sequences and the arrangement of macronuclear genes relative to eliminated DNA. This allows us to compare the E. crassus micronuclear genome organization to that of another distantly related hypotrichous ciliate, Oxytricha nova. The clones from E. crassus segregate into three prevalent classes: those containing primarily eliminated repetitive DNA (Class I); those containing macronuclear genes in addition to repetitive sequences (Class II); and those containing only eliminated unique sequence DNA (Class III). All of the repetitive sequences in these clones belong to the same highly abundant repetitive element family. Our results demonstrate that the sequence organization of the E. crassus and O. nova micronuclear genomes is related in that the macronuclear genes are clustered together in the micronuclear genome and the eliminated unique sequences occur in long stretches that are uninterrupted by repetitive sequences. In both organisms a single repetitive element family comprises the majority of the eliminated interspersed middle repetitive DNA and appears to be preferentially associated with the macronuclear sequence clusters. The similarities in the sequence organization in these two organisms suggest that clustering of macronuclear genes plays a role in the chromosome fragmentation process.

Multiple sequence versions of the Oxytricha fallax 81-MAC alternate processing family

The 81-MAC family consists of three sizes of macronuclear chromosomes in Oxytricha fallax. Clones of these and of micronuclear homologs have been classified according to DNA sequence into three highly homologous (95.9-97.9%), but distinct versions. Version A is represented by a micronuclear clone and by clones of two different-sized macronuclear chromosomes, showing that alternate processing of micronuclear DNA is responsible for the variety of sizes of macronuclear chromosomes. Three Internal Eliminated Sequences (IES's) are demonstrated in Version A micronuclear DNA. Two have been sequenced and show short, flanking direct repeats but no inverted terminal repeats. Version C micronuclear DNA has interruptions in the macronuclear homology which correspond closely to the Version A IES's. Whether they are true IES's is unknown because no Version C macronuclear DNA has been demonstrated. Version C micronuclear DNA may be "macronuclear-homologous" but "micronucleus-limited" and not "macronucleus-destined." Version B is represented by macronuclear DNA clones, but no micronuclear clones. Vegetative micronuclear aneuploidy is suggested. The possible role of micronuclear defects in somatic karyonidal senescence is discussed in light of the precise macronuclear chromosome copy controls demonstrated within the 81-MAC family. These controls apparently operate throughout karyonidal life to maintain 1) a constant absolute amount of 81-MAC sequences in the macronucleus and 2) a constant stoichiometry within the family, both according to version and chromosome size.

Alternative processing of sequences during macronuclear development in Tetrahymena thermophila

DNA is eliminated during development of the somatic MACronucleus from the germinal MICronucleus in the ciliated protozoan, Tetrahymena thermophila. Facultatively persistent sequences are a class of sequences that persist in the MAC DNA of some cell lines but are eliminated from the MAC DNA of other cell lines. One cloned MAC fragment contains a persistent sequence as well as sequences normally retained in the MAC. When this cloned fragment was used to construct MAC restriction maps of this region in cell lines whose MAC DNAs do, or do not, contain the persistent sequence, extensive variation in the map flanking this region was observed. The different DNA rearrangements of this MIC segment are epigenetically determined during or soon after MAC development. Moreover, different rearrangements may occur among the 45 copies of this MIC segment as a MAC is formed, resulting in polymorphisms that are later resolved by phenotypic assortment.

Mobile elements bounded by C4A4 telomeric repeats in Oxytricha fallax

A novel family of micronuclear elements termed telomere-bearing elements (TBEs) is described. All 1900 family members are eliminated during macronuclear development. We conclude that they are transposons, first because the members are moderately conserved in sequence and probably dispersed in the genome. Second, in two cases, sequence comparison of the termini and flanks of the element with the corresponding empty site indicate that elements cause 3 bp target duplications (AAT) upon insertion; the 3 bp are part of the 5 bp target sequence, AATGA. Lastly, both elements carry 77 or 78 bp inverted terminal repeats. The tip of each inverted terminal repeat is the 17 bp telomere-like sequence 5' C1A4C4A4C4. At least half of the elements have these 17 bp or an extremely similar sequence. One possible pathway for transposition into new micronuclear sites starts in the developing macronucleus with excision to create a free linear form to which telomeres are added, followed by a low frequency of movement to the micronucleus, and insertion into the germ-line micronuclear DNA.