The build up of mutation-selection- drift balance in laboratory Drosophila populations

The build up of an equilibrium between mutation, selection, and drift in populations of moderate size is an important evolutionary issue, and can be critical in the conservation of endangered populations. We studied this process in two Drosophila melanogaster populations initially lacking genetic variability (C1 and C2) that were subsequently maintained during 431 or 165 generations with effective population size N(e) approximately 500 (estimated by lethal complementation analysis). Each population originated synchronously to a companion set of full-sib mutation accumulation (MA) lines, C1 and MA1 were derived from an isogenic origin and C2 and MA2 from a single MA1 line at generation 265. The results suggest that both C1 and C2 populations were close to the mutation-selection-drift balance for viability and bristle traits, and are consistent with a 2.5-fold increase of the mutation rate in C2 and MA2. Despite this increase, the average panmictic viability in C2 was only slightly below that of C1, indicating that the expressed loads due to segregating deleterious mutation were small, in agreement with the low deleterious mutation rate (0.015-0.045) previously reported for the MA1 lines. In C1, the nonlethal inbreeding depression rate for viability was 30% of that usually estimated in segregating populations. The genetic variance for bristles regenerated in C1 and C2 was moderately smaller than the average value reported for natural populations, implying that they have accumulated a substantial adaptive potential. In light of neutral and selective predictions, these results suggest that bristle additive variance was predominantly due to segregation of mutations with deleterious effects of the order of 10(-3), and is consistent with relatively weak causal stabilizing selection (V(s) approximately 30).     

Accumulation of deleterious mutations: additional Drosophila melanogaster estimates and a simulation of the effects of selection

We report an assay of egg-to-adult viability in full-sibling mutation accumulation (MA) lines derived from a completely homozygous population of Drosophila melanogaster and maintained for 210 generations. A simultaneous evaluation was also made of a large population derived from the same origin and maintained as a control for the same period. We also present computer simulations to explore the possible decline in viability of the control population due to mutation accumulation and the possible effect of selection within and between MA lines. For this purpose, we used two mutational models independent from the data analyzed and based on radically different assumptions. The first model implies a large number of mutations of small effect, whereas the second implies a much smaller number of mutations with much larger effects. The observed rate of decline in mean viability was very small but significant (0.077%). The rate of increase in among line variance (0.189 x 10(-3)) was similar to those obtained previously in the same lines. The simulation results indicated that a model of many mutations of small effect is incompatible with the evolution of the mean viability of the control and MA lines over generations, the distribution of line means after 210 generations of mutation accumulation, and the pattern of line extinction over generations. Basically, this model predicted a large drop in viability, both in the control and particularly the MA lines, that is not observed empirically. It also predicted a rate of line extinction too low in the early generations and too high in the later ones. In contrast, the model based on few mutations of large effect was generally consistent with all the observations.     

The rate of mutation and the homozygous and heterozygous mutational effects for competitive viability: a long-term experiment with Drosophila melanogaster

The effect of 250 generations of mutation accumulation (MA) on the second chromosome competitive viability of Drosophila melanogaster was analyzed both in homozygous and heterozygous conditions. We used full-sib MA lines, where selection hampers the accumulation of severely deleterious mutations but is ineffective against mildly deleterious ones. A large control population was simultaneously evaluated. Competitive viability scores, unaffected by the expression of mutations in heterozygosis, were obtained relative to a Cy/L(2) genotype. The rate of decline in mean DeltaM approximately 0.1% was small. However, that of increase in variance DeltaV approximately 0.08 x 10(-3) was similar to the values obtained in previous experiments when severely deleterious mutations were excluded. The corresponding estimates of the mutation rate lambda > or = 0.01 and the average effect of mutations E(s) < or = 0.08 are in good agreement with Bateman-Mukai and minimum distance estimates for noncompetitive viability obtained from the same MA lines after 105 generations. Thus, competitive and noncompetitive viability show similar mutational properties. The regression estimate of the degree of dominance for mild-to-moderate deleterious mutations was approximately 0.3, suggesting that the pertinent value for new unselected mutations should be somewhat smaller.     

