Genome architecture drives protein evolution in ciliates

Studies of microbial eukaryotes have been pivotal in the discovery of biological phenomena, including RNA editing, self-splicing RNA, and telomere addition. Here we extend this list by demonstrating that genome architecture, namely the extensive processing of somatic (macronuclear) genomes in some ciliate lineages, is associated with elevated rates of protein evolution. Using newly developed likelihood-based procedures for studying molecular evolution, we investigate 6 genes to compare 1) ciliate protein evolution to that of 3 other clades of eukaryotes (plants, animals, and fungi) and 2) protein evolution in ciliates with extensively processed macronuclear genomes to that of other ciliate lineages. In 5 of the 6 genes, ciliates are estimated to have a higher ratio of nonsynonymous/synonymous substitution rates, consistent with an increase in the rate of protein diversification in ciliates relative to other eukaryotes. Even more striking, there is a significant effect of genome architecture within ciliates as the most divergent proteins are consistently found in those lineages with the most highly processed macronuclear genomes. We propose a model whereby genome architecture-specifically chromosomal processing, amitosis within macronuclei, and epigenetics-allows ciliates to explore protein space in a novel manner. Further, we predict that examination of diverse eukaryotes will reveal additional evidence of the impact of genome architecture on molecular evolution.     

Chromosome Copy Number Variation and Control in the Ciliate Chilodonella uncinata

Copy number variations are widespread in eukaryotes. The unusual genome architecture of ciliates, in particular, with its process of amitosis in macronuclear division, provides a valuable model in which to study copy number variation. The current model of amitosis envisions stochastic distribution of macronuclear chromosomes during asexual reproduction. This suggests that amitosis is likely to result in high levels of copy number variation in ciliates, as dividing daughter cells can have variable copy numbers of chromosomes if chromosomal distribution during amitosis is a stochastic process. We examined chromosomal distribution during amitosis in Chilodonella uncinata, a ciliate with gene-size macronuclear chromosomes. We quantified 4 chromosomes in evolving populations of C. uncinata and modeled the amitotic distribution process. We found that macronuclear chromosomes differ in copy number from one another but that copy number does not change as expected under a stochastic process. The chromosome carrying SSU increased in copy number, which is consistent with selection to increase abundance; however, two other studied chromosomes displayed much lower than expected among-line variance. Our models suggest that balancing selection is sufficient to explain the observed maintenance of chromosome copy during asexual reproduction.     

A Simple and Rapid PCR-based Method to Isolate Complete Small Macronuclear Minichromosomes from Hypotrich Ciliates: 5S rDNA and S26 Ribosomal Protein Gene of Oxytricha (Sterkiella) nova

Hypotrich ciliates present a macronuclear genome consisting of gene-sized instead of chromosome-sized DNA molecules. Exploiting this unique eukaryotic genome feature, we introduce, for the first time in ciliates, a rapid and easy PCR method using telomeric primers to isolate small complete macronuclear DNA molecules or minichromosomes. Two presumably abundant macronuclear DNA molecules, containing ribosomal genes, were amplified from the Oxytricha (Sterkiella) nova complete genome after using this method, and then were cloned and sequenced. The 5S rDNA sequence of O. (S.) nova is the third one reported among hypotrich ciliates; its primary and secondary structure is compared with other eukaryotic 5S rRNAs. The ribosomal protein S26 gene is the first one reported among ciliates. This “End-End-PCR” method might be useful to obtain similar gene-sized macronuclear molecules from other hypotrich ciliates, and, therefore, to increase our knowledge on ribosomal genes in these eukaryotic microorganisms.     

Macronuclear gene-sized molecules of hypotrichs.

