Allocation of resources to male and female functions in hermaphrodites

The question of"how a self-fertile hermaphrodite will distribute the resources that it allocates to reproduction is studied by means of the ESS approach. Different models of the relations between allocation to male function, the male and female fertilities, and the selfing rate, yield different conclusions about how much resource should be allocated to male function. Values below a half are obtained with one model, while another can give values greater than a half. Even with no selfing, values other than a half are usually obtained; with both models studied, the values decrease with increasing selling. If the selfing rate is assumed to be independent of the fraction of resources allocated to male function, it can be shown that the ESS allocation to male function always decreases as selling increases. The types of relations that might be expected in species with different types of breeding biology, and some data on allocation to male function, are reviewed.
The implications for the fitness of male- and female-sterility mutations are discussed. It is argued that the concavity or convexity of the curve relating female fertility to male fertility is not a good guide to when hermaphroditism should exist when there is some selfing. Even with a concave relation, male-sterility mutants can have a higher fitness than hermaphrodites, if there is some selling and inbreeding depression. Also, when the selfing rate depends on allocation to male I unction, an hermaphrodite ESS does not always exist when the function is concave (as it does when there is no selfing), and such an ESS may exist when the relation is convex. The fitness of male- or female-sterility mutants may also depend on the existence of 'fixed costs'. It is shown that these do not ailed the ESS allocation of resources.     

Evolutionary consequences of gender plasticity in genetically dimorphic breeding systems

A functional view of gender helps evolutionary biologists evaluate the mechanisms underlying breeding-system evolution. Evolutionary pathways from hermaphroditism to dioecy include the intermediate breeding systems of gynodioecy and androdioecy. These pathways start with the invasion of unisexual mutants, females or males, respectively, followed by alteration of the hermaphrodites to allocate more to the sexual function that the unisexuals lack. Eventually, hermaphrodites become unisexual and dioecy has evolved. Some species evolving along these pathways stop short of completing this second step, or even revert back from dioecy. We evaluate the hypothesis that gender plasticity is involved in these transitions to and from dioecy. Evidence from studies of subdioecious species that have evolved along the gynodioecy pathway suggests that gender plasticity occurs and stabilizes subdioecy by lowering the cost of producing seed. Factors influencing species evolving toward androdioecy, or reverting to androdioecy from dioecy, appear to be more varied and include reproductive assurance, herbivory and gender plasticity. In general, gender specialization appears to be favored in resource-poor environments regardless of which pathway is taken to dioecy.     

A pooling approach to detect signatures of selective sweeps in genome scans using microsatellites

We have evaluated a pooling approach that can reduce the number of polymerase chain reactions in a screen for selective sweeps by more than an order of magnitude. We show that the complex peak pattern that results from pooling of all samples from a given population is a faithful reflection of the composite pattern of the individual alleles, although with an under‐representation of the larger alleles. Candidate loci for selective sweeps can be identified by visual inspection of the pool patterns. We have also implemented a software tool, which can find suitable microsatellite loci in the vicinity of annotated genes.     

Working at the interface of phylogenetics and population genetics: a biogeographical analysis of Triaenops spp. (Chiroptera: Hipposideridae)

New applications of genetic data to questions of historical biogeography have revolutionized our understanding of how organisms have come to occupy their present distributions. Phylogenetic methods in combination with divergence time estimation can reveal biogeographical centres of origin, differentiate between hypotheses of vicariance and dispersal, and reveal the directionality of dispersal events. Despite their power, however, phylogenetic methods can sometimes yield patterns that are compatible with multiple, equally well-supported biogeographical hypotheses. In such cases, additional approaches must be integrated to differentiate among conflicting dispersal hypotheses. Here, we use a synthetic approach that draws upon the analytical strengths of coalescent and population genetic methods to augment phylogenetic analyses in order to assess the biogeographical history of Madagascar's Triaenops bats (Chiroptera: Hipposideridae). Phylogenetic analyses of mitochondrial DNA sequence data for Malagasy and east African Triaenops reveal a pattern that equally supports two competing hypotheses. While the phylogeny cannot determine whether Africa or Madagascar was the centre of origin for the species investigated, it serves as the essential backbone for the application of coalescent and population genetic methods. From the application of these methods, we conclude that a hypothesis of two independent but unidirectional dispersal events from Africa to Madagascar is best supported by the data.     

Monte Carlo experiments and the defense of diffusion models in molecular population genetics