Inferences on the role of insertion in a mutation accumulation experiment with Drosophila melanogaster using RAPDs

The genetic variability for RAPDs band pattern was studied in a set of 157 mutation accumulation (MA) lines of Drosophila melanogaster. These MA lines were derived from the same isogenic base population and subsequently maintained by full-sib mating during 132 generations. The ancestral pattern of the original isogenic base can be unambiguously established as the consensus pattern of the MA lines and, because these lines are expected to be homozygous, dominance for band pattern is not a concern. Only repeatable changes in band pattern were considered. The number of ancestral bands detected implies that nine-nucleotide targets are enough for repeatable PCR amplification. Compared with the ancestral pattern, one MA line lost one band and two MA lines gained a new one. These results can be accounted for by the insertion of transposable elements occurring at a rate 0.07 < i < 0.21 per whole haploid genome and generation. This range is typical for Drosophila and consistent with the previously observed mobility for the roo family, supporting the generality of previous estimates of spontaneous mutation rates for morphological and fitness traits based on these MA lines. The sequence of one of the new bands suggests that the Idefix family is also active in the lines.     

The effects of spontaneous mutation on competitive fitness in Drosophila melanogaster

Abstract We have analysed the effect of 288 generations of mutation accumulation (MA) on chromosome II competitive fitness in 21 full‐sib lines of Drosophila melanogaster and in a large control population, all derived from the same isogenic base. The rate of mean log‐fitness decline and that of increase of the between‐line variance were consistent with a low rate (λ ≈ 0.03 per gamete and generation), and moderate average fitness effect [E(s) ≈ 0.1] of deleterious mutation. Subsequently, crosses were made between pairs of MA lines, and these were maintained with effective size on the order of a few tens. In these crosses, MA recombinant chromosomes quickly recovered to about the average fitness level of control chromosomes. Thus, deleterious mutations responsible for the fitness decline were efficiently selected against in relatively small populations, confirming that their effects were larger than a few percent.     

Mutation accumulation and the effect of copia insertions in Drosophila melanogaster

Repeated efforts to estimate the genomic deleterious mutation rate per generation (U) in Drosophila melanogaster have yielded inconsistent estimates ranging from 0.01 to nearly 1. We carried out a mutation-accumulation experiment with a cryopreserved control population in hopes of resolving some of the uncertainties raised by these estimates. Mutation accumulation (MA) was carried out by brother sister mating of 150 sublines derived from two inbred lines. Fitness was measured under conditions chosen to mimic the ancestral laboratory environment of these genotypes. We monitored the insertions of a transposable element, copia, that proved to accumulate at the unusually high rate of 0.24 per genome per generation in one of our MA lines. Mutational variance in fitness increased at a rate consistent with previous studies, yielding a mutational coefficient of variation greater than 3%. The performance of the cryopreserved control relative to the MA lines was inconsistent, so estimates of mutation rate by the Bateman-Mukai method are suspect. Taken at face value, these data suggest a modest decline in fitness of about 0.3% per generation. The element number of copia was a significant predictor of fitness within generations; on average, insertions caused a 0.76% loss in fitness, although the confidence limits on this estimate are wide.     

Effect of rate of inbreeding on inbreeding depression in Drosophila melanogaster

Summary This experiment was designed to study the relationship between rate of inbreeding and observed inbreeding depression of larval viability, adult fecundity and cold shock mortality in Drosophila melanogaster. Rates of inbreeding used were full-sib mating and closed lines of N=4 and N=20. Eight generations of mating in the N=20 lines, three generations in the N=4 lines and one generation of full-sib mating were synchronised to simultaneously produce individuals with an expected level of inbreeding coefficient (F) of approximately 0.25. Inbreeding depression for the three traits was significant at F=0.25. N=20 lines showed significantly less inbreeding depression than full-sib mated lines for larval viability at approximately the same level of F. A similar trend was observed for fecundity. No effect of rate of inbreeding depression was found for cold shock mortality, but this trait was measured with less precision than the other two. Natural selection acting on loci influencing larval viability and fecundity during the process of inbreeding could explain these results. Selection is expected to be more effective with slow rates of inbreeding because there are more generations and greater opportunity for selection to act before F=0.25 is reached. Selection intensities seem to have been different in the three traits measured. Selection was most intense for larval viability, less intense for fecundity and, perhaps, negligible at loci influencing cold shock mortality.     