The macronuclear genome of hypotrichous ciliates consists of DNA molecules of gene-sized length. A macronuclear DNA molecule contains a single coding region. We have analyzed the many hypotrich macronuclear DNA sequences sequenced by us and others. No highly conserved promoter sequences nor replication initiation sequences have been identified in the 5' nor in the 3' non-translated regions, suggesting that promoter function in hypotrichs may differ from other eukaryotes. The macronuclear genes are intron-poor; approximately 19% of the genes sequenced to date have one to three introns. Not all macronuclear DNA molecules may be transcribed; some macronuclear molecules may not have any coding function. Codon bias in hypotrichs is different in many respects from other ciliates and from other eukaryotes.     

Multiple genes of apparent algal origin suggest ciliates may once have been photosynthetic

Plantae (as defined by Cavalier-Smith, 1981) plastids evolved via primary endosymbiosis whereby a heterotrophic protist enslaved a photosynthetic cyanobacterium. This "primary" plastid spread into other eukaryotes via secondary endosymbiosis. An important but contentious theory in algal evolution is the chromalveolate hypothesis that posits chromists (cryptophytes, haptophytes, and stramenopiles) and alveolates (ciliates, apicomplexans, and dinoflagellates) share a common ancestor that contained a red-algal-derived "secondary" plastid. Under this view, the existence of several later-diverging plastid-lacking chromalveolates such as ciliates and oomycetes would be explained by plastid loss in these lineages. To test the idea of a photosynthetic ancestry for ciliates, we used the 27,446 predicted proteins from the macronuclear genome of Tetrahymena thermophila to query prokaryotic and eukaryotic genomes. We identified 16 proteins of possible algal origin in the ciliates Tetrahymena and Paramecium tetraurelia. Fourteen of these are present in other chromalveolates. Here we compare and contrast the likely scenarios for algal-gene origin in ciliates either via multiple rounds of horizontal gene transfer (HGT) from algal prey or symbionts, or through endosymbiotic gene transfer (EGT) during a putative photosynthetic phase in their evolution.     

Chromosome organization and gene expansion in the highly fragmented genome of the ciliate Strombidium stylifer

Chromosomes are well-organized carriers of genetic information in eukaryotes and are usually quite long, carrying hundreds and thousands of genes. Intriguingly, a clade of single-celled ciliates, Spirotrichea, feature nanochromosomes—also called “gene-sized chromosomes”. These chromosomes predominantly carry only one gene, flanked by short telomere sequences. However, the organization and copy number variation of the chromosomes in these highly fragmented genomes remain unexplored in many groups of Spirotrichea, including the marine Strombidium. Using deep genome sequencing, we assembled the macronuclear genome of Strombidium stylifer into more than 18,000 nanochromosomes (~2.4 Kb long on average). Our results show that S. stylifer occupies an intermediate position during the evolutionary history of Strombidium lineage and experienced significant expansions in several gene families related to guanyl ribonucleotide binding. Based on the nucleotide distribution bias analysis and conserved motifs search in non-genic regions, we found that the subtelomeric regions have a conserved adenine-thymine (AT)-rich sequence motif. We also found that the copy number of nanochromosomes lacks precise regulation. This work sheds light on the unique features of chromosome structure in eukaryotes with highly fragmented genomes and reveals that a rather specialized evolutionary strategy at the genomic level has resulted in great diversity within the ciliated lineages.     

The Oxytricha trifallax Macronuclear Genome: A Complex Eukaryotic Genome with 16,000 Tiny 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.     

A simple and rapid PCR-based method to isolate complete small macronuclear minichromosomes from hypotrich ciliates: 5S rDNA and S26 ribosomal protein gene of Oxytricha (Sterkiella) nova

Hypotrich ciliates present a macronuclear genome consisting of gene-sized instead of chromosome-sized DNA molecules. Exploiting this unique eukaryotic genome feature, we introduce, for the first time in ciliates, a rapid and easy PCR method using telomeric primers to isolate small complete macronuclear DNA molecules or minichromosomes. Two presumably abundant macronuclear DNA molecules, containing ribosomal genes, were amplified from the Oxytricha (Sterkiella) nova complete genome after using this method, and then were cloned and sequenced. The 5S rDNA sequence of O. (S.) nova is the third one reported among hypotrich ciliates; its primary and secondary structure is compared with other eukaryotic 5S rRNAs. The ribosomal protein S26 gene is the first one reported among ciliates. This "End-End-PCR" method might be useful to obtain similar gene-sized macronuclear molecules from other hypotrich ciliates, and, therefore, to increase our knowledge on ribosomal genes in these eukaryotic microorganisms.     

