A phylogenetically informed search for an alternative Macrostomum model species, with notes on taxonomy, mating behavior, karyology, and genome size

Contrasting the form and strength of pre- and postcopulatory sexual selection in a flatworm

Sexual traits may be selected during multiple consecutive episodes of selection, occurring before, during, or after copulation. The overall strength and form of selection acting on traits may thus be determined by how selection (co-)varies along different episodes. However, it is challenging to measure pre- and postcopulatory phenotypic traits alongside variation in fitness components at each different episode. Here, we used a transgenic line of the transparent flatworm Macrostomum lignano expressing green fluorescent protein (GFP) in all cell types, including sperm cells, enabling in vivo sperm tracking. We assessed the mating success, sperm-transfer efficiency, and sperm fertilizing efficiency of GFP(+) focal worms in which we measured 13 morphological traits. We found linear selection on sperm production rate arising from pre- and postcopulatory components and on copulatory organ shape arising from sperm fertilizing efficiency. We further found nonlinear (mostly concave) selection on combinations of copulatory organ and sperm morphology traits arising mostly from sperm-transfer efficiency and sperm fertilizing efficiency. Our study provides a fine-scale quantification of sexual selection, showing that both the form and strength of selection can change across fitness components. Quantifying how sexual selection builds up along episodes of selection allows us to better understand the evolution of sexually selected traits.

Unraveling the Unusual Subgenomic Organization in the Neopolyploid Free-Living Flatworm Macrostomum lignano

Whole genome duplication (WGD) is an evolutionary event resulting in a redundancy of genetic material. Different mechanisms of WGD, allo- or autopolyploidization, lead to distinct evolutionary trajectories of newly formed polyploids. Genome studies on such species are important for understanding the early stages of genome evolution. However, assembling neopolyploid is a challenging task due to the presence of 2 homologous (or homeologous) chromosome sets and therefore the existence of the extended paralogous regions in its genome. Post-WGD evolution of polyploids includes cytogenetic diploidization leading to the formation of species, whose polyploid origin might be hidden by disomic inheritance. Earlier we uncovered the hidden polyploid origin of the free-living flatworms of the genus Macrostomum (Macrostomum lignano, M. janickei, and M. mirumnovem). Cytogenetic diploidization in these species is accompanied by intensive chromosomal rearrangements including chromosomes fusions. In this study, we unravel the M. lignano genome organization through generation and sequencing of 2 sublines of the commonly used inbred line of M. lignano (called DV1) differing only in a copy number of the largest chromosome (MLI1). Using nontrivial assembly free comparative analysis of their genomes, we deciphered DNA sequences belonging to MLI1 and validated them by sequencing the pool of microdissected MLI1. Here we presented the uncommon mechanism of genome rediplodization of M. lignano, which consists of (i) presence of 3 subgenomes, which emerged via formation of large fused chromosomes and its variants, and (ii) sustaining their heterozygosity through inter- and intrachromosomal rearrangements.

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.

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.

Three new species of Macrostomum (Platyhelminthes, Macrostomorpha) from China and Australia, with notes on taxonomy and phylogenetics

In this paper, three species of the macrostomid free-living flatworm genus Macrostomum are described. Two species, Macrostomumlittorale Wang & Shi, sp. nov. and M.shekouense Wang & Shi, sp. nov., were collected from coastal water at Shenzhen, Guangdong Province, China. One species, M.brandi Wang & Shi, sp. nov., was collected from Perth, Western Australia and Queenscliff, Victoria, Australia. Macrostomumlittorale sp. nov. differs from congeneric species within the genus in the length of the stylet, diameter of stylet opening, and the bend of the stylet. Macrostomumshekouense sp. nov. and M.brandi sp. nov. differ from similar species within the genus in the stylet morphology, position of the female antrum and developing eggs, or presence or absence of the false seminal vesicle. Phylogenetic analysis based on cytochrome c oxidase subunit I (COI) gene shows that M.littorale sp. nov. and M.hystrix are sister clades on two well-separated branch, M.shekouense sp. nov. and M.brandi sp. nov. are sister clades on two well-separated branches. Accordingly, both morphological and molecular evidence support M.littorale sp. nov., M.shekouense sp. nov., and M.brandi sp. nov. as three new species.

Mating behavior and reproductive morphology predict macroevolution of sex allocation in hermaphroditic flatworms

Background

