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.     

Plasticity in sex allocation in the plant Mercurialis annua is greater for hermaphrodites sampled from dimorphic than from monomorphic populations

Plants are notoriously variable in gender, ranging in sex allocation from purely male through hermaphrodite to purely female. This variation can have both a genetic and an adaptive plastic component. In gynodioecious species, where females co‐occur with hermaphrodites, hermaphrodites tend to shift their allocation towards greater maleness when growing under low‐resource conditions, either as a result of hermaphrodites shifting away from an expensive female function, or because of enhanced siring advantages in the presence of females. Similarly, in the androdioecious plant Mercurialis annua, where hermaphrodites co‐exist with males, hermaphrodites also tend to enhance their relative male allocation under low‐resource conditions. Here, we ask whether this response differs between hermaphrodites that have been evolving in the presence of males, in a situation analogous to that supposed for gynodioecious populations, vs. those that have been evolving in their absence. We grew hermaphrodites of M. annua from populations in which males were either present or absent under different levels of nutrient availability and compared their reaction norms. We found that, overall, hermaphrodites from populations with males tended to be more female than those from populations lacking males. Importantly, hermaphrodites' investment in pollen and seed production was more plastic when they came from populations with males than without them, reducing their pollen production at low resource availability and increasing their seed production at high resource availability. These results are consistent with the hypothesis that plasticity in sex allocation is enhanced in hermaphrodites that have likely been exposed to variation in mating opportunities due to fluctuations in the frequency of co‐occurring males.     

Sweeps in time: leveraging the joint distribution of branch lengths

Current methods of identifying positively selected regions in the genome are limited in two key ways: the underlying models cannot account for the timing of adaptive events and the comparison between models of selective sweeps and sequence data is generally made via simple summaries of genetic diversity. Here, we develop a tractable method of describing the effect of positive selection on the genealogical histories in the surrounding genome, explicitly modeling both the timing and context of an adaptive event. In addition, our framework allows us to go beyond analyzing polymorphism data via the site frequency spectrum or summaries thereof and instead leverage information contained in patterns of linked variants. Tests on both simulations and a human data example, as well as a comparison to SweepFinder2, show that even with very small sample sizes, our analytic framework has higher power to identify old selective sweeps and to correctly infer both the time and strength of selection. Finally, we derived the marginal distribution of genealogical branch lengths at a locus affected by selection acting at a linked site. This provides a much-needed link between our analytic understanding of the effects of sweeps on sequence variation and recent advances in simulation and heuristic inference procedures that allow researchers to examine the sequence of genealogical histories along the genome.     

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.     

Tree‐sequence recording in SLiM opens new horizons for forward‐time simulation of whole genomes

There is an increasing demand for evolutionary models to incorporate relatively realistic dynamics, ranging from selection at many genomic sites to complex demography, population structure, and ecological interactions. Such models can generally be implemented as individual‐based forward simulations, but the large computational overhead of these models often makes simulation of whole chromosome sequences in large populations infeasible. This situation presents an important obstacle to the field that requires conceptual advances to overcome. The recently developed tree‐sequence recording method (Kelleher, Thornton, Ashander, & Ralph, 2018), which stores the genealogical history of all genomes in the simulated population, could provide such an advance. This method has several benefits: (1) it allows neutral mutations to be omitted entirely from forward‐time simulations and added later, thereby dramatically improving computational efficiency; (2) it allows neutral burn‐in to be constructed extremely efficiently after the fact, using “recapitation”; (3) it allows direct examination and analysis of the genealogical trees along the genome; and (4) it provides a compact representation of a population's genealogy that can be analysed in Python using the msprime package. We have implemented the tree‐sequence recording method in SLiM 3 (a free, open‐source evolutionary simulation software package) and extended it to allow the recording of non‐neutral mutations, greatly broadening the utility of this method. To demonstrate the versatility and performance of this approach, we showcase several practical applications that would have been beyond the reach of previously existing methods, opening up new horizons for the modelling and exploration of evolutionary processes.     

Detecting Selection in Multiple Populations by Modeling Ancestral Admixture Components

One of the most powerful and commonly used approaches for detecting local adaptation in the genome is the identification of extreme allele frequency differences between populations. In this article, we present a new maximum likelihood method for finding regions under positive selection. It is based on a Gaussian approximation to allele frequency changes and it incorporates admixture between populations. The method can analyze multiple populations simultaneously and retains power to detect selection signatures specific to ancestry components that are not representative of any extant populations. Using simulated data, we compare our method to related approaches, and show that it is orders of magnitude faster than the state-of-the-art, while retaining similar or higher power for most simulation scenarios. We also apply it to human genomic data and identify loci with extreme genetic differentiation between major geographic groups. Many of the genes identified are previously known selected loci relating to hair pigmentation and morphology, skin, and eye pigmentation. We also identify new candidate regions, including various selected loci in the Native American component of admixed Mexican-Americans. These involve diverse biological functions, such as immunity, fat distribution, food intake, vision, and hair development.     

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.     

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.     

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.     

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.     

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.     

The population genetics of adaptation on correlated fitness landscapes: the block model

Several recent theoretical studies of the genetics of adaptation have focused on the mutational landscape model, which considers evolution on rugged fitness landscapes (i.e., ones having many local optima). Adaptation in this model is characterized by several simple results. Here I ask whether these results also hold on correlated fitness landscapes, which are smoother than those considered in the mutational landscape model. In particular, I study the genetics of adaptation in the block model, a tunably rugged model of fitness landscapes. Considering the scenario in which adaptation begins from a high fitness wild-type DNA sequence, I use extreme value theory and computer simulations to study both single adaptive steps and entire adaptive walks. I show that all previous results characterizing single steps in adaptation in the mutational landscape model hold at least approximately on correlated landscapes in the block model; many entire-walk results, however, do not.     

