Mutation-selection balance accounting for genetic variation for viability inDrosophila melanogasteras deduced from an inbreeding and artificial selection experiment

Stage- and thermal-specific genetic architecture for preadult viability in natural populations of Drosophila melanogaster

Studying the processes affecting variation for preadult viability is essential to understand the evolutionary trajectories followed by natural populations. This task requires focusing on the complex nature of the phenotype-genotype relationship by taking into account usually neglected aspects of the phenotype and recognizing the modularity between different ontogenetic stages. Here, we describe phenotypic variability for viability during the larval and pupal stages in lines derived from three natural populations of Drosophila melanogaster, as well as the variability for phenotypic plasticity and canalization at two different rearing temperatures. The results indicate that the three populations present significant phenotypic differences for preadult viability. Furthermore, distinct aspects of the phenotype (means, plasticity, canalization, plasticity of canalization) are affected by different genetic bases underlying changes in viability in a stage- and environment-specific manner. These findings explain the generalized maintenance of genetic variability for this fitness trait.

Effects of genetic drift and gene flow on the selective maintenance of genetic variation

Explanations for the genetic variation ubiquitous in natural populations are often classified by the population-genetic processes they emphasize: natural selection or mutation and genetic drift. Here we investigate models that incorporate all three processes in a spatially structured population, using what we call a construction approach, simulating finite populations under selection that are bombarded with a steady stream of novel mutations. As expected, the amount of genetic variation compared to previous models that ignored the stochastic effects of drift was reduced, especially for smaller populations and when spatial structure was most profound. By contrast, however, for higher levels of gene flow and larger population sizes, the amount of genetic variation found after many generations was greater than that in simulations without drift. This increased amount of genetic variation is due to the introduction of slightly deleterious alleles by genetic drift and this process is more efficient when migration load is higher. The incorporation of genetic drift also selects for fitness sets that exhibit allele-frequency equilibria with larger domains of attraction: they are "more stable." Moreover, the finiteness of populations strongly influences levels of local adaptation, selection strength, and the proportion of allele-frequency vectors that can be distinguished from the neutral expectation.

Simulating the maintenance of a rare fish morph experiencing negative frequency dependent selection

Empirical work assessing the maintenance of rare genotypes in natural populations is difficult over very long time scales. Skirting this problematic issue is possible with theory and simulations. Major theoretical constructs, including mutation-selection balance and balancing selection, explain the theoretical maintenance of rare genotypes, and the occurrence of multiple, rare genotypes over time. Additionally, numerical simulations are valuable tools for assessing evolving biological systems because they allow for monitoring systems over long time scales, as well as for controlling model parameters, thus contributing to the exploration of system dynamics that cannot be assessed in nature. Here we employed numerical simulations to explore the importance of several biological factors that contribute to the maintenance of a fish color-pattern polymorphism. We present a numerical model of a two-morph fish polymorphism that allowed us to test the sensitivity of the rare morph's persistence and the population's stability to multiple parameters. Our simulations ran over 10,000 years (where one year is approximately one generation) and demonstrated the maintenance of a stable polymorphism with a rare morph which persisted at a frequency of ~10(-2), which is in-fact the frequency of the rare, mottled black mosquitofish morph in natural populations. This pigmentation polymorphism is stable, independent of changes in population size, but can be destabilized with very high predation when coupled with very low birth rates. Employing models with empirical fitness estimates is a valuable tool for monitoring rare vertebrate morphs in nature, however few studies exist that have accomplished this task. Our approach can be adapted for modeling rare morphs (particularly in additional live-bearing fishes like sailfin mollies) that also harbor rare, pigmentation morphs within large populations.

The maintenance of genetic variation for oviposition rate in two-spotted spider mites: inferences from artificial selection

