Spontaneous Mutational Variances and Covariances for Fitness-Related Traits in Drosophila Melanogaster

Starting from a completely homozygous population of Drosophila melanogaster, 176 lines were derived and independently maintained by a single brother-sister mating per generation. Three fitness-related traits were considered (fecundity, egg-to-pupa and pupa-to-adult viabilities). Mutational heritabilities of these traits and genetic correlations between all possible pairs were calculated from the between line divergence (codivergence), after 104-106 generations of mutation accumulation. Mutational heritabilities ranged from 0.60 X 10(-3) to 0.82 X 10(-3) and correlations from -0.11 to 0.25. These values are likely to be underestimates due to selection against deleterious mutations. The distribution of the means of the lines was asymmetric, positive for fecundity and negative for both viability components. The coefficients of asymmetry are also likely to be biased, again due to selection. Extreme lines from the two tails of the distribution were examined in detail. Homozygous line effects were all negative for viability traits but predominantly positive for fecundity, indicating the fixation of mutations with positive effects on the latter. Corresponding heterozygous line effects showed a variable degree of dominance.     

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.     

The Effects of Spontaneous Mutation on Quantitative Traits. I. Variances and Covariances of Life History Traits

We have accumulated spontaneous mutations in the absence of natural selection in Drosophila melanogaster by backcrossing 200 heterozygous replicates of a single high fitness second chromosome to a balancer stock for 44 generations. At generations 33 and 44 of accumulation, we extracted samples of chromosomes and assayed their homozygous performance for female fecundity early and late in adult life, male and female longevity, male mating ability early and late in adult life, productivity (a measure of fecundity times viability) and body weight. The variance among lines increased significantly for all traits except male mating ability and weight. The rate of increase in variance was similar to that found in previous studies of egg-to-adult viability, when calculated relative to trait means. The mutational correlations among traits were all strongly positive. Many correlations were significantly different from 0, while none was significantly different from 1. These data suggest that the mutation-accumulation hypothesis is not a sufficient explanation for the evolution of senescence in D. melanogaster. Mutation-selection balance does seem adequate to explain a substantial proportion of the additive genetic variance for fecundity and longevity.     

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.     

Spontaneous mutational correlations for life-history, morphological and behavioral characters in Caenorhabditis elegans

The pattern of mutational covariance among traits plays a central, but largely untested, role in many theories in evolutionary genetics. Here we estimate the pattern of phenotypic, environmental, and mutational correlations for a set of life-history, behavioral, and morphological traits using 67 self-fertilizing lines of Caenorhabditis elegans, each having independently experienced an average of 370 generations of spontaneous mutation accumulation. Bivariate relationships of mutational effects indicate the existence of extensive pleiotropy. We find that mutations may tend to produce manifold effects on suites of functionally related traits; however, our data do not support the idea of completely parcelated pleiotropy, in which functional units are separately affected by mutations. Positive net phenotypic and mutational correlations are common for life-history traits, with environmental correlations being comparatively smaller and of the same sign for most pairs of traits. Observed mutational correlations are shown to be higher than those produced by the chance accumulation of nonpleiotropic mutations in the same lines.     

Selection for distinct gene expression properties favours the evolution of mutational robustness in gene regulatory networks

Mutational robustness is a genotype's tendency to keep a phenotypic trait with little and few changes in the face of mutations. Mutational robustness is both ubiquitous and evolutionarily important as it affects in different ways the probability that new phenotypic variation arises. Understanding the origins of robustness is specially relevant for systems of development that are phylogenetically widespread and that construct phenotypic traits with a strong impact on fitness. Gene regulatory networks are examples of this class of systems. They comprise sets of genes that, through cross‐regulation, build the gene activity patterns that define cellular responses, different tissues or distinct cell types. Several empirical observations, such as a greater robustness of wild‐type phenotypes, suggest that stabilizing selection underlies the evolution of mutational robustness. However, the role of selection in the evolution of robustness is still under debate. Computer simulations of the dynamics and evolution of gene regulatory networks have shown that selection for any gene activity pattern that is steady and self‐sustaining is sufficient to promote the evolution of mutational robustness. Here, I generalize this scenario using a computational model to show that selection for different aspects of a gene activity phenotype increases mutational robustness. Mutational robustness evolves even when selection favours properties that conflict with the stationarity of a gene activity pattern. The results that I present support an important role for stabilizing selection in the evolution of robustness in gene regulatory networks.     

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.     

Quantitative genetics of ovariole number in Drosophila melanogaster. II. Mutational variation and genotype-environment interaction.

