The Evolutionary Interplay between Adaptation and Self-Fertilization

Characteristics of duplicated gene expression and DNA methylation regulation in different tissues of allopolyploid Brassica napus

Plant polyploidization increases the complexity of epigenomes and transcriptional regulation, resulting in genome evolution and enhanced adaptability. However, few studies have been conducted on the relationship between gene expression and epigenetic modification in different plant tissues after allopolyploidization. In this study, we studied gene expression and DNA methylation modification patterns in four tissues (stems, leaves, flowers and siliques) of Brassica napusand its diploid progenitors. On this basis, the alternative splicing patterns and cis-trans regulation patterns of four tissues in B. napus and its diploid progenitors were also analyzed. It can be seen that the number of alternative splicing occurs in the B. napus is higher than that in the diploid progenitors, and the IR type increases the most during allopolyploidy. In addition, we studied the fate changes of duplicated genes after allopolyploidization in B. napus. We found that the fate of most duplicated genes is conserved, but the number of neofunctionalization and specialization is also large. The genetic fate of B. napus was classified according to five replication types (WGD, PD, DSD, TD, TRD). This study also analyzed generational transmission analysis of expression and DNA methylation patterns. Our study provides a reference for the fate differentiation of duplicated genes during allopolyploidization.

Weak response to selection on stigma-anther distance in a primarily selfing population of yellow monkeyflower

Stebbins hypothesized that selfing lineages are evolutionary dead ends because they lack adaptive potential. While selfing populations often possess limited nucleotide variability compared with closely related outcrossers, reductions in the genetic variability of quantitative characters remain unclear, especially for key traits determining selfing rates. Yellow monkeyflower (Mimulus guttatus) populations generally outcross and maintain extensive quantitative genetic variation in floral traits. Here, we study the Joy Road population (Bodega Bay, CA, USA) of M. guttatus, where individuals exhibit stigma-anther distances (SAD) typical of primarily selfing monkeyflowers. We show that this population is closely related to nearby conspecifics on the Pacific Coast with a modest 33% reduction in genome-wide variation compared with a more highly outcrossing population. A five-generation artificial selection experiment challenged the hypothesis that the Joy Road population harbours comparatively low evolutionary potential in stigma-anther distance, a critical determinant of selfing rate in Mimulus. Artificial selection generated a weak phenotypic response, with low realized heritabilities (0.020-0.028) falling 84% below those measured for floral characters in more highly outcrossing M. guttatus. These results demonstrate substantial declines in evolutionary potential with a transition toward selfing. Whether these findings explain infrequent reversals to outcrossing or general limits on adaptation in selfers requires further investigation.

Complex parental effects impact variation in larval thermal tolerance in a vertically transmitting coral

Coral populations must be able to adapt to changing environmental conditions for coral reefs to persist under climate change. The adaptive potential of these organisms is difficult to forecast due to complex interactions between the host animal, dinoflagellate symbionts and the environment. Here we created 26 larval families from six Montipora capitata colonies from a single reef, showing significant, heritable variation in thermal tolerance. Our results indicate that 9.1% of larvae are expected to exhibit four times the thermal tolerance of the general population. Differences in larval thermotolerance were driven mainly by maternal contributions, but we found no evidence that these effects were driven by symbiont identity despite vertical transmission from the dam. We also document no evidence of reproductive incompatibility attributable to symbiont identity. These data demonstrate significant genetic variation within this population which provides the raw material upon which natural selection can act.

Genetic Causes and Genomic Consequences of Breakdown of Distyly in Linum trigynum

Distyly is an iconic floral polymorphism governed by a supergene, which promotes efficient pollen transfer and outcrossing through reciprocal differences in the position of sexual organs in flowers, often coupled with heteromorphic self-incompatibility. Distyly has evolved convergently in multiple flowering plant lineages, but has also broken down repeatedly, often resulting in homostylous, self-compatible populations with elevated rates of self-fertilization. Here, we aimed to study the genetic causes and genomic consequences of the shift to homostyly in Linum trigynum, which is closely related to distylous Linum tenue. Building on a high-quality genome assembly, we show that L. trigynum harbors a genomic region homologous to the dominant haplotype of the distyly supergene conferring long stamens and short styles in L. tenue, suggesting that loss of distyly first occurred in a short-styled individual. In contrast to homostylous Primula and Fagopyrum, L. trigynum harbors no fixed loss-of-function mutations in coding sequences of S-linked distyly candidate genes. Instead, floral gene expression analyses and controlled crosses suggest that mutations downregulating the S-linked LtWDR-44 candidate gene for male self-incompatibility and/or anther height could underlie homostyly and self-compatibility in L. trigynum. Population genomic analyses of 224 whole-genome sequences further demonstrate that L. trigynum is highly self-fertilizing, exhibits significantly lower genetic diversity genome-wide, and is experiencing relaxed purifying selection and less frequent positive selection on nonsynonymous mutations relative to L. tenue. Our analyses shed light on the loss of distyly in L. trigynum, and advance our understanding of a common evolutionary transition in flowering plants.

Genomic Insights into High-Altitude Adaptation: A Comparative Analysis of Roscoea alpina and R. purpurea in the Himalayas

Environmental stress at high altitudes drives the development of distinct adaptive mechanisms in plants. However, studies exploring the genetic adaptive mechanisms of high-altitude plant species are scarce. In the present study, we explored the high-altitude adaptive mechanisms of plants in the Himalayas through whole-genome resequencing. We studied two widespread members of the Himalayan endemic alpine genus Roscoea (Zingiberaceae): R. alpina (a selfing species) and R. purpurea (an outcrossing species). These species are distributed widely in the Himalayas with distinct non-overlapping altitude distributions; R. alpina is distributed at higher elevations, and R. purpurea occurs at lower elevations. Compared to R. purpurea, R. alpina exhibited higher levels of linkage disequilibrium, Tajima's D, and inbreeding coefficient, as well as lower recombination rates and genetic diversity. Approximately 96.3% of the genes in the reference genome underwent significant genetic divergence (FST ≥ 0.25). We reported 58 completely divergent genes (FST = 1), of which only 17 genes were annotated with specific functions. The functions of these genes were primarily related to adapting to the specific characteristics of high-altitude environments. Our findings provide novel insights into how evolutionary innovations promote the adaptation of mountain alpine species to high altitudes and harsh habitats.

Single-cell eQTL mapping in yeast reveals a tradeoff between growth and reproduction

