PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
 
Sci Rep. 2017; 7: 11674.
Published online 2017 September 15. doi:  10.1038/s41598-017-12021-2
PMCID: PMC5601457

Differential regulation of vaginal lipocalins (OBP, MUP) during the estrous cycle of the house mouse

Abstract

Female house mice produce pheromone-carrying major urinary proteins (MUPs) in a cycling manner, thus reaching the maximum urinary production just before ovulation. This is thought to occur to advertise the time of ovulation via deposited urine marks. This study aimed to characterize the protein content from the house mouse vaginal flushes to detect putative vaginal-advertising molecules for a direct identification of reproductive states. Here we show that the mouse vaginal discharge contains lipocalins including those from the odorant binding (OBP) and major urinary (MUP) protein families. OBPs were highly expressed but only slightly varied throughout the cycle, whilst several MUPs were differentially abundant. MUP20 or ‘darcin’, was thought to be expressed only by males. However, in females it was significantly up-regulated during estrus similarly as the recently duplicated central/group-B MUPs (sMUP17 and highly expressed sMUP9), which in the mouse urine are male biased. MUPs rise between proestrus and estrus, remain steady throughout metestrus, and are co-expressed with antimicrobial proteins. Thus, we suggest that MUPs and potentially also OBPs are important components of female vaginal advertising of the house mouse.

Introduction

During an estrous cycle, most mammalian females will pass through the four consecutive phases including proestrus, estrus, metestrus, and diestrus (e.g. Mus 1, 2, Apodemus 3, Rattus 4). The estrous cycle in mice is strongly affected by pheromones. Suppression of estrus in single-sex populations5, rapid return to estrus when paired with males6 or exposed to male urine7 are well known effects and they were also observed in other mammals3, 8. The house mouse (Mus musculus) uses a system of lipocalin transporters of active volatile organic compounds–VOCs912, which together serve as signals. In mice, these signals are manifested via expression of large quantities of male-biased13 and highly homologous1416 major urinary proteins (MUPs) in the liver. MUPs are products of a gene cluster on the chromosome 4, that contains ~21 coding genes (and a similar number of pseudogenes), and can be divided into two groups, the group-A (ancestral), containing Mup3, Mup4, Mup5, Mup6, Mup20 and Mup21, and the group-B, consisting of 15 other Mups (i.e. Mup1, Mup2, Mup7-Mup19) sharing almost 99% sequence identity17, 18, reviewed in citations19, 20. They bind some volatile organic compounds (VOCs) in their eight-stranded beta barrel and transport them to the outside environment with urine10, 2123, where they act as an honest, cheat-proof display of an individual’s health and condition24. Such scent-marking signals have strong effects on the reproductive success of the signaller25. Interestingly, group-A and group-B MUPs are not just in the urine, they were recently identified also in tears26 and saliva27. However, most studies focused on the male components of chemical signalling and, thus, there is very little evidence on the roles of female MUPs.

A group of lipocalins that may also be involved in the transport of chemical signals and potentially in the mouse vaginal advertising are members of the odorant binding protein family (OBPs), because in other mammals, they were detected in their vaginal secretions (e.g. hamsters28 and potentially bank voles29) and/or urine (e.g. mole rats30 and bank voles29). The X-linked Obp genes of the mouse were thought to involve just the two nasal members – Obp1a, and Obp1b 31. However, Obp genes have also undergone a series of duplications in mice, and they occur in a cluster of seven genes (i.e. including Prb, Probasin) and two pseudogenes on the X chromosome32, 33. In wild M. m. musculus mice, Obps were predicted from the genome data19, detected as transcripts (Obp1 - KJ605390, Obp2 - KJ605391, Obp5 (Obp1a) - KJ605392, Obp6 - KJ605393, and Obp7 (Obp1b) - KJ605394)33 and their expression was corroborated on the level of protein in the mouse tears26 and saliva27. Mouse OBPs are produced by various oro-facial tissues including olfactory (OE), vomeronasal (VNO), nasal-lymphoid, and lacrimal glands/tissues26, 33 and they are naturally transported with mucosal secretions to the oral cavity, where they are highly abundant27. It has also been demonstrated that they are involved in the process of rapid ligand internalization by means of removing small organic molecules including odorants from nasal mucosa, this has been demonstrated for OBP5 (i.e. OBP1a in olfactory epithelia)34. However, data on these proteins from the mouse vaginal secretions were missing.

In mice, chemical signalling was thought to be manifested via urine marks, and so many experiments have been performed to study in detail this component of chemical signalling. For example, female house mice use major urinary proteins to advertise their estrus with urine marks. It has been shown in the laboratory mice35 and in wild house mice (M. m. musculus)36 that (i) males up-regulate urinary MUPs during social contacts, and that (ii) females use MUPs to advertise their reproductive states by varying the concentration of MUPs in the urine during the estrous cycle35. The urine volatiles during estrus are recognized by a synergistic action of the VNO-OE system of the mouse, and it has been demonstrated that the main olfactory system (MOS) is involved primarily in the attraction from a distance, while the VNO plays a major role in close proximity during pre-copulatory behaviour37. Moreover, there is also a study on rats (Rattus rattus) also showing that a lipocalin with chemosignalling functions was highly expressed during estrus and metestrus phases in the urine38. Thus, it is possible that the close proximity and chemical signals on bodies of interacting individuals could also be important during sexual interactions because this component of signalling is directly linked to an interacting individual, i.e. the signal owner. Thus, the aim of our study was to analyse the relative protein content of vaginal flushes and its variation throughout the cycle to detect potential signalling proteins, and to evaluate the role of particular estrous phases in regulating the expression of these proteins. Furthermore, we also aimed to determine the level of correlation between the house mouse estrous cycle and the host regulation and degradation of bacteria to further the understanding of the emergence of complex odour mixtures.

Materials and Methods

Ethical Standards

All animal procedures were carried out in strict accordance with the law of the Czech Republic paragraph 17 no. 246/1992 and the local ethics committee of the Faculty of Science, Charles University in Prague chaired by Dr. Stanislav Vybíral specifically approved this study in accordance with accreditation no. 27335/2013–17214 valid until 2019.

