- Research article
- Open Access
Mutation accumulation in Tetrahymena
© Brito et al; licensee BioMed Central Ltd. 2010
- Received: 21 May 2010
- Accepted: 15 November 2010
- Published: 15 November 2010
The rate and fitness effects of mutations are key in understanding the evolution of every species. Traditionally, these parameters are estimated in mutation accumulation experiments where replicate lines are propagated in conditions that allow mutations to randomly accumulate without the purging effect of natural selection. These experiments have been performed with many model organisms but we still lack empirical estimates of the rate and effects of mutation in the protists.
We performed a mutation accumulation (MA) experiment in Tetrahymena thermophila, a species that can reproduce sexually and asexually in nature, and measured both the mean decline and variance increase in fitness of 20 lines. The results obtained with T. thermophila were compared with T. pyriformis that is an obligate asexual species. We show that MA lines of T. thermophila go to extinction at a rate of 1.25 clonal extinctions per bottleneck. In contrast, populations of T. pyriformis show a much higher resistance to extinction. Variation in gene copy number is likely to be a key factor in explaining these results, and indeed we show that T. pyriformis has a higher mean copy number per cell than T. thermophila. From fitness measurements during the MA experiment, we infer a rate of mutation to copy number variation of 0.0333 per haploid MAC genome of T. thermophila and a mean effect against copy number variation of 0.16. A strong effect of population size in the rate of fitness decline was also found, consistent with the increased power of natural selection.
The rate of clonal extinction measured for T. thermophila is characteristic of a mutational degradation and suggests that this species must undergo sexual reproduction to avoid the deleterious effects detected in the laboratory experiments. We also suggest that an increase in chromosomal copy number associated with the phenotypic assortment of amitotic divisions can provide an alternative mechanism to escape the deleterious effect of random chromosomal copy number variation in species like T. pyriformis that lack the resetting mechanism of sexual reproduction. Our results are relevant to the understanding of cell line longevity and senescence in ciliates.
- Copy Number Variation
- Deleterious Mutation
- Mutation Accumulation
- Ferric Sodium
- Chromosomal Copy Number
Understanding the rate at which new mutations arise and their effects on fitness has been central in the development of many important theories in evolutionary genetics such as those that attempt to explain the evolution of sex and recombination , patterns of DNA sequence variation and molecular evolution in chromosomal regions with different recombination rates , the evolution of diploidy  and polyploidy , and the risk of small populations to become extinct . Mutation accumulation (MA) is a classical experiment to estimate the rate at which deleterious mutations arise in populations and their effects on fitness [recently reviewed in ]. This approach is characterized by the propagation of replicates of highly inbreed and clonal populations over a long period of time in an environment where mutations over the whole genome are allowed to randomly accumulate without the purging effect of selection. Because different replicate lines will accumulate distinct mutations, the variance in fitness among lines is expected to increase over time. Fitness assays carried out with contemporaneous and ancestral populations may thus explain the phenotypic effects of the new accumulated mutations and can be used to provide estimates of the haploid genomic mutation rate per generation (U), the mean mutational effect (E(s)) [7, 8], and the distribution of effects of new mutations (f(s)) . Several model organisms, including some microbes, have been subject to such MA experiments and have permitted estimates of the rates and effects of mutations. MA experiments are expected to increase the genetic load and potentially lead to extinction due to fitness decline , as it has been observed, for example, in strains of E. coli  or in yeast .
Mutation accumulation studies have already been done with a wide variety of organisms, including viruses , bacteria [14, 15], fungi , plants , and animals [e.g. [7, 18]]. To our knowledge, however, experiments to estimate the rate and effects of mutations have never been performed with any of the diverse lineages of protists. Although Lynch and Gabriel  had previously estimated the mutation load of Paramecium caudatum this estimate was performed on data collected to determine the species maximum lifespan and therefore replicate lines were obtained from a diverse population and initial phenotypic variability was not set to zero as it is typical of MA studies.
Ciliates have a unique form of genomic information. As typical ciliate, Tetrahymena thermophila presents a marked nuclear dimorphism characterized by a diploid micronucleus (Mic) with five chromosomes that functions primarily as the germ line, and a somatic polyploid macronucleus (Mac) thought to be responsible for all gene expression during vegetative growth [20–23]. The Mac originates from the Mic and goes through several developmentally programmed DNA rearrangements such as site-specific fragmentation of the five chromosomes leading to 225 autonomous replicating pieces (ARP) that comprises a total of 104 Mb in length and more than 27000 predicted protein-coding genes [24–26]. Each of these ARP undergoes several rounds of replication until each locus is present at an average level of 45 copies per Mac genome [27, 28].
Another important consequence of the random segregation of the homologous chromosome copies at each cell division is the production of homozygous phenotypes from sexually-produced heterozygotes by the process of phenotypic assortment . The general explanation for phenotypic assortment simply requires a randomly distribution of alleles between the two daughter cells [28, 36]. The rate of assortment is dependent on the initial input ratio and on the number of assorting units and proceeds to homozygosity independently of the dominance relationships between alleles. The allele assortment rate has been estimated in T. thermophila for several loci and leads to a value of about 45 assorting units typical of the Mac of this species [27, 28, 36].
In this study we carry out a mutation accumulation experiment with the ciliate protozoan Tetrahymena thermophila with the aim to determine its somatic macronuclear (Mac) genomic mutation rate and the effect of new deleterious mutations on the mean population fitness. Even though experiments similar to this MA have been performed in the past [32, 37, 38], mainly to determine the possible immortality of ciliates , this is the first attempt to estimate the genomic mutation rate and effects of deleterious mutations in a protist. Due to the particular mode that characterizes the genetic segregation in the macronucleous of Tetrahymena that defies the rules of diploid Mendelian inheritance, we developed a simple simulation model to estimate both U - the rate at which copy number variation in the MAC occurs- and s - the mean effect on fitness caused by deviations of the optimum copy number. The MA experiment carried out in this study with T. thermophila was replicated with cell lines of Tetrahymena pyriformis with the aim to compare their rates of clonal extinction. The macronucleous of T. pyriformis is also fragmented into many ARP, each one highly amplified, and cell divisions are characterized by the amitotic segregation of alleles between the two daughter cells [31, 40]. In contrast to T. thermophila, T. pyriformis organisms are amicronucleated, i.e. they lack the germinal micronucleus and are therefore unable to reproduce sexually with both germinal and somatic functions being carried out by the macronucleus alone. Most importantly, these organisms lack the ability to purge deleterious mutation through meiosis and conjugation, and therefore they provide an important contrast to the somatic MA occurring the in Mac of T. thermophila. Although during the course of the experiment T. thermophila accumulates mutations in both Mic and Mac nuclei, our phenotypic measures of fitness only detect the mutations accumulated in the macronucleus. Hence the results are comparable between species and they show that the mutation load in T. thermophila macronucleus is much greater than in T. pyriformis. Tetrahymena offers the opportunity to investigate (1) the establishment of the somatic vs. germ line during evolution process, (2) the advantage of having either sexual or asexual life cycles, and (3) the dynamics of and diploidy vs. polyploidy coexisting in the same cell/organism. Finally, T. thermophila is an important model organism in biomedical research and understanding the magnitude of its genomic mutation rate is critical to any evolutionary interpretation of its phenotypes.
Mutation accumulation experiment and fitness assays
Additional file 2: Video of an abnormal phenotype of Tetrahymena thermophila due to the MA experiment. Live images from T. thermophila cells with DNA stained with Hoechst 33342 and directly recorded with a Leica DMRA2 microscope. This video shows an example of a cell with an abnormal longitudinal dividing plan and an altered circular swimming behavior. (MOV 6 MB)
Relative proportion of T. pyriformis gene copy number per ng of DNA relative to T. thermophila
The rate at which T. thermophila clones become extinct through continuous reduction in population size imposed by the strong bottlenecks is characteristic of a mutational meltdown . During the course of the experiment, a synergistic interaction is created between population size and the random effect of genetic drift, and this promotes the build up of the population mutation load which, in its extreme, causes a mutation meltdown and drives populations to extinction . In fact, although some extinctions have been previously observed in mutator strains of bacteria  and yeasts , no such high rate of extinction was observed in MA experiments with other organisms. In our experiments T. thermophila showed a decline rate of 1.25 clonal extinctions per bottleneck, a value that is at least two and three orders of magnitude faster than the rates of 0.01 and 0.06 obtained in the yeast and bacteria experiments. It is unlikely that the fast rate observed in T. thermophila is due solely to the accumulation of random deleterious mutations, such as nucleotide changes, transpositions and/or small insertions-deletions as observed in other organisms. More plausible is that this rate of extinction reflects the process of gain and loss of chromosomes that is characteristic of the amitotic division of the macronucleous of T. thermophila. Theory predicts that the accumulation of deleterious mutations in MA lines is accompanied by a steady decrease in the mean population fitness . However, in our study we observed occasional recoveries of the mean population fitness (Figure 3 and Additional file 3: Fitness measurements for each population across bottlenecks), which can only be obtained when the rate of back mutations is sufficiently high for offspring to be fitter than the parental generation . In our simulations we mimic these recoveries when we assume no bias for the probability of gain and loss of chromosomal copies in the amitotic division, a condition that still lacks empirical confirmation.
Lynch and Gabriel  analyzed the data from Takagi and Yoshida  and, by regressing the fraction of extinct or non-dividing cells on the number of cell divisions estimated the macronuclear mutation load (μs) of 0.0009 for a clone of P. caudatum - a ciliate that shares with Tetrahymena many of the genetic features associated with the presence of two nuclei and the amitotic division of the macronucleus. There are several differences both in the experimental design used by Takagi and Yoshida and the model used by Lynch and Gabriel to reanalyzed the data from that experiment. Takagi and Yoshida followed the evolution of 10 clones, replicated 3 times each, derived from a large population of P caudatum. Each of these lines were propagated with daily transfers of one cell, and occasionally replaced by members of sister sublines, a procedure we have not done and is not typical in MA experiments. Lynch and Gabriel applied a model to the data obtained by measuring the fraction of clonal populations with dividing cells during the experiment. Lynch and Gabriel model assumes that mutations occur following a Poisson process, just as we do in our model, and that each mutation causes a negative effect on fitness, s, that is irreversible. This was not assumed in our model. In our model each mutation causes a deviation from optimal copy number (W(Z) = exp(-s|Z-Zm|), see simulation model in Methods section), which causes a negative fitness effect that can be reverted. Under Lynch and Gabriel model, no selection occurs in the experiment and a linear decline in log mean fitness is expected over time, which translates into an estimate of Us, the mutational load. Lynch and Gabriel point out that since several simplifying assumptions are done in the model, including lack of back-mutation and lack of selection, this estimated is likely downwardly biased. Importantly not only Takagi and Yoshida did not perform a classical MA experiment but also they did not measure competitive fitness. Lynch and Gabriel performed a regression of the fraction of clones with dividing cells on the number of cell divisions to get a rough estimate of the mutational load in P. caudatum of 0.0009. If we applied the same model used by Lynch and Gabriel on our data of the extinction pattern (as used for P. caudatum) we obtain an equivalent value for the macronuclear mutation load in T. thermophila of 0.004. This value is four times larger than the one estimated for P. caudatum, which may reflect differences in macronuclear gene content between these two species and differences in experimental design. This load measures the proportional reduction in fitness due to deleterious mutations and lies within the range of μs (0.02 to 0.0002) estimated for eukaryotes . Interestingly, a 0.002 rate of clonal death was independently estimated by Merriam and Bruns  for T. thermophila. These authors demonstrated that this rate is consistent with the rate of phenotypic assortment of a single gene, and suggested that the lower the copy number in the macronucleus, the higher the probability of expressing a deleterious allele in the homozygous form, and thus kill the cell. Our proposed mechanism also assumes that an increase in copy number in the macronuclear genome would provide higher robustness to clonal extinction. In fact, the selection coefficient (s) against a change in copy number away from an optimum value (Zm) may depend upon the absolute number of copies. It is easy to imagine that a cell carrying 44 copies of a gene whose optimum number of copies is 45 will have a lower fitness relative to a cell carrying 89 copies whose optimum value is 90. If this is the case then we can easily understand the difference in results of the MA of T. thermophila and T. pyriformis. For the same rate U a smaller s will lead to a much lower level of extinctions in the same time period (Figure 6). Lynch and Gabriel (1990) also speculate that the high longevities attributed to the amicronucleated Tetrahymena could be associated with an increase in mean extinction time when the effects of mutations are highly variable relative to their expectations - exactly what is expected when the gain and loss of chromosomal copies is achieved through a mechanism that is random relatively to the DNA composition of what is transferred.
The ecology associated with natural populations of Tetrahymena imply that they periodically suffer drastic fluctuations in population size  and hence suggests that T. thermophila species must recurrently undergo sexual reproduction in order to avoid the deleterious effects seen in our laboratory experiments due to extreme and periodic size reductions. The robustness of T. pyriformis to the accumulation of deleterious mutations in the MA lines, on the other hand, requires that species that lack the resetting mechanism of sex must take advantage of another mechanism that would confer protection from the deleterious effects of randomly accumulated mutations. The dependence of the rate of phenotypic assortment on the level of copy number suggests that longevity would benefit from (i) high copy number that allows for longer assortment times [27, 39], and (ii) shorter chromosomes sizes so that the linkage between deleterious and beneficial mutations can be broken down thus promoting more efficient selection. Our computer simulations show that the observed rate of clonal death can be obtained simply as a consequence of the evolutionary dynamics of an unfaithful mechanism of copy number segregation that, in the face of strong genetic drift, is unable to maintain the cell line at an optimum number of copies.
One important result from our study is the suggestion that the mechanism that regulates copy number in the macronucleus is natural selection acting on the diversity generated by a stochastic process of gain and loss of chromosomes. It has always been a puzzle how the continuous amitotic mechanism of macronuclear division does not lead to massive chromosomal imbalances and death from aneuploidy For this reason, a regulatory mechanism of copy number control through replication was favored [45, 46]. It was observed that a cell-cycle mechanism such as the one proposed by Preer and Preer  can rescue the cells from extreme chromosomal imbalances. But the fact that once started replication would duplicate all DNA content present in the macronucleus  strongly suggests this mechanism would be a last resort solution.
This is the not the first time that a stochastic model is advocated to explain the chromosomal segregation in ciliate macronucleus . However, unlike previous models we postulate that changes in chromosome copy number, even the smaller changes compatible to what is expected from each cell fission , carry fitness costs that are seen by selection. We can therefore hypothesize that under large effective population sizes the probability of population extinction through genetic imbalances should be small. In fact, once we allow for selection to be effective (by increasing Ne) no such rate of extinction is observed, suggesting that the rate of extinction was due to ineffective selection in the face of intense genetic drift. Our model does not explain assortment of rDNA given these genes are present in such high copy numbers unless these copies occur in synteny. Interestingly, Merrian and Burns  estimated a rate of assortment for rdnA2 [a mutation in the gene for ribosomal DNA that is present in 10000 copies, ] that is similar to the assortment rate of other genes maintained at 45 copies. It was suggested an epigenetic regulation independent of the assortment of other genes could explain the regulation of these highly amplified gene copies [45, 49].
The aim of this study was to estimate the somatic mutation rate of T. thermophila by performing a classic experiment of mutation accumulation. We observed a rate of clonal death and a fitness dynamics that are consistent with a mechanism of gain/loss of chromosomes, as is characteristic of the amitotic division of Tetrahymena macronucleus. The high similarity of estimates available for the mutation load in T. thermophila also strengthens the hypothesis that the mechanism responsible for the accumulation of deleterious effects is dependent on a species-specific feature such as chromosomal copy number. We suggest that copy number and the rate of chromosomal gain/loss at each amitotic division are crucial to understand macronuclear degeneration, cell line longevity, and senescence in ciliates. It also renders the hypothesis of longevity controlled by programmed cell death unlikely.
Tetrahymena thermophila, B2086, mating type II, and Tetrahymena pyriformis GL-C; these strains were obtained from the Tetrahymena Stock Center at Cornell University http://tetrahymena.vet.cornell.edu.
Mutation accumulation experiment
At time T0 twenty lines were established from single cells all originated from the same founding clonal culture. Because we were interested in estimating the somatic mutation rate, our initial cell was randomly selected from a culture of mature cells that were grown at large population sizes in standard conditions in the lab. This allowed us to start the experiment with clonal populations that had already gone through phenotypic assortment and experienced the effect of natural selection. The twenty cells were incubated and let to grow in individual drops of 30 μl of rich medium (SPP: 1% proteose peptone, 0.1% yeast extract, 0.22% glucose, and 0.003% EDTA ferric sodium salt) for 48 hours that corresponds to approximately 16 generations in the ancestral cell population. Drops were maintained on sterile Petri dishes placed in humid boxes at 30° and 28° for T. thermophila and T. pyriformis, respectively, which are the temperatures normally used in laboratory to grow cultures of each species. These lines were periodically bottlenecked (every 48 hours) by randomly sampling one individual cell and transferred it to a fresh medium drop. This procedure was repeated for 20 bottlenecks or until extinction of all lines. Every two bottlenecks populations were frozen and fitness was measured through competition with the ancestral type. Whenever a drop was detected with no more remaining live cells, a new cell was taken from the previous drop in an attempted to maintain the culture line, in most cases these attempts were not successful and a the culture line would be considered extinct.
To evaluate the effect of bottleneck size in the extinction rate obtained for T. thermophila an additional experiment was performed where twenty lines of T. thermophila all derived from the same founding clonal population, were grown in the same experimental conditions as above but with daily transferences of 1 μl that corresponded to approximately 1000 individual cells.
Mean and variance in population fitness were measured in the ancestral population and in the MA population replicates every two bottlenecks until extinction of all lines. Fitness was measured through competition of the experimental samples of T. thermophila with a T. thermophila CU522 strain expressing GFP under the promoter MTT1 (metallothionein 1 protein) dependent of cadmium chloride (Seixas et al unpublished data) and relative fitness of the evolved strains relative to the ancestral ones was interpolated from the competition assays done with the GFP-expressing cultures. Competitions were performed by putting one cell of each strain in a single medium drop and letting it incubate at 30°C for 48 hours. Initial cell number of each genotype (Ni) was 1, and final densities (Nf) were estimated from three sets of 10 μl fluorescent induced cultures collected from each competition drop and mounted in slides that were observed with a Leica DMRA2 microscope. GFP expression was induced by 2 μl of CdCl2 (0, 5 μg/ml) to each drop for one hour, and slides were prepared by adding 10 μl of extraction buffer and 14 μl of 2% paraformaldehyde to the 10 μl culture sample, that were let to dry for 30 min at 30°C, and later mounted in a slide with DABCO. All competitions were performed in triplicate. Live images from Tetrahymena cells were obtained with Leica DMRA2 microscope, and nuclei were stained with Hoechst 33342 for 15 minutes at 30°, at a final concentration of 2 μg/μl.
DNA extractions and Q-PCR assays
Sequence Genbank accession numbers and primers used to assay gene copy number in T. thermophila and T. pyriformis
GenBank Acession number
Q-PCR assay was done twice for each gene and in each experiment each sample was run in quadruplicate. The threshold cycle (CT) for each well was calculated following the instrument instructions. Standard curves were prepared for DNA templates of both species and the relative differences in copy number (level of polyploidization) was estimated for each gene from the ratio of the averaged CT obtained for T. pyriformis and T. thermophila at equal DNA concentrations (1X).
We built a model that simulates the experimental conditions of the mutation accumulation experiment as well as the variation in the chromosome copy number that occurs at each cell division due to the amitotical division of the Mac genome. Cells in our model are characterized by their mean chromosomal copy number (Z), and their fitness (viability) is solely dependent on how Z deviates from an optimal value (Zm). A population starting from a single cell would have Z=Zm = 45 copies in the macronucleus genome, and for simplicity this initial cell does not carry any recessive deleterious mutations. At every generation the number of copies either increases or decreases, with equal probability, at a rate U according to a Poisson distribution. Note that U parameterizes variation in copy number and not rate of acquisition of point mutations, as typically assumed in other organisms where mutation accumulation experiments have been performed. The relative fitness of an individual is determined by W(Z) = exp(-s|Z-Zm|), where s is the selection strength against changes in the copy number, and Z-Zm is the distance between that individuals' chromosomal copy number and the species' optimal copy number. This fitness landscape is the simplest function reflecting stabilizing selection for optimal copy number, and s is the effect on fitness of a change in one copy from the optimum number. When s is small this is the same as the effect of a single deleterious point mutation in classical models of MA . Every generation population size fluctuates with reproduction (by duplication), mutation and viability selection. We let the population evolve until stationarity followed by a bottleneck of a randomly selected individual, or until extinction. The simulation ends after 16 bottlenecks or until all replicate lines are extinct. We explored our model to determine the parameter space of U (mutation rate towards increase or decrease in chromosomal copy number) and s (selection strength against changes in copy number) that recovers results compatible with the empirical observations, namely the number of extinctions observed in the 20 populations over the 16 bottlenecks, the mean and variance in extinction time observed experimentally. Specifically we initially searched a large parameter space defined by U in the interval [0.1; 12.0] and s in the interval [0.01; 1.0], in spaces of 0.05 for each of the parameters. This allowed us to eliminate combinations of U and s that lead to simulation results where the pattern of extinction in the experiment was not produced. From this search we restricted our parameter space to new intervals for U and s that were now searched in smaller spaces, for compatibility with the extinction pattern of the data. This involved and extensive search for U between 0.1 and 10.0, every 0.1, and s between 0.02 and 0.20, every 0.005. For the same parameter values of U and s, the simulation run was replicated 20 times. With this approach we can delimit the range of values of U and s that in our model produce results for number of extinctions, mean and variance of extinction time, that are compatible with the empirical results. The dynamics of fitness decline obtained for each set (U, s) compatible with the empirical results for the rate of extinction can then be determined in the simulations and compared with the fitness dynamics measured experimentally.
This research was supported by Fundação para a Ciência e Tecnologia: PTDC/BIA-BDE/73163/2006 to I.G. and SFRH/BPD/26852/2006 to PHB. EG was supported by IEFP. We also thank Miguel Godinho Ferreira for helpful discussions, Michael Parkhouse for reviewing the manuscript, and Inês Conceição for technical assistance.
- Otto SP, Lenormand T: Resolving the paradox of sex and recombination. Nat Rev Gen. 2002, 3 (4): 252-261. 10.1038/nrg761.View ArticleGoogle Scholar
- Gordo I, Charlesworth B: Genetic linkage and molecular evolution. Curr Biol. 2001, 11 (17): 684-686. 10.1016/S0960-9822(01)00408-0.View ArticleGoogle Scholar
- Kondrashov AS, Crow JF: Haploidy or diploidy: which is better?. Nature. 1991, 351 (6324): 314-315. 10.1038/351314a0.View ArticlePubMedGoogle Scholar
- Gerstein AC, Otto SP: Ploidy and the causes of genomic evolution. J Hered. 2009, 100 (5): 571-581. 10.1093/jhered/esp057.View ArticlePubMedGoogle Scholar
- Lynch M, Conery J, Burger R: Mutation accumulation and the extinction of small populations. American Naturalist. 1995, 146 (4): 489-518. 10.1086/285812.View ArticleGoogle Scholar
- Halligan DL, Keightley PD: Spontaneous mutation accumulation studies in evolutionary genetics. Ann Rev Ecol Syst. 2009, 40: 151-172. 10.1146/annurev.ecolsys.39.110707.173437.View ArticleGoogle Scholar
- Mukai T, Chigusa SI, Mettler LE, Crow JF: Mutation rate and dominance of genes affecting viability in Drosophila melanogaster. Genetics. 1972, 72 (2): 335-355.PubMed CentralPubMedGoogle Scholar
- Drake JW, Charlesworth B, Charlesworth D, Crow JF: Rates of spontaneous mutation. Genetics. 1998, 148 (4): 1667-1686.PubMed CentralPubMedGoogle Scholar
- Keightley PD: The distribution of mutation effects on viability in Drosophila melanogaster. Genetics. 1994, 138 (4): 1315-1322.PubMed CentralPubMedGoogle Scholar
- Lynch M, Bürger R, Butcher D, Gabriel W: The mutational meltdown in asexual populations. J Hered. 1993, 84 (5): 339-344.PubMedGoogle Scholar
- Funchain P, Yeung A, Stewart JL, Lin R, Slupska MM, Miller JH: The consequences of growth of a mutator strain of Escherichia coli as measured by loss of function among multiple gene targets and loss of fitness. Genetics. 2000, 154 (3): 959-970.PubMed CentralPubMedGoogle Scholar
- Zeyl C, Mizesko M, de Visser JA: Mutational meltdown in laboratory yeast populations. Evolution. 2001, 55 (5): 909-917. 10.1554/0014-3820(2001)055[0909:MMILYP]2.0.CO;2.View ArticlePubMedGoogle Scholar
- Clarke DK, Duarte EA, Moya A, Elena SF, Domingo E, Holland J: Genetic bottlenecks and population passages cause profound fitness differences in RNA viruses. J Virol. 1993, 67 (1): 222-228.PubMed CentralPubMedGoogle Scholar
- Kibota T, Lynch M: Estimate of the genomic mutation rate deleterious to overall fitness in E. coli. Nature. 1996, 381 (6584): 694-696. 10.1038/381694a0.View ArticlePubMedGoogle Scholar
- Trindade S, Perfeito L, Gordo I: Rate and effects of spontaneous mutations that effect fitness in mutator Escherichia coli. Phil Trans R Soc Lond B. 2010, 00: 1-10.Google Scholar
- Zeyl C, DeVisser JA: Estimates of the rate and distribution of fitness effects of spontaneous mutation in Saccharomyces cerevisiae. Genetics. 2001, 157 (1): 53-61.PubMed CentralPubMedGoogle Scholar
- Schultz ST, Lynch M, Willis JH: Spontaneous deleterious mutation in Arabidopsis thaliana. PNAS. 1999, 96 (20): 11393-11398. 10.1073/pnas.96.20.11393.PubMed CentralView ArticlePubMedGoogle Scholar
- Keightley PD, Caballero A: Genomic mutation rates for lifetime reproductive output and lifespan in Caenorhabditis elegans. PNAS. 1997, 94 (8): 3823-3827. 10.1073/pnas.94.8.3823.PubMed CentralView ArticlePubMedGoogle Scholar
- Lynch M, Gabriel W: Mutation load and the survival of small populations. Evolution. 1990, 44 (7): 1725-1737. 10.2307/2409502.View ArticleGoogle Scholar
- Hausmann K, Bradbury PC: Ciliates: cells as organisms. 1996, Stuttgart: Gustav FisherGoogle Scholar
- Katz LA: Evolution of nuclear dualism in ciliates: a reanalysis in light of recent molecular data. Int J Syst Evol Microbiol. 2001, 51 (Pt 4): 1587-1592.View ArticlePubMedGoogle Scholar
- Prescott DM: The DNA of ciliated protozoa. Microbiol Rev. 1994, 58 (2): 233-267.PubMed CentralPubMedGoogle Scholar
- Raikov IB: Nuclei of ciliates. Ciliates: cells as organisms. Edited by: Hausmann K, Bradbury PC. 1996, Stuttgart: Gustav Fisher, 221-242.Google Scholar
- Eisen JA, Coyne RS, Wu M, Wu D, Thiagarajan M, Wortman JR, Badger JH, Ren Q, Amedeo P, Jones KM, et al: Macronuclear genome sequence of the ciliate Tetrahymena thermophila, a model eukaryote. PLoS Biology. 2006, 4 (9): e286-10.1371/journal.pbio.0040286.PubMed CentralView ArticlePubMedGoogle Scholar
- Conover RK, Brunk CF: Macronuclear DNA molecules of Tetrahymena thermophila. Mol Cell Biol. 1986, 6 (3): 900-905.PubMed CentralView ArticlePubMedGoogle Scholar
- Yao MC, Choi J, Yokoyama S, Austerberry CF, Yao CH: DNA elimination in Tetrahymena: a developmental process involving extensive breakage and rejoining of DNA at defined sites. Cell. 1984, 36: 433-440. 10.1016/0092-8674(84)90236-8.View ArticlePubMedGoogle Scholar
- Merriam EV, Bruns PJ: Phenotypic assortment in Tetrahymena thermophila: assortment kinetics of antibiotic-resistance markers, tsA, death, and the highly amplified rDNA locus. Genetics. 1988, 120 (2): 389-395.PubMed CentralPubMedGoogle Scholar
- Doerder FP, Deak JC, Lief JH: Rate of phenotypic assortment in Tetrahymena thermophila. Dev Genet. 1992, 13 (2): 126-132. 10.1002/dvg.1020130206.View ArticlePubMedGoogle Scholar
- Berger JD, Schmidt HJ: Regulation of macronuclear DNA content in Paramecium tetraurelia. J Cell Biol. 1978, 76 (1): 116-126. 10.1083/jcb.76.1.116.View ArticlePubMedGoogle Scholar
- Doerder FP, DeBault LE: Life cycle variation and regulation of macronuclear DNA content in Tetrahymena thermophila. Chromosoma. 1978, 69 (1): 1-19. 10.1007/BF00327377.View ArticlePubMedGoogle Scholar
- Doerder FP, Ditaranto J, DeBault LE: Evolutionary constraints on quantitative variation and regulation of macronuclear dna content in the genus Tetrahymena. J Cell Sci. 1981, 49: 177-193.PubMedGoogle Scholar
- Takagi Y, Yoshida M: Clonal death associated with the number of fissions in Paramecium caudatum. J Cell Sci. 1980, 41: 177-191.PubMedGoogle Scholar
- Cleffmann G: Amount of DNA produced during extra S phases in Tetrahymena. J Cell Biol. 1975, 66 (1): 204-209. 10.1083/jcb.66.1.204.View ArticlePubMedGoogle Scholar
- Doerder FP: Regulation of Macronuclear DNA Content in Tetrahymena thermophila. Journal of Protozoology. 1979, 26 (1): 28-35.View ArticleGoogle Scholar
- Nanney DL, Preparata RM: Genetic evidence concerning the structure of the Tetrahymena thermophila macronucleus. Journal of Protozoology. 1979, 26: 2-9.View ArticleGoogle Scholar
- Orias E, Flacks M: Macronuclear genetics of Tetrahymena. I. Random distribution of macronuclear gene copies in T. pyriformis, SYNGEN 1. Genetics. 1975, 79: 187-206.PubMed CentralPubMedGoogle Scholar
- Smith-Sonneborn J: DNA repair and longevity assurance in Paramecium tetraurelia. Science. 1979, 203 (4385): 1115-1117. 10.1126/science.424739.View ArticlePubMedGoogle Scholar
- Sonneborn T: The relation of autogamy to senescence and rejuvenescence in Paramecium aurelia. Journal of Eukaryotic Microbiology. 1954, 1 (1): 38-53. 10.1111/j.1550-7408.1954.tb00792.x.Google Scholar
- Bell G: Sex and death in protozoa. The history of an obsession. 1988, Cambridge: Cambridge University PressGoogle Scholar
- Nanney DL: Aging and long-term temporal regulation in ciliated protozoa. A critical review. Mech Ageing Dev. 1974, 3 (2): 81-105. 10.1016/0047-6374(74)90008-6.View ArticlePubMedGoogle Scholar
- Estes S, Phillips P, Denver D, Thomas W, Lynch M: Mutation accumulation in populations of varying size: the distribution of mutational effects for fitness correlates in Caenorhabditis elegans. Genetics. 2004, 166 (3): 1269-10.1534/genetics.166.3.1269.PubMed CentralView ArticlePubMedGoogle Scholar
- Bodenbender J, Prohaska A, Jauker F, Hipke H, Cleffmann G: DNA elimination and its relation to quantities in the macronucleus of Tetrahymena. Developmental Genetics. 1992, 13 (2): 103-110. 10.1002/dvg.1020130203.View ArticlePubMedGoogle Scholar
- Lynch M, Conery J, Burger R: Mutational meltdowns in sexual populations. Evolution. 1995, 49 (6): 1067-1080. 10.2307/2410432.View ArticleGoogle Scholar
- Doerder FP, Arslanyolu M, Saad Y, Kaczmarek M, Mendoza M, Mita B: Ecological genetics of Tetrahymena thermophila: mating types, i-antigens, multiple alleles and epistasis. J Euk Microbiol. 1996, 43 (2): 95-100. 10.1111/j.1550-7408.1996.tb04487.x.View ArticleGoogle Scholar
- Larson DD, Umthun AR, Shaiu WL: Copy number control in the Tetrahymena macronuclear genome. J Protozool. 1991, 38 (3): 258-263.View ArticlePubMedGoogle Scholar
- Preer JR, Preer LB: The size of macronuclear DNA and Its relationship to models for maintaining genic balance. Journal of Protozoology. 1979, 26: 14-18.View ArticleGoogle Scholar
- Duerr HP, Eichner M, Ammermann D: Modeling senescence in hypotrichous ciliates. Protist. 2004, 155 (1): 45-52. 10.1078/1434461000163.View ArticlePubMedGoogle Scholar
- Gorovsky MA: Genome organization and reorganization in Tetrahymena. Annu Rev Genet. 1980, 14: 203-239. 10.1146/annurev.ge.14.120180.001223.View ArticlePubMedGoogle Scholar
- Juranek SA, Lipps HJ: New insights into the macronuclear development in ciliates. Int Rev Cytol. 2007, 262: 219-251. 10.1016/S0074-7696(07)62005-1.View ArticlePubMedGoogle Scholar
- Gaertig J, Thatcher T, Gu L, Gorovsky M: Electroporation-mediated replacement of a positively and negatively selectable-tubulin gene in Tetrahymena thermophila. PNAS. 1994, 91: 4549-4553. 10.1073/pnas.91.10.4549.PubMed CentralView ArticlePubMedGoogle Scholar
- Kuo C-H, Wares JP, Kissinger C: The apicomplexan whole-genome phylogeny: an analysis of incongruence among gene trees. Molecular Biology and Evolution. 2008, 25 (12): 2689-2698. 10.1093/molbev/msn213.PubMed CentralView ArticlePubMedGoogle Scholar
- McGrath KE, Yu SM, Heruth DP, Kelly AA, Gorovsky MA: Regulation and evolution of the single alpha-tubulin gene of the ciliate Tetrahymena thermophila. Cell Motil Cytoskeleton. 1994, 27 (3): 272-283. 10.1002/cm.970270308.View ArticlePubMedGoogle Scholar
- Soares H, Cyrne L, Barahona I, Rodrigues-Pousada C: Different patterns of expression of beta-tubulin genes in Tetrahymena pyriformis during reciliation. Eur J Biochem. 1991, 197 (2): 291-299. 10.1111/j.1432-1033.1991.tb15910.x.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.