Expression quantitative trait loci (eQTLs) provide a key bridge between noncoding DNA sequence variants and organismal traits. The effects of eQTLs can differ among tissues, cell types, and cellular states, but these differences are obscured by gene expression measurements in bulk populations. We developed a one-pot approach to map eQTLs in Saccharomyces cerevisiae by single-cell RNA sequencing (scRNA-seq) and applied it to over 100,000 single cells from three crosses. We used scRNA-seq data to genotype each cell, measure gene expression, and classify the cells by cell-cycle stage. We mapped thousands of local and distant eQTLs and identified interactions between eQTL effects and cell-cycle stages. We took advantage of single-cell expression information to identify hundreds of genes with allele-specific effects on expression noise. We used cell-cycle stage classification to map 20 loci that influence cell-cycle progression. One of these loci influenced the expression of genes involved in the mating response. We showed that the effects of this locus arise from a common variant (W82R) in the gene GPA1, which encodes a signaling protein that negatively regulates the mating pathway. The 82R allele increases mating efficiency at the cost of slower cell-cycle progression and is associated with a higher rate of outcrossing in nature. Our results provide a more granular picture of the effects of genetic variants on gene expression and downstream traits.

Can heterosis and inbreeding depression explain the maintenance of outcrossing in a cleistogamous perennial?

PREMISE

What maintains mixed mating is an evolutionary enigma. Cleistogamy-the production of both potentially outcrossing chasmogamous and obligately selfing cleistogamous flowers on the same individual plant-is an excellent system to study the costs of selfing. Inbreeding depression can prevent the evolution of greater selfing within populations, and heterosis in crosses between populations may further tip the balance in favor of outcrossing. Few empirical estimates of inbreeding depression and heterosis in the same system exist for cleistogamous species.

METHODS

We investigate the potential costs of selfing by quantifying inbreeding depression and heterosis in three populations of the cleistogamous perennial Ruellia humilis Nutt (Acanthaceae). We performed three types of hand-pollinations-self, outcross-within, and outcross-between populations-and measured seed number, germination, total flower production, and estimated cumulative fitness for the resulting progeny in a greenhouse experiment.

RESULTS

We found moderate inbreeding depression for cumulative fitness (<30%) in two populations, but outbreeding depression for crosses within a third population (-26%). For between-population crosses, there was weak to modest heterosis (11-47%) in two of the population combinations, but modest to strong outbreeding depression (-21 to -71%) in the other four combinations.

CONCLUSIONS

Neither inbreeding depression nor heterosis was of sufficient magnitude to explain the continued production of chasmogamous flowers given the relative energetic advantage of cleistogamous flowers previously estimated for these populations. Outbreeding depression either within or between populations makes the maintenance of chasmogamous flowers even harder to explain. More information is needed on the genetic basis of cleistogamy to resolve this conundrum.

Influence of the mutation load on the genomic composition of hybrids between outcrossing and self-fertilizing species

Hybridization is a natural process whereby two diverging evolutionary lineages reproduce and create offspring of mixed ancestry. Differences in mating systems (e.g., self-fertilization and outcrossing) are expected to affect the direction and extent of hybridization and introgression in hybrid zones. Among other factors, selfers and outcrossers are expected to differ in their mutation loads. This has been studied both theoretically and empirically; however, conflicting predictions have been made on the effects mutation loads of parental species with different mating systems can have on the genomic composition of hybrids. Here, we develop a multi-locus, selective model to study how the different mutation load built up in selfers and outcrossers as a result of selective interference and homozygosity impact the long-term genetic composition of hybrid populations. Notably, our results emphasize that genes from the parental population with lesser mutation load get rapidly overrepresented in hybrid genomes, regardless of the hybrids own mating system. When recombination tends to be more important than mutation, outcrossers' genomes tend to be of higher quality and prevail. When recombination rates are low, however, selfers' genomes may reach higher quality than outcrossers' genomes and prevail in the hybrids. Taken together, these results provide concrete insights into one of the multiple factors influencing hybrid genome ancestry and introgression patterns in hybrid zones containing species with different mating systems.

Genomic diversity landscapes in outcrossing and selfing Caenorhabditis nematodes

Caenorhabditis nematodes form an excellent model for studying how the mode of reproduction affects genetic diversity, as some species reproduce via outcrossing whereas others can self-fertilize. Currently, chromosome-level patterns of diversity and recombination are only available for self-reproducing Caenorhabditis, making the generality of genomic patterns across the genus unclear given the profound potential influence of reproductive mode. Here we present a whole-genome diversity landscape, coupled with a new genetic map, for the outcrossing nematode C. remanei. We demonstrate that the genomic distribution of recombination in C. remanei, like the model nematode C. elegans, shows high recombination rates on chromosome arms and low rates toward the central regions. Patterns of genetic variation across the genome are also similar between these species, but differ dramatically in scale, being tenfold greater for C. remanei. Historical reconstructions of variation in effective population size over the past million generations echo this difference in polymorphism. Evolutionary simulations demonstrate how selection, recombination, mutation, and selfing shape variation along the genome, and that multiple drivers can produce patterns similar to those observed in natural populations. The results illustrate how genome organization and selection play a crucial role in shaping the genomic pattern of diversity whereas demographic processes scale the level of diversity across the genome as a whole.

Strong Positive Selection Biases Identity-By-Descent-Based Inferences of Recent Demography and Population Structure in Plasmodium falciparum

Malaria genomic surveillance often estimates parasite genetic relatedness using metrics such as Identity-By-Decent (IBD). Yet, strong positive selection stemming from antimalarial drug resistance or other interventions may bias IBD-based estimates. In this study, we utilized simulations, a true IBD inference algorithm, and empirical datasets from different malaria transmission settings to investigate the extent of such bias and explore potential correction strategies. We analyzed whole genome sequence data generated from 640 new and 4,026 publicly available Plasmodium falciparum clinical isolates. Our findings demonstrated that positive selection distorts IBD distributions, leading to underestimated effective population size and blurred population structure. Additionally, we discovered that the removal of IBD peak regions partially restored the accuracy of IBD-based inferences, with this effect contingent on the population's background genetic relatedness. Consequently, we advocate for selection correction for parasite populations undergoing strong, recent positive selection, particularly in high malaria transmission settings.

Genetic diversity and selection signatures in a gene bank panel of maize inbred lines from Southeast Europe compared with two West European panels

Southeast Europe (SEE) is a very important maize-growing region, comparable to the Corn belt region of the United States, with similar dent germplasm (dent by dent hybrids). Historically, this region has undergone several genetic material swaps, following the trends in the US, with one of the most significant swaps related to US aid programs after WWII. The imported accessions used to make double-cross hybrids were also mixed with previously adapted germplasm originating from several more distant OPVs, supporting the transition to single cross-breeding. Many of these materials were deposited at the Maize Gene Bank of the Maize Research Institute Zemun Polje (MRIZP) between the 1960s and 1980s. A part of this Gene Bank (572 inbreds) was genotyped with Affymetrix Axiom Maize Genotyping Array with 616,201 polymorphic variants. Data were merged with two other genotyping datasets with mostly European flint (TUM dataset) and dent (DROPS dataset) germplasm. The final pan-European dataset consisted of 974 inbreds and 460,243 markers. Admixture analysis showed seven ancestral populations representing European flint, B73/B14, Lancaster, B37, Wf9/Oh07, A374, and Iodent pools. Subpanel of inbreds with SEE origin showed a lack of Iodent germplasm, marking its historical context. Several signatures of selection were identified at chromosomes 1, 3, 6, 7, 8, 9, and 10. The regions under selection were mined for protein-coding genes and were used for gene ontology (GO) analysis, showing a highly significant overrepresentation of genes involved in response to stress. Our results suggest the accumulation of favorable allelic diversity, especially in the context of changing climate in the genetic resources of SEE.