Sex allocation is the distribution of resources to male or female reproduction. In hermaphrodites, this concerns an individual’s resource allocation to, for example, the production of male or female gametes. Macroevolutionary studies across hermaphroditic plants have revealed that the self-pollination rate and the pollination mode are strong predictors of sex allocation. Consequently, we expect similar factors such as the selfing rate and aspects of the reproductive biology, like the mating behaviour and the intensity of postcopulatory sexual selection, to predict sex allocation in hermaphroditic animals. However, comparative work on hermaphroditic animals is limited. Here, we study sex allocation in 120 species of the hermaphroditic free-living flatworm genus Macrostomum. We ask how hypodermic insemination, a convergently evolved mating behaviour where sperm are traumatically injected through the partner’s epidermis, affects the evolution of sex allocation. We also test the commonly-made assumption that investment into male and female reproduction should trade-off. Finally, we ask if morphological indicators of the intensity of postcopulatory sexual selection (female genital complexity, male copulatory organ length, and sperm length) can predict sex allocation.

Results

We find that the repeated evolution of hypodermic insemination predicts a more female-biased sex allocation (i.e., a relative shift towards female allocation). Moreover, transcriptome-based estimates of heterozygosity reveal reduced heterozygosity in hypodermically mating species, indicating that this mating behavior is linked to increased selfing or biparental inbreeding. Therefore, hypodermic insemination could represent a selfing syndrome. Furthermore, across the genus, allocation to male and female gametes is negatively related, and larger species have a more female-biased sex allocation. Finally, increased female genital complexity, longer sperm, and a longer male copulatory organ predict a more male-biased sex allocation.

Conclusions

Selfing syndromes have repeatedly originated in plants. Remarkably, this macroevolutionary pattern is replicated in Macrostomum flatworms and linked to repeated shifts in reproductive behavior. We also find a trade-off between male and female reproduction, a fundamental assumption of most theories of sex allocation. Beyond that, no theory predicts a more female-biased allocation in larger species, suggesting avenues for future work. Finally, morphological indicators of more intense postcopulatory sexual selection appear to predict more intense sperm competition.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01234-1.

Frequent origins of traumatic insemination involve convergent shifts in sperm and genital morphology

Traumatic insemination is a mating behavior during which the (sperm) donor uses a traumatic intromittent organ to inject an ejaculate through the epidermis of the (sperm) recipient, thereby frequently circumventing the female genitalia. Traumatic insemination occurs widely across animals, but the frequency of its evolution, the intermediate stages via which it originates, and the morphological changes that such shifts involve remain poorly understood. Based on observations in 145 species of the free-living flatworm genus Macrostomum, we identify at least nine independent evolutionary origins of traumatic insemination from reciprocal copulation, but no clear indication of reversals. These origins involve convergent shifts in multivariate morphospace of male and female reproductive traits, suggesting that traumatic insemination has a canalizing effect on morphology. We also observed sperm in both the sperm receiving organ and within the body tissue of two species. These species had intermediate trait values indicating that traumatic insemination evolves through initial internal wounding during copulation. Finally, signatures of male-female coevolution of genitalia across the genus indicate that sexual selection and sexual conflict drive the evolution of traumatic insemination, because it allows donors to bypass postcopulatory control mechanisms of recipients.

B Chromosomes in Free-Living Flatworms of the Genus Macrostomum (Platyhelminthes, Macrostomorpha)

B chromosomes (Bs) or supernumerary chromosomes are extra chromosomes in the species karyotype that can vary in its copy number. Bs are widespread in eukaryotes. Usually, the Bs of specimens collected from natural populations are the object of the B chromosome studies. We applied another approach analyzing the Bs in animals maintained under the laboratory conditions as lines and cultures. In this study, three species of the Macrostomum genus that underwent a recent whole-genome duplication (WGD) were involved. In laboratory lines of M. lignano and M. janickei, the frequency of Bs was less than 1%, while in the laboratory culture of M. mirumnovem, it was nearer 30%. Their number in specimens of the culture varied from 1 to 14. Mosaicism on Bs was discovered in parts of these animals. We analyzed the distribution of Bs among the worms of the laboratory cultures during long-term cultivation, the transmission rates of Bs in the progeny obtained from crosses of worms with different numbers of Bs, and from self-fertilized isolated worms. The DNA content of the Bs in M. mirumnovem was analyzed with the chromosomal in situ suppression (CISS) hybridization of microdissected DNA probes derived from A chromosomes (As). Bs mainly consisted of repetitive DNA. The cytogenetic analysis also revealed the divergence and high variation in large metacentric chromosomes (LMs) containing numerous regions enriched for repeats. The possible mechanisms of the appearance and evolution of Bs and LMs in species of the Macrostomum genus were also discussed.

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.