Whole‐genome resequencing uncovers molecular signatures of natural and sexual selection in wild bighorn sheep

The identification of genes influencing fitness is central to our understanding of the genetic basis of adaptation and how it shapes phenotypic variation in wild populations. Here, we used whole‐genome resequencing of wild Rocky Mountain bighorn sheep (Ovis canadensis) to >50‐fold coverage to identify 2.8 million single nucleotide polymorphisms (SNPs) and genomic regions bearing signatures of directional selection (i.e. selective sweeps). A comparison of SNP diversity between the X chromosome and the autosomes indicated that bighorn males had a dramatically reduced long‐term effective population size compared to females. This probably reflects a long history of intense sexual selection mediated by male–male competition for mates. Selective sweep scans based on heterozygosity and nucleotide diversity revealed evidence for a selective sweep shared across multiple populations at RXFP2, a gene that strongly affects horn size in domestic ungulates. The massive horns carried by bighorn rams appear to have evolved in part via strong positive selection at RXFP2. We identified evidence for selection within individual populations at genes affecting early body growth and cellular response to hypoxia; however, these must be interpreted more cautiously as genetic drift is strong within local populations and may have caused false positives. These results represent a rare example of strong genomic signatures of selection identified at genes with known function in wild populations of a nonmodel species. Our results also showcase the value of reference genome assemblies from agricultural or model species for studies of the genomic basis of adaptation in closely related wild taxa.     

Comparative and population genomics approaches reveal the basis of adaptation to deserts in a small rodent

Organisms that live in deserts offer the opportunity to investigate how species adapt to environmental conditions that are lethal to most plants and animals. In the hot deserts of North America, high temperatures and lack of water are conspicuous challenges for organisms living there. The cactus mouse (Peromyscus eremicus) displays several adaptations to these conditions, including low metabolic rate, heat tolerance, and the ability to maintain homeostasis under extreme dehydration. To investigate the genomic basis of desert adaptation in cactus mice, we built a chromosome‐level genome assembly and resequenced 26 additional cactus mouse genomes from two locations in southern California (USA). Using these data, we integrated comparative, population, and functional genomic approaches. We identified 16 gene families exhibiting significant contractions or expansions in the cactus mouse compared to 17 other Myodontine rodent genomes, and found 232 sites across the genome associated with selective sweeps. Functional annotations of candidate gene families and selective sweeps revealed a pervasive signature of selection at genes involved in the synthesis and degradation of proteins, consistent with the evolution of cellular mechanisms to cope with protein denaturation caused by thermal and hyperosmotic stress. Other strong candidate genes included receptors for bitter taste, suggesting a dietary shift towards chemically defended desert plants and insects, and a growth factor involved in lipid metabolism, potentially involved in prevention of dehydration. Understanding how species adapted to deserts will provide an important foundation for predicting future evolutionary responses to increasing temperatures, droughts and desertification in the cactus mouse and other species.     

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.     

Mechanisms of population differentiation in marbled murrelets: historical versus contemporary processes

Abstract.— Mechanisms of population differentiation in highly vagile species such as seabirds are poorly understood. Previous studies of marbled murrelets ( Brachyramphus marmoratus ; Charadriiformes: Alcidae) found significant population genetic structure, but could not determine whether this structure is due to historical vicariance (e.g., due to Pleistocene glaciers), isolation by distance, drift or selection in peripheral populations, or nesting habitat selection. To discriminate among these possibilities, we analyzed sequence variation in nine nuclear introns from 120 marbled murrelets sampled from British Columbia to the western Aleutian Islands. Mismatch distributions indicated that murrelets underwent at least one population expansion during the Pleistocene and probably are not in genetic equilibrium. Maximum-likelihood analysis of allele frequencies suggested that murrelets from "mainland" sites (from the Alaskan Peninsula east) are genetically different from those in the Aleutians and that these two lineages diverged prior to the last glaciation. Analyses of molecular variance, as well as estimates of gene flow derived using coalescent theory, indicate that population genetic structure is best explained by peripheral isolation of murrelets in the Aleutian Islands, rather than by selection associated with different nesting habitats. No isolation-by-distance effects could be detected. Our results are consistent with a rapid expansion of murrelets from a single refugium during the early–mid Pleistocene, subsequent isolation and divergence in two or more refugia during the final Pleistocene glacial advance, and secondary contact following retreat of the ice sheets. Population genetic structure now appears to be maintained by distance effects combined with small populations and a highly fragmented habitat in the Aleutian Islands.     

The population genetics of phenotypic deterioration in experimental populations of Bacillus subtilis

Although many examples of trait loss exist in nature, the underlying population genetic mechanism responsible for the loss is usually unknown. Selective or neutral processes can result in the deterioration of a trait, and often one of these is inferred based on indirect evidence. Furthermore, selective pressures that are unique to particular environments and the effect these might have on the population genetic cause of trait loss are not well understood. Here we describe an experimental evolution system where two different environments were used for addressing the population genetic cause of trait loss throughout evolutionary time. We found that growth in minimal medium (i.e., prototrophy) was lost in all populations regardless of the experimental environment and that the pattern of trait loss in one environment was due to selection, whereas in the other environment the cause remains inconclusive.