Sex allocation plasticity on a transcriptome scale: Socially sensitive gene expression in a simultaneous hermaphrodite

Transcriptomic changes in the posterior pallium of male zebra finches associated with social niche conformance

Animals plastically adjust their physiological and behavioural phenotypes to conform to their social environment-social niche conformance. The degree of sexual competition is a critical part of the social environment to which animals adjust their phenotypes, but the underlying genetic mechanisms are poorly understood. We conducted a study to investigate how differences in sperm competition risk affect the gene expression profiles of the testes and two brain areas (posterior pallium and optic tectum) in breeding male zebra finches (Taeniopygia castanotis). In this pre-registered study, we investigated a large sample of 59 individual transcriptomes. We compared two experimental groups: males held in single breeding pairs (low sexual competition) versus those held in two pairs (elevated sexual competition) per breeding cage. Using weighted gene co-expression network analysis (WGCNA), we observed significant effects of the social treatment in all three tissues. However, only the treatment effects found in the pallium were confirmed by an additional randomisation test for statistical robustness. Likewise, the differential gene expression analysis revealed treatment effects only in the posterior pallium (ten genes) and optic tectum (six genes). No treatment effects were found in the testis at the single gene level. Thus, our experiments do not provide strong evidence for transcriptomic adjustment specific to manipulated sperm competition risk. However, we did observe transcriptomic adjustments to the manipulated social environment in the posterior pallium. These effects were polygenic rather than based on few individual genes with strong effects. Our findings are discussed in relation to an accompanying paper using the same animals, which reports behavioural results consistent with the results presented here.

Optimized protocols for RNA interference in Macrostomum lignano

Macrostomum lignano, a marine free-living flatworm, has emerged as a potent invertebrate model in developmental biology for studying stem cells, germline, and regeneration processes. In recent years, many tools have been developed to manipulate this worm and to facilitate genetic modification. RNA interference is currently the most accessible and direct technique to investigate gene functions. It is obtained by soaking worms in artificial seawater containing dsRNA targeting the gene of interest. Although easy to perform, the original protocol calls for daily exchange of dsRNA solutions, usually until phenotypes are observed, which is both time- and cost-consuming. In this work, we have evaluated alternative dsRNA delivery techniques, such as electroporation and osmotic shock, to facilitate the experiments with improved time and cost efficiency. During our investigation to optimize RNAi, we demonstrated that, in the absence of diatoms, regular single soaking in artificial seawater containing dsRNA directly produced in bacteria or synthesized in vitro is, in most cases, sufficient to induce a potent gene knockdown for several days with a single soaking step. Therefore, this new and highly simplified method allows a very significant reduction of dsRNA consumption and lab work. In addition, it enables performing experiments on a larger number of worms at minimal cost.

Sex-limited experimental evolution drives transcriptomic divergence in a hermaphrodite

The evolution of gonochorism from hermaphroditism is linked with the formation of sex chromosomes, as well as the evolution of sex-biased and sex-specific gene expression to allow both sexes to reach their fitness optimum. There is evidence that sexual selection drives the evolution of male-biased gene expression in particular. However, previous research in this area in animals comes from either theoretical models or comparative studies of already old sex chromosomes. We therefore investigated changes in gene expression under 3 different selection regimes for the simultaneous hermaphrodite Macrostomum lignano subjected to sex-limited experimental evolution (i.e. selection for fitness via eggs, sperm, or a control regime allowing both). After 21 and 22 generations of selection for male-specific or female-specific fitness, we characterized changes in whole-organism gene expression. We found that female-selected lines had changed the most in their gene expression. Although annotation for this species is limited, gene ontology term and Kyoto Encyclopedia of Genes and Genomes pathway analyses suggest that metabolic changes (e.g. biosynthesis of amino acids and carbon metabolism) are an important adaptive component. As predicted, we found that the expression of genes previously identified as testis-biased candidates tended to be downregulated in the female-selected lines. We did not find any significant expression differences for previously identified candidates of other sex-specific organs, but this may simply reflect that few transcripts have been characterized in this way. In conclusion, our experiment suggests that changes in testis-biased gene expression are important in the early evolution of sex chromosomes and gonochorism.

Dynamics of seminal fluid production after mating