(Un)expected Similarity of the Temporary Adhesive Systems of Marine, Brackish, and Freshwater Flatworms

Many free-living flatworms have evolved a temporary adhesion system, which allows them to quickly attach to and release from diverse substrates. In the marine Macrostomum lignano, the morphology of the adhesive system and the adhesion-related proteins have been characterised. However, little is known about how temporary adhesion is performed in other aquatic environments. Here, we performed a 3D reconstruction of the M. lignano adhesive organ and compared it to the morphology of five selected Macrostomum, representing two marine, one brackish, and two freshwater species. We compared the protein domains of the two adhesive proteins, as well as an anchor cell-specific intermediate filament. We analysed the gene expression of these proteins by in situ hybridisation and performed functional knockdowns with RNA interference. Remarkably, there are almost no differences in terms of morphology, protein regions, and gene expression based on marine, brackish, and freshwater habitats. This implies that glue components produced by macrostomids are conserved among species, and this set of two-component glue functions from low to high salinity. These findings could contribute to the development of novel reversible biomimetic glues that work in all wet environments and could have applications in drug delivery systems, tissue adhesives, or wound dressings.

Large-scale phylogenomics of the genus Macrostomum (Platyhelminthes) reveals cryptic diversity and novel sexual traits

Free-living flatworms of the genus Macrostomum are small and transparent animals, representing attractive study organisms for a broad range of topics in evolutionary, developmental, and molecular biology. The genus includes the model organism M. lignano for which extensive molecular resources are available, and recently there is a growing interest in extending work to additional species in the genus. These endeavours are currently hindered because, even though >200 Macrostomum species have been taxonomically described, molecular phylogenetic information and geographic sampling remain limited. We report on a global sampling campaign aimed at increasing taxon sampling and geographic representation of the genus. Specifically, we use extensive transcriptome and single-locus data to generate phylogenomic hypotheses including 145 species. Across different phylogenetic methods and alignments used, we identify several consistent clades, while their exact grouping is less clear, possibly due to a radiation early in Macrostomum evolution. Moreover, we uncover a large undescribed diversity, with 94 of the studied species likely being new to science, and we identify multiple novel morphological traits. Furthermore, we identify cryptic speciation in a taxonomically challenging assemblage of species, suggesting that the use of molecular markers is a prerequisite for future work, and we describe the distribution of putative synapomorphies and suggest taxonomic revisions based on our finding. Our large-scale phylogenomic dataset now provides a robust foundation for comparative analyses of morphological, behavioural and molecular evolution in this genus.

Two New Brackish-Water Species of Macrostomum (Platyhelminthes: Macrostomorpha) from China and Their Phylogenetic Positions

In this paper, two new brackish-water species of the macrostomid turbellarian genus Macrostomum, Macrostomum pseudosinense sp. nov. and Macrostomum taurinum sp. nov., collected from coastal water at Shenzhen, Guangdong Province, China, are described based on morphological, histological, and molecular phylogenetic analyses. Macrostomum pseudosinense sp. nov. differs from similar species within the genus in the length of the stylet (152 ± 15.0 µm), diameter of stylet opening (20 ± 4.0 µm proximally; 7 ± 0.5 µm distally), two bends of the stylet, and the non-spiral end of the stylet. Macrostomum taurinum sp. nov. differs from its congeners in the length of the stylet (81 ± 7.4 µm), the stylet bending position and angle (50% and 60°), diameter of stylet proximal opening (15 ± 3.0 µm), sperm with bristles and brush, and the smooth-walled ovaries. Phylogenetic analyses inferred from nuclear 18S and 28S rRNA genes support the establishments of these two new species. In addition, reciprocal mating behavior of M. pseudosinense sp. nov. was observed and documented.

The repeatable opportunity for selection differs between pre- and postcopulatory fitness components