Genome Scan of Rice Landrace Populations Collected Across Time Revealed Climate Changes' Selective Footprints in the Genes Network Regulating Flowering Time

Analyses of the genetic bases of plant adaptation to climate changes, using genome-scan approaches, are often conducted on natural populations, under hypothesis of out-crossing reproductive regime. We report here on a study based on diachronic sampling (1980 and 2011) of the autogamous crop species, Oryza sativa and Oryza glaberrima, in the tropical forest and the Sudanian savannah of West Africa. First, using historical meteorological data we confirmed changes in temperatures (+ 1 °C on average) and rainfall regime (less predictable and reduced amount) in the target areas. Second, phenotyping the populations for phenology, we observed significantly earlier heading time in the 2010 samples. Third, implementing two genome-scan methods (one of which specially developed for selfing species) on genotyping by sequencing genotypic data of the two populations, we detected 31 independent selection footprints. Gene ontology analysis detected significant enrichment of these selection footprints in genes involved in reproductive processes. Some of them bore known heading time QTLs and genes, including OsGI, Hd1 and OsphyB. This rapid adaptive evolution, originated from subtle changes in the standing variation in genetic network regulating heading time, did not translate into predominance of multilocus genotypes, as it is often the case in selfing plants, and into notable selective sweeps. The high adaptive potential observed results from the multiline genetic structure of the rice landraces, and the rather large and imbricated genetic diversity of the rice meta-population at the farm, the village and the region levels, that hosted the adaptive variants in multiple genetic backgrounds before the advent of the environmental selective pressure. Our results illustrate the evolution of in situ diversity through processes of human and natural selection, and provide a model for rice breeding and cultivars deployment strategies aiming resilience to climate changes. It also calls for further development of population genetic models for adaptation of plant populations to environmental changes. To our best knowledge, this is the first study dealing with climate-changes' selective footprint in crops.

Pollinator loss causes rapid adaptive evolution of selfing and dramatically reduces genome-wide genetic variability

Although selfing populations harbor little genetic variation limiting evolutionary potential, the causes are unclear. We experimentally evolved large, replicate populations of Mimulus guttatus for nine generations in greenhouses with or without pollinating bees and studied DNA polymorphism in descendants. Populations without bees adapted to produce more selfed seed yet exhibited striking reductions in DNA polymorphism despite large population sizes. Importantly, the genome-wide pattern of variation cannot be explained by a simple reduction in effective population size, but instead reflects the complicated interaction between selection, linkage, and inbreeding. Simulations demonstrate that the spread of favored alleles at few loci depresses neutral variation genome wide in large populations containing fully selfing lineages. It also generates greater heterogeneity among chromosomes than expected with neutral evolution in small populations. Genome-wide deviations from neutrality were documented in populations with bees, suggesting widespread influences of background selection. After applying outlier tests to detect loci under selection, two genome regions were found in populations with bees, yet no adaptive loci were otherwise mapped. Large amounts of stochastic change in selfing populations compromise evolutionary potential and undermine outlier tests for selection. This occurs because genetic draft in highly selfing populations makes even the largest changes in allele frequency unremarkable.

Genomic regions and candidate genes selected during the breeding of rice in Vietnam

Vietnam harnesses a rich diversity of rice landraces adapted to a range of conditions, which constitute a largely untapped source of diversity for the continuous improvement of cultivars. We previously identified a strong population structure in Vietnamese rice, which is captured in five Indica and four Japonica subpopulations, including an outlying Indica‐5 group. Here, we leveraged that strong differentiation and 672 native rice genomes to identify genomic regions and genes putatively selected during the breeding of rice in Vietnam. We identified significant distorted patterns in allele frequency (XP‐CLR) and population differentiation scores (F_ST) resulting from differential selective pressures between native subpopulations, and later annotated them with QTLs previously identified by GWAS in the same panel. We particularly focussed on the outlying Indica‐5 subpopulation because of its likely novelty and differential evolution, where we annotated 52 selected regions, which represented 8.1% of the rice genome. We annotated the 4576 genes in these regions and selected 65 candidate genes as promising breeding targets, several of which harboured alleles with nonsynonymous substitutions. Our results highlight genomic differences between traditional Vietnamese landraces, which are likely the product of adaption to multiple environmental conditions and regional culinary preferences in a very diverse country. We also verified the applicability of this genome scanning approach to identify potential regions harbouring novel loci and alleles to breed a new generation of sustainable and resilient rice.

Detecting selection using extended haplotype homozygosity (EHH)-based statistics in unphased or unpolarized data

Analysis of population genetic data often includes a search for genomic regions with signs of recent positive selection. One of such approaches involves the concept of extended haplotype homozygosity (EHH) and its associated statistics. These statistics typically require phased haplotypes, and some of them necessitate polarized variants. Here, we unify and extend previously proposed modifications to loosen these requirements. We compare the modified versions with the original ones by measuring the false discovery rate in simulated whole-genome scans and by quantifying the overlap of inferred candidate regions in empirical data. We find that phasing information is indispensable for accurate estimation of within-population statistics (for all but very large samples) and of cross-population statistics for small samples. Ancestry information, in contrast, is of lesser importance for both types of statistic. Our publicly available R package rehh incorporates the modified statistics presented here.

Physical geography, isolation by distance and environmental variables shape genomic variation of wild barley (Hordeum vulgare L. ssp. spontaneum) in the Southern Levant

Determining the extent of genetic variation that reflects local adaptation in crop-wild relatives is of interest for the purpose of identifying useful genetic diversity for plant breeding. We investigated the association of genomic variation with geographical and environmental factors in wild barley (Hordeum vulgare L. ssp. spontaneum) populations of the Southern Levant using genotyping by sequencing (GBS) of 244 accessions in the Barley 1K+ collection. The inference of population structure resulted in four genetic clusters that corresponded to eco-geographical habitats and a significant association between lower gene flow rates and geographical barriers, e.g. the Judaean Mountains and the Sea of Galilee. Redundancy analysis (RDA) revealed that spatial autocorrelation explained 45% and environmental variables explained 15% of total genomic variation. Only 4.5% of genomic variation was solely attributed to environmental variation if the component confounded with spatial autocorrelation was excluded. A synthetic environmental variable combining latitude, solar radiation, and accumulated precipitation explained the highest proportion of genomic variation (3.9%). When conditioned on population structure, soil water capacity was the most important environmental variable explaining 1.18% of genomic variation. Genome scans with outlier analysis and genome-environment association studies were conducted to identify adaptation signatures. RDA and outlier methods jointly detected selection signatures in the pericentromeric regions, which have reduced recombination, of the chromosomes 3H, 4H, and 5H. However, selection signatures mostly disappeared after correction for population structure. In conclusion, adaptation to the highly diverse environments of the Southern Levant over short geographical ranges had a limited effect on the genomic diversity of wild barley. This highlighted the importance of nonselective forces in genetic differentiation.