Seminal fluid proteins (SFPs) play vital roles for optimizing reproductive success in diverse animals. Underlining their significance, SFP production and transfer are highly plastic, e.g., depending on the presence of rivals or mating status of partners. However, surprisingly little is known about replenishing SFPs after mating. This is especially relevant in species that mate multiple times, as they continuously produce and use SFPs throughout their reproductive life. Here we examined the expression pattern of SFP genes after mating in the great pond snail, Lymnaea stagnalis. Our results show that two out of the six SFP genes investigated here were upregulated 1 week after mating. Surprisingly, most SFP genes did not change their expression immediately after mating. Even after 1 week, when supposedly seminal fluid is fully replenished, the expression of SFP genes is rather high. In addition, the difference with previous studies hints at the possibility that SFP production after mating is plastic and depends on the mating history of female-acting snails. Our results shed light on unexplored aspects of SFP production, thereby expanding the understanding of reproductive strategies in animals.

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.

Precocious Sperm Exchange in the Simultaneously Hermaphroditic Nudibranch, Berghia stephanieae

Sexual systems vary greatly across molluscs. This diversity includes simultaneous hermaphroditism, with both sexes functional at the same time. Most nudibranch molluscs are thought to be simultaneous hermaphrodites, but detailed studies of reproductive development and timing remain rare as most species cannot be cultured in the lab. The aeolid nudibranch, Berghia stephanieae, is one such species that can be cultured through multiple generations on the benchtop. We studied B. stephanieae reproductive timing to establish when animals first exchange sperm and how long sperm can be stored. We isolated age- and size-matched individuals at sequential timepoints to learn how early individuals can exchange sperm. Individuals isolated at 10 weeks post initial feeding (wpf; ∼13 weeks postlaying [wpl]) can produce fertilized eggs. This is 6 weeks before animals first lay egg masses, indicating that sperm exchange occurs well before individuals are capable of laying eggs. Our results indicate that male gonads become functional for animals between 6 mm (∼6 wpf, ∼9 wpl) and 9 mm (∼12 wpf, ∼15 wpl) in length. That is much smaller (and sooner) than the size (and age) of individuals at first laying (12–19 mm; ∼16 wpf, ∼19 wpl), indicating that male and female functions do not develop simultaneously. We also tracked the number of fertilized eggs in each egg mass, which remained steady for the first 10–15 egg masses, followed by a decline to near-to-no fertilization. This dataset provides insights into the precise timing of the onset of functionality of the male and female reproductive systems in B. stephanieae. These data contribute to a broader understanding of reproductive development and the potential for understanding the evolution of diverse sexual systems in molluscs.

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.

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.

Larger males facilitate population expansion in Ophraella communa

One of the most intriguing concepts in animal ecology is the reproductive advantages offered by larger body size, and the females prefer to mate with larger males to gain reproductive advantage. Currently, it is not clear how females recognize signs of male 'quality' and what mechanisms are involved in producing offspring with direct or indirect benefits. Our study aims to assess the preferences of females for males in Ophraella communa, determine the reproductive benefits and reveal the underlying mechanism behind this advantage. We demonstrate that male body size is an important determinant in the evolutionary process of O. communa, affecting female mate choice. Moreover, our study establishes that females prefer males with a larger body size, and this could further improve the developmental and reproductive fitness of their offspring. Finally, we focus on the seminal fluid proteins (SFPs) in O. communa, determine differentially expressed genes (i.e. OcACE, OcCBP and OcSFP) by analysing their proteomes and transcriptomes, and define the role of these SFPs-related genes through RNAi. Our study proved that the reproductive benefit of large males may be regulated by biased expression of crucial SFPs genes. The present study advances our understanding of the biological significance of preferential mating.

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.

Seminal fluid and accessory male investment in sperm competition

Sperm production and allocation strategies have been a central concern of sperm competition research for the past 50 years. But during the 'sexual cascade' there may be strong selection for alternative routes to maximizing male fitness. Especially with the evolution of internal fertilization, a common and by now well-studied example is the accessory ejaculate investment represented by seminal fluid, the complex mixture of proteins, peptides and other components transferred to females together with sperm. How seminal fluid investment should covary with sperm investment probably depends on the mechanism of seminal fluid action. If seminal fluid components boost male paternity success by directly enhancing sperm function or use, we might often expect a positive correlation between the two forms of male investment, whereas trade-offs seem more likely if seminal fluid acts independently of sperm. This is largely borne out by a broad taxonomic survey to establish the prevailing patterns of seminal fluid production and allocation during animal evolution, in light of which I discuss the gaps that remain in our understanding of this key ejaculate component and its relationship to sperm investment, before outlining promising approaches for examining seminal fluid-mediated sperm competitiveness in the post-genomic era. This article is part of the theme issue 'Fifty years of sperm competition'.

