Population genetic structure and mating system evolution in freshwater pulmonates

Bulinus snails in the Lake Victoria Basin in Kenya: Systematics and their role as hosts for schistosomes

The planorbid gastropod genus Bulinus consists of 38 species that vary in their ability to vector Schistosoma haematobium (the causative agent of human urogenital schistosomiasis), other Schistosoma species, and non-schistosome trematodes. Relying on sequence-based identifications of bulinids (partial cox1 and 16S) and Schistosoma (cox1 and ITS), we examined Bulinus species in the Lake Victoria Basin in Kenya for naturally acquired infections with Schistosoma species. We collected 6,133 bulinids from 11 sites between 2014-2021, 226 (3.7%) of which harbored Schistosoma infections. We found 4 Bulinus taxa from Lake Victoria (B. truncatus, B. tropicus, B. ugandae, and B. cf. transversalis), and an additional 4 from other habitats (B. globosus, B. productus, B. forskalii, and B. scalaris). S. haematobium infections were found in B. globosus and B. productus (with infections in the former predominating) whereas S. bovis infections were identified in B. globosus, B. productus, B. forskalii, and B. ugandae. No nuclear/mitochondrial discordance potentially indicative of S. haematobium/S. bovis hybridization was detected. We highlight the presence of Bulinus ugandae as a distinct lake-dwelling taxon closely related to B. globosus yet, unlike all other members of the B. africanus species group, is likely not a vector for S. haematobium, though it does exhibit susceptibility to S. bovis. Other lake-dwelling bulinids also lacked S. haematobium infections, supporting the possibility that they all lack compatibility with local S. haematobium, thereby preventing widespread transmission of urogenital schistosomiasis in the lake's waters. We support B. productus as a distinct species from B. nasutus, B. scalaris as distinct from B. forskalii, and add further evidence for a B. globosus species complex with three lineages represented in Kenya alone. This study serves as an essential prelude for investigating why these patterns in compatibility exist and whether the underlying biological mechanisms may be exploited for the purpose of limiting schistosome transmission.

House mouse subspecies do differ in their social structure

It is widely acknowledged that population structure can have a substantial impact on evolutionary trajectories. In social animals, this structure is strongly influenced by relationships among the population members, so studies of differences in social structure between diverging populations or nascent species are of prime interest. Ideal models for such a study are two house mouse subspecies, Mus musculus musculus and M. m. domesticus, meeting in Europe along a secondary contact zone. Though the latter subspecies has usually been supposed to form tighter and more isolated social units than the former, the evidence is still inconclusive. Here, we carried out a series of radiofrequency identification experiments in semi-natural enclosures to gather large longitudinal data sets on individual mouse movements. The data were summarized in the form of uni- and multi-layer social networks. Within them, we could delimit and describe the social units ("modules"). While the number of estimated units was similar in both subspecies, domesticus revealed a more "modular" structure. This subspecies also showed more intramodular social interactions, higher spatial module separation, higher intramodular persistence of parent-offspring contacts, and lower multiple paternity, suggesting more effective control of dominant males over reproduction. We also demonstrate that long-lasting modules can be identified with basic reproductive units or demes. We thus provide the first robust evidence that the two subspecies differ in their social structure and dynamics of the structure formation.

Effects of the anxiolytic benzodiazepine oxazepam on freshwater gastropod reproduction: a prospective study

Psychoactive drugs have emerged as contaminants over the last few decades. These drugs are frequently prescribed and poorly eliminated by wastewater treatment plants, and many are present at non-negligible concentrations in surface waters. Several studies have investigated the non-target organism toxicity of one such drug, oxazepam, a benzodiazepine anxiolytic frequently detected in rivers. However, very little is known about the impact of this drug on reproduction. We investigated the effects of environmentally relevant concentrations of oxazepam on Radix balthica, a freshwater gastropod widespread in Europe. We identified the reproductive organs of Radix balthica. We then exposed this gastropod to oxazepam for two months and assessed several reproductive parameters, from reproductive organ status to behavioral parameters. We found that adults exposed to 10 µg/L oxazepam display an increase in the density of spermatozoa, and that adults exposed to 0.8 µg/L oxazepam displayed a decrease in the number of eggs per egg mass over time. By contrast, oxazepam had no effect on shell length, the size of male reproductive organs or social interactions. Finally, a locomotor activity analysis showed the distance covered over time decreased in all conditions of exposure to oxazepam, potentially reflecting a disturbance of exploratory activity. These results shed light on the effects of oxazepam on the reproduction of a non-target freshwater mollusk.