In the 1960s molecular population geneticists used Monte Carlo experiments to evaluate particular diffusion equation models. In this paper I examine the nature of this comparative evaluation and argue for three claims: first, Monte Carlo experiments are genuine experiments: second, Monte Carlo experiments can provide an important meansfor evaluating the adequacy of highly idealized theoretical models; and, third, the evaluation of the computational adequacy of a diffusion model with Monte Carlo experiments is significantlydifferent from the evaluation of the emperical adequacy of the same diffusion model.     

Seed fitness of hermaphrodites in areas with females and anther smut disease: Silene acaulis and Microbotryum violaceum

Sex-dependent infection rates could change the effective sex ratio of a population. Here, I tested whether females and hermaphrodites of Silene acaulis were equally likely to be infected by Microbotryum violaceum, a fungus that sterilizes the host, and whether sex allocation in hermaphrodites differed between low and high disease plots. Sex ratios of healthy and diseased plants were estimated in five natural plots. Fitness gained through seed production was estimated by measuring seed quantity and quality for each sex morph in eight plots for 2 yr; four plots had 1-5% disease frequency and four plots had 18-25% disease frequency. Sex ratios of healthy and diseased plants did not differ in five plots. The proportion of fitness hermaphrodites gained through ovules varied from 25 to 48%, indicating that this population is near the cosexual end of gynodioecy. Variation in functional gender of hermaphrodites was not explained by sex-dependent infection rates. Spatial heterogeneity in resources and microclimate seems to be important in explaining both disease frequency and variation in seed production by females and hermaphrodites.     

Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows

We present here a new version of the Arlequin program available under three different forms: a Windows graphical version (Winarl35), a console version of Arlequin (arlecore), and a specific console version to compute summary statistics (arlsumstat). The command-line versions run under both Linux and Windows. The main innovations of the new version include enhanced outputs in XML format, the possibility to embed graphics displaying computation results directly into output files, and the implementation of a new method to detect loci under selection from genome scans. Command-line versions are designed to handle large series of files, and arlsumstat can be used to generate summary statistics from simulated data sets within an Approximate Bayesian Computation framework.     

Relative fitness of females and hermaphrodites in a natural gynodioecious population of wild radish,Raphanus sativus L. (Brassicaceae): comparison based on molecular genotyping

In many gynodioecious species, sex determination involves both cytoplasmic male‐sterility (CMS) genes and nuclear genes that restore male function. Differences in fitness among genotypes affect the dynamics of those genes, and thus that of gynodioecy. We used a molecular marker to discriminate between hermaphrodites with and without a CMS gene in gynodioecious Raphanus sativus. We compared fitness through female function among the three genotypes: females, hermaphrodites with the CMS gene and those without it. Although there was no significant difference among the genotypes in seed size, hermaphrodites without the CMS gene produced significantly more seeds, and seeds with a higher germination rate than the other genotypes, suggesting no fitness advantage for females and no benefit to bearing the CMS gene. Despite the lack of fitness advantage for females in the parameter values we estimated, a theoretical model of gynodioecy shows it can be maintained if restorer genes impose a cost paid in pollen production. In addition, we found that females invest more resources into female reproduction than hermaphrodites when they become larger. If environmental conditions enable females to grow larger this would facilitate the dynamics of CMS genes.     

Estimation of mating system parameters in an evolving gynodioecous population of cultivated sunflower (Helianthus annuus L.)

Cultivated plants have been molded by human-induced selection, including manipulations of the mating system in the twentieth century. How these manipulations have affected realized parameters of the mating system in freely evolving cultivated populations is of interest for optimizing the management of breeding populations, predicting the fate of escaped populations and providing material for experimental evolution studies. To produce modern varieties of sunflower (Helianthus annuus L.), self-incompatibility has been broken, recurrent generations of selfing have been performed and male sterility has been introduced. Populations deriving from hybrid-F1 varieties are gynodioecious because of the segregation of a nuclear restorer of male fertility. Using both phenotypic and genotypic data at 11 microsatellite loci, we analyzed the consanguinity status of plants of the first three generations of such a population and estimated parameters related to the mating system. We showed that the resource reallocation to seed in male-sterile individuals was not significant, that inbreeding depression on seed production averaged 15-20% and that cultivated sunflower had acquired a mixed-mating system, with ～50% of selfing among the hermaphrodites. According to theoretical models, the female advantage and the inbreeding depression at the seed production stage were too low to allow the persistence of male sterility. We discuss our methods of parameter estimation and the potential of such study system in evolutionary biology.     

