Accumulation of Spontaneous Mutations in the Ciliate Tetrahymena thermophila

Knowledge of the rate and fitness effects of mutations is essential for understanding the process of evolution. Mutations are inherently difficult to study because they are rare and are frequently eliminated by natural selection. In the ciliate Tetrahymena thermophila, mutations can accumulate in the germline genome without being exposed to selection. We have conducted a mutation accumulation (MA) experiment in this species. Assuming that all mutations are deleterious and have the same effect, we estimate that the deleterious mutation rate per haploid germline genome per generation is U = 0.0047 (95% credible interval: 0.0015, 0.0125), and that germline mutations decrease fitness by s = 11% when expressed in a homozygous state (95% CI: 4.4%, 27%). We also estimate that deleterious mutations are partially recessive on average (h = 0.26; 95% CI: –0.022, 0.62) and that the rate of lethal mutations is <10% of the deleterious mutation rate. Comparisons between the observed evolutionary responses in the germline and somatic genomes and the results from individual-based simulations of MA suggest that the two genomes have similar mutational parameters. These are the first estimates of the deleterious mutation rate and fitness effects from the eukaryotic supergroup Chromalveolata and are within the range of those of other eukaryotes.     

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.     

Intergenomic epistasis for fitness: within-population interactions between cytoplasmic and nuclear genes in Drosophila melanogaster

The symbiotic relationship between the mitochondrial and nuclear genomes coordinates metabolic energy production and is fundamental to life among eukaryotes. Consequently, there is potential for strong selection to shape interactions between these two genomes. Substantial research attention has focused on the possibility that within-population sequence polymorphism in mitochondrial DNA (mtDNA) is maintained by mitonuclear fitness interactions. Early theory predicted that selection will often eliminate mitochondrial polymorphisms. However, recent models demonstrate that intergenomic interactions can promote the maintenance of polymorphism, especially if the nuclear genes involved are linked to the X chromosome. Most empirical studies to date that have assessed cytonuclear fitness interactions have studied variation across populations and it is still unclear how general and strong such interactions are within populations. We experimentally tested for cytonuclear interactions within a laboratory population of Drosophila melanogaster using 25 randomly sampled cytoplasmic genomes, expressed in three different haploid nuclear genetic backgrounds, while eliminating confounding effects of intracellular bacteria (e.g., Wolbachia). We found sizable cytonuclear fitness interactions within this population and present limited evidence suggesting that these effects were sex specific. Moreover, the relative fitness of cytonuclear genotypes was environment specific. Sequencing of mtDNA (2752 bp) revealed polymorphism within the population, suggesting that the observed cytoplasmic genetic effects may be mitochondrial in origin.     

Massive colonization of protein-coding exons by selfish genetic elements in Paramecium germline genomes

Ciliates are unicellular eukaryotes with both a germline genome and a somatic genome in the same cytoplasm. The somatic macronucleus (MAC), responsible for gene expression, is not sexually transmitted but develops from a copy of the germline micronucleus (MIC) at each sexual generation. In the MIC genome of Paramecium tetraurelia, genes are interrupted by tens of thousands of unique intervening sequences called internal eliminated sequences (IESs), which have to be precisely excised during the development of the new MAC to restore functional genes. To understand the evolutionary origin of this peculiar genomic architecture, we sequenced the MIC genomes of 9 Paramecium species (from approximately 100 Mb in Paramecium aurelia species to >1.5 Gb in Paramecium caudatum). We detected several waves of IES gains, both in ancestral and in more recent lineages. While the vast majority of IESs are single copy in present-day genomes, we identified several families of mobile IESs, including nonautonomous elements acquired via horizontal transfer, which generated tens to thousands of new copies. These observations provide the first direct evidence that transposable elements can account for the massive proliferation of IESs in Paramecium. The comparison of IESs of different evolutionary ages indicates that, over time, IESs shorten and diverge rapidly in sequence while they acquire features that allow them to be more efficiently excised. We nevertheless identified rare cases of IESs that are under strong purifying selection across the aurelia clade. The cases examined contain or overlap cellular genes that are inactivated by excision during development, suggesting conserved regulatory mechanisms. Similar to the evolution of introns in eukaryotes, the evolution of Paramecium IESs highlights the major role played by selfish genetic elements in shaping the complexity of genome architecture and gene expression.

A comparative genomics study of nine Paramecium species reveals successful invasion of genes by transposable elements in their germline genomes, showing that the internal eliminated sequences (IESs) followed an evolutionary trajectory remarkably similar to that of spliceosomal introns.     

Genome diversity in microbial eukaryotes