Subjects, housing conditions and experimental design

In this experiment, we used a total of 9 G1 wild-derived Mus musculus musculus females (120 days old) and 9 unfamiliar males (120 days old) for a period of social stimulation during which vaginal fluid samples were collected on a daily basis. We opted for biological replicates instead of methodological duplicates simply for the fact that in similar studies the duplicates were highly correlated26. Females were housed in pairs with a male and in cages divided by a wire mesh from individual males with square openings (diameter was 1 cm) allowing communication but suppressing direct contact (13:11 hrs, D:N, temperature t = 23 °C). The experimental cages were supplied with fresh bedding at the beginning of the experiment and provided with water and food ad libitum. Vaginal lavage was performed daily between 8:00 and 9:00am by gentle flushing with a pipette using 20 μl of the 0.9% saline solution. Samples were centrifuged at 300 rcf for 10 minutes, 4 °C. Supernatant was used for the protein analysis and the cell pellet was used to prepare slides on a concentrator (StatSpin) for further cytological analyses. These slides were differentially stained to visualize nucleated and cornified cells with May-Gruenwald (3 min) and Giemsa (10 min) staining solutions.

Protein Digestion

Prior to our experiments, we tested the efficiency of ethanol and acetone precipitation of proteins from the vaginal lavage samples. Some acetone precipitated samples did not enter the gel (PAGE) whilst all proteins entered the gel after cold ethanol precipitation. Thus, all protein samples were precipitated with the ice-cold ethanol (20 minutes) and centrifuged at 10 000 rcf for 10 minutes, 0 °C. This was followed by a re-suspension of dried pellets in the digestion buffer (1% SDC, 100 mM TEAB – pH = 8.5). Protein concentration of each lysate was determined using the BCA assay kit (Fisher Scientific). Cysteines in 20 μg of proteins were reduced with a final concentration of 5 mM TCEP (60 °C for 60 min) and blocked with10mM MMTS (i.e. S-methyl methanethiosulfonate, 10 min Room Temperature). Samples were cleaved with trypsin (1 ug of trypsin per sample) in 37 °C overnight. Peptides were desalted on a Michrom C18 column.

nLC-MS2 Analysis

Nano Reversed phase columns were used (EASY-Spray column, 50 cm × 75 µm ID, PepMap C18, 2 µm particles, 100 Å pore size). Mobile phase buffer A was composed of water, and 0.1% formic acid. Mobile phase B contained acetonitrile, and 0.1% formic acid. Samples were loaded onto a trap column (Acclaim PepMap300, C18, 5 µm, 300 Å Wide Pore, 300 µm × 5 mm, 5 Cartridges) for 4 min at 15 μl/min loading buffer was composed of water, 2% acetonitrile and 0.1% trifluoroacetic acid. After 4 minutes ventile was switched and Mobile phase B increased from 4% to 35% B at 60 min, 75% B at 61 min, hold for 8 minutes, and 4% B at 70 min, hold for 15 minutes until the end of run.

Eluting peptide cations were converted to gas-phase ions by electrospray ionization and analysed on a Thermo Orbitrap Fusion (Q-OT-qIT, Thermo). Survey scans of peptide precursors from 350 to 1450 m/z were performed at 120 K resolution (at 200 m/z) with a 5 × 105 ion count target. Tandem MS was performed by isolation at 1.5 Th with the quadrupole, HCD fragmentation with normalized collision energy of 30 and rapid scan MS analysis in the ion trap. The MS2 ion count target was set to 104 and the max injection time was 35 ms. Only those precursors with charge state 2–6 were sampled for MS2. The dynamic exclusion duration was set to 45 s with a 10 ppm tolerance around the selected precursor and its isotopes. Monoisotopic precursor selection was turned on. The instrument was run in top speed mode with 2 s cycles.

Protein analysis

LC-MS data were analysed and quantified with MaxQuant software (version 1.5.3.8)39. The false discovery rate (FDR) was set to 1% for both proteins and peptides and we specified a minimum peptide length of seven amino acids. The Andromeda search engine was used for the MS/MS spectra search against our modified Uniprot Mus musculus database (downloaded on June, 2015), containing 44,900 entries. We modified our databases such that all MUP, OBP sequences were removed and instead of them we have added a complete list of MUPs from Ensembl database, and OBPs from NCBI (sensu - citation33). Next we added some Tremble sequences that were missing in Uniprot, for example KLKs, BPIs, SPINKs, SCGB/ABPs, and LCNs. Enzyme specificity was set as C-terminal to Arg and Lys, also allowing cleavage at proline bonds40 and a maximum of two missed cleavages. Dithiomethylation of cysteine was selected as fixed modification and N-terminal protein acetylation and methionine oxidation as variable modifications. The “match between runs” feature of MaxQuant was used to transfer identifications to other LC-MS/MS runs based on their masses and retention time (maximum deviation 0.7 min) and this was also used in all quantification experiments. Quantifications were performed with the label-free algorithms described recently39 using a combination of unique and razor peptides. All statistical analyses were performed in R software41. First, the dataset was normalized to diminish potential differences due to differential protein extractability and also due to potential differences caused by different signal intensity between samples. We used a normalization based upon quantiles, which normalizes a matrix of peak areas/intensities with the function normalize.quantiles from ‘preprocessCore’ routines under the Bioconductor package42. This method is based upon the concept of a quantile-quantile plot extended to n dimensions. To check that the data distribution conforms to the same type of distribution after normalization, we used ‘mixtools’43. Second, we used the Power Law Global Error Model (PLGEM)44 to detect differentially expressed/abundant proteins using the functions plgem.fit and plgem-stn43. Original and normalized LC-MS/MS data are provided in Dataset 1.

Data validation with targeted quantification of selected proteins