Spontaneous and Ethyl Methanesulfonate-Induced Mutations Controlling Viability in DROSOPHILA MELANOGASTER. II. Homozygous Effect of Polygenic Mutations

Polygenic mutations affecting viability were accumulated on the second chromosome of Drosophila melanogaster by treating flies with EMS in successive generations. The treated chromosomes were later made homozygous and tested for their effects on viability by comparison of the frequency of such homozygotes with that of other genotypes in the same culture. The treated wild-type chromosomes were kept heterozygous in Pm/+ males by mating individual males in successive generations to Cy/Pm females. The number of generations of accumulation was 1 to 30 generations, depending on the concentration of EMS. A similar experiment for spontaneous polygenic mutations was also conducted by accumulating mutations for 40 generations. The lower limit of the spontaneous mutation rate of viability polygenes is estimated to be 0.06 per second chromosome per generation, which is about 12 times as high as the spontaneous recessive lethal mutation rate, 0.005. EMS-induced polygenic mutations increase linearly with the number of treated generations and with the concentration of EMS. The minimum mutation rate of viability polygenes is about 0.017 per 10(-4)m, which is only slightly larger than the lethal rate of 0.013 per 10(-4) m. The maximum estimate of the viability reduction of a single mutant is about 6 to 10 percent of the normal viability. The data are consistent with a constant average effect per mutant at all concentrations, but this is about three times as high as that for spontaneous mutants. It is obvious that one can obtain only a lower limit for the mutation rate, since some mutants may have effects so near to zero that they cannot be detected. The possibility of measuring something other than the lower limit is discussed. The ratio of the load due to detrimental mutants to that caused by lethals, the D/L ratio, is about 0.2 to 0.3 for EMS-induced mutants, as compared to about 0.5 for spontaneous mutants. This is to be expected if EMS treatment produces a large fraction of small deletions and other chromosome rearrangements which are more likely to be lethal.     

Spontaneous mutational variation for body size in Caenorhabditis elegans

We measured the impact of new mutations on genetic variation for body size in two independent sets of C. elegans spontaneous mutation-accumulation (MA) lines, derived from the N2 strain, that had been maintained by selfing for 60 or 152 generations. The two sets of lines gave broadly consistent results. The change of among-line genetic variation between cryopreserved controls and the MA lines implied that broad sense heritability increased by 0.4% per generation. Overall, MA reduced mean body size by approximately 0.1% per generation. The genome-wide rate for mutations with detectable effects on size was estimated to be approximately 0.0025 per haploid genome per generation, and their mean effects were approximately 20%. The proportion of mutations that increase body size was estimated by maximum likelihood to be no more than 20%, suggesting that the amount of mutational variation available for selection for increased size could be quite small. This hypothesis was supported by an artificial selection experiment on adult body size, started from a single highly inbred N2 individual. We observed a strongly asymmetrical response to selection of a magnitude consistent with the input of mutational variance observed in the MA experiment.     

Effects of inbreeding on economic traits of channel catfish

Summary Inbred channel catfish (Ictalurus punctatus) were produced from two generations of full-sib matings to study the effect of inbreeding on reproduction, growth and survival. A randomly mated control line was propagated from the same base population to be used for the evaluation of the inbred fish. First generation inbred (I1) and control (C1) lines comprised five full-sib families each. Second generation inbred (I2) and control (C2) lines were produced by mating each male catfish from the I1 or C1 line to two females in sequence, one from the I1 and one from the C1 line. The design also produced two reciprocal outcross lines to be compared to their contemporary inbred and control lines. The coefficient of inbreeding for the inbred line increased from 0.25 in generation 1 to 0.375 in generation 2. The inbreeding coefficient was zero for all other lines. The resulting fish were performance tested in two locations, Tifton, Georgia and Auburn, Alabama and no genotype-environment interactions occurred. Results indicated that one generation of inbreeding increased number of days required for eggs to hatch by 21%, but did not significantly influence spawn weight or hatchability score. However, inbred females produced more eggs/kg body weight than control females. Two generations of full-sib mating in Georgia did not depress weight when expressed as a deviation to random controls but was depressed 13–16% when expressed as a deviation to half-sib out-crosses. Second generation inbreds produced in Alabama exhibited a 19% depression for growth rate when compared to either random or half-sib outcross controls. Survival rates at various age intervals was not decreased by inbreeding. The amount of inbreeding depression varied among families and between sexes.This study was supported by State and Hatch Funds allocated to the Georgia and Alabama Agricultural Experiment Stations     