Reintroduction of freshwater macroinvertebrates: challenges and opportunities

Species reintroductions - the translocation of individuals to areas in which a species has been extirpated with the aim of re-establishing a self-sustaining population - have become a widespread practice in conservation biology. Reintroduction projects have tended to focus on terrestrial vertebrates and, to a lesser extent, fishes. Much less effort has been devoted to the reintroduction of invertebrates into restored freshwater habitats. Yet, reintroductions may improve restoration outcomes in regions where impoverished regional species pools limit the self-recolonisation of restored freshwaters. We review the available literature on macroinvertebrate reintroductions, focusing on identifying the intrinsic and extrinsic factors that determine their success or failure. Our study reveals that freshwater macroinvertebrate reintroductions remain rare, are often published in the grey literature and, of the attempts made, approximately one-third fail. We identify life-cycle complexity and remaining stressors as the two factors most likely to affect reintroduction success, illustrating the unique challenges of freshwater macroinvertebrate reintroductions. Consideration of these factors by managers during the planning process and proper documentation - even if a project fails - may increase the likelihood of successful outcomes in future reintroduction attempts of freshwater macroinvertebrates.

Mixed breeding system in the hermaphroditic land slug Arion intermedius (Stylommatophora, Arionidae)

Theory suggests that hermaphroditic plants and animals should be either entirely outcrossing or entirely selfing. As such, very few hermaphroditic plants and basommatophoran snails have a mixed breeding system. However, reliable estimates of selfing rates are lacking for most hermaphroditic animals. This partly prevents to delineate the relative contributions of the selective factors that determine selfing and outcrossing rates in hermaphroditic animal taxa. Here, we studied the population genetic structure of, and breeding system in, 11 populations of the hermaphroditic land slug Arion intermedius using five polymorphic microsatellite loci. Moreover, genotype frequencies deviated significantly from Hardy-Weinberg equilibrium expectations for most of the loci in all populations suggesting some level of selfing. Estimates of the selfing level s, suggest moderate levels of outcrossing (mean s based on FIS = 0.84; mean s based on the two-locus heterozygosity disequilibrium = 0.20, or with a ML approach = 0.22). Our study therefore suggests that A. intermedius has a mixed breeding system. A re-analysis of allozyme data from another arionid slug ( subgenus Carinarion) indicates that mixed breeding may be more common in arionid slugs than hitherto was assumed. These results seem therefore at variance with current theoretical and empirical predictions and opens perspectives for the study on the evolutionary factors driving mixed breeding systems in animals.

Life history traits variation in heterogeneous environment: The case of a freshwater snail resistance to pond drying

Ecologists and population geneticists have long suspected that the diversity of living organisms was connected to the structure of their environment. In heterogeneous environments, diversifying selection combined to restricted gene flow may indeed lead to locally adapted populations. The freshwater snail, Galba truncatula, is a good model to address this question because it is present in a heterogeneous environment composed of temporary and permanent waters. In order to test the selective importance of those environments, we proposed here to measure survival of lineages from both habitats during drought episodes. To this purpose, we experimentally submitted adults and juveniles individuals from both habitats to drought. We found a difference in desiccation resistance between temporary and permanents waters only for adults. Adults from temporary habitats were found more resistant to drought. This divergence in desiccation resistance seems to explain the unexpected life history traits differences between habitats observed.

Factors and processes shaping the population structure and distribution of genetic variation across the species range of the freshwater snail radix balthica (Pulmonata, Basommatophora)

BACKGROUND

Factors and processes shaping the population structure and spatial distribution of genetic diversity across a species' distribution range are important in determining the range limits. We comprehensively analysed the influence of recurrent and historic factors and processes on the population genetic structure, mating system and the distribution of genetic variability of the pulmonate freshwater snail Radix balthica. This analysis was based on microsatellite variation and mitochondrial haplotypes using Generalised Linear Statistical Modelling in a Model Selection framework.

RESULTS