Surprising Differences in the Variability of Y Chromosomes in African and Cosmopolitan Populations of Drosophila melanogaster

The nonrecombining Drosophila melanogaster Y chromosome is heterochromatic and has few genes. Despite these limitations, there remains ample opportunity for natural selection to act on the genes that are vital for male fertility and on Y factors that modulate gene expression elsewhere in the genome. Y chromosomes of many organisms have low levels of nucleotide variability, but a formal survey of D. melanogaster Y chromosome variation had yet to be performed. Here we surveyed Y-linked variation in six populations of D. melanogaster spread across the globe. We find surprisingly low levels of variability in African relative to Cosmopolitan (i.e., non-African) populations. While the low levels of Cosmopolitan Y chromosome polymorphism can be explained by the demographic histories of these populations, the staggeringly low polymorphism of African Y chromosomes cannot be explained by demographic history. An explanation that is entirely consistent with the data is that the Y chromosomes of Zimbabwe and Uganda populations have experienced recent selective sweeps. Interestingly, the Zimbabwe and Uganda Y chromosomes differ: in Zimbabwe, a European Y chromosome appears to have swept through the population.     

Molecular population genetics and the search for adaptive evolution in plants

The first papers on plant molecular population genetics were published approximately 10 years ago. Since that time, well over 50 additional studies of plant nucleotide polymorphism have been published, and many of these studies focused on detecting the signature of balancing or positive selection at a locus. In this review, we discuss some of the theoretical and statistical issues surrounding the detection of selection, with focus on plant populations, and we also summarize the empirical plant molecular population genetics literature. At face value, the literature suggests that a history of balancing or positive selection in plant genes is rampant. In two well-studied taxa (maize and Arabidopsis) over 20% of studied genes have been interpreted as containing the signature of selection. We argue that this is probably an overstatement of the prevalence of natural selection in plant genomes, for two reasons. First, demographic effects are difficult to incorporate and have generally not been well integrated into the plant population genetics literature. Second, the genes studied to date are not a random sample, so selected genes may be overrepresented. The next generation of studies in plant molecular population genetics requires additional sampling of local populations, explicit comparisons among loci, and improved theoretical methods to control for demography. Eventually, candidate loci should be confirmed by explicit consideration of phenotypic effects.     

Evolutionary game dynamics in a finite asymmetric two-deme population and emergence of cooperation

We consider evolutionary game dynamics in a finite population subdivided into two demes with both unequal deme sizes and different migration rates. Assuming viability differences in the population according to a linear game within each deme as a result of pairwise interactions, we specify conditions for weak selection favoring a mutant strategy to go to fixation, under the structured-coalescent assumptions, and their connections with evolutionary stability concepts. In the framework of the Iterated Prisoner's Dilemma with strategy ‘tit-for-tat’ as mutant strategy and ‘always defect’ as resident strategy, we deduce a condition under which the emergence of cooperation is favored by selection, when the game matrix is the same in both demes. We show how this condition extends the one-third law for a panmictic population and when an asymmetry in the spatial structure of a two-deme population facilitates the emergence of the cooperative tit-for-tat strategy in comparison with both its symmetric and panmictic population structure counterparts. We find that the condition is less stringent in the asymmetric scenario versus the symmetric scenario if both the fraction of the population in the deme where the mutant was initially introduced, and the expected proportion of migrant offspring in this deme among all migrant offspring after population regulation, are smaller than, or equal to, , provided they are not too small. On the other hand, the condition is less stringent than the one-third law, which holds in the panmictic case, if the latter proportion remains not too close to 1.     

Combining models of protein translation and population genetics to predict protein production rates from codon usage patterns

