Hybridization of mouse lemurs: different patterns under different ecological conditions
© Hapke et al; licensee BioMed Central Ltd. 2011
Received: 24 May 2011
Accepted: 11 October 2011
Published: 11 October 2011
Skip to main content
© Hapke et al; licensee BioMed Central Ltd. 2011
Received: 24 May 2011
Accepted: 11 October 2011
Published: 11 October 2011
Several mechanistic models aim to explain the diversification of the multitude of endemic species on Madagascar. The island's biogeographic history probably offered numerous opportunities for secondary contact and subsequent hybridization. Existing diversification models do not consider a possible role of these processes. One key question for a better understanding of their potential importance is how they are influenced by different environmental settings. Here, we characterized a contact zone between two species of mouse lemurs, Microcebus griseorufus and M. murinus, in dry spiny bush and mesic gallery forest that border each other sharply without intermediate habitats between them. We performed population genetic analyses based on mtDNA sequences and nine nuclear microsatellites and compared the results to a known hybrid zone of the same species in a nearby wide gradient from dry spiny bush over transitional forest to humid littoral forest.
In the spiny-gallery system, Microcebus griseorufus is restricted to the spiny bush; Microcebus murinus occurs in gallery forest and locally invades the dryer habitat of its congener. We found evidence for bidirectional introgressive hybridization, which is closely linked to increased spatial overlap within the spiny bush. Within 159 individuals, we observed 18 hybrids with mitochondrial haplotypes of both species. Analyses of simulated microsatellite data indicate that we identified hybrids with great accuracy and that we probably underestimated their true number. We discuss short-term climatic fluctuations as potential trigger for the dynamic of invasion and subsequent hybridization. In the gradient hybrid zone in turn, long-term aridification could have favored unidirectional nuclear introgression from Microcebus griseorufus into M. murinus in transitional forest.
Madagascar's southeastern transitional zone harbors two very different hybrid zones of mouse lemurs in different environmental settings. This sheds light on the multitude of opportunities for the formation of hybrid zones and indicates an important influence of environmental factors on secondary contact and hybridization. Our findings suggest that hybridization could enhance the adaptability of mouse lemurs without necessarily leading to a loss of distinctiveness. They point to a potential role of hybridization in Madagascar's diversification history that requires further investigation.
There is increasing evidence that natural hybridization between animal species occurs more frequently than previously appreciated [1, 2]. Moreover, the use of genetic techniques has facilitated the detection of many cases of introgressive hybridization between animals . Hybridization can constitute a threat for biodiversity when rare endemic species come into contact with widespread invaders and are hybridized out of existence [4, 5], but there are also examples of hybridizing species that maintain distinctiveness in the face of interspecific gene flow (e.g. [6, 7]). Introgressive hybridization can even allow for the transfer of beneficial adaptations between species, facilitate rapid adaptation to changing environmental conditions and thus play an important role for diversification [2, 3, 8]. According to different models, environmental factors can influence hybrid zones to various degrees. Tension zones are maintained by a balance of immigration of the parental species and endogenous selection against hybrids, which is independent of environmental factors . In other kinds of hybrid zones, environmental selection in mosaics or gradients of different habitats influences hybridization [1, 10, 11]. According to the tension zone and mosaic zone models, hybrids are generally less fit than the parental species [9, 11]. According to the bounded hybrid superiority model , they are more fit than the parental species in ecotonal habitats. Under the evolutionary novelty model , endogenous and environmental selection act in concert, and certain hybrid genotypes can be as fit as or more fit than the parental species in ecotonal and parental habitats. Temporal change of environmental conditions can also influence hybridization. Prominent examples for hybrid superiority due to climatic fluctuations are Darwin's finches on Daphne Major Island .
Madagascar has an extremely high level of endemism both at the species level and at higher taxonomic levels and is among the world's eight hottest biodiversity hotspots . Due to its long isolation from other landmasses, the island harbors many endemic radiations that gave rise to numerous microendemic species with very restricted ranges [14–16]. In various groups of organisms, there are also macroendemic species with larger ranges that overlap with those of their microendemic congeners. Examples are leaf chameleons [17, 18], cophyline frogs [17, 19], tufted-tailed rats  and mouse lemurs . Patterns of microendemism in Madagascar are just beginning to emerge in the course of a recent and ongoing wave of species detections [14, 22]. In recent years, a variety of models have been proposed that aim to explain the evolution of this diversity of microendemics [15, 19, 23, 24]. Most of these models focus on allopatric or vicariant speciation. Consequently, a potential role of hybridization for the diversification of Malagasy endemics has so far widely been neglected although the present biogeography offers numerous opportunities for the formation of hybrid zones along environmental gradients ( but see ). There is a controversial discussion about a potential role of Pleistocene climatic changes for the evolution of microendemics in Madagascar [15, 18, 24]. At least it is clear that Madagascar's current biogeography arose after dramatic changes in the course of climatic fluctuations, which are best documented for the Holocene . Most endemic organisms must have undergone range shifts, retractions and expansions in the course of climatic fluctuations . The study of microevolutionary processes at species boundaries is thus important for a better understanding of how species adapted to changing environmental conditions and how the diversity of microendemics was maintained throughout the vicissitudes of the Pleistocene. It appears probable that there were even more opportunities for the formation of hybrid zones in the dynamic past of Madagascar's biogeography than is evident from current patterns. At present, there is a belt of humid forest in eastern Madagascar, seasonally dry deciduous forest in western Madagascar and dry spiny bush in the South. Dry western and humid eastern formations are largely separated from each other by mostly forestless central highlands. Dry-humid ecotones are restricted to the North and the South, but there is evidence for the existence of forest corridors that connected western dry and eastern humid forest through the central highlands during the Holocene [23, 26]. Heckman et al.  discuss secondary contact through such corridors and hybridization as potential alternative explanation to incomplete lineage sorting for phylogeographic patterns in western and eastern mouse lemurs.
There is an ongoing discussion about the question how many species of mouse lemurs should be distinguished (e.g. [28, 29]). Weisrock et al.  delineated 16 population-level lineages based on sequence data from two mitochondrial and four nuclear loci from a large sample of mouse lemurs comprising localities all across Madagascar. According to the authors , it depends on the species concept applied how many of the 16 lineages deserve species rank. The currently recognized species Microcebus murinus comprises several of these lineages. Weisrock et al.  did not formally subdivide Microcebus murinus into different species. In the study presented here, we use the name "Microcebus murinus" in the sense of "M. murinus sensu latu" in Weisrock et al. .
Potential contact zones between species adapted to dry and mesic or humid conditions are not limited to such wide ecotonal gradients in Madagascar. For example, gallery forests with mesic conditions exist throughout the southern dry spiny bush, thus providing a virtually endless line of contact between different kinds of habitat.
It is unclear, which role hybridization could have played in Madagascar's diversification history. One key question towards a better understanding of its possible role is how far hybridization is influenced by ecological conditions. The aim of our study presented here was to contribute to the answering of the latter question by comparing contact zones of the same model species in different ecological settings. We questioned how far hybridization between Microcebus griseorufus and M. murinus would be restricted to the presence of an intermediate habitat where hybrids might have selective advantages. Would the same species also hybridize where they come into contact under different ecological conditions and would such a hybrid zone have a different structure? In order to answer these questions, we performed a population genetic study of the two species within the western rain shadow of the Anosy Mountains where spiny bush and gallery forest border each other sharply without intermediate habitats between them. We refer to this system as the spiny-gallery contact zone and to the hybrid zone in transitional forest as the gradient contact zone in the following (Figures 1 and 2).
Our genetic data yield evidence for hybridization in the spiny-gallery contact zone. Hybridization patterns in the two zones differ in several aspects. We discuss these differences and possible environment-related processes that could influence hybridization.
We trapped mouse lemurs with banana-baited Sherman traps set up approximately 1 m above ground level. After anesthetizing the animals, we took small tissue samples from the ear for genetic analyses. We stored the samples in 90% ethanol at ambient temperature. We released all animals in the late afternoon of the same day at the respective sites of capture. We performed the trapping of mouse lemurs and the collection of tissue samples in compliance with respective research authorizations by the Malagasy Ministère de l'Environnement, des Eaux et Forêts (No 179/06/MINENV.EF/SG/DGEF/DPB/SCBLF/RECH, No 0174/07 -MINENV.EF/SG/DGEEF/DVRN/SPE) and by the Malagasy Ministère de l'Environnement, des Eaux et Forêts et du Tourisme (No 091/08/MEEFT/SG/DGEF/DSAP/SSE).
We isolated DNA from tissue samples using the DNeasy Blood & Tissue Kit (Qiagen) following the protocol for the purification of total DNA from animal tissues. We genotyped all individuals using the same genetic markers as Gligor et al. , which comprised the mitochondrial hypervariable region 1 (HV1) and nine nuclear microsatellite loci.
We amplified the HV1 with one primer binding in the constant region of the mitochondrial D-loop (mih1coau: 5'-GTTATAGTTTCAGGTTAGTCA-3') and one of the following primers binding in the cytochrome b gene: mih1cbau (5'-GATCTACTTATCCTTACATGA-3'), Mcytbf (5'-CTAGTAGAATGRATCTGAGG-3'), MrFTDcytbinfw58 (5'-GATTCTTCGCATTCCACTTC-3') or TsimMgCytbfw2 (5'-TCGGACAAGTGGCCTCTAT-3'). Typical PCR conditions comprised an initial denaturation step of 2 min at 92°C, 35 to 40 cycles of 40 s denaturation at 92°C, 60 s annealing at 55°C and 60 s elongation at 72°C, and one final elongation step of 5 min at 72°C. When using primer mih1cbau, we changed the annealing temperature to 54°C. When using primers Mcytbf or MrFTDcytbinfw58, we changed the elongation time to 70 s. We performed wax-mediated hot-start-PCR using the Qiagen Core Kit. For sequencing of the HV1 on both strands, we used the BigDye version 3.1 kit (Applied Biosystems) on a 3130 Genetic Analyzer (Applied Biosystems). We used mih1coau combined with either mih1cbau, TsimMgCytbfw2 or mih1cbin2Mz (5'-TTATACCWACYGTAAGYCTT-3') as primers for sequencing.
We applied the microsatellite markers 33104, Mm21, Mm22, Mm39, Mm51, Mm30, Mm42, Mm43b and Mm60  with the following modifications: We performed a first wax-mediated hot-start PCR with non-labeled primers and then reamplified 1-4 μl of the resulting product in repeated unidirectional extensions of one fluorescently labeled primer both with the Qiagen Core Kit. We changed the annealing temperature for locus Mm42 to 54°C. We used a 3130 Genetic Analyzer (Applied Biosystems) for electrophoresis and the GENEMAPPER 3.0 software (Applied Biosystems) for raw data analysis. We used 6-FAM and HEX as fluorescent dyes and GeneScan™-350 ROX™ (Applied Biosystems) as length standard.
We aligned the sequences of the HV1 using the CLUSTAL W module implemented in BIOEDIT version 7.0.1  and corrected the resulting alignment visually. We collapsed identical sequences to haplotypes using FABOX . We performed phylogenetic tree reconstructions in order to assign haplotypes to the two species of mouse lemurs. To this aim, we added reference sequences of Microcebus griseorufus ([GenBank: EU109652], ) and M. murinus ([GenBank: DQ865143], ) to the dataset. For tree reconstructions based on maximum parsimony and Bayesian inference, we used a sequence of Microcebus ravelobensis ([GenBank: AF285455], ) as outgroup.
We used PAUP* version 4.0b10  for a maximum parsimony tree reconstruction. We performed a heuristic search with 100 random addition replicates and TBR branch swapping and a bootstrap analysis with 100 replicates for the evaluation of relative levels of support for internal nodes.
For the following tree reconstructions, we selected most appropriate substitution models for our data based on the Akaike Information Criterion as implemented in JMODELTEST version 0.1.1 [41, 42]. We performed likelihood calculations for 24 models, which included 3 substitution schemes, equal or unequal base frequencies, a proportion of invariable sites and rate variation among sites with 4 rate categories on maximum likelihood optimized trees.
We used the MRBAYES software version 3.1.2 [43, 44] for a tree reconstruction via Bayesian inference. We applied the HKY+G model in 2 independent analyses with 4 Markov Chain Monte Carlo chains and sampled the resulting trees every 100th generation. We stopped the analysis after 2,000,000 generations, when the standard deviation of split frequencies was 0.006532. We discarded the first 5,000 samples as burnin and summed the parameters from the remaining 15,000 samples. There was no increasing or decreasing trend in the log probabilities over generations. The potential scale reduction factor was 1.000 for all parameters. Based on these observations, we concluded that the analysis had converged. We discarded the first 5,000 trees as burnin and computed a consensus tree and posterior probabilities for internal nodes from the remaining 15,000 trees.
We used PHYML version 3.0  for a maximum likelihood tree reconstruction, where we did not include an outgroup. For this dataset, JMODELTEST selected the HKY+G substitution model with a transition-transversion-parameter kappa of 80.0366 and a gamma shape parameter alpha of 0.0740. We applied the HKY+G model with 4 rate categories and fixed kappa and alpha to the values estimated by JMODELTEST. We performed tree searches using the SPR method with 5 random starting trees and the simultaneous NNI method with a BioNJ starting tree and selected the best tree overall with the aid of the BEST method. We evaluated support of internal branches by a bootstrap analysis with 100 replicates.
We compared the degree of spatial overlap between individuals with mitochondrial haplotypes of Microcebus griseorufus and M. murinus at Mangatsiaka and Tsimelahy as follows. For each individual with a griseorufus-like haplotype, we determined the 6 nearest neighbors based on pairwise distances between individual trapping positions rounded to the nearest 10 m (Additional file 1: Individual trap positions at Mangatsiaka and Tsimelahy). We then calculated the average proportion of murinus-like haplotypes in the nearest neighborhood of griseorufus-like individuals for each sampling site. In order to evaluate a potential impact of unequal numbers of individuals with the two species' haplotypes, we calculated averages from 100 datasets, where we reduced the greater group to the same size as the smaller group by random resampling. We wrote the program SOA  for these analyses.
We further tested for significant differences of the proportions of murinus-like and griseorufus-like mitochondrial haplotypes among individuals captured at Mangatsiaka and Tsimelahy with the aid of Fisher's exact test as implemented in PASW Statistics 17.0.
For the following tests, we divided the individuals into local samples according to sampling sites and griseorufus-like (-Mg) or murinus-like (-Mm) mitochondrial haplotypes. We excluded some resulting small samples with sample sizes between 1 and 7 and included the following ones: Hzf-Mg (n = 16), Abt-Mm (n = 19), Mtk-Mg (n = 12), Mtk-Mm (n = 63), Tml-Mg (n = 17) and Tml-Mm (n = 19). We used FSTAT version 22.214.171.124  for the following analyses. We estimated global F-statistics and FST between pairs of samples according to Weir and Cockerham . We tested for deviations from Hardy-Weinberg equilibrium for each locus and over all loci in each sample with 54,000 permutations of the original data. We tested for linkage disequilibria for all pairs of loci in each sample with the aid of LINKDOS  as implemented in GENETIX version 4.05.2 .
We used two Bayesian methods implemented in the programs STRUCTURE version 2.1  and NEWHYBRIDS version 1.1  to identify hybrids based on individual microsatellite genotypes in comparison to mitochondrial haplotypes. We used three datasets of microsatellite genotypes: One large dataset comprised all individuals in the study; two smaller ones comprised all individuals at Mangatsiaka and all individuals at Tsimelahy.
With NEWHYBRIDS, we estimated posterior probabilities to belong to one of six predefined categories of purebreds and hybrids for each individual's genotype. We used the default genotype frequency class file with the following categories: Mg (purebred Microcebus griseorufus), Mm (purebred M. murinus), F1 (Mg × Mm), F2 (F1 × F1) and two classes of first generation backcrosses Mg-Bx1 (F1 × Mg) and Mm-Bx1 (F1 × Mm). We performed Markov chain Monte Carlo computations with a burnin period of 100,000 steps and a post-burnin period of 1,000,000 steps with Jeffreys-like priors for the mixing proportions and allele frequencies and without prior population or allele frequency information. As recommended by Vähä and Primmer , we used a threshold of 0.5 as criterion for the detection of hybrids and the distinction between hybrid categories. We regarded an individual as purebred when the posterior probability to be purebred from the species corresponding to its mitochondrial haplotype was > 0.5, as a hybrid when it was ≤0.5 and as a specific category of hybrid when the corresponding posterior probability was > 0.5.
STRUCTURE performs a Bayesian clustering of genotypes into a number of clusters K predefined by the user and estimates membership coefficients as posterior probabilities for each genotype to belong to each of the K clusters. We used the admixture ancestry model with independent allele frequencies and without prior population information. We fixed the allele frequency parameter lambda to one and let the program estimate a uniform value for the model parameter alpha. We performed 10 independent runs with a burnin of 10,000 generations and a post-burnin period of 40,000 generations for different values of K from 1 to 9 for the large dataset and from 1 to 5 for the two smaller datasets. We then used the ad hoc statistic ΔK  to determine the most appropriate number of clusters for our data. ΔK unambiguously indicated K = 2 as most appropriate for all datasets. We then performed optimal alignments of the results from the 10 independent runs with K = 2 for each dataset and calculated average membership coefficients using the full-search algorithm of CLUMPP version 1.1.1 . We applied a threshold of 0.1 to the membership-coefficients as criterion for hybrid detection. We regarded an individual as purebred when the average membership coefficient for the cluster corresponding to its mitochondrial haplotype was > 0.9 and as a hybrid when it was ≤ 0.9.
In order to evaluate the power of our microsatellite data to detect hybrids and to distinguish between different hybrid classes, we analyzed simulated datasets with STRUCTURE and NEWHYBRIDS. Within our original microsatellite genotypes, we identified all individuals observed as purebred from the species corresponding to their mitochondrial haplotypes with probabilities > = 0.9 both with STRUCTURE and NEWHYBRIDS. We then used these individuals' genotypes to create simulated datasets with the software HYBRIDLAB version 1.0 . We analyzed three sets of simulated data:
Simulation A: We simulated 100 datasets, which each contained genotypes of 100 Mg (purebred Microcebus griseorufus), 100 Mm (purebred M. murinus), 10 F1 (Mg × Mm) and 10 F2 (F1 × F1). Each file further contained genotypes that we call first generation backcrosses, 10 Mg-Bx1 (F1 × Mg) and 10 Mm-Bx1 (F1 × Mm) and second generation backcrosses, 10 Mg-Bx2 (Mg-Bx1 × Mg) and 10 Mm-Bx2 (Mm-Bx1 × Mm). We analyzed all simulated datasets with STRUCTURE and NEWHYBRIDS under the same settings as our real data. With STRUCTURE, we performed 10 independent runs with K = 2 for each dataset and summed up the results with CLUMPP. We evaluated different threshold values of 0.3, 0.2 and 0.1 for hybrid detection. With NEWHYBRIDS, we used a threshold value of 0.5 as with the real data. We then calculated efficiency and accuracy similar to Vähä and Primmer . Efficiency is the proportion of individuals of a certain category correctly identified (e.g. true F1-hybrids identified as hybrids), and accuracy is the proportion of individuals assigned to a certain category that was correctly identified (e.g. individuals identified as purebred M. griseorufus that are truly purebred of that species).
Simulation B: At Mangatsiaka, we had sampled 12 individuals with griseorufus-like and 63 with murinus-like mitochondrial haplotypes. We questioned if the number of hybrids that we observed in the real data at Mangatsiaka could be an artifact due to these unequal sample sizes. We simulated a scenario where all individuals at Mangatsiaka are purebred with the aid of 100 datasets, which comprised each 12 purebred Microcebus griseorufus and 63 purebred M. murinus. We analyzed these datasets under the same settings as the real data with both programs and calculated the proportions of datasets where the number of false-positive hybrids was as great as or greater than the number of hybrids detected in the real data.
Simulation C: In the real data, we had identified several individuals as hybrids with only one of both programs or in only one of two datasets. We questioned if such discrepancies more probably indicated false-positive or false-negative hybrids. We used real genotypes of purebred Microcebus griseorufus and M. murinus from Tsimelahy to simulate 100 local datasets with the same size and composition as in simulation A. We then enlarged each dataset by adding 100 purebred individuals of each species, which we simulated based on purebred genotypes from Hazofotsy and Ambatoabo. We analyzed all local and enlarged datasets with both programs under the same settings as in simulations A and B and evaluated the proportions of true hybrids among classifications that were discrepant between programs or datasets. We then combined the evidence from both programs and corresponding datasets by accepting all individuals as hybrids that we identified as such at least once and calculated efficiencies and accuracies.
Pairwise FST between samples
Tests of heterozygote deficiency and linkage disequilibria
FIS All Loci
Simulation A, STRUCTURE, efficiency and accuracy with different threshold values
Simulation A, NEWHYBRIDS, efficiency and accuracy
Mg (n = 10,000)
Mm (n = 10,000)
F1 (n = 1000)
F2 (n = 1000)
Mg-Bx1 (n = 1000)
Mm-Bx1 (n = 1000)
Mg-Bx2 (n = 1000)
Mm-Bx2 (n = 1000)
The simulated scenario of no hybridization at Mangatsiaka appears as highly improbable with respect to the results of simulation B (Additional file 6: Simulation B: Scenario of no hybrids at Mangatsiaka). Most simulated purebred individuals were correctly identified with both programs. The number of purebred Microcebus murinus misclassified as hybrids was smaller than the number of hybrids with murinus-like haplotypes detected in the real data with both programs in all simulated datasets. The respective number for Microcebus griseorufus was reached in only 4% of the simulated datasets with NEWHYBRIDS.
Simulation C, discrepancies between datasets and programs
Discrepant between sets
Disc. true Hyb
Discrepant between programs
Disc. true Hyb
N Hyb obs. per set/pair sets
Our results demonstrate that Microcebus griseorufus and M. murinus hybridize within the spiny-gallery contact zone. In the following, we discuss our findings from this area and compare them to those from the gradient contact zone. We refer to Gligor et al.  for a more detailed characterization of the latter.
The environmental settings are very different in both contact zones. In the spiny-gallery zone, a network of narrow gallery forests penetrates the dry spiny bush. Both habitats border each other sharply without intermediate zones between them. In the gradient contact zone, there is a large climatic gradient with an approximately 10 km wide intermediate zone with transitional forests between dry spiny bush west of it and humid littoral forests to the east.
According to the distributions of mitochondrial haplotypes in our study area, Microcebus griseorufus appears as strictly bound to the spiny bush and M. murinus as closely associated with the gallery forest. The fact that the latter also invades adjacent habitats of different types leads to local sympatric admixture within the spiny bush at Mangatsiaka. This is in contrast to observations at other places. At Beza Mahafaly in southwestern Madagascar, Heckman et al.  observed only griseorufus-like mitochondrial haplotypes despite morphological evidence for the presence of both species [21, 56]. At Berenty, west of our study area, Yoder et al.  observed mitochondrial haplotypes of Microcebus murinus in gallery forests and of M. griseorufus in adjacent spiny bush without admixture. In the gradient contact zone investigated by Gligor et al. , griseorufus-like haplotypes where restricted to spiny bush and murinus-like haplotypes to transitional and littoral forest.
Other than in the gradient contact zone, we observed hybrids with mitochondrial haplotypes of both species in the spiny-gallery contact zone. The number of detected hybrids indicates a limited degree of hybridization although we might have underestimated it. F1-hybrids are apparently rare. It appears thus as probable that F2-hybrids are rare as well and that most hybrids are backcrosses or crosses between backcrosses. This would indicate that the formation of F1-hybrids is more difficult than the interbreeding of hybrids with purebreds or with other hybrids, which would provide opportunities for introgression as soon as a few hybrids exist [2, 57]. Our results point to a limited degree of bidirectional introgressive hybridization, which would also be in line with the observation of heterozygote deficiencies and linkage disequilibria in subsamples of both species separated according to their mitochondrial haplotypes. The directionality of introgressive hybridization is not to be confounded with the directionality of habitat-invasion. The latter is apparently unidirectional because Microcebus murinus invades the spiny bush, whereas we have no evidence for an invasion of the gallery forest by M. griseorufus.
Gligor et al.  did not apply a threshold value for hybrid detection when using STRUCTURE. In their study, ΔK  indicated K = 2 as the most probable uppermost structure but there was also some support for K = 3. We re-evaluated the original results of Gligor et al.  with K = 2 and applied a threshold of 0.1 as in our study presented here. Among 38 individuals from the transition zone, we identified 28 hybrids in a large dataset and 35 hybrids in a smaller dataset (Additional file 7: Re-evaluated identification of hybrids in the gradient contact zone). All individuals in the spiny bush and littoral forest were classified as purebred. With K = 3, most individuals from the transition zone formed a third cluster of genotypes, while those from the spiny bush and littoral forest were assigned to the other two clusters . None of the murinus-like mitochondrial haplotypes in the transition zone occurred in the littoral forests further east. Based on this observation and their results from microsatellite data, Gligor et al.  concluded on a unidirectional introgression of griseorufus-like nuclear alleles into autochthonous populations of Microcebus murinus in the transition zone.
Despite bidirectionality, hybridization in the spiny-gallery contact zone is apparently asymmetric since the majority of hybrids carry murinus-like mitochondrial haplotypes, which is similar to the gradient contact zone, where all hybrids carry murinus-like haplotypes. In the spiny-gallery zone, secondary contact occurs most probably within the habitat of Microcebus griseorufus, which is invaded by M. murinus. In the gradient contact zone, the respective mechanism is unknown since Gligor et al.  did not detect recent purebred immigrants in either direction. They discussed cyto-nuclear incompatibility [58, 59] versus environmental selection, potentially on cyto-nuclear gene complexes (e.g. ), as possible causes for asymmetric introgression. Our findings in the spiny-gallery contact zone now contradict strict incompatibility of either species' mitochondrial genome with a hybridized nuclear background. Locally unequal abundances of both species in the spiny-gallery contact zone could play a role for asymmetric hybridization. This is most obvious at Mangatsiaka, where we observed most hybrids and where Microcebus murinus appears to be more abundant than its congener.
Microcebus griseorufus appears as a specialist restricted to the southern spiny bush and M. murinus as a generalist that colonizes a wider range of habitats within an extremely large range from northern to southeastern Madagascar. In general, widespread species tend to be locally more abundant than species with restricted ranges. This has been explained by metapopulation dynamics  and by Brown's  niche-breadth hypothesis, which predicts that among closely related species, those with the broader niche will be both locally more abundant and more widespread. It remains, however, questionable if greater ecological plasticity alone can explain the apparent success of Microcebus murinus at Mangatsiaka. According to Ganzhorn and Schmid , the species is extremely sensitive to altered environmental conditions. They compared populations in primary and secondary dry deciduous forest in western Madagascar. In secondary forest, Microcebus murinus displayed lower population densities and extremely small year-to-year survival rates. The authors explained this by higher ambient temperatures and lower availability of tree holes, which decreased the possibilities for energy-saving daily torpor and hibernation.
The significant linkage disequilibria and heterozygote deficiencies in the spiny-gallery zone could most likely point to selection against or in favor of certain allelic combinations in the course of introgressive hybridization. Theoretically, they could also be produced by recent intraspecific admixture, e.g. in the course of colonization-recolonization dynamics, but, with respect to very small pairwise FST values between samples of the same mitochondrial species, such a scenario appears not very plausible. As well, in the gradient-contact zone, significant heterozygote deficiencies and linkage disequilibria, which were mostly restricted to the transition zone, pointed to selection . Environmental selection appears as probable in both zones because genetic patterns are congruent with habitat patterns. The distribution of purebreds and hybrids in the spiny-gallery zone closely reflects the mosaic-like distribution of spiny bush and gallery forest. In the gradient zone, three clusters of genotypes observed with STRUCTURE with K = 3 were largely congruent with the three vegetation zones along the west-eastern climatic gradient . With K = 2, membership coefficients were predominantly griseorufus-like in the western part of the transition zone and predominantly murinus-like in the eastern part (Additional file 7: Re-evaluated identification of hybrids in the gradient contact zone). Genetic distances between sites were significantly correlated with vegetation zone differences under control of geographic distances .
A remaining open question is why we did observe heterozygote deficiencies and linkage disequilibria at several sampling sites in the spiny-gallery zone, where we identified no or few hybrids. One could explain this by the difficult detection of hybrids, which might be more pronounced in smaller local samples. Another potential explanation could be temporal change of the frequency of hybridization, which could be due to stochastically fluctuating environmental conditions. Indeed, there is some evidence for strong environmental stochasticity that could have an important impact on hybridization in the spiny-gallery contact zone. According to Dewar and Richard , Madagascar has significantly less predictable rainfall than continental Africa. In the North and the South, this unpredictability takes the form of high interannual variation of total precipitation. The weather station at Behara, which is situated 18 km west of Mangatsiaka in our study area, displayed a mean annual rainfall of 532 mm and the least predictable climate among 15 weather stations throughout Madagascar included in the study of Dewar and Richard . The climate in our study area is characterized by a short rainy season and a long dry season (e.g. ). The region is irregularly struck by severe droughts, which can last one or several subsequent years [66–68]. According to Elmqvist et al. , the frequency of droughts in southern Madagascar is increasing since the 1970's (but see ). Without long-term data from our study area, we can only speculate how these dramatic fluctuations of rainfall might influence hybridization. Neaves et al.  discuss density fluctuations of sympatric grey kangaroos due to fluctuations in rainfall as potential cause of occasional hybridization. A hypothetical scenario for our study area could be that density fluctuations due to irregular rainfall are most pronounced for Microcebus murinus, which is adapted to more mesic conditions. Such fluctuations could even follow different rhythms or have different amplitudes near gallery forests with headwaters in the spiny bush or in the rainforest. Occasional high densities would lead to increased overlap with Microcebus griseorufus and subsequent hybridization. Intermittent droughts could reduce the density of Microcebus murinus, reduce invasive pressure and shift the adaptive landscape towards more divergent selection, which could help to maintain species integrity.
Pollen data from southern Madagascar  and subfossil remains at Andrahomana cave [71–73] indicate that the vegetation zones in the gradient contact zone could have shifted in eastward direction in the course of long-term aridification during the last 3000 years approximately. Accordingly, Gligor et al.  propose that such a shift of the adaptive landscape would have provided a growing advantage for introgressing griseorufus-like alleles. The apparently strong degree of hybridization within the transition zone could thus result from long-term past introgression. Without long-term climate data from this area, it remains unclear if selection in the course of ongoing aridification could be responsible for heterozygote deficiencies and linkage disequilibria within the transition zone. It could even be possible that short-term climatic fluctuations superimpose a long-term trend of aridification and contribute to ongoing change of environmental selective pressures.
The environment-related process proposed by Gligor et al.  implies an adaptive value of hybridization in the gradient zone, which facilitates adaptation to changing environmental conditions. In the spiny-gallery zone, so far, we have regarded the apparently great relative abundance of Microcebus murinus at Mangatsiaka as a prerequisite that possibly facilitated hybridization. In fact, we were astonished to see that the great majority of individuals at this site carried murinus-like haplotypes because it appeared as suboptimal habitat for the species. The headwaters of the temporary watercourses at Mangatsiaka are situated within dry spiny bush to the north and northeast, while those of all other gallery forests in our study are situated in the eastern rainforest. Accordingly, there is only very sparse gallery vegetation at Mangatsiaka. It appears now as possible that Microcebus murinus locally adapted to dryer conditions at Mangatsiaka and that hybridization could have facilitated this adaptation.
There is increasing evidence for potentially beneficial consequences of hybridization such as facilitation of adaptability and diversification [2, 3, 8]. Hybridization can enhance the invasibility of intruding species . For Microcebus griseorufus in turn, one could expect rather negative consequences of hybridization. Hybridization with a more common, widespread invader is potentially deleterious and often seen as a threat for rare specialized endemics [4, 5]. On the other hand, hybridizing species can remain distinct when hybrid zones are narrow  or when genotypes display bimodal distributions [6, 75]. Even rare endemic specialists can maintain distinctiveness despite hybridization with a more common invader .
Microcebus griseorufus and M. murinus maintain their distinctiveness in the two contact zones in different ways. In the gradient contact zone, introgressed populations in transitional forest cannot be assigned to either parental species, but the two species remain distinct in adjacent spiny bush and littoral forest. In the spiny-gallery contact zone, the strongly bimodal microsatellite genotypes and their concordance with mitochondrial haplotypes indicate that selection is apparently strong enough to preserve the distinctiveness of both species in the face of local sympatric admixture and bidirectional introgressive hybridization. This does not preclude the potential acquisition of single beneficial adaptations. The identification of respective candidate loci, however, would require a genome wide scan with numerous loci.
Microcebus griseorufus and M. murinus are sister species, but they are not among the most closely related species of mouse lemurs detected so far (e.g. [27, 29, 39]). The fact that they hybridize generates the expectation that hybrid zones between further species of mouse lemurs could exist or could have existed in the past. Our study demonstrates that the two species display very different patterns of hybridization under different ecological conditions. This highlights the importance of environmental factors for the formation of different kinds of contact zones. Our results show that the transition zone of southeastern Madagascar harbors different hybrid zones in very different ecological settings within a small geographic area. This exemplifies that a multitude of opportunities for the formation of different contact and hybrid zones might exist or have existed within the complex biogeography of Madagascar that still need to be explored. Our study points to interesting perspectives on the potential role of hybridization in the evolution of Madagascar's endemics that require further investigation: It appears as possible that macroendemic species acquire beneficial adaptations through hybridization with microendemic congeners that allow them to extend their ranges into novel kinds of habitat. At the same time, the two hybrid zones of mouse lemurs exemplify that species can maintain distinctiveness despite introgressive hybridization.
Madagascar's biogeography has been repeatedly reshuffled during the Pleistocene and Holocene climatic fluctuations. These processes appear as rather recent and rapid when compared to phylogenetic timescales. Consequently, their importance for the diversification of endemics in Madagascar is controversial [15, 16, 18, 24]. It might be possible to reconcile these conflicting views partially when we pay more attention to the rapid evolutionary process of hybridization.
The fieldwork for this study was carried out under the "Accord de Collaboration" between Madagascar National Parks (MNP), the University of Antananarivo and the University of Hamburg. We thank the Ministère de l'Environnement, des Eaux et Forêts et du Tourisme, MNP and the University of Antananarivo for the authorization and support of this study. We are indebted to Chantal Andrianarivo, Jocelyn Rakotomalala, Domoina Rakotomalala of MNP, the late Olga Ramilijaona, Daniel Rakotondravony and Jean-Baptiste Ramanamanjato for their collaboration and support. We gratefully acknowledge Ernest Refaly, Kai Schütte, Ole Theisinger, Susanne Kobbe, Roberto Rakoto, Fiadana Ijo, Tahiana Andriaharimalala, Miha Mahatanty, Mosa Vonjoloke, Soja Fahamaro, Maharavo, Mandiavato and Bertin Tsimihanta for their help in the field and Irene Tomaschewski for help in the lab. We are indebted to Jörg Ganzhorn who conceived of the field study and supervised it and who provided helpful comments on the manuscript. We thank Hans Zischler for important critical discussions about the analyses of genetic data. We thank two anonymous reviewers for constructive criticism on an earlier version of this manuscript.
The DFG (GA 342/14-1, HA 4762/1-1, HA 4762/2-1) and DAAD provided financial support for this study.
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.