Evidence for inter-population variation in waiting times in a self-fertilizing flatworm

Propensity for selfing varies within a population of hermaphroditic snails: coexistence of selfers, outcrossers and mixed-mating individuals

To understand mating-system evolution in self-compatible hermaphrodites, variation in selfing rates is highly relevant. Empirical studies are rarely designed to capture variation between individuals, instead often comparing species and populations. Yet, evolution primarily occurs within populations, rendering among-individual variation essential. Observed individual selfing rates depend on the environment (e.g. differences in mate availability) and individuals' propensity for selfing. We quantified individual variation in selfing propensity in the snail Radix balthica by conducting laboratory mating trials that manipulated mate availability (low versus moderate) and estimating selfing rates from progeny arrays. We also measured female lifetime fitness. We found substantial among-individual variation in selfing propensity, including pure selfers (32%), pure outcrossers (31%) and mixed-mating individuals that selfed and outcrossed (37%). Experimental levels of mate availability did not significantly affect selfing rates. Selfers had reduced female liftetime fitness. Our results show that the propensity for selfing can differ considerably among individuals, with similar proportions of selfers, outcrossers and mixed maters. As mate availability did not affect selfing, our 'moderate' experimental level of mate availability might still have been too low to prompt selfers to outcross. This and the observed fitness differences also cautiously suggest that investigating the heritability of selfing propensities might be worthwhile in this population.

Genome assemblies of the simultaneously hermaphroditic flatworms Macrostomum cliftonense and Macrostomum hystrix

The free-living, simultaneously hermaphroditic flatworms of the genus Macrostomum are increasingly used as model systems in various contexts. In particular, Macrostomum lignano, the only species of this group with a published genome assembly, has emerged as a model for the study of regeneration, reproduction, and stem-cell function. However, challenges have emerged due to M. lignano being a hidden polyploid, having recently undergone whole-genome duplication and chromosome fusion events. This complex genome architecture presents a significant roadblock to the application of many modern genetic tools. Hence, additional genomic resources for this genus are needed. Here, we present such resources for Macrostomum cliftonense and Macrostomum hystrix, which represent the contrasting mating behaviors of reciprocal copulation and hypodermic insemination found in the genus. We use a combination of PacBio long-read sequencing and Illumina shot-gun sequencing, along with several RNA-Seq data sets, to assemble and annotate highly contiguous genomes for both species. The assemblies span ∼227 and ∼220 Mb and are represented by 399 and 42 contigs for M. cliftonense and M. hystrix, respectively. Furthermore, high BUSCO completeness (∼84-85%), low BUSCO duplication rates (8.3-6.2%), and low k-mer multiplicity indicate that these assemblies do not suffer from the same assembly ambiguities of the M. lignano genome assembly, which can be attributed to the complex karyology of this species. We also show that these resources, in combination with the prior resources from M. lignano, offer an excellent foundation for comparative genomic research in this group of organisms.

Waiting for love but not forever: Modeling the evolution of waiting time to selfing in hermaphrodites

Although mixed mating systems involving both selfing and outcrossing are fairly common in hermaphrodites, the mechanisms maintaining mixed mating are still unknown in many cases. In some species, individuals that have not yet found a mating partner delay self-fertilization for some time. This “waiting time” to selfing (WT) can exhibit heritable variation between individuals and is subject to two opposing selection pressures: waiting longer increases the density-dependent probability to encounter a mate within that time and thereby the chance to avoid inbreeding depression (ID) in offspring, but also increases the risk of dying before reproduction. It has long been hypothesized that fluctuations in population density and thus mate availability can lead to stable intermediate WTs, but to our knowledge there are so far no quantitative models that also take into account the joint evolutionary dynamics of ID. We use an individual-based model and a mathematical approximation to explore how delayed selfing evolves in response to density and density fluctuations. We find that at high density, when individuals meet often, WT evolution is dominated by genetic drift; at intermediate densities, strong ID causes WT to increase; and at low densities, ID is purged and WT approaches zero. Positive feedback loops drive the system to either complete selfing or complete outcrossing. Fluctuating density can slow down convergence to these alternative stable states. However, mixed mating, in the sense of either a stable polymorphism in WT, or stable intermediate waiting times, was never observed. Thus, additional factors need to be explored to explain the persistence of delayed selfing.

Evolution of sex allocation plasticity in a hermaphroditic flatworm genus

Sex allocation theory in simultaneous hermaphrodites predicts that optimal sex allocation is influenced by local sperm competition, which occurs when related sperm compete to fertilize a given set of eggs. Different factors, including the mating strategy and the ability to self‐fertilize, are predicted to affect local sperm competition and hence the optimal SA. Moreover, since the local sperm competition experienced by an individual can vary temporally and spatially, this can favour the evolution of sex allocation plasticity. Here, using seven species of the free‐living flatworm genus Macrostomum, we document interspecific variation in sex allocation, but neither their mating strategy nor their ability to self‐fertilize significantly predicted sex allocation among these species. Since we also found interspecific variation in sex allocation plasticity, we further estimated standardized effect sizes for plasticity in response to (i) the presence of mating partners (i.e. in isolation vs. with partners) and (ii) the strength of local sperm competition (i.e. in small vs. large groups). We found that self‐fertilization predicted sex allocation plasticity with respect to the presence of mating partners, with plasticity being lower for self‐fertilizing species. Finally, we showed that interspecific variation in sex allocation is higher than intraspecific variation due to sex allocation plasticity. Our study suggests that both sex allocation and sex allocation plasticity are evolutionarily labile, with self‐fertilization predicting the latter in Macrostomum.

Faster Rates of Molecular Sequence Evolution in Reproduction-Related Genes and in Species with Hypodermic Sperm Morphologies

Sexual selection drives the evolution of many striking behaviors and morphologies and should leave signatures of selection at loci underlying these phenotypes. However, although loci thought to be under sexual selection often evolve rapidly, few studies have contrasted rates of molecular sequence evolution at such loci across lineages with different sexual selection contexts. Furthermore, work has focused on separate sexed animals, neglecting alternative sexual systems. We investigate rates of molecular sequence evolution in hermaphroditic flatworms of the genus Macrostomum. Specifically, we compare species that exhibit contrasting sperm morphologies, strongly associated with multiple convergent shifts in the mating strategy, reflecting different sexual selection contexts. Species donating and receiving sperm in every mating have sperm with bristles, likely to prevent sperm removal. Meanwhile, species that hypodermically inject sperm lack bristles, potentially as an adaptation to the environment experienced by hypodermic sperm. Combining functional annotations from the model, Macrostomum lignano, with transcriptomes from 93 congeners, we find genus-wide faster sequence evolution in reproduction-related versus ubiquitously expressed genes, consistent with stronger sexual selection on the former. Additionally, species with hypodermic sperm morphologies had elevated molecular sequence evolution, regardless of a gene's functional annotation. These genome-wide patterns suggest reduced selection efficiency following shifts to hypodermic mating, possibly due to higher selfing rates in these species. Moreover, we find little evidence for convergent amino acid changes across species. Our work not only shows that reproduction-related genes evolve rapidly also in hermaphroditic animals, but also that well-replicated contrasts of different sexual selection contexts can reveal underappreciated genome-wide effects.