Migration without interbreeding: Evolutionary history of a highly selfing Mediterranean grass inferred from whole genomes

Wild plant populations show extensive genetic subdivision and are far from the ideal of panmixia which permeates population genetic theory. Understanding the spatial and temporal scale of population structure is therefore fundamental for empirical population genetics - and of interest in itself, as it yields insights into the history and biology of a species. In this study we extend the genomic resources for the wild Mediterranean grass Brachypodium distachyon to investigate the scale of population structure and its underlying history at whole-genome resolution. A total of 86 accessions were sampled at local and regional scales in Italy and France, which closes a conspicuous gap in the collection for this model organism. The analysis of 196 accessions, spanning the Mediterranean from Spain to Iraq, suggests that the interplay of high selfing and seed dispersal rates has shaped genetic structure in B. distachyon. At the continental scale, the evolution in B. distachyon is characterized by the independent expansion of three lineages during the Upper Pleistocene. Today, these lineages may occur on the same meadow yet do not interbreed. At the regional scale, dispersal and selfing interact and maintain high genotypic diversity, thus challenging the textbook notion that selfing in finite populations implies reduced diversity. Our study extends the population genomic resources for B. distachyon and suggests that an important use of this wild plant model is to investigate how selfing and dispersal, two processes typically studied separately, interact in colonizing plant species.

Genomic Signatures of Sexual Selection on Pollen-Expressed Genes in Arabis alpina

Fertilization in angiosperms involves the germination of pollen on the stigma, followed by the extrusion of a pollen tube that elongates through the style and delivers two sperm cells to the embryo sac. Sexual selection could occur throughout this process when male gametophytes compete for fertilization. The strength of sexual selection during pollen competition should be affected by the number of genotypes deposited on the stigma. As increased self-fertilization reduces the number of mating partners, and the genetic diversity and heterozygosity of populations, it should thereby reduce the intensity of sexual selection during pollen competition. Despite the prevalence of mating system shifts, few studies have directly compared the molecular signatures of sexual selection during pollen competition in populations with different mating systems. Here we analyzed whole-genome sequences from natural populations of Arabis alpina, a species showing mating system variation across its distribution, to test whether shifts from cross- to self-fertilization result in molecular signatures consistent with sexual selection on genes involved in pollen competition. We found evidence for efficient purifying selection on genes expressed in vegetative pollen, and overall weaker selection on sperm-expressed genes. This pattern was robust when controlling for gene expression level and specificity. In agreement with the expectation that sexual selection intensifies under cross-fertilization, we found that the efficacy of purifying selection on male gametophyte-expressed genes was significantly stronger in genetically more diverse and outbred populations. Our results show that intra-sexual competition shapes the evolution of pollen-expressed genes, and that its strength fades with increasing self-fertilization rates.

Understanding how an amphicarpic species with a mixed mating system responds to fire: a population genetic approach

Amphicarpic plants produce both above-ground and below-ground seeds. Because below-ground seeds are protected in the soil and may maintain viability when above-ground conditions are stressful, they were proposed as an adaptation to recolonize a site after disturbance. However, whether below-ground seeds are the main colonizers after a disturbance remains unknown. Our goal was to understand whether recolonization by an amphicarpic species after fire was accomplished primarily through germination of seeds produced above-ground or below-ground. We investigated Polygala lewtonii, an amphicarpic, perennial species endemic to fire-prone Florida sandhill and scrub, where fire kills plants but subsequently increases recruitment and population sizes. Polygala lewtonii produces three flower types: above-ground chasmogamous flowers and above-ground and below-ground cleistogamous flowers, with previous research demonstrating chasmogamous flowers produce a much greater proportion of seeds than cleistogamous flowers. We quantified outcrossing in seeds produced by chasmogamous flowers to determine whether it differed from the 100 % self-fertilized below-ground seeds. Approximately 25 % of seeds from chasmogamous flowers showed evidence of cross-pollination. Assuming that chasmogamous flowers produce the majority of the above-ground seeds, as was shown previously, this indicates it is possible to differentiate between germination by above-ground versus below-ground seeds in post-fire colonization. We next compared genetic diversity, admixture, inbreeding and population genetic structure pre- and post-fire. If fire promoted germination of chasmogamous seeds, heterozygosity and admixture would increase, and genetic structure and inbreeding would decrease. Instead, inbreeding and genetic structure increased and admixture decreased, suggesting that the below-ground selfed seeds (with limited dispersal ability) increased their contribution to the population after fire, possibly because fire reduced above-ground seed viability. Additionally, new alleles not found previously in range-wide analyses emerged from the seed bank post-fire. These results suggest that amphicarpy is a powerful adaptation to preserve genetic variation, maintain adaptive potential and promote rapid post-fire colonization.

Interplay between extreme drift and selection intensities favors the fixation of beneficial mutations in selfing maize populations

Population and quantitative genetic models provide useful approximations to predict long-term selection responses sustaining phenotypic shifts, and underlying multilocus adaptive dynamics. Valid across a broad range of parameters, their use for understanding the adaptive dynamics of small selfing populations undergoing strong selection intensity (thereafter High Drift-High selection regime, HDHS) remains to be explored. Saclay Divergent Selection Experiments (DSEs) on maize flowering time provide an interesting example of populations evolving under HDHS, with significant selection responses over 20 generations in two directions. We combined experimental data from Saclay DSEs, forward individual-based simulations, and theoretical predictions to dissect the evolutionary mechanisms at play in the observed selection responses. We asked two main questions: How do mutations arise, spread, and reach fixation in populations evolving under HDHS? How does the interplay between drift and selection influence observed phenotypic shifts? We showed that the long-lasting response to selection in small populations is due to the rapid fixation of mutations occurring during the generations of selection. Among fixed mutations, we also found a clear signal of enrichment for beneficial mutations revealing a limited cost of selection. Both environmental stochasticity and variation in selection coefficients likely contributed to exacerbate mutational effects, thereby facilitating selection grasp and fixation of small-effect mutations. Together our results highlight that despite a small number of polymorphic loci expected under HDHS, adaptive variation is continuously fueled by a vast mutational target. We discuss our results in the context of breeding and long-term survival of small selfing populations.

Arabis alpina: A perennial model plant for ecological genomics and life-history evolution

