Sperm competition risk drives plasticity in seminal fluid composition
- Steven A. Ramm†1, 2,
- Dominic A. Edward†1,
- Amy J. Claydon1, 3,
- Dean E. Hammond3,
- Philip Brownridge3,
- Jane L. Hurst1,
- Robert J. Beynon3 and
- Paula Stockley1Email author
© Ramm et al. 2015
Received: 20 August 2015
Accepted: 7 October 2015
Published: 27 October 2015
Ejaculates contain a diverse mixture of sperm and seminal fluid proteins, the combination of which is crucial to male reproductive success under competitive conditions. Males should therefore tailor the production of different ejaculate components according to their social environment, with particular sensitivity to cues of sperm competition risk (i.e. how likely it is that females will mate promiscuously). Here we test this hypothesis using an established vertebrate model system, the house mouse (Mus musculus domesticus), combining experimental data with a quantitative proteomics analysis of seminal fluid composition. Our study tests for the first time how both sperm and seminal fluid components of the ejaculate are tailored to the social environment.
Our quantitative proteomics analysis reveals that the relative production of different proteins found in seminal fluid – i.e. seminal fluid proteome composition – differs significantly according to cues of sperm competition risk. Using a conservative analytical approach to identify differential expression of individual seminal fluid components, at least seven of 31 secreted seminal fluid proteins examined showed consistent differences in relative abundance under high versus low sperm competition conditions. Notably three important proteins with potential roles in sperm competition – SVS 6, SVS 5 and CEACAM 10 – were more abundant in the high competition treatment groups. Total investment in both sperm and seminal fluid production also increased with cues of heightened sperm competition risk in the social environment. By contrast, relative investment in different ejaculate components was unaffected by cues of mating opportunities.
Our study reveals significant plasticity in different ejaculate components, with the production of both sperm and non-sperm fractions of the ejaculate strongly influenced by the social environment. Sperm competition risk is thus shown to be a key factor in male ejaculate production decisions, including driving plasticity in seminal fluid composition.
KeywordsEvolution Phenotypic plasticity Proteomics Random Forest analysis Seminal fluid Sperm competition
Ejaculates are a complex mixture of sperm and seminal fluid components, many features of which are thought to have been shaped by sperm competition [1, 2]. For example, much recent research has focused on the optimal allocation of limited sperm reserves to maximise fertilisation success [3–5]. A related but separate question concerns the extent to which overall investment in ejaculate production varies with sperm competition. There is widespread support for the prediction that higher levels of sperm competition select for greater investment in sperm production, both across [6–9] and within [10–14] species, including the ability to plastically adjust sperm production to match the prevailing competitive conditions [15–19]. However, seminal fluid contains a diverse array of proteins that exert various influences on female reproductive biology, and it is likely that varying levels of sperm competition could also impact on optimal seminal fluid production and/or composition [20–24]. A further important consideration is that investment in ejaculate production might also be influenced by the number of females available for mating [25, 26]. This complicates interpretation of previous studies because sperm competition and mating rate will often be confounded [27, 28] and under certain conditions, mating rate might predominate in determining optimal male sperm production strategies .
Here we use an established vertebrate model system, the house mouse (Mus musculus domesticus), to investigate plasticity in production of different ejaculate components in relation to cues of sperm competition risk and mating rate. House mice have a relatively low overall investment in sperm production compared to many other rodents . However, males typically experience a significant but variable risk of sperm competition [30, 31], which previously has been linked to phenotypic plasticity in sperm production [17, 32]. Moreover, we have shown that there is a rapid turnover of seminal fluid proteins in this species, notably more rapid than spermatogenesis , suggesting potential for short-term plasticity in responses of non-sperm ejaculate components. Mating opportunities are also likely to vary substantially for male house mice in natural populations, depending on the number of resident females in the territory of a particular male and the proximity of neighbouring territories [30, 34].
To determine the relative importance of sperm competition risk and anticipated mating rates on phenotypic plasticity in ejaculate production requires a controlled experiment with independent manipulation of appropriate cues used to assess these two factors. In rodents, scent represents a crucial sensory modality used in a variety of reproductive contexts . For example, male rodents are sensitive to cues of sperm competition conveyed by conspecific male odours [17, 21, 36] and are capable of sophisticated discrimination of odours from female conspecifics [37, 38]. The perception of sperm competition risk and mating opportunities (hereafter ‘mating rate’) should therefore be influenced by exposure to the odours of adult male and female conspecifics, respectively, in addition to direct encounters with them. That is, a high encounter rate with multiple males and/or their odours should indicate a high risk of sperm competition, whereas a high encounter rate with multiple females and/or their odours should indicate a high potential mating rate.
To test whether male house mice vary investment in the production of different ejaculate components according to social conditions, we used a 2 × 2 factorial experimental design and independently manipulated cues of sperm competition risk and potential mating rate. Combined with detailed quantitative proteomics analysis, our experiment reveals that the relative abundance of different proteins found in seminal fluid (i.e. seminal fluid composition) and sperm production both differ according to cues of sperm competition risk. By contrast, relative investment in different ejaculate components was unaffected by cues of potential mating rate. Sperm competition risk is thus shown to be a key factor influencing male ejaculate production decisions that drives plasticity in seminal fluid composition.
The experiment simultaneously manipulated subjects’ social experience of both male and female conspecifics in a 2 × 2 factorial design to simulate both high or low sperm competition risk (HSC/LSC; subjects were provided with cues from three or one male conspecifics, respectively) and high or low anticipated mating rates (HMR/LMR; subjects were provided with cues from four or two female conspecifics, respectively – see detailed description in the Methods).
Mean sperm production parameters under contrasting cues of sperm competition risk and mating opportunities
Sperm competition (SC)
Mating rate (MR)
Paired testes mass (g)
Epididymal sperm count (x106)
Daily sperm production (x106)
Seminal vesicles mass (g)
0.196 ± 0.007
6.72 ± 0.32
3.04 ± 0.24
0.151 ± 0.015
0.188 ± 0.009
6.69 ± 0.47
3.50 ± 0.33
0.139 ± 0.013
0.169 ± 0.006
5.36 ± 0.27
2.82 ± 0.14
0.117 ± 0.012
0.174 ± 0.012
5.78 ± 0.46
2.60 ± 0.16
0.115 ± 0.012
Sperm competition risk but not potential mating rate impacts on ejaculate production
Paired testes mass
Epididymal sperm count
Daily sperm production
Seminal vesicles mass
SC × MR
These results are supported by mixed model analyses taking into account potential non-independence of males from cages grouped within the same enclosure. The effects of competition remained either significant (epididymal sperm numbers: χ2 = 4.07, d.f. = 1, P = 0.044) or very nearly so (daily sperm production: χ2 = 3.57, d.f. = 1, P = 0.059; testis mass: χ2 = 3.44, d.f. = 1, P = 0.063). Taken together with our earlier findings , these data strongly support the conclusion that sperm competition cues significantly impact upon sperm production parameters (Fig. 1), and that this effect is most pronounced on epididymal sperm numbers and daily sperm production rates, and somewhat weaker on the size of the testes.
Neither sperm competition nor mating rate manipulations affected measures of ejaculate quality, such as percentage sperm motility (treatment group effect: F 3,28 = 2.19, P = 0.1), curvilinear sperm velocity (F 3,28 = 1.11, P = 0.4) or straight-line sperm velocity (F 3,28 = 1.72, P = 0.2).
Seminal fluid production
Seminal fluid composition
Progenesis QI analysis of individual proteins
Random Forest analysis of proteome composition
Next, we employed a more stringent approach to analysing the seminal vesicle proteome, performing a single analysis of proteome composition and focusing solely on that proportion of the proteome that is secreted and that could thus have functions relevant to sperm competition outcomes within the female reproductive tract. To do so, we implemented a Random Forest (RF) classification model for proteins secreted by the seminal vesicles (see Methods). In contrast to the standard Progenesis QI analysis reported above (see Additional file 1), for this analysis we used only known ejaculated proteins (from ) to normalise the samples, again using a minimum of three unique peptides per protein for identification and quantification (see Additional file 2). Normalising to ejaculated proteins is a more conservative approach that controls for any potential differences between treatments in the proportion of secreted proteins present in each sample (see Methods). This might occur, for example, if larger seminal vesicles associated with high competition were to contain a higher ratio of lumen to tissue, and hence a higher proportion of ejaculated to non-ejaculated proteins. We used this approach for analysis of all four treatment groups, to investigate effects of both sperm competition risk (number of males encountered) and potential mating rate (number of females encountered).
Random forest analysis of plasticity in seminal fluid composition
95 % CI
45.14 - 47.63
28.99 - 31.76
26.60 - 29.61
10.64 - 14.67
10.16 - 13.86
7.95 - 12.80
7.75 - 12.02
PATE 4 / SVS 7
4.52 - 8.51
3.56 - 8.09
Serine protease inhibitor A3K
3.21 - 7.92
2.65 - 7.28
2.19 - 6.88
Metalloproteinase inhibitor 1
−0.17 - 4.98
−0.16 - 4.01
−0.15 - 4.67
−1.16 - 3.92
−1.48 - 2.25
−2.11 - 2.15
−2.44 - 2.50
−2.70 - 2.07
−2.24 - 0.99
−3.77 - 0.16
Nucleoside diphosphate kinase B
−3.71 - 0.30
−4.80 - 0.14
78 kDa glucose-regulated protein
−5.36 - -1.08
−5.83 - -1.85
−6.64 - -2.85
−7.48 - -3.23
−7.88 - -3.58
−8.10 - -3.49
−8.7 - -5.00
To determine whether males from the same enclosure tended to be more similar, and thus whether our results may be affected by a common environment effect, we also generated a RF classification model that analysed differences in protein abundance that might have arisen between enclosures. The model performed poorly at classifying males to the correct enclosure, having a classification error rate of 100 %. Additionally, when analysing the differences between males described in the proximity matrix of this model, males were not more similar to their enclosure partners than individuals from other enclosures.
Finally, we investigated the protein composition of seminal fluid using an RF classification model that analysed differences in protein abundances between samples according to the potential mating rate treatment (number of females encountered). This RF model had a classification error rate of 63.3 % (11 correct: 19 incorrect) that could easily have been achieved by chance (968 of 1,000 bootstrap replicates achieved this error rate or lower; P = 0.97). This shows that cues of potential mating rate did not influence the relative abundance of seminal vesicle proteins, as can be visually inferred from the absence of separation of samples with respect to mating rate treatment in the multidimensional scaling plot summarising the RF analysis (Fig. 6b).
Our findings reveal that social cues of sperm competition risk are a critical factor in determining male investment in ejaculate traits. Whilst previous studies have largely focused on the impact of social conditions on ejaculate allocation [4, 5, 20, 43, 44] (but see, e.g. [15–17]), we demonstrate here the importance of the social environment on male ejaculate production strategies. We show that: (1) sperm competition risk cues strongly influence production of both seminal fluid proteins and sperm; but (2) cues of potential mating opportunities do not appear to influence ejaculate production in this species.
Sperm competition risk drives plasticity in seminal fluid composition and investment
Using a proteomics approach to identify and quantify seminal vesicle proteins, we found strong evidence that cues of sperm competition risk influence seminal fluid proteome composition. The response of SVS proteins to cues of sperm competition risk are of particular interest. Collectively, the SVS proteins make up the vast majority of the seminal vesicle secretion; they are likely to be important to male competitive fertilisation success due to established roles in copulatory plug formation, sperm transport and sperm motility (reviewed in [42, 45]), and many exhibit interspecific evidence of rapid evolution [46–48]. Although caution is required in interpreting results of our Progenesis QI analyses, overall these suggest a consistent trend of up-regulated SVS proteins under high competition conditions. By contrast, the 24 other ejaculated proteins that we investigated showed no such directional bias (15 were more abundant in the high competition treatment, and 9 in the low competition treatment). This pattern is supported by results of our Random Forest analyses showing significant differences in secreted seminal fluid proteome composition, with SVS 5 and SVS 6 among the seven proteins identified as important for classifying samples according to sperm competition treatment group. Although a role in sperm competition for SVS 5 and SVS 6 appears likely, elucidating their precise function requires further investigation. Roles in relation to copulatory plugs could be one possibility. However, neither SVS 5 nor SVS 6 contain transglutaminase substrate domains, and so are unlikely to be actively involved in copulatory plug formation . Interestingly, these proteins have also been suggested to be serine protease inhibitors , based on their proximity in the genome to other WFDC proteins, some of which have this activity. However, SVS 5 and SVS 6 do not possess the classic four disulphide bond potential, and homology modelling to predict a three dimensional protein structure  leads to low confidence models that do not map to known serine protease inhibitors. There has been no direct demonstration of anti-protease activity of these proteins, and although the role of SVS 5 in copulatory plug survival has recently been questioned , this does not preclude other functions in sperm competition unrelated to plug formation or dissolution. Also, in relation to copulatory plug formation, we note that the coagulum transglutaminase TGM4 – which catalyses the cross-linking of SVS proteins – was not itself significantly up-regulated under high sperm competition risk according to our initial Progenesis QI analysis.
Although we found limited evidence for an influence of the social environment on proteins related to copulatory plug formation, there is some evidence for an up-regulation of proteins linked to sperm motility under conditions of high sperm competition risk. For example, according to our initial Progenesis QI analysis, SVS 7 (caltrin, PATE4), a protein that enhances sperm motility , was significantly up-regulated in our high competition treatment groups. Further support for this idea comes from more robust evidence of increased CEACAM 10 production under high competition conditions. CEACAM 10 is a glycoprotein that enhances sperm motility , and its increased production under high sperm competition risk is also confirmed by our Random Forest analysis.
Our Random Forest analysis indicates that two other proteins – cathepsin D, an aspartic protease, and annexin A5 (function unknown) – were down-regulated under conditions of high sperm competition (in addition DPP3 was down-regulated under high competition according to our initial Progenesis QI analysis). These responses are less easily explained, although it is possible that reduction of the potential proteolytic environment of the seminal vesicle may exert other effects influencing reproductive success under sperm competition .
In summary, these findings based on quantitative protein analysis in a model vertebrate provide direct and broad-scale evidence of plasticity in seminal fluid protein production linked to sperm competition risk. We are aware of only one previous experimental study that attempted to measure plasticity in seminal fluid production in response to sperm competition cues, and this measured gene rather than protein expression levels. Fedorka et al.  manipulated the social environment of male fruit flies (Drosophila melanogaster), and found evidence for plastic expression of two out of three seminal fluid protein-coding genes investigated – Acp26Aa and Acp62f (but not Acp70A) – both of which were down-regulated 72 h after eclosion under high competition conditions.
A further novel finding of the present study is that male house mice appear to adjust investment in the overall size of the seminal vesicles according to cues of average sperm competition risk. The marginal significance of this response means the result should be interpreted with caution, especially because no similar response was found in our previous study . Nevertheless, it seems likely that altered male investment in the size of the seminal vesicles, and thus presumably in the total amount of seminal fluid production, could be adaptive under varying sperm competition risk. Indeed, recent evidence points to an influence of seminal vesicle mass on competitive male reproductive success in semi-natural house mouse populations , as would be expected if variation in seminal fluid production is relevant to sperm competition outcomes. Socially-mediated plasticity in seminal vesicle size in a second rodent species, the bank vole Myodes glareolus , the correlated evolution of sperm and seminal fluid investment across rodent species , and the rapid evolutionary dynamics of specific seminal fluid components [42, 46–48], all serve to further emphasize the likely importance of non-sperm ejaculate components in mammalian sperm competition. Future experimental work should now focus on the fitness consequences of variation in total and individual seminal fluid protein production (cf. ).
Sperm competition risk drives sperm production plasticity
Male house mice in this study adjusted overall investment in sperm production according to cues of average sperm competition risk experienced during sexual development. Specifically, those males that were repeatedly exposed to social cues from three rival males had higher daily sperm production rates and higher epididymal sperm counts compared to those that had encountered social cues from just one rival male. Importantly, these changes occurred irrespective of whether subject males were regularly also exposed to social cues from two or four potential mates, suggesting that cues of likely mating rate at this level do not affect sperm production plasticity in this species (see below).
Our findings are generally consistent with results from two previous studies in which social experience of male house mice with male (but not female) conspecifics was manipulated to vary cues of average sperm competition level, and where sperm production was also higher following regular encounters with more male conspecifics [17, 32]. In the present study, and contrary to previous findings [17, 32], we additionally report a moderate but significant influence of male social experience on the testis mass of male house mice, consistent with widely held assumptions linking testis size to sperm production rates (see  and references therein). Nonetheless, we also found a disproportionate increase in the other sperm production parameters (epididymal sperm numbers and daily sperm production based on spermatid head counts from testicular homogenates) relative to testis size, a pattern that is highly consistent across studies . This implies that an increased testis size alone cannot explain all of the sperm production plasticity response seen in male house mice. Similar evidence for invertebrates suggests that this may be a common feature of phenotypic plasticity in sperm production (e.g. ), although the precise adjustments within the testis that enable males to plastically increase investment in sperm production under heightened sperm competition remain to be elucidated (reviewed in ; see also [58–62] for evolutionary responses to sperm competition besides gross testis size).
In addition to sperm production parameters, we also investigated a measure of sperm quality. Despite recent evidence that mouse sperm motility parameters respond to different post-copulatory sexual selection regimes under experimental evolution , we found no evidence of plasticity in sperm motility parameters in response to cues of either sperm competition risk or likely mating rate, at least under in vitro conditions, consistent with earlier findings in this species .
Female-mediated cues do not impact on ejaculate production
Although male house mice in our study adjusted overall investment in ejaculate production according to their encounters with male conspecifics, we found no evidence that ejaculate investment increased as predicted in response to encounters with female conspecifics. Here we predicted a potential response to likely mating rates: males that regularly encounter more females might perceive more potential mating opportunities and increase their sperm and seminal fluid production accordingly. Since no such effect was observed, our findings suggest that variation in mating rate may not be as important as sperm competition in determining optimal ejaculate investment decisions for male house mice, at least at the relatively modest but naturalistic levels of variation in potential mating opportunities presented within our experimental treatment groups. In general, male house mice appear relatively constrained in the number of ejaculates they are able to produce within a given 24 hour period [63–65], such that some degree of sperm limitation may be expected to occur where two or more females within their territory are simultaneously receptive . Cues encountered from four females might, therefore, reasonably be regarded as offering a potentially high mating rate for male house mice, although we cannot rule out the possibility that more extreme variation in female encounter rates might reveal a significant response in sperm production traits. For example, an experimental evolution study in Drosophila melanogaster revealed that mating rate can significantly impact on sperm production parameters under relatively extreme female-biased sex ratios  (cf. ). Similar effects could also explain why some other experimental evolution studies have not observed an evolutionary response in testis size in male-biased populations, where sperm competition is elevated but mating rate is simultaneously reduced . In contrast to these previous studies however, the cues of likely sperm competition and mating rate that we manipulated in the present study were indirect; hence, it is possible that repeated emptying and re-filling of sperm reserves under a realised elevation in mating rate might affect sperm production parameters differently (see also [16, 55]).
In conclusion, our study reveals significant plasticity in different ejaculate components, with both the relative abundances of different seminal fluid proteins and sperm production strongly influenced by the social environment. Moreover, we show that experience with male conspecifics – and thus the perceived level of sperm competition – predominates over experience with female conspecifics in determining variable investment in ejaculate components. This confirms sperm competition risk as a key factor in male ejaculate production decisions. We also show here for the first time that this risk drives plasticity in seminal fluid composition as well as sperm numbers. Our results, therefore, emphasise the importance of considering the whole ejaculate when seeking to understand male responses to sperm competition [2, 69].
To distinguish the potential effects of sperm competition risk and potential mating rate cues on sperm production plasticity, we exposed recently-weaned experimental subjects to long-term treatments that manipulated their social experience with both male conspecifics (as a cue of sperm competition risk) and female conspecifics (as a cue of anticipated mating rate). The manipulations were achieved by closely controlling the social experience of subjects, including controlled physical encounters with both male and female conspecifics and regular exposure to conspecific odours.
Subject males (n = 32) were from a colony of wild house mice that had been outbred for six or fewer generations in captivity and originally derived from local populations in Cheshire. Each male was individually housed in a 48 cm × 11.5 cm × 12 cm cage for the duration of the experiment (M3, North Kent Plastic Cages Ltd., UK), with Corn Cob Absorb 10/14 substrate and paper-wool nest material, with ad libitum access to food (LabDiet 5002) and water. Subjects were maintained under controlled environmental conditions: temperature 20–21 °C, relative humidity 45–65 % and a reversed 12:12 h light cycle (lights off at 08.00). Females used for the experiments (n = 36) were unrelated, sexually mature individuals from the same colony, housed in (unrelated) pairs in M3 cages distributed on a rack in the same room as the experimental males. All social experience manipulations were performed during the dark period.
Cues of sperm competition risk and potential mating rate
The experiment was designed to simultaneously manipulate cues of sperm competition risk and potential mating rate. This was achieved by a 2 × 2 factorial design incorporating contrasting levels of both sperm competition (high sperm competition (HSC) vs. low sperm competition (LSC)) and mating rate (high mating rate (HMR) vs. low mating rate (LMR)), with n = 8 males per block (n = 32 in total) and social experience manipulated over a period of ten weeks. Recently weaned subject male mice (age at start of experiment: mean 32.9 days, range 31-33 days) were randomly allocated to treatment groups, whilst ensuring brothers were not grouped together and with minor modifications to ensure body masses were balanced at the start of the experiment (mean ± SEM; HSC, LMR: 15.76 ± 0.74 g HSC, HMR: 16.31 ± 0.78 g; LSC, LMR: 15.88 ± 0.56 g; LSC, HMR: 16.18 ± 0.42 g). Males did not differ in body mass between treatment groups at the end of the experiment (mean ± SEM; HSC, LMR: 20.14 ± 0.83 g; HSC, HMR: 20.11 ± 0.89 g; LSC, LMR: 19.83 ± 0.68 g; LSC, HMR: 18.99 ± 0.78 g; F 3,28 = 0.44, P = 0.7).
Each experimental male was placed (within its home cage) inside a high-sided enclosure (each 1.25 m × 0.6 m × 0.8 m high) with either one (LSC treatment groups) or three (HSC treatment groups) other males present in the same enclosure (also within their home cages). All males had the opportunity to interact with the other males in their enclosure during weeks 1, 3 and 5 of the experiment. To do this we released each male in turn from their home cage for a period of 30 min each. This ensured that males were always separated by cage bars to prevent escalation of any aggressive interactions. In HSC treatment groups, this contact experience was supplemented by the transfer of odour cues (approximately 25 g soiled bedding) among experimental cages three times a week for all ten weeks of the experiment. This followed a weekly schedule whereby each experimental male received one odour cue from each of the three other experimental males housed in the same enclosure. A similar schedule was followed for the LSC treatment groups except that odour cues from the competitor male were transferred once a week. On the other two transfer days each male received a sham odour transfer from their own cage to control for cage handling between treatments.
Cues of likely mating rate were manipulated by transferring female odour cues (approximately 25 g soiled bedding) to subject male cages. Odour cues came from either one cage containing a pair of females (LMR) or from two cages each containing a pair of females (HMR). Transfers were carried out on two days each week, each time with odour cues from two females transferred in the LMR treatment, and from four females in the HMR treatment. Female odour cues were transferred on different days to the transfer of male odour cues, for the full ten weeks of the experiment. Each pair of females contributed odour cues consistently to only one male enclosure (containing either two or four males, depending on the SC treatment). We also allowed males to interact with females by placing the female cages within the high-sided enclosures and releasing each male in turn from its home cage, during which time the other male cages were removed. This was done for a period of 30 min per male in weeks 2, 4 and 6 of the experiment. Male and female cages were cleaned (including replacement of all soiled substrate and nest material with fresh substrate and nest material) in weeks 2, 4, 6, 8, and 10. Ejaculate production parameters were assessed by removing subjects from the experiment in a randomised order at the end of week 11 and during week 12.
Seminal fluid investment
The seminal vesicles are the largest of the accessory reproductive glands in rodents and the major site of seminal fluid production . At the end of the experiment, we therefore recorded the mass of the left and right seminal vesicles, as a measure of investment in seminal fluid production. Samples were then frozen at −20 °C for later analysis.
Proteomics and label-free quantification of seminal fluid
Whole seminal vesicles were thawed then homogenised on ice in 20x volume 50 mM ammonium bicarbonate (based on fresh weights). Each homogenate was assayed for protein concentration; 50 μg of protein per sample was proteolysed with trypsin, after treatment with dithiothreitol and iodoacetamide to reduce and alkylate cysteine residues. The final digestion volume was 200 μL. Digests were further diluted to ensure an approximate concentration of 0.8 μg/μL. Samples were analysed for 30 of the 32 experimental subjects (samples for two subjects suffered degradation as a result of being defrosted prematurely and refrozen).
To quantify the abundance of specific seminal fluid proteins, mass spectrometry was performed as described previously . Briefly, the samples, as prepared above, were analysed as tryptic peptides, with 1 μL of the final dilution resolved by reversed-phase (C18) ultra performance liquid chromatography (Waters nanoAcquity) over a 30 min linear organic gradient of 3–40 % buffer B (0.1 % formic acid in acetonitrile), prior to tandem mass spectrometry using a LTQ-Orbitrap Velos (Thermo Scientific, Waltham, MA, USA). High resolution, accurate mass data were acquired in a data-dependent manner, with the top 20 most intense peptides in each MS scan selected for fragmentation.
The raw data were processed using Progenesis QI (v2, Nonlinear Dynamics) to determine protein abundances. The data from all raw files were automatically aligned according to retention time to produce an aggregate spectrum, from which charge states +1 and > +4 were excluded. This aggregated spectrum contains tandem MS fragmentation data from all aligned runs enabling maximal protein identifications across all samples. Data were then separated into four experimental treatment groups, with eight males acting as biological replicates for both LSC groups and seven males for the HSC conditions. The aggregate peak list file (.mgf) was created using standard default settings and specified to contain only data relating to peptides ranked 1–4. This was then searched against a protein database containing reviewed UniProt entries for Mus musculus plus a small number of unreviewed entries relating to proteins identified elsewhere , using Mascot (v 2.5.1). A 10 ppm peptide tolerance and 0.5 Da MS/MS tolerance were set, with modifications of fixed cysteinyl carbamidomethylation and variable oxidation of methionine. Trypsin was the specified enzyme, allowing for one miscleavage. The .xml file generated by Mascot was imported back into Progenesis QI for feature assignment, normalisation and relative quantification. Data were normalised 1) to all proteins identified in the seminal vesicle samples (see Additional file 1 – data are split by treatment groups according to high or low sperm competition risk), and 2) to ejaculated protein(s) in the list identified by Dean et al.  (see Additional file 2 – data are split by treatment groups according to high or low sperm competition risk, and high or low potential mating rate). Normalisation to all proteins is the default approach; however, in the latter case we normalised to ejaculated proteins because there was a trend for males in our high competition treatment groups to develop larger seminal vesicles than those in the low competition groups. This is important to take into account because even though we analysed an equivalent total amount of protein for each male, larger seminal vesicles might contain a different ratio of lumen to tissue, and hence a different ratio of ejaculated to non-ejaculated proteins. Protein quantification was based on averaging the individual abundances of the top three most abundant peptides for each protein and then comparing them relatively across runs. Summaries of the quantitative proteomics data and Progenesis QI analyses are provided in Fig. 1 and Additional files 1 and 2, respectively.
Statistical analyses of seminal fluid composition
A total of 383 proteins were identified in the seminal vesicles from 30 experimental subjects, although our analysis considers only 147 of these that were quantified using at least three peptides (see Additional files 1 and 2). The seminal vesicle samples analysed contain a complex mixture of proteins at widely varying concentrations. In addition to proteins found in the lumen of the seminal vesicle, these 147 proteins also include seminal vesicle tissue proteins. So that only ejaculated proteins were considered, we compared our 147 proteins against the list of 69 male-derived proteins identified by Dean et al.  that are contained in the male ejaculate, and limited further analysis to the 31 proteins common to both datasets. Note that we do not expect to identify all 69 ejaculated proteins identified by Dean et al.  as these will also be sourced from other reproductive tissues, such as the prostate. By limiting our analysis to 31 proteins known to be transferred to the female at ejaculation we therefore focus here on proteins that could potentially influence male success in sperm competition.
Abundances were analysed using Random Forest (RF) classification models . RF is a classification algorithm highly suited to the analysis of high-dimensional proteomics and other -omics datasets . In brief, RF creates an ensemble of decision trees (a forest) in which each tree is trained to classify a set of samples to different classes (e.g. experimental treatment groups) using a set of variables (e.g. protein abundances). If the abundance of a protein, or proteins, differs between experimental treatment groups, the RF model will be more accurate at classifying samples to the correct class. Each tree in the forest is trained using a different random subset of samples and proteins from the entire dataset. Remaining samples not used in the training stage are termed ‘out of bag’ (OOB) samples. After training, each tree is used to predict the experimental treatment of each OOB sample and these are subsequently averaged across all trees to give an overall classification for each sample. The predicted classifications are then compared to the actual experimental treatment of each sample to determine the overall accuracy of the RF model. If the RF model can accurately predict the experimental treatment of each sample, this indicates differences in the abundance of proteins between treatments.
The results of the RF model can be used to identify which proteins differ most between classes or experimental treatments. This is done by extracting variable importance scores for each protein. These importance scores describe the contribution of each protein to the accurate classification of samples across all the trees in the forest. Proteins that differ most between treatments will make a greater contribution to the accurate classification of samples and consequently have higher variable importance scores.
RF classification models were implemented using the randomForest package  in R v2.15.1 (R Foundation for Statistical Computing, Vienna, Austria). Each RF model contained 10,000 trees. The number of variables sampled in each tree was determined using the tuneRF function in randomForest. Outliers were defined as samples that consistently achieved an outlier score greater than 10 .
Two RF models were generated to test for differences based on treatment group. The first model classified samples according to the number of males encountered by subjects in each group whilst the second model classified samples according to the number of females encountered by subjects in each group. As RF is a stochastic modelling approach, that is, each tree is based on a random subset of samples and variables, each RF model was generated 1,000 times. OOB classification errors are reported as the median error rate across all 1,000 replicates of each model. The significance of OOB classification errors was determined by performing a further 1,000 bootstrap replicates of each model with random permutations of the classification variable.
Variable importance (VI) scores are reported as the mean decrease in classification accuracy attributable to each protein. VI scores were averaged and 95 % confidence intervals estimated from the 1,000 replicates of each RF model. Variables are defined as being important for differentiating samples only when the VI score is greater than the absolute value of the lowest VI score. VI scores were only reported for a model if the OOB classification error rate of the model was significantly greater than expected by chance.
RF models are highly suited to the analysis of proteomics datasets . Advantages of this approach include the ability to deal with a larger number of variables relative to sample size, the test is non-parametric, models are not sensitive to outliers, overfitting is not a problem, and variable importances and sample proximities are easily generated. An alternative method of analysis to identify significant differences in individual protein abundances between treatments would be to perform a series of statistical tests, for example t-tests, with subsequent correction for multiple comparisons. However, there is considerable debate over the correct method to adopt to control for multiple comparisons. Broadly speaking, methods control either for the familywise error rate (e.g. Bonferroni, Holm-Bonferroni, Hochberg-Bonferonni) or false discovery rate (FDR; e.g. Benjamini-Hochberg, Benjamini-Hochberg-Yekutieli). The chance of type I error, that is, falsely rejecting the null hypothesis, will depend on the method adopted. To validate our analysis and provide a comparison of different methods available to analyse our data, we have also conducted multiple t-tests to analyse differences according to the sperm competition treatment (number of males). P-values resulting from t-tests were corrected for multiple comparisons by the FDR using the Benjamini-Hochberg method to yield q-values.
At the end of the experiment we determined three sperm production parameters: testes mass (immediately, combined mass of both testes), epididymal sperm numbers (immediately, based on the caput of the right epididymis), and daily sperm production rate (using frozen samples of the right testis), as previously described . The latter measure is based on spermatid head counts from testicular homogenates which, because the timing of spermatogenesis in mice is known, can be converted into a dynamic measure of sperm output per testis per day (see [17, 72]). One male was excluded from the analysis of daily sperm production rate due to problems processing frozen material. We also recorded male body mass at the beginning and end of the experiment. All measurements were taken blind to the treatment group of the experimental subject being measured.
We estimated sperm motility parameters for subject males based on sperm recovered from the left epididymis, as described by Lemaître et al. . Briefly, at the end of the experiment sperm were collected by making ten incisions in the cauda epididymis with a scalpel blade into 200 μL BWW medium . This was incubated for 15 min at 37 °C and a 20 μl sample then transferred to a glass slide and a cover slip applied. For each sample we obtained two video recordings of sperm, each of 2 s duration and captured at 75 frames/s using a Leica DM1000 microscope with heated stage set to 37 °C coupled to a Point Grey Flea2 (FL2-03S2M-C) 1394b camera and FlyCapture® software (Point Grey Research Inc., Richmond, BC, Canada). Videos were analysed in ImageJ  with the CASA plugin  using parameters optimised for mouse sperm.
Statistical analyses of sperm data
Effects of the two experimental factors (sperm competition risk, potential mating rate) and their interaction were assessed using general linear models, with data log-transformed as appropriate. We also conducted mixed model analyses to control for potential non-independence of subjects whose home cages were grouped within the same enclosure (i.e. with enclosure ID fitted as a random effect). The significance of individual fixed terms was assessed using likelihood ratio tests comparing models with and without the term of interest. Analyses were performed in JMP v10 and using the lme4 package in R v2.15.1 (R Foundation for Statistical Computing, Vienna, Austria).
This research adhered to the Association for the Study of Animal Behaviour/Animal Behaviour Society Guidelines for the Use of Animals in Research, the legal requirements of the country in which the work was carried out and all institutional guidelines.
Availability of data and materials
The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium  via the PRIDE partner repository with the dataset identifier PXD002900. Morphological data are provided as supplementary information in Additional file 3.
We thank J. Fick, L. Hogg, R. Humphries, S. Jopson, J. F. Lemaître, J. Waters, and A. Davidson for help conducting the experiments, and members of the Mammalian Behaviour & Evolution Group for constructive feedback. This study was funded by research grants from the Leverhulme Trust (grant number F/0025/W) and the Natural Environment Research Council, UK (grant number NE/I013008/1).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Birkhead TR, Hosken DJ, Pitnick S. Sperm biology: an evolutionary perspective. Burlington, MA: Academic Press; 2009.Google Scholar
- Perry JC, Sirot L, Wigby S. The seminal symphony: how to compose an ejaculate. Trends Ecol Evol. 2013;28:414–22.View ArticlePubMedGoogle Scholar
- Dewsbury D. Ejaculate cost and male choice. Am Nat. 1982;119:601–10.View ArticleGoogle Scholar
- Wedell N, Gage MJ, Parker GA. Sperm competition, male prudence and sperm-limited females. Trends Ecol Evol. 2002;17:313–20.View ArticleGoogle Scholar
- Parker GA, Pizzari T. Sperm competition and ejaculate economics. Biol Rev. 2010;85:897–934.View ArticlePubMedGoogle Scholar
- Harcourt AH, Harvey PH, Larson SG, Short RV. Testis weight, body weight and breeding system in primates. Nature. 1981;293:55–7.View ArticlePubMedGoogle Scholar
- Parker GA, Ball MA, Stockley P, Gage MJ. Sperm competition games: a prospective analysis of risk assessment. Proc R Soc B. 1997;264:1793–802.PubMed CentralView ArticlePubMedGoogle Scholar
- Stockley P, Gage MJ, Parker GA, Moller AP. Sperm competition in fishes: the evolution of testis size and ejaculate characteristics. Am Nat. 1997;149:933–54.View ArticlePubMedGoogle Scholar
- Ramm SA, Parker GA, Stockley P. Sperm competition and the evolution of male reproductive anatomy in rodents. Proc R Soc B. 2005;272:949–55.PubMed CentralView ArticlePubMedGoogle Scholar
- Hosken DJ, Ward PI. Experimental evidence for testis size evolution via sperm competition. Ecol Lett. 2001;4:10–3.View ArticleGoogle Scholar
- Pitnick S, Miller GT, Reagan J, Holland B. Males’ evolutionary responses to experimental removal of sexual selection. Proc R Soc B. 2001;268:1071–80.PubMed CentralView ArticlePubMedGoogle Scholar
- Simmons LW, Garcia-Gonzalez F. Evolutionary reduction in testes size and competitive fertilization success in response to the experimental removal of sexual selection in dung beetles. Evolution. 2008;62:2580–91.View ArticlePubMedGoogle Scholar
- Gay L, Hosken D, Vasudev R, Tregenza T, Eady P. Sperm competition and maternal effects differentially influence testis and sperm size in Callosobruchus maculatus. J Evol Biol. 2009;22:1143–50.View ArticlePubMedGoogle Scholar
- Firman RC, Simmons LW. Experimental evolution of sperm competitiveness in a mammal. BMC Evol Biol. 2011;11:19.PubMed CentralView ArticlePubMedGoogle Scholar
- Gage MJ. Continuous variation in reproductive strategy as an adaptive response to population density in the moth Plodia interpunctella. Proc R Soc B. 1995;261:25–30.View ArticleGoogle Scholar
- Bjork A, Dallai R, Pitnick S. Adaptive modulation of sperm production rate in Drosophila bifurca, a species with giant sperm. Biol Lett. 2007;3:517–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Ramm SA, Stockley P. Adaptive plasticity of mammalian sperm production in response to social experience. Proc R Soc B. 2009;276:745–51.PubMed CentralView ArticlePubMedGoogle Scholar
- Janicke T, Marie-Orleach L, De Mulder K, Berezikov E, Ladurner P, Vizoso DB, et al. Sex allocation adjustment to mating group size in a simultaneous hermaphrodite. Evolution. 2013;67:3233–42.View ArticlePubMedGoogle Scholar
- Moatt JP, Dytham C, Thom MD. Sperm production responds to perceived sperm competition risk in male Drosophila melanogaster. Physiol Behav. 2014;131:111–4.View ArticlePubMedGoogle Scholar
- Wigby S, Sirot L, Linklater J, Buehner N, Calboli F, Bretman A, et al. Seminal fluid protein allocation and male reproductive success. Curr Biol. 2009;19:751–7.PubMed CentralView ArticlePubMedGoogle Scholar
- Lemaître JF, Ramm SA, Hurst JL, Stockley P. Social cues of sperm competition influence accessory reproductive gland size in a promiscuous mammal. Proc R Soc B. 2011;278:1171–6.PubMed CentralView ArticlePubMedGoogle Scholar
- Fedorka KM, Winterhalter WE, Ware B. Perceived sperm competition intensity influences seminal fluid protein production prior to courtship and mating. Evolution. 2011;65:584–90.View ArticlePubMedGoogle Scholar
- Harris WE, Moore PJ. Sperm competition and male ejaculate investment in Nauphoeta cinerea: effects of social environment during development. J Evol Biol. 2004;18:474–80.View ArticleGoogle Scholar
- Perry JC, Rowe L. Condition-dependent ejaculate size and composition in a ladybird beetle. Proc R Soc B. 2010;277:3639–47.PubMed CentralView ArticlePubMedGoogle Scholar
- Short RV. Sexual selection and its component parts, somatic and genital selection, as illustrated by Man and the Great Apes. Adv Stud Behav. 1979;9:131–58.View ArticleGoogle Scholar
- Parker GA, Ball MA. Sperm competition, mating rate and the evolution of testis and ejaculate sizes: a population model. Biol Lett. 2005;1:235–8.PubMed CentralView ArticlePubMedGoogle Scholar
- Reuter M, Linklater JR, Lehmann L, Fowler K, Chapman T, Hurst GD. Adaptation to experimental alterations of the operational sex ratio in populations of Drosophila melanogaster. Evolution. 2008;62:401–12.View ArticlePubMedGoogle Scholar
- Vahed K, Parker DJ. The evolution of large testes: sperm competition or male mating rate? Z Tierpsychol. 2012;118:107–17.Google Scholar
- Kenagy GJ, Trombulak SC. Size and function of mammalian testes in relation to body size. J Mammal. 1986;67:1–22.View ArticleGoogle Scholar
- Dean MD, Ardlie KG, Nachman MW. The frequency of multiple paternity suggests that sperm competition is common in house mice (Mus domesticus). Mol Ecol. 2006;15:4141–51.PubMed CentralView ArticlePubMedGoogle Scholar
- Firman RC, Simmons LW. The frequency of multiple paternity predicts variation in testes size among island populations of house mice. J Evol Biol. 2008;21:1524–33.View ArticlePubMedGoogle Scholar
- Firman RC, Klemme I, Simmons LW. Strategic adjustments in sperm production within and between two island populations of house mice. Evolution. 2013;67:3061–70.PubMedGoogle Scholar
- Claydon AJ, Ramm SA, Pennington A, Hurst JL, Stockley P, Beynon RJ. Heterogenous turnover of sperm and seminal vesicle proteins in the mouse revealed by dynamic metabolic labelling. Mol Cell Proteomics. 2012;11:M111.014993. doi:10.1074/mcp.M111.014993.PubMed CentralView ArticlePubMedGoogle Scholar
- Bronson FH. The reproductive ecology of the house mouse. Quart Rev Biol. 1979;54:265–99.View ArticlePubMedGoogle Scholar
- Hurst JL, Beynon RJ. Chemical communication in societies of rodents. In: d’Ettorre P, Hughes D, editors. The sociobiology of communication. Oxford: Oxford University Press; 2008. p. 97–118.View ArticleGoogle Scholar
- delBarco-Trillo J, Ferkin MH. Male mammals respond to a risk of sperm competition conveyed by odours of conspecific males. Nature. 2004;431:446–9.View ArticlePubMedGoogle Scholar
- Kruczek M, Golas A. Behavioural development of conspecific odour preferences in bank voles, Clethrionomys glareolus. Behav Process. 2003;64:31–9.View ArticleGoogle Scholar
- delBarco-Trillo J, LaVenture AB, Johnston RE. Male hamsters discriminate estrous state from vaginal secretions and individuals from flank marks. Behav Process. 2009;82:18–24.View ArticleGoogle Scholar
- Dean MD, Findlay GD, Hoopmann MR, Wu CC, MacCoss MJ, Swanson WJ, et al. Identification of ejaculated proteins in the house mouse (Mus domesticus) via isotopic labeling. BMC Genomics. 2011;12:306.PubMed CentralView ArticlePubMedGoogle Scholar
- McKee CM, Xu D, Kessler BM, Muschel RJ. Proteomic analysis reveals a proteolytic feedback loop in murine seminal fluid. Prostate. 2013;73:1397–440.View ArticleGoogle Scholar
- Mangels R, Young B, Keeble S, Ardekani R, Meslin C, et al. Genetic and phenotypic influences on copulatory plug survival in mice. Heredity. 2015. doi: 10.1038/hdy.2015.50Google Scholar
- Ramm SA, McDonald L, Hurst JL, Beynon RJ, Stockley P. Comparative proteomics reveals evidence for evolutionary diversification of rodent seminal fluid and its functional significance in sperm competition. Mol Biol Evol. 2009;26:189–98.View ArticlePubMedGoogle Scholar
- Sirot LK, Wolfner MF, Wigby S. Protein-specific manipulation of ejaculate composition in response to female mating status in Drosophila melanogaster. Proc Natl Acad Sci U S A. 2011;108:9922–6.PubMed CentralView ArticlePubMedGoogle Scholar
- Kelly CD, Jennions MD. Sexual selection and sperm quantity: meta-analyses of strategic ejaculation. Biol Rev. 2011;86:863–84.View ArticlePubMedGoogle Scholar
- Dean MD. Genetic disruption of the copulatory plug in mice leads to severely reduced fertility. PLoS Genet. 2013;9:e1003185.PubMed CentralView ArticlePubMedGoogle Scholar
- Ramm SA, Oliver PL, Ponting CP, Stockley P, Emes RD. Sexual selection and the adaptive evolution of mammalian ejaculate proteins. Mol Biol Evol. 2008;25:207–19.View ArticlePubMedGoogle Scholar
- Karn RC, Clark NL, Nguyen ED, Swanson WJ. Adaptive evolution in rodent seminal vesicle secretion proteins. Mol Biol Evol. 2008;25:2301–10.PubMed CentralView ArticlePubMedGoogle Scholar
- Ramm SA. Sperm competition and its evolutionary consequences in rodents. UK: University of Liverpool; 2006.Google Scholar
- Hurle B, Swanson W. NISC Comparative Sequencing Program, Green ED. Comparative sequence analyses reveal rapid and divergent evolutionary changes of the WFDC locus in the primate lineage. Genome Res. 2007;17:276–86.PubMed CentralView ArticlePubMedGoogle Scholar
- Kelley LA, Mezulis S, Yates CM, Wass MN, Sternberg MJE. The Phyre2 web portal for protein modeling, prediction and analysis. Nat Protoc. 2015;10:845–58.View ArticlePubMedGoogle Scholar
- Luo CW, Lin HJ, Chen YH. A novel heat-labile phospholipid-binding protein, SVS VII, in mouse seminal vesicle as a sperm motility enhancer. J Biol Chem. 2001;276:6913–21.View ArticlePubMedGoogle Scholar
- Li SH, Lee RKK, Hsiao YL, Chen YH. Demonstration of a glycoprotein derived from the Ceacam10 gene in mouse seminal vesicle secretions. Biol Reprod. 2005;73:546–53.View ArticlePubMedGoogle Scholar
- Morikawa W, Yamamoto K, Ishikawa S, Takemoto S, Ono M, Fukushi JI, et al. Angiostatin generation by cathepsin D secreted by human prostate carcinoma cells. J Biol Chem. 2000;275:38912–20.View ArticlePubMedGoogle Scholar
- Stockley P, Ramm SA, Sherborne AL, Thom MD, Paterson S, Hurst JL. Baculum morphology predicts reproductive success of male house mice under sexual selection. BMC Biol. 2013;11:66.PubMed CentralView ArticlePubMedGoogle Scholar
- Schärer L, Ladurner P, Rieger RM. Bigger testes do work more: experimental evidence that testis size reflects testicular cell proliferation activity in the marine invertebrate, the free-living flatworm Macrostomum sp. Behav Ecol Sociobiol. 2004;56:420–5.View ArticleGoogle Scholar
- Schärer L, Vizoso D. Phenotypic plasticity in sperm production rate: there’s more to it than testis size. Evol Ecol. 2007;21:295–306.View ArticleGoogle Scholar
- Ramm SA, Schärer L. The evolutionary ecology of testicular function: size isn’t everything. Biol Rev. 2014;89:874–88.View ArticlePubMedGoogle Scholar
- Schärer L, Da Lage JL, Joly D. Evolution of testicular architecture in the Drosophilidae: a role for sperm length. BMC Evol Biol. 2008;8:143.PubMed CentralView ArticlePubMedGoogle Scholar
- Lüpold S, Linz GM, Rivers JW, Westneat DF, Birkhead TR. Sperm competition selects beyond relative testes size in birds. Evolution. 2009;63:391–402.View ArticlePubMedGoogle Scholar
- Firman RC, Simmons LW. Experimental evolution of sperm quality via postcopulatory sexual selection in house mice. Evolution. 2010;64:1245–56.PubMedGoogle Scholar
- Firman RC, GarcIa-GonzAlez F, Thyer E, Wheeler S, Yamin Z, Yuan M, et al. Evolutionary change in testes tissue composition among experimental populations of house mice. Evolution. 2015;69:848–55.View ArticlePubMedGoogle Scholar
- Ramm SA, Stockley P. Sperm competition and sperm length influence the rate of mammalian spermatogenesis. Biol Lett. 2010;6:219–21.PubMed CentralView ArticlePubMedGoogle Scholar
- Huber M, Bronson F, Desjardins C. Sexual activity of aged male mice: correlation with level of arousal, physical endurance, pathological status, and ejaculatory capacity. Biol Reprod. 1980;23:305–16.View ArticlePubMedGoogle Scholar
- Preston B, Stockley P. The prospect of sexual competition stimulates premature and repeated ejaculation in a mammal. Curr Biol. 2006;16:R239–41.View ArticlePubMedGoogle Scholar
- Ramm SA, Stockley P. Sequential male mate choice under sperm competition risk. Behav Ecol. 2014;25:660–7.PubMed CentralView ArticlePubMedGoogle Scholar
- Rusu AS, Krackow S. Kin-preferential cooperation, dominance-dependent reproductive skew, and competition for mates in communally nesting female house mice. Behav Ecol Sociobiol. 2004;56:298–305.View ArticleGoogle Scholar
- Wigby S, Chapman T. Female resistance to male harm evolves in response to manipulation of sexual conflict. Evolution. 2004;58:1028–37.View ArticlePubMedGoogle Scholar
- Crudgington HS, Fellows S, Badcock NS, Snook RR. Experimental manipulation of sexual selection promotes greater male mating capacity but does not alter sperm investment. Evolution. 2009;63:926–38.View ArticlePubMedGoogle Scholar
- Simmons LW, Fitzpatrick JL. Sperm wars and the evolution of male fertility. Reproduction. 2012;144:519–34.View ArticlePubMedGoogle Scholar
- Touw WG, Bayjanov JR, Overmars L, Backus L, Boekhorst J, Wels M, et al. Data mining in the Life Sciences with Random Forest: a walk in the park or lost in the jungle? Brief Bioinform. 2013;14:315–26.PubMed CentralView ArticlePubMedGoogle Scholar
- Liaw A, Wiener M. Classification and regression by randomForest. R News. 2002;2:18–22.Google Scholar
- Seung H, Wolfe G, Rocca M. Performing a testicular spermatid head count. Curr Protoc Toxicol. 2003;16:16.7.Google Scholar
- Koyama S, Kamimura S. Lowered sperm motility in subordinate social status of mice. Physiol Behav. 1999;65:665–9.View ArticlePubMedGoogle Scholar
- ImageJ. http://rsbweb.nih.gov/ij
- Wilson-Leedy JG, Ingermann RL. Development of a novel CASA system based on open source software for characterization of zebrafish sperm motility parameters. Theriogenology. 2007;67:661–72.View ArticlePubMedGoogle Scholar
- Vizcaíno JA, Deutsch EW, Wang R, Csordas A, Reisinger F, Ríos D, et al. ProteomeXchange provides globally co-ordinated proteomics data submission and dissemination. Nature Biotechnol. 2014;30:223–6.View ArticleGoogle Scholar