BACKGROUND
in recent years, there has been growing interest in using genealogical tree structure to reconstruct the demographic and temporal context of diverging populations and species  <cit> .
this tree-based approach permits the examination of lineages as they have diverged in the past using principles derived from coalescent theory, since ancestral polymorphisms that are shared among lineages are sorted during population segregation and speciation.
a comparative demographic approach utilizing tree-based methods has been implemented in a variety of biological sub-fields, including systematics, phylogeography, conservation, and life history.
further, estimates of divergence times, migration rates, and effective population sizes among groups of organisms are now inferred with consideration of the genealogical structure of multilocus data sets  <cit> .
these estimates may be synthesized and compared among sister taxa and used to infer the impact of the biogeographic or climatic context or other historical processes on patterns of genetic differences among groups  <cit> .
in general, it is assumed that a complex suite of geologic, climatic, and population genetic forces have led to divergence and speciation in allopatry, where a reproductive barrier effectively divides an ancestral species into two populations of roughly equal size.
with the continual progress of next-generation dna sequencing and the increasing abundance of sequence data from across the genome, researchers are beginning to test alternate models of speciation including parapatric divergence along ecological clines  <cit> .
few studies, however, have used recently developed multilocus coalescent methods to fully understand the geography of speciation and the frequency of peripatric events in nature.
peripatric speciation can be seen as a subset of allopatric speciation, with a founder effect leading to the formation of a new species as a small population becomes physically separated from a broadly distributed ancestor.
by definition then, the effective population size  of the diverging population will be smaller than that of the common ancestor in a peripatric scenario.
although evidence of peripatric divergence is commonly inferred throughout archipelagos using traditional phylogenetic approaches , examples from strictly mainland taxa are relatively rare.
new algorithms can effectively estimate demographic parameters on a species tree to help disentangle the historical context of divergence and speciation and to test the plausibility of a peripatric scenario in comparison with the generally accepted null model of allopatric divergence  <cit> .
mitochondrial dna  has been commonly used in studies of speciation given its small effective population size relative to nuclear dna  loci.
accordingly, mitochondrial gene trees undergo lineage sorting much more rapidly, and thus have a higher probability of yielding gene trees that are reciprocally monophyletic for the hypothesized species.
however, the stochasticity inherent in single-locus gene trees makes this a controversial approach for estimating species histories, as previous authors have discussed .
therefore, the implementation of a comparative demographic approach incorporating multiple unlinked loci is now regarded as a more powerful approach.
mouse lemurs  are small, nocturnal primates that inhabit a wide range of habitats throughout madagascar.
as many as  <dig> species have been described within the genus  <cit> , thus a remarkable increase from the two species taxonomy that was recognized as recently as  <dig>  <cit> .
the geography of the island has historically been considered to play a significant role in explaining divergence patterns among lemurs and other vertebrate species  and specifically microcebus species  <cit> .
within the genus, however, m. murinus is unique in having a relatively vast geographic distribution.
whereas all other microcebus species have geographically limited ranges --- some restricted to isolated forest fragments --- m. murinus is distributed along most of the western half and along the south of the island .
over the extensive range of m. murinus, the species has been found to overlap with multiple congeners  <cit> , many of which exhibit a high degree of endemicity.
one example of this pattern of sympatry couples m. murinus with m. griseorufus, with multiple lines of evidence for distribution overlap and potential hybridization between the two species in southern madagascar  <cit> .
within their respective ranges, m. murinus and m. griseorufus also exhibit marked differences in abundance.
microcebus murinus is a common, widespread species throughout the west  that has been studied extensively since its description  <cit> .
conversely, m. griseorufus has been recorded from relatively few locations within a narrow range in the southwest of the island in the recent past  <cit> .
despite their partially overlapping distributions and large genetic distances  <cit> , m. murinus and m. griseorufus are sister taxa  <cit> .
given the current ranges of m. murinus and m. griseorufus and their status as sister species, we consider two possible scenarios for the speciation event within this small clade.
one scenario is that m. griseorufus was derived from the most recent common ancestor  as a small founding population, with divergence the result of random genetic drift .
alternatively, we consider that the mrca was a small, narrowly distributed population located in southwestern madagascar.
under this alternative hypothesis, the two modern species resulted from an event that partitioned that ancestral species into comparably sized populations prior to a substantial range expansion in m. murinus .
these two hypotheses differ most critically in their assumptions regarding the ancestral range and population size: was the mrca more like modern m. murinus  or more like modern m. griseorufus ?
to test these alternative scenarios, we analyze multilocus dna sequence data with coalescent-based methods to estimate relevant demographic parameters.
if m. griseorufus diverged as a small founding population from a large ancestral population , then we expect 1) that m. murinus has a greater effective population size and higher levels of genetic diversity than m. griseorufus, 2) that estimates of effective population size for the mrca are more similar to those for m. murinus than to those for m. griseorufus, 3) that historical changes in population size in m. murinus were minimal, and finally, 4) that m. griseorufus is more likely to segregate as a clade in individual gene trees.
under the alternative hypothesis of recent allopatric speciation , we expect that 1) that m. murinus and m. griseorufus have similar effective population sizes and levels of genetic diversity, 2) that ne estimates for the mrca are smaller than those for contemporary species or similar to m. griseorufus, 3) that significant population expansion has occurred in m. murinus, and finally, 4) that individual gene trees are more likely to show incomplete lineage sorting between the two species.
methods
data assembly
to estimate relevant demographic parameters for both species, we mined genbank for all available high-quality sequences as both m. murinus and m. griseorufus have been used in a number of recent molecular studies .
we sought to utilize genes with sufficient sequence variation to reliably estimate coalescent parameters.
we also chose data sets with samples covering a large portion of the known geographic range of both species to assemble the most geographically complete data set possible .
although recent studies have suggested that m. murinus may possibly contain at least three undocumented species,  <cit> , we treat these populations as m. murinus sensu lato for the purpose of this study.
the final data set consisted of the following nuclear loci: alpha enolase intron , alpha fibrinogen intron , von willebrand factor intron .
although the adora  <dig> gene has been sequenced numerous times within microcebus, levels of variation were too low to provide any meaningful information for this study.
we also assembled and concatenated information from two mitochondrial loci including cytochrome b  and cytochrome c oxidase subunit ii  for a total of  <dig>  bp.
the total number of sequences per locus was as follows: enol =  <dig> sequences; fib =  <dig> sequences; vwf =  <dig> sequences; mtdna =  <dig> concatenated sequences.
all of the nuclear sequences represented phased haplotypes, which were included for all subsequent analyses.
cheirogaleus major was used as an outgroup for all phylogenetic and demographic analyses .
although cheirogaleus can be considered a relatively distant outgroup, we were interested in adding a temporal component to our analyses  and sought to maintain consistency with taxon sampling across analyses.
all sequence manipulation was performed in geneious v.  <dig> .
<dig>  <cit> .
multiple sequence alignments were conducted using mafft v.  <dig>   <cit>  within geneious.
sequence diversity and effective population sizes
we first tested each marker for signs of recombination using the program rdp v.  <dig>   <cit> .
each marker was tested for recombination events using the geneconv  <cit> , maxchi  <cit> , and rdp methods  <cit> .
default settings were used for all analyses.
because no signs of recombination were detected, all subsequent analyses utilized the entire read of each locus.
general sequence diversity statistics for each locus and species including number of haplotypes, haplotype diversity, number of segregating sites, and nucleotide diversity were calculated using dnasp v.  <dig> .
<dig>  <cit> .
using the species tree topology , cheirogaleus major) prescribed by both nuclear and mitochondrial phylogenies  <cit> , we used the program bayesian phylogenetics and phylogeography  to estimate coalescent-scaled population sizes  and time of divergence .
this method accommodates the species tree divergence patterns as well as gene tree lineage sorting processes.
a gamma prior g, with mean 2/1000 =  <dig>  was used for the population size parameters .
the age of the root in the species tree  was also assigned a gamma prior g, while the other divergence time parameters were assigned the dirichlet prior .
locus-specific mutation rates were allowed to vary, and we specified a heredity multiplier value g to account for the combined mtdna and ndna data.
we also performed analyses using larger priors for both τ <dig> and θ ) to evaluate the sensitivity of our results to the choice of prior.
for this study, the parameters of interest were θm, θg, θmg, and τmg, where m = m. murinus, g = m. griseorufus, and mg = the mrca of m. murinus and m. griseorufus.
the default numbers of generations and sampling intervals were used for all analyses.
due to computational issues with the full data set, all bpp analyses were implemented on a reduced data set consisting of approximately  <dig> sequences per locus per microcebus species and one to six cheirogaleus sequences per locus.
sequences were randomly sampled to encompass a broad geographic range for each species.
previous simulation-based studies have shown that similar sample sizes can be sufficient to infer speciation processes  <cit> .
additional analyses on further reduced data sets  yielded similar results.
all analyses were run at least twice to check for consistency among runs.
demographic expansion
to further test between the alternative speciation hypotheses, and to determine if the relatively large geographic range of m. murinus was due to population expansion following the split with m. griseorufus, we tested for signs of demographic change through time for both species.
under a peripatric speciation model, an ancestral species with a large range similar to that of modern m. murinus would have diverged into m. murinus and m. griseorufus populations, with the former experiencing little to no demographic change associated with the speciation event and the latter experiencing a significant population bottleneck event.
conversely, under the allopatric speciation model, a narrowly-distributed ancestral species would have diverged into m. murinus and m. griseorufus populations, with the former subsequently experiencing a significant population growth event and the latter experiencing a relatively smaller but possibly detectable demographic change after the speciation event.
to test for signs of demographic expansion or contraction we implemented multilocus bayesian methods.
specifically, we constructed extended bayesian skyline plots  for both species using beast v.  <dig> .
<dig>  <cit> .
studies have shown that doubling the number of independent loci can reduce error and 95% credible intervals in demographic reconstruction by √ <dig>  <cit> .
thus, compared to single locus estimates, our multilocus data provided a powerful approach for estimating demographic trends in mouse lemurs.
first, we used jmodeltest v.  <dig> .
<dig>  <cit>  to calculate model likelihood scores for each locus and to estimate optimal models using bic .
because mtdna is inherited as a single linked unit, and to minimize the computational burden for beast, the concatenated mtdna data were treated as a single partition.
the three nuclear loci were each specified as a separate partition for model fitting.
we tested the likelihood of  <dig> commonly used models in beast.
the first two blocks represent per locus models selected for extended bayesian skyline plot analyses.
‘*beast’ represents a reduced data set encompassing both species and the outgroup .
we were also interested in adding a temporal component to the demographic analyses.
although there is no fossil record for lemurs, recent studies have utilized fossil information from more distantly related groups to date divergence times within the lemurs .
using multiple calibration points outside the clade and relaxed clock methods, these authors estimated the split between cheirogaleus and microcebus to be approximately  <dig> million years ago .
we first used this information to estimate the substitution rate for each of the test loci in beast using a reduced data set of both species and cheirogaleus .
best-fitting models were calculated and used for all divergence dating and rate estimation.
because alignments contained representative alleles from both multiple species and multiple individuals within species, we used *beast  <cit>  to estimate the posterior distribution of substitution rates.
we defined three species  and grouped alleles accordingly.
for each analysis, the root node of the species tree was calibrated with a normal distribution around a mean of 25 ma and standard deviation of 5 ma, which encompassed the 95% hpd estimates from previous studies  <cit> .
to increase computational efficiency, we ran four independent *beast analyses  to estimate the posterior distribution of rates.
both strict and relaxed clock  models were tested.
all analyses were run between 10– <dig> million generations; sampling was chosen at intervals to utilize  <dig>  draws from the posterior.
following analyses, the program tracer v <dig>   <cit>  was used to examine effective sample size  values  and examine the posterior distribution of relevant parameters.
we used the 95% hpd of the substitution rate for each locus as a uniform prior to add a temporal component to all ebsp analyses.
all ebsp analyses used a strict clock.
operators were modified according to author recommendations and analyses were run for  <dig> million generations  or  <dig> million generations  to obtain adequate ess values.
all ebsp and *beast analyses were implemented via the duke shared cluster resource .
phylogenetic analysis
under a scenario of peripatric speciation, we expected to detect differing signals in the degree of incomplete lineage sorting in gene trees for m. griseorufus and m. murinus.
for example, some coalescence times within m. murinus would predate speciation, whereas m. griseorufus would show a higher degree of reciprocal monophyly among different loci due to its smaller effective population size, particularly during the bottleneck that was hypothesized to be associated with the peripatric speciation event.
conversely, under an allopatric speciation model we would expect the degree of incomplete lineage sorting among gene trees to be similar for both species.
thus, we performed maximum likelihood  phylogenetic analysis of each locus using raxml v.  <dig> .
<dig>  <cit> .
because all mitochondrial genes are linked, we performed a single ml analysis for the concatenated mtdna loci.
for each gene we ran a full ml analysis followed by rapid bootstrapping  <cit>  using the automre bootstopping criterion.
RESULTS
sequence diversity and effective population sizes
in general, sequence diversity characteristics showed moderate values for each gene for both species .
although average values were slightly higher, m. murinus did not exhibit consistently higher diversity values than m. griseorufus based on haplotype diversity, nucleotide diversity, or average number of nucleotide differences.
multiple runs of bpp gave similar results indicating adequate sampling of the posterior.
ess values were also high for all parameters .
using the prior g for θ and τ <dig>  mean effective population size was slightly greater for m. murinus  versus m. griseorufus .
however, 95% hpds for the two species almost completely overlapped .
the mrca had a significantly smaller mean population size than either species .
bpp results using a gamma prior of g for θ and τ <dig> resulted in slightly larger parameter estimates .
with these priors, the mean effective population size of the mrca was intermediate between m. murinus and m. griseorufus.
.
however, confidence intervals again overlapped substantially.
in all analyses, effective population size estimates of the mrca were more similar to m. griseorufus than they were to m. murinus.
s
h
hd
pi
k
demographic expansion
for all *beast analyses, the 95% hpd for the coefficient of variation parameter included zero for all relaxed clock analyses, indicating that a strict clock was sufficient to explain the data.
mean estimated rates of nucleotide substitution  and 95% hpds for each locus were as follows: mtdna =  <dig>  ; enol =  <dig>  ; fib =  <dig>  ; vwf =  <dig>  ).
ebsp results for m. murinus indicated an increase in effective population size at about 1 ma, with a subsequent rapid decline in size starting approximately  <dig> ka and continuing to the present .
the mean number of population size changes throughout the history of m. murinus was estimated as  <dig>  .
conversely, results for m. griseorufus showed signs of relatively constant population size through time  with an estimated mean of  <dig>  size changes .
phylogenetic analysis
maximum likelihood phylogenetic analyses of the mtdna data revealed reciprocal monophyly for both species with strong bootstrap support .
gene copies within microcebus griseorufus were also monophyletic based on ml analysis of the enol locus .
conversely, none of the nuclear markers showed reciprocal monophyly for m. murinus.
moderate geographic structure among populations was detected from the mtdna analysis, particularly for m. murinus .
however, nuclear gene tree analyses suggested little to no signal of geographic population structure within species as haplotypes were shared among many localities .
species tree analyses suggested that divergence of m. griseorufus and m. murinus occurred approximately 3–6 ma.
discussion
speciation models
our motivation for this study was to differentiate between competing models of speciation that can potentially explain the historical divergence between a sister species pair of mouse lemurs.
these primates are of particular interest for such a study given their cryptic species diversity, highly threatened status, and their endemic distribution in madagascar, one of earth’s hottest biodiversity hotspots.
given the difference in the sizes of their geographic distributions , we initially hypothesized that m. griseorufus diverged peripatrically from a geographically widespread common ancestor .
the observation that two sister species occur with different yet overlapping ranges in madagascar raises obvious questions regarding the driving mechanisms behind their divergence.
given that one species, m. griseorufus, shows a more limited though contiguous range with m. murinus is on its face entirely congruent with a peripatric model of speciation.
to test this hypothesis, we formed a set of predictions that rely on a series of historical demographic variables including comparisons of effective population size in both the modern species and their common ancestor, as well as patterns of population size change in the history of the focal species.
the majority of our results reject a model of peripatric speciation in favor of a model of allopatric divergence with subsequent range expansion for m. murinus.
contemporary and ancestral ne
our bpp results using both large and small priors for divergence times and population sizes strongly suggest that contemporary ne values are similar for both species.
for example, although mean estimates of ne for m. murinus are slightly larger than those for m. griseorufus, 95% hpds overlapped significantly.
in analyses using the priors g, ne estimates of the mrca are substantially smaller than for either contemporary species, although using larger priors results in ne estimates of the mrca that are intermediate to the two contemporary species.
however, regardless of the prior used, ne estimates for the mrca were more similar to m. griseorufus than to the more widespread m. murinus.
these results support a model in which the mrca was a species with a relatively small geographic distribution likely situated within southern madagascar.
these conclusions are congruent with our phylogenetic results and with estimates of genetic diversity for both species.
although genetic diversity is slightly higher in m. murinus, this may be an artifact of both our sampling regime and/or recent demographic trends for this species .
under a peripatric scenario we would expect to find substantially larger diversity values in m. murinus, a pattern that was not recovered with any marker.
the reciprocal monophyly and population structure we observed with the mtdna and incomplete lineage sorting in nuclear markers is also congruent with other recent studies of these species .
these results strongly support the allopatric rather than the peripatric speciation model.
of particular importance for this conclusion is the estimation of ne for the mrca of m. murinus and m. griseorufus.
the model employed in bpp allows for the combined analysis of multiple genetic markers in a coalescent framework, a necessary approach as individual loci may suffer from rate heterogeneities and idiosyncratic gene genealogies  <cit> .
because a single genetic locus provides a limited and highly stochastic perspective, we find that examining four independent loci is effective for estimates when combined, but may remain insufficient for reconciling historical demographic processes when analyzed individually given the limited power of single locus analyses  <cit> .
a rejection of a peripatric model for mouse lemurs is also similar to a recent study on mantellid frogs based on range overlap analysis, where the authors found that range size differences among sister species increased with evolutionary age  <cit> .
for example, many sister species of frogs were composed of microendemics encompassing similarly-sized geographic ranges.
under a peripatric scenario, range asymmetry would be high in younger species.
demographic changes
our results regarding population size changes provide further support for the allopatric model.
although only four independent loci are used to infer demographic changes, four loci are predicted to reduce the error by one-half as compared to single-locus estimates  <cit> .
ebsp results suggest demographic expansion for m. murinus in the quaternary around 1 ma.
demographic expansion during the quaternary has also been documented previously for populations of m. murinus in northwest madagascar  <cit>  as well as for malagasy rodents  <cit> .
however, our results differ from previous single-locus studies of m. murinus that have utilized different analytical methods and found evidence for more recent episodes of expansion during the pleistocene and holocene  <cit> .
evidence from palynological records indicates that the pleistocene climate and vegetation of madagascar, like most of the world, was quite different than the climate of today  <cit> .
for example, during the last glacial maximum  ~40– <dig>  years ago, humid forest was likely restricted to isolated refugia scattered throughout the island, whereas dry, xeric vegetation was allowed to expand.
indeed, a wide body of evidence is available that suggests that a large portion of madagascar experienced substantially drier conditions during the lgm than the present .
both m. murinus and m. griseorufus are common in dry environments and both were likely affected by climate conditions associated with quaternary madagascar.
following the demographic expansion of m. murinus during the quaternary, our results suggest recent and substantial population decline of this species beginning approximately  <dig> ka and continuing to the present.
it has been proposed that vegetation shifts associated with pleistocene climate change were more substantial in western dry forest versus the arid spiny forest to the south  <cit> , which may partly explain why no evidence of recent population decline is indicated for m. griseorufus.
furthermore, evidence suggests that humans first colonized madagascar ~ <dig>  years ago and subsequently had a rapid and profound impact on the native biota and their habitats  <cit> .
a variety of hypotheses have been put forth to explain the decline of malagasy flora and fauna subsequent to human colonization, including increased frequency of fire  <cit> , drought  <cit> , hunting  <cit> , invasive species  <cit> , disease  <cit> , and synergistic anthropogenic influences  <cit> .
regardless of the exact mechanism, it is highly probable that the recent and rapid population decline inferred from our data for m. murinus has been exacerbated by subsequent anthropogenic influences beginning around  <dig>  years ago.
the relatively constant population size of m. griseorufus suggests that human impacts and habitat fragmentation throughout the southern spiny forests may have been less severe than impacts throughout western dry forests.
alternative speciation hypotheses
our results indicate that the large geographic range of m. murinus seems to be a uniquely derived feature of this species.
this phenomenon begs investigation.
even so, the evidence presented in this study fails to specifically explain the mechanism either promoting the geographic range expansion of m. murinus or limiting the range of m. griseorufus.
the initial divergence between the two species may have resulted from a geographically-based vicariant event or an ecological niche separation.
there is no discernible extrinsic barrier separating the two species, but there is evidence that suggests ecological segregation  <cit> .
although some studies suggest that m. murinus preferentially inhabits dry forest habitat in northwestern madagascar  <cit> , it could be argued that m. murinus is a generalist species as it is often found in both dry deciduous forest and wet, gallery forest habitats in the southeast  <cit> .
these habitat associations contrast with m. griseorufus, which is more common in xeric, spiny forest  <cit> , though recently m. griseorufus has been shown to inhabit both spiny and gallery forests at the beza mahafaly private reserve  <cit> .
therefore, while these habitat preferences may have been important in the original subdivision of the ancestral species, there is no direct evidence supporting the hypothesis that this drove their initial divergence.
the use of next-generation dna sequencing methods, in particular rad-seq or whole transcriptomes, may be useful for investigating and quantifying genomic islands of divergence and adaptation in this system.
additionally, multilocus demographic methods such as those used here should be combined with future projections of species distribution models to better inform conservation practices.
this is of utmost importance for highly threatened taxa inhabiting areas that are experiencing high rates of habitat loss, such as the case with madagascar’s lemur fauna  <cit> .
CONCLUSIONS
we estimated historical demographic parameters in a multilocus coalescent framework to test the predictions associated with two models of speciation that may have driven the divergence of m. griseorufus and m. murinus.
the majority of our results reject the hypothesis of peripatric speciation.
our results instead favor a model of allopatric divergence from a range-restricted common ancestor in southwestern madagascar, with subsequent range expansions for m. murinus.
whether due to ecological constraint or interspecific competition, m. griseorufus is presently restricted to the arid spiny forest in the south, whereas m. murinus has successfully expanded throughout much of western madagascar and limited areas to the southeast.
the methods used here can be easily applied to address similar evolutionary questions in other systems to help elucidate the geographic context of divergence and speciation.
in turn, these approaches can help guide conservation priorities when synthesized with complex geospatial methods and species distribution models.
abbreviations
ne: effective population size; mtdna: mitochondrial dna; ndna: nuclear dna; mrca: most recent common ancestor; enol: alpha enolase intron; fib: alpha fibrinogen intron; vwf: von willebrand factor intron; cytb: cytochrome b; coii: cytochrome c oxidase subunit ii; bpp: bayesian phylogenetics and phylogeography; ma: million years ago; ess: effective sample size; dscr: duke shared cluster resource; ml: maximum likelihood; hpd: highest posterior density; θm: effective population size of microcebus murinus; θg: effective population size of microcebus griseorufus; θmg: effective population size of the most recent common ancestor; τmg: divergence time of the most recent common ancestor; ka: thousand years ago; lgm: last glacial maximum.
competing interests
the authors declare that they have no competing interests.
authors’ contributions
cb, kh, ar, ady designed the study.
cb collected and analyzed the data.
cb, ar, ady wrote the paper.
all authors read an approved the final manuscript.
authors’ information
cb is presently a postdoctoral associate in the laboratory of anne yoder at duke university.
his research interests include combining next-generation dna sequencing data with spatially explicit modeling to understand the contemporary and historical processes responsible for shaping patterns of genomic variation in natural populations.
he is also interested in using next-generation sequencing of environmental samples to quantify patterns of biodiversity.
supplementary material
additional file 1: table s1
database of all individuals and haplotypes used for this study.
the last four columns indicate which haplotypes were used for each analysis.
click here for file
additional file 2: table s2
results from bpp analysis using priors g for both τ <dig> and θ. parameters were estimated on a fixed species tree ).
click here for file
additional file 3: table s3
results from bpp analysis using priors g for both τ <dig> and θ. parameters were estimated on a fixed species tree ).
click here for file
additional file 4: figure s1
maximum likelihood mtdna gene tree  for all microcebus griseorufus  and m. murinus  sequences used for this study.
values at nodes represent bootstrap support values > <dig> calculated using the automre function in raxml.
for ease of visualization the outgroup taxon  is removed.
click here for file
additional file 5: figure s2
maximum likelihood alpha enolase gene tree for all microcebus griseorufus  and m. murinus  sequences used for this study.
values at nodes represent bootstrap support values > <dig> calculated using the automre function in raxml.
for ease of visualization the outgroup taxon  is removed.
click here for file
additional file 6: figure s3
maximum likelihood alpha fibrinogen gene tree for all microcebus griseorufus  and m. murinus  sequences used for this study.
values at nodes represent bootstrap support values > <dig> calculated using the automre function in raxml.
for ease of visualization the outgroup taxon  is removed.
click here for file
additional file 7: figure s4
maximum likelihood von willebrand factor gene tree for all microcebus griseorufus  and m. murinus  sequences used for this study.
values at nodes represent bootstrap support values > <dig> calculated using the automre function in raxml.
for ease of visualization the outgroup taxon  is removed.
click here for file
acknowledgments
this research was supported by nsf award deb- <dig> and a biodiversity leadership award from the bay & paul foundations to a.d.y.
the manuscript was significantly improved thanks to comments from m. noor and discussion with j. wall.
the data sets supporting the results of this article are available in the dryad repository, .

 in recent years, there has been growing interest in using genealogical tree structure to reconstruct the demographic and temporal context of diverging populations and species. this tree-based approach permits the examination of lineages as they have diverged in the past using principles derived from coalescent theory, since ancestral polymorphisms that are shared among lineages are sorted during population segregation and speciation. a comparative demographic approach utilizing tree-based methods has been implemented in a variety of biological sub-fields, including systematics, phylogeography, conservation, and life history. further, estimates of divergence times, migration rates, and effective population sizes among groups of organisms are now inferred with consideration of the genealogical structure of multilocus data sets. these estimates may be synthesized and compared among sister taxa and used to infer the impact of the biogeographic or climatic context or other historical processes on patterns of genetic differences among groups. in general, it is assumed that a complex suite of geologic, climatic, and population genetic forces have led to divergence and speciation in allopatry, where a reproductive barrier effectively divides an ancestral species into two populations of roughly equal size. with the continual progress of next-generation dna sequencing and the increasing abundance of sequence data from across the genome, researchers are beginning to test alternate models of speciation including parapatric divergence along ecological clines. few studies, however, have used recently developed multilocus coalescent methods to fully understand the geography of speciation and the frequency of peripatric events in nature. peripatric speciation can be seen as a subset of allopatric speciation, with a founder effect leading to the formation of a new species as a small population becomes physically separated from a broadly distributed ancestor. by definition then, the effective population size of the diverging population will be smaller than that of the common ancestor in a peripatric scenario. although evidence of peripatric divergence is commonly inferred throughout archipelagos using traditional phylogenetic approaches, examples from strictly mainland taxa are relatively rare. new algorithms can effectively estimate demographic parameters on a species tree to help disentangle the historical context of divergence and speciation and to test the plausibility of a peripatric scenario in comparison with the generally accepted null model of allopatric divergence. mitochondrial dna has been commonly used in studies of speciation given its small effective population size relative to nuclear dna loci. accordingly, mitochondrial gene trees undergo lineage sorting much more rapidly, and thus have a higher probability of yielding gene trees that are reciprocally monophyletic for the hypothesized species. however, the stochasticity inherent in single-locus gene trees makes this a controversial approach for estimating species histories, as previous authors have discussed. therefore, the implementation of a comparative demographic approach incorporating multiple unlinked loci is now regarded as a more powerful approach. mouse lemurs are small, nocturnal primates that inhabit a wide range of habitats throughout madagascar. as many as species have been described within the genus, thus a remarkable increase from the two species taxonomy that was recognized as recently as. the geography of the island has historically been considered to play a significant role in explaining divergence patterns among lemurs and other vertebrate species and specifically microcebus species. within the genus, however, m. murinus is unique in having a relatively vast geographic distribution. whereas all other microcebus species have geographically limited ranges --- some restricted to isolated forest fragments --- m. murinus is distributed along most of the western half and along the south of the island. over the extensive range of m. murinus, the species has been found to overlap with multiple congeners, many of which exhibit a high degree of endemicity. one example of this pattern of sympatry couples m. murinus with m. griseorufus, with multiple lines of evidence for distribution overlap and potential hybridization between the two species in southern madagascar. within their respective ranges, m. murinus and m. griseorufus also exhibit marked differences in abundance. microcebus murinus is a common, widespread species throughout the west that has been studied extensively since its description. conversely, m. griseorufus has been recorded from relatively few locations within a narrow range in the southwest of the island in the recent past. despite their partially overlapping distributions and large genetic distances, m. murinus and m. griseorufus are sister taxa. given the current ranges of m. murinus and m. griseorufus and their status as sister species, we consider two possible scenarios for the speciation event within this small clade. one scenario is that m. griseorufus was derived from the most recent common ancestor as a small founding population, with divergence the result of random genetic drift. alternatively, we consider that the mrca was a small, narrowly distributed population located in southwestern madagascar. under this alternative hypothesis, the two modern species resulted from an event that partitioned that ancestral species into comparably sized populations prior to a substantial range expansion in m. murinus. these two hypotheses differ most critically in their assumptions regarding the ancestral range and population size: was the mrca more like modern m. murinus or more like modern m. griseorufus? to test these alternative scenarios, we analyze multilocus dna sequence data with coalescent-based methods to estimate relevant demographic parameters. if m. griseorufus diverged as a small founding population from a large ancestral population, then we expect 1) that m. murinus has a greater effective population size and higher levels of genetic diversity than m. griseorufus, 2) that estimates of effective population size for the mrca are more similar to those for m. murinus than to those for m. griseorufus, 3) that historical changes in population size in m. murinus were minimal, and finally, 4) that m. griseorufus is more likely to segregate as a clade in individual gene trees. under the alternative hypothesis of recent allopatric speciation, we expect that 1) that m. murinus and m. griseorufus have similar effective population sizes and levels of genetic diversity, 2) that ne estimates for the mrca are smaller than those for contemporary species or similar to m. griseorufus, 3) that significant population expansion has occurred in m. murinus, and finally, 4) that individual gene trees are more likely to show incomplete lineage sorting between the two species. methods data assembly to estimate relevant demographic parameters for both species, we mined genbank for all available high-quality sequences as both m. murinus and m. griseorufus have been used in a number of recent molecular studies. we sought to utilize genes with sufficient sequence variation to reliably estimate coalescent parameters. we also chose data sets with samples covering a large portion of the known geographic range of both species to assemble the most geographically complete data set possible. although recent studies have suggested that m. murinus may possibly contain at least three undocumented species,, we treat these populations as m. murinus sensu lato for the purpose of this study. the final data set consisted of the following nuclear loci: alpha enolase intron, alpha fibrinogen intron, von willebrand factor intron. although the adora gene has been sequenced numerous times within microcebus, levels of variation were too low to provide any meaningful information for this study. we also assembled and concatenated information from two mitochondrial loci including cytochrome b and cytochrome c oxidase subunit ii for a total of  bp. the total number of sequences per locus was as follows: enol =  sequences; fib =  sequences; vwf =  sequences; mtdna =  concatenated sequences. all of the nuclear sequences represented phased haplotypes, which were included for all subsequent analyses. cheirogaleus major was used as an outgroup for all phylogenetic and demographic analyses. although cheirogaleus can be considered a relatively distant outgroup, we were interested in adding a temporal component to our analyses and sought to maintain consistency with taxon sampling across analyses. all sequence manipulation was performed in geneious v... multiple sequence alignments were conducted using mafft v. within geneious. sequence diversity and effective population sizes we first tested each marker for signs of recombination using the program rdp v.. each marker was tested for recombination events using the geneconv, maxchi, and rdp methods. default settings were used for all analyses. because no signs of recombination were detected, all subsequent analyses utilized the entire read of each locus. general sequence diversity statistics for each locus and species including number of haplotypes, haplotype diversity, number of segregating sites, and nucleotide diversity were calculated using dnasp v... using the species tree topology, cheirogaleus major) prescribed by both nuclear and mitochondrial phylogenies, we used the program bayesian phylogenetics and phylogeography to estimate coalescent-scaled population sizes and time of divergence. this method accommodates the species tree divergence patterns as well as gene tree lineage sorting processes. a gamma prior g, with mean 2/1000 =  was used for the population size parameters. the age of the root in the species tree was also assigned a gamma prior g, while the other divergence time parameters were assigned the dirichlet prior. locus-specific mutation rates were allowed to vary, and we specified a heredity multiplier value g to account for the combined mtdna and ndna data. we also performed analyses using larger priors for both τ and θ ) to evaluate the sensitivity of our results to the choice of prior. for this study, the parameters of interest were θm, θg, θmg, and τmg, where m = m. murinus, g = m. griseorufus, and mg = the mrca of m. murinus and m. griseorufus. the default numbers of generations and sampling intervals were used for all analyses. due to computational issues with the full data set, all bpp analyses were implemented on a reduced data set consisting of approximately sequences per locus per microcebus species and one to six cheirogaleus sequences per locus. sequences were randomly sampled to encompass a broad geographic range for each species. previous simulation-based studies have shown that similar sample sizes can be sufficient to infer speciation processes. additional analyses on further reduced data sets yielded similar results. all analyses were run at least twice to check for consistency among runs. demographic expansion to further test between the alternative speciation hypotheses, and to determine if the relatively large geographic range of m. murinus was due to population expansion following the split with m. griseorufus, we tested for signs of demographic change through time for both species. under a peripatric speciation model, an ancestral species with a large range similar to that of modern m. murinus would have diverged into m. murinus and m. griseorufus populations, with the former experiencing little to no demographic change associated with the speciation event and the latter experiencing a significant population bottleneck event. conversely, under the allopatric speciation model, a narrowly-distributed ancestral species would have diverged into m. murinus and m. griseorufus populations, with the former subsequently experiencing a significant population growth event and the latter experiencing a relatively smaller but possibly detectable demographic change after the speciation event. to test for signs of demographic expansion or contraction we implemented multilocus bayesian methods. specifically, we constructed extended bayesian skyline plots for both species using beast v... studies have shown that doubling the number of independent loci can reduce error and 95% credible intervals in demographic reconstruction by √. thus, compared to single locus estimates, our multilocus data provided a powerful approach for estimating demographic trends in mouse lemurs. first, we used jmodeltest v.. to calculate model likelihood scores for each locus and to estimate optimal models using bic. because mtdna is inherited as a single linked unit, and to minimize the computational burden for beast, the concatenated mtdna data were treated as a single partition. the three nuclear loci were each specified as a separate partition for model fitting. we tested the likelihood of commonly used models in beast. the first two blocks represent per locus models selected for extended bayesian skyline plot analyses. ‘*beast’ represents a reduced data set encompassing both species and the outgroup. we were also interested in adding a temporal component to the demographic analyses. although there is no fossil record for lemurs, recent studies have utilized fossil information from more distantly related groups to date divergence times within the lemurs. using multiple calibration points outside the clade and relaxed clock methods, these authors estimated the split between cheirogaleus and microcebus to be approximately million years ago. we first used this information to estimate the substitution rate for each of the test loci in beast using a reduced data set of both species and cheirogaleus. best-fitting models were calculated and used for all divergence dating and rate estimation. because alignments contained representative alleles from both multiple species and multiple individuals within species, we used *beast to estimate the posterior distribution of substitution rates. we defined three species and grouped alleles accordingly. for each analysis, the root node of the species tree was calibrated with a normal distribution around a mean of 25 ma and standard deviation of 5 ma, which encompassed the 95% hpd estimates from previous studies. to increase computational efficiency, we ran four independent *beast analyses to estimate the posterior distribution of rates. both strict and relaxed clock models were tested. all analyses were run between 10– million generations; sampling was chosen at intervals to utilize draws from the posterior. following analyses, the program tracer v was used to examine effective sample size values and examine the posterior distribution of relevant parameters. we used the 95% hpd of the substitution rate for each locus as a uniform prior to add a temporal component to all ebsp analyses. all ebsp analyses used a strict clock. operators were modified according to author recommendations and analyses were run for million generations or million generations to obtain adequate ess values. all ebsp and *beast analyses were implemented via the duke shared cluster resource. phylogenetic analysis under a scenario of peripatric speciation, we expected to detect differing signals in the degree of incomplete lineage sorting in gene trees for m. griseorufus and m. murinus. for example, some coalescence times within m. murinus would predate speciation, whereas m. griseorufus would show a higher degree of reciprocal monophyly among different loci due to its smaller effective population size, particularly during the bottleneck that was hypothesized to be associated with the peripatric speciation event. conversely, under an allopatric speciation model we would expect the degree of incomplete lineage sorting among gene trees to be similar for both species. thus, we performed maximum likelihood phylogenetic analysis of each locus using raxml v... because all mitochondrial genes are linked, we performed a single ml analysis for the concatenated mtdna loci. for each gene we ran a full ml analysis followed by rapid bootstrapping using the automre bootstopping criterion.  sequence diversity and effective population sizes in general, sequence diversity characteristics showed moderate values for each gene for both species. although average values were slightly higher, m. murinus did not exhibit consistently higher diversity values than m. griseorufus based on haplotype diversity, nucleotide diversity, or average number of nucleotide differences. multiple runs of bpp gave similar results indicating adequate sampling of the posterior. ess values were also high for all parameters. using the prior g for θ and τ mean effective population size was slightly greater for m. murinus versus m. griseorufus. however, 95% hpds for the two species almost completely overlapped. the mrca had a significantly smaller mean population size than either species. bpp results using a gamma prior of g for θ and τ resulted in slightly larger parameter estimates. with these priors, the mean effective population size of the mrca was intermediate between m. murinus and m. griseorufus.. however, confidence intervals again overlapped substantially. in all analyses, effective population size estimates of the mrca were more similar to m. griseorufus than they were to m. murinus. s h hd pi k demographic expansion for all *beast analyses, the 95% hpd for the coefficient of variation parameter included zero for all relaxed clock analyses, indicating that a strict clock was sufficient to explain the data. mean estimated rates of nucleotide substitution and 95% hpds for each locus were as follows: mtdna = ; enol = ; fib = ; vwf =  ). ebsp results for m. murinus indicated an increase in effective population size at about 1 ma, with a subsequent rapid decline in size starting approximately ka and continuing to the present. the mean number of population size changes throughout the history of m. murinus was estimated as. conversely, results for m. griseorufus showed signs of relatively constant population size through time with an estimated mean of size changes. phylogenetic analysis maximum likelihood phylogenetic analyses of the mtdna data revealed reciprocal monophyly for both species with strong bootstrap support. gene copies within microcebus griseorufus were also monophyletic based on ml analysis of the enol locus. conversely, none of the nuclear markers showed reciprocal monophyly for m. murinus. moderate geographic structure among populations was detected from the mtdna analysis, particularly for m. murinus. however, nuclear gene tree analyses suggested little to no signal of geographic population structure within species as haplotypes were shared among many localities. species tree analyses suggested that divergence of m. griseorufus and m. murinus occurred approximately 3–6 ma. discussion speciation models our motivation for this study was to differentiate between competing models of speciation that can potentially explain the historical divergence between a sister species pair of mouse lemurs. these primates are of particular interest for such a study given their cryptic species diversity, highly threatened status, and their endemic distribution in madagascar, one of earth’s hottest biodiversity hotspots. given the difference in the sizes of their geographic distributions, we initially hypothesized that m. griseorufus diverged peripatrically from a geographically widespread common ancestor. the observation that two sister species occur with different yet overlapping ranges in madagascar raises obvious questions regarding the driving mechanisms behind their divergence. given that one species, m. griseorufus, shows a more limited though contiguous range with m. murinus is on its face entirely congruent with a peripatric model of speciation. to test this hypothesis, we formed a set of predictions that rely on a series of historical demographic variables including comparisons of effective population size in both the modern species and their common ancestor, as well as patterns of population size change in the history of the focal species. the majority of our results reject a model of peripatric speciation in favor of a model of allopatric divergence with subsequent range expansion for m. murinus. contemporary and ancestral ne our bpp results using both large and small priors for divergence times and population sizes strongly suggest that contemporary ne values are similar for both species. for example, although mean estimates of ne for m. murinus are slightly larger than those for m. griseorufus, 95% hpds overlapped significantly. in analyses using the priors g, ne estimates of the mrca are substantially smaller than for either contemporary species, although using larger priors results in ne estimates of the mrca that are intermediate to the two contemporary species. however, regardless of the prior used, ne estimates for the mrca were more similar to m. griseorufus than to the more widespread m. murinus. these results support a model in which the mrca was a species with a relatively small geographic distribution likely situated within southern madagascar. these conclusions are congruent with our phylogenetic results and with estimates of genetic diversity for both species. although genetic diversity is slightly higher in m. murinus, this may be an artifact of both our sampling regime and/or recent demographic trends for this species. under a peripatric scenario we would expect to find substantially larger diversity values in m. murinus, a pattern that was not recovered with any marker. the reciprocal monophyly and population structure we observed with the mtdna and incomplete lineage sorting in nuclear markers is also congruent with other recent studies of these species. these results strongly support the allopatric rather than the peripatric speciation model. of particular importance for this conclusion is the estimation of ne for the mrca of m. murinus and m. griseorufus. the model employed in bpp allows for the combined analysis of multiple genetic markers in a coalescent framework, a necessary approach as individual loci may suffer from rate heterogeneities and idiosyncratic gene genealogies. because a single genetic locus provides a limited and highly stochastic perspective, we find that examining four independent loci is effective for estimates when combined, but may remain insufficient for reconciling historical demographic processes when analyzed individually given the limited power of single locus analyses. a rejection of a peripatric model for mouse lemurs is also similar to a recent study on mantellid frogs based on range overlap analysis, where the authors found that range size differences among sister species increased with evolutionary age. for example, many sister species of frogs were composed of microendemics encompassing similarly-sized geographic ranges. under a peripatric scenario, range asymmetry would be high in younger species. demographic changes our results regarding population size changes provide further support for the allopatric model. although only four independent loci are used to infer demographic changes, four loci are predicted to reduce the error by one-half as compared to single-locus estimates. ebsp results suggest demographic expansion for m. murinus in the quaternary around 1 ma. demographic expansion during the quaternary has also been documented previously for populations of m. murinus in northwest madagascar as well as for malagasy rodents. however, our results differ from previous single-locus studies of m. murinus that have utilized different analytical methods and found evidence for more recent episodes of expansion during the pleistocene and holocene. evidence from palynological records indicates that the pleistocene climate and vegetation of madagascar, like most of the world, was quite different than the climate of today. for example, during the last glacial maximum ~40–  years ago, humid forest was likely restricted to isolated refugia scattered throughout the island, whereas dry, xeric vegetation was allowed to expand. indeed, a wide body of evidence is available that suggests that a large portion of madagascar experienced substantially drier conditions during the lgm than the present. both m. murinus and m. griseorufus are common in dry environments and both were likely affected by climate conditions associated with quaternary madagascar. following the demographic expansion of m. murinus during the quaternary, our results suggest recent and substantial population decline of this species beginning approximately ka and continuing to the present. it has been proposed that vegetation shifts associated with pleistocene climate change were more substantial in western dry forest versus the arid spiny forest to the south, which may partly explain why no evidence of recent population decline is indicated for m. griseorufus. furthermore, evidence suggests that humans first colonized madagascar ~  years ago and subsequently had a rapid and profound impact on the native biota and their habitats. a variety of hypotheses have been put forth to explain the decline of malagasy flora and fauna subsequent to human colonization, including increased frequency of fire, drought, hunting, invasive species, disease, and synergistic anthropogenic influences. regardless of the exact mechanism, it is highly probable that the recent and rapid population decline inferred from our data for m. murinus has been exacerbated by subsequent anthropogenic influences beginning around  years ago. the relatively constant population size of m. griseorufus suggests that human impacts and habitat fragmentation throughout the southern spiny forests may have been less severe than impacts throughout western dry forests. alternative speciation hypotheses our results indicate that the large geographic range of m. murinus seems to be a uniquely derived feature of this species. this phenomenon begs investigation. even so, the evidence presented in this study fails to specifically explain the mechanism either promoting the geographic range expansion of m. murinus or limiting the range of m. griseorufus. the initial divergence between the two species may have resulted from a geographically-based vicariant event or an ecological niche separation. there is no discernible extrinsic barrier separating the two species, but there is evidence that suggests ecological segregation. although some studies suggest that m. murinus preferentially inhabits dry forest habitat in northwestern madagascar, it could be argued that m. murinus is a generalist species as it is often found in both dry deciduous forest and wet, gallery forest habitats in the southeast. these habitat associations contrast with m. griseorufus, which is more common in xeric, spiny forest, though recently m. griseorufus has been shown to inhabit both spiny and gallery forests at the beza mahafaly private reserve. therefore, while these habitat preferences may have been important in the original subdivision of the ancestral species, there is no direct evidence supporting the hypothesis that this drove their initial divergence. the use of next-generation dna sequencing methods, in particular rad-seq or whole transcriptomes, may be useful for investigating and quantifying genomic islands of divergence and adaptation in this system. additionally, multilocus demographic methods such as those used here should be combined with future projections of species distribution models to better inform conservation practices. this is of utmost importance for highly threatened taxa inhabiting areas that are experiencing high rates of habitat loss, such as the case with madagascar’s lemur fauna.  we estimated historical demographic parameters in a multilocus coalescent framework to test the predictions associated with two models of speciation that may have driven the divergence of m. griseorufus and m. murinus. the majority of our results reject the hypothesis of peripatric speciation. our results instead favor a model of allopatric divergence from a range-restricted common ancestor in southwestern madagascar, with subsequent range expansions for m. murinus. whether due to ecological constraint or interspecific competition, m. griseorufus is presently restricted to the arid spiny forest in the south, whereas m. murinus has successfully expanded throughout much of western madagascar and limited areas to the southeast. the methods used here can be easily applied to address similar evolutionary questions in other systems to help elucidate the geographic context of divergence and speciation. in turn, these approaches can help guide conservation priorities when synthesized with complex geospatial methods and species distribution models. abbreviations ne: effective population size; mtdna: mitochondrial dna; ndna: nuclear dna; mrca: most recent common ancestor; enol: alpha enolase intron; fib: alpha fibrinogen intron; vwf: von willebrand factor intron; cytb: cytochrome b; coii: cytochrome c oxidase subunit ii; bpp: bayesian phylogenetics and phylogeography; ma: million years ago; ess: effective sample size; dscr: duke shared cluster resource; ml: maximum likelihood; hpd: highest posterior density; θm: effective population size of microcebus murinus; θg: effective population size of microcebus griseorufus; θmg: effective population size of the most recent common ancestor; τmg: divergence time of the most recent common ancestor; ka: thousand years ago; lgm: last glacial maximum. competing interests the authors declare that they have no competing interests. authors’ contributions cb, kh, ar, ady designed the study. cb collected and analyzed the data. cb, ar, ady wrote the paper. all authors read an approved the final manuscript. authors’ information cb is presently a postdoctoral associate in the laboratory of anne yoder at duke university. his research interests include combining next-generation dna sequencing data with spatially explicit modeling to understand the contemporary and historical processes responsible for shaping patterns of genomic variation in natural populations. he is also interested in using next-generation sequencing of environmental samples to quantify patterns of biodiversity. supplementary material additional file 1: table s1 database of all individuals and haplotypes used for this study. the last four columns indicate which haplotypes were used for each analysis. click here for file additional file 2: table s2 results from bpp analysis using priors g for both τ and θ. parameters were estimated on a fixed species tree ). click here for file additional file 3: table s3 results from bpp analysis using priors g for both τ and θ. parameters were estimated on a fixed species tree ). click here for file additional file 4: figure s1 maximum likelihood mtdna gene tree for all microcebus griseorufus and m. murinus sequences used for this study. values at nodes represent bootstrap support values > calculated using the automre function in raxml. for ease of visualization the outgroup taxon is removed. click here for file additional file 5: figure s2 maximum likelihood alpha enolase gene tree for all microcebus griseorufus and m. murinus sequences used for this study. values at nodes represent bootstrap support values > calculated using the automre function in raxml. for ease of visualization the outgroup taxon is removed. click here for file additional file 6: figure s3 maximum likelihood alpha fibrinogen gene tree for all microcebus griseorufus and m. murinus sequences used for this study. values at nodes represent bootstrap support values > calculated using the automre function in raxml. for ease of visualization the outgroup taxon is removed. click here for file additional file 7: figure s4 maximum likelihood von willebrand factor gene tree for all microcebus griseorufus and m. murinus sequences used for this study. values at nodes represent bootstrap support values > calculated using the automre function in raxml. for ease of visualization the outgroup taxon is removed. click here for file acknowledgments this research was supported by nsf award deb- and a biodiversity leadership award from the bay & paul foundations to a.d.y. the manuscript was significantly improved thanks to comments from m. noor and discussion with j. wall. the data sets supporting the results of this article are available in the dryad repository,.