Populations of R. balthica were found throughout North-Western Europe with range margins marked either by dispersal barriers or the presence of other Radix taxa. Overall, the population structure was characterised by distance independent passive dispersal mainly along a Southwest-Northeast axis, the absence of isolation-by-distance together with rather isolated and genetically depauperated populations compared to the variation present in the entire species due to strong local drift. A recent, climate driven range expansion explained most of the variance in genetic variation, reducing at least temporarily the genetic variability in this area. Other factors such as geographic marginality and dispersal barriers play only a minor role.

CONCLUSIONS

To our knowledge, such a population structure has rarely been reported before. It might nevertheless be typical for passively dispersed, patchily distributed taxa (e.g. freshwater invertebrates). The strong local drift implied in such a structure is expected to erode genetic variation at both neutral and coding loci and thus probably diminish evolutionary potential. This study shows that the analysis of multiple factors is crucial for the inference of the processes shaping the distribution of genetic variation throughout species ranges.

Genetic diversity and structure in two species of Leavenworthia with self-incompatible and self-compatible populations

Self-fertilization is a common mating system in plants and is known to reduce genetic diversity, increase genetic structure and potentially put populations at greater risk of extinction. In this study, we measured the genetic diversity and structure of two cedar glade endemic species, Leavenworthia alabamica and L. crassa. These species have self-incompatible (SI) and self-compatible (SC) populations and are therefore ideal for understanding how the mating system affects genetic diversity and structure. We found that L. alabamica and L. crassa had high species-level genetic diversity (H(e)=0.229 and 0.183, respectively) and high genetic structure among their populations (F(ST)=0.45 and 0.36, respectively), but that mean genetic diversity was significantly lower in SC compared with SI populations (SC vs SI, H(e) for L. alabamica was 0.065 vs 0.206 and for L. crassa was 0.084 vs 0.189). We also found significant genetic structure using maximum-likelihood clustering methods. These data indicate that the loss of SI leads to the loss of genetic diversity within populations. In addition, we examined genetic distance relationships between SI and SC populations to analyze possible population history and origins of self-compatibility. We find there may have been multiple origins of self-compatibility in L. alabamica and L. crassa. However, further work is required to test this hypothesis. Finally, given their high genetic structure and that individual populations harbor unique alleles, conservation strategies seeking to maximize species-level genetic diversity for these or similar species should protect multiple populations.

Low genetic diversity in a snail intermediate host (Biomphalaria pfeifferi Krass, 1848) and schistosomiasis transmission in the Senegal River Basin

Population genetic perturbations of intermediate hosts, often a consequence of human pressure on environmental resources, can precipitate unexpectedly severe disease outbreaks. Such disturbances are set to become increasingly common following range changes concomitant with climate shifts, dwindling natural resources and major infrastructure changes such as hydroprojects. Construction of the Diama dam in the Senegal River Basin (SRB) reduced river salinity, enabling the freshwater snail intermediate host Biomphalaria pfeifferi to rapidly expand its distribution. A serious public health problem ensued, with an epidemic of intestinal schistosomiasis occurring in the previously schistosome-free Richard-Toll region within 2 years. The current study aimed to assess the population variability of B. pfeifferi in the SRB, and speculate upon its subsequent impact on host-parasite interactions following such engineered ecological change. Genetic variation at nine polymorphic microsatellite loci revealed little population differentiation in SRB snails compared with those from natural habitats in Zimbabwe, where Schistosoma mansoni transmission is much lower. 'Open' SRB habitats are associated with greater water contact, smaller population sizes and less genetic diversity, with sites downstream of Richard-Toll showing greater inter- and intrapopulation variation, concomitant with less frequent human contact. These observations may be explained by rapid expansion into pristine habitat selecting for high fecundity genotypes at the expense of schistosome resistance, presenting S. mansoni with genetically homogenous highly fecund susceptible populations around the focal point, promoting development of a highly compatible host-parasite relationship. Longitudinal study of such systems may prove important in predicting public health risks engendered by future environmental engineering projects.

Highly polymorphic microsatellite markers for Radix balthica (Linnaeus 1758)

We present data for eight polymorphic microsatellite markers isolated from a microsatellite-enriched DNA library for the freshwater snail Radix balthica. Three of them were specific for R. balthica while five also amplified polymorphic products in two congeneric species. Test application on populations from all over the species range has shown that these loci are highly informative for analysing population structure and estimating migration rates. Observed deviations from Hardy-Weinberg equilibrium are attributed to a mixed mating system.