Despite the directional selection acting on life-history traits, substantial amounts of standing variation for these traits have frequently been found. This variation may result from balancing selection (e.g., through genetic trade-offs) or from mutation-selection balance. These mechanisms affect allele frequencies in different ways: Under balancing selection alleles are maintained at intermediate frequencies, whereas under mutation-selection balance variation is generated by deleterious mutations and removed by directional selection, which leads to asymmetry in the distribution of allele frequencies. To investigate the importance of these two mechanisms in maintaining heritable variation in oviposition rate of the two-spotted spider mite, we analyzed the response to artificial selection. In three replicate experiments, we selected for higher and lower oviposition rate, compared to control lines. A response to selection only occurred in the downward direction. Selection for lower oviposition rate did not lead to an increase in any other component of fitness, but led to a decline in female juvenile survival. The results suggest standing variation for oviposition rate in this population consists largely of deleterious alleles, as in a mutation-selection balance. Consequently, the standing variation for this trait does not appear to be indicative of its adaptive potential.

Genetic adaptation to captivity in species conservation programs

As wild environments are often inhospitable, many species have to be captive-bred to save them from extinction. In captivity, species adapt genetically to the captive environment and these genetic adaptations are overwhelmingly deleterious when populations are returned to wild environments. I review empirical evidence on (i) the genetic basis of adaptive changes in captivity, (ii) factors affecting the extent of genetic adaptation to captivity, and (iii) means for minimizing its deleterious impacts. Genetic adaptation to captivity is primarily due to rare alleles that in the wild were deleterious and partially recessive. The extent of adaptation to captivity depends upon selection intensity, genetic diversity, effective population size and number of generation in captivity, as predicted by quantitative genetic theory. Minimizing generations in captivity provides a highly effective means for minimizing genetic adaptation to captivity, but is not a practical option for most animal species. Population fragmentation and crossing replicate captive populations provide practical means for minimizing the deleterious effects of genetic adaptation to captivity upon populations reintroduced into the wild. Surprisingly, equalization of family sizes reduces the rate of genetic adaptation, but not the deleterious impacts upon reintroduced populations. Genetic adaptation to captivity is expected to have major effects on reintroduction success for species that have spent many generations in captivity. This issue deserves a much higher priority than it is currently receiving.

Testing alternative methods for purging genetic load using the housefly (Musca domestica L.)

When a population faces long-term inbreeding, artificial selection, in principle, can enhance natural selection processes for purging the exposed genetic load. However, strong purge pressures might actually decrease fitness through the inadvertent fixation of deleterious alleles and allelic combinations. We tested lines of the housefly (Musca domestica L.) for the effectiveness of artificial selection to promote the adaptation to small population size. Specifically, replicate populations were held at average census sizes of 54 for nine generations or 30 for 14 generations while being subjected to artificial selection pressure for increased fitness in overall mating propensity (i.e., the proportion of virgin male-female pairs initiating copulation within 30 min), while also undergoing selection to create differences among lines in multivariate components of courtship performance. In the 14-generation experiment, a subset of the lines were derived from a founder-flush population (i.e., derived from three male-female pairs). In both experiments, we also maintained parallel non-selection lines to assess the potential for natural purging through serial inbreeding alone. Sub-populations derived from a stock newly derived from the wild responded to artificial selection for increased mating propensity, but only in the short-term, with eventual rebounds back to the original levels. Serial inbreeding in these lines simply reduced mating propensity. In sub-populations derived from the same base population, but 36 generations later, both artificial selection and serial inbreeding increased mating propensity, but mainly to restore the level found upon establishment in the laboratory. Founder-flush lines responded as well as the non-bottlenecked controls, so we base our major conclusions on the comparisons between fresh-caught and long-term laboratory stocks. We suggest that the effectiveness of the alternative purge protocols depended upon the amount of genetic load already exposed, such that prolonged periods of relaxed or altered selection pressures of the laboratory rendered a population more responsive to purging protocols.

Reversed selection responses in small populations of the housefly (Musca domestica L.)