The hidden ageing costs of sperm competition

Ageing and sexual selection are intimately linked. There is by now compelling evidence from studies performed across diverse organisms that males allocating resources to mating competition incur substantial physiological costs, ultimately increasing ageing. However, although insightful, we argue here that to date these studies cover only part of the relationship linking sexual selection and ageing. Crucially, allocation to traits important in post-copulatory sexual selection, that is sperm competition, has been largely ignored. As we demonstrate, such allocation could potentially explain much diversity in male and female ageing patterns observed both within and among species. We first review how allocation to sperm competition traits such as sperm and seminal fluid production depends on the quality of resources available to males and can be associated with a wide range of deleterious effects affecting both somatic tissues and the germline, and thus modulate ageing in both survival and reproductive terms. We further hypothesise that common biological features such as plasticity, prudent sperm allocation and seasonality of ejaculate traits might have evolved as counter-adaptations to limit the ageing costs of sperm competition. Finally, we discuss the implications of these emerging ageing costs of sperm competition for current research on the evolutionary ecology of ageing.

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.

RNA-Seq of three free-living flatworm species suggests rapid evolution of reproduction-related genes

BACKGROUND

The genus Macrostomum consists of small free-living flatworms and contains Macrostomum lignano, which has been used in investigations of ageing, stem cell biology, bioadhesion, karyology, and sexual selection in hermaphrodites. Two types of mating behaviour occur within this genus. Some species, including M. lignano, mate via reciprocal copulation, where, in a single mating, both partners insert their male copulatory organ into the female storage organ and simultaneously donate and receive sperm. Other species mate via hypodermic insemination, where worms use a needle-like copulatory organ to inject sperm into the tissue of the partner. These contrasting mating behaviours are associated with striking differences in sperm and copulatory organ morphology. Here we expand the genomic resources within the genus to representatives of both behaviour types and investigate whether genes vary in their rate of evolution depending on their putative function.

RESULTS

We present de novo assembled transcriptomes of three Macrostomum species, namely M. hystrix, a close relative of M. lignano that mates via hypodermic insemination, M. spirale, a more distantly related species that mates via reciprocal copulation, and finally M. pusillum, which represents a clade that is only distantly related to the other three species and also mates via hypodermic insemination. We infer 23,764 sets of homologous genes and annotate them using experimental evidence from M. lignano. Across the genus, we identify 521 gene families with conserved patterns of differential expression between juvenile vs. adult worms and 185 gene families with a putative expression in the testes that are restricted to the two reciprocally mating species. Further, we show that homologs of putative reproduction-related genes have a higher protein divergence across the four species than genes lacking such annotations and that they are more difficult to identify across the four species, indicating that these genes evolve more rapidly, while genes involved in neoblast function are more conserved.

CONCLUSIONS

This study improves the genus Macrostomum as a model system, by providing resources for the targeted investigation of gene function in a broad range of species. And we, for the first time, show that reproduction-related genes evolve at an accelerated rate in flatworms.

Effects of two seminal fluid transcripts on post-mating behaviour in the simultaneously hermaphroditic flatworm Macrostomum lignano

The seminal fluid proteins (SFPs) transferred to mating partners along with sperm often play crucial roles in mediating post-mating sexual selection. One way in which sperm donors can maximize their own reproductive success is by modifying the partner's (sperm recipient's) post-copulatory behaviour to prevent or delay re-mating, thereby decreasing the likelihood or intensity of sperm competition. Here, we adopted a quantitative genetic approach combining gene expression and behavioural data to identify candidates that could mediate such a response in the simultaneously hermaphroditic flatworm Macrostomum lignano. We identified two putative SFPs-Mlig-pro46 and Mlig-pro63-linked to both mating frequency and 'suck' frequency, a distinctive behaviour, in which, upon ejaculate receipt, the worm places its pharynx over its female genital opening and apparently attempts to remove the received ejaculate. We, therefore, performed a manipulative experiment using RNA interference-induced knockdown to ask how the loss of Mlig-pro46 and Mlig-pro63 expression, singly and in combination, affects mating frequency, partner suck propensity and sperm competitive ability. None of the knockdown treatments impacted strongly on the mating frequency or sperm competitive ability, but knockdown of Mlig-pro63 resulted in a significantly decreased suck propensity of mating partners. This suggests that Mlig-pro63 may normally act as a cue in the ejaculate to trigger recipient suck behaviour and-given that other proteins in the ejaculate have the opposite effect-could be one component of an ongoing arms race between donors and recipients over the control of ejaculate fate. However, the adaptive significance of Mlig-pro46 and Mlig-pro63 from a donor perspective remains enigmatic.