The causes of mutation accumulation in mitochondrial genomes

A fundamental observation across eukaryotic taxa is that mitochondrial genomes have a higher load of deleterious mutations than nuclear genomes. Identifying the evolutionary forces that drive this difference is important to understanding the rates and patterns of sequence evolution, the efficacy of natural selection, the maintenance of sex and recombination and the mechanisms underlying human ageing and many diseases. Recent studies have implicated the presumed asexuality of mitochondrial genomes as responsible for their high mutational load. We review the current body of knowledge on mitochondrial mutation accumulation and recombination, and conclude that asexuality, per se, may not be the primary determinant of the high mutation load in mitochondrial DNA (mtDNA). Very little recombination is required to counter mutation accumulation, and recent evidence suggests that mitochondrial genomes do experience occasional recombination. Instead, a high rate of accumulation of mildly deleterious mutations in mtDNA may result from the small effective population size associated with effectively haploid inheritance. This type of transmission is nearly ubiquitous among mitochondrial genomes. We also describe an experimental framework using variation in mating system between closely related species to disentangle the root causes of mutation accumulation in mitochondrial genomes.

Predominance of outcrossing in Lymnaea stagnalis despite low apparent fitness costs of self-fertilization

We have quantified the natural mating system in eight populations of the simultaneously hermaphroditic aquatic snail Lymnaea stagnalis, and studied the ecological and genetic forces that may be directing mating system evolution in this species. We investigated whether the natural mating system can be explained by the availability of mates, by the differential survival of self- and cross-fertilized snails in nature, and by the effects of mating system on parental fecundity and early survival. The natural mating system of L. stagnalis was found to be predominantly cross-fertilizing. Density of snails in the populations had no relationship with the mating system, suggesting that outcrossing rates are not limited by mate availability at the population densities observed. Contrary to expectations for outcrossing species, we detected no evidence for inbreeding depression in survival in nature with inferential population genetic methods. Further, experimental manipulations of mating system in the laboratory revealed that self-fertilization had no effect on parental fecundity, and only minor effects on offspring survival. Predominance of cross-fertilization despite low apparent fitness costs of self-fertilization is at odds with the paradigm that high self-fertilization depression is necessary for maintenance of cross-fertilization in self-compatible hermaphrodites.

Endocrine disruption in aquatic pulmonate molluscs: few evidences, many challenges

As compared to other groups of aquatic gastropods, documented examples of endocrine disruption in pulmonates are rather limited. This is quite surprising because the endocrine control of physiological functions has been extensively studied in these animals. In the model-species Lymnaea stagnalis, the neurohormonal regulation of reproduction has been thoroughly investigated, and the primary structure of several peptides and receptors involved in endocrine processes has been established. However, the use of this knowledge has been fairly limited in the context of ecotoxicology, to investigate the effects of endocrine-disrupting chemicals. The present review summarizes the main and more recent findings on the neuroendocrine control of reproduction in aquatic pulmonate snails (Basommatophora). It then comprehensively describes selected in vivo laboratory and semi-field studies which provide evidence for possible endocrine disrupting effects of estrogenic and androgenic test compounds [e.g., ethynylestradiol, methyltestosterone (MT)], and of environmental contaminants [e.g., cadmium (Cd), tributyltin (TBT), and nonylphenol (NP), pesticides]. Finally, challenging perspectives for future research are discussed.

Breeding system and demography shape population genetic structure across ecological and climatic zones in the African freshwater snail, Bulinus forskalii (Gastropoda, Pulmonata), intermediate host for schistosomes

The role of breeding system and population bottlenecks in shaping the distribution of neutral genetic variation among populations inhabiting patchily distributed, ephemeral water bodies was examined for the hermaphroditic freshwater snail Bulinus forskalii, intermediate host for the medically important trematode Schistosoma guineensis. Levels of genetic variation at 11 microsatellite loci were assessed for 600 individuals sampled from 19 populations that span three ecological and climatic zones (ecozones) in Cameroon, West Africa. Significant heterozygote deficiencies and linkage disequilibria indicated very high selfing rates in these populations. Despite this and the large genetic differentiation detected between populations, high levels of genetic variation were harboured within these populations. The high level of gene flow inferred from assignment tests may be responsible for this pattern. Indeed, metapopulation dynamics, including high levels of gene flow as well as extinction/contraction and recolonization events, are invoked to account for the observed population structuring, which was not a consequence of isolation-by-distance. Because B. forskalii populations inhabiting the northern, Sahelian area are subject to more pronounced annual cycles of drought and flood than the southern equatorial ones, they were expected to be subject to population bottlenecks of increased frequency and severity and, therefore, show reduced genetic variability and elevated population differentiation. Contrary to predictions, the populations inhabiting the most northerly ecozone exhibited higher genetic diversity and lower genetic differentiation than those in the most southerly one, suggesting that elevated gene flow in this region is counteracting genetic drift.