Many model organisms were chosen and achieved prominence because of an advantageous combination of their life‐history characteristics, genetic properties and also practical considerations. Discoveries made in Arabidopsis thaliana, the most renowned noncrop plant model species, have markedly stimulated studies in other species with different biology. Within the family Brassicaceae, the arctic–alpine Arabis alpina has become a model complementary to Arabidopsis thaliana to study the evolution of life‐history traits, such as perenniality, and ecological genomics in harsh environments. In this review, we provide an overview of the properties that facilitated the rapid emergence of A. alpina as a plant model. We summarize the evolutionary history of A. alpina, including genomic aspects, the diversification of its mating system and demographic properties, and we discuss recent progress in the molecular dissection of developmental traits that are related to its perennial life history and environmental adaptation. From this published knowledge, we derive open questions that might inspire future research in A. alpina, other Brassicaceae species or more distantly related plant families.

The integrative biology of genetic dominance

Dominance is a basic property of inheritance systems describing the link between a diploid genotype at a single locus and the resulting phenotype. Models for the evolution of dominance have long been framed as an opposition between the irreconcilable views of Fisher in 1928 supporting the role of largely elusive dominance modifiers and Wright in 1929, who viewed dominance as an emerging property of the structure of enzymatic pathways. Recent theoretical and empirical advances however suggest that these opposing views can be reconciled, notably using models investigating the regulation of gene expression and developmental processes. In this more comprehensive framework, phenotypic dominance emerges from departures from linearity between any levels of integration in the genotype‐to‐phenotype map. Here, we review how these different models illuminate the emergence and evolution of dominance. We then detail recent empirical studies shedding new light on the diversity of molecular and physiological mechanisms underlying dominance and its evolution. By reconciling population genetics and functional biology, we hope our review will facilitate cross‐talk among research fields in the integrative study of dominance evolution.

Genetic diversity and population structure of a global invader Mayweed chamomile (Anthemis cotula): management implications

Mayweed chamomile (Anthemis cotula) is a globally invasive, troublesome annual weed but knowledge of its genetic diversity, population structure in invaded regions and invasion patterns remains unstudied. Therefore, germplasm from 19 A. cotula populations (sites) from three geographically distinct invaded regions: the Walla Walla Basin (located in southern Washington) and the Palouse (located in both northern Idaho and eastern Washington), Pacific Northwest, USA and Kashmir Valley, India were grown in the greenhouse for DNA extraction and sequencing. A total of 18 829 single-nucleotide polymorphisms were called and filtered for each of 89 samples. Pairwise F ST, Nei's genetic distance, heterozygosity, Wright's inbreeding coefficient (F) and self-fertilization rates were estimated for populations within and among the three regions with a total of 19 populations comprised of 89 individuals. Overall measurements of genetic variation were low but significant among regions, populations and individuals. Despite the weak genetic structure, two main genetic clusters were evident, one comprised of populations from Palouse and Kashmir Valley, the other comprised of populations from the Walla Walla Basin. Significant selfing was observed in populations from the Walla Walla Basin and Palouse but not from Kashmir Valley, indicating that Mayweed chamomile in the Pacific Northwest, USA could persist with low pollinator or pollen donor densities. Although F ST values between the regions indicate Palouse populations are more closely related to Kashmir Valley than to Walla Walla Basin populations, based on Migrate-n analysis, panmixis was the most likely model, suggesting an unrestricted gene flow among all three regions. Our study indicated that Kashmir Valley populations either originated from or shared the origin with the Palouse populations, suggesting human-mediated migration of A. cotula between regions.

Evidence of spontaneous selfing and disomic inheritance in Geranium robertianum

Knowing species' breeding system and mating processes occurring in populations is important not only for understanding population dynamics, gene flow processes, and species' response to climate change, but also for designing control plans of invasive species. Geranium robertianum, a widespread biennial herbaceous species showing high morphological variation and wide ecological amplitude, can become invasive outside its distribution range. A mixed-mating system may be expected given the species' floral traits. However, autonomous selfing is considered as a common feature. Genetic variation and structure, and so population mating processes, have not been investigated in wild populations. We developed 15 polymorphic microsatellite markers to quantify genetic variation and structure in G. robertianum. To investigate whether selfing might be the main mating process in natural conditions, we sampled three generations of plants (adult, F1, and F2) for populations from the UK, Spain, Belgium, Germany, and Sweden, and compared open-pollinated with outcrossed hand-pollinated F2 progeny. The highly positive Wright's inbreeding coefficient (F IS) values in adults, F1, and open-pollinated F2 progeny and the low F IS values in outcross F2 progeny supported autonomous selfing as the main mating process for G. robertianum in wild conditions, despite the presence of attractive signals for insect pollination. Genetic differentiation among samples was found, showing some western-eastern longitudinal trend. Long-distance seed dispersal might have contributed to the low geographic structure. Local genetic differentiation may have resulted not only from genetic drift effects favored by spontaneous selfing, but also from ecological adaptation. The presence of duplicate loci with disomic inheritance is consistent with the hypothesis of allotetraploid origin of G. robertianum. The fact that most microsatellite markers behave as diploid loci with no evidence of duplication supports the hypothesis of ancient polyploidization. The differences in locus duplication and the relatively high genetic diversity across G. robertianum range despite spontaneous autonomous selfing suggest multiple events of polyploidization.

Model design artificial intelligence and research of adaptive network intrusion detection and defense system using fuzzy logic

Along with improvement of technology in network and continuous expansion of network economy and network applications, the Internet has gradually become an indispensable part of the modern society. However, an endless stream of hacker attacks and network virus events make network security issues stand out. Therefore, network security has become a hot spot in computer network research and development. This paper aims at establishing a real-time detection and dynamic defense security system and makes an in-depth study of intrusion detection technology and defense decision-making technology. The strategy involved in finding the intrusion behavior since the fuzzy base contains the better group of rules. We have utilized an automated fuzzy rule generation strategy. An adaptive network intrusion detection and defense system model is established, and the architecture of the model is discussed in detail. The platform independence, good self-adaptability, expansibility, multi-level data analysis and dynamic defense decision-making are expounded. The experiment proves that the model proposed in this article has a good self-adaptability and open construction, and effectively combines the functions of intrusion detection and defense decision-making.

Drivers of linkage disequilibrium across a species' geographic range

While linkage disequilibrium (LD) is an important parameter in genetics and evolutionary biology, the drivers of LD remain elusive. Using whole-genome sequences from across a species’ range, we assessed the impact of demographic history and mating system on LD. Both range expansion and a shift from outcrossing to selfing in North American Arabidopsis lyrata were associated with increased average genome-wide LD. Our results indicate that range expansion increases short-distance LD at the farthest range edges by about the same amount as a shift to selfing. However, the extent over which LD in genic regions unfolds was shorter for range expansion compared to selfing. Linkage among putatively neutral variants and between neutral and deleterious variants increased to a similar degree with range expansion, providing support that genome-wide LD was positively associated with mutational load. As a consequence, LD combined with mutational load may decelerate range expansions and set range limits. Finally, a small number of genes were identified as LD outliers, suggesting that they experience selection by either of the two demographic processes. These included genes involved in flowering and photoperiod for range expansion, and the self-incompatibility locus for mating system.