The genomic peculiarities among microbial eukaryotes challenge the conventional wisdom of genome evolution. Currently, many studies and textbooks explore principles of genome evolution from a limited number of eukaryotic lineages, focusing often on only a few representative species of plants, animals and fungi. Increasing emphasis on studies of genomes in microbial eukaryotes has and will continue to uncover features that are either not present in the representative species (e.g. hypervariable karyotypes or highly fragmented mitochondrial genomes) or are exaggerated in microbial groups (e.g. chromosomal processing between germline and somatic nuclei). Data for microbial eukaryotes have emerged from recent genome sequencing projects, enabling comparisons of the genomes from diverse lineages across the eukaryotic phylogenetic tree. Some of these features, including amplified rDNAs, subtelomeric rDNAs and reduced genomes, appear to have evolved multiple times within eukaryotes, whereas other features, such as absolute strand polarity, are found only within single lineages.     

Molecular Population Genetics of Rice Domestication

Domestication is a selection process that genetically modifies species to meet human needs. A most intriguing feature of domestication is the extreme phenotypic diversification among breeds. What could be the ultimate source of such genetic variations? Another notable outcome of artificial selection is the reduction in the fitness of domesticated species when they live in the wild without human assistance. The complete sequences of the two subspecies of rice cultivars provide an opportunity to address these questions. Between the two subspecies, we found much higher rates of non‐synonymous (N) than synonymous (S) substitutions and the N/S ratios are higher between cultivars than between wild species. Most interestingly, substitutions of highly dissimilar amino acids that are deleterious and uncommon between natural species are disproportionately common between the two subspecies of rice. We suggest strong selection in the absence of effective recombination may be the driving force, which we called the domestication‐associated Hill‐Robertson effect. These hitchhiking mutations may contribute to some fitness reduction in cultivars. Comparisons of the two genomes also reveal the existence of highly divergent regions in the genomes. Haplotypes in these regions often form highly polymorphic linkage blocks that are much older than speciation between wild species. Genes from such regions could contribute to the differences between indica and japonica and are likely to be involved in the diversifying selection under domestication. Their existence suggests that the amount of genetic variation within the single progenitor species Oryza rufipogon may be insufficient to account for the variation among rice cultivars, which may come from a more inclusive gene pool comprising most of the A‐genome wild species. Genes from the highly polymorphic regions also provide strong support for the independent domestication of the two subspecies. The genomic variation in rice has revealing implications for studying the genetic basis of indica‐japonica differentiation under rice domestication and subsequent improvement.     

Temperature-specific outcomes of cytoplasmic-nuclear interactions on egg-to-adult development time in seed beetles

The integration of the mitochondrial and nuclear genomes coordinates cellular energy production and is fundamental to life among eukaryotes. Therefore, there is potential for strong selection to shape the interactions between the two genomes. Several studies have now demonstrated that epistatic interactions between cytoplasmic and nuclear genes for fitness can occur both at a "within" and "across" population level. Genotype-by-environment interactions are common for traits that are encoded by nuclear genes, but the effects of environmental heterogeneity on traits that are partly encoded by cytoplasmic genes have received little attention despite the fact that there are reasons to believe that phenotypic effects of cytoplasmic genetic variation may often be environment specific. Consequently, the importance of environmental heterogeneity to the outcomes of cyto-nuclear fitness interactions and to the maintenance of mitochondrial polymorphism is unclear. Here, we assess the influence of temperature on cyto-nuclear effects on egg-to-adult development time in seed beetles (Callosobruchus maculatus). We employed an "across-population" design, sourcing beetles from five distinct populations and using backcrossing to create orthogonal combinations of distinct introgression lines, fixed for their cytoplasmic and nuclear lineages. We then assayed development times at two different temperatures and found sizeable cyto-nuclear effects in general, as well as temperature- and block-specific cyto-nuclear effects. These results demonstrate that environmental factors such as temperature do exert selection on cytoplasmic genes by favoring specific cyto-nuclear genetic combinations, and are consistent with the suggestion that complex genotype-by-environment interactions may promote the maintenance of polymorphism in mitochondrial genes.     

Applying the genetic theories of ageing to the cytoplasm: cytoplasmic genetic covariation for fitness and lifespan

