Genomic variation in the vomeronasal receptor gene repertoires of inbred mice
© Wynn et al.; licensee BioMed Central Ltd. 2012
Received: 2 February 2012
Accepted: 18 July 2012
Published: 21 August 2012
Vomeronasal receptors (VRs), expressed in sensory neurons of the vomeronasal organ, are thought to bind pheromones and mediate innate behaviours. The mouse reference genome has over 360 functional VRs arranged in highly homologous clusters, but the vast majority are of unknown function. Differences in these receptors within and between closely related species of mice are likely to underpin a range of behavioural responses. To investigate these differences, we interrogated the VR gene repertoire from 17 inbred strains of mice using massively parallel sequencing.
Approximately half of the 6222 VR genes that we investigated could be successfully resolved, and those that were unambiguously mapped resulted in an extremely accurate dataset. Collectively VRs have over twice the coding sequence variation of the genome average; but we identify striking non-random distribution of these variants within and between genes, clusters, clades and functional classes of VRs. We show that functional VR gene repertoires differ considerably between different Mus subspecies and species, suggesting these receptors may play a role in mediating behavioural adaptations. Finally, we provide evidence that widely-used, highly inbred laboratory-derived strains have a greatly reduced, but not entirely redundant capacity for differential pheromone-mediated behaviours.
Together our results suggest that the unusually variable VR repertoires of mice have a significant role in encoding differences in olfactory-mediated responses and behaviours. Our dataset has expanded over nine fold the known number of mouse VR alleles, and will enable mechanistic analyses into the genetics of innate behavioural differences in mice.
KeywordsVomeronasal Receptor Olfaction Pheromone Behaviour Genome sequencing Single nucleotide polymorphism Mouse
The clone-based sequencing and assembly of the first high quality mouse genome, from the inbred C57BL/6J strain, revealed over 20,000 protein coding genes [1, 2]. Approximately three quarters of these are direct 1:1 orthologues of genes identified in the human genome, but a number of gene families show striking rodent-specific expansions. These are enriched in genes associated with reproduction, reflecting an influence of sexual competition on the evolution of the mouse genome . Mouse sexual behaviour and physiology is strongly regulated by pheromones: biochemical signals emitted by conspecifics that directly influence the behaviour or physiology of a receiving animal . Accordingly, families of genes encoding protein pheromones are among those expanded in rodents [4–6]. This is mirrored by extensive species-specific expansions (with some selective losses and pseudogenisations) in gene families encoding rodent vomeronasal receptors (VRs) [1, 7–12], which are putative receptors for pheromones [13, 14].
The rodent VR repertoire consists of three distantly related families of G-protein coupled receptor, each expressed in sensory neurons of the vomeronasal organ in the nose: V1Rs, V2Rs and Formyl-peptide receptors (FPRs) [15–20]. Each vomeronasal sensory neuron (VSN) expresses a very restricted sub-set of VRs, typically one or a few of the same sub-family, thereby patterning each neuron to detect a limited number of ligands [15, 19–23]. With a few exceptions [14, 24, 25], direct relationships between specific VRs and their pheromone ligands remain unresolved. However some general trends have emerged: neurons expressing V1Rs appear tuned to detect small volatile chemicals [26, 27], those expressing V2Rs have been shown to detect proteinaceous ligands [13, 28, 29], while those expressing FPRs are stimulated by disease and inflammation-related peptides [19, 20]. In addition, a significant subset VSNs detect chemical signals emitted by sympatric and predatory species, suggesting some VRs may be specifically tuned to allo-specific, instead of conspecific, cues [30–32].
The number of VRs identified in the C57BL/6J genome steadily increased with each draft genome assembly release, but has stabilised as the genome nears completion [10, 11]. Recent analyses find 239 functional V1R genes , 121 functional V2R genes  and 7 FPR genes (of which 5 are expressed in the VNO [19, 20]) distributed across the mouse genome in tightly clustered arrays. When their sequences are compared, the VRs within each cluster tend to fall within phylogenetically related clades, indicating that they expanded by localised duplication events [8, 9]. Comparisons of VR clades across a range of mammals show that each species has a “semi-private” repertoire, consistent with a functional role as receptors for species-specific signals such as pheromones [7–11, 23].
Although the variation in VR repertoire between divergent species is well documented , much less is known about the microevolution of VR repertoires between very closely related species. An analysis of 18 V1Rs between M. m. musculus and M. spretus revealed dynamic modulation in evolutionary pressures, including examples of positive selection and lineage-specific pseudogenisation . Similarly, a recent study of 44 V1Rs in M. m. musculus and M. m. domesticus found evidence of genes shaped by negative selection and random drift, with a small proportion having evidence of positive selection . Thus different rodent VR genes may be evolving under very different selective pressures, perhaps depending on the functional nature of their ligands; though whether these correlate within VR families, clades or even clusters is unknown.
Our understanding of intra-specific VR variation is even more limited. On one hand VRs may be expected to be under stabilising selection within a species, favouring low variation in receptor characteristics to ensure essential chemo signals (predator cues or sex-specific signals, for example [13, 30, 31]) are accurately detected and transduced. On the other hand, sexual selection driven by individual recognition may favour variation between individuals in both chemo signals and receptors . Genome wide studies of copy number variation have reported significant variation in the size of VR clusters within inbred laboratory strains of mice [36, 37], though accurate quantification of VR content and sequence level analyses of specific receptors are lacking. However, the value of using highly inbred, domesticated lab mice to study pheromone signalling has also been questioned . Pheromone signalling may be largely redundant in laboratory mice releasing VR genes from selective constraint; alternatively these strains may offer a snapshot of the VR variation within the founder population, frozen by generations of inbreeding . It is also possible that the very process of laboratory domestication has artificially shaped pheromone receptor content or variation, as was recently demonstrated in C. elegans.
The development of massively parallel sequencing (MPS) now offers the opportunity to rapidly and affordably resequence genomes from multiple individuals of the same species . Recently the Mouse Genomes Project produced high coverage genome sequence of 17 inbred strains of mouse, including several laboratory-derived strains and wild-derived M. musculus musculus, M. musculus domesticus and M. musculus castaneus strains [41–43]. These are widely considered to be three subspecies of Mus musculus. M. m. musculus ranges across eastern Europe and northern Asia, domesticus is found throughout western Europe, the Middle East and north Africa, whereas castaneus extends throughout south and east Asia (reviewed in ). The laboratory-derived mouse strains are overwhelmingly domesticus in origin, with minor genomic contributions from castaneus x(~5%) and musculus (0.3%) . The Mouse Genomes Project also sequenced the genome of a wild-derived Mus spretus strain [41–43]. M. spretus is a distinct species of short-tailed mouse with a sympatric range to domesticus, and serves as a genetic outgroup to the Mus musculus subspecies.
Here we interrogate the full complement of coding VR gene sequences in these 17 strains, compared to the C57BL/6J reference. We show that MPS can accurately resolve over half of mouse VRs, but that copy number variation (CNV) and non-specific short read mapping confounds complete repertoire analysis. We identify divergent patterns in the number and distribution of single-nucleotide polymorphisms (SNPs) within genes, clades and clusters of VRs, as well as between strains, and present evidence that the functional vomeronasal repertoire may vary significantly between Mus species and subspecies.
Generating an accurate dataset of VR SNP variation
Variation in vomeronasal receptor genes identified by massively parallel sequencing
VRs after parsing
Resolved sequence (in Kb)
Unresolved sequence %
Non-synonymous SNPs (%)
Private SNPs (%)
To assess the accuracy of our dataset, we compared VR sequences identified by MPS with those generated by traditional long-read, capillary sequencing. Recently Kurzweil et al. reported the sequencing of a cluster of V1R genes and pseudogenes from a SPRET/EiJ BAC library . We identified nine BAC sequenced V1R genes (8.1 kb of sequence) that was also in our parsed dataset. Assuming no error in the BAC sequencing, we calculate a false positive SNP call rate of 1.3%, with a false negative rate of 0.5%. This equates to 99.95% accuracy in base-pair calling. We were also able to identify four genes (3.6 kb of sequence) that we had excluded for failing to meet our quality threshold. In these genes the false-positive and false-negative SNP rates are significantly higher, 5% and 17.5% respectively, which supports the use of our parsing strategy to generate an accurate, albeit conservative, SNP catalogue for this study.
VR SNP distribution among mouse strains
Functional variation in VRs
V1Rs and FPRs are members of the class A rhodopsin-like receptor superfamily . They are encoded by small, typically single exon, genes and dominated by the core seventrans-membrane domains with only a small extracellular domain . V2Rs belong to the class C glutamate-like receptor superfamily; they are multi-exonic and are characterised by a long N-terminal extracellular domain. This region is predicted to be involved in ligand recognition and, in contrast to the highly conserved seven trans-membrane domains, is highly variable between different members of the V2R family (Figure2B) [9, 49]. We divided 1220 V2R genes from 16 strains into N-terminal and trans-membrane domains, and calculated the ratio of non-synonymous to synonymous SNPs for each. In both laboratory-derived and wild-derived strains, we found a statistically significant bias in functional variation towards the N-terminal domains (Figure2C). This is consistent with an enhanced selective constraint on the conserved seven trans-membrane domains and suggests that V2Rs, at least, are not entirely redundant in domesticated mice.
SNP distribution between and within VR families
As we found no evidence of differences between entire VR sub-families we next tested whether there was statistical correlation in variation within sub-families, based on the hypothesis that closely related receptors are likely to be activated by closely related ligands that may, in turn, regulate similar types of behaviours. Considering wild-derived and lab-derived strains independently, principal component analyses (PCA, an unsupervised linear feature extraction method that discovers the directions of maximal variances in data) found highly significant correlations between non-synonymous SNP distribution and phylogeny within lab strain V2Rs (Two-way ANOVA, variance by V2R clade F12,858 = 17.99, P < 0.0001). For example, over 90% of the V2R variation in lab-derived strains can be represented by four principal components (eigenvalues <2). Grouping of VRs by cluster is apparent by the first two components in PCA, as shown by V2RB and V2RC, which are highly conserved among lab-derived strains, and V2RA5, which is highly variable (Figure3, Figure4B). Similar correlations can be observed within V1Rs and wild-derived strains (Two-way ANOVA, variance by V2R clade F12,264 = 5.52, P < 0.0001) and even within closely related receptor clades, such as the two non-contiguous VR clusters that together generate V2RA8, Vmn2r57-61 on chromosome 7and Vmn2r90-110 on chromosome 17 (Figure3, Figure4C). Physically clustered VRs may therefore act as “units” of variation, with constituent receptor genes sharing characteristics of SNP accumulation or resistance. This implies that at least some clusters of VRs were shaped by contrasting selective pressures before or during strain derivation, but can these differences be reconciled with what is currently known about VR function?
SNP distribution between functional classes of V2R
Recently Isogai et al. undertook a comprehensive survey of immediate early gene upregulation in VSNs after exposure to olfactory stimuli derived from conspecifics and a diverse range of other animals. They then identified the VRs expressed in these VSNs and showed that individual V2Rs appear specifically tuned to detect the nature of the signal emitter, whereas V1Rs are activated by cues from multiple species . Although both the specific chemosignals and the behavioural consequences of detection remain unresolved, these groups of V2Rs represent the first functional classification of VRs in mice. The receptors can be classed as those that detect conspecific odour cues (soiled bedding from the same species of mouse), predators (including birds of prey, snakes and mammals) and sympatric non-predatory rodents (including M. spicilegus and M. m. musculus) . We reasoned that these functional classes of V2R would likely be under divergent selective pressures within and between species/subspecies of mouse. Indeed, within the strains that represent different species/subspecies, the V2Rs that detect sympatric cues differ significantly more from the C57BL/6J reference than those that detect either conspecific (Two-way ANOVA, Bonferroni post hoc test, P < 0.001) or predator cues (P < 0.01) (Figure4D). However, within domesticus strains there is no significant difference in variation by function (P > 0.05).
Massively parallel sequencing of VRs
To our knowledge, this is the first application of MPS to specifically resequence and analyse vomeronasal receptor genes. Therefore we first assessed whether current technology was capable of accurately resolving the sequences of these highly homologous, clustered sequences. Across entire genomes, 77-87% of reference sequence can be resolved, depending on the strain . In contrast we were able to resolve between 50-60% of known coding VR sequence with high accuracy. While some of the unresolved sequence is certainly a consequence of strain-specific VR deletions or duplications, the majority is likely due to the challenge of uniquely mapping short sequence reads to highly repetitive genes . Two strands of evidence support this conclusion: firstly the VR clades with the highest sequence homology are enriched in the unresolved sequence, and these are also typically unresolved across all strains. Secondly in C57BL/6NJ, a sub-line of the C57BL/6J reference strain that has just one SNP difference over 300 kb of VR sequence, only 60% resolution was achieved. It is improbable that such a large proportion of the VR repertoire would differ between such closely related sub-strains while virtually no SNPs were accumulated. Instead it is likely that 60% resolution is a baseline value for mapping reads to VRs in this dataset. In comparison, we were able to resolve only 50% of the spretus VRs, which may indicate that as much as one tenth of the domesticus functional VR repertoire is missing in spretus (for an example, see Additional file 3: Figure S1). MPS is therefore currently only partially effective for accurately resequencing VRs in mouse. However, as mapping success is largely influenced by read length, the rapid development of MPS technology suggests that more of the VR repertoire will be resolved in future resequencing efforts .
Considerations in using inbred strains for behavioural research
Domestication by hundreds of generations of inbreeding has resulted in strains of laboratory mouse that display constrained innate behavioural responses when compared to wild (or recently wild-derived) relatives [52–54]. We find that variation within VRs to be similarly restricted within lab-derived strains, reinforcing the genomic limitation of these strains in studies of olfactory-mediated mouse behaviour, especially those that involve individual recognition or the detection of genetic relatedness by olfaction [35, 38, 55]. Nevertheless, inbred lab-derived strains do differ from one another in some olfactory and/or pheromone mediated behaviours [14, 56, 57]. Strain-specific variance in genetically-encoded pheromone production will certainly account for some of these differences [5, 14, 28, 58], but we also find evidence of unusual strain-specific differences in VRs that may influence behaviour. For example, in some cases unusual clustering of rare privates SNP scan be observed. 129P2/OlaHsd has no private SNPs in V1Rs but 8% of the V2R SNPs in that strain are private and almost all are clustered within a single VR clade, V2RA5. Furthermore 129P2/OlaHsd, which serves as the background strain for many knockout mice lines , also has over 5% of the C57BL/6J functional VR repertoire either truncated, frame-shifted or deleted. It is important to stress that the truncating SNPs and frame-shifting indels described here do not definitively denote pseudogenisation, as alternative splicing or termination codons near the 3′ end of genes may still result in fully functional receptors being produced. Nor can we confidently report this the upper limit of VR pseudogenisation, as other SNPs that result in non-conservative amino-acid substitutions may also result in functionally disabled genes and are not recorded here. Nevertheless, these data will assist strain selection and control when conducting or comparing experiments on olfactory-mediated behaviour in lab mice. The VR variation reported here may also explain conflicting reports in the literature about the effectiveness of some putative pheromones .
Selection on VRs in inbred mice
We found that the distribution of SNPs in lab-derived strains is not random, but correlated with VR clade and genomic clusters, suggesting at least some VRs remain under selective constraint during the process of domestication by humans. It is tempting to speculate that those clades with the least variation mediate critical pheromone-mediated behaviours, such as reproduction or maternal care, which are essential even in a laboratory environment. However the paucity of VRs with defined ligands or functions makes it impossible to draw firm conclusions at this time. Of the VRs with known ligands [14, 20, 24, 25], Vmn2r116 (also known as V2rp5) was not resolved in this study. Two others, Vmn2r26 (also known as V2r1b) and Vmn1r49 (also known as V1rb2), are highly conserved among lab-derived strains. The FPRs display a range of variation, but it is presently unclear what behaviours these receptors mediate. It is apparent, however, that the receptors in the atypical V2RC clade, which may act as functional co-receptors with other V2Rs , are among the most conserved VRs.
As a corollary to selection for essential VRs, receptors that are widely inactivated may be either disadvantageous (or inconsequential) to domesticated environments. Recently a study demonstrated examples of acquired pheromone resistance during domestication of C. elegans, by deletion of dauer pheromone receptor genes . Similar pressures may have acted on VRs during mouse domestication, especially as a number of pheromone-mediated behaviours, such as aggression or inbreeding avoidance, are likely to be disadvantageous in a laboratory environment. Unfortunately, because the reference genome is also from a laboratory-derived inbred strain, this study is not well suited to detect such deletions. However mapping the VR sequences described here to a de novo genome assembly of a wild-derived domesticus mouse genome would permit the identification of novel VRs missing in domesticated mice.
VRs between species and subspecies of mice
We used four wild-derived strains as representatives of “wild” species and subspecies of Mus: M. musculus musculus, M. musculus domesticus M. musculus castaneus and M. spretus. These strains are also inbred to homozygosity, albeit through fewer generations, and therefore under-represent the true genomic variation in wild-caught mice. Nevertheless here we have shown that they faithfully capture many of the pheromone receptor gene sequences found in wild individuals , supporting their use as genetically controlled, “pseudo-wild” animals to investigate VR function. However, some wild-derived inbred strains have genomic contamination with classical strains  and our data reveals specific examples of this when compared with sub-sets of VRs from wild-caught mice . Most strikingly, the entire V1RA and V1RB receptor cluster in the musculus wild-derived strain (PWK/PhJ) are identical to the C57BL/6J reference (see Figure3A) and do not match any of the alleles from wild musculus individuals (see Figure 1D). Even with such cases of introgression between subspecies, we were able to identify significant differences in VR variation depending on the nature of their ligand. Most species/subspecies of Mus will encounter olfactory signals from similar predatory species that activate subsets of VRs resulting in defensive behaviours [30–32]. There is relatively little coding variation in those receptors between domesticus and non-domesticus derived strains, consistent with strong selective pressure. In contrast, some VRs are tuned to detect cues from other sympatric mouse subspecies [30, 32]. The functions of these receptors are unknown, but they may mediate a behavioural barrier that contributes to reproductive isolation between subspecies . We found that, collectively, these receptors are more variable in non-domesticus derived strains than those detecting predators or conspecifics. Vmn1r67 (also known as V1re10) has been proposed as a candidate for mediating subspecific recognition between these subspecies . The gene was not included in the analysis described above, but we did identify the same 24 SNP differences between musculus and domesticus alleles reported previously. Unexpectedly, four laboratory-derived strains (A/J, BALB/cJ, NOD/ShiLtJ and NZO/HlLtJ) also share the presumably introgressed musculus allele while the others have the domesticus allele (with one additional SNP). Thus it may now be possible to exploit these introgressions to delineate a function for this unusually variable receptor in subspecific recognition, using more commonly available laboratory-derived mice.
Here we have resolved, analysed and compared the full coding sequences of over 3300 pheromone receptors of unknown function, from 17 strains of inbred mice using massively parallel sequencing data of whole genomes. We describe complex patterns of non-random sequence variation that indicate these receptors are under divergent selective pressures that correlate with proposed ligand, phylogeny, chromosomal clustering, and protein domain. Moreover, we demonstrate that VR protein sequences are unusually variable between mice, and that the repertoire of functional genes differs significantly, mirroring the inter-individual diversity in pheromone ligands. Together this supports a two-component model in which differences in both receptor and pheromone ligand sequence (and expression) may together genetically encode a diverse range of innate responses. We anticipate, as VRs are increasingly deorphanised in mice, the data described here will be invaluable for investigation into the genetic basis of behavioural differences.
We retrieved VR cDNA sequences from 239 V1R , 121 V2R , and 7 FPR sequences [19, 20], and use the nomenclature of those authors. We used BLAT to determine the genomic co-ordinates of each gene exon in the NCBIM37 mouse reference genome. One gene, Vmn2r119, was excluded because it mapped to the same location as Vmn2r118. We then used Mouse Genomes Project SNP query tool to compare the genome sequences of 17 strains (European Nucleotide Archive: [ENA: ERP000034–ERP000050]) across each exon of each VR gene, relative to the NCBIM37 reference . Full detail of how each strain was sequenced and mapped, and how the SNPs were called is described by Keane et al.. Briefly, frame shifts between 1 and 50 bp were called by Dindel . SNPs were called by agreement of two or more of four different SNP calling algorithms: Samtools varFilter 8, Genome Analysis Toolkit 10, iMR and QCALL 11. The SNP query tool parsed each SNP with a quality score >10, and then calculated the predicted transcriptional consequence from its location in an open reading frame: reporting synonymous, non-synonymous or ambiguous SNPs, and additionally when a termination codon is gained or lost. We then used Look-seq with a MAPQ cutoff >30 to identify and exclude genes lacking sequence coverage across the entire coding region . These included clusters of genes with extremely high homology (e.g. entire V1RD and V2RA4 clades), and genes that appear to have specific deletions in only some lineages (e.g. see Additional file 3: Figure S1). Only the latter are described as deletions in Figures 2E and 3. Genes that had > 50% ambiguous SNP calls were also excluded. Although these ambiguous calls are experimentally unresolved, they are typically the result of multiple haplotypes in the sequenced strain mapping to the same location in the reference strain (e.g. see Additional file 4: Figure S2), which is consistent with a duplication of that locus. These are described as duplications in Figures 2E and 3. Finally, if an orthologous gene was excluded in more than three strains for any of the above reasons, SNPs from that gene were then excluded from all other strains for statistical purposes. Additional file 1: Table S1 shows the number of genes excluded for each reason. All the SNPs in the dataset used in this study are listed in Additional files 56 and 7. The mouse strains used in this study are: 129P2/OlaHsd, 129S1/SvImJ, 129S5SvEvBrd, A/J, AKR/J, BALB/cJ, C3H/HeJ, C57BL/6NJ, CAST/EiJ, CBA/J, DBA/2J, LP/J, NOD/ShiLtJ, NZO/HlLtJ, PWK/PhJ, SPRET/EiJ and WSB/EiJ.
SNP density heat maps were generated using Cluster and Java Treeview [66, 67], and arranged using the phylogeny reconstructed using MEGA from aligned C57BL/6J VR cDNA sequences using the neighbour-joining method with the Kimura-2 parameter model of substitution . The heat-mapping of wild-caught isolates was arranged by hierarchical clustering with complete linkage.
TMHMM was used to predict the location of the beginning of the first trans-membrane domain in each V2R . The number of synonymous and non-synonymous SNPs either side of this location was calculated and normalised to the domain size. SNPs located in the short loops between each trans-membrane span were included in the seven trans-membrane domain count. The relative location of truncation codons were also identified using TMHMM.
Comparisons with published sequence data
We aligned the V1R sequences reported by Kurzweil et al. to a de novo scaffold of the SPRET/EiJ genome  using CLUSTAL W2  and identified each gene by sequence identity and synteny. The nine gene pairs compared were YUA. 5 (Vmn1r40), YUA.3 (Vmn1r42), YUB.2 (Vmn1r43), YUA.4 (Vmn1r44), YUB.1ps (Vmn1r46), YUC.1 (Vmn1r49), YUC.3 (Vmn1r50), YUC.5 (Vmn1r51) and YUD.4 (Vmn1r54). V1R sequences from wild isolates [GenBank: JF782602–JF783819] were compared to the orthologous gene in each strain genome. Five genes were excluded because either no sequence could be generated in some isolates  or they did not meet the sequence quality threshold in our dataset. The wild mouse dataset contains two VR sequences per isolate (Park, S.H. and Zhang, J., personal communication), these were paired together and a match was called if at least one was identical to a strain haplotype. When multiple strains had an identical haplotype, the wild-derived strain of the same subspecies was called ahead of the laboratory derived strains or the wild-derived strain from a difference subspecies.
Total counts of synonymous and non-synonymous SNPs across the entire 17 genomes are from the data described in  and further personal communication from those authors. The number of non-synonymous SNPs across all 17 genomes total 136,968 in lab-derived strains and 178,126 in wild-derived strains, and the number of synonymous SNPs total 253,181 in lab-derived strains and 361,993 in wild-derived strains. The total number of SNPs in VRs is detailed in Table1. The SNPs per codon densities were calculated from 599,688,770 bp of coding sequence across the entire 17 genomes (Goodstadt, L., personal communication), and 5, 236, 930 bp of VR coding sequence (Table1).
V2R functional categorisation
V2Rs that respond to conspecific, predator or sympatric cues were identified from . Briefly, Isogai et al. exposed domesticus-derived mice to olfactory stimuli from each source and then quantified the overlap of Egr1 positive cells in the VNO (indicating neuronal activation) with in situ hybridization of probes against VRs. Specific probes to 29 V2R genes in our parsed dataset co localized with cells responsive to conspecific, predator or sympatric cues. Two V2Rs were excluded because they were responsive to more than one class of stimulus, leaving 11 conspecific, 12 predator and 4 sympatric responsive V2Rs. The coding variation in these genes was compared within 13 domesticus-derived strains and the non-domesticus strains: CAST/EiJ, PWK/PhJ and SPRET/EiJ.
Two-way ANOVA with repeated measures for strain followed by Bonferroni post hoc tests, or two-tailed t-tests, were applied as appropriate using GraphPad Prism 5. Statistical significance was considered when P < 0.05. PCA was performed using non-synonymous SNPs / kb values from 13 lab-derived strains using the publicly available R package for statistical computing, version 2.14.0 . Ten missing values were encountered due to lineage specific deletions/duplications. The average of the two neighbouring genes was used in these cases.
Copy number variation
Formyl-peptide receptor family
Massively parallel sequencing
Principal component analysis
Single nucleotide polymorphism
Vomeronasal receptor family 1
Vomeronasal receptor family 2
Vomeronasal sensory neuron.
We are grateful to David Adams, Thomas Keane and Leo Goodstadt of the Mouse Genomes Project for sharing pre-publication sequence data and unpublished analysis of the 17 genomes, Carl Anderson for providing statistical advice, Ian Jackson for advice on mouse strains, Maria Levitin for technical support, and Rachel Larder and Ximena Ibarra-Soria for helpful discussion. This work was supported by the Wellcome Trust, grant number 098051. DWL is a Wellcome Trust Career Development Fellow and KJC is the recipient of a Wellcome Trust PhD studentship.
- Church DM, Goodstadt L, Hillier LW, Zody MC, Goldstein S, She X, Bult CJ, Agarwala R, Cherry JL, DiCuccio M, et al: Lineage-specific biology revealed by a finished genome assembly of the mouse. PLoS Biol. 2009, 7: e1000112-10.1371/journal.pbio.1000112.PubMed CentralView ArticlePubMedGoogle Scholar
- Waterston RH, Lindblad-Toh K, Birney E, Rogers J, Abril JF, Agarwal P, Agarwala R, Ainscough R, Alexandersson M, An P, et al: Initial sequencing and comparative analysis of the mouse genome. Nature. 2002, 420: 520-562. 10.1038/nature01262.View ArticlePubMedGoogle Scholar
- Tirindelli R, Dibattista M, Pifferi S, Menini A: From pheromones to behavior. Physiol Rev. 2009, 89: 921-956. 10.1152/physrev.00037.2008.View ArticlePubMedGoogle Scholar
- Logan DW, Marton TF, Stowers L: Species specificity in major urinary proteins by parallel evolution. PLoS One. 2008, 3: e3280-10.1371/journal.pone.0003280.PubMed CentralView ArticlePubMedGoogle Scholar
- Kimoto H, Sato K, Nodari F, Haga S, Holy TE, Touhara K: Sex- and strain-specific expression and vomeronasal activity of mouse ESP family peptides. Curr Biol. 2007, 17: 1879-1884. 10.1016/j.cub.2007.09.042.View ArticlePubMedGoogle Scholar
- Karn RC, Laukaitis CM: The mechanism of expansion and the volatility it created in three pheromone gene clusters in the mouse (Mus musculus) genome. Genome Biol Evol. 2009, 1: 494-503.PubMed CentralView ArticlePubMedGoogle Scholar
- Lane RP, Young J, Newman T, Trask BJ: Species specificity in rodent pheromone receptor repertoires. Genome Res. 2004, 14: 603-608. 10.1101/gr.2117004.PubMed CentralView ArticlePubMedGoogle Scholar
- Shi P, Bielawski JP, Yang H, Zhang YP: Adaptive diversification of vomeronasal receptor 1 genes in rodents. J Mol Evol. 2005, 60: 566-576. 10.1007/s00239-004-0172-y.View ArticlePubMedGoogle Scholar
- Yang H, Shi P, Zhang YP, Zhang J: Composition and evolution of the V2r vomeronasal receptor gene repertoire in mice and rats. Genomics. 2005, 86: 306-315. 10.1016/j.ygeno.2005.05.012.View ArticlePubMedGoogle Scholar
- Young JM, Massa HF, Hsu L, Trask BJ: Extreme variability among mammalian V1R gene families. Genome Res. 2010, 20: 10-18. 10.1101/gr.098913.109.PubMed CentralView ArticlePubMedGoogle Scholar
- Young JM, Trask BJ: V2R gene families degenerated in primates, dog and cow, but expanded in opossum. Trends Genet. 2007, 23: 212-215. 10.1016/j.tig.2007.03.004.View ArticlePubMedGoogle Scholar
- Yang H, Shi P: Molecular and evolutionary analyses of formyl peptide receptors suggest the absence of VNO-specific FPRs in primates. J Genet Genomics. 2010, 37: 771-778. 10.1016/S1673-8527(09)60094-1.View ArticlePubMedGoogle Scholar
- Chamero P, Marton TF, Logan DW, Flanagan K, Cruz JR, Saghatelian A, Cravatt BF, Stowers L: Identification of protein pheromones that promote aggressive behaviour. Nature. 2007, 450: 899-902. 10.1038/nature05997.View ArticlePubMedGoogle Scholar
- Haga S, Hattori T, Sato T, Sato K, Matsuda S, Kobayakawa R, Sakano H, Yoshihara Y, Kikusui T, Touhara K: The male mouse pheromone ESP1 enhances female sexual receptive behaviour through a specific vomeronasal receptor. Nature. 2010, 466: 118-122. 10.1038/nature09142.View ArticlePubMedGoogle Scholar
- Dulac C, Axel R: A novel family of genes encoding putative pheromone receptors in mammals. Cell. 1995, 83: 195-206. 10.1016/0092-8674(95)90161-2.View ArticlePubMedGoogle Scholar
- Herrada G, Dulac C: A novel family of putative pheromone receptors in mammals with a topographically organized and sexually dimorphic distribution. Cell. 1997, 90: 763-773. 10.1016/S0092-8674(00)80536-X.View ArticlePubMedGoogle Scholar
- Matsunami H, Buck LB: A multigene family encoding a diverse array of putative pheromone receptors in mammals. Cell. 1997, 90: 775-784. 10.1016/S0092-8674(00)80537-1.View ArticlePubMedGoogle Scholar
- Ryba NJ, Tirindelli R: A new multigene family of putative pheromone receptors. Neuron. 1997, 19: 371-379. 10.1016/S0896-6273(00)80946-0.View ArticlePubMedGoogle Scholar
- Liberles SD, Horowitz LF, Kuang D, Contos JJ, Wilson KL, Siltberg-Liberles J, Liberles DA, Buck LB: Formyl peptide receptors are candidate chemosensory receptors in the vomeronasal organ. Proc Natl Acad Sci USA. 2009, 106: 9842-9847. 10.1073/pnas.0904464106.PubMed CentralView ArticlePubMedGoogle Scholar
- Riviere S, Challet L, Fluegge D, Spehr M, Rodriguez I: Formyl peptide receptor-like proteins are a novel family of vomeronasal chemosensors. Nature. 2009, 459: 574-577. 10.1038/nature08029.View ArticlePubMedGoogle Scholar
- Martini S, Silvotti L, Shirazi A, Ryba NJ, Tirindelli R: Co-expression of putative pheromone receptors in the sensory neurons of the vomeronasal organ. J Neurosci. 2001, 21: 843-848.PubMedGoogle Scholar
- Silvotti L, Moiani A, Gatti R, Tirindelli R: Combinatorial co-expression of pheromone receptors, V2Rs. J Neurochem. 2007, 103: 1753-1763. 10.1111/j.1471-4159.2007.04877.x.View ArticlePubMedGoogle Scholar
- Silvotti L, Cavalca E, Gatti R, Percudani R, Tirindelli R: A recent class of chemosensory neurons developed in mouse and rat. PLoS One. 2011, 6: e24462-10.1371/journal.pone.0024462.PubMed CentralView ArticlePubMedGoogle Scholar
- Leinders-Zufall T, Ishii T, Mombaerts P, Zufall F, Boehm T: Structural requirements for the activation of vomeronasal sensory neurons by MHC peptides. Nat Neurosci. 2009, 12: 1551-1558. 10.1038/nn.2452.View ArticlePubMedGoogle Scholar
- Boschat C, Pelofi C, Randin O, Roppolo D, Luscher C, Broillet MC, Rodriguez I: Pheromone detection mediated by a V1r vomeronasal receptor. Nat Neurosci. 2002, 5: 1261-1262. 10.1038/nn978.View ArticlePubMedGoogle Scholar
- Del Punta K, Leinders-Zufall T, Rodriguez I, Jukam D, Wysocki CJ, Ogawa S, Zufall F, Mombaerts P: Deficient pheromone responses in mice lacking a cluster of vomeronasal receptor genes. Nature. 2002, 419: 70-74. 10.1038/nature00955.View ArticlePubMedGoogle Scholar
- Leinders-Zufall T, Lane AP, Puche AC, Ma W, Novotny MV, Shipley MT, Zufall F: Ultrasensitive pheromone detection by mammalian vomeronasal neurons. Nature. 2000, 405: 792-796. 10.1038/35015572.View ArticlePubMedGoogle Scholar
- Leinders-Zufall T, Brennan P, Widmayer P, PC S, Maul-Pavicic A, Jager M, Li XH, Breer H, Zufall F, Boehm T: MHC class I peptides as chemosensory signals in the vomeronasal organ. Science. 2004, 306: 1033-1037. 10.1126/science.1102818.View ArticlePubMedGoogle Scholar
- Kimoto H, Haga S, Sato K, Touhara K: Sex-specific peptides from exocrine glands stimulate mouse vomeronasal sensory neurons. Nature. 2005, 437: 898-901. 10.1038/nature04033.View ArticlePubMedGoogle Scholar
- Ben-Shaul Y, Katz LC, Mooney R, Dulac C: In vivo vomeronasal stimulation reveals sensory encoding of non specific and allospecific cues by the mouse accessory olfactory bulb. Proc Natl Acad Sci USA. 2010, 107: 5172-5177. 10.1073/pnas.0915147107.PubMed CentralView ArticlePubMedGoogle Scholar
- Papes F, Logan DW, Stowers L: The vomeronasal organ mediates interspecies defensive behaviors through detection of protein pheromone homologs. Cell. 2010, 141: 692-703. 10.1016/j.cell.2010.03.037.PubMed CentralView ArticlePubMedGoogle Scholar
- Isogai Y, Si S, Pont-Lezica L, Tan T, Kapoor V, Murthy VN, Dulac C: Molecular organization of vomeronasal chemoreception. Nature. 2011Google Scholar
- Kurzweil VC, Getman M, Green ED, Lane RP: Dynamic evolution of V1R putative pheromone receptors between Mus musculus and Mus spretus. BMC Genomics. 2009, 10: 74-10.1186/1471-2164-10-74.PubMed CentralView ArticlePubMedGoogle Scholar
- Park SH, Podlaha O, Grus WE, Zhang J: The microevolution of V1r vomeronasal receptor genes in mice. Genome Biol Evol. 2011, 3: 401-412. 10.1093/gbe/evr039.PubMed CentralView ArticlePubMedGoogle Scholar
- Hurst JL, Payne CE, Nevison CM, Marie AD, Humphries RE, Robertson DH, Cavaggioni A, Beynon RJ: Individual recognition in mice mediated by major urinary proteins. Nature. 2001, 414: 631-634. 10.1038/414631a.View ArticlePubMedGoogle Scholar
- Cutler G, Marshall LA, Chin N, Baribault H, Kassner PD: Significant gene content variation characterizes the genomes of inbred mouse strains. Genome Res. 2007, 17: 1743-1754. 10.1101/gr.6754607.PubMed CentralView ArticlePubMedGoogle Scholar
- She X, Cheng Z, Zollner S, Church DM, Eichler EE: Mouse segmental duplication and copy number variation. Nat Genet. 2008, 40: 909-914. 10.1038/ng.172.PubMed CentralView ArticlePubMedGoogle Scholar
- Cheetham SA, Smith AL, Armstrong SD, Beynon RJ, Hurst JL: Limited variation in the major urinary proteins of laboratory mice. Physiol Behav. 2009, 96: 253-261. 10.1016/j.physbeh.2008.10.005.View ArticlePubMedGoogle Scholar
- McGrath PT, Xu Y, Ailion M, Garrison JL, Butcher RA, Bargmann CI: Parallel evolution of domesticated Caenorhabditis species targets pheromone receptor genes. Nature. 2011, 477: 321-325. 10.1038/nature10378.PubMed CentralView ArticlePubMedGoogle Scholar
- van der Weyden L, White JK, Adams DJ, Logan DW: The mouse genetics toolkit: revealing function and mechanism. Genome Biol. 2011, 12: 224-10.1186/gb-2011-12-6-224.PubMed CentralView ArticlePubMedGoogle Scholar
- Mouse Genomes Project. http://www.sanger.ac.uk/mousegenomes/
- Keane TM, Goodstadt L, Danecek P, White MA, Wong K, Yalcin B, Heger A, Agam A, Slater G, Goodson M, et al: Mouse genomic variation and its effect on phenotypes and gene regulation. Nature. 2011, 477: 289-294. 10.1038/nature10413.PubMed CentralView ArticlePubMedGoogle Scholar
- Yalcin B, Wong K, Agam A, Goodson M, Keane TM, Gan X, Nellaker C, Goodstadt L, Nicod J, Bhomra A, et al: Sequence-based characterization of structural variation in the mouse genome. Nature. 2011, 477: 326-329. 10.1038/nature10432.PubMed CentralView ArticlePubMedGoogle Scholar
- Guenet JL, Bonhomme F: Wild mice: an ever-increasing contribution to a popular mammalian model. Trends Genet. 2003, 19: 24-31. 10.1016/S0168-9525(02)00007-0.View ArticlePubMedGoogle Scholar
- Yang H, Wang JR, Didion JP, Buus RJ, Bell TA, Welsh CE, Bonhomme F, Yu AH, Nachman MW, Pialek J, et al: Subspecific origin and haplotype diversity in the laboratory mouse. Nat Genet. 2011, 43: 648-655. 10.1038/ng.847.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang H, Bell TA, Churchill GA, Pardo-Manuel De Villena F: On the subspecific origin of the laboratory mouse. Nat Genet. 2007, 39: 1100-1107. 10.1038/ng2087.View ArticlePubMedGoogle Scholar
- Nordstrom KJ, Sallman Almen M, Edstam MM, Fredriksson R, Schioth HB: Independent hhsearch, needleman-wunsch-based, and motif analyses reveal the overall hierarchy for most of the g protein-coupled receptor families. Mol Biol Evol. 2011, 28: 2471-2480. 10.1093/molbev/msr061.View ArticlePubMedGoogle Scholar
- Bryson-Richardson RJ, Logan DW, Currie PD, Jackson IJ: Large-scale analysis of gene structure in rhodopsin-like GPCRs: evidence for widespread loss of an ancient intron. Gene. 2004, 338: 15-23. 10.1016/j.gene.2004.05.001.View ArticlePubMedGoogle Scholar
- Rondard P, Goudet C, Kniazeff J, Pin JP, Prezeau L: The complexity of their activation mechanism opens new possibilities for the modulation of mGlu and GABAB class C G protein-coupled receptors. Neuropharmacology. 2011, 60: 82-92. 10.1016/j.neuropharm.2010.08.009.View ArticlePubMedGoogle Scholar
- Trapnell C, Salzberg SL: How to map billions of short reads onto genomes. Nat Biotechnol. 2009, 27: 455-457. 10.1038/nbt0509-455.PubMed CentralView ArticlePubMedGoogle Scholar
- Schatz MC, Delcher AL, Salzberg SL: Assembly of large genomes using second-generation sequencing. Genome Res. 2010, 20: 1165-1173. 10.1101/gr.101360.109.PubMed CentralView ArticlePubMedGoogle Scholar
- Koide T, Ikeda K, Ogasawara M, Shiroishi T, Moriwaki K, Takahashi A: A new twist on behavioral genetics by incorporating wild-derived mouse strains. Exp Anim. 2011, 60: 347-354. 10.1538/expanim.60.347.View ArticlePubMedGoogle Scholar
- Le Roy I, Roubertoux PL, Jamot L, Maarouf F, Tordjman S, Mortaud S, Blanchard C, Martin B, Guillot PV, Duquenne V: Neuronal and behavioral differences between Mus musculus domesticus (C57BL/6JBy) and Mus musculus castaneus (CAST/Ei). Behav Brain Res. 1998, 95: 135-142. 10.1016/S0166-4328(97)00218-0.View ArticlePubMedGoogle Scholar
- Blanchard DC, Griebel G, Blanchard RJ: Mouse defensive behaviors: pharmacological and behavioral assays for anxiety and panic. Neurosci Biobehav Rev. 2001, 25: 205-218. 10.1016/S0149-7634(01)00009-4.View ArticlePubMedGoogle Scholar
- Sherborne AL, Thom MD, Paterson S, Jury F, Ollier WE, Stockley P, Beynon RJ, Hurst JL: The genetic basis of inbreeding avoidance in house mice. Curr Biol. 2007, 17: 2061-2066. 10.1016/j.cub.2007.10.041.PubMed CentralView ArticlePubMedGoogle Scholar
- Priebe K, Romeo RD, Francis DD, Sisti HM, Mueller A, McEwen BS, Brake WG: Maternal influences on adult stress and anxiety-like behavior in C57BL/6J and BALB/cJ mice: a cross-fostering study. Dev Psychobiol. 2005, 47: 398-407. 10.1002/dev.20098.View ArticlePubMedGoogle Scholar
- Jones SE, Brain PF: Performances of inbred and outbred laboratory mice in putative tests of aggression. Behav Genet. 1987, 17: 87-96. 10.1007/BF01066013.View ArticlePubMedGoogle Scholar
- Roberts SA, Simpson DM, Armstrong SD, Davidson AJ, Robertson DH, McLean L, Beynon RJ, Hurst JL: Darcin: a male pheromone that stimulates female memory and sexual attraction to an individual male’s odour. BMC Biol. 2010, 8: 75-10.1186/1741-7007-8-75.PubMed CentralView ArticlePubMedGoogle Scholar
- Flanagan KA, Webb W, Stowers L: Analysis of male pheromones that accelerate female reproductive organ development. PLoS One. 2011, 6: e16660-10.1371/journal.pone.0016660.PubMed CentralView ArticlePubMedGoogle Scholar
- Ishii T, Mombaerts P: Coordinated coexpression of two vomeronasal receptor V2R genes per neuron in the mouse. Mol Cell Neurosci. 2011, 46: 397-408. 10.1016/j.mcn.2010.11.002.View ArticlePubMedGoogle Scholar
- Vošlajerová BB, Macholan M, Baird SJ, Munclinger P, Dufkova P, Laukaitis CM, Karn RC, Luzynski K, Tucker PK, Pialek J: Reinforcement selection acting on the European house mouse hybrid zone. Mol Ecol. 2011, 20: 2403-2424. 10.1111/j.1365-294X.2011.05106.x.View ArticleGoogle Scholar
- Karn RC, Young JM, Laukaitis CM: A candidate subspecies discrimination system involving a vomeronasal receptor gene with different alleles fixed in M. m. domesticus and M. m. musculus. PLoS One. 2010, 5: e12638-10.1371/journal.pone.0012638.PubMed CentralView ArticlePubMedGoogle Scholar
- Mouse Genomes Project SNP query tool. http://www.sanger.ac.uk/cgi-bin/modelorgs/mousegenomes/snps.pl
- Albers CA, Lunter G, MacArthur DG, McVean G, Ouwehand WH, Durbin R: Dindel: accurate indel calls from short-read data. Genome Res. 2011, 21: 961-973. 10.1101/gr.112326.110.PubMed CentralView ArticlePubMedGoogle Scholar
- Manske HM, Kwiatkowski DP: LookSeq: a browser-based viewer for deep sequencing data. Genome Res. 2009, 19: 2125-2132. 10.1101/gr.093443.109.PubMed CentralView ArticlePubMedGoogle Scholar
- Saldanha AJ: Java Treeview–extensible visualization of microarray data. Bioinformatics. 2004, 20: 3246-3248. 10.1093/bioinformatics/bth349.View ArticlePubMedGoogle Scholar
- de Hoon MJ, Imoto S, Nolan J, Miyano S: Open source clustering software. Bioinformatics. 2004, 20: 1453-1454. 10.1093/bioinformatics/bth078.View ArticlePubMedGoogle Scholar
- Kumar S, Nei M, Dudley J, Tamura K: MEGA: a biologist-centric software for evolutionary analysis of DNA and protein sequences. Brief Bioinform. 2008, 9: 299-306. 10.1093/bib/bbn017.PubMed CentralView ArticlePubMedGoogle Scholar
- Krogh A, Larsson B, von Heijne G, Sonnhammer EL: Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol. 2001, 305: 567-580. 10.1006/jmbi.2000.4315.View ArticlePubMedGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, et al: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23: 2947-2948. 10.1093/bioinformatics/btm404.View ArticlePubMedGoogle Scholar
- The R Project for Statistical Computing. http://www.r-project.org/
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.