Variation in macronuclear genome content of three ciliates with extensive chromosomal fragmentation: a preliminary analysis

The genome architecture of ciliates, including features such as nuclear dualism and large-scale genome rearrangements, impacts gene and genome evolution in these organisms. To better understand the structure of macronuclear chromosomes in ciliates with extensively processed chromosomes, a sample of complete macronuclear chromosomes was sequenced from three ciliate species: Metopus es (Class [Cl]: Armophorea), Nyctotherus ovalis (Cl: Armophorea), and Chilodonella uncinata (Cl: Phyllopharyngea). By cloning whole macronuclear chromosomes into a plasmid vector, we generated nine clones from each of M. es and C. uncinata, and 37 clones from N. ovalis. Analysis of these macronuclear chromosomes provides insight into the evolution of genome features such as chromosome content, gene structure, and genetic code. Phylogenetic patterns can be found in telomere structure and codon usage, which are both more similar in M. es and N. ovalis than C. uncinata. In addition, we provide evidence of lateral transfer of a bacterial endo-beta-mannanase gene onto a M. es chromosome and report the discovery of a 42-bp conserved sequence motif within N. ovalis untranslated regions.     

Evolution of Germline-Limited Sequences in Two Populations of the Ciliate Chilodonella uncinata

Ciliates are microbial eukaryotes that separate their nuclear functions into a germline micronucleus and a somatic macronucleus. During development of the macronucleus the genome undergoes a series of reorganization events that includes the precise excision of intervening DNA. Here, we determine the architecture of four loci in the micronuclear and macronuclear genomes of the ciliate Chilodonella uncinata and compare the levels of variation in micronuclear-limited sequences to macronuclear destined sequences at two of these loci. We find that within a population, germline-limited sequences are evolving at the same rate as other putatively neutral sites, but between populations germline-limited sequences are accumulating mutations at a much faster rate than other sites. We also find evidence of macronuclear recombination and incomplete elimination of intervening DNA, which result in increased diversity in the macronuclear genome. Our results support the assertion that the unusual genomic features of ciliates can result in rapid and unpredicted patterns of diversification.     

Genetic and evolutionary consequences of symbiosis

Permanent stable symbioses, primarily microbial, are analyzed as parasexual phenomena from the evolutionary point of view. Such associations bring together in single individuals heritable traits of high selective advantage in certain environments. By convergent evolution several types of associations have repeatedly arisen: motile photosynthetic forms, nitrogen fixing and wood digesting complexes, and so forth. Many examples are discussed from the point of view of the number of originally independent genomes that comprise the recognizable individuals.The level of partner integration in many associations is analyzed. Examples of many levels: genic, gene product, metabolite, behavioral, and the methods by which they can be distinguished are discussed.The literature concerning a large number of associations is reviewed: Methanobacillus; predatory and consortia bacteria; blue green algal sheath dwelling bacteria; anaerobic worm-bacterial; algal, and foreign chloroplast retention by heterotrophic eukaryotes (ciliates, coelenterates, mollusks); the double nucleated photosynthetic dinoflagellate (Peridinium balticum); hindgut microbes in termites and woodroaches (Pyrsonympha, Barbulanympha and their associated spirochetes and other bacteria); sand dwelling and other ciliates and their associated bacteria; and so forth. The status of observations and artificial systems claiming evidence for transfer of genes between very distantly related organisms is critically discussed.A continuum from nearly completely autonomous partners (e.g., zoochlorellae in invertebrate animals) to nearly unrecognizable merged components (e.g., gamma particles in Blastocladiella) is found to exist among examples of extant organisms. The diversity and prevalence of such associations support the concept that there are many precedents for the steps hypothesized in the serial endosymbiotic theory of the origin of eukaryotic cells.     