The free-living flatworm Macrostomum lignano

Macrostomum lignano is a free-living flatworm that is emerging as an attractive experimental animal for research on a broad range of biological questions. One feature setting it apart from other flatworms is the successful establishment of transgenesis methods, facilitated by a steady supply of eggs in the form of single-cell zygotes that can be readily manipulated. This, in combination with the transparency of the animal and its small size, creates practical advantages for imaging and fluorescence-activated cell sorting in studies related to stem cell biology and regeneration. M. lignano can regenerate most of its body parts, including the germline, thanks to the neoblasts, which represent the flatworm stem cell system. Interestingly, neoblasts seem to have a high capacity of cellular maintenance, as M. lignano can survive up to 210 Gy of γ-irradiation, and partially offset the negative consequence of ageing. As a non-self-fertilizing simultaneous hermaphrodite that reproduces in a sexual manner, M. lignano is also used to study sexual selection and other evolutionary aspects of sexual reproduction. Work over the past several years has led to the development of molecular resources and tools, including high-quality genome and transcriptome assemblies, transcriptional profiling of the germline and somatic neoblasts, gene knockdown, and in situ hybridization. The increasingly detailed characterization of this animal has also resulted in novel research questions, such as bio-adhesion based on its adhesion-release glands and genome evolution due to its recent whole-genome duplication.

Seminal fluid-mediated fitness effects in the simultaneously hermaphroditic flatworm Macrostomum lignano

As a class, seminal fluid proteins are expected to exert strong effects on mating partners due to the selection pressures of sperm competition and sexual conflict. But because of the complexity of this secretion, linking specific proteins to downstream effects on own fitness-via manipulating the reproductive behavior, physiology, and ultimately the sperm utilization of mating partners-is not straightforward. Here, we adopted a systematic gene knockdown approach to screen for seminal fluid-mediated fitness effects in the simultaneously hermaphroditic flatworm Macrostomum lignano. We focused on 18 transcripts in M. lignano seminal fluid, testing how their RNA interference-induced knockdown impacted on three aspects of donor (male) reproductive success: (a) fertility (offspring production of the partner); (b) defensive sperm competitive ability, P 1; and (c) offensive sperm competitive ability, P 2. In general, the knockdown of most individual transcripts appeared to have only a minor impact on male reproductive success, though we found evidence that the knockdown of up to five different transcripts impacted on fertility; the knockdown of two other transcripts resulted in reduced P 2; and knockdown of a further transcript actually increased P 2. We thus identify a number of candidate seminal fluid transcripts that appear to modulate offspring production and sperm competitiveness in M. lignano. That only a minority of transcripts exhibit such a pattern likely reflects both the difficulty of accurately estimating sperm competitiveness and the functional redundancy of seminal fluid.

Genetic and environmental variation in transcriptional expression of seminal fluid proteins

Seminal fluid proteins (SFPs) are crucial mediators of sexual selection and sexual conflict. Recent studies have chiefly focused on environmentally induced plasticity as one source of variation in SFP expression, particularly in response to differing sperm competition levels. However, understanding the evolution of a trait in heterogenous environments requires estimates of both environmental and genetic sources of variation, as well as their interaction. Therefore, we investigated how environment (specifically mating group size, a good predictor of sperm competition intensity), genotype and genotype-by-environment interactions affect seminal fluid expression. To do so, we reared 12 inbred lines of a simultaneously hermaphroditic flatworm Macrostomum lignano in groups of either two or eight worms and measured the expression levels of 58 putative SFP transcripts. We then examined the source of variation in the expression of each transcript individually and for multivariate axes extracted from a principal component analysis. We found that mating group size did not affect expression levels according to the single transcript analyses, nor did it affect the first principal component (presumably representing overall investment in seminal fluid production). However, mating group size did affect the relative expression of different transcripts captured by the second principal component (presumably reflecting variation in seminal fluid composition). Most transcripts were genetically variable in their expression level and several exhibited genotype-by-environment interactions; relative composition also showed high genetic variation. Collectively, our results reveal the tightly integrated nature of the seminal fluid transcriptome and provide new insights into the quantitative genetic basis of seminal fluid investment and composition.