Nearby genomic variants are often co-inherited because of limited recombination. The extent of non-random association of alleles at different loci is called linkage disequilibrium (LD) and is commonly used in genomic analyses, for example to detect regions under selection or to determine effective population size. Here we reversed testing and addressed how demographic history may affect LD within a species. Using genomic data from more than a thousand individuals of North American Arabidopsis lyrata from across the entire species’ range, we quantified the effect of postglacial range expansion and a shift in mating system from outcrossing to selfing on LD. We show that both factors lead to increased LD, and that the maximal effect of range expansion is comparable with a shift in mating system to selfing. Heightened LD involves deleterious mutations, and therefore, LD can also serve as an indicator of mutation accumulation. Furthermore, we provide evidence that some genes experienced stronger increases in LD possibly due to selection associated with the two demographic changes. Our results provide a novel and broad view on the evolutionary factors shaping LD that may also apply to the very many species that underwent postglacial range expansion.

Genomic characterization of selfing in the cyclic parthenogen Daphnia magna

Inbreeding refers to the fusion of related individuals' gametes, with self-fertilization (selfing) being an extreme form of inbreeding-involving gametes produced by the same individual. Selfing is expected to reduce heterozygosity by an average of 50% in one generation; however, little is known about the empirical variation on a genome level surrounding this figure and the factors that affect variation. We selfed genotypes of the cyclic parthenogen Daphnia magna and analysed whole genomes of mothers and selfed offspring, observing the predicted 50% heterozygosity reduction on average. We also saw substantial variation around this value and significant differences among mother-offspring pairs. Crossover analysis confirmed the known trend of recombination occurring more often towards the telomeres. This effect was shown, through simulations, to increase the variance of heterozygosity reduction compared to when a uniform distribution of crossovers was used. Similarly, we simulated inbred line production after several generations of selfing and we observed higher variance in achieved homozygosity when we consider a higher recombination rate towards the telomeres. Our empirical and simulation study highlights that the expected mean values of heterozygosity reduction show remarkable variation, which can help understand, for example, differences among inbred individuals.

Flower evolution in the presence of heterospecific gene flow and its contribution to lineage divergence

The success of species depends on their ability to exploit ecological resources in order to optimize their reproduction. However, species are not usually found within single-species ecosystems but in complex communities. Because of their genetic relatedness, closely related lineages tend to cluster within the same ecosystem, rely on the same resources, and be phenotypically similar. In sympatry, they will therefore compete for the same resources and, in the case of flowering plants, exchange their genes through heterospecific pollen transfer. These interactions, nevertheless, pose significant challenges to species co-existence because they can lead to resource limitation and reproductive interference. In such cases, divergent selective pressures on floral traits will favour genotypes that isolate or desynchronize the reproduction of sympatric lineages. The resulting displacement of reproductive characters will, in turn, lead to pre-mating isolation and promote intraspecific divergence, thus initiating or reinforcing the speciation process. In this review, we discuss the current theoretical and empirical knowledge on the influence of heterospecific pollen transfer on flower evolution, highlighting its potential to uncover the ecological and genomic constraints shaping the speciation process.

The patterns of deleterious mutations during the domestication of soybean

Globally, soybean is a major protein and oil crop. Enhancing our understanding of the soybean domestication and improvement process helps boost genomics-assisted breeding efforts. Here we present a genome-wide variation map of 10.6 million single-nucleotide polymorphisms and 1.4 million indels for 781 soybean individuals which includes 418 domesticated (Glycine max), 345 wild (Glycine soja), and 18 natural hybrid (G. max/G. soja) accessions. We describe the enhanced detection of 183 domestication-selective sweeps and the patterns of putative deleterious mutations during domestication and improvement. This predominantly selfing species shows 7.1% reduction of overall deleterious mutations in domesticated soybean relative to wild soybean and a further 1.4% reduction from landrace to improved accessions. The detected domestication-selective sweeps also show reduced levels of deleterious alleles. Importantly, genotype imputation with this resource increases the mapping resolution of genome-wide association studies for seed protein and oil traits in a soybean diversity panel.

Identifying the fitness consequences of sex in complex natural environments

In the natural world, sex prevails, despite its costs. Although much effort has been dedicated to identifying the intrinsic costs of sex (e.g., the cost of males), few studies have identified the ecological fitness consequences of sex. Furthermore, correlated biological traits that differ between sexuals and asexuals may alter these costs, or even render the typical costs of sex irrelevant. We conducted a large-scale, multisite, reciprocal transplant using multiple sexual and asexual genotypes of a native North American wildflower to show that sexual genotypes have reduced lifetime fitness, despite lower herbivory. We separated the effects of sex from those of hybridity, finding that overwinter survival is elevated in asexuals regardless of hybridity, but herbivores target hybrid asexuals more than nonhybrid asexual or sexual genotypes. Survival is lowest in homozygous sexual lineages, implicating inbreeding depression as a cost of sex. Our results show that the consequences of sex are shaped not just by sex itself, but by complex natural environments, correlated traits, and the identity and availability of mates.

Expressed mutational load increases toward the edge of a species' geographic range

There is no general explanation for why species have restricted geographic distributions. One hypothesis posits that range expansion or increasing scarcity of suitable habitat results in accumulation of mutational load due to enhanced genetic drift, which constrains population performance toward range limits and further expansion. We tested this hypothesis in the North American plant, Arabidopsis lyrata. We experimentally assessed mutational load by crossing plants of 20 populations from across the entire species range and by raising the offspring of within- and between-population crosses at five common garden sites within and beyond the range. Offspring performance was tracked over three growing seasons. The heterosis effect, depicting expressed mutational load, was increased in populations with heightened genomic estimates of load, longer expansion distance or long-term isolation, and a selfing mating system. The decline in performance of within-population crosses amounted to 80%. Mutation accumulation due to past range expansion and long-term isolation of populations in the area of range margins is therefore a strong determinant of population-mean performance, and the magnitude of effect may be sufficient to cause range limits.

Selective Sweeps Under Dominance and Inbreeding

A major research goal in evolutionary genetics is to uncover loci experiencing positive selection. One approach involves finding 'selective sweeps' patterns, which can either be 'hard sweeps' formed by de novo mutation, or 'soft sweeps' arising from recurrent mutation or existing standing variation. Existing theory generally assumes outcrossing populations, and it is unclear how dominance affects soft sweeps. We consider how arbitrary dominance and inbreeding via self-fertilization affect hard and soft sweep signatures. With increased self-fertilization, they are maintained over longer map distances due to reduced effective recombination and faster beneficial allele fixation times. Dominance can affect sweep patterns in outcrossers if the derived variant originates from either a single novel allele, or from recurrent mutation. These models highlight the challenges in distinguishing hard and soft sweeps, and propose methods to differentiate between scenarios.

Postglacial ecotype formation under outcrossing and self-fertilization in Arabidopsis lyrata