Two genetic models exist to explain the evolution of ageing – mutation accumulation (MA) and antagonistic pleiotropy (AP). Under MA, a reduced intensity of selection with age results in accumulation of late‐acting deleterious mutations. Under AP, late‐acting deleterious mutations accumulate because they confer beneficial effects early in life. Recent studies suggest that the mitochondrial genome is a major player in ageing. It therefore seems plausible that the MA and AP models will be relevant to genomes within the cytoplasm. This possibility has not been considered previously. We explore whether patterns of covariation between fitness and ageing across 25 cytoplasmic lines, sampled from a population of Drosophila melanogaster, are consistent with the genetic associations predicted under MA or AP. We find negative covariation for fitness and the rate of ageing, and positive covariation for fitness and lifespan. Notably, the direction of these associations is opposite to that typically predicted under AP.     

MITOCHONDRIAL–NUCLEAR EPISTASIS AFFECTS FITNESS WITHIN SPECIES BUT DOES NOT CONTRIBUTE TO FIXED INCOMPATIBILITIES BETWEEN SPECIES OF DROSOPHILA

Efficient mitochondrial function requires physical interactions between the proteins encoded by the mitochondrial and nuclear genomes. Coevolution between these genomes may result in the accumulation of incompatibilities between divergent lineages. We test whether mitochondrial–nuclear incompatibilities have accumulated within the Drosophila melanogaster species subgroup by combining divergent mitochondrial and nuclear lineages and quantifying the effects on relative fitness. Precise placement of nine mtDNAs from D. melanogaster, D. simulans, and D. mauritiana into two D. melanogaster nuclear genetic backgrounds reveals significant mitochondrial–nuclear epistasis affecting fitness in females. Combining the mitochondrial genomes with three different D. melanogaster X chromosomes reveals significant epistasis for male fitness between X‐linked and mitochondrial variation. However, we find no evidence that the more than 500 fixed differences between the mitochondrial genomes of D. melanogaster and the D. simulans species complex are incompatible with the D. melanogaster nuclear genome. Rather, the interactions of largest effect occur between mitochondrial and nuclear polymorphisms that segregate within species of the D. melanogaster species subgroup. We propose that a low mitochondrial substitution rate, resulting from a low mutation rate and/or efficient purifying selection, precludes the accumulation of mitochondrial–nuclear incompatibilities among these Drosophila species.     

Evolution of the cancer genome

Human tumors result from an evolutionary process operating on somatic cells within tissues, whereby natural selection operates on the phenotypic variability generated by the accumulation of genetic, genomic and epigenetic alterations. This somatic evolution leads to adaptations such as increased proliferative, angiogenic, and invasive phenotypes. In this review we outline how cancer genomes are beginning to be investigated from an evolutionary perspective. We describe recent progress in the cataloging of somatic genetic and genomic alterations, and investigate the contributions of germline as well as epigenetic factors to cancer genome evolution. Finally, we outline the challenges facing researchers who investigate the processes driving the evolution of the cancer genome.     

Fitness change in relation to mutation number in spontaneous mutation accumulation lines of Chlamydomonas reinhardtii

Although all genetic variation ultimately stems from mutations, their properties are difficult to study directly. Here, we used multiple mutation accumulation (MA) lines derived from five genetic backgrounds of the green algae Chlamydomonas reinhardtii that have been previously subjected to whole genome sequencing to investigate the relationship between the number of spontaneous mutations and change in fitness from a nonevolved ancestor. MA lines were on average less fit than their ancestors and we detected a significantly negative correlation between the change in fitness and the total number of accumulated mutations in the genome. Likewise, the number of mutations located within coding regions significantly and negatively impacted MA line fitness. We used the fitness data to parameterize a maximum likelihood model to estimate discrete categories of mutational effects, and found that models containing one to two mutational effect categories (one neutral and one deleterious category) fitted the data best. However, the best‐fitting mutational effects models were highly dependent on the genetic background of the ancestral strain.     

Exploration of Genetic Variations through Single‐cell Whole‐genome Sequencing in the Model Ciliate Tetrahymena thermophila

Ciliates are unicellular eukaryotes with separate germline and somatic genomes and diverse life cycles, which make them a unique model to improve our understanding of population genetics through the detection of genetic variations. However, traditional sequencing methods cannot be directly applied to ciliates because the majority are uncultivated. Single‐cell whole‐genome sequencing (WGS) is a powerful tool for studying genetic variation in microbes, but no studies have been performed in ciliates. We compared the use of single‐cell WGS and bulk DNA WGS to detect genetic variation, specifically single nucleotide polymorphisms (SNPs), in the model ciliate Tetrahymena thermophila. Our analyses showed that (i) single‐cell WGS has excellent performance regarding mapping rate and genome coverage but lower sequencing uniformity compared with bulk DNA WGS due to amplification bias (which was reproducible); (ii) false‐positive SNP sites detected by single‐cell WGS tend to occur in genomic regions with particularly high sequencing depth and high rate of C:G to T:A base changes; (iii) SNPs detected in three or more cells should be reliable (an detection efficiency of 83.4–97.4% was obtained for combined data from three cells). This analytical method could be adapted to measure genetic variation in other ciliates and broaden research into ciliate population genetics.     