Effects of inbreeding on the genetic diversity of populations

The study of variability within species is important to all biologists who use genetic markers. Since the discovery of molecular variability among normal individuals, data have been collected from a wide range of organisms, and it is important to understand the major factors affecting diversity levels and patterns. Comparisons of inbreeding and outcrossing populations can contribute to this understanding, and therefore studying plant populations is important, because related species often have different breeding systems. DNA sequence data are now starting to become available from suitable plant and animal populations, to measure and compare variability levels and test predictions.

Levels of DNA polymorphism vary with mating system in the nematode genus caenorhabditis

Self-fertilizing species often harbor less genetic variation than cross-fertilizing species, and at least four different models have been proposed to explain this trend. To investigate further the relationship between mating system and genetic variation, levels of DNA sequence polymorphism were compared among three closely related species in the genus Caenorhabditis: two self-fertilizing species, Caenorhabditis elegans and C. briggsae, and one cross-fertilizing species, C. remanei. As expected, estimates of silent site nucleotide diversity were lower in the two self-fertilizing species. For the mitochondrial genome, diversity in the selfing species averaged 42% of diversity in C. remanei. Interestingly, the reduction in genetic variation was much greater for the nuclear than for the mitochondrial genome. For two nuclear genes, diversity in the selfing species averaged 6 and 13% of diversity in C. remanei. We argue that either population bottlenecks or the repeated action of natural selection, coupled with high levels of selfing, are likely to explain the observed reductions in species-wide genetic diversity.

The effect of isolation on the life-history traits of Pseudosuccinea columella (Pulmonata: Lymnaeidae)

A population of Pseudosuccinea columella was raised under laboratory conditions and its life tables were determined in isolated and paired snails. Isolated snails were significantly larger in shell size than paired snails from five weeks of age onward. Also, statistically significant differences were found for the number of eggs per mass per individual from week 5 to 9, isolated snails exhibiting the highest values. The intrinsic and finite rates of increase were greater in isolated than in paired snails. Either an inhibition of the reproductive output between individuals or the advantage of selfing may be the cause of the differences in this species, acting as a possible mechanism that increase the fitness of isolated snails.

Is there a geographical pattern in the breeding system of a complex of hermaphroditic slugs (Mollusca: Gastropoda: Carinarion)?

Allozyme analyses of the hermaphroditic slugs Arion (Carinarion) fasciatus, A. (C.) circumscriptus and A. (C.) silvaticus have suggested that the three species in North America and north-west Europe predominantly reproduce uniparentally, most probably by selfing. We used allozyme electrophoresis to investigate the population genetic structure of these species throughout a larger part of their native European distribution. Our results show that the previously reported "species" specific allozyme markers are no longer valid if populations from central Europe are investigated, and A. fasciatus and A. silvaticus appear to be "paraphyletic" taxa. In contrast to the general belief that selfing organisms show low gene diversities, the high selfing rates in N-NE European Carinarion do not necessarily result in low gene diversities. Moreover, our data suggest a geographical pattern in the prevalence of outcrossing, at least in A. fasciatus, with selfing in N-NE Europe and a mixed breeding system (i.e. selfing and outcrossing) in central Europe. Possible scenarios for the disjunct distribution of breeding systems in Carinarion are discussed.

Allozyme diversity in Leavenworthia populations with different inbreeding levels