In species with multiple mating, intense sexual selection may occur both before and after copulation. However, comparing the strength of pre- and postcopulatory selection is challenging, because (i) postcopulatory processes are generally difficult to observe and (ii) the often-used opportunity for selection (I) metric contains both deterministic and stochastic components. Here, we quantified pre- and postcopulatory male fitness components of the simultaneously hermaphroditic flatworm, Macrostomum lignano. We did this by tracking fluorescent sperm-using transgenics-through the transparent body of sperm recipients, enabling to observe postcopulatory processes in vivo. Moreover, we sequentially exposed focal worms to three independent mating groups, and in each assessed their mating success, sperm-transfer efficiency, sperm fertilizing efficiency, and partner fecundity. Based on these multiple measures, we could, for each fitness component, combine the variance (I) with the repeatability (R) in individual success to assess the amount of repeatable variance in individual success-a measure we call the repeatable opportunity for selection (IR ). We found higher repeatable opportunity for selection in sperm-transfer efficiency and sperm fertilizing efficiency compared to mating success, which clearly suggests that postcopulatory selection is stronger than precopulatory selection. Our study demonstrates that the opportunity for selection contains a repeatable deterministic component, which can be assessed and disentangled from the often large stochastic component, to provide a better estimate of the strength of selection.

Regions enriched for DNA repeats in chromosomes of Macrostomum mirumnovem, a species with a recent Whole Genome Duplication

The free-living flatworm Macrostomum mirumnovem is a neopolyploid species whose genome underwent a recent Whole Genome Duplication (WGD). In the result of chromosome fusions of the ancient haploid chromosome set, large metacentric chromosomes were formed. In addition to three pairs of small metacentrics, the current karyotype of M. mirumnovem contains two pairs of large metacentric chromosomes, MMI1 and MMI2. The generation of microdissected DNA libraries enriched for DNA repeats followed by DNA probe preparation and fluorescent in situ hybridization (FISH) were performed. The DNA probes obtained marked chromosome regions enriched for different DNA repeats in the M. mirumnovem chromosomes. The size and localization of these regions varied in different copies of large chromosomes. They varied even in homologous chromosomes, suggesting their divergence due to genome re-diploidization after a WGD. Besides the newly formed chromosome regions enriched for DNA repeats, B chromosomes were found in the karyotypes of the studied specimens of M. mirumnovem. These B chromosomes varied in size and morphology. FISH with microdissected DNA probes revealed that some Bs had a distinct DNA content. FISH could paint differently B chromosomes in different worms and even in the same sample. B chromosomes could carry a bright specific fluorescent signal or could show no fluorescent signal at all. In latter cases, the specific FISH signal could be absent even in the pericentromeric region of the B chromosome. Possible mechanisms of B chromosome formation and their further evolution are discussed. The results obtained indicate an important role that repetitive DNAs play in genome re-diploidization initiating a rapid differentiation of large chromosome copies. Taking together, karyotype peculiarities (a high level of intraspecific karyotypic diversity associated with chromosome number variation, structural chromosomal rearrangements, and the formation of new regions enriched for DNA repeats) and some phenotypic features of M. mirumnovem (small body size, short lifecycle, easy maintenance in the laboratory) make this species a perspective model in the studies of genomic and karyotypic evolution in species passed through a recent WGD event.

Generation of microdissected DNA probes from metaphase chromosomes when chromosome identification by routine staining is impossible

Application of microdissected DNA libraries and DNA probes in numerous and various modern molecular cytogenetic studies showed them as an efficient and reliable tool in the analysis of chromosome reorganization during karyotypic evolution and in the diagnosis of human chromosome pathology. An important advantage of DNA probe generation by metaphase chromosome microdissection followed by sequence-independent polymerase chain reaction in comparison with the method of DNA probe generation using chromosome sorting is the possibility of DNA probe preparation from chromosomes of an individual sample without cell line establishment for the production of a large number of metaphase chromosomes. One of the main requirements for successful application of this technique is a possibility for identification of the chromosome of interest during its dissection and collection of its material from metaphase plates spread on the coverslip. In the present study, we developed and applied a technique for generation of microdissected DNA probes in the case when chromosome identification during microdissection appeared to be impossible. The technique was used for generation of two sets of Whole Chromosome Paints (WCPs) from all chromosomes of two species of free-living flatworms in the genus Macrostomum, M. mirumnovem and M. cliftonensis. The single-copy chromosome technique including separate collection of all chromosomes from one metaphase plate allowed us to generate WCPs that painted specifically the original chromosome by Chromosome In Situ Suppression Hybridization (CISS-Hybridization). CISS-Hybridization allowed identifying the original chromosome(s) used for DNA probe generation. Pooled WCPs derived from homologous chromosomes increased the intensity and specificity of chromosome painting provided by CISS-Hybridization. In the result, the obtained DNA probes appeared to be good enough for application in our studies devoted to analysis of karyotypic evolution in the genus Macrostomum and for analysis of chromosome rearrangements among the worms of laboratory cultures of M. mirumnovem.