The rare alleles model of mutation-selection balance (MSB) hypothesis for the maintenance of genetic variation was evaluated for two quantitative traits, ovariole number and body size. Mutational variances (VM) for these traits, estimated from mutation accumulation lines, were 4.75 and 1.97 x 10(-4) times the environmental variance (VE), respectively. The mutation accumulation lines were studied in three environments to test for genotype x environment interaction (GEI) of new mutations; significant mutational GEI was found for both traits. Mutations for ovariole number have a quadratic relationship with competitive fitness, suggesting stabilizing selection for the trait; there is no significant correlation between mutations for body size and competitive fitness. Under MSB, the ratio of segregating genetic variance, VG, to mutational variance, VM, estimates the inverse of the selection coefficient against a heterozygote for a new mutation. Estimates of VG/VM for ovariole number and body size were both approximately 1.1 x 10(4). Thus, MSB can explain the level of variation, if mutations affecting these traits are under very weak selection, which is inconsistent with the empirical observation of stabilizing selection, or if the estimate of VM is biased downward by two orders of magnitude. GEI is a possible alternative explanation.     

The mutational structure of metabolism in Caenorhabditis elegans

A properly functioning organism must maintain metabolic homeostasis. Deleterious mutations degrade organismal function, presumably at least in part via effects on metabolic function. Here we present an initial investigation into the mutational structure of the Caenorhabditis elegans metabolome by means of a mutation accumulation experiment. We find that pool sizes of 29 metabolites vary greatly in their vulnerability to mutation, both in terms of the rate of accumulation of genetic variance (the mutational variance, VM) and the rate of change of the trait mean (the mutational bias, ΔM). Strikingly, some metabolites are much more vulnerable to mutation than any other trait previously studied in the same way. Although we cannot statistically assess the strength of mutational correlations between individual metabolites, principal component analysis provides strong evidence that some metabolite pools are genetically correlated, but also that there is substantial scope for independent evolution of different groups of metabolites. Averaged over mutation accumulation lines, PC3 is positively correlated with relative fitness, but a model in which metabolites are uncorrelated with fitness is nearly as good by Akaike's Information Criterion.     

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.     

The evolution of age-specific mortality rates in Drosophila melanogaster: genetic divergence among unselected lines.

Age-specific effects of spontaneous mutations on mortality rates in Drosophila are inferred from three large demographic experiments. Data were collected from inbred lines that were allowed to accumulate spontaneous mutations for 10, 19, and 47 generations. Estimates of age-specific mutational variance for mortality were based on data from all three experiments, totalling approximately 225,000 flies, using a model developed for genetic analysis of age-dependent traits (the character process model). Both within- and among-generation analyses suggest that the input of genetic variance is greater for early life mortality rates than for mortality at older ages. In females, age-specific mutational variances ranged over an order of magnitude from 5.96 x 10(-3) at 2 wk posteclosion to 0.02 x 10(-3) at 7 wk. The male data show a similar pattern. Age-specific genetic variances were substantially less at generation 47 than at generation 19-an unexplained observation that is likely due to block effects. Mutational correlations among mortality rates at different ages tend to increase with the accumulation of new mutations. Comparison of the mutation-accumulation lines at generations 19 and 47 with their respective control lines suggests little age-specific mutational bias.     

The short-term and long-term effects of parental age in the bean weevil (Acanthoscelides obtectus)

We examined the influence of parental age on life history traits of their offspring in the lines of bean weevil that have evolved different rates of senescence. Measurements included preadult traits (egg size, embryonic developmental time, total preadult developmental time, preadult viability) and adult traits (body weight, total realized fecundity of females, first day of egg laying, early fecundity, late fecundity and longevity). The negative parental age effects were observed for all traits except for the early and total realized fecundity. We did not detect statistically significant line×parental age interactions for either preadult- or adult-survival, so offspring survival did not change with parental age after selection for early vs. late reproduction. It seems that selection acting on the quality of offspring produced by parents of different ages has not been responsible for the evolution of senescence in bean weevil. This revised version was published online in August 2006 with corrections to the Cover Date.     

Environment dependence of mutational parameters for viability in Drosophila melanogaster