We compared the efficacy of artificial and natural selection processes in purging the genetic load of perpetually small populations. We subjected replicate lines of the housefly (Musca domestica L.), recently derived from the wild, to artificial selection for increased mating propensity (i.e., the proportion of male-female pairs initiating copulation within 30 min) in efforts to cull out the inbreeding depression effects of long-term small population size (as determined by a selection protocol for increased assortative mating). We also maintained parallel non-selection lines for assessing the spontaneous purge of genetic load due to inbreeding alone. We thus evaluated the fitness of artificially and 'naturally' purging populations held at census sizes of 40 individuals over the course of 18 generations. We found that the artificially selected lines had significant increases in mating propensity (up to 46% higher from the beginning of the protocol) followed by reversed selection responses back to the initial levels, resulting in non-significant heritabilities. Nevertheless, the 'naturally' selected lines had significantly lower fitness overall (a 28% reduction from the beginning of the protocol), although lower effective population sizes could have contributed to this effect. We conclude that artificial selection bolstered fitness, but only in the short-term, because the inadvertent fixation of extant genetic load later resulted in pleiotropic fitness declines. Still, the short-term advantage of the selection protocol likely contributed to the success of the speciation experiment since our recently-derived housefly populations are particularly vulnerable to inbreeding depression effects on mating behavior.

Genetic load in subdivided populations: interactions between the migration rate, the size and the number of subpopulations

We assess the relative importance of migration rate, size and number of subpopulations on the genetic load of subdivided populations. Using diffusion approximations, we show that in most cases subdivision has detrimental effects on fitness. Moreover, our results suggest that fitness increases with subpopulation size, so that for the same total population size, genetic load is relatively lower when there are a small number of large subpopulations. Using elasticity analysis, we show that the size of the subpopulations appears to be the parameter that most strongly determines genetic load. interconnecting subpopulations via migration would also be of importance for population fitness when subpopulations are small and gene flow is low. Interestingly, the number of subpopulations has minor influence on genetic load except for the case of both very slightly deleterious mutations and small subpopulations. Elasticities decrease as the magnitude of deleterious effects increases. In other words, population structure does not matter for very deleterious alleles, but strongly affects fitness for slightly deleterious alleles.

Persistence time of loss-of-function mutations at nonessential loci affecting eye color in Drosophila melanogaster

Persistence time of a mutant allele, the expected number of generations before its elimination from the population, can be estimated as the ratio of the number of segregating mutations per individual over the mutation rate per generation. We screened two natural populations of Drosophila melanogaster for mutations causing clear-cut eye phenotypes and detected 25 mutant alleles, falling into 19 complementation groups, in 1164 haploid genomes, which implies 0.021 eye mutations/genome. The de novo haploid mutation rate for the same set of loci was estimated as 2 x 10(-4) in a 10-generation mutation-accumulation experiment. Thus, the average persistence time of all mutations causing clear-cut eye phenotypes is approximately 100 generations (95% confidence interval: 61-219). This estimate shows that the strength of selection against phenotypically drastic alleles of nonessential loci is close to that against recessive lethals. In both cases, deleterious alleles are apparently eliminated by selection against heterozygous individuals, which show no visible phenotypic differences from wild type.

Relaxation of selection with equalization of parental contributions in conservation programs: an experimental test with Drosophila melanogaster

Equalization of parental contributions is one of the most simple and widely recognized methods to maintain genetic diversity in conservation programs, as it halves the rate of increase in inbreeding and genetic drift. It has, however, the negative side effect of implying a reduced intensity of natural selection so that deleterious genes are less efficiently removed from the population with possible negative consequences on the reproductive capacity of the individuals. Theoretical results suggest that the lower fitness resulting from equalization of family sizes relative to that for free contribution schemes is expected to be substantial only for relatively large population sizes and after many generations. We present a long-term experiment with Drosophila melanogaster, comparing the fitness performance of lines maintained with equalization of contributions (EC) and others maintained with no management (NM), allowing for free matings and contributions from parents. Two (five) replicates of size N = 100 (20) individuals of each type of line were maintained for 38 generations. As expected, EC lines retained higher gene diversity and allelic richness for four microsatellite markers and a higher heritability for sternopleural bristle number. Measures of life-history traits, such as egg-to-adult viability, mating success, and global fitness declined with generations, but no significant differences were observed between EC and NM lines. Our results, therefore, provide no evidence to suggest that equalization of family sizes entails a disadvantage on the reproductive capacity of conserved populations in comparison with no management procedures, even after long periods of captivity.

The effect of antagonistic pleiotropy on the estimation of the average coefficient of dominance of deleterious mutations

We investigate the impact of antagonistic pleiotropy on the most widely used methods of estimation of the average coefficient of dominance of deleterious mutations from segregating populations. A proportion of the deleterious mutations affecting a given studied fitness component are assumed to have an advantageous effect on another one, generating overdominance on global fitness. Using diffusion approximations and transition matrix methods, we obtain the distribution of gene frequencies for nonpleiotropic and pleiotropic mutations in populations at the mutation-selection-drift balance. From these distributions we build homozygous and heterozygous chromosomes and assess the behavior of the estimators of dominance. A very small number of deleterious mutations with antagonistic pleiotropy produces substantial increases on the estimate of the average degree of dominance of mutations affecting the fitness component under study. For example, estimates are increased three- to fivefold when 2% of segregating loci are over-dominant for fitness. In contrast, strengthening pleiotropy, where pleiotropic effects are assumed to be also deleterious, has little effect on the estimates of the average degree of dominance, supporting previous results. The antagonistic pleiotropy model considered, applied under mutational parameters described in the literature, produces patterns for the distribution of chromosomal viabilities, levels of genetic variance, and homozygous mutation load generally consistent with those observed empirically for viability in Drosophila melanogaster.

Analysis of the estimators of the average coefficient of dominance of deleterious mutations

We investigate the sources of bias that affect the most commonly used methods of estimation of the average degree of dominance (h) of deleterious mutations, focusing on estimates from segregating populations. The main emphasis is on the effect of the finite size of the populations, but other sources of bias are also considered. Using diffusion approximations to the distribution of gene frequencies in finite populations as well as stochastic simulations, we assess the behavior of the estimators obtained from populations at mutation-selection-drift balance under different mutational scenarios and compare averages of h for newly arisen and segregating mutations. Because of genetic drift, the inferences concerning newly arisen mutations based on the mutation-selection balance theory can have substantial upward bias depending upon the distribution of h. In addition, estimates usually refer to h weighted by the homozygous deleterious effect in different ways, so that inferences are complicated when these two variables are negatively correlated. Due to both sources of bias, the widely used regression of heterozygous on homozygous means underestimates the arithmetic mean of h for segregating mutations, in contrast to their repeatedly assumed equality in the literature. We conclude that none of the estimators from segregating populations provides, under general conditions, a useful tool to ascertain the properties of the degree of dominance, either for segregating or for newly arisen deleterious mutations. Direct estimates of the average h from mutation-accumulation experiments are shown to suffer some bias caused by purging selection but, because they do not require assumptions on the causes maintaining segregating variation, they appear to give a more reliable average dominance for newly arisen mutations.

Quantitative genetic analyses of complex behaviours in Drosophila

Behaviours are exceptionally complex quantitative traits. Sensitivity to environmental variation and genetic background, the presence of sexual dimorphism, and the widespread functional pleiotropy that is inherent in behavioural phenotypes pose daunting challenges for unravelling their underlying genetics. Drosophila melanogaster provides an attractive system for elucidating the unifying principles of the genetic architectures that drive behaviours, as genetically identical individuals can be reared rapidly in controlled environments and extensive publicly accessible genetic resources are available. Recent advances in quantitative genetic and functional genomic approaches now enable the extensive characterization of complex genetic networks that mediate behaviours in this important model organism.