Deleterious Mutations, Apparent Stabilizing Selection and the Maintenance of Quantitative Variation

Apparent stabilizing selection on a quantitative trait that is not causally connected to fitness can result from the pleiotropic effects of unconditionally deleterious mutations, because as N. Barton noted, "...individuals with extreme values of the trait will tend to carry more deleterious alleles...." We use a simple model to investigate the dependence of this apparent selection on the genomic deleterious mutation rate, U; the equilibrium distribution of K, the number of deleterious mutations per genome; and the parameters describing directional selection against deleterious mutations. Unlike previous analyses, we allow for epistatic selection against deleterious alleles. For various selection functions and realistic parameter values, the distribution of K, the distribution of breeding values for a pleiotropically affected trait, and the apparent stabilizing selection function are all nearly Gaussian. The additive genetic variance for the quantitative trait is kQa2, where k is the average number of deleterious mutations per genome, Q is the proportion of deleterious mutations that affect the trait, and a2 is the variance of pleiotropic effects for individual mutations that do affect the trait. In contrast, when the trait is measured in units of its additive standard deviation, the apparent fitness function is essentially independent of Q and a2; and beta, the intensity of selection, measured as the ratio of additive genetic variance to the "variance" of the fitness curve, is very close to s = U/k, the selection coefficient against individual deleterious mutations at equilibrium. Therefore, this model predicts appreciable apparent stabilizing selection if s exceeds about 0.03, which is consistent with various data. However, the model also predicts that beta must equal Vm/VG, the ratio of new additive variance for the trait introduced each generation by mutation to the standing additive variance. Most, although not all, estimates of this ratio imply apparent stabilizing selection weaker than generally observed. A qualitative argument suggests that even when direct selection is responsible for most of the selection observed on a character, it may be essentially irrelevant to the maintenance of variation for the character by mutation-selection balance. Simple experiments can indicate the fraction of observed stabilizing selection attributable to the pleiotropic effects of deleterious mutations.     

Experimental estimate of the abundance and effects of nearly neutral mutations in the RNA virus phi 6

Although the frequency and effects of neutral and nearly neutral mutations are critical to evolutionary patterns and processes governed by genetic drift, the small effects of such mutations make them difficult to study empirically. Here we present the results of a mutation-accumulation experiment designed to assess the frequencies of deleterious mutations with undetectable effects. We promoted the accumulation of spontaneous mutations by subjecting independent lineages of the RNA virus 6 to repeated population bottlenecks of a single individual. We measured fitness following every bottleneck to obtain a complete picture of the timing and effects of the accumulated mutations with detectable effects and sequenced complete genomes to determine the number of mutations that were undetected by the fitness assays. To estimate the effects of the undetected mutations, we implemented a likelihood model developed for quantitative trait locus (QTL) data (Otto and Jones 2000) to estimate the number and effects of the undetected mutations from the measured number and effects of the detected mutations. Using this method we estimated a deleterious mutation rate of U = 0.03 and a gamma effects distribution with mean s=0.093 and coefficient of variation = 0.204. Although our estimates of U and s fall within the range of recent mutation rate and effect estimates in eukaryotes, the fraction of mutations with detectable effects on laboratory fitness (39%) appears to be far higher in 6 than in eukaryotes.     

What cost mitochondria? The maintenance of functional mitochondrial DNA within and across generations

The peculiar biology of mitochondrial DNA (mtDNA) potentially has detrimental consequences for organismal health and lifespan. Typically, eukaryotic cells contain multiple mitochondria, each with multiple mtDNA genomes. The high copy number of mtDNA implies that selection on mtDNA functionality is relaxed. Furthermore, because mtDNA replication is not strictly regulated, within-cell selection may favour mtDNA variants with a replication advantage, but a deleterious effect on cell fitness. The opportunities for selfish mtDNA mutations to spread are restricted by various organism-level adaptations, such as uniparental transmission, germline mtDNA bottlenecks, germline selection and, during somatic growth, regular alternation between fusion and fission of mitochondria. These mechanisms are all hypothesized to maintain functional mtDNA. However, the strength of selection for maintenance of functional mtDNA progressively declines with age, resulting in age-related diseases. Furthermore, organismal adaptations that most probably evolved to restrict the opportunities for selfish mtDNA create secondary problems. Owing to predominantly maternal mtDNA transmission, recombination among mtDNA from different individuals is highly restricted or absent, reducing the scope for repair. Moreover, maternal inheritance precludes selection against mtDNA variants with male-specific effects. We finish by discussing the consequences of life-history differences among taxa with respect to mtDNA evolution and make a case for the use of microorganisms to experimentally manipulate levels of selection.     