The genomic rate of mildly deleterious mutations (U) figures prominently in much evolutionary and ecological theory. In Drosophila melanogaster, estimates of U have varied widely, from <0.1 to nearly 1 per zygote. The source of this variation is unknown, but could include differences in the conditions used for assaying fitness traits. We examined how assay conditions affect estimates of the rates and effects of viability-depressing mutations in two sets of lines with accumulated spontaneous mutations on the second chromosome. In each set, the among-line variance in egg-to-adult viability was significantly greater when viability was assayed using a high parental density than when it was assayed using a low density. In contrast, the proportional decline in viability due to new mutations did not differ between densities. Two other manipulations, lowering the temperature and adding ethanol to the medium, had no significant effects on either the mean decline or among-line variance. Cross-environment genetic correlations in viability were generally close to one, implying that most mutations reduced viability in all environments. Using data from the low-density, lower-bound estimates of U approached the classic, high values of Mukai and Ohnishi; at the high density, U estimates were similar to recently reported low values. The difference in estimated mutation rates, taken at face value, would imply that many mutations affected fitness at low density but not at high density, but this is shown to be incompatible with the observed high cross-environment correlations. Possible reasons for this discrepancy are discussed. Regardless of the interpretation, the results show that assay conditions can have a large effect on estimates of mutational parameters for fitness traits.     

Quantifying natural seasonal variation in mutation parameters with mutation accumulation lines

Mutations create novel genetic variants, but their contribution to variation in fitness and other phenotypes may depend on environmental conditions. Furthermore, natural environments may be highly heterogeneous. We assessed phenotypes associated with survival and reproductive success in over 30,000 plants representing 100 mutation accumulation lines of Arabidopsis thaliana across four temporal environments at a single field site. In each of the four assays, environmental variance was substantially larger than mutational variance. For some traits, whether mutational variance was significantly varied between seasons. The founder genotype had mean trait values near the mean of the distribution of the mutation accumulation lines in all field experiments. New mutations also contributed more phenotypic variation than would be predicted, given phenotypic and sequence‐level divergence among natural populations of A. thaliana. The combination of large environmental variance with a mean effect of mutation near zero suggests that mutations could contribute substantially to standing genetic variation.     

Simultaneous Estimation of Additive and Mutational Genetic Variance in an Outbred Population of Drosophila serrata

How new mutations contribute to genetic variation is a key question in biology. Although the evolutionary fate of an allele is largely determined by its heterozygous effect, most estimates of mutational variance and mutational effects derive from highly inbred lines, where new mutations are present in homozygous form. In an attempt to overcome this limitation, middle-class neighborhood (MCN) experiments have been used to assess the fitness effect of new mutations in heterozygous form. However, because MCN populations harbor substantial standing genetic variance, estimates of mutational variance have not typically been available from such experiments. Here we employ a modification of the animal model to analyze data from 22 generations of Drosophila serrata bred in an MCN design. Mutational heritability, measured for eight cuticular hydrocarbons, 10 wing-shape traits, and wing size in this outbred genetic background, ranged from 0.0006 to 0.006 (with one exception), a similar range to that reported from studies employing inbred lines. Simultaneously partitioning the additive and mutational variance in the same outbred population allowed us to quantitatively test the ability of mutation-selection balance models to explain the observed levels of additive and mutational genetic variance. The Gaussian allelic approximation and house-of-cards models, which assume real stabilizing selection on single traits, both overestimated the genetic variance maintained at equilibrium, but the house-of-cards model was a closer fit to the data. This analytical approach has the potential to be broadly applied, expanding our understanding of the dynamics of genetic variance in natural populations.     

Whole-genome effects of ethyl methanesulfonate-induced mutation on nine quantitative traits in outbred Drosophila melanogaster

We induced mutations in Drosophila melanogaster males by treating them with 21.2 mm ethyl methanesulfonate (EMS). Nine quantitative traits (developmental time, viability, fecundity, longevity, metabolic rate, motility, body weight, and abdominal and sternopleural bristle numbers) were measured in outbred heterozygous F3 (viability) or F2 (all other traits) offspring from the treated males. The mean values of the first four traits, which are all directly related to the life history, were substantially affected by EMS mutagenesis: the developmental time increased while viability, fecundity, and longevity declined. In contrast, the mean values of the other five traits were not significantly affected. Rates of recessive X-linked lethals and of recessive mutations at several loci affecting eye color imply that our EMS treatment was equivalent to approximately 100 generations of spontaneous mutation. If so, our data imply that one generation of spontaneous mutation increases the developmental time by 0.09% at 20 degrees and by 0.04% at 25 degrees, and reduces viability under harsh conditions, fecundity, and longevity by 1.35, 0.21, and 0.08%, respectively. Comparison of flies with none, one, and two grandfathers (or greatgrandfathers, in the case of viability) treated with EMS did not reveal any significant epistasis among the induced mutations.     

Properties of ethylmethane sulfonate-induced mutations affecting life-history traits in Caenorhabditis elegans and inferences about bivariate distributions of mutation effects