Packaging of Subchromosomal-Size DNA Molecules in Chromatin Bodies in the Ciliate Macronucleus

Molecular Biology - A fundamental difference between somatic nuclei (macronuclei) of ciliates and cell nuclei of higher eukaryotes is that the macronuclear genome is a huge number (up to tens or...     

The Kinetic and Analytical Complexities of the DNA Genomes of Certain Marine and Fresh-Water Ciliates1

We investigated the macronuclear DNA genomes of several marine and fresh-water ciliates. The marine forms studied were: Uronema nigricans, Parauronema virginianum, Parauronema acutum, and two strains of Miamiensis avidus; the fresh-water ciliates included: Tetrahymena pyriformis, Paramecium octaurelia, and P. caudatum. The organisms were cultured axenically and the DNA extracted from isolated and purified macronuclear preparations. Reassociation rate constants of purified DNA preparations used to calculate kinetic complexity were determined both optically and by hydroxyapatite chromatography. Analytical complexity was determined chemically. Ciliate macronuclear DNA appeared to reassociate as a single unique sequence, except for a small fraction (4% of the total DNA) that was repetitive and renatured rapidly. Values for the kinetic complexities of macronuclear DNA in these forms varied over a relatively narrow range, from 1.5 to 3.8 times 10¹⁰ daltons, and were only 7–15x larger than that of the bacterium Escherichia coli. On the other hand, values for analytical complexities of macronuclear DNA of marine and fresh-water ciliates varied over two orders of magnitude and were related to the size of the animals. It is suggested that ploidy levels of macronuclear DNA in these ciliates may represent a functionally permanent amplification of the genome.     

Macronuclear DNA organization and transcription in Paramecium primaurelia

Macronuclear DNA was isolated from Paramecium primaurelia, stock 168. Although the macronucleus is polyploid to the extent of 840C, in other respect the DNA appears to be simply organized, having neither satellite sequences nor substantial amounts of intermediately repetitive sequence. The sequence complexity of macronuclear DNA is quite low for a eukaryote cell, being approximately 19 times more complex than the genome of Escherichia coli. In addition, the GC content is low (25%) and the isolated DNA molecules have lengths mostly in the range 0.2–5 μm. In these various respects, the macronuclear DNA of Paramecium is similar to that of other ciliates. A clone of Paramecium cultured under controlled conditions contains polyadenylated RNA sequences which are homologous to 5–8% of the macronuclear DNA. Sequence complexity analysis indicates that the polyadenylated RNA contains two abundance classes of molecules, one present at low frequency and transcribed from approximately 10⁴ genes, the other at 100 times greater concentration and transcribed from about 100 genes. The relevance of these results to the control of gene expression in Paramecium is discussed.     

Evolution of the Germline Actin Gene in Hypotrichous Ciliates: Multiple Nonscrambled IESs at Extremely Conserved Locations in Two Urostylids

In hypotrichous ciliates, macronuclear chromosomes are gene‐sized, and micronuclear genes contain short, noncoding internal eliminated segments (IESs) as well as macronuclear‐destined segments (MDSs). In the present study, we characterized the complete macronuclear gene and two to three types of micronuclear actin genes of two urostylid species, i.e. Pseudokeronopsis rubra and Uroleptopsis citrina. Our results show that (1) the gain/loss of IES happens frequently in the subclass Hypotrichia (formerly Stichotrichia), and high fragmentation of germline genes does not imply for gene scrambling; and (2) the micronuclear actin gene is scrambled in the order Sporadotrichida but nonscrambled in the orders Urostylida and Stichotrichida, indicating the independent evolution of MIC‐actin gene patterns in different orders of hypotrichs; (3) locations of MDS–IES junctions of micronuclear actin gene in coding regions are conserved among closely related species.     

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.