PERSPECTIVE: SPONTANEOUS DELETERIOUS MUTATION

Mildly deleterious mutation has been invoked as a leading explanation for a diverse array of observations in evolutionary genetics and molecular evolution and is thought to be a significant risk of extinction for small populations. However, much of the empirical evidence for the deleterious-mutation process derives from studies of Drosophila melanogaster, some of which have been called into question. We review a broad array of data that collectively support the hypothesis that deleterious mutations arise in flies at rate of about one per individual per generation, with the average mutation decreasing fitness by about only 2% in the heterozygous state. Empirical evidence from microbes, plants, and several other animal species provide further support for the idea that most mutations have only mildly deleterious effects on fitness, and several other species appear to have genomic mutation rates that are of the order of magnitude observed in Drosophila. However, there is mounting evidence that some organisms have genomic deleterious mutation rates that are substantially lower than one per individual per generation. These lower rates may be at least partially reconciled with the Drosophila data by taking into consideration the number of germline cell divisions per generation. To fully resolve the existing controversy over the properties of spontaneous mutations, a number of issues need to be clarified. These include the form of the distribution of mutational effects and the extent to which this is modified by the environmental and genetic background and the contribution of basic biological features such as generation length and genome size to interspecific differences in the genomic mutation rate. Once such information is available, it should be possible to make a refined statement about the long-term impact of mutation on the genetic integrity of human populations subject to relaxed selection resulting from modern medical procedures.     

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.     

MALE‐BIASED FITNESS EFFECTS OF SPONTANEOUS MUTATIONS IN DROSOPHILA MELANOGASTER

In populations with males and females, sexual selection may often represent a major component of overall selection. Sexual selection could act to eliminate deleterious alleles in concert with other forms of selection, thereby improving the fitness of sexual populations. Alternatively, the divergent reproductive strategies of the sexes could promote the maintenance of sexually antagonistic variation, causing sexual populations to be less fit. The net impact of sexual selection on fitness is not well understood, due in part to limited data on the sex‐specific effects of spontaneous mutations on total fitness. Using a set of mutation accumulation lines of Drosophila melanogaster, we found that mutations were deleterious in both sexes and had larger effects on fitness in males than in females. This pattern is expected to reduce the mutation load of sexual females and promote the maintenance of sexual reproduction.     

Genomic parasites or symbionts? Modeling the effects of environmental pressure on transposition activity in asexual populations

Transposable elements are DNA segments capable of persisting in host genomes by self-replication in spite of deleterious mutagenic effects. The theoretical dynamics of these elements within genomes has been studied extensively, and population genetic models predict that they can invade and maintain as a result of both intra-genomic and inter-individual selection in sexual species. In asexuals, the success of selfish DNA is more difficult to explain. However, most theoretical work assumes constant environment. Here, we analyze the impact of environmental change on the dynamics of transposition activity when horizontal DNA exchange is absent, based on a stochastic computational model of transposable element proliferation. We argue that repeated changes in the phenotypic optimum in a multidimensional fitness landscape may induce explosive bursts of transposition activity associated with faster adaptation. However, long-term maintenance of transposition activity is unlikely. This could contribute to the significant variation in the transposable element copy number among closely related species.     

How do mammalian transposons induce genetic variation? A conceptual framework

In this essay, we discuss new insights into the wide‐ranging impacts of mammalian transposable elements (TE) on gene expression and function. Nearly half of each mammalian genome is comprised of these mobile, repetitive elements. While most TEs are ancient relics, certain classes can move from one chromosomal location to another even now. Indeed, striking recent data show that extensive transposition occurs not only in the germline over evolutionary time, but also in developing somatic tissues and particular human cancers. While occasional germline TE insertions may contribute to genetic variation, many other, similar TEs appear to have little or no impact on neighboring genes. However, the effects of somatic insertions on gene expression and function remain almost completely unknown. We present a conceptual framework to understand how the ages, allele frequencies, molecular structures, and especially the genomic context of mammalian TEs each can influence their various possible functional consequences. Editor's suggested further reading in BioEssays Evolution of eukaryotic genome architecture: Insights from the study of a rapidly evolving metazoan, Oikopleura dioica Abstract