BACKGROUND
adaptation leaves its mark in the genome.
these genomic signatures of adaptation can be investigated using phylogenetic and population genetic approaches, which address the different evolutionary timescales of interspecific and intraspecific variation.
a rapidly growing body of literature is documenting the loci involved in adaptation at one of these levels  <cit> .
however, there are few studies which directly evaluate the relationship between the two, or separate studies examining the same loci at different timescales .
hence it is an open question whether positive selection that occurs at loci among species is generally reflected in positive selection at the same loci within species.
different selection regimes make different predictions for the relationship between past and ongoing selection.
in cases where there is an ongoing co-evolutionary dynamic, such as in host-parasite or red queen systems, selection on key loci may be more or less continuous and a signal of selection is expected both in the past and present.
under balancing selection, for example, multiple alleles may be under on-going frequency-dependent selection across speciation events, and similar signals of selection can therefore be expected to act from the past to the present.
another possibility is that directional selection during a certain evolutionary period leads to the fixation of novel and advantageous variants.
after this first and potentially episodic period of adaptation, purifying selection may then follow to stabilize and maintain the acquired adaptations.
under these conditions, there will be discordance between inferences of past and ongoing selection.
many other patterns are also possible: for example, changes in the effective population size  will affect fixation of mildly deleterious alleles by drift and selection on compensatory alleles, so changes in ne alone may lead to different signatures of selection in different timeframes  <cit> .
in mammals, g-protein coupled chemoreceptors form one of the gene families that are most consistently found to be under positive selection in phylogenetic analyses  <cit> .
one of the most diverse classes of chemoreceptors in terrestrial mammals is vomeronasal receptors  that function in intraspecific communication and predator detection.
there is evidence for positive selection on both types of vr  in mammals  <cit> .
this is consistent with the possibility of arms races between the scent of predators and the vrs of their prey, or the action of sexual selection on vr evolution for intraspecific communication.
mouse lemurs  are nocturnal strepsirrhine primates with an elaborated olfactory behavioural repertoire .
olfactory stimuli play an important role in their intraspecific communication as well as in predator detection .
they have the largest repertoire of vr genes  of any primate.
we previously showed a high level of positive selection acting on multiple families of v1r genes, as well as v2r genes, during mouse lemur evolution  <cit> .
we hypothesised that this positive selection may relate to predator detection and/or intraspecific communication.
if this is the case then one would also expect current selection on mouse lemur vr loci.
to the best of our knowledge, population genetics of vrs has not previously been investigated in the wild.
functional loci have been relatively neglected in studies of demography, which favour neutral markers .
the phylogeography of many taxa has been influenced by pleistocene glaciation cycles, and this has been demonstrated in temperate regions as well as in the tropics .
for madagascar with its highly complex ecogeography of central-eastern highlands, a dry western and a humid eastern zone and various seasonal forest types, various biogeographic hypotheses have been proposed that take into account pleistocene climatic changes, topographic barriers such as mountains and large rivers  <cit> .
it has been suggested that various recent endemic radiations in lemurs and other vertebrates are causally linked to these climatic fluctuations that led to cyclic vegetation changes within and across isolated centres of endemism  <cit> .
mouse lemurs  are a highly diverse genus with currently  <dig> described species  <cit>  that are mostly confined to rather local or regional geographic ranges that are separated from each other by large rivers.
allopatric speciation has therefore been proposed as predominant mechanism of diversification within this genus  <cit> .
however, one species, the grey mouse lemur , occurs from southern madagascar up to northwestern madagascar and can be found in partial sympatry with at least five other local or regional mouse lemur species, m. griseorufus, m. berthae
, m. myoxinus, m. ravelobensis and m. bongolavensis .
two recent molecular studies suggest that this large geographic range is the result of a rather recent expansion from southern to northwestern madagascar at some time point in the late pleistocene  <cit> .
schneider et al.
<cit>  inferred this expansion based on modelling and simulating the diversity and gene genealogy of mitochondrial sequence data obtained from several populations in northwestern madagascar which also included samples from the same study site as chosen for this study.
the haplotype diversity sampled in that study could be best explained by a succession of two spatial expansions that could be dated to the late pleistocene .
it was hypothesised that an ancestral population of m. murinus colonized this region  before the last glacial maximum  and then contracted into riverine refugia during the dry period that coincided with the lgm and presumably with the preceding glacial periods in madagascar  <cit> .
from there it may subsequently have expanded again in parallel with forest expansion  <cit> .
in the second study, blair et al.
<cit>  inferred the rather recent expansion of m. murinus by modelling the split between m. murinus and its sister species, m. griseorufus, in southwestern madagascar.
coalescent methods were employed on 55– <dig> sequences for four molecular markers , respectively, stemming from localities in southern to western madagascar.
results were concordant with allopatric speciation from a narrowly distributed common ancestor that lived in southwest madagascar.
only m. murinus underwent a subsequent range expansion to the north and experienced severe population dynamics during the pleistocene  <cit> .
given these scenarios, the demographic history of m. murinus populations in northwestern madagascar, the region which it last colonized, should be substantially different from those of their sympatric congeners that most likely evolved in those areas over longer timescales  <cit> .
here we assay sequence variation in  <dig> v1r and two v2r loci in populations of the two sympatric mouse lemur species  from northwest madagascar.
we perform various tests of neutrality and investigate the correlation between past and present selective pressures.
in addition, we perform demographic analyses to assess the suitability of vr loci for making demographic inferences.
methods
data collection
we extracted dna from  <dig> grey mouse lemurs, m. murinus , and  <dig> golden-brown mouse lemurs, m. ravelobensis , using a phenol-chloroform protocol and a repli-g wga kit .
small ear biopsies were collected between may and october  <dig> in the ankarafantsika national park in northwestern madagascar by s. thorén .
all animals were live-trapped in the  <dig>  ha study site jba , and we selected individuals from different trapping locations to minimise the risk of sampling several individuals from the same sleeping group who are most likely related  <cit> .
all sampled animals of the same species are defined to belong to the same population.
we designed locus-specific primer pairs for  <dig> v1r and both v2r loci using the online software primer3plus  <cit> .
all loci  are expressed in the vno of the grey mouse lemur  <cit> .
previous work on these loci in mouse lemurs provided no evidence for paralogues or recent duplicates  <cit> .
we used twelve loci from seven of the nine monophyletic v1r clusters and three unclustered loci  <cit> .
clusters i and viii were not analysed because locus-specific primers could not be successfully designed for these.
the pcr amplicons of intronless v1r loci covered the whole locus as all primers bind outside of the coding sequence.
vn2r <dig> consists of six exons and because of long intron sequences between the exons we designed exon-specific external primers.
in contrast to our first description of v2rs  <cit> , vn2r <dig> consists of only five exons.
previously, we had used cdna of extracted rna without introns and sequenced a fragment spanning exon  <dig> to  <dig> according to closely related v2rs of family d in mice.
a short fragment  between these exons was assigned to exon  <dig>  however, the present study which uses external primers reveals that the 22 bp of exon  <dig> rather belongs to exon  <dig> which is longer in mouse lemurs than in mice.
the analogous exon  <dig> of mice is missing in mouse lemurs, which is in concordance with genomic data of the two strepsirrhines daubentonia madagascariensis and otolemur garnettii where no “exon 4” was detected previously , but a similar 22 bp insertion at the end of exon  <dig> was in fact present.
in the present study we could not design external primers for exon  <dig> of vn2r <dig> because of missing genomic data.
due to internal primers we are missing 108 bp  of the sequence of exon  <dig> and therefore results on vn2r <dig> are based on  <dig> % of the whole coding sequence.table  <dig>
v1r and v2r loci analysed with corresponding gene cluster according to hohenbrink et al.
<cit>  and total length
vn1r mmur001
vn1r mmur011
vn1r mmur031
vn1r mmur033*
vn1r mmur040
vn1r mmur041
vn1r mmur043
vn1r mmur048*
vn1r mmur049
vn1r mmur060*
vn1r mmur065
vn1r mmur066*
vn1r mmur067
vn1r mmur074*
vn1r mmur075
vn2r1
vn2r2*
a = locus 3 bp longer in m. ravelobensis, *: loci and clusters under significant positive selection  <cit> , uncl unclustered
we used mytaq dna polymerase for amplification  with the following pcr conditions: 94 °c for 2 min,  <dig> times , 72 °c for 5 min.
pcr products were sequenced on both strands using bigdye terminator  <dig>   under standard conditions and run on an applied biosystems  <dig> capillary sequencing machine.
consensus sequences were built with seqman  <dig>  .
sequences were aligned and analysed using mega  <dig>  <cit> .
before data analyses we concatenated the exons of the v2rs.
data analyses
we used dnasp  <dig>   <cit>  to unphase the two alleles of each diploid sequence and to identify the different haplotypes.
nucleotide diversity and haplotype diversity  were calculated with dnasp  <dig>  to estimate the genetic variation within the population.
dnasp  <dig>  was used to estimate the number of polymorphic synonymous  and nonsynonymous substitutions  per site  <cit>  and p
n/p
s ratios were calculated.
mcdonald-kreitman tests  were conducted online  <cit>  to calculate the ratio of fixed differences to polymorphic differences for synonymous and nonsynonymous snps.
here, all haplotypes of m. murinus and the most closely related haplotype sequence of m. ravelobensis were entered to test for positive selection in m. murinus, and vice versa to test in m. ravelobensis.
arlequin  <dig>   <cit>  was used to calculate tajima’s d and fu’s fs with  <dig> simulated samples.
population contractions or balancing selection can result in significant positive values of tajima’s d, whereas negative values indicate population growth  <cit> .
fu’s fs show negative values if the data contains an excess of rare haplotypes also indicating population growth and/or positive selection  <cit> .
in combination, positive d and positive fs values indicate an excess of intermediate-frequency alleles after population subdivision or balancing selection, whereas negative d and negative fs values reflect relative excess of rare variants and reveal population growth  <cit> .
we used false discovery rate  to correct for multiple testing across both species for each analysis, setting fdr to  <dig>   <cit> .
we analysed the distribution of nonsynonymous snps along the v1r protein .
observed vs. expected χ2-tests were used to compare the observed distribution of nonsynonymous snps in the v1r protein with the expected distribution using statistica  <dig>  .
a previous study on v1rs in strepsirrhines reported that the ligand binding site of the v1r protein is potentially formed by about half of the  <dig> and 5th transmembrane region and the in-between 2nd extracellular loop   <cit> .
this estimation was based on vr sequence data of cluster i only, but assuming structural similarities between clusters, we also tested if nonsynonymous substitutions were concentrated on the binding site proposed by yoder et al.
<cit> .
all statistical comparisons of dependent data between the two species were conducted with the wilcoxon matched pairs test in statistica.
here, the sample size was large enough to ignore p-value corrections that would have been necessary for smaller sample sizes  <cit> .
no such analyses were performed for v2rs since there is relatively little structural information available for them.
arlequin  <dig>  was also used to test two expansion models  on the mismatch distributions of the haplotypes within this single population of each species that span the same spatial scale.
a mismatch distribution is a distribution of the number of nucleotide mismatches between all pairs of nucleotide sequences of one locus within a given sample.
for this study each individual  always entered two sequences into the data pool.
mismatch distributions can be directly compared between loci or species in this study, because of the same sampling regime , same spatial spread of the samples, and the same sample size of  <dig> nucleotide sequences per locus and species.
the shape of a mismatch distribution has been shown to be influenced by demographic events like past expansions or population bottlenecks  <cit> .
mismatch distributions are bell-shape  in populations having increased in the past as a consequence of one demographic  <cit>  or spatial expansion  <cit>  or l-shaped in case of a very recent size reduction  <cit> .
a previous modelling approach  <cit>  revealed that two successive expansions can generate a tri-modal mismatch distribution and the position of these modes corresponds to the time of the expansion.
in contrast, populations at demographic equilibrium show a more ragged distribution  <cit>  .
demographic expansions usually result from past genetic bottlenecks, whereas spatial expansions usually follow a colonisation event by relatively few founder individuals.
we tested the expansion models available in arlequin  <dig>  to evaluate the evidence for a preceding colonisation event  <cit> .
the models also calculate τ- values that reflect the time of the expansion , although exact time points are difficult to infer as reliable mutation rates are often not known.
however, higher values of τ indicate that the expansion happened further in the past.
RESULTS
comparison of genetic diversity
for v1r, almost all comparisons showed substantially higher genetic diversity in m. ravelobensis than m. murinus .
across all loci, m. ravelobensis possessed significantly more polymorphic sites than m. murinus , more haplotypes , higher nucleotide diversity , higher haplotype diversity  and higher numbers of protein alleles .
all alleles at all loci had open reading frames of the expected length except locus mmur <dig> in m. ravelobensis, where an allele with an internal stop codon, indicating pseudogenization, occurred at a frequency of 6/40 =  <dig> .table  <dig> measures of genetic diversity for each locus and species
001
011
031
033
040
041
043
048
049
060
065
066
067
074
075
vn2r1
vn2r2
diff.
aa sequences: number of different amino acid sequences, mmur: m. murinus, mrav: m. ravelobensis, Ø v1r/v2r: mean
for both species, there were more haplotypes at v2r loci than v1r loci .
however, the number of polymorphic sites in v2rs was increased more dramatically in m. murinus than in m. ravelobensis.
the number of unique amino acid sequences was considerably lower than the number of haplotypes in vn2r <dig>  but this was not the case in vn2r <dig> indicating that here most haplotypes differed by at least one nonsynonymous substitution.
overall, the genetic diversity in v1rs and vn2r <dig> differs between the species, whereas it was equally high in vn2r <dig> of both species.
tests of neutrality
several loci in both species had significantly negative neutrality tests using uncorrected p-values .
however, using fdr, only three loci remained significant, all for fu’s fs in m. ravelobensis .
no significantly positive values were found and test values did not differ significantly between the species .table  <dig> results of the neutrality tests, and p
n/p
s ratios for each locus and species
p
n/p
s
001
011
031
033
040
041
043
048
049
060
065
066
067
074
075
vn2r1
vn2r2
*: p <  <dig>  with fdr of q =  <dig> ; mmur: m. murinus, mrav: m. ravelobensis, —: p
s was zero
tests for selection and distribution of mutations across vr proteins
the p
s of v1rs was significantly higher than p
n in both species .
the p
n/p
s ratios of v1rs did not differ significantly between species .
using uncorrected p-values, mcdonald-kreitman tests  were only significant for vn1r mmur <dig> in m. ravelobensis , but this result was not robust to multiple testing correction using fdr.
there was no significant correlation between the d
n/d
s of loci estimated across microcebus species  <cit>  and the p
n/p
s ratios determined intraspecifically in the current study .table  <dig> results of mk tests and mismatch distributions for each locus and species
001
011
031
033
040
041
043
048
049
060
065
066
067
074
075
vn2r1
vn2r2
*: p <  <dig>  with fdr of q =  <dig> ; mmur: m. murinus, mrav: m. ravelobensis, τ
demo: τ-values for demographic expansion model, τ
spat: τ-values for spatial expansion model, —: incalculable because of absence of variation
the distribution of nonsynonymous snps in the domains of the v1r protein did not differ significantly from the expected distribution in any species when looking at the data of all v1r loci combined .
similarly, no single locus showed a significant deviation from the expected distributions .
a proposed odorant binding site  <cit>  comprised about 20% of the v1r protein.
there was no evidence for concentration of ns mutations in this region: across all loci,  <dig> % of m. murinus ns substitutions and  <dig> % of m. ravelobensis ns substitutions occurred here, and no single locus showed an excess of ns substitutions in this region.
mismatch distributions
the mismatch distributions showed huge variation between loci and species .
in m. murinus most v1r loci  had half-bell shaped distributions close to zero pairwise differences or the peak was at zero .
furthermore, m. murinus had four loci with unimodal distributions  but no locus with multimodal or ragged distributions .
in contrast, three loci in m. ravelobensis showed ragged mismatch distributions  and six had mostly broad unimodal distributions .
the remaining six v1rs showed half-bell shaped distributions similar to m. murinus.
notably, unimodal distributions in m. murinus were still close to zero pairwise differences with a high peak, whereas they were generally broader and flat in m. ravelobensis .fig.
<dig> observed and simulated mismatch distributions of vn1r mmur <dig> and  <dig> in m. murinus  and m. ravelobensis ; the two loci were selected to show the three observed types of distributions: half-bell shaped , unimodal  and ragged ; mismatch distributions of the remaining loci are shown in additional file 2; simulated  = simulated under demographic expansion model , simulated  = simulated under spatial expansion model
m. murinus
m. ravelobensis
hb
um
rg
Σ
hb: half-bell  shape distribution, um: unimodal distribution, rg: ragged distribution, Σ: sum
the patterns of occurrence of half-bell shaped, unimodal and ragged distributions in the two species are summarized in table  <dig>  about half of the loci  had similar types of distributions in both species, but six of the remaining mismatch distributions showed higher variation in m. ravelobensis than in m. murinus.
in v2rs , the mismatch distributions in m. murinus showed a ragged distribution for vn2r <dig> and a unimodal distribution for vn2r <dig>  the distributions in m. ravelobensis were unimodal for both v2rs.
using uncorrected p-values, tests of demographic expansion were significant in one locus of m. murinus  and three v1r loci of m. ravelobensis , and tests of spatial expansion were significant in one locus of m. murinus .
however, only one of these results remained significant after fdr correction – demographic expansion for mmur <dig> in m. murinus.
the τ-values were significantly larger in m. ravelobensis compared to m. murinus using both models  indicating an older starting point of the putative expansion in m. ravelobensis compared to m. murinus.
discussion
selection in the recent evolutionary history of vomeronasal receptor genes
it was previously shown that the majority of v1r gene clusters in mouse lemurs evolved under strong positive selection and repeated gene duplication led to the evolution of a large v1r repertoire  <cit> .
positive selection still acted on v1rs during the diversification of mouse lemurs as indicated by analyses of single v1r loci across different mouse lemur species  <cit>  or of a v1r subfamily across different lemur species  <cit> .
positive selection, which was probably involved in generating this high v1r diversity in evolutionary timescales, could in principle be ongoing in present-day populations.
in contrast to this expectation, the vrs we studied seem to be mostly evolving under purifying selection.
several results support the presence of purifying selection at the population level: 1) mcdonald-kreitman tests were mostly non-significant .
2) neutrality tests were non-significant in most cases .
3) the p
n/p
s ratio was less than one in most loci of both species.
4) nonsynonymous substitutions were randomly distributed within the v1r protein indicating neutral evolution rather than positive selection.
this contrasts with our previous phylogenetic study where replacement substitutions were significantly concentrated at particular parts of the protein, consistent with odorant binding sites  <cit> .
for one locus, mmur <dig> in m. ravelobensis, there is even evidence for ongoing loss of function, with an allele encoding a pseudogene segregating at appreciable frequency.
only few loci showed patterns consistent with current positive selection.
vn1r mmur <dig> in m. ravelobensis was the only example of a significant mk test with an excess of non-synonymous snps segregating in the population, although this was not robust to fdr correction.
however, the high haplotype diversity and p
n/p
s ratio greater than one suggest that the possibility of balancing selection at this locus warrants further investigation.
three loci in m. ravelobensis showed significant evidence for departure from neutrality, with a negative fu’s fs that was robust to fdr correction.
the cause of this non-neutrality however is unclear – it may be due to directional selection or population expansion, which was not rejected for these loci from the mismatch distribution tests.
in a previous study we argued that diversification and positive selection on vr loci may be involved in reproductive isolation and speciation of mouse lemurs  <cit> .
if indeed functionally linked to reproduction, it would also be possible that the loci that are potentially under positive selection in the present study could be involved in olfactory mate choice or pheromonal communication in the context of finding a suitable mate .
selection may thus act upon the vrs in the vno and partly the main olfactory epithelium   <cit> .
although some information is available on the function of certain v1r clusters in mice , knowledge of the biological function of certain v1rs or v2rs are lacking completely for primates  <cit> .
in view of the rich and functional v1r repertoires of nocturnal strepsirrhines, future studies are urgently needed that shed light on their biological functions.
fruitful future experimental approaches may include behavioural assays involving individuals with polymorphisms at individual loci, and tests of the effect of odorants on the activity of vrs expressed in vitro.
by contrast, the use of vno tissue slices for immunohistochemistry or of anesthetised individuals for electrophysiological recordings, which has been a useful method in rodents  <cit> , is not a viable option for these endangered species.
the higher haplotype diversity for v2r genes than v1r loci is interesting.
the repertoire of v1r genes is much more diverse in mouse lemurs than their repertoire of v2r genes, with ~ <dig> v1r loci  <cit>  and only  <dig> v2r loci currently known in mouse lemurs  <cit> .
it is possible that the paucity of v2r genes in the vno is partly compensated by a high allelic diversity which also translates into a relatively high number of amino acid sequences and could therefore lead to a further increase of the olfactory sensory resolution in the vno.
the presence of only weak evidence for ongoing selection at vr loci in these two mouse lemur species is in strong contrast to the prevalence of positive selection at the same loci in a comparative study  <cit> .
this is one of the first studies to explicitly compare patterns of positive selection at these two levels, although there have been a few studies in humans  <cit> .
this obvious discordance between the results of the two studies suggests that fixation of beneficial mutations at vrs during mouse lemur evolution may have been highly episodic, i.e.
occurring over short periods of time.
however, one issue for further consideration is whether the power to detect positive selection is the same in the two approaches.
future studies would be needed to address this issue, for example using simulations.
demographic history of two sympatric mouse lemur species
under the assumption of a divergent phylogeographic history of both mouse lemur species in northwestern madagascar, we predicted to find differences in the genetic diversity of the vrs between the two species.
this prediction was confirmed by several datasets.
the species with the longer phylogeographic history in the region, m. ravelobensis, possessed significantly more polymorphic sites, a higher number of haplotypes, a higher haplotype and nucleotide diversity, as well as a higher number of different amino acid sequences in its v1r loci than its congener with the supposedly shorter phylogeographic history in the region, m. murinus.
despite some degree of variability between loci, the results from these functional loci therefore support the conclusion of previous studies about the relatively recent expansion of m. murinus into northwestern madagascar  <cit>  and suggest a distinct founder effect in this species.
a similar interspecific difference was visible in the diversity of the two v2r loci, although these were generally more diverse in m. murinus than the v1rs.
based on the scenario that m. murinus colonized northwestern madagascar only sometime in the late pleistocene  <cit> , we predicted to find signals of a stronger and/or more recent bottleneck in m. murinus than in m. ravelobensis which most likely evolved somewhere in this region earlier on  <cit>  and may have maintained a larger ancient effective population size than the expanding founder population of m. murinus.
as already noted above, although the majority of loci showed negative values in the summary statistics  of both species, which would indicate population expansion under neutrality and/or positive selection, only very few loci deviated significantly from mutation-drift equilibrium.
on the other hand, only m. ravelobensis showed significantly negative values of fu’s fs in both v2r loci.
however, the overall similarity in the patterns of these summary statistics does not support the hypotheses of a largely different demographic history of both species.
there may be several reasons for these findings: first, it is possible that these results mostly reflect the most recent demographic history of both species that may have been rather similar in the late pleistocene forests of northwestern madagascar.
the extent of the forest surface most likely contracted in all western lowland areas of madagascar towards the last glacial maximum  and expanded again only afterwards .
during the lgm, all species, independent on whether they had a long or short phylogeographic history in these lowland forests, probably underwent population contractions and expanded again afterwards together with the forests  <cit> .
in addition, both tested mouse lemur species have most likely been equally affected by the anthropogenic habitat loss that started in large parts of madagascar after the arrival of man within the last few thousand years and continues until today .
second, it is possible that similar selection regimes in the two species acting across multiple loci could affect the allelic diversity of vrs.
however, this is unlikely, since there was little evidence for widespread directional selection across vr loci .
in addition to the summary statistics discussed above, mismatch distributions were used to compare the distribution of haplotype diversity within both species.
three types of mismatch distributions were identified in both mouse lemur species, which differed in relative frequency: 1) half-bell shaped or l-shaped distributions showing a lack of variation which indicates purifying selection or recent bottlenecks  <cit> ; 2) unimodal distributions that are seen after one demographic or spatial expansion  <cit> ; 3) ragged distributions that are typical for populations at demographic equilibrium when colonisation events are very old and diversification is not constrained.
although the two species shared the half-bell shape distribution in five loci, m. ravelobensis showed the more diverse types of distributions  more often than m. murinus .
in accordance with the results on genetic diversity presented above, these findings are in agreement with the hypothesis of a larger ancient effective population size in m. ravelobensis that was able to maintain a larger degree of genetic variability across time than its congener m. murinus.
the mismatch distributions of both species did not differ significantly from the simulated distribution after one demographic or spatial expansion in most loci .
however, the τ-values were significantly higher in m. ravelobensis indicating that the putative expansion of m. ravelobensis is older than that of m. murinus.
a more precise estimation of the time since expansion is not possible, since reasonable mutation rates for vr loci are not available and evolutionary rates were shown to vary across the entire gene family  <cit> .
therefore, these species differences cannot be easily reconciled with the present knowledge on certain historic vegetation changes.
CONCLUSIONS
the current vr diversity of m. murinus and m. ravelobensis in northwestern madagascar appears to be shaped by various processes such as divergent scenarios of population expansions, purifying selection, loss of function and a potential contribution from positive selection.
whereas strong positive selection, found in the whole vr repertoire  and within individual gene clusters occurred in the past  <cit> , ongoing selection may have shifted towards purifying selection in the majority of v1r loci to maintain the adaptive function of individual receptors, e.g.
in the context of olfactory reproductive isolation between species as well as sex or kairomone recognition.
this study only analysed a small subset of the large vr repertoire of mouse lemurs but gives important insights into the recent evolution of vrs and suggests a previously unknown shift in selection pressures acting on these functional genes that are probably of highest relevance for nocturnal solitary foragers that rely heavily on olfactory communication.
functional vr loci may not be best-suited for demographic modelling considering the difficulty of differentiating between signals of purifying selection and recent population bottlenecks.
however, the simultaneous analysis of synonymous and nonsynonymous substitutions can help to disentangle these different processes.
future studies on the functional diversity and molecular evolution of olfactory receptor genes will certainly add substantially to our understanding of adaptive radiations, local adaptation and reproductive strategies in various mammalian clades whose life styles, social systems and reproductive strategies rely on pheromonal communication.
additional files
additional file 1: ids of the sampled individuals; the table contains the ids of all males and females that were used for sequencing  <dig> vr loci in the two mouse lemur species microcebus murinus  and m. ravelobensis  in the study site jba.
additional file 2: observed and simulated mismatch distributions of m. murinus and m. ravelobensis under the demographic and the spatial expansion model; graphic representation of the  <dig> mismatch distributions per species that were not included as examples in the main manuscript.
the two species-specific distributions for each locus are displayed side by side .
simulated  = simulated after demographic expansion model , simulated  = simulated after spatial expansion model .

 adaptation leaves its mark in the genome. these genomic signatures of adaptation can be investigated using phylogenetic and population genetic approaches, which address the different evolutionary timescales of interspecific and intraspecific variation. a rapidly growing body of literature is documenting the loci involved in adaptation at one of these levels. however, there are few studies which directly evaluate the relationship between the two, or separate studies examining the same loci at different timescales. hence it is an open question whether positive selection that occurs at loci among species is generally reflected in positive selection at the same loci within species. different selection regimes make different predictions for the relationship between past and ongoing selection. in cases where there is an ongoing co-evolutionary dynamic, such as in host-parasite or red queen systems, selection on key loci may be more or less continuous and a signal of selection is expected both in the past and present. under balancing selection, for example, multiple alleles may be under on-going frequency-dependent selection across speciation events, and similar signals of selection can therefore be expected to act from the past to the present. another possibility is that directional selection during a certain evolutionary period leads to the fixation of novel and advantageous variants. after this first and potentially episodic period of adaptation, purifying selection may then follow to stabilize and maintain the acquired adaptations. under these conditions, there will be discordance between inferences of past and ongoing selection. many other patterns are also possible: for example, changes in the effective population size will affect fixation of mildly deleterious alleles by drift and selection on compensatory alleles, so changes in ne alone may lead to different signatures of selection in different timeframes. in mammals, g-protein coupled chemoreceptors form one of the gene families that are most consistently found to be under positive selection in phylogenetic analyses. one of the most diverse classes of chemoreceptors in terrestrial mammals is vomeronasal receptors that function in intraspecific communication and predator detection. there is evidence for positive selection on both types of vr in mammals. this is consistent with the possibility of arms races between the scent of predators and the vrs of their prey, or the action of sexual selection on vr evolution for intraspecific communication. mouse lemurs are nocturnal strepsirrhine primates with an elaborated olfactory behavioural repertoire. olfactory stimuli play an important role in their intraspecific communication as well as in predator detection. they have the largest repertoire of vr genes of any primate. we previously showed a high level of positive selection acting on multiple families of v1r genes, as well as v2r genes, during mouse lemur evolution. we hypothesised that this positive selection may relate to predator detection and/or intraspecific communication. if this is the case then one would also expect current selection on mouse lemur vr loci. to the best of our knowledge, population genetics of vrs has not previously been investigated in the wild. functional loci have been relatively neglected in studies of demography, which favour neutral markers. the phylogeography of many taxa has been influenced by pleistocene glaciation cycles, and this has been demonstrated in temperate regions as well as in the tropics. for madagascar with its highly complex ecogeography of central-eastern highlands, a dry western and a humid eastern zone and various seasonal forest types, various biogeographic hypotheses have been proposed that take into account pleistocene climatic changes, topographic barriers such as mountains and large rivers. it has been suggested that various recent endemic radiations in lemurs and other vertebrates are causally linked to these climatic fluctuations that led to cyclic vegetation changes within and across isolated centres of endemism. mouse lemurs are a highly diverse genus with currently described species that are mostly confined to rather local or regional geographic ranges that are separated from each other by large rivers. allopatric speciation has therefore been proposed as predominant mechanism of diversification within this genus. however, one species, the grey mouse lemur, occurs from southern madagascar up to northwestern madagascar and can be found in partial sympatry with at least five other local or regional mouse lemur species, m. griseorufus, m. berthae, m. myoxinus, m. ravelobensis and m. bongolavensis. two recent molecular studies suggest that this large geographic range is the result of a rather recent expansion from southern to northwestern madagascar at some time point in the late pleistocene. schneider et al. inferred this expansion based on modelling and simulating the diversity and gene genealogy of mitochondrial sequence data obtained from several populations in northwestern madagascar which also included samples from the same study site as chosen for this study. the haplotype diversity sampled in that study could be best explained by a succession of two spatial expansions that could be dated to the late pleistocene. it was hypothesised that an ancestral population of m. murinus colonized this region before the last glacial maximum and then contracted into riverine refugia during the dry period that coincided with the lgm and presumably with the preceding glacial periods in madagascar. from there it may subsequently have expanded again in parallel with forest expansion. in the second study, blair et al. inferred the rather recent expansion of m. murinus by modelling the split between m. murinus and its sister species, m. griseorufus, in southwestern madagascar. coalescent methods were employed on 55– sequences for four molecular markers, respectively, stemming from localities in southern to western madagascar. results were concordant with allopatric speciation from a narrowly distributed common ancestor that lived in southwest madagascar. only m. murinus underwent a subsequent range expansion to the north and experienced severe population dynamics during the pleistocene. given these scenarios, the demographic history of m. murinus populations in northwestern madagascar, the region which it last colonized, should be substantially different from those of their sympatric congeners that most likely evolved in those areas over longer timescales. here we assay sequence variation in v1r and two v2r loci in populations of the two sympatric mouse lemur species from northwest madagascar. we perform various tests of neutrality and investigate the correlation between past and present selective pressures. in addition, we perform demographic analyses to assess the suitability of vr loci for making demographic inferences. methods data collection we extracted dna from grey mouse lemurs, m. murinus, and golden-brown mouse lemurs, m. ravelobensis, using a phenol-chloroform protocol and a repli-g wga kit. small ear biopsies were collected between may and october in the ankarafantsika national park in northwestern madagascar by s. thorén. all animals were live-trapped in the  ha study site jba, and we selected individuals from different trapping locations to minimise the risk of sampling several individuals from the same sleeping group who are most likely related. all sampled animals of the same species are defined to belong to the same population. we designed locus-specific primer pairs for v1r and both v2r loci using the online software primer3plus. all loci are expressed in the vno of the grey mouse lemur. previous work on these loci in mouse lemurs provided no evidence for paralogues or recent duplicates. we used twelve loci from seven of the nine monophyletic v1r clusters and three unclustered loci. clusters i and viii were not analysed because locus-specific primers could not be successfully designed for these. the pcr amplicons of intronless v1r loci covered the whole locus as all primers bind outside of the coding sequence. vn2r consists of six exons and because of long intron sequences between the exons we designed exon-specific external primers. in contrast to our first description of v2rs, vn2r consists of only five exons. previously, we had used cdna of extracted rna without introns and sequenced a fragment spanning exon to according to closely related v2rs of family d in mice. a short fragment between these exons was assigned to exon however, the present study which uses external primers reveals that the 22 bp of exon rather belongs to exon which is longer in mouse lemurs than in mice. the analogous exon of mice is missing in mouse lemurs, which is in concordance with genomic data of the two strepsirrhines daubentonia madagascariensis and otolemur garnettii where no “exon 4” was detected previously, but a similar 22 bp insertion at the end of exon was in fact present. in the present study we could not design external primers for exon of vn2r because of missing genomic data. due to internal primers we are missing 108 bp of the sequence of exon and therefore results on vn2r are based on % of the whole coding sequence.table v1r and v2r loci analysed with corresponding gene cluster according to hohenbrink et al. and total length vn1r mmur001 vn1r mmur011 vn1r mmur031 vn1r mmur033* vn1r mmur040 vn1r mmur041 vn1r mmur043 vn1r mmur048* vn1r mmur049 vn1r mmur060* vn1r mmur065 vn1r mmur066* vn1r mmur067 vn1r mmur074* vn1r mmur075 vn2r1 vn2r2* a = locus 3 bp longer in m. ravelobensis, *: loci and clusters under significant positive selection, uncl unclustered we used mytaq dna polymerase for amplification with the following pcr conditions: 94 °c for 2 min, times, 72 °c for 5 min. pcr products were sequenced on both strands using bigdye terminator under standard conditions and run on an applied biosystems capillary sequencing machine. consensus sequences were built with seqman. sequences were aligned and analysed using mega. before data analyses we concatenated the exons of the v2rs. data analyses we used dnasp to unphase the two alleles of each diploid sequence and to identify the different haplotypes. nucleotide diversity and haplotype diversity were calculated with dnasp to estimate the genetic variation within the population. dnasp was used to estimate the number of polymorphic synonymous and nonsynonymous substitutions per site and p n/p s ratios were calculated. mcdonald-kreitman tests were conducted online to calculate the ratio of fixed differences to polymorphic differences for synonymous and nonsynonymous snps. here, all haplotypes of m. murinus and the most closely related haplotype sequence of m. ravelobensis were entered to test for positive selection in m. murinus, and vice versa to test in m. ravelobensis. arlequin was used to calculate tajima’s d and fu’s fs with simulated samples. population contractions or balancing selection can result in significant positive values of tajima’s d, whereas negative values indicate population growth. fu’s fs show negative values if the data contains an excess of rare haplotypes also indicating population growth and/or positive selection. in combination, positive d and positive fs values indicate an excess of intermediate-frequency alleles after population subdivision or balancing selection, whereas negative d and negative fs values reflect relative excess of rare variants and reveal population growth. we used false discovery rate to correct for multiple testing across both species for each analysis, setting fdr to. we analysed the distribution of nonsynonymous snps along the v1r protein. observed vs. expected χ2-tests were used to compare the observed distribution of nonsynonymous snps in the v1r protein with the expected distribution using statistica. a previous study on v1rs in strepsirrhines reported that the ligand binding site of the v1r protein is potentially formed by about half of the and 5th transmembrane region and the in-between 2nd extracellular loop. this estimation was based on vr sequence data of cluster i only, but assuming structural similarities between clusters, we also tested if nonsynonymous substitutions were concentrated on the binding site proposed by yoder et al.. all statistical comparisons of dependent data between the two species were conducted with the wilcoxon matched pairs test in statistica. here, the sample size was large enough to ignore p-value corrections that would have been necessary for smaller sample sizes. no such analyses were performed for v2rs since there is relatively little structural information available for them. arlequin was also used to test two expansion models on the mismatch distributions of the haplotypes within this single population of each species that span the same spatial scale. a mismatch distribution is a distribution of the number of nucleotide mismatches between all pairs of nucleotide sequences of one locus within a given sample. for this study each individual always entered two sequences into the data pool. mismatch distributions can be directly compared between loci or species in this study, because of the same sampling regime, same spatial spread of the samples, and the same sample size of nucleotide sequences per locus and species. the shape of a mismatch distribution has been shown to be influenced by demographic events like past expansions or population bottlenecks. mismatch distributions are bell-shape in populations having increased in the past as a consequence of one demographic or spatial expansion or l-shaped in case of a very recent size reduction. a previous modelling approach revealed that two successive expansions can generate a tri-modal mismatch distribution and the position of these modes corresponds to the time of the expansion. in contrast, populations at demographic equilibrium show a more ragged distribution. demographic expansions usually result from past genetic bottlenecks, whereas spatial expansions usually follow a colonisation event by relatively few founder individuals. we tested the expansion models available in arlequin to evaluate the evidence for a preceding colonisation event. the models also calculate τ- values that reflect the time of the expansion, although exact time points are difficult to infer as reliable mutation rates are often not known. however, higher values of τ indicate that the expansion happened further in the past.  comparison of genetic diversity for v1r, almost all comparisons showed substantially higher genetic diversity in m. ravelobensis than m. murinus. across all loci, m. ravelobensis possessed significantly more polymorphic sites than m. murinus, more haplotypes, higher nucleotide diversity, higher haplotype diversity and higher numbers of protein alleles. all alleles at all loci had open reading frames of the expected length except locus mmur in m. ravelobensis, where an allele with an internal stop codon, indicating pseudogenization, occurred at a frequency of 6/40 = .table measures of genetic diversity for each locus and species 001 011 031 033 040 041 043 048 049 060 065 066 067 074 075 vn2r1 vn2r2 diff. aa sequences: number of different amino acid sequences, mmur: m. murinus, mrav: m. ravelobensis, Ø v1r/v2r: mean for both species, there were more haplotypes at v2r loci than v1r loci. however, the number of polymorphic sites in v2rs was increased more dramatically in m. murinus than in m. ravelobensis. the number of unique amino acid sequences was considerably lower than the number of haplotypes in vn2r but this was not the case in vn2r indicating that here most haplotypes differed by at least one nonsynonymous substitution. overall, the genetic diversity in v1rs and vn2r differs between the species, whereas it was equally high in vn2r of both species. tests of neutrality several loci in both species had significantly negative neutrality tests using uncorrected p-values. however, using fdr, only three loci remained significant, all for fu’s fs in m. ravelobensis. no significantly positive values were found and test values did not differ significantly between the species.table results of the neutrality tests, and p n/p s ratios for each locus and species p n/p s 001 011 031 033 040 041 043 048 049 060 065 066 067 074 075 vn2r1 vn2r2 *: p <  with fdr of q = ; mmur: m. murinus, mrav: m. ravelobensis, —: p s was zero tests for selection and distribution of mutations across vr proteins the p s of v1rs was significantly higher than p n in both species. the p n/p s ratios of v1rs did not differ significantly between species. using uncorrected p-values, mcdonald-kreitman tests were only significant for vn1r mmur in m. ravelobensis, but this result was not robust to multiple testing correction using fdr. there was no significant correlation between the d n/d s of loci estimated across microcebus species and the p n/p s ratios determined intraspecifically in the current study.table results of mk tests and mismatch distributions for each locus and species 001 011 031 033 040 041 043 048 049 060 065 066 067 074 075 vn2r1 vn2r2 *: p <  with fdr of q = ; mmur: m. murinus, mrav: m. ravelobensis, τ demo: τ-values for demographic expansion model, τ spat: τ-values for spatial expansion model, —: incalculable because of absence of variation the distribution of nonsynonymous snps in the domains of the v1r protein did not differ significantly from the expected distribution in any species when looking at the data of all v1r loci combined. similarly, no single locus showed a significant deviation from the expected distributions. a proposed odorant binding site comprised about 20% of the v1r protein. there was no evidence for concentration of ns mutations in this region: across all loci, % of m. murinus ns substitutions and % of m. ravelobensis ns substitutions occurred here, and no single locus showed an excess of ns substitutions in this region. mismatch distributions the mismatch distributions showed huge variation between loci and species. in m. murinus most v1r loci had half-bell shaped distributions close to zero pairwise differences or the peak was at zero. furthermore, m. murinus had four loci with unimodal distributions but no locus with multimodal or ragged distributions. in contrast, three loci in m. ravelobensis showed ragged mismatch distributions and six had mostly broad unimodal distributions. the remaining six v1rs showed half-bell shaped distributions similar to m. murinus. notably, unimodal distributions in m. murinus were still close to zero pairwise differences with a high peak, whereas they were generally broader and flat in m. ravelobensis.fig. observed and simulated mismatch distributions of vn1r mmur and in m. murinus and m. ravelobensis; the two loci were selected to show the three observed types of distributions: half-bell shaped, unimodal and ragged; mismatch distributions of the remaining loci are shown in additional file 2; simulated  = simulated under demographic expansion model, simulated  = simulated under spatial expansion model m. murinus m. ravelobensis hb um rg Σ hb: half-bell shape distribution, um: unimodal distribution, rg: ragged distribution, Σ: sum the patterns of occurrence of half-bell shaped, unimodal and ragged distributions in the two species are summarized in table  about half of the loci had similar types of distributions in both species, but six of the remaining mismatch distributions showed higher variation in m. ravelobensis than in m. murinus. in v2rs, the mismatch distributions in m. murinus showed a ragged distribution for vn2r and a unimodal distribution for vn2r the distributions in m. ravelobensis were unimodal for both v2rs. using uncorrected p-values, tests of demographic expansion were significant in one locus of m. murinus and three v1r loci of m. ravelobensis, and tests of spatial expansion were significant in one locus of m. murinus. however, only one of these results remained significant after fdr correction – demographic expansion for mmur in m. murinus. the τ-values were significantly larger in m. ravelobensis compared to m. murinus using both models indicating an older starting point of the putative expansion in m. ravelobensis compared to m. murinus. discussion selection in the recent evolutionary history of vomeronasal receptor genes it was previously shown that the majority of v1r gene clusters in mouse lemurs evolved under strong positive selection and repeated gene duplication led to the evolution of a large v1r repertoire. positive selection still acted on v1rs during the diversification of mouse lemurs as indicated by analyses of single v1r loci across different mouse lemur species or of a v1r subfamily across different lemur species. positive selection, which was probably involved in generating this high v1r diversity in evolutionary timescales, could in principle be ongoing in present-day populations. in contrast to this expectation, the vrs we studied seem to be mostly evolving under purifying selection. several results support the presence of purifying selection at the population level: 1) mcdonald-kreitman tests were mostly non-significant. 2) neutrality tests were non-significant in most cases. 3) the p n/p s ratio was less than one in most loci of both species. 4) nonsynonymous substitutions were randomly distributed within the v1r protein indicating neutral evolution rather than positive selection. this contrasts with our previous phylogenetic study where replacement substitutions were significantly concentrated at particular parts of the protein, consistent with odorant binding sites. for one locus, mmur in m. ravelobensis, there is even evidence for ongoing loss of function, with an allele encoding a pseudogene segregating at appreciable frequency. only few loci showed patterns consistent with current positive selection. vn1r mmur in m. ravelobensis was the only example of a significant mk test with an excess of non-synonymous snps segregating in the population, although this was not robust to fdr correction. however, the high haplotype diversity and p n/p s ratio greater than one suggest that the possibility of balancing selection at this locus warrants further investigation. three loci in m. ravelobensis showed significant evidence for departure from neutrality, with a negative fu’s fs that was robust to fdr correction. the cause of this non-neutrality however is unclear – it may be due to directional selection or population expansion, which was not rejected for these loci from the mismatch distribution tests. in a previous study we argued that diversification and positive selection on vr loci may be involved in reproductive isolation and speciation of mouse lemurs. if indeed functionally linked to reproduction, it would also be possible that the loci that are potentially under positive selection in the present study could be involved in olfactory mate choice or pheromonal communication in the context of finding a suitable mate. selection may thus act upon the vrs in the vno and partly the main olfactory epithelium. although some information is available on the function of certain v1r clusters in mice, knowledge of the biological function of certain v1rs or v2rs are lacking completely for primates. in view of the rich and functional v1r repertoires of nocturnal strepsirrhines, future studies are urgently needed that shed light on their biological functions. fruitful future experimental approaches may include behavioural assays involving individuals with polymorphisms at individual loci, and tests of the effect of odorants on the activity of vrs expressed in vitro. by contrast, the use of vno tissue slices for immunohistochemistry or of anesthetised individuals for electrophysiological recordings, which has been a useful method in rodents, is not a viable option for these endangered species. the higher haplotype diversity for v2r genes than v1r loci is interesting. the repertoire of v1r genes is much more diverse in mouse lemurs than their repertoire of v2r genes, with ~ v1r loci and only v2r loci currently known in mouse lemurs. it is possible that the paucity of v2r genes in the vno is partly compensated by a high allelic diversity which also translates into a relatively high number of amino acid sequences and could therefore lead to a further increase of the olfactory sensory resolution in the vno. the presence of only weak evidence for ongoing selection at vr loci in these two mouse lemur species is in strong contrast to the prevalence of positive selection at the same loci in a comparative study. this is one of the first studies to explicitly compare patterns of positive selection at these two levels, although there have been a few studies in humans. this obvious discordance between the results of the two studies suggests that fixation of beneficial mutations at vrs during mouse lemur evolution may have been highly episodic, i.e. occurring over short periods of time. however, one issue for further consideration is whether the power to detect positive selection is the same in the two approaches. future studies would be needed to address this issue, for example using simulations. demographic history of two sympatric mouse lemur species under the assumption of a divergent phylogeographic history of both mouse lemur species in northwestern madagascar, we predicted to find differences in the genetic diversity of the vrs between the two species. this prediction was confirmed by several datasets. the species with the longer phylogeographic history in the region, m. ravelobensis, possessed significantly more polymorphic sites, a higher number of haplotypes, a higher haplotype and nucleotide diversity, as well as a higher number of different amino acid sequences in its v1r loci than its congener with the supposedly shorter phylogeographic history in the region, m. murinus. despite some degree of variability between loci, the results from these functional loci therefore support the conclusion of previous studies about the relatively recent expansion of m. murinus into northwestern madagascar and suggest a distinct founder effect in this species. a similar interspecific difference was visible in the diversity of the two v2r loci, although these were generally more diverse in m. murinus than the v1rs. based on the scenario that m. murinus colonized northwestern madagascar only sometime in the late pleistocene, we predicted to find signals of a stronger and/or more recent bottleneck in m. murinus than in m. ravelobensis which most likely evolved somewhere in this region earlier on and may have maintained a larger ancient effective population size than the expanding founder population of m. murinus. as already noted above, although the majority of loci showed negative values in the summary statistics of both species, which would indicate population expansion under neutrality and/or positive selection, only very few loci deviated significantly from mutation-drift equilibrium. on the other hand, only m. ravelobensis showed significantly negative values of fu’s fs in both v2r loci. however, the overall similarity in the patterns of these summary statistics does not support the hypotheses of a largely different demographic history of both species. there may be several reasons for these findings: first, it is possible that these results mostly reflect the most recent demographic history of both species that may have been rather similar in the late pleistocene forests of northwestern madagascar. the extent of the forest surface most likely contracted in all western lowland areas of madagascar towards the last glacial maximum and expanded again only afterwards. during the lgm, all species, independent on whether they had a long or short phylogeographic history in these lowland forests, probably underwent population contractions and expanded again afterwards together with the forests. in addition, both tested mouse lemur species have most likely been equally affected by the anthropogenic habitat loss that started in large parts of madagascar after the arrival of man within the last few thousand years and continues until today. second, it is possible that similar selection regimes in the two species acting across multiple loci could affect the allelic diversity of vrs. however, this is unlikely, since there was little evidence for widespread directional selection across vr loci. in addition to the summary statistics discussed above, mismatch distributions were used to compare the distribution of haplotype diversity within both species. three types of mismatch distributions were identified in both mouse lemur species, which differed in relative frequency: 1) half-bell shaped or l-shaped distributions showing a lack of variation which indicates purifying selection or recent bottlenecks; 2) unimodal distributions that are seen after one demographic or spatial expansion; 3) ragged distributions that are typical for populations at demographic equilibrium when colonisation events are very old and diversification is not constrained. although the two species shared the half-bell shape distribution in five loci, m. ravelobensis showed the more diverse types of distributions more often than m. murinus. in accordance with the results on genetic diversity presented above, these findings are in agreement with the hypothesis of a larger ancient effective population size in m. ravelobensis that was able to maintain a larger degree of genetic variability across time than its congener m. murinus. the mismatch distributions of both species did not differ significantly from the simulated distribution after one demographic or spatial expansion in most loci. however, the τ-values were significantly higher in m. ravelobensis indicating that the putative expansion of m. ravelobensis is older than that of m. murinus. a more precise estimation of the time since expansion is not possible, since reasonable mutation rates for vr loci are not available and evolutionary rates were shown to vary across the entire gene family. therefore, these species differences cannot be easily reconciled with the present knowledge on certain historic vegetation changes.  the current vr diversity of m. murinus and m. ravelobensis in northwestern madagascar appears to be shaped by various processes such as divergent scenarios of population expansions, purifying selection, loss of function and a potential contribution from positive selection. whereas strong positive selection, found in the whole vr repertoire and within individual gene clusters occurred in the past, ongoing selection may have shifted towards purifying selection in the majority of v1r loci to maintain the adaptive function of individual receptors, e.g. in the context of olfactory reproductive isolation between species as well as sex or kairomone recognition. this study only analysed a small subset of the large vr repertoire of mouse lemurs but gives important insights into the recent evolution of vrs and suggests a previously unknown shift in selection pressures acting on these functional genes that are probably of highest relevance for nocturnal solitary foragers that rely heavily on olfactory communication. functional vr loci may not be best-suited for demographic modelling considering the difficulty of differentiating between signals of purifying selection and recent population bottlenecks. however, the simultaneous analysis of synonymous and nonsynonymous substitutions can help to disentangle these different processes. future studies on the functional diversity and molecular evolution of olfactory receptor genes will certainly add substantially to our understanding of adaptive radiations, local adaptation and reproductive strategies in various mammalian clades whose life styles, social systems and reproductive strategies rely on pheromonal communication. additional files additional file 1: ids of the sampled individuals; the table contains the ids of all males and females that were used for sequencing vr loci in the two mouse lemur species microcebus murinus and m. ravelobensis in the study site jba. additional file 2: observed and simulated mismatch distributions of m. murinus and m. ravelobensis under the demographic and the spatial expansion model; graphic representation of the mismatch distributions per species that were not included as examples in the main manuscript. the two species-specific distributions for each locus are displayed side by side. simulated  = simulated after demographic expansion model, simulated  = simulated after spatial expansion model. 