Analysis of response to 20 generations of selection for body composition in mice: fit to infinitesimal model assumptions

Data were analysed from a divergent selection experiment for an indicator of body composition in the mouse, the ratio of gonadal fat pad to body weight (GFPR). Lines were selected for 20 generations for fat (F), lean (L) or were unselected (C), with three replicates of each. Selection was within full-sib families, 16 families per replicate for the first seven generations, eight subsequently. At generation 20, GFPR in the F lines was twice and in the L lines half that of C. A log transformation removed both asymmetry of response and heterogeneity of variance among lines, and so was used throughout. Estimates of genetic variance and heritability (approximately 50%) obtained using REML with an animal model were very similar, whether estimated from the first few generations of selection, or from all 20 generations, or from late generations having fitted pedigree. The estimates were also similar when estimated from selected or control lines. Estimates from REML also agreed with estimates of realised heritability. The results all accord with expectations under the infinitesimal model, despite the four-fold changes in mean. Relaxed selection lines, derived from generation 20, showed little regression in fatness after 40 generations without selection.     

Temporal uniformity of the spontaneous mutational variance of quantitative traits in Drosophila melanogaster

Spontaneous mutations were allowed to accumulate over 209 generations in more than 100 lines, all of them independently derived from a completely homozygous population of Drosophila melanogaster and subsequently maintained under strong inbreeding (equivalent to full-sib mating). Traits scored were: abdominal (AB) and sternopleural (ST) bristle number, wing length (WL) and egg-to-adult viability (V). On two occasions--early (generations 93-122) and late (generations 169-209)--ANOVA estimates of the mutational variance and the mutational line x generation interaction variance were obtained. Mutational heritabilities of morphological traits ranged from 2 x 10(-4) to 2 x 10(-3) and the mutational coefficient of variation of viability was 0.01. For AB, WL and V, temporal uniformity of the mutational variance was observed. However, a fluctuation of the mutational heritability of ST was detected and could be ascribed to random genotype x environment interaction.     

Fitness decline in spontaneous mutation accumulation lines of Caenorhabditis elegans with varying effective population sizes

The rate and fitness effects of new mutations have been investigated by mutation accumulation (MA) experiments in which organisms are maintained at a constant minimal population size to facilitate the accumulation of mutations with minimal efficacy of selection. We evolved 35 MA lines of Caenorhabditis elegans in parallel for 409 generations at three population sizes (N = 1, 10, and 100), representing the first spontaneous long-term MA experiment at varying population sizes with corresponding differences in the efficacy of selection. Productivity and survivorship in the N = 1 lines declined by 44% and 12%, respectively. The average effects of deleterious mutations in N = 1 lines are estimated to be 16.4% for productivity and 11.8% for survivorship. Larger populations (N = 10 and 100) did not suffer a significant decline in fitness traits despite a lengthy and sustained regime of consecutive bottlenecks exceeding 400 generations. Together, these results suggest that fitness decline in very small populations is dominated by mutations with large deleterious effects. It is possible that the MA lines at larger population sizes contain a load of cryptic deleterious mutations of small to moderate effects that would be revealed in more challenging environments.     

THE RATE AND EFFECTS DISTRIBUTION OF VIABILITY MUTATION IN DROSOPHILA: MINIMUM DISTANCE ESTIMATION