To corroborate evidence on the variation of protein abundances obtained by LC-MS approach, we have selected several ‘key’ proteins (i.e. significant lipocalins that were discussed in the manuscript) and quantified them by Selected Reaction Monitoring (SRM) using triple quadrupole MS/MS instrument TSQ Quantiva (Thermo Scientific). Targeted SRM assays are the mass spectrometry equivalent of a Western blot and are intended to complement discovery-based analysis45. Based on the data from untargeted LC-MS experiment we chose 2 proteotypic peptides per protein for SRM experiment. For each peptide 3 transitions were selected (Dataset 1). LC separation was done with the same setup and conditions as for untargeted LC-MS experiment. Samples were analysed in positive mode with resolution for both Q1 and Q3 set to 0.7, collision gas pressure was set to 1.5 mTorr, cycle time was 1.5 s. Collision energy for each peptide was predicted and results were obtained using Skyline Daily 3.6.1.46. Datasets were normalized the same way as LC-MS data – see above.

Results

The vaginal fluid proteome of the house mouse

We have generated the vaginal fluid proteome of the house mouse, M. m. musculus from 9 females, each in three estrous phases (i.e. proestrus, estrus, metestrus) and detected a total of 2507 proteins at 0.01 FDR (i.e. False Discovery Rate for all peptides and proteins). However, some proteins were detected only in one or two instances per individual. Thus, we reduced our dataset such that all the rows that had just two counts per row were deleted as well as those rows where the median expression in all three groups was lower than 1 (i.e. number of proteins decreased to N = 986). To reduce the influence of false positive abundances due to differences in signal intensities between individuals we quantile-normalized a matrix of protein abundances with ‘preprocessCore’ routine within the Bioconductor package in R software43; this step ensures that differential expression (i.e. abundances or ‘expression patterns’) is measured instead of differential extractability of proteins from complex mucosal secretions. This normalization strategy has resulted in highly similar datasets with similar data distribution, almost the same mean values, and lower standard deviation (i.e. RAW DATA: 19.2 ± 4.9, NORMALIZED DATA: 26.4 ± 0.1, see Fig. 1A) thus decreasing the potential of obtaining false positive values.

Figure 1
Graphical representation of differentially expressed/abundant proteins. Before normalization, the data revealed some variation between individuals (A–yellow bars). However, after the quantile-normalization procedure (A–green bars), the ...

After normalization of our data, many of those proteins that were invariant throughout the cycle and most expressed belong to a group of “housekeeping” genes (HK) which are used in many experiments for normalization. In our data, these included: beta actin (ACTB), serum albumin (ALB), phosphoglycerate kinase 1 (PGK1), fructose-bisphosphate aldolase (ALDOA), Glyceraldehyde 3-phosphate dehydrogenase (G3P, gapdh), and for example all of the 16 detected proteasomal subunits that are provided in Dataset 1. Thus, the lack of any variation in these highly expressed proteins signifies that the quantile-normalization of our data was sufficiently robust and does not require further HK normalization.

Next, we searched for differentially expressed proteins throughout the cycle using the Power Law Global Error Model (PLGEM)44. This model was first developed to quantify microarray data44, however, due to similar statistical properties – namely the n-binomial distributions of signal values (i.e. deviating from normality) – it has proved to be an amenable model for the quantification of label-free MS-based proteomics data47. We calculated the signal-to-noise ratio – STN (equation provided in citation47), because it explicitly takes unequal variances into account and because it penalizes proteins that have higher variance in each class more than those proteins that have a high variance in one class and a low variance in another44. To create statistical baseline, PLGEM can only be fitted on a set of replicates from the same experimental condition, so we have done this for estrus phase. Correlation between the mean values and standard deviations was high (r2 = 0.985, Pearson = 0.977, Supplementary Fig. 1) so we continued with the resampled STNs and calculated differences with corresponding p-values between proestrus and estrus, and between metestrus and estrus. Mean value differences between proestrus and estrus (Fig. 1B), and metestrus and estrus (Fig. 1C) are visualized with MA plots (only protein names with P < 0.05 and fold change FC > 2 or FC < −2 are shown), where in both the plots Fig. 1B,C it is obvious that prevailing number of proteins that are differentially expressed are those that are abundant in estrus (i.e. lower parts of Fig. 1B,C and left parts of the volcano plots in Fig. 1D,E). The plot also demonstrates that the most expressed and at the same time differentially regulated protein is sMUP9 (proestrus→estrus 2 fold, see below).

Differentially expressed proteins – markers of estrus