Levels of allozyme diversity were studied in five species of the annual plant genus Leavenworthia. The breeding systems of the populations were also characterized in terms of measures of autogamous seed-set and self-fertility. The populations that appeared more inbreeding by these criteria also tended to have lower within-population gene diversity (Hs) values, but there was nevertheless considerable allozyme variability in two of the four sets of inbred populations studied. Comparing the outcrosser L. stylosa and the related inbreeders L. uniflora and L. torulosa the reduction in diversity associated with inbreeding was very high, as no variants were seen within populations of the inbreeders. In L. crassa and L. alabamica, the reduction was lower, but correcting for the fact that this comparison is over less extreme selfing rate differences, it is estimated that the difference between extreme inbreeding and outcrossing populations of L. crassa considerably exceeds twofold. These results are compared with published studies within closely related taxa with different inbreeding levels. The effect of inbreeding on allozyme diversity is consistently larger than has been estimated from comparisons of unrelated species.

Genetic diversity in Leavenworthia populations with different inbreeding levels

Levels of neutral genetic diversity within and between populations were compared between outcrossing (self-incompatible) and inbreeding populations in the annual plant genus Leavenworthia. Two taxonomically independent comparisons are possible, since self-incompatibility has been lost twice in the group of species studied. Within inbred populations of L.uniflora and L.crassa, no DNA sequence variants were seen among the alleles sampled, but high diversity was seen in alleles from populations of the outcrosser L. stylosa, and in self-incompatible L. crassa populations. Diversity between populations was seen in all species. Although total diversity values were lower in the sets of inbreeding populations, between-population values were as high or higher, than those in the outcrossing taxa. Possible reasons for these diversity patterns are discussed. As the effect of inbreeding appears to be a greater than twofold reduction in diversity, we argue that some process such as selection for advantageous mutations, or against deleterious mutations, or bottlenecks occurring predominantly in the inbreeders, appears necessary to account for the findings. If selection for advantageous mutations is responsible, it appears that it must be some form of local adaptive selection, rather than substitution of alleles that are advantageous throughout the species. This is consistent with the finding of high between-population diversity in the inbreeding taxa.

Phylogeography and genetic divergence of some lymnaeid snails, intermediate hosts of human and animal fascioliasis with special reference to lymnaeids from the Bolivian Altiplano

A population genetic study using starch gel electrophoresis was performed on populations of several species of lymnaeid snails acting as intermediate hosts for Fasciola hepatica (Trematoda, Plathyhelminth). Lymnaea viatrix was collected in 16 sites from the Bolivian Northern Altiplano. L. cubensis were obtained in one site from Venezuela, one site from Guadeloupe, three sites from Cuba and one site from the Dominican Republic. L. truncatula were collected in one site from France, one from Portugal and one from Morocco. Multilocus enzyme electrophoresis (MEE) were determined for 282 snails at 18 loci. A complete monomorphism was encountered at each geographic site. However, among these 18 loci, 13 are polymorphic and low and high levels of genetic divergence were observed between samples. Two genotypic groups can be differentiated by their multilocus genotypes. The western genotypic group associates together samples from Venezuela, Guadeloupe, Cuba and Dominican Republic (L. cubensis) while samples from France, Portugal and Morocco (L. truncatula) belong to the eastern genotypic group. Surprisingly, the Northern Bolivian Altiplano populations (L. viatrix) do not present any genetic divergence with the Portuguese sample. Therefore, the Bolivian snails belong entirely to the eastern genetic group. Within each group slight genetic divergences were observed. These results strongly support the European origin of the lymnaeid snails from the Northern Bolivian Altiplano.

Outcrossing rates and correlated matings in a predominantly selfing freshwater snail

We estimated outcrossing rates and correlated matings in a natural population of the hermaphroditic freshwater snail Ancylus fluviatilis, using progeny arrays and polymorphic allozyme markers. This design circumvents drawbacks of earlier studies, and yields maximum likelihood estimates of outcrossing rates at both the individual and population level. Adult snails were sampled during the spring breeding period, allowing for prior copulations and sperm transfer under natural conditions. Once in the laboratory, they were kept isolated for 21 days and laid egg capsules, and emerging hatchlings from 42 families were scored for five polymorphic loci. 87 of 848 offspring exhibited non-maternal alleles, representing a minimum population-level outcrossing of 10.3%. Maximum likelihood estimates are in the range 13-15%. Importantly, there is significant among-family heterogeneity, as only 16 of 42 families are estimated to have outcrossed, in accord with high estimates of correlated selfing within progeny arrays. Moreover, the high proportion of full sibs among outcrossed sibs suggests a limited extent of multiple paternity. Our data reveal important heterogeneity in the mating system of an animal hermaphrodite, and indicate the potential for evolutionary change in the breeding system.