The formation of ecotypes has been invoked as an important driver of postglacial biodiversity, because many species colonized heterogeneous habitats and experienced divergent selection. Ecotype formation has been predominantly studied in outcrossing taxa, while far less attention has been paid to the implications of mating system shifts. Here, we addressed whether substrate-related ecotypes exist in selfing and outcrossing populations of Arabidopsis lyrata subsp. lyrata and whether the genomic footprint differs between mating systems. The North American subspecies colonized both rocky and sandy habitats during postglacial range expansion and shifted the mating system from predominantly outcrossing to predominantly selfing in a number of regions. We performed an association study on pooled whole-genome sequence data of 20 selfing or outcrossing populations, which suggested genes involved in adaptation to substrate. Motivated by enriched gene ontology terms, we compared root growth between plants from the two substrates in a common environment and found that plants originating from sand grew roots faster and produced more side roots, independent of mating system. Furthermore, single nucleotide polymorphisms associated with substrate-related ecotypes were more clustered among selfing populations. Our study provides evidence for substrate-related ecotypes in A. lyrata and divergence in the genomic footprint between mating systems. The latter is the likely result of selfing populations having experienced divergent selection on larger genomic regions due to higher genome-wide linkage disequilibrium.

Mating system impacts the genetic architecture of adaptation to heterogeneous environments

Self-fertilisation has consequences for variation across the genome as it reduces effective population size, effect recombination rates and pollen flow, with implications for local adaptation. We conducted simulations of divergent stabilising selection on a quantitative trait with drift, pollen flow, mutation, recombination and different outcrossing rates. We quantified trait divergence and the genetic architecture of adaptation. We conducted an F_ST outlier analysis to identify candidate loci and quantified the impact of mating system on detectability. Selfing promoted trait divergence mainly through reductions in pollen flow. Moreover, trait architecture became more diffuse with selfing. Average effect size of trait loci was lower, while the number of loci, and their clustering distance increased. The genetic architecture of selfers was also more diffuse than outcrossers for equivalent migration rates. However, when deleterious alleles were included, architectures became more concentrated in selfers, likely to be because of reductions in population size caused by mutational meltdown and impacts of background selection on N_e . Our simulations demonstrate that mating system has important impacts on adaptive divergence of traits and the genetic landscape underlying that divergence. Selfing has a significant effect on detectability of regions of the genome important for adaptation because of neutral divergence and diffuse trait architecture.

How relatedness between mates influences reproductive success: An experimental analysis of self-fertilization and biparental inbreeding in a marine bryozoan

Kin associations increase the potential for inbreeding. The potential for inbreeding does not, however, make inbreeding inevitable. Numerous factors influence whether inbreeding preference, avoidance, or tolerance evolves, and, in hermaphrodites where both self-fertilization and biparental inbreeding are possible, it remains particularly difficult to predict how selection acts on the overall inbreeding strategy, and to distinguish the type of inbreeding when making inferences from genetic markers. Therefore, we undertook an empirical analysis on an understudied type of mating system (spermcast mating in the marine bryozoan, Bugula neritina) that provides numerous opportunities for inbreeding preference, avoidance, and tolerance. We created experimental crosses, containing three generations from two populations to estimate how parental reproductive success varies across parental relatedness, ranging from self, siblings, and nonsiblings from within the same population. We found that the production of viable selfed offspring was extremely rare (only one colony produced three selfed offspring) and biparental inbreeding more common. Paternity analysis using 16 microsatellite markers confirmed outcrossing. The production of juveniles was lower for sib mating compared with nonsib mating. We found little evidence for consistent inbreeding, in terms of nonrandom mating, in adult samples collected from three populations, using multiple population genetic inferences. Our results suggest several testable hypotheses that potentially explain the overall mating and dispersal strategy in this species, including early inbreeding depression, inbreeding avoidance through cryptic mate choice, and differential dispersal distances of sperm and larvae.

A high-throughput amplicon-based method for estimating outcrossing rates

BACKGROUND

The outcrossing rate is a key determinant of the population-genetic structure of species and their long-term evolutionary trajectories. However, determining the outcrossing rate using current methods based on PCR-genotyping individual offspring of focal plants for multiple polymorphic markers is laborious and time-consuming.

RESULTS

