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.     

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.     

Recombination and Assortment in the Macronucleus of TETRAHYMENA THERMOPHILA: A Theoretical Study by Computer Simulation

The compound nature of the macronucleus of Tetrahymena thermophila presents multiple opportunities for recombination between genes on the same macronuclear chromosome. Such recombinants should be detectable through their assortment at subsequent amitotic macronuclear divisions. Thus, a macronucleus that is initially AB/ab should produce recombinant assortees of the genotypes Ab/aB. Computer simulation shows that, when the recombination frequency is two or fewer times per cell cycle, recombinant assortees are produced at experimentally measurable frequencies of less than 40%. At higher recombination frequencies, linked genes appear to assort independently. The simulations also show that recombination during macronuclear development can be distinguished from recombination in subsequent cell cycles only if the first appearance of recombinant assortees is 100 or more fissions after conjugation. The use of macronuclear recombination and assortment as a means of mapping macronuclear genes is severely constrained by the large variances in assortment outcomes; with experimentally small sample sizes, such mapping is impossible.     

Organization of the Euplotes crassus Micronuclear Genome1

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.     

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.     

Conserved linkage of two genes on the same macronuclear chromosome in spirotrichous ciliates

Macronuclear chromosomes of spirotrichous ciliates are mainly “nanochromosomes” containing only a single gene. We identified a two-gene chromosome in the spirotrich Sterkiella histriomuscorum (formerly Oxytricha trifallax) which, unlike other characterized two-gene molecules, contains reading frames oriented tail to tail. These are homologs of ribosomal protein L29 (RPL29) and cyclophilin. We found that both genes are transcribed, with their polyadenylation sites on opposite strands separated by only 135 bp. Furthermore, both genes in S. histriomuscorum are present only on one macronuclear chromosome and do not occur alone or linked to other genes. The corresponding micronuclear locus is fragmented into three nonscrambled gene segments (MDSs), separated by two noncoding segments (IESs). We also found that these two genes are linked on a macronuclear chromosome, similarly arranged tail to tail, in the three spirotrichs Stylonychia lemnae, Uroleptus sp., and Holosticha sp.. In addition, single-gene macronuclear chromosomes containing only the RPL29 gene were detected in the earlier diverged Holosticha and Uroleptus. These observations suggest a possible evolutionary trend towards loss of chromosomal breakage between these two genes. This study is the first to examine gene linkage in the macronucleus of several spirotrichs and may provide insight into the evolution of multi-gene macronuclear chromosomes and chromosomal fragmentation in spirotrichs. Electronic Supplementary Material Supplementary material is available for this article at     

Does Ribosomal DNA Get Out of the Micronuclear Chromosome in Paramecium tetraurelia by Means of a Rolling Circle?

The macronuclear genes coding for rRNA (ribosomal DNA [rDNA]) of Paramecium tetraurelia, stock 51, are arranged in polymers consisting of units made up of a transcribed coding region and a nontranscribed spacer region. The whole macronuclear polymer ends with a portion of the spacer on either end followed by a telomere. Six kinds of macronuclear units, or genes, were mapped. Spacers were different, and transcribed regions were the same. These genes are found in markedly different numbers in the macronucleus. The most common gene shows two regions in the spacer where a sequence is followed by a direct repeat. The next most common gene is similar but shows a deletion plus a number of base pair substitutions. Although most cosmid clones contain only a single kind of gene, many contain more than one. These are thought to be produced by somatic crossing over. The four micronuclear genes that have been isolated consist of a single central transcribed region and portions of the spacer on either end. Sequencing indicates that the two ends of the molecule are partially redundant. While the spacer region at the right end of the macronuclear polymer is derived from the micronuclear spacer on the right, the spacer at the left end of the macronuclear polymer is derived from regions of the micronuclear spacer on both the right and the left. To account for this situation, a rolling-circle model for generation of the macronuclear rDNA from the micronuclear DNA is proposed.     

Studies on macronuclear DNA from Paramecium aurelia