The homozygous effects of ethylmethane sulfonate (EMS)-induced mutations in Caenorhabditis elegans are compared across life-history traits. Mutagenesis has a greater effect on early than late reproductive output, since EMS-induced mutations tend to cause delayed reproduction. Mutagenesis changes the mean and variance of longevity much less than reproductive output traits. Mutations that increase total or early productivity are not detected, but the net effect of mutations is to increase and decrease late productivity to approximately equal extents. Although most mutations decrease longevity, a mutant line with increased longevity was found. A flattening of mortality curves with age is noted, particularly in EMS lines. We infer that less than one-tenth of mutations that have fitness effects in natural conditions are detected in the laboratory, and such mutations have moderately large effects ( approximately 20% of the mean). Mutational correlations for life-history traits are strong and positive. Correlations between early or late productivity and longevity are of similar magnitude. We develop a maximum-likelihood procedure to infer bivariate distributions of mutation effects. We show that strong mutation-induced genetic correlations do not necessarily imply strong directional correlations between mutational effects, since correlation is also generated by lines carrying different numbers of mutations.     

Susceptibility of the male fitness phenotype to spontaneous mutation

Adult reproductive success can account for a large fraction of male fitness, however, we know relatively little about the susceptibility of reproductive traits to mutation-accumulation (MA). Estimates of the mutational rate of decline for adult fitness and its components are controversial in Drosophila melanogaster, and post-copulatory performance has not been examined. We therefore separately measured the consequences of MA for total male reproductive success and its major pre-copulatory and post-copulatory components: mating success and sperm competitive success. We also measured juvenile viability, an important fitness component that has been well studied in MA experiments. MA had strongly deleterious effects on both male viability and adult fitness, but the latter declined at a much greater rate. Mutational pressure on total fitness is thus much greater than would be predicted by viability alone. We also noted a significant and positive correlation between all adult traits and viability in the MA lines, suggesting pleiotropy of mutational effect as required by 'good genes' models of sexual selection.     

What affects the predictability of evolutionary constraints using a G‐matrix? The relative effects of modular pleiotropy and mutational correlation

Phenotypic traits do not always respond to selection independently from each other and often show correlated responses to selection. The structure of a genotype‐phenotype map (GP map) determines trait covariation, which involves variation in the degree and strength of the pleiotropic effects of the underlying genes. It is still unclear, and debated, how much of that structure can be deduced from variational properties of quantitative traits that are inferred from their genetic (co) variance matrix ( G ‐matrix). Here we aim to clarify how the extent of pleiotropy and the correlation among the pleiotropic effects of mutations differentially affect the structure of a G ‐matrix and our ability to detect genetic constraints from its eigen decomposition. We show that the eigenvectors of a G ‐matrix can be predictive of evolutionary constraints when they map to underlying pleiotropic modules with correlated mutational effects. Without mutational correlation, evolutionary constraints caused by the fitness costs associated with increased pleiotropy are harder to infer from evolutionary metrics based on a G ‐matrix's geometric properties because uncorrelated pleiotropic effects do not affect traits' genetic correlations. Correlational selection induces much weaker modular partitioning of traits' genetic correlations in absence then in presence of underlying modular pleiotropy.     

Perturbations in the impact of mutational activity on Vλ genes in systemic lupus erythematosus

To assess the impact of somatic hypermutation and selective influences on the Vλ light chain repertoire in systemic lupus erythematosus (SLE), the frequency and pattern of mutations were analyzed in individual CD19⁺ B cells from a patient with previously undiagnosed SLE. The mutational frequency of nonproductive and productive rearrangements in the SLE patient was greater (3.1 × 10^-2 vs 3.4 × 10^-2, respectively) than that in normal B cells (1.2 × 10^-2 vs 2.0 × 10^-2, both P < 0.001). The frequencies of mutated rearrangements in both the nonproductive and productive repertoires were significantly higher in the patient with SLE than in normal subjects. Notably, there were no differences in the ratio of replacement to silent (R/S) mutations in the productive and nonproductive repertoires of the SLE patient, whereas the R/S ratio in the framework regions of productive rearrangements of normal subjects was reduced, consistent with active elimination of replacement mutations in this region. The pattern of mutations was abnormal in the SLE patient, with a significant increase in the frequency of G mutations in both the productive and nonproductive repertoires. As in normal subjects, however, mutations were found frequently in specific nucleotide motifs, the RGYW/WRCY sequences, accounting for 34% (nonproductive) and 46% (productive) of all mutations. These data are most consistent with the conclusion that in this SLE patient, the mutational activity was markedly greater than in normal subjects and exhibited some abnormal features. In addition, there was decreased subsequent positive or negative selection of mutations. The enhanced and abnormal mutational activity along with disturbances in selection may play a role in the emergence of autoreactivity in this patient with SLE.