We have developed an amplicon-based, high-throughput enabled method for estimating the outcrossing rate and have applied this to an example of scented versus non-scented Capsella (Shepherd's Purse) genotypes. Our results show that the method is able to robustly capture differences in outcrossing rates. They also highlight potential biases in the estimates resulting from differential haplotype sharing of the focal plants with the pollen-donor population at individual amplicons.

CONCLUSIONS

This novel method for estimating outcrossing rates will allow determining this key population-genetic parameter with high-throughput across many genotypes in a population, enabling studies into the genetic determinants of successful pollinator attraction and outcrossing.

Rapid functional divergence after small-scale gene duplication in grasses

BACKGROUND

Gene duplication has played an important role in the evolution and domestication of flowering plants. Yet little is known about how plant duplicate genes evolve and are retained over long timescales, particularly those arising from small-scale duplication (SSD) rather than whole-genome duplication (WGD) events.

RESULTS

We address this question in the Poaceae (grass) family by analyzing gene expression data from nine tissues of Brachypodium distachyon, Oryza sativa japonica (rice), and Sorghum bicolor (sorghum). Consistent with theoretical predictions, expression profiles of most grass genes are conserved after SSD, suggesting that functional conservation is the primary outcome of SSD in grasses. However, we also uncover support for widespread functional divergence, much of which occurs asymmetrically via the process of neofunctionalization. Moreover, neofunctionalization preferentially targets younger (child) duplicate gene copies, is associated with RNA-mediated duplication, and occurs quickly after duplication. Further analysis reveals that functional divergence of SSD-derived genes is positively correlated with both sequence divergence and tissue specificity in all three grass species, and particularly with anther expression in B. distachyon.

CONCLUSIONS

Our results suggest that SSD-derived grass genes often undergo rapid functional divergence that may be driven by natural selection on male-specific phenotypes. These observations are consistent with those in several animal species, suggesting that duplicate genes take similar evolutionary trajectories in plants and animals.

Parental legacy, demography, and admixture influenced the evolution of the two subgenomes of the tetraploid Capsella bursa-pastoris (Brassicaceae)

Allopolyploidy is generally perceived as a major source of evolutionary novelties and as an instantaneous way to create isolation barriers. However, we do not have a clear understanding of how two subgenomes evolve and interact once they have fused in an allopolyploid species nor how isolated they are from their relatives. Here, we address these questions by analyzing genomic and transcriptomic data of allotetraploid Capsella bursa-pastoris in three differentiated populations, Asia, Europe, and the Middle East. We phased the two subgenomes, one descended from the outcrossing and highly diverse Capsella grandiflora (CbpCg) and the other one from the selfing and genetically depauperate Capsella orientalis (CbpCo). For each subgenome, we assessed its relationship with the diploid relatives, temporal changes of effective population size (Ne), signatures of positive and negative selection, and gene expression patterns. In all three regions, Ne of the two subgenomes decreased gradually over time and the CbpCo subgenome accumulated more deleterious changes than CbpCg. There were signs of widespread admixture between C. bursa-pastoris and its diploid relatives. The two subgenomes were impacted differentially depending on geographic region suggesting either strong interploidy gene flow or multiple origins of C. bursa-pastoris. Selective sweeps were more common on the CbpCg subgenome in Europe and the Middle East, and on the CbpCo subgenome in Asia. In contrast, differences in expression were limited with the CbpCg subgenome slightly more expressed than CbpCo in Europe and the Middle-East. In summary, after more than 100,000 generations of co-existence, the two subgenomes of C. bursa-pastoris still retained a strong signature of parental legacy but their evolutionary trajectory strongly varied across geographic regions.

Structure of multilocus genetic diversity in predominantly selfing populations

Predominantly selfing populations are expected to have reduced effective population sizes due to nonrandom sampling of gametes, demographic stochasticity (bottlenecks or extinction-recolonization), and large scale hitchhiking (reduced effective recombination). Thus, they are expected to display low genetic diversity, which was confirmed by empirical studies. The structure of genetic diversity in predominantly selfing species is dramatically different from outcrossing ones, with populations often dominated by one or a few multilocus genotypes (MLGs) coexisting with several rare genotypes. Therefore, multilocus diversity indices are relevant to describe diversity in selfing populations. Here, we use simulations to provide analytical expectations for multilocus indices and examine whether selfing alone can be responsible for the high-frequency MLGs persistent through time in the absence of selection. We then examine how combining single and multilocus indices of diversity may be insightful to distinguish the effects of selfing, population size, and more complex demographic events (bottlenecks, migration, admixture, or extinction-recolonization). Finally, we examine how temporal changes in MLG frequencies can be insightful to understand the evolutionary trajectory of a given population. We show that combinations of selfing and small demographic sizes can result in high-frequency MLGs, as observed in natural populations. We also show how different demographic scenarios can be distinguished by the parallel analysis of single and multilocus indices of diversity, and we emphasize the importance of temporal data for the study of predominantly selfing populations. Finally, the comparison of our simulations with empirical data on populations of Medicago truncatula confirms the pertinence of our simulation framework.

polyDFEv2.0: testing for invariance of the distribution of fitness effects within and across species

Summary

Distribution of fitness effects (DFE) of mutations can be inferred from site frequency spectrum (SFS) data. There is mounting interest to determine whether distinct genomic regions and/or species share a common DFE, or whether evidence exists for differences among them. polyDFEv2.0 fits multiple SFS datasets at once and provides likelihood ratio tests for DFE invariance across datasets. Simulations show that testing for DFE invariance across genomic regions within a species requires models accounting for distinct sources of heterogeneity (chance and genuine difference in DFE) underlying differences in SFS data in these regions. Not accounting for this will result in the spurious detection of DFE differences.

Availability and Implementation

polyDFEv2.0 is implemented in C and is accompanied by a series of R functions that facilitate post-processing of the output. It is available as source code and compiled binaries under a GNU General Public License v3.0 from https://github.com/paula-tataru/polyDFE.

Supplementary information

Supplementary data are available at Bioinformatics online.

Sexual selection and inbreeding: Two efficient ways to limit the accumulation of deleterious mutations

Theory and empirical data showed that two processes can boost selection against deleterious mutations, thus facilitating the purging of the mutation load: inbreeding, by exposing recessive deleterious alleles to selection in homozygous form, and sexual selection, by enhancing the relative reproductive success of males with small mutation loads. These processes tend to be mutually exclusive because sexual selection is reduced under mating systems that promote inbreeding, such as self‐fertilization in hermaphrodites. We estimated the relative efficiency of inbreeding and sexual selection at purging the genetic load, using 50 generations of experimental evolution, in a hermaphroditic snail (Physa acuta). To this end, we generated lines that were exposed to various intensities of inbreeding, sexual selection (on the male function), and nonsexual selection (on the female function). We measured how these regimes affected the mutation load, quantified through the survival of outcrossed and selfed juveniles. We found that juvenile survival strongly decreased in outbred lines with reduced male selection, but not when female selection was relaxed, showing that male‐specific sexual selection does purge deleterious mutations. However, in lines exposed to inbreeding, where sexual selection was also relaxed, survival did not decrease, and even increased for self‐fertilized juveniles, showing that purging through inbreeding can compensate for the absence of sexual selection. Our results point to the further question of whether a mixed strategy combining the advantages of both mechanisms of genetic purging could be evolutionary stable.

The interaction between sex-specific selection and local adaptation in species without separate sexes

Local adaptation in hermaphrodite species can be based on a variety of fitness components, including survival, as well as both female and male sex-functions within individuals. When selection via female and male fitness components varies spatially (e.g. due to environmental heterogeneity), local adaptation will depend, in part, on variation in selection through each fitness component, and the extent to which genetic trade-offs between sex-functions maintain genetic variation necessary for adaptation. Local adaptation will also depend on the hermaphrodite mating system because self-fertilization alters several key factors influencing selection and the maintenance of genetic variance underlying trade-offs between the sex-functions (sexually antagonistic polymorphism). As a first step to guide intuition regarding sex-specific adaptation in hermaphrodites, we develop a simple theoretical model incorporating the essential features of hermaphrodite mating and adaptation in a spatially heterogeneous environment, and explore the interaction between sex-specific selection, self-fertilization and local adaptation. Our results suggest that opportunities for sex-specific local adaptation in hermaphrodites depend strongly on the extent of self-fertilization and inbreeding depression. Using our model as a conceptual framework, we provide a broad overview of the literature on sex-specific selection and local adaptation in hermaphroditic plants and animals, emphasizing promising future directions in light of our theoretical predictions.This article is part of the theme issue 'Linking local adaptation with the evolution of sex differences'.

Evolution of the selfing syndrome: Anther orientation and herkogamy together determine reproductive assurance in a self-compatible plant

Capacity for autonomous self-fertilization provides reproductive assurance, has evolved repeatedly in the plant kingdom, and typically involves several changes in flower morphology and development (the selfing syndrome). Yet, the relative importance of different traits and trait combinations for efficient selfing and reproductive success in pollinator-poor environments is poorly known. In a series of experiments, we tested the importance of anther-stigma distance and the less studied trait anther orientation for efficiency of selfing in the perennial herb Arabis alpina. Variation in flower morphology among eight self-compatible European populations was correlated with efficiency of self-pollination and with pollen limitation in a common-garden experiment. To examine whether anther-stigma distance and anther orientation are subject to directional and/or correlational selection, and whether this is because these traits affect pollination success, we planted a segregating F2 population at two native field sites. Selection strongly favored a combination of introrse anthers and reduced anther-stigma distance at a site where pollinator activity was low, and supplemental hand-pollination demonstrated that this was largely because of their effect on securing self-pollination. The results suggest that concurrent shifts in more than one trait can be crucial for the evolution of efficient self-pollination and reproductive assurance in pollinator-poor habitats.