PLGEM analysis of protein abundances on the level of FC > 2 revealed that a total of 69 proteins (i.e. 7%) were up-regulated in estrus. Of these 69 proteins, a total of 17 proteins (1.7%) were consequently down-regulated in metestrus (Fig. 2A), whilst 52 (5.3%) of them remained unchanged on the level of P < 0.05 and FC > 2. This means that only 1.7% of all proteins were detected as significant markers of estrus. These are: SPA9, LYG1, PHLD, HGFA, ASAH1, HORN, DAG1, CP089, AOC1, CTL4, PPBT, CATC, MA2B2, PPGB, DPP2, FUCO, TGM3. Gene ontology analysis with STRING (https://string-db.org/) revealed that – on the level of FDR = 0.007 – these 17 genes are concurrently involved in catalytic (N = 12 proteins, p = 0.0047) or hydrolase (10 proteins, p = 0.0005) activities (GO functions), and that their functions are localized in extracellular regions (N = 14, p = 8.11e-7) or extracellular exosome (N = 12, p = 8.11e-7) and or lysosome (N = 6, p = 3.67e-5). Thus, the proteins involved in the detected degradation processes – typically occurring during cornification and bacterial invasions during estrus (i.e. a phase, which is characteristic of producing enucleated keratinized cells) – can be treated as markers of this particular phase. For example, Hornerin (HORN) is a component of the epidermal cornified cell envelopes48 and a marker of cell differentiation49. In our data, Hornerin was upregulated in estrus (3.9 fold) and significantly down-regulated in metestrus (3.4 fold). Hornerin is involved in the process of cell differentiation and with mucin/episialin (MUC1)50 forms anti-inflammatory epidermal barriers during estrus. The most differentially expressed protein in estrus (i.e. proestrus→estrus 5.5 fold, estrus→metestrus −5.1 fold) was Serpin A9 (SPA9). SPA9 is a member of the serpin family of serine protease inhibitors.

Figure 2
Graphical representation of individual changes in protein abundances (FC > 2 and P < 0.05) with heat maps. (A)–shows significant proteins abundant only in estrus (i.e. estrus markers, green bar next ...

Differentially expressed anti-microbial proteins

Our analysis (above) revealed that a biological process that is tightly linked to cell proliferation during the estrous cycle involves mechanisms of anti-microbial defence, because the total number of bacteria is highest during estrus and declines during metesturs51. Mucins, particularly MUC1, are effective barriers in the protection from microbial infection50. In our data, MUC1 was significantly up-regulated in estrus (4.5 fold). Similarly, MUC9 (i.e. OVGP1 - Oviduct-specific glycoprotein also known as mucin-9) was also up-regulated during estrus (4 fold) and metestrus. Another protein, which forms viscous barriers in oviducts under low pH during estrus, and which is up-regulated during estrus is Uromodulin (UROM, 2.3 fold). Significantly up-regulated anti-microbial proteins during estrus also involved alpha-macroglobulin or OVOS (Ovostatin, 3 fold) which is slightly down-regulated but not significantly during metestrus. Similar pattern is also typical for the Lysozyme g-like protein 1 (LYG1) which is up-regulated during estrus (2.5 fold) and down-regulated (4.8 fold) during metestrus. Other secretory antimicrobial proteins expressed in polarized epithelial cells involve those that bind to particular structures of bacterial membranes and those that physically break the membranes due to their amphipathic electrostatics. We have detected a member of the PLUNC family – the Lipopolysaccharide-binding protein (LBP, BPIFD2). The saliva proteome of the mouse contains seven members of the bactericidal/permeability-increasing proteins (i.e. BPI52, 53) which are male biased and include BPIA1, BPIB1, BPIB2, BPIB3, BPIFA2, BPIFB5, BPIFB9B27. BPI/PLUNC proteins have an antibacterial activity against gram-negative bacteria52. The vaginal fluid proteome, however, contains only LBP (BPIFD2), which is up-regulated during estrus (4.4 fold) and remains expressed during metestrus.

Anti-microbial proteins, which are significantly up-regulated during metestrus include the Cathelicidin antimicrobial peptide or CAMP/CRAMP (3 fold). A natural antibiotics CAMP forms an amphipathic alpha-helix similar to other antimicrobial peptides, and functional studies showed that CAMP is a potent antibiotics against gram-negative bacteria by inhibiting the growth of a variety of bacterial strains and is expressed by neutrophils and macrophages54. NGP (Neutrophilic granule protein) or ‘bectenecin’ – also belongs to cathelicidins, has a cathelicidin protein domain, and in our data, it is also significantly up-regulated in metestrus (2.8 fold), likely because it is co-expressed with CAMP in neutrophils which invade vaginal environment during metestrus. CAMP is regulated by the serine-proteases Kallikreins 5 and 755 and thus belongs to an extended KLK/LEKTI network members that are crucial for homeostasis of stratified epithelia56. We did not detect KLK5 and KLK7 in the vaginal secretion. However, we have detected KLK1, KLK10, KLK11, KLK12, KLK13, and KLK14 invariantly expressed throughout the cycle, whilst KLK8 varied such that it was up-regulated during proestrus and estrus and down-regulated during metestrus (4.9 fold). SPINK5, which negatively regulates KLK5 expression was down-regulated (2.5 fold) during estrus.

Lipocalins

Previously, it has been shown with microarrays that the mouse (CD 1) uterine transcriptome contains transcripts coding several major urinary proteins, namely Mup3, Mup4, Mup5, Mup6, Mup9, and Mup20 57. In their experiment, all detected Mup transcipts were significantly up-regulated in estrus (i.e. when compared to proestrus). Our data revealed a support for this trend on the level of protein. We have detected MUP20-darcin (unique peptides: VFVEYIHVLENSLALK, FAQLSEEHGIVR), MUP3, MUP5, and also a group of MUP proteins with no unique peptides. These peptides, however, are shared between the recently duplicated group-B MUPs and are, therefore, vizualized in Fig. 2B as the most likely hits though they are mixtures of several MUPs: sMUP17 (i.e. either MUP13 and/or MUP17), sMUP9 (i.e. MUPs 1, 2, 7, 8, 9, 10, 11, 12, 14, 15, 16, 18, 19). From the biological function perspective, it is plausible to group them simply for the fact that most group-B MUPs are highly homologous in M. m. musculus, have similar beta barrel structures, and there is almost no individual variability in the production of these proteins (i.e. at least in the liver/urine) in this sub-species of the house mouse1416. Our proteomic data included also metestrus (i.e. when compared to citation57 where they concentrated on the proestrus/estrus analyses), thus we were interested whether MUPs and other lipocalins are predominantly expressed during estrus and declining in metestrus. A total of three MUPs were up-regulated in estrus: MUP20 (4.2 fold), the highly expressed sMUP9 (2.1 fold), and sMUP17 (3.2 fold) but all of them remained abundant in metestrus.

To provide further evidence on the variation of significant MUPs and LCN11 (see below) and particularly MUP20-darcin, we have selected specific unique-peptide peaks and analysed their variation with Selected Reaction Monitoring (SRM) using triple quadrupole MS/MS instrument. This method is a MS-based equivalent to Western blotting but without antibodies45. We have chosen last three females of the dataset (F7, F7, F9) where two of them (F7, F8) followed the trend from PLGEM analysis whilst one female (F9) had a delay in the rising phase of protein abundance (F9 - green colour in Fig. 3). On the level of p < 0.05, the correlation between LC-MS and SRM data for particular females (differentiated by unique colours in Fig. 3) was high in all three females (note the insets with the correlation coefficients in Fig. 3). Thus, our SRM approach corroborated some of the trends detected with untargeted LC-MS/MS based approach.

Figure 3
Selected Reaction Monitoring (SRM) of unique peptides for MUP20, sMUP9, sMUP17 and LCN11 using the triple quadrupole MS/MS instrument. We used this method as an alternative to Western blotting and revealed that untargeted LC-MS/MS data (first/grey column) ...

For the first time (i.e. to our knowledge), we have detected the four members of the recently described33 family of odorant binding proteins OBP1 (KJ605390), OBP2 (KJ605391), OBP5 (KJ605392), and OBP7 (KJ605394) in vaginal fluids of the house mouse as highly abundant proteins. Our PLGEM analysis revealed significant up-regulations during estrus (e.g. OBP5 1.6 fold, p = 0.0027). However, the fold change is below 2 and it is in question whether these differences are biologically relevant. We have also detected lipocalins that have not been linked to chemical communication before. Lipocalin 11 (LCN11) (i.e. proestrus→estrus 2.3 fold) follows similar and significant trend as MUPs but its total expression is lower. Retinol-binding protein 1(RBP1) had an opposite trend (i.e. proestrus→estrus −2.13, estrus→metestrus 2.8 fold) with the lowest abundances during estrus. Siderocalin (LCN2) belons to a group of the most abundant proteins and was invariant throughout the cycle.

Discussion

Vaginal fluids represent an important source of markers of reproductive states. Proteins obtained from the vaginal secretions by non-invasive lavage methods, and identified with the state-of-the-arts label-free proteomics enabled us to ask questions related to reproduction, chemical communication and evolution but with a detail that may help to support particular hypotheses dealing with the function of signalling proteins. The samples, however, are a mixture of proteins stemming out of multiple processes that take place in ovaries, uterine walls and horns, cervix, and various internal vaginal glands. In this paper, we aimed to characterize the expression pattern of proteins, which are known for their involvement in chemical communication with a particular focus on their variation throughout the estrous cycle of the wild derived house mice Mus musculus musculus. We used wild derived mice for the fact that the laboratory mice are influenced by the differential contribution of blocks of genes from the two house-mouse subspecies M. m. domesticus and M. m. musculus to current laboratory strains58 that may mask their natural behaviour.

We used the cytology screening of the cells from vaginal lavage and selected samples according to predominace of the cell types typical for each phase13. Consequent analyses of the vaginal fluid proteome served to corroborate our cytology determination. We have detected a total of 1.7% of proteins that are prevailingly expressed during estrus, Fig. 2B. They involve markers of cell keratinization as well as markers of catalytic activity and anti-microbial defence. Anti-microbial proteins function as natural innate-immune responses to an outburst of bacteria typically occurring during estrus51 such as Lysozyme G1 (LYG1) which has a hydrolase activity and is involved in the degradation of peptidoglycans from bacterial membranes (i.e. GO:0009253 peptidoglycan catabolic process). Bacterial peptidoglycans are continuously recognized by Peptidoglycan recognition protein 1 (PGRP1) which activates bacterial tool-component systems59 and is logically invariant throughout the cycle in our data.

A marker of interactions between symbiotic bacteria and the host is the highly expressed lipocalin 2/siderocalin (LCN2). In our data, the highly abundant LCN2 was invariantly expressed throughout the cycle, thus demonstrating that bacteria are present in all phases of the cycle. It is known that during metabolic degradation in most mucosal tissues, bacteria attempt to acquire ‘free’ iron by a secretion of high-affinity iron sequestrating siderophores. The mammalian host, however, limits this process by the production of Lipocalin 260 which efficiently scavenges for catecholate-type siderophores61. Thus, high amounts of produced LCN2 represent an efficient regulatory element that prevents uncontrolled bacterial growth also in the mouse vaginal environment. WFDC proteins (i.e. ‘Whey acidic proteins four disulphide core’) were also shown to have anti-microbial properties62 and the two members WFDC12 and WFDC18 are present in the house mouse saliva as proteins encoded by the submandibular gland transcripts Wfdc12 and Wfdc18 27. In this study, we have detected WFDC2 also invariantly expressed throughout the cycle. To add, an antimicrobial system of defence in vaginal environment, thus, seems to have two components, one that is static and invariant over the cycle, and another (dynamic) that reacts to bacterial growth and dynamically functions during estrus (e.g. LYG1, LBP, OVOS) and metestrus (e.g. CAMP, NGP).

Lipocalins of the laboratory mouse were already reported to vary between proestrus and estrus with the use of microarrays of the uterine horns57. They have provided evidence that all the detected lipocalin-coding transcripts, namely Mup3, Mup4, Mup5, Mup6, Mup9, and Mup20 rise from proestrus to estrus. Our protein analysis corroborates their findings to some extent and extends their view because we included also metestrus. We show for the first time in wild-derived M. m. musculus mice that, MUP20 or ‘darcin’ is the most up-regulated lipocalin when passing from proestrus to estrus. This is interesting because some studies originally thought that urinary MUP20 is a male unique lipocalin with pheromonal effects even without ligands and which stimulates an inherent attraction for particular males63. Here we provide evidence, that MUP20 is not male-unique because it is produced by females in their vaginal secretions during estrous cycle (i.e. M. m. musculus) and also by females of the laboratory mice in their uterine horns57. Although, we have no data on what MUP20 binds in the vaginal secretions, we show that MUP20 production significantly rises and reaches its peak in estrus. Interestingly, major urinary proteins were also detected in metestrus which – if they were supposed to signal receptivity – would decrease a male’s chances to correctly evaluate the best time for mating. Furthermore, although male-biased, MUP20 and other MUPs are also produced by salivary27 and lacrimal26 glands of males and females. It is also possible that MUP20 plays different roles in different tissues which may be further explored with behavioural experiments and/or with GC-MS techniques.

There is as yet another alternative explanation to why are some lipocalins expressed also during metestrus. It is possible that MUPs play other roles besides sexual signalling similarly as olfactory receptors play other roles besides the detection of chemical signals64. We already suggested this alternative to a common view that MUPs serve only to chemical communication by providing evidence that MUPs are expressed in a sexually dimorphic manner in one tissue (e.g. liver,13) but are non-dimorphic and with lower expression levels elsewhere (e.g. saliva27 and tears26). Here we suggest a likely hypothesis that lipocalins (i.e. including major urinary proteins and odorant binding proteins) may have the potential to detoxify mucosal tissues by removing various organic compounds. We have postulated this explanation as the ‘Toxic waste hypothesis’19 which states that the same products of metabolic degradation might have been an ideal source of signals that have driven the evolution of chemical communication by means of a consequent toxic-waste perception and recognition. Some of these organic compounds might have come from bacterial degradation because we found various antimicrobial proteins differentially acting during estrous cycle. This hypothesis, however, needs to be further tested with GC-MS techniques. Furthermore, lipocalin involvement in detoxification has been experimentally evidenced, so for example MUPs are already known to transport toxic substances (i.e. other than known pheromones) out of the body65.

To conclude, MUPs and OBPs may function as the transporting devices of volatiles to signal female receptivity in the house mouse, M. m. musculus. The most abundant lipocalins during estrus were MUP20, and the group-B MUPs - sMUP9 and sMUP17. However, they remain expressed during metestrus. Thus, it is likely that these proteins may have other roles besides sexual signalling. We have also detected variation in bactericidal proteins, thus further supporting an idea that females may have adopted alternative strategies in controlling microbiota thus yielding different odor profiles throughout the cycle.

Electronic supplementary material

Dataset 1(183K, pdf)

Acknowledgements

This research was supported by the project BIOCEV (CZ.1.05/1.1.00/02.0109) and National Program for Sustainability II (LQ1604) from the Ministry of Education, Youth and Sports. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. We are also grateful to Helena Uhlířová for her careful and patient lab assistance and to Karel Harant from the Mass Spectrometry and Proteomics Service Laboratory, Faculty of Science, Charles University for performing the LC-MS/MS run.

Author Contributions

Author Contributions

P.S., M.C., and R.S. wrote the main manuscript text, B.K. helped to annotate proteomic datasets. P.T. performed experiments on MS instruments. M.C. had responsibilities over the animal housing and performed animal experiments. All authors have reviewed the manuscript.

Notes

Competing Interests

The authors declare that they have no competing interests.

Footnotes

Electronic supplementary material

Supplementary information accompanies this paper at 10.1038/s41598-017-12021-2.

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

References

1. Cora MC, Kooistra L, Travlos G. Vaginal Cytology of the Laboratory Rat and Mouse: Review and Criteria for the Staging of the Estrous Cycle Using Stained Vaginal Smears. Toxicol Pathol. 2015;43:776–793. doi: 10.1177/0192623315570339. [PubMed] [Cross Ref]
2. Byers SL, Wiles MV, Dunn SL, Taft RA. Mouse estrous cycle identification tool and images. PLoS One. 2012;7:e35538. doi: 10.1371/journal.pone.0035538. [PMC free article] [PubMed] [Cross Ref]
3. Stopka P, Macdonald DW. Signal interchange during mating in the Wood mouse (Apodemus sylvaticus): the concept of active and passive signalling. Behaviour. 1998;135:231–249. doi: 10.1163/156853998793066339. [Cross Ref]
4. Goldman JM, Murr AS, Cooper RL. The rodent estrous cycle: characterization of vaginal cytology and its utility in toxicological studies. Birth Defects Res B Dev Reprod Toxicol. 2007;80:84–97. doi: 10.1002/bdrb.20106. [PubMed] [Cross Ref]
5. Whitten WK. Occurence of anoestrus in mice caged in groups. Journal of Endocrinology. 1959;18:102–107. doi: 10.1677/joe.0.0180102. [PubMed] [Cross Ref]
6. Whitten WK. Modification of the oestrous cycle of the mouse by external stimuli associated with the male. Changes in the oestrous cycle determined by vaginal smears. J. Endocrinol. 1958;17:307–313. doi: 10.1677/joe.0.0170307. [PubMed] [Cross Ref]
7. Marsden HM, Bronson FH. Estrus Synchrony in Mice: Alteration by Exposure to Male Urine. Science. 1965;144:1469. doi: 10.1126/science.144.3625.1469. [PubMed] [Cross Ref]
8. Bronson FH, Marsden HM. Male-Induced synchrony of estrus in deermice. General and comparative endocrinology. 1964;4:634–637. doi: 10.1016/0016-6480(64)90073-5. [PubMed] [Cross Ref]
9. Sharrow SD, Novotny MV, Stone MJ. Thermodynamic analysis of binding between mouse major urinary protein-I and the pheromone 2-sec-butyl-4,5-dihydrothiazole. Biochemistry. 2003;42:6302–6309. doi: 10.1021/bi026423q. [PubMed] [Cross Ref]
10. Sharrow SD, Vaughn JL, Žídek L, Novotny MV, Stone MJ. Pheromone binding by polymorphic mouse major urinary proteins. Protein Science. 2002;11:2247–2256. doi: 10.1110/ps.0204202. [PubMed] [Cross Ref]
11. Novotny MV, et al. A unique urinary constituent, 6-hydroxy-6-methyl-3-heptanone, is a pheromone that accelerates puberty in female mice. Chem. Biol. 1999;6:377–383. doi: 10.1016/S1074-5521(99)80049-0. [PubMed] [Cross Ref]
12. Zidek L, et al. NMR Mapping of the Recombinant Mouse Major Urinary Protein I Binding site Occupied by the Pheromone 2-sec-Butyl-4,5-dihydrothiazole. Biochemistry. 1999;38:9850–9861. doi: 10.1021/bi990497t. [PubMed] [Cross Ref]
13. Stopková R, Stopka P, Janotová K, Jedelsky PL. Species-specific expression of major urinary proteins in the house mice (Mus musculus musculus and Mus musculus domesticus) J Chem Ecol. 2007;33:861–869. doi: 10.1007/s10886-007-9262-9. [PubMed] [Cross Ref]
14. Enk VM, et al. Regulation of highly homologous major urinary proteins in house mice quantified with label-free proteomic methods. Mol Biosyst. 2016;12:3005–3016. doi: 10.1039/C6MB00278A. [PMC free article] [PubMed] [Cross Ref]
15. Thoss M, et al. Diversity of major urinary proteins (MUPs) in wild house mice. Sci Rep. 2016;6:38378. doi: 10.1038/srep38378. [PMC free article] [PubMed] [Cross Ref]
16. Thoß, M., Luzynski, K., Ante, M., Miller, I. & Penn, D. J. Major urinary protein (MUP) profiles show dynamic changes rather than individual ‘barcode’ signatures. Frontiers in Ecology and Evolution3, doi:10.3389/fevo.2015.00071 (2015). [PMC free article] [PubMed]
17. Logan, D. W., Marton, T. F. & Stowers, L. Species Specificity in Major Urinary Proteins by Parallel Evolution. PLoS ONE3, doi: 10.1371/journal.pone.0003280 (2008). [PMC free article] [PubMed]
18. Mudge JM, et al. Dynamic instability of the major urinary protein gene family revealed by genomic and phenotypic comparisons between C57 and 129 strain mice. Genome Biol. 2008;9:R91. doi: 10.1186/gb-2008-9-5-r91. [PMC free article] [PubMed] [Cross Ref]
19. Stopková, R., Hladovcová, D., J., K., Vyoral, D. & Stopka, P. Multiple roles of secretory lipocalins (MUP, OBP) in mice. Folia Zool. 58, 29–40 (2009).
20. Janotová K, Stopka P. Mechanisms of chemical communication: the role of Major Urinary Proteins. Folia Zool. 2009;58:41–55.
21. Kwak J, et al. Changes in volatile compounds of mouse urine as it ages: their interactions with water and urinary proteins. Physiol Behav. 2013;120:211–219. doi: 10.1016/j.physbeh.2013.08.011. [PubMed] [Cross Ref]
22. Timm DE, Baker LJ, Mueller H, Zidek L, Novotny MV. Structural basis of pheromone binding to mouse major urinary protein (MUP-I) Protein Science. 2001;10:997–1004. doi: 10.1110/ps.52201. [PubMed] [Cross Ref]
23. Phelan MM, et al. The structure, stability and pheromone binding of the male mouse protein sex pheromone darcin. PLoS One. 2014;9:e108415. doi: 10.1371/journal.pone.0108415. [PMC free article] [PubMed] [Cross Ref]
24. Zala SM, Potts WK, Penn DJ. Scent-marking displays provide honest signals of health and infection. Behavioral Ecology. 2004;15:338–344. doi: 10.1093/beheco/arh022. [Cross Ref]
25. Thonhauser KE, Raveh S, Hettyey A, Beissmann H, Penn DJ. Scent marking increases male reproductive success in wild house mice. Anim Behav. 2013;86:1013–1021. doi: 10.1016/j.anbehav.2013.09.004. [PMC free article] [PubMed] [Cross Ref]
26. Stopkova, R., Klempt, P., Kuntova, B. & Stopka, P. On the tear proteome of the house mouse (Mus musculus musculus) in relation to chemical signalling. PeerJ 6, e3541, doi:10.7717/peerj.3541 (2017). [PMC free article] [PubMed]
27. Stopka P, et al. On the saliva proteome of the Eastern European house mouse (Mus musculus musculus) focusing on sexual signalling and immunity. Sci Rep. 2016;6:32481. doi: 10.1038/srep32481. [PMC free article] [PubMed] [Cross Ref]
28. Singer A, Macrides F, clancy AN, Agosta WC. Purification and Analysis of a Proteinaceous Aphrodisiac Pheromone from Hamster Vaginal Discharge. Journal of Biological Chemistry. 1986;261:13323–13326. [PubMed]
29. Stopkova R, et al. Novel OBP genes similar to hamster Aphrodisin in the bank vole, Myodes glareolus. BMC Genomics. 2010;11:45. doi: 10.1186/1471-2164-11-45. [PMC free article] [PubMed] [Cross Ref]
30. Hagemeyer P, et al. Searching for major urinary proteins (MUPs) as chemosignals in urine of subterranean rodents. J Chem Ecol. 2011;37:687–694. doi: 10.1007/s10886-011-9971-y. [PubMed] [Cross Ref]
31. Pes D, Dal Monte M, Ganni M, Pelosi P. Isolation of two odorant-binding proteins from mouse nasal tissue. Comp. Biochem. Physiol. 1992;103B:1011–1017. [PubMed]
32. Stopkova R, Dudkova B, Hajkova P, Stopka P. Complementary roles of mouse lipocalins in chemical communication and immunity. Biochem Soc T. 2014;42:893–898. doi: 10.1042/BST20140053. [PubMed] [Cross Ref]
33. Stopkova, R. et al. Mouse lipocalins (MUP, OBP, LCN) are co-expressed in tissues involved in chemical communication. Frontiers in Ecology and Evolution4, doi:10.3389/fevo.2016.00047 (2016).
34. Strotmann J, Breer H. Internalization of odorant-binding proteins into the mouse olfactory epithelium. Histochemistry and cell biology. 2011;136:357–369. doi: 10.1007/s00418-011-0850-y. [PubMed] [Cross Ref]
35. Stopka P, Janotova K, Heyrovsky D. The advertisement role of major urinary proteins in mice. Physiology & Behavior. 2007;91:667–670. doi: 10.1016/j.physbeh.2007.03.030. [PubMed] [Cross Ref]
36. Janotova K, Stopka P. The level of major urinary proteins is socially regulated in wild Mus musculus musculus. J Chem Ecol. 2011;37:647–656. doi: 10.1007/s10886-011-9966-8. [PubMed] [Cross Ref]
37. Achiraman S, Ponmanickam P, Ganesh DS, Archunan G. Detection of estrus by male mice: synergistic role of olfactory-vomeronasal system. Neurosci Lett. 2010;477:144–148. doi: 10.1016/j.neulet.2010.04.051. [PubMed] [Cross Ref]
38. Muthukumar S, et al. Urinary lipocalin protein in a female rodent with correlation to phases in the estrous cycle: an experimental study accompanied by in silico analysis. PLoS One. 2013;8:e71357. doi: 10.1371/journal.pone.0071357. [PMC free article] [PubMed] [Cross Ref]
39. Cox J, et al. Accurate proteome-wide label-free quantification by delayed normalization and maximal peptide ratio extraction, termed MaxLFQ. Mol Cell Proteomics. 2014;13:2513–2526. doi: 10.1074/mcp.M113.031591. [PMC free article] [PubMed] [Cross Ref]
40. Rodriguez J, Gupta N, Smith RD, Pevzner PA. Does trypsin cut before proline? J Proteome Res. 2008;7:300–305. doi: 10.1021/pr0705035. [PubMed] [Cross Ref]
41. Crawley, M. J. The R Book. (Wiley Publishing, 2007).
42. Bolstad BM, Irizarry RA, Astrand M, Speed TP. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003;19:185–193. doi: 10.1093/bioinformatics/19.2.185. [PubMed] [Cross Ref]
43. Gentleman RC, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5:R80. doi: 10.1186/gb-2004-5-10-r80. [PMC free article] [PubMed] [Cross Ref]
44. Pavelka N, et al. A power law global error model for the identification of differentially expressed genes in microarray data. BMC Bioinformatics. 2004;5:203. doi: 10.1186/1471-2105-5-203. [PMC free article] [PubMed] [Cross Ref]
45. Arnott D, et al. Selective detection of membrane proteins without antibodies: a mass spectrometric version of the Western blot. Mol Cell Proteomics. 2002;1:148–156. doi: 10.1074/mcp.M100027-MCP200. [PubMed] [Cross Ref]
46. MacLean B, et al. Skyline: an open source document editor for creating and analyzing targeted proteomics experiments. Bioinformatics. 2010;26:966–968. doi: 10.1093/bioinformatics/btq054. [PMC free article] [PubMed] [Cross Ref]
47. Pavelka N, et al. Statistical similarities between transcriptomics and quantitative shotgun proteomics data. Mol Cell Proteomics. 2008;7:631–644. doi: 10.1074/mcp.M700240-MCP200. [PubMed] [Cross Ref]
48. Henry J, et al. Hornerin is a component of the epidermal cornified cell envelopes. FASEB J. 2011;25:1567–1576. doi: 10.1096/fj.10-168658. [PubMed] [Cross Ref]
49. Makino T, Takaishi M, Morohashi M, Huh NH. Hornerin, a novel profilaggrin-like protein and differentiation-specific marker isolated from mouse skin. J Biol Chem. 2001;276:47445–47452. doi: 10.1074/jbc.M107512200. [PubMed] [Cross Ref]
50. DeSouza MM, et al. MUC1/episialin: a critical barrier in the female reproductive tract. Journal of Reproductive Immunology. 1999;45:127–158. doi: 10.1016/S0165-0378(99)00046-7. [PubMed] [Cross Ref]
51. Noguchi K, Tsukumi K, Urano T. Qualitative and quantitative differences in normal vaginal flora of conventionally reared mice, rats, hamsters, rabbits, and dogs. Comp Med. 2003;53:404–412. [PubMed]
52. LeClair EE. Four reasons to consider a novel class of innate immune molecules in the oral epithelium. J Dent Res. 2003;82:944–950. doi: 10.1177/154405910308201202. [PubMed] [Cross Ref]
53. Leclair, E. E. Four BPI (bactericidal/permeability-increasing protein)-like genes expressed in the mouse nasal, oral, airway and digestive epithelia. Biochem Soc Trans31, 801–805, 10.1042/ (2003). [PubMed]
54. Gallo RL, et al. Identification of CRAMP, a cathelin-related antimicrobial peptide expressed in the embryonic and adult mouse. J Biol Chem. 1997;272:13088–13093. doi: 10.1074/jbc.272.20.13088. [PubMed] [Cross Ref]
55. Yamasaki K, et al. Kallikrein-mediated proteolysis regulates the antimicrobial effects of cathelicidins in skin. FASEB J. 2006;20:2068–2080. doi: 10.1096/fj.06-6075com. [PubMed] [Cross Ref]
56. Kasparek P, et al. KLK5 and KLK7 Ablation Fully Rescues Lethality of Netherton Syndrome-Like Phenotype. PLoS Genet. 2017;13:e1006566. doi: 10.1371/journal.pgen.1006566. [PMC free article] [PubMed] [Cross Ref]
57. Yip KS, Suvorov A, Connerney J, Lodato NJ, Waxman DJ. Changes in mouse uterine transcriptome in estrus and proestrus. Biol Reprod. 2013;89:13. doi: 10.1095/biolreprod.112.107334. [PMC free article] [PubMed] [Cross Ref]
58. Mouse Genome Sequencing C, et al. Initial sequencing and comparative analysis of the mouse genome. Nature. 2002;420:520–562. doi: 10.1038/nature01262. [PubMed] [Cross Ref]
59. Royet J, Gupta D, Dziarski R. Peptidoglycan recognition proteins: modulators of the microbiome and inflammation. Nat Rev Immunol. 2011;11:837–851. [PubMed]
60. Goetz DH, et al. The Neutrophil Lipocalin NGAL Is a Bacteriostatic Agent that Interferes with Siderophore-Mediated Iron Acquisition. Molecular Cell. 2002;10:1033–1043. doi: 10.1016/S1097-2765(02)00708-6. [PubMed] [Cross Ref]
61. Flo TH, et al. Lipocalin 2 mediates an innate immune response to bacterial infection by sequestrating iron. Nature. 2004;432:917–921. doi: 10.1038/nature03104. [PubMed] [Cross Ref]
62. Scott A, Weldon S, Taggart CC. SLPI and elafin: multifunctional antiproteases of the WFDC family. Biochem Soc Trans. 2011;39:1437–1440. doi: 10.1042/BST0391437. [PubMed] [Cross Ref]
63. Roberts, S. A. et al. Darcin: a male pheromone that stimulates female memory and sexual attraction to an individual male’s odour. BMC Biol. 8, 10.1186/1741-7007-1188-1175. (2010). [PMC free article] [PubMed]
64. Ferrer I, et al. Olfactory Receptors in Non-Chemosensory Organs: The Nervous System in Health and Disease. Front Aging Neurosci. 2016;8:163. doi: 10.3389/fnagi.2016.00163. [PMC free article] [PubMed] [Cross Ref]
65. Kwak J, Strasser E, Luzynski K, Thoss M, Penn DJ. Are MUPs a Toxic Waste Disposal System? PLoS One. 2016;11:e0151474. doi: 10.1371/journal.pone.0151474. [PMC free article] [PubMed] [Cross Ref]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group