Signals of recent spatial expansions in the grey mouse lemur (Microcebus murinus)
© Schneider et al; licensee BioMed Central Ltd. 2010
Received: 4 December 2009
Accepted: 22 April 2010
Published: 22 April 2010
Pleistocene events have shaped the phylogeography of many taxa worldwide. Their genetic signatures in tropical species have been much less explored than in those living in temperate regions. We analysed the genetic structure of a Malagasy primate species, a mouse lemur with a wide distribution (M. murinus), in order to investigate such phylogeographic processes on a large tropical island. We also evaluated the effects of anthropogenic pressures (fragmentation/deforestation) and natural features (geographic distance, rivers) on genetic structure in order to complement our understanding of past and present processes of genetic differentiation.
The analysis of the mitochondrial D-loop sequences of 195 samples from 15 study sites (10 from a continuous forest and five from isolated forest fragments) from two adjacent Inter-River-Systems (IRSs) revealed that forest fragmentation and the river restrict gene flow, thereby leading to an increased genetic differentiation between populations beyond the effect of isolation-by-distance. Demographic simulations detected signals of two successive spatial expansions that could be preliminarily dated to the late Pleistocene and early Holocene. The haplotype network revealed geographic structure and showed deep molecular divergences within and between the IRSs that would be congruent with a two-step colonization scenario.
This study supports the hypothesis of a relatively recent spatial expansion of the grey mouse lemur in northwestern Madagascar, which may also explain why this taxon, in contrast to its congeners, has not yet undergone allopatric speciation in the studied area and possibly across its presently wide range.
Pleistocene events have shaped the phylogeography of many taxa worldwide . Particular emphasis has been put on these dynamics in temperate regions . Tropical biomes have also been shown to have undergone complex ecological dynamics following worldwide Pleistocene climate changes , but their genetic signatures in present-day natural populations are much less explored. The genetic structure of natural populations is not only shaped by these ancient processes but also by ongoing natural (e.g. migration, genetic drift) and more recent anthropogenic factors .
One model region for the study of the combined effects of ancient and recent factors is Madagascar. This island has been isolated from all other landmasses during the last 80 million years (My) and different vertebrate groups have undergone complex phylogeographic histories . Moreover, Madagascar has been strongly affected by anthropogenic disturbances during the last 2000 years following the arrival of humans on the island [6, 7]. Approximately 90% of the original vegetation is believed to have already disappeared (, but see  for less extreme estimates) and most forests are now heavily fragmented [10, 11]. This poses a particular threat to forest-dependent species, since barriers such as novel savannahs can significantly modify their genetic structure . Such barriers can notably reduce gene flow between populations which can lead to a loss of alleles, an increase in homozygosity and inbreeding through isolation .
Lemurs are an endemic mammal group on Madagascar. Within lemurs, mouse lemurs (Microcebus spp.) form an exceptionally diverse genus. The number of described Microcebus species has increased from four to 18 in recent years (overview in [14, 15]) due to largely increased sampling efforts and to the application of phylogenetic analyses of DNA sequences. Mouse lemurs are small, nocturnal, solitary foragers that inhabit a large variety of forest habitats with no more than two species co-occurring in a given area (see  for a review).
The high species diversity of lemurs in general and of mouse lemurs in particular has been explained as a joint effect of topographic barriers (large rivers, mountains) and/or Pleistocene climatic and vegetation changes [17–19]. Recent studies indicate that most large rivers act as genetic barriers for mouse lemurs [19, 20] and it has been suggested that a long-term separation between adjacent Inter-River-Systems (IRSs) promoted speciation in this clade . This process may explain the very limited distributions of many species.
In contrast to most other lemur species, the grey mouse lemur (M. murinus) has a wide distribution as it inhabits the dry deciduous forests from southern to northwestern Madagascar . The apparent lack of speciation events within this taxon could be the result of a lower speciation rate, a higher migratory potential, or a relatively recent expansion of this species into the western IRSs of Madagascar. The latter hypothesis is supported by three lines of evidence and shall be further investigated in this study: i) In a recent study Kappeler et al.  estimated the age of the M. murinus clade to about 1.4 My, which sets an upper Pleistocene limit for its expansion; ii) M. murinus possesses an ecological preference for dry habitats  which are typical for southern Madagascar; iii) M. murinus is the closest relative of M. griseorufus from southern Madagascar [19, 23] which suggests that this expansion may have originated in southern Madagascar. This geographic setting would indeed suggest a late Pleistocene or even postglacial colonization of western to northwestern Madagascar by the grey mouse lemur.
This study aims to assess the effects of phylogeographic processes, anthropogenic pressures (fragmentation/deforestation), and natural features (geographic distance, rivers) on the genetic structure of the grey mouse lemur at its northern distribution edge. Populations in this area are predicted to show signs of a Pleistocene spatial expansion from the likely southern origin towards the northwestern IRSs. Spatial simulations were used to identify the expansion model which fitted best to the patterns of molecular diversity observed in extant M. murinus populations. The results are compared to previously published data on M. ravelobensis, which may have arrived before. By means of spatial simulations, we could detect signals of two successive spatial expansions of M. murinus in northwestern Madagascar which could be preliminarily dated to the late Pleistocene and the early Holocene.
A total of 367 monomorphic (81.9%) and 81 polymorphic positions (18.1%) were identified in the D-loop sequences (455 bp) of the 195 individuals. A total of 47 different haplotypes were found which differed in 1-37 positions in pairwise comparisons.
Fragment size, sample size, genetic diversity and results of neutrality tests for the different sites
(abbreviation in bold)
Fragment size (km2)
S te Marie
In general, the ANP samples showed the highest genetic diversity, in particular sites with larger sample sizes (n > 5). The four samples with most haplotypes (n = 7) were located within the ANP. The three ANP samples with only one or two haplotypes were also the smallest samples with three or five individuals, respectively. Among the isolated fragments, Tanambao possessed the highest number of haplotypes (n = 6), while the other sites were less diverse. When comparing haplotype diversity, allelic richness or nucleotide diversity statistically between the sites in the ANP and the isolated forest fragments, however, no significant difference could be revealed (Mann-Whitney-U: Hd: U = 19, n1 = 10, n2 = 5, n.s.; AR(r): U = 14.5, n1 = 8, n2 = 5, n.s.; Nd: U = 17; n1 = 10, n2 = 5, n.s.). This lack of statistical evidence could be due to a lack of power given the small overall number of study sites and some variability among fragments and among ANP sites.
Table 1 shows the results of Tajima's D and Fu's Fs tests. Six out of 12 values of Tajima's D were negative, two of them significantly (Beronono, Mangatelo), while none of the populations showed significantly positive D values. Similarly, no Fu's Fs values were significant (Table 1). No tests were carried out for Ambodimanga, Komandria and Maroakata due to their lack of polymorphism.
Results of the mismatch distributions, tested for a single spatial expansion (Model 1, 1000 bootstraps)
Mismatch observed mean
Mismatch observed variance
S te Marie
Genetic differentiation between populations
Pairwise ΦST values ranged widely from 0 to 1, but most values (83 out of 91) were above 0.2 and only seven were not significant (Additional file 1: Table S1). Therefore, the sites can be generally regarded as genetically differentiated from each other. Genetic (ΦST values) and geographic distances (km) were significantly and positively correlated in the IRS1 (Mantel test, rp = 0.487, n = 13, p < 0.001) suggesting the existence of isolation-by-distance within IRS1.
Influence of the Mahajamba River
The importance of the Mahajamba River as a barrier to gene flow between IRS1 and IRS2 was tested with an AMOVA and with permutation tests. All variance components and fixation indices were large and significantly different from zero. The highest proportion of the molecular variance (41.2%) was explained by the variations between sites within each IRS while 32.8% is explained by the river (the remaining 26.0% of the molecular variance is explained by the variation within sites). These values should be interpreted with care, however, since we only had two samples in IRS2 which are almost fixed for very different haplotypes. The permutation tests showed that average ΦST values across the Mahajamba river were significantly higher than those between samples from the same side of the river (averages: 0.87 vs. 0.55, p < 0.001), but we note that some pairwise ΦST values between populations from the same side of the river were also very high (seven out of 34 were greater than 0.87, the across river average).
Influence of forest fragmentation
Genetic structure by haplotype network
This study provides evidence for two successive spatial expansions of M. murinus in northwestern Madagascar (see below) but could also demonstrate some effects of forest fragmentation, geographic distance and the large river Mahajamba on the genetic structure of this widely distributed mouse lemur species. These anthropogenic and environmental effects shall be discussed first, before we move to the phylogeographic processes.
Effects of forest fragmentation and environmental features on genetic diversity and genetic structure
Our results suggest that forest fragmentation significantly increased genetic drift in the isolated fragments. However, even though the genetic diversity of samples from the ANP was generally high, the difference with isolated fragments was not significant. This suggests that a larger dataset (both in terms of number of study sites and within-site sample sizes) would be required to draw stronger conclusions about the effect of forest fragmentation on the diversity of M. murinus. When compared to the sympatric M. ravelobensis, the number of haplotypes and the haplotype diversity per site were not significantly different (HANP: 3-8, HIFF: 1-8, HdANP: 0.542-0.880, HdIFF: 0-0.656, , four comparisons with MWU-test: n.s.). Therefore, the effect of fragmentation on genetic diversity can be assumed to be similar in both species.
Besides affecting the genetic diversity within populations, fragmentation also increased the genetic differentiation between them beyond the effects of isolation-by-distance which was also detected. The average ΦST values between sites isolated by savannah were significantly higher than those between sites from the ANP (i.e. continuous forest) despite comparable geographic distances. Furthermore, haplotype sharing was typically not observed between isolated fragments but only between sites in the ANP. The only exception was Ste Marie that shared its four haplotypes with different Park populations. It must be emphasized, though, that Ste Marie is separated from the ANP only by a 3 km stretch of savannah. Surprisingly, a previous study on M. ravelobensis in Ste Marie did not find haplotype sharing with ANP populations . These differences in their genetic diversity may be the result of divergent colonization histories or of different migration abilities during the fragmentation process.
Our study provides evidence for a significant influence of the Mahajamba River on genetic differentiation but not as high as in their reddish sister species. Guschanski et al.  found that 82.7% of the molecular variance in the reddish mouse lemurs of northwestern Madagascar was explained by the Sofia and Mahajamba Rivers, while we found that the Mahajamba River explained only 32% of the molecular variation in M murinus. Note that at that time M. ravelobensis and M. bongolavensis specimen were still believed to form two clades within the same species. D-loop haplotypes from the western or eastern side of the river Mahajamba differed in more than 59 base pairs , compared to the minimum of 24 bp differences found in this study for M. murinus. A later study showed that the two reddish clades had already diverged into two separate species and that the river had probably promoted this speciation event . This comparative evidence suggests that either the common ancestors of the reddish mouse lemur species colonized the IRS1 and IRS2 earlier than M. murinus or that M. murinus possesses larger migratory abilities that sustain gene flow between the IRSs (around the headwaters of the large rivers), counteracting genetic drift on both sides of the river. The latter explanation seems to be less likely, since the haplotypes of either side of the river have diverged completely, indicating that the river effectively prevents gene flow between the IRSs. A later arrival of M. murinus in the northwest compared to the ancestors of M. ravelobensis and M. bongolavensis would also explain why the grey mouse lemur clade did not yet undergo speciation in these two IRSs. Ultimately, this barrier function of the river may be sufficient so that populations in the two IRSs may eventually diverge into distinct species.
Taken together, the AMOVA and the network analyses suggest that the genetic differentiation in M. murinus within IRS1 and IRS2 reached levels comparable to the ones existing between the IRSs. This may indicate that the genetic lineages within the IRSs have diverged very early during the colonization process or that they result from several colonization events. These hypotheses are discussed in the next section.
Signatures of a recent colonization history
Assuming a Pleistocene colonization scenario which is suggested by its recent date of divergence (about 1.4 Mya) from its southern sister species , we predicted that populations might still show signs of the initial spatial expansion into the IRS1. The occurrence of a range expansion was confirmed by our spatial analyses, as none of the samples significantly rejected the hypothesis of spatial expansion (Model 1). Surprisingly, most samples pointed to only one of the two different periods of expansion. Our modelling approach (Model 2) showed that two successive spatial expansions in small-sized populations generate a trimodal distribution similar to the observed distributions. It may be surprising that Model 1 is not significantly rejected even when there are three modes in the mismatch distribution. This can be explained by the fact that the test is always picking up the bigger of the two last modes as a signal of spatial expansion, without considering the smaller one (see Figure 2). The estimations for the age of the expansion (Tau values, Table 2) also correspond to this largest mode. The fact that only very few Tajima's D and no Fu's Fs values were significantly negative, may seem contradictory to the mismatch distribution results. This apparent contradiction may be explained by the lower statistical power of these statistics compared to mismatch distributions to detect spatial expansions when local demes exchange only few migrants (low Nm, which seems to be the case for M. murinus) as demonstrated theoretically  and empirically in the eastern tiger salamander . This is also confirmed by our simulations of the double expansion where only 55.1% and 6.9% of the respective Tajima's D and Fu's Fs were negative and only 3.7% and 0.0% of those values were significant. Furthermore, a very recent simulation study by Städler et al.  confirmed that even when populations are undergoing a very important spatial expansion, Tajima's D can exhibit positive values (apparently indicating a bottleneck). Finally, these tests may also be affected by the recent fragmentation as noted above. Altogether, these analyses show the limitations of single statistics such as Tajima's D and Fu's Fs to detect demographic events when more than one such event has taken place [31, 32].
The fact that we pick up the signals of spatial expansions in the mismatch analyses suggests that they took place relatively recently. It has already been suggested that the Malagasy forest habitats underwent severe changes during the Pleistocene due to climatic changes associated with the glacial and interglacial periods [3, 18]. Together with the estimate of 1.4 My for the last common ancestor of M. murinus , it seems likely that M. murinus may have colonized the IRSs at its northern distribution range in the second half of the Pleistocene.
Based on the modes of the observed mismatch distributions expected under a spatial expansion model, it is possible to calculate rough estimates for the dates of these two successive spatial expansions. Assuming a constant mutation rate of 10-6 per site per generation (one estimate used for a closely related lemur species ) and using the range of modes obtained for the samples, we can calculate the time of the expansions as tau*/(2×mu), where tau* is the mode of the mismatch distribution and mu is the mutation rate for the whole sequence. The first expansion was thereby estimated to have occurred between 26,500 to 33,500 years ago (modes between 24 and 30, Figure 2). The more recent expansion would then have taken place between 3,300 and 14,000 years ago (modes between three and 13). These estimations suggest that two recent expansions took place rather recently, one before and one after the last glacial maximum (LGM). Apparently, some populations have kept memory of the older and others of the more recent expansion, as was predicted by our simulations (Figure 4).
Given the results of this study, the following scenario for the colonization of IRS1 by M. murinus is proposed: One ancestral population of M. murinus colonized the IRS1 probably before the last glacial maximum (LGM). M. murinus may then have retreated when forests contracted during the dry period associated with the LGM towards the persisting riverine forests along the Mahajamba and Betsiboka. It subsequently re-colonized the area (second wave of expansion) in a period of recent forest expansion following the LGM. This second spatial expansion is supported by the detection of a contact zone with high haplotype diversity, which parallels the events that were described in several other species as a result of postglacial expansions with subsequent admixture [34, 35]. This contact zone is running through some central sites which exhibit haplotypes from three different haplotypes clusters. This interesting and new hypothesis deserves further investigations and needs to be formally tested in the future. In particular, the results obtained with our mtDNA dataset must be confirmed and complemented by other markers, such as nuclear STRs. Moreover, more refined simulation work taking into account different factors such as environmental heterogeneity or uncertainties about the mutation rate is needed to better define and date these processes. Indeed, no reliable estimate for the mutation rate is currently available for mouse lemurs. Using a mutation rate of 10-7 , for example, would increase the age of the expansions by one order of magnitude (265,000 to 335,000 yrs for the older expansion and 33,000-140,000 yrs for the younger expansion). Although such a shift would still indicate a late Pleistocene dynamics, a reliable dating of these colonization processes would enable us to understand in more detail the evolution of Madagascar's extraordinary biodiversity and the underlying driving forces of diversification. So far, most studies have emphasized the ancient processes that triggered the many radiations that occurred on this island [36–38]. This is the first study detecting a potential genetic signature of late Pleistocene and maybe even postglacial biogeographic dynamics on this mini-continent. The further exploration of modern molecular modelling approaches may ultimately allow us to judge upon the relative importance of certain Pleistocene climate changes for the development of the biogeographic pattern that we can observe today. Mouse lemurs are an extremely suitable model to explore these processes across all Madagascar.
Study sites and samples
The study sites are located in northwestern Madagascar between two major rivers, the Betsiboka in the south and the Sofia in the north. The Mahajamba River divides this area into two IRSs (IRS1 and IRS2, Figure 1). IRS1 includes the Ankarafantsika National Park (ANP) as well as the isolated forest fragments Ste Marie, Tanambao and Mangatelo. IRS2 contains the fragments Tsinjomitondraka and Maroakata. The ANP (S 16°19', E 46°48') is one of the largest remaining forest areas in western Madagascar  and contains ten study sites (Figure 1). Two sites (Beronono, Bealana) are still connected with the ANP (Radespiel & Rakotondravony, unpublished data) and are thus considered hereafter as being part of the ANP group, even though they do not lie within the official borders of the ANP.
The forest fragments differed largely in size (2 km2 - 40.8 km2) and were isolated from each other and the ANP by variable stretches of savannah. Euclidean distances between study sites ranged from 13 to 109 km with a maximum distance of 68 km between sites within the ANP. A total of 195 individual samples (♂: 106, ♀: 88, sex undetermined: 1), collected between the years 2000-2004, were available for this study. The number of individuals per site varied between 3 and 27 (Tab. 1). Samples were collected in form of ear biopsies (2-4 mm2) after capturing mouse lemurs systematically with Sherman Live traps as described in Olivieri et al. . Double sampling of single individuals was prevented by marking all captured animals with an individual ear-cut pattern. Capturing and sampling protocols adhered to the legal requirements of Madagascar and were approved by CAFF/CORE, the "Direction Générale des Eaux et Forêts", and by ANGAP (required for the sites in the Ankarafantsika National Park).
DNA was isolated from the ear biopsies using a standard phenol/chloroform extraction . A specified PCR for mitochondrial D-loop sequences was conducted using the universal mammalian control region primers L15997 (5'-CACCATTAGCACCCAAAGCT-3') and H16498 (5'-CCTGAAGTAGGAACCAGATG-3') . The amplification of mtDNA followed the routines described in Guschanski et al. . Purification of amplified PCR-products was performed with the Invitek purification kit (MSB® Spin PCRapace) following the manufacturer's protocol. Sequencing was carried out by the company Macrogen (Seoul, South Korea, http://www.macrogen.com/eng/macrogen/macrogen_main.jsp). Sequences were analyzed with SeqMan II 6.00 (© 1989-2004 DNASTAR) and subsequently cut to a uniform length of 455 bp. The alignment of all sequences was constructed with the program Mega 4.0.
Analysis of the genetic diversity and demographic history of the populations
The number of haplotypes and the haplotype and nucleotide diversities (Nd, ) were calculated with the program DnaSP 4.0 in order to determine the genetic variability within each population. Missing data or alignment gaps were not taken into account. Allelic richness (AR(r)), corrected for sample size variations, was calculated for each population with a sample size >3 from haplotype frequencies using the rarefaction method proposed by Petit et al.  with the software RAREFAC. We used a simple regression (GLM Type II, sum of squares) to test whether the number of haplotypes or haplotype diversity could be explained by sample size. Furthermore, a Mann-Whitney U-Test was used to determine if there was a difference in genetic diversity (Hd, AR(r) and Nd) between the ANP samples and the samples from isolated forest fragments. All tests were done using the STATISTICA 5.5 software.
Tajima's D  und Fu's Fs  are known to be sensitive to departures from mutation-drift equilibrium due to population size changes (expansions, bottlenecks) and selection. They were thus computed for all samples using the Arlequin v.3 software . Following the recommendations in the manual, the Fu's Fs-statistic was regarded as significant with a p-value lower than 0.02 instead of 0.05. Negative values of both statistics point towards population growth and/or positive selection, whereas positive Tajima's D values indicate bottlenecks and/or balancing selection.
In addition to these tests, mismatch distributions were constructed for polymorphic samples with a sample size of at least 10 individuals. A mismatch distribution is the distribution of the number of nucleotide mismatches between all pairs of DNA sequences belonging to a population sample. The shape of the mismatch distribution has been shown to be influenced by past demographic events such as expansions and bottlenecks. The distribution is usually bell-shape in populations having increased demographically in the past , while it shows one or two modes in populations having passed through a range expansion, depending on the population density and the amount of migrants exchanged with neighbouring populations .
Using the software Arlequin v.3, it is possible to formally test if a mismatch distribution rejects statistically the null hypothesis of one spatial expansion, which we call Model 1 hereafter. In this model it is assumed that a population expanded spatially from one deme into a larger area divided in same-sized demes interconnected by gene flow. We estimated also estimated Tau, the time of the expansion (in mutation units, Tau = 2 Tμ).
Analysis of genetic structure
Genetic differentiation between populations was measured using ΦST. The model of Tamura & Nei  was chosen with the gamma correction factor of 0.6354, as determined with the program Modeltest 3.5 mac  on the basis of the Akaike Information Criterion (AIC). The significance of these values was estimated with 1000 permutations as implemented in Arlequin v.3. A Mantel test  was performed with XLSTAT v. 2007.7 in order to investigate the relationship between geographic and genetic distances (ΦST-values). The statistical significance was determined by means of 10,000 permutations.
An Analysis of Molecular Variance (AMOVA) was carried out in order to detect genetic structure generated by the river Mahajamba, separating the two IRSs. For this analysis, samples were divided into two groups corresponding to IRS1 und IRS2, which makes the results directly comparable to study of Guschanski et al.  on the other species of mouse lemurs in the area. The input file was constructed with Seqtrans 1.1 (Dr. S. M. Funk) and the AMOVA was performed with Arlequin v.3 . Two populations with very small sample sizes (n < 5: Bealana, Komandria) were excluded from this analysis.
The influence of the forest fragmentation on genetic differentiation was tested for the populations in IRS1 by comparing the pairwise ΦST values obtained between populations separated by stretches of savannah (ANP group vs. fragments and between fragments) and those obtained for pairs of populations separated by a "continuous" forest habitat, i.e. within the ANP group including Beronono and Bealana. Since geographic distance is a confounding factor that may inflate "savannah" distances, we only compared pairwise ΦST values for geographic distances lower than the highest one found within the ANP. Consequently, of the three fragments available, Tanambao was excluded from the analysis due to its large distance from the ANP. We regressed ΦST values against geographical distances within the ANP and then compared the residuals obtained with those corresponding to pairs involving a savannah crossing. If the savannah generated an increase in genetic differentiation we expect that the savannah residuals should be more positive than the ANP residuals (i.e. the ΦST values would be higher for a similar geographic distance). The difference between the two sets of residuals was tested using a permutation approach (10,000 permutations). We also tested for a correlation between the residuals and geographical distances. These computations were performed with R  partly with available and partly with self-written scripts.
Finally, for the analysis of genetic structure and haplotype sharing, a haplotype network was constructed with the programme TCS 1.21 . Clusters of closely related haplotypes (connected by less than 10 mutation steps) were identified and mapped geographically.
NS was an undergraduate Biology student in the lab of UR when conducting this population genetics study. LC is interested in understanding how genetic data can be used to reconstruct the demographic history of species with an emphasis on endangered species and on humans. MC leads a research program on population genetics, with a focus on the genetic consequences of past demographic events and natural selection. UR is co-leading a long-term project on the evolution and behavioral ecology of nocturnal lemurs in northwestern Madagascar which is conducted by members of the Institute of Zoology, University of Veterinary Medicine Hannover.
We would like to thank CAFF/CORE, the "Direction Générale des Eaux et Forêts" and ANGAP for giving us permission to conduct the fieldwork in Madagascar. For their long lasting collaboration, we thank the late B. Rakotosamimanana, D. Rakotondravony and the late O. Ramilijaona. We are grateful to G. Olivieri, M. Craul, K. Guschanski, K. Marquart, R. Rakotondravony, B. Randrianambinina, and S. Rasoloharijaona for their help in collecting the samples. We thank H. Breitrück for her assistance in the lab and R. Brüning for technical help. We thank E. Zimmermann, B. Crouau-Roy, A. Coutinho, and A. Sanchez-Mazas for constant support. Many thanks go to E.S. Poloni and to two anonymous reviewers for critical comments on the manuscript. The study was made possible through the funding provided by the German Research Council (DFG Ra 502/7 and Ra 502/9 to U.R.), the PROCOPE programme (D/0707573 to U.R., # 17831 to L.C.), the Institut Français de la Biodiversité, Programme Biodiversité de l'Océan Indien (# CD-AOOI-07-003 to L.C.), the FCT Project PTDC/BIA-BDE/71299/2006 to L.C., and the Swiss National Foundation (SNF) grant #3100A0-112651 to A. Sanchez-Mazas. The European Science Foundation is also thanked for funding the ConGen Workshop "Population genetics modelling and habitat fragmentation: separating recent and ancient events for efficient conservation" (ref. 1469), during which the collaboration on the spatial simulations started.
- Hewitt GM: Post-glacial re-colonization of European biota. Biol J Linn Soc. 1999, 68: 87-112. 10.1111/j.1095-8312.1999.tb01160.x.View ArticleGoogle Scholar
- Taberlet P, Fumagalli L, Wust-Saucy AG, Cosson JF: Comparative phylogeography and postglacial colonization routes in Europe. Mol Ecol. 1998, 7: 453-464. 10.1046/j.1365-294x.1998.00289.x.View ArticlePubMedGoogle Scholar
- Burney DA: Madagascar's prehistoric ecosystem. The Natural History of Madagascar. Edited by: Goodman SM, Benstead JP. 2003, Chicago: University of Chicago Press, 47-51.Google Scholar
- Goossens B, Chikhi L, Ancrenaz M, Lackman-Ancrenaz I, Andau P, Bruford MW: Genetic signature of anthropogenic population collapse in orang-utans. PLOS Biology. 2006, 4: 285-291. 10.1371/journal.pbio.0040025.View ArticleGoogle Scholar
- Poux C, Madsen O, Marquard E, Vieites DR, de Jong WW, Vences M: Asynchronous colonization of Madagascar by the four endemic clades of primates, tenrecs, carnivores, and rodents as inferred from nuclear genes. Syst Biol. 2005, 54: 719-730. 10.1080/10635150500234534.View ArticlePubMedGoogle Scholar
- Green GM, Sussman RW: Deforestation history of the eastern rain forests of Madagascar from satellite images. Science. 1990, 248: 212-215. 10.1126/science.248.4952.212.View ArticlePubMedGoogle Scholar
- Harper GJ, Steininger MK, Tucker CJ, Juhn D, Hawkins F: Fifty years of deforestation and forest fragmentation in Madagascar. Environ Cons. 2007, 34: 1-9. 10.1017/S0376892907003554.View ArticleGoogle Scholar
- Smith AP: Deforestation, fragmentation, and reserve design in western Madagascar. Tropical Forest Remnants: Ecology, Management, and Conservation of fragmented Communities. Edited by: Laurance WF, Bierregaard RO Jr. 1997, Chicago: University of Chicago Press, 415-441.Google Scholar
- Bond WJ, Silander JA, Ranivonasy J, Ratsirarson J: The antiquity of Madagascar's grasslands and the rise of C4 grassy biomes. J Biogeogr. 2008, 35: 1743-1758. 10.1111/j.1365-2699.2008.01923.x.View ArticleGoogle Scholar
- Nicoll ME: Forests outside protected areas. The Natural History of Madagascar. Edited by: Goodman SM, Benstead JP. 2003, Chicago: University of Chicago Press, 1432-1437.Google Scholar
- Quéméré E, Louis EE, Ribéron A, Chikhi L, Crouau-Roy B: Non-invasive conservation genetics of the critically endangered golden-crowned sifaka (Propithecus tattersalli): high diversity and significant genetic differentiation over a small range. Conserv Genet. DOI 10.1007/s10592-009-9837-9, Google Scholar
- Radespiel U, Rakotondravony R, Chikhi L: Natural and anthropogenic determinants of genetic structure in the largest remaining metapopulation of the endangered golden-brown mouse lemur, Microcebus ravelobensis. Am J Primatol. 2008, 70: 1-11.Google Scholar
- Frankham R: Conservation genetics. Ann Rev Genet. 1995, 29: 305-527. 10.1146/annurev.ge.29.120195.001513.View ArticlePubMedGoogle Scholar
- Mittermeier RA, Ganzhorn J, Konstant W, Glander K, Tattersall I, Groves CP, Rylands AB, Hapke A, Ratsimbazafy J, Mayor MI, Louis EE, Rumpler Y, Schwitzer C, Rasoloarison RM: Lemur diversity in Madagascar. Int J Prim. 2008, 29: 1607-1656. 10.1007/s10764-008-9317-y.View ArticleGoogle Scholar
- Louis EE, Engberg SE, McGuire SM, McCormick MJ, Randriamampionona R, Ranaivoarisoa JF, Bailey CA, Mittermeier RA, Lei R: Revision of the mouse lemurs, Microcebus (Primates, Lemuriformes), of northern and northwestern Madagascar with descriptions of two new species at Montagne d'Ambre National Park and Antafondro Classified Forest. Primate Conservation. 2008, 23: 19-38. 10.1896/052.023.0103.View ArticleGoogle Scholar
- Radespiel U: Ecological diversity and seasonal adaptations of mouse lemurs (Microcebus spp.). Lemurs: Ecology and Adaptation. Edited by: Gould L, Sauther ML. 2006, New York: Springer, 211-233.View ArticleGoogle Scholar
- Martin RD: Prosimians: from obscurity to extinction?. Creatures of the Dark: The Nocturnal Prosimians. Edited by: Alterman L, Doyle GA, Izard MK. 1995, New York: Plenum Press, 535-563.View ArticleGoogle Scholar
- Wilmé L, Goodman SM, Ganzhorn JU: Biogeographic evolution of Madagascar's microendemic biota. Science. 2006, 312: 1063-1065. 10.1126/science.1122806.View ArticlePubMedGoogle Scholar
- Olivieri G, Zimmermann E, Randrianambinina B, Rasoloharijaona S, Rakotondravony D, Guschanski K, Radespiel U: The ever-increasing diversity in mouse lemurs: Three new species in north and northwestern Madagascar. Mol Phyl Evol. 2007, 4: 309-327. 10.1016/j.ympev.2006.10.026.View ArticleGoogle Scholar
- Radespiel U, Olivieri G, Rasolofoson DW, Rakotondratsimba G, Rakotonirainy O, Rasoloharijaona S, Randrianambinina B, Ratsimbazafy JH, Ratelolahy F, Randriamboavonjy T, Rasolofoharivelo T, Craul M, Rakotozafy L, Randrianarison RM: Exceptional diversity of mouse lemurs (Microcebus spp.) in the Makira region with the description of one new species. Am J Primatol. 2008, 70: 1033-1046. 10.1002/ajp.20592.View ArticlePubMedGoogle Scholar
- Kappeler PM, Rasoloarison RM, Razafimanantsoa L, Walter L, Roos C: Morphology, behaviour and molecular evolution of giant mouse lemurs (Mirza spp.) Gray, with description of a new species. Primate Report. 1870, 71: 3-26.Google Scholar
- Rakotondravony R, Radespiel U: Varying patterns of coexistence of two mouse lemur species (Microcebus ravelobensis and M. murinus) in a heterogeneous landscape. Am J Primatol. 2009, 71: 928-938. 10.1002/ajp.20732.View ArticlePubMedGoogle Scholar
- Yoder AD, Rasoloarison RM, Goodman SM, Irwin JA, Atsalis S, Ravosa MJ, Ganzhorn JU: Remarkable species diversity in Malagasy mouse lemurs (Primates, Microcebus). PNAS. 2000, 97: 11325-11330. 10.1073/pnas.200121897.PubMed CentralView ArticlePubMedGoogle Scholar
- Guschanski K, Olivieri G, Funk SM, Radespiel U: MtDNA reveals strong genetic differentiation among geographically isolated populations of the golden brown mouse lemur, Microcebus ravelobensis. Conserv Genet. 2007, 8: 809-821. 10.1007/s10592-006-9228-4.View ArticleGoogle Scholar
- Olivieri GL, Sousa V, Chikhi L, Radespiel U: From genetic diversity and structure to conservation: genetic signature of recent population declines in three mouse lemur species (Microcebus spp.). Biol Conserv. 2008, 141: 1257-1271. 10.1016/j.biocon.2008.02.025.View ArticleGoogle Scholar
- Wakeley J: Nonequilibrium migration in human history. Genetics. 1999, 153: 1863-1871.PubMed CentralPubMedGoogle Scholar
- Wakeley J: The effects of subdivision on the genetic divergence of populations and species. Evolution. 2000, 54: 1092-1101.View ArticlePubMedGoogle Scholar
- Olivieri G: Molecular phylogenetics and conservation genetics of mouse lemurs (Microcebus spp.) in northwestern Madagascar. PhD thesis. 2007, University of Veterinary Medicine Hannover, Institute of ZoologyGoogle Scholar
- Ray N, Currat M, Excoffier L: Intra-deme molecular diversity in spatially expanding populations. Mol Biol Evol. 2003, 20: 76-86. 10.1093/molbev/msg009.View ArticlePubMedGoogle Scholar
- Bos DH, Gopurenko D, Williams RN, Dewoody JA: Inferring population history and demography using microsatellites, mitochondrial DNA, and major histocompatibility complex (MHC) genes. Evolution. 2008, 62: 1458-1468. 10.1111/j.1558-5646.2008.00364.x.View ArticlePubMedGoogle Scholar
- Städler T, Haubold B, Merino C, Stephan W, Pfaffelhuber P: The impact of sampling schemes on the site frequency spectrum in nonequilibrium subdivided populations. Genetics. 2009, 182: 205-216. 10.1534/genetics.108.094904.PubMed CentralView ArticlePubMedGoogle Scholar
- Chikhi L, Bruford MW: Mammalian population genetics and genomics. Mammalian Genomics. Edited by: Ruvinsky A, Marshall Graves J. 2005, UK: CABI Publishers, 539-583. full_text.View ArticleGoogle Scholar
- Fredsted T, Schierup MH, Groeneveld LH, Kappeler PM: Genetic structure, lack of sex-biased dispersal and behavioral flexibility in the pair-living fat-tailed dwarf lemur, Cheirogaleus medius. Behav Ecol Sociobiol. 2007, 61: 943-954. 10.1007/s00265-006-0323-9.View ArticleGoogle Scholar
- Pennington RT, Dick CW: Diversification of the Amazonian flora and its relation to key geological and environmental events: a molecular perspective. Neogene history of Western Amazonia and its significance for modern biodiversity. Edited by: Hoorn C, Vonhof H, Wesselingh F. Oxford: Blackwell PublishingGoogle Scholar
- Petit RJ, Aguinagalde I, de Beaulieu JL, Bittkau C, Brewer S, Cheddadi R, Ennos R, Fineschi S, Grivet D, Lascoux M, Mohanty A, Muller-Starck G, Demesure-Musch B, Palme A, Martin JP, Rendell S, Vendramin GG: Glacial refugia: hotspots but not melting pots of genetic diversity. Science. 2003, 300: 1563-1565. 10.1126/science.1083264.View ArticlePubMedGoogle Scholar
- Yoder AD, Yang ZH: Divergence dates for Malagasy lemurs estimated from multiple gene loci: geological and evolutionary context. Mol Ecol. 2004, 13: 757-773. 10.1046/j.1365-294X.2004.02106.x.View ArticlePubMedGoogle Scholar
- Poux C, Madsen O, Glos J, de Jong WW, Vences M: Molecular phylogeny and divergence times of Malagasy tenrecs: Influence of data partitioning and taxon sampling on dating analyses. BMC Evol Biol. 2008, 8: 102-10.1186/1471-2148-8-102.PubMed CentralView ArticlePubMedGoogle Scholar
- Townsend TM, Vieites DR, Glaw F, Vences M: Testing species-level diversification hypotheses in Madagascar: The case of microendemic Brookesia Leaf Chamaeleons. Syst Biol. 2009, 58: 641-656. 10.1093/sysbio/syp073.View ArticlePubMedGoogle Scholar
- Maniatis T, Fritsch EF, Sambrook J: Molecular cloning: a laboratory manual. 1982, New York: Cold Spring Harbor Laboratory PressGoogle Scholar
- Gerloff U, Hartung B, Fruth B, Hohmann G, Tautz D: Intracommunity relationships, dispersal pattern and paternity success in a wild living community of Bonobos (Pan paniscus) determined from DNA analysis of faecal samples. Proc R Soc B. 1999, 266: 1189-1195. 10.1098/rspb.1999.0762.PubMed CentralView ArticlePubMedGoogle Scholar
- Nei M: Molecular Evolutionary Genetics. 1987, New York: Columbia University PressGoogle Scholar
- Petit RJ, El Mousadik A, Pons O: Identifying populations for conservation on the basis of genetic markers. Conserv Biol. 1998, 12: 844-855. 10.1046/j.1523-1739.1998.96489.x.View ArticleGoogle Scholar
- Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989, 123: 585-595.PubMed CentralPubMedGoogle Scholar
- Fu YX: Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997, 147: 915-925.PubMed CentralPubMedGoogle Scholar
- Excoffier L, Laval G, Schneider S: Arlequin ver. 3.0: An integrated software package for population genetics data analysis. Evolutionary Bioinformatics. 2005, 1: 47-50.Google Scholar
- Rogers AR, Harpending H: Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol. 1992, 9: 552-569.PubMedGoogle Scholar
- Currat M, Ray N, Excoffier L: SPLATCHE: A program to simulate genetic diversity taking into account environmental heterogeneity. Mol Ecol Notes. 2004, 4: 139-142. 10.1046/j.1471-8286.2003.00582.x.View ArticleGoogle Scholar
- Currat M, Excoffier L: The effect of the Neolithic expansion on European molecular diversity. Proc R Soc B. 2005, 272: 679-688. 10.1098/rspb.2004.2999.PubMed CentralView ArticlePubMedGoogle Scholar
- Currat M: Effets des expansions des populations humaines en Europe sur leur diversité génétique. PhD thesis. 2004, University of Geneva, Department of AnthropologyGoogle Scholar
- Tamura K, Nei M: Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol. 1993, 10: 512-526.PubMedGoogle Scholar
- Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.View ArticlePubMedGoogle Scholar
- Mantel N: The detection of disease clustering and a generalized regression approach. Cancer Res. 1967, 27: 209-220.PubMedGoogle Scholar
- Ihaka R, Gentleman R: R: A language for data analysis and graphics. Journal of Computational and Graphical Statistics. 1996, 5: 299-314. 10.2307/1390807.Google Scholar
- Clement M, Posada D, Crandall KA: TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000, 9: 1657-1659. 10.1046/j.1365-294x.2000.01020.x.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.