here, we re-sequence vrs and perform a suite of selection and demographic analyses in sympatric populations of two species of mouse lemurs in northwestern madagascar. m. ravelobensis had consistently higher genetic diversity at vrs than m. murinus. in general, we find little evidence for positive selection, with most loci evolving under purifying selection and one locus even showing evidence of functional loss in m. ravelobensis. however, a few loci in m. ravelobensis show potential evidence of positive selection. using mismatch distributions and expansion models, we infer a more recent colonisation of the habitat by m. murinus than by m. ravelobensis, which most likely speciated in this region earlier on. we address this issue using vomeronasal receptors, a diverse gene family in mammals involved in intraspecific communication and predator detection. a major effort is underway to use population genetic approaches to identify loci involved in adaptation. in conclusion, this study reveals a substantial heterogeneity over time in selection on vr loci, suggesting that vr evolution is episodic. vnoevolutionmadagascarmicrocebus murinusmicrocebus ravelobensisv1rv2rgenetic diversitydemographyselectionhttp://dx.doi.org/ /501100001663volkswagen foundationi/84798hohenbrink philipp issue-copyright-statement© the author 2017 in mouse lemurs, we previously demonstrated that both subfamilies of vrs show a strong signal of directional selection in interspecific analyses. these findings suggest that the analysis of vr variation is useful in inferring demographic and phylogeographic history of mouse lemurs. one issue that has so far received limited attention is whether loci that show a phylogenetic signal of positive selection in the past also show evidence of ongoing positive selection at the population level.
 a major effort is underway to use population genetic approaches to identify loci involved in adaptation. one issue that has so far received limited attention is whether loci that show a phylogenetic signal of positive selection in the past also show evidence of ongoing positive selection at the population level. we address this issue using vomeronasal receptors, a diverse gene family in mammals involved in intraspecific communication and predator detection. in mouse lemurs, we previously demonstrated that both subfamilies of vrs show a strong signal of directional selection in interspecific analyses. we predicted that ongoing sexual selection and/or co-evolution with predators may lead to current directional or balancing selection on vrs. here, we re-sequence vrs and perform a suite of selection and demographic analyses in sympatric populations of two species of mouse lemurs in northwestern madagascar.  m. ravelobensis had consistently higher genetic diversity at vrs than m. murinus. in general, we find little evidence for positive selection, with most loci evolving under purifying selection and one locus even showing evidence of functional loss in m. ravelobensis. however, a few loci in m. ravelobensis show potential evidence of positive selection. using mismatch distributions and expansion models, we infer a more recent colonisation of the habitat by m. murinus than by m. ravelobensis, which most likely speciated in this region earlier on.  these findings suggest that the analysis of vr variation is useful in inferring demographic and phylogeographic history of mouse lemurs. in conclusion, this study reveals a substantial heterogeneity over time in selection on vr loci, suggesting that vr evolution is episodic. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords vnoevolutionmadagascarmicrocebus murinusmicrocebus ravelobensisv1rv2rgenetic diversitydemographyselectionhttp://dx.doi.org/ /501100001663volkswagen foundationi/84798hohenbrink philipp issue-copyright-statement© the author 2017 