Successful mating and hybridisation in two closely related flatworm species despite significant differences in reproductive morphology and behaviour

Reproductive traits are some of the fastest diverging characters and can serve as reproductive barriers. The free-living flatworm Macrostomum lignano, and its congener M. janickei are closely related, but differ substantially in their male intromittent organ (stylet) morphology. Here, we examine whether these morphological differences are accompanied by differences in behavioural traits, and whether these could represent barriers to successful mating and hybridization between the two species. Our data shows that the two species differ in many aspects of their mating behaviour. Despite these differences, the species mate readily with each other in heterospecific pairings. Although both species have similar fecundity in conspecific pairings, the heterospecific pairings revealed clear postmating barriers, as few heterospecific pairings produced F1 hybrids. These hybrids had a stylet morphology that was intermediate between that of the parental species, and they were fertile. Finally, using a mate choice experiment, we show that the nearly two-fold higher mating rate of M. lignano caused it to mate more with conspecifics, leading to assortative mating, while M. janickei ended up mating more with heterospecifics. Thus, while the two species can hybridize, the mating rate differences could possibly lead to higher fitness costs for M. janickei compared to M. lignano.

Genome and Karyotype Reorganization after Whole Genome Duplication in Free-Living Flatworms of the Genus Macrostomum

The genus Macrostomum represents a diverse group of rhabditophoran flatworms with >200 species occurring around the world. Earlier we uncovered karyotype instability linked to hidden polyploidy in both M. lignano (2n = 8) and its sibling species M. janickei (2n = 10), prompting interest in the karyotype organization of close relatives. In this study, we investigated chromosome organization in two recently described and closely related Macrostomum species, M. mirumnovem and M. cliftonensis, and explored karyotype instability in laboratory lines and cultures of M. lignano (DV1/10, 2n = 10) and M. janickei in more detail. We revealed that three of the four studied species are characterized by karyotype instability, while M. cliftonensis showed a stable 2n = 6 karyotype. Next, we performed comparative cytogenetics of these species using fluorescent in situ hybridization (FISH) with a set of DNA probes (including microdissected DNA probes generated from M. lignano chromosomes, rDNA, and telomeric DNA). To explore the chromosome organization of the unusual 2n = 9 karyotype discovered in M. mirumnovem, we then generated chromosome-specific DNA probes for all chromosomes of this species. Similar to M. lignano and M. janickei, our findings suggest that M. mirumnovem arose via whole genome duplication (WGD) followed by considerable chromosome reshuffling. We discuss possible evolutionary scenarios for the emergence and reorganization of the karyotypes of these Macrostomum species and consider their suitability as promising animal models for studying the mechanisms and regularities of karyotype and genome evolution after a recent WGD.

Variation in sex allocation plasticity in three closely related flatworm species

Sex allocation (SA) theory for simultaneous hermaphrodites predicts an influence of group size on SA. Since group size can vary within an individual's lifetime, this can favor the evolution of phenotypically plastic SA. In an emerging comparative context, we here report on SA plasticity in three closely related Macrostomum flatworm species, namely Macrostomum janickei, Macrostomum cliftonensis, and Macrostomum mirumnovem. For each species, we experimentally raised worms in three group sizes (isolated, pairs, and octets) and two enclosure sizes (small and large) in all factorial combinations and studied the effects of these factors on different estimates of SA. In addition, we also evaluated whether isolated worms engage in self-fertilization. We found that all species have plastic SA, with M. cliftonensis being more plastic than the other two species, as assessed by comparing standardized effect sizes of (a) the presence/absence of mating partners and (b) the strength of sexual competition. Moreover, we found that sperm production rate-but not sperm morphology-is plastic in M. cliftonensis, and that only M. mirumnovem self-fertilized during our observation period. Our study suggests that both SA and SA plasticity can diverge even between closely related species.