coalescent effective population size estimates were similar for both extant species and larger than that estimated for their most recent common ancestor. gene tree results show similar levels of incomplete lineage sorting within species with respect to the species tree, and locus-specific estimates of genetic diversity are concordant for both species. if the peripatric speciation model is supported, we predict substantially greater genetic diversity in m. murinus, relative to m. griseorufus. further, we expect a larger effective population size in m. murinus and in the common ancestor of the two species than in m. griseorufus, with a concomitant decrease in gene tree/species tree incongruence in the latter and weak signs of demographic expansion in m. murinus. given that these two species are sister taxa, the asymmetric but overlapping geographic ranges are consistent with a model of peripatric speciation. our results reject a model of peripatric divergence. results suggest that speciation occurred in allopatry from a common ancestor narrowly distributed throughout southwest madagascar, with subsequent range expansion for m. murinus. here, we examine speciation patterns in a pair of sister species of mouse lemur, microcebus murinus and m. griseorufus. the mouse lemurs of madagascar are a system in which evolution has generated a large number of species over a relatively recent time frame. these two species have ranges that are disparately proportioned in size, with m. murinus showing a much more extensive range that marginally overlaps that of m. griseorufus. multilocus demographic analyses suggest range expansions for m. murinus, with this species also experiencing more recent population declines over the past thousand years. population decline in m. murinus is likely related to patterns of climate change in madagascar throughout the pleistocene, potentially exacerbated by continual anthropogenic perturbation. debate continues as to whether allopatric speciation or peripatric speciation through a founder effect is the predominant force driving evolution in vertebrates. coalescent methodsmodes of speciationhistorical demographylemur evolutionmicrocebusmultilocusperipatric speciation
 debate continues as to whether allopatric speciation or peripatric speciation through a founder effect is the predominant force driving evolution in vertebrates. the mouse lemurs of madagascar are a system in which evolution has generated a large number of species over a relatively recent time frame. here, we examine speciation patterns in a pair of sister species of mouse lemur, microcebus murinus and m. griseorufus. these two species have ranges that are disparately proportioned in size, with m. murinus showing a much more extensive range that marginally overlaps that of m. griseorufus. given that these two species are sister taxa, the asymmetric but overlapping geographic ranges are consistent with a model of peripatric speciation. to test this hypothesis, we analyze dna sequence data from four molecular markers using coalescent methods. if the peripatric speciation model is supported, we predict substantially greater genetic diversity in m. murinus, relative to m. griseorufus. further, we expect a larger effective population size in m. murinus and in the common ancestor of the two species than in m. griseorufus, with a concomitant decrease in gene tree/species tree incongruence in the latter and weak signs of demographic expansion in m. murinus.  our results reject a model of peripatric divergence. coalescent effective population size estimates were similar for both extant species and larger than that estimated for their most recent common ancestor. gene tree results show similar levels of incomplete lineage sorting within species with respect to the species tree, and locus-specific estimates of genetic diversity are concordant for both species. multilocus demographic analyses suggest range expansions for m. murinus, with this species also experiencing more recent population declines over the past thousand years.  results suggest that speciation occurred in allopatry from a common ancestor narrowly distributed throughout southwest madagascar, with subsequent range expansion for m. murinus. population decline in m. murinus is likely related to patterns of climate change in madagascar throughout the pleistocene, potentially exacerbated by continual anthropogenic perturbation. genome-level data are needed to quantify the role of niche specialization and adaptation in shaping the current ranges of these species. coalescent methodsmodes of speciationhistorical demographylemur evolutionmicrocebusmultilocusperipatric speciation 