The empirical distribution of the mean viability of mutation accumulation lines, obtained from three published experiments, was analyzed using minimum-distance estimation. In two cases (Mukai et al. 1972; Ohnishi 1977), mutations were allowed to accumulate in copies of chromosome II protected from natural selection and recombination. In the other one (Fernández and López-Fanjul 1996), they accumulated in inbred lines derived from an isogenic stock. In contrast with currently accepted hypotheses, we consistently estimated low (about 0.01) genomic viability mutation rates, λ, and a small kurtosis of the distribution of mutational effects on viability (a) in the three datasets. Minimum-distance estimates of the per-generation mean viability change due to mutation (λE[a]) were also obtained. These were very similar for both chromosomal datasets, their absolute values being about five times smaller than estimates obtained from the observed change in mean viability during the mutation process. It must be noted that, in both experiments, viability was measured relative to the Cy chromosome of a Cy/Pm stock. Thus, an unnoticed viability increase in this Cy chromosome may have resulted in overestimation of the mean viability reduction in the lines. In parallel, minimum-distance estimation of λE(a) from inbred lines data (where the selective pressure during the accumulation process was larger) was even somewhat smaller, in absolute value, and very close to the estimate obtained by comparing the mean viability of the lines with that of the control isogenic line. The evolutionary importance of these results, as well as their relevance to the solution of the mutational load paradox, is discussed.     

Increased Selection Response in Larger Populations. II. Selection for Ethanol Vapor Resistance in Drosophila Melanogaster at Two Population Sizes

The effect of large population size on selection response was investigated using Drosophila melanogaster, with four "small" lines of 160 selected parents/generation compared to two "large" lines of 1,600 selected parents/generation. All lines were selected under similar conditions at a selection intensity of approximately 0.55 standard deviations, for 65 generations, for increased ethanol vapor resistance (measured in minutes required to become anesthetized). Two unselected control lines of 320 parents/generation were also maintained. A significant effect of population size was found. The final treatment means and standard errors were: 27.91 +/- 1.28 min (two "large" lines); 19.40 +/- 1.54 min (four "small" lines); and 4.98 +/- 0.35 min (two control lines). To estimate the mutation rate for the trait, two isogenic lines of about 400 selected parents were selected for 29 generations. The mean increase in additive genetic variance per generation was 0.0009 times the initial environmental variance of the outbred lines. This is comparable to other reported mutation rates. Mutation can explain part of the difference in evolved resistance between treatments, but it appears that even at rather large population sizes, a large difference in long-term response can be obtained in larger outbred lines, from more complete utilization of the initial genetic variation.     

A Case of High Rate of Spontaneous Mutation Affecting Viability in DROSOPHILA MELANOGASTER

A spontaneous lethal mutation rate approximately twenty to thirty times greater than normal has been discovered in second and third chromosomes derived from an irradiated isogenic line and paired with marked inversion chromosomes. Mutations resulting in reductions of viability of varying magnitude short of complete lethality apparently also occur at a very high rate in the third but not in the second chromosome. The pattern of accumulation of lethal mutations over several generations and viability frequency distributions within generations have been studied in a number of independent experiments. High mutation rate occurs in heterozygous isogenic-derived second and third chromosomes, either together or apart, irrespective of the genetic constitution of nonhomologous chromosomes. High mutation rates were not observed using the same methods with chromosomes of an inbred line from a different source. The possible mechanisms responsible for these results are discussed.     

The effect of spontaneous mutations on competitive ability

Understanding the impact of spontaneous mutations on fitness has many theoretical and practical applications in biology. Although mutational effects on individual morphological or life‐history characters have been measured in several classic genetic model systems, there are few estimates of the rate of decline due to mutation for complex fitness traits. Here, we estimate the effects of mutation on competitive ability, an important complex fitness trait, in a model system for ecological and evolutionary genomics, Daphnia. Competition assays were performed to compare fitness between mutation‐accumulation (MA) lines and control lines from eight different genotypes from two populations of Daphnia pulicaria after 30 and 65 generations of mutation accumulation. Our results show a fitness decline among MA lines relative to controls as expected, but highlight the influence of genomic background on this effect. In addition, in some assays, MA lines outperform controls providing insight into the frequency of beneficial mutations.     