Macronuclear DNA was isolated from purified macronuclei of Paramecium aurelia and the size distribution was determined with regard to growth phase and method of extraction. DNA molecules as long as 105 microns and as short as 0.2 microns were observed. It was concluded that the method of extraction affected the observed length of DNA extracted and that macronuclear DNA isolated from cells in balanced growth was less susceptible to nuclease degradation than was DNA isolated from cells in stationary phase. Renaturation studies were performed on macronuclear DNA and a kinetic complexity of 22-times E. coli DNA was determined. This value was similar to those values reported for Tetrahymena and Stylonychia macronuclear DNA. Correcting for GC base content yielded a kinetic complexity for Paramecium macronuclear DNA of 11-times E. coli DNA which corresponded to 3 X 10(10) daltons. There would be about 1400 copies of a unit genome of this complexity within each newly replicated macronucleus. Density gradient analysis indicated that the genes coding for ribosomal RNA had a greater density in CsCl than the bulk DNA. Molecular hybridization studies indicated that the genes coding for 25 S RNA represented 0.14 percent of the total macronuclear DNA. Correcting for GC base content, this corresponded to 30-35 25 S RNA genes per unit genome. These results on Paramecium are discussed in relationship to other ciliate macronuclear DNA.     

Automated image analysis for array hybridization experiments.

MOTIVATION: Image analysis is a major part of data evaluation for array hybridization experiments in molecular biology. The program presented here is designed to analyze automatically images from hybridization experiments with various arrangements: different kinds of probes (oligonucleotides or complex probes), different supports (nylon filters or glass slides), different labeling of probes (radioactively or fluorescently). The program is currently applied to oligonucleotide fingerprinting projects and complex hybridizations. The only precondition for the use of the program is that the targets are arrayed in a grid, which can be approximately transformed to an orthogonal equidistant grid by a projective mapping. RESULTS: We demonstrate that our program can cope with the following problems: global distortion of the grid, missing of grid nodes, local deviation of the spot from its specified grid position. This is checked by different quality measures. The image analysis of oligonucleotide fingerprint experiments on an entire genetic library is used, in clustering procedures, to group related clones together. The results show that the program yields automatically generated high quality input data for follow up analysis such as clustering procedures. AVAILABILITY: The executable files will be available upon request for academics.     

Organization of gene and non-gene sequences in micronuclear DNA of Oxytricha nova.

In order to study the derivation of the macronuclear genome from the micronuclear genome in Oxytricha nova micronuclear DNA was partially digested with EcoRI, size fractionated, and then cloned in the lambda phage Charon 8. Clones were selected a) at random b) by hybridization with macronuclear DNA or c) by hybridization with clones of macronuclear DNA. One group of these clones contains only unique sequence DNA, and all of these had sequences that were homologous to macronuclear sequences. The number of macronuclear genes with sequences homologous to these micronuclear clones indicates that macronuclear sequences are clustered in the micronuclear genome. Many micronuclear clones contain repetitive DNA sequences and hybridize to numerous EcoRI fragments of total micronuclear DNA, yielding similar but non-identical patterns. Some micronuclear clones containing these repetitive sequences also contained unique sequence DNA that hybridized to a macronuclear sequence. These clones define a major interspersed repetitive sequence family in the micronuclear genome that is eliminated during formation of the macronuclear genome.     

Variable copy number of macronuclear DNA molecules in Tetrahymena.

In Tetrahymena, the DNA of the macronucleus exists as very large (100 to 4,000-kb) linear molecules that are randomly partitioned to the daughter cells during cell division. This genetic system leads directly to an assortment of alleles such that all loci become homozygous during vegetative growth. Apparently, there is a copy number control mechanism operative that adjusts the number of each macronuclear DNA molecule so that macronuclear DNA molecules (with their loci) are not lost and aneuploid death is a rare event. In comparing Southern analyses of the DNA from various species of Tetrahymena using histone H4 genes as a probe, we find different band intensities in many species. These differences in band intensities primarily reflect differences in the copy number of macronuclear DNA molecules. The variation in copy number of macronuclear DNA molecules in some species is greater than an order of magnitude. These observations are consistent with a developmental control mechanism that operates by increasing the macronuclear copy number of specific DNA molecules (and the genes located on these molecules) to provide the relatively high gene copy number required for highly expressed proteins.     

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.     