Genes are often biased in their codon usage. The degree of bias displayed often changes with expression level and intragenic position. Numerous indices, such as the codon adaptation index, have been developed to measure this bias. Although the expression level of a gene and index values are correlated, the heuristic nature of these metrics limits their ability to explain this relationship. As an alternative approach, this study integrates mechanistic models of cellular and population processes in a nested manner to develop a stochastic evolutionary model of a protein's production rate (SEMPPR). SEMPPR assumes that the evolution of codon bias is driven by selection to reduce the cost of nonsense errors and that this selection is counteracted by mutation and drift. Through the application of Bayes' theorem, SEMPPR generates a posterior probability distribution for the protein production rate of a given gene. Conceptually, SEMPPR's predictions are based on the degree of adaptation to reduce the cost of nonsense errors observed in the codon usage pattern of the gene. As an illustration, SEMPPR was parameterized using the Saccharomyces cerevisiae genome and its predictions tested using available empirical data. The results indicate that SEMPPR's predictions are as reliable index based ones. In addition, SEMPPR's output is more easily interpreted and its predictions could be improved through refinements of the models upon which it is built.     

Molecular evolution and population genetics of duplicated accessory gland protein genes in Drosophila

To investigate the potential importance of gene duplication in D. melanogaster accessory gland protein (Acp) gene evolution we carried out a computational analysis comparing annotated D. melanogaster Acp genes to the entire D. melanogaster genome. We found that two known Acp genes are actually members of small multigene families. Polymorphism and divergence data from these duplicated genes suggest that in at least four cases, protein divergence between D. melanogaster and D. simulans is a result of directional selection. One putative Acp revealed by our computational analysis shows evidence of a recent selective sweep in a non-African population (but not in an African population). These data support the idea that selection on reproduction-related genes may drive divergence of populations within species, and strengthen the conclusion that Acps may often be under directional selection in Drosophila.     

Recurrent gene amplification and soft selective sweeps during evolution of multidrug resistance in malaria parasites

When selection is strong and beneficial alleles have a single origin, local reductions in genetic diversity are expected. However, when beneficial alleles have multiple origins or were segregating in the population prior to a change in selection regime, the impact on genetic diversity may be less clear. We describe an example of such a "soft" selective sweep in the malaria parasite Plasmodium falciparum that involves adaptive genome rearrangements. Amplification in copy number of genome regions containing the pfmdr1 gene on chromosome 5 confer resistance to mefloquine and spread rapidly in the 1990s. Using flanking microsatellite data and real-time polymerase chain reaction determination of copy number, we show that 5-15 independent amplification events have occurred in parasites on the Thailand/Burma border. The amplified genome regions (amplicons) range in size from 14.7 to 49 kb and contain 2-11 genes, with 2-4 copies arranged in tandem. To examine the impact of drug selection on flanking variation, we genotyped 48 microsatellites on chromosome 5 in 326 parasites from a single Thai location. Diversity was reduced in a 170- to 250-kb (10-15 cM) region of chromosomes containing multiple copies of pfmdr1, consistent with hitchhiking resulting from the rapid recent spread of selected chromosomes. However, diversity immediately flanking pfmdr1 is reduced by only 42% on chromosomes bearing multiple amplicons relative to chromosomes carrying a single copy. We highlight 2 features of these results: 1) All amplicon break points occur in monomeric A/T tracts (9-45 bp). Given the abundance of these tracts in P. falciparum, we expect that duplications will occur frequently at multiple genomic locations and have been underestimated as drivers of phenotypic evolution in this pathogen. 2) The signature left by the spread of amplified genome segments is broad, but results in only limited reduction in diversity. If such "soft" sweeps are common in nature, statistical methods based on diversity reduction may be inefficient at detecting evidence for selection in genome-wide marker screens. This may be particularly likely when mutation rate is high, as appears to be the case for gene duplications, and in pathogen populations where effective population sizes are typically very large.     

Existence and nonexistence of steady-state solutions for a selection-migration model in population genetics

We discuss a selection-migration model in population genetics, involving two alleles A₁ and A₂ such that A₁ is at an advantage over A₂ in certain subregions and at a disadvantage in others. It is shown that if A₁ is at an overall disadvantage to A₂ and the rate of gene flow is sufficiently large than A₁ must die out; on the other hand, if the two alleles are in some sense equally advantaged overall, then A₁ and A₂ can coexist no matter how great the rate of gene flow.     

The evolution from females to hermaphrodites results in a sexual conflict over mating in androdioecious nematode worms and clam shrimp