Heterosis for viability,fecundity, and male fertility in Drosophila melanogaster: comparison of mutational and standing variation.

If genetic variation for fitness traits in natural populations ("standing" variation) is maintained by recurrent mutation, then quantitative-genetic properties of standing variation should resemble those of newly arisen mutations. One well-known property of standing variation for fitness traits is inbreeding depression, with its converse of heterosis or hybrid vigor. We measured heterosis for three fitness traits, pre-adult viability, female fecundity, and male fertility, among a set of inbred Drosophilia melanogaster lines recently derived from the wild, and also among a set of lines that had been allowed to accumulate spontaneous mutations for over 200 generations. The inbred lines but not the mutation-accumulation (MA) lines showed heterosis for pre-adult viability. Both sets of lines showed heterosis for female fecundity, but heterosis for male fertility was weak or absent. Crosses among a subset of the MA lines showed that they were strongly differentiated for male fertility, with the differences inherited in autosomal fashion; the absence of heterosis for male fertility among the MA lines was therefore not caused by an absence of mutations affecting this trait. Crosses among the inbred lines also gave some, albeit equivocal, evidence for male fertility variation. The contrast between the results for female fecundity and those for male fertility suggests that mutations affecting different fitness traits may differ in their average dominance properties, and that such differences may be reflected in properties of standing variation. The strong differentiation among the MA lines in male fertility further suggests that mutations affecting this trait occur at a high rate.     

Effects of Population Size and Selection Intensity on Short-Term Response to Selection for Postweaning Gain in Mice

The effects of population size and selection intensity on the mean response was examined after 14 generations of within full-sib family selection for postweaning gain in mice. Population sizes of 1, 2, 4, 8 and 16 pair matings were each evaluated at selection intensities of 100% (control), 50% and 25% in a replicated experiment. Selection response per generation increased as selection intensity increased. Selection response and realized heritability tended to increase with increasing population size. Replicate variability in realized heritability was large at population sizes of 1, 2 and 4 pairs. Genetic drift was implicated as the primary factor causing the reduced response and lowered repeatability at the smaller population sizes. Lines with intended effective population sizes of 62 yielded larger selection responses per unit selection differential than lines with effective population sizes of 30 or less.     

Reduced Genetic Load Revealed by Slow Inbreeding in Drosophila Melanogaster

The rate of decline in reproductive fitness in populations of Drosophilia melanogaster inbred at an initial rate of ~1% per generation has been investigated under both competitive and noncompetitive conditions. Breeding population size was variable in the inbred lines with an estimated harmonic mean of 66.7 +/- 2.2. Of the 60 lines maintained without reserves, 75% survived a period of 210 generations of slow inbreeding and were then rapidly inbred by full-sib mating to near-homozygosity. The initial rate of inbreeding was estimated to be 0.96 +/- 0.16% per generation, corresponding to an effective population size of ~50. However, the rate of inbreeding declined significantly with time to average only 0.52 +/- 0.08% per generation over the 210 generation period, most likely due to associative overdominance built up by genetic sampling and selection in the small populations. The total inbreeding depression in fitness was estimated to be 87 +/- 3% for competitive ability and 27 +/- 5% for fitness under uncrowded conditions, corresponding to rates of decline of 2.0 +/- 0.3 and 0.32 +/- 0.07%, respectively, per 1% increase in the inbreeding coefficient. The frequency of lethal second chromosomes in the resultant near-homozygous lines was of the order of 5%, lethal free second chromosomes showed a mean viability under both crowded and uncrowded conditions of ~95%, and their population cage fitness was 60% that of Cy/+ heterozygotes. It can be concluded that homozygous genotypes from which deleterious genes of major effect have been eliminated during slow inbreeding may show far less depression in reproductive fitness than suggested by earlier studies of wild chromosome homozygotes. The loss in fitness due to homozygosity throughout the entire genome may be as little as 85-90% under competitive conditions, and 25-30% in an optimal environment.