EC_oligos: automated and whole-genome primer design for exons within one or between two genomes

SUMMARY: EC_oligos designs oligonucleotides (oligos) from exons of annotated genomic sequence information. It can automatically and rapidly select oligos that are conserved between two sets of sequence data, and can pair up oligos for use as PCR primers. It can do this on a whole-genome scale and according to user-defined criteria. AVAILABILITY: The source code, executable program and user manual are available at ftp://ftp.ebi.ac.uk/pub/software/dos/EC_oligos/.     

Intergenic sequence inspector: searching and identifying bacterial RNAs

Intergenic Sequence Inspector (ISI) is a program that helps in identifying bacterial ribonucleic acids (RNAs).It automatically extracts, selects and visualizes candidate intergenic regions (IGRs) that bear conservations between phylogenetically related species, displaying sequence and structural signatures of RNA genes encompassing putative promoters, terminators, RNA secondary structure predictions and the G+C percent of the selected sequences. ISI is intentionally modular to analyze various bacterial genomes according to their intrinsic specificities and to the available data. ISI provides a sum of information to the user who can evaluate whether or not an IGR is susceptible to express RNAs.     

DNA elimination and its relation to quantities in the macronucleus of Tetrahymena

The macronucleus of Tetrahymena contains a large number of DNA molecules of subchromosomal size. They belong to about 270 species each one occurring at an average number of 45 copies Macronuclei divide unequally and nothing is known of segregation control. This and the elimination and degradation of DNA during macronuclear amitosis make the clonal stability of macronuclei a problem of qualitative and quantitative control on a subchromosomal level. We studied the contribution of DNA elimination to the quantitative composition of the macronucleus cytophotometrically in single cells of different strains. This was done under standard conditions and under conditions known to influence the amount of macronuclear DNA. The following results were found: Elimination of DNA occurs at almost every division. The size of the elimination body is highly variable but still positively correlated with the macronuclear DNA content. In T. thermophila the amount of eliminated DNA is 2.5% of the G2 content and is not dependent on the growth state. It varies with species, amounting to as much as 8% in T pigmentosa. During conditions which increase the macronuclear DNA content, very little DNA is eliminate. On the other hand, large amounts are eliminated under other conditions causing the macronuclear DNA content to decrease. DNA to be eliminated at division is synthesized at the same time as bulk DNA. We developed a computer program which helps us study the effects of DNA elimination and unequal divisions upon the copy numbers of subchromosomal DNA classes. The result indicates that in a given cell line at least one of the DNA molecules becoms extinct after 60 generations which we expect would cause the cell's extinction and restrict a clone's life to 60 generations. As this does not happen in nature, there must be some control of the copy numbers preventing their extinction during vegetative multiplication. Whether elimination increases or decreases the imbalance of genes remains to be investigated. © 1992 Wiley-Liss, Inc.     

Automatic clustering of orthologs and inparalogs shared by multiple proteomes

MOTIVATION: The complete sequencing of many genomes has made it possible to identify orthologous genes descending from a common ancestor. However, reconstruction of evolutionary history over long time periods faces many challenges due to gene duplications and losses. Identification of orthologous groups shared by multiple proteomes therefore becomes a clustering problem in which an optimal compromise between conflicting evidences needs to be found. RESULTS: Here we present a new proteome-scale analysis program called MultiParanoid that can automatically find orthology relationships between proteins in multiple proteomes. The software is an extension of the InParanoid program that identifies orthologs and inparalogs in pairwise proteome comparisons. MultiParanoid applies a clustering algorithm to merge multiple pairwise ortholog groups from InParanoid into multi-species ortholog groups. To avoid outparalogs in the same cluster, MultiParanoid only combines species that share the same last ancestor. To validate the clustering technique, we compared the results to a reference set obtained by manual phylogenetic analysis. We further compared the results to ortholog groups in KOGs and OrthoMCL, which revealed that MultiParanoid produces substantially fewer outparalogs than these resources. AVAILABILITY: MultiParanoid is a freely available standalone program that enables efficient orthology analysis much needed in the post-genomic era. A web-based service providing access to the original datasets, the resulting groups of orthologs, and the source code of the program can be found at http://multiparanoid.cgb.ki.se.