The nematode worm Caenorhabditis elegans and the clam shrimp Eulimnadia texana are two well‐studied androdioecious species consisting mostly of self‐fertilizing hermaphrodites and few males. To understand how androdioecy can evolve, a simple two‐step mathematical model of the evolutionary pathway from a male–female species to a selfing‐hermaphrodite species is constructed. First, the frequency of mutant females capable of facultative self‐fertilization increases if the benefits of reproductive assurance exceed the cost. Second, hermaphrodites become obligate self‐fertilizers if the fitness of selfed offspring exceeds one‐half the fitness of outcrossed offspring. Genetic considerations specific to C. elegans and E. texana show that males may endure as descendants of the ancestral male–female species. These models combined with an extensive literature review suggest a sexual conflict over mating in these androdioecious species: selection favours hermaphrodites that self and males that outcross. The strength of selection on hermaphrodites and males differs, however. Males that fail to outcross suffer a genetic death. Hermaphrodites may never encounter a rare male, and those that do and outcross only bear less fecund offspring. This asymmetric sexual conflict results in an evolutionary stand‐off: rare, but persistent males occasionally fertilize common, but reluctant hermaphrodites. A consequence of this stand‐off may be an increase in the longevity of the androdioecious mating system.     

Molecular population genetics of herbivore-induced protease inhibitor genes in European aspen (Populus tremula L., Salicaceae)

Plants defend themselves against the attack of natural enemies by using an array of both constitutively expressed and induced defenses. Long-lived woody perennials are overrepresented among plant species that show strong induced defense responses, whereas annual plants and crop species are underrepresented. However, most studies of plant defense genes have been performed on annual or short-lived perennial weeds or crop species. Here I use molecular population genetic methods to survey six wound-inducible protease inhibitors (PIs) in a long-lived woody, perennial plant species, the European aspen (Populus tremula), to evaluate the likelihood of either recurrent selective sweeps or balancing selection maintaining amino acid polymorphisms in these genes. The results show that none of the six PI genes have reduced diversities at synonymous sites, as would be expected in the presence of recurrent selective sweeps. However, several genes show some evidence of nonneutral evolution such as enhanced linkage disequilibrium and a large number of high-frequency-derived mutations. A group of at least four Kunitz trypsin inhibitor genes appear to have experienced elevated levels of nonsynonymous substitutions, indicating allelic turnover on an evolutionary timescale. One gene, TI1, has enhanced levels of intraspecific polymorphism at nonsynonymous sites and also has an unusual haplotype structure characterized by two divergent haplotypes occurring at roughly equal frequencies in the sample. One haplotype has very low levels of intraallelic nucleotide diversity, whereas the other haplotype has levels of diversity comparable to other genes in P. tremula. Patterns of sequence diversity at TI1 do not fit a simple model of either balancing selection or recurrent selective sweeps. This suggests that selection at TI1 is more complex, possibly involving allelic cycling.     

Strategic analysis in evolutionary genetics and the theory of games

This paper is written in memory of John Maynard Smith. In a brief survey it discusses essential aspects of how game theory in biology relates to its counterpart in economics, the major transition in game theory initiated by Maynard Smith, the discrepancies between genetic and phenotypic models in evolutionary biology, and a balanced way of reconciling these models. In addition, the paper discusses modern problems in understanding games at the genetic level using the examples of conflict between endosymbionts and their hosts, and the molecular interactions between parasites and the mammalian immune system.     

How to track and assess genotyping errors in population genetics studies

Genotyping errors occur when the genotype determined after molecular analysis does not correspond to the real genotype of the individual under consideration. Virtually every genetic data set includes some erroneous genotypes, but genotyping errors remain a taboo subject in population genetics, even though they might greatly bias the final conclusions, especially for studies based on individual identification. Here, we consider four case studies representing a large variety of population genetics investigations differing in their sampling strategies (noninvasive or traditional), in the type of organism studied (plant or animal) and the molecular markers used [microsatellites or amplified fragment length polymorphisms (AFLPs)]. In these data sets, the estimated genotyping error rate ranges from 0.8% for microsatellite loci from bear tissues to 2.6% for AFLP loci from dwarf birch leaves. Main sources of errors were allelic dropouts for microsatellites and differences in peak intensities for AFLPs, but in both cases human factors were non-negligible error generators. Therefore, tracking genotyping errors and identifying their causes are necessary to clean up the data sets and validate the final results according to the precision required. In addition, we propose the outline of a protocol designed to limit and quantify genotyping errors at each step of the genotyping process. In particular, we recommend (i) several efficient precautions to prevent contaminations and technical artefacts; (ii) systematic use of blind samples and automation; (iii) experience and rigor for laboratory work and scoring; and (iv) systematic reporting of the error rate in population genetics studies.