Delimiting species without nuclear monophyly in Madagascar's mouse lemurs

The genomic diversity of the Eliurus genus in northern Madagascar with a putative new species

Madagascar exhibits extraordinarily high level of species richness and endemism, while being severely threatened by habitat loss and fragmentation (HL&F). In front of these threats to biodiversity, conservation effort can be directed, for instance, in the documentation of species that are still unknown to science, or in investigating how species respond to HL&F. The tufted-tail rats genus (Eliurus spp.) is the most speciose genus of endemic rodents in Madagascar, with 13 described species, which occupy two major habitat types: dry or humid forests. The large species diversity and association to specific habitat types make the Eliurus genus a suitable model for investigating species adaptation to new environments, as well as response to HL&F (dry vs humid). In the present study, we investigated Eliurus spp. genomic diversity across northern Madagascar, a region covered by both dry and humid fragmented forests. From the mitochondrial DNA (mtDNA) and nuclear genomic (RAD-seq) data of 124 Eliurus individuals sampled in poorly studied forests of northern Madagascar, we identified an undescribed Eliurus taxon (Eliurus sp. nova). We tested the hypothesis of a new Eliurus species using several approaches: i) DNA barcoding; ii) phylogenetic inferences; iii) species delimitation tests based on the Multi-Species Coalescent (MSC) model, iv) genealogical divergence index (gdi); v) an ad-hoc test of isolation-by-distance within versus between sister-taxa, vi) comparisons of %GC content patterns and vii) morphological analyses. All analyses support the recognition of the undescribed lineage as a putative distinct species. In addition, we show that Eliurus myoxinus, a species known from the dry forests of western Madagascar, is, surprisingly, found mostly in humid forests in northern Madagascar. In conclusion, we discuss the implications of such findings in the context of Eliurus species evolution and diversification, and use the distribution of northern Eliurus species as a proxy for reconstructing past changes in forest cover and vegetation type in northern Madagascar.

Molecular phylogenetics of the sucking louse genus Lemurpediculus (Insecta: Phthiraptera), ectoparasites of lemurs, with descriptions of three new species

Sucking lice live in intimate association with their hosts and often display a high degree of host specificity. The present study investigated sucking lice of the genus Lemurpediculus from six mouse lemur (Microcebus) and two dwarf lemur (Cheirogaleus) species endemic to the island of Madagascar, considered a biodiversity hotspot. Louse phylogenetic trees were created based on cytochrome C oxidase subunit I (COI), elongation factor 1α (EF1α) and internal transcribed spacer 1 (ITS1) sequences. While clustering according to host species was generally observed for COI and ITS1, suggesting high host specificity of the examined lice, EF1α sequences alone did not distinguish between lice of different Microcebus species, possibly due to rather recent divergence. As bootstrap support for basal tree structure was rather low, further data are necessary to resolve the evolutionary history of louse-mouse lemur associations. Three new species of sucking lice are described: Lemurpediculus zimmermanni sp. Nov. From Microcebus ravelobensis, Lemurpediculus gerpi sp.nov. from Microcebus gerpi, and Lemurpediculus tsimanampesotsae sp. nov. from Microcebus griseorufus. These new species are compared with all known congeneric species and identifying features are illustrated for all known species of Lemurpediculus.

Long-term host-pathogen evolution of endogenous beta- and gammaretroviruses in mouse lemurs with little evidence of recent retroviral introgression

Madagascar's flora and fauna have evolved in relative isolation since the island split from the African and Indian continents. When the last common ancestors of lemurs left Africa between 40 and 70 million years ago, they carried a subset of the viral diversity of the mainland population within them, which continued to evolve throughout the lemur radiation. Relative to other primate radiations, we know very little about the past or present viral diversity of lemurs, particularly mouse lemurs. Using high-throughput sequencing, we identified two gammaretroviruses and three betaretroviruses in the genomes of four species of wild mouse lemurs. The two gammaretroviruses and two betaretroviruses have not previously been described. One betaretrovirus was previously identified. All identified viruses are present in both Lorisiformes and Lemuriformes but absent from haplorrhine primates. The estimated ages of these viruses are consistent with the estimated divergence dates of the host lineages, suggesting they colonized the lemur genome after the Haplorrhine-Strepsirrhine split, but before the Lorisiformes-Lemuriformes split and before the colonization of Madagascar. The viral phylogenies connect multiple lineages of retroviruses from non-lemur and non-Madagascar-native species, suggesting substantial cross-species transmission occurred deep in the primate clade prior to its geographic dispersal. These phylogenies provide novel insights into known retroviral clades. They suggest that the origin of gammaretroviruses in rodents or bats may be premature and that the Jaagsiekte sheep virus clade may be older and more broadly distributed among mammals than previously thought.

Skull morphological evolution in Malagasy endemic Nesomyinae rodents

Madagascar is a large island to the south-east of Africa and in many ways continental in size and ecological complexity. Here we aim to define how skull morphology of an endemic and monophyletic clade of rodents (sub-family Nesomyinae), that show considerable morphological variation, have evolved and how their disparity is characterized in context of the geographical and ecological complexity of the island. We performed a two-dimensional geometric morphometric analysis on 370 dorsal and 399 ventral skull images of 19 species (comprising all nine extant endemic genera) and tested the influence of three ecological parameters (climate, locomotor habitat and nychthemeral cycle) in a phylogenetic context on size and shape. The results indicate that skull shape appears to importantly reflect phylogeny, whereas skull size does not carry a significant phylogenetic signal. Skull shape is significantly influenced by climate while, skull size is not impacted by any of the ecological factors tested, which is controversial to expectations in an insular context. In conclusion, Nesomyinae must have evolved under unusual types of local constraints, preventing this radiation from demonstrating strong ecological release.

Past environmental changes affected lemur population dynamics prior to human impact in Madagascar

Quaternary climatic changes have been invoked as important drivers of species diversification worldwide. However, the impact of such changes on vegetation and animal population dynamics in tropical regions remains debated. To overcome this uncertainty, we integrated high-resolution paleoenvironmental reconstructions from a sedimentary record covering the past 25,000 years with demographic inferences of a forest-dwelling primate species (Microcebus arnholdi), in northern Madagascar. Result comparisons suggest that climate changes through the African Humid Period (15.2 - 5.5 kyr) strongly affected the demographic dynamics of M. arnholdi. We further inferred a population decline in the last millennium which was likely shaped by the combination of climatic and anthropogenic impacts. Our findings demonstrate that population fluctuations in Malagasy wildlife were substantial prior to a significant human impact. This provides a critical knowledge of climatically driven, environmental and ecological changes in the past, which is essential to better understand the dynamics and resilience of current biodiversity.

SPART: A versatile and standardized data exchange format for species partition information

A wide range of data types can be used to delimit species and various computer-based tools dedicated to this task are now available. Although these formalized approaches have significantly contributed to increase the objectivity of species delimitation (SD) under different assumptions, they are not routinely used by alpha-taxonomists. One obvious shortcoming is the lack of interoperability among the various independently developed SD programs. Given the frequent incongruences between species partitions inferred by different SD approaches, researchers applying these methods often seek to compare these alternative species partitions to evaluate the robustness of the species boundaries. This procedure is excessively time consuming at present, and the lack of a standard format for species partitions is a major obstacle. Here, we propose a standardized format, SPART, to enable compatibility between different SD tools exporting or importing partitions. This format reports the partitions and describes, for each of them, the assignment of individuals to the "inferred species". The syntax also allows support values to be optionally reported, as well as original trees and the full command lines used in the respective SD analyses. Two variants of this format are proposed, overall using the same terminology but presenting the data either optimized for human readability (matricial SPART) or in a format in which each partition forms a separate block (SPART.XML). ABGD, DELINEATE, GMYC, PTP and TR2 have already been adapted to output SPART files and a new version of LIMES has been developed to import, export, merge and split them.

Mouse Lemurs in an Assemblage of Cheirogaleid Primates in Menabe Central, Western Madagascar – Three Reasons to Coexist

Ecological communities are structured by interactions between coexisting species that mutually influence their distribution and abundance. Ecologically similar species are expected to exclude one another from suitable habitat, so the coexistence of two mouse lemur species in an assemblage of several closely related cheirogaleid primates in the central Menabe region of Madagascar requires explanation. We assessed the occurrence of Madame Berthe’s mouse lemurs (Microcebus berthae) and Gray mouse lemurs (Microcebus murinus), and of two larger cheirogaleids, Coquerel’s giant mouse lemur (Mirza coquereli) and the western fat-tailed dwarf lemur (Cheirogaleus medius), by nocturnal line transect walks between 2003 and 2007. We explored interspecific interactions for four different scenarios with varying resource availability (degraded and non-degraded habitat in the wet and dry season), both on the regional spatial scale and on a finer local (transect) scale. We tested whether the interspecific distribution of mouse lemur individuals indicates interspecific competition and whether their regional coexistence might be stabilized by interactions with M. coquereli or C. medius. We developed the “Inter-Species Index of Attraction” (ISIA) to quantify the observed interspecific interactions within transects and determined if these were significantly different from a null model generated by a combination of randomization and bootstrapping to control for intraspecific aggregation. For the two mouse lemurs, interspecific spatial exclusion was most pronounced during the resource-poor dry season, consistent with the hypothesis of feeding competition. Seasonally varying distribution patterns indicated resource tracking in a spatio-temporally heterogeneous environment. The interspecific distribution of individuals suggested that the larger cheirogaleids benefit M. berthae at the expense of the more abundant M. murinus: spatial associations of both, M. coquereli and C. medius, with M. murinus were negative in most scenarios and across spatial scales, but neutral or even positive with M. berthae. Thus, our study revealed that coexistence among ecologically similar heterospecifics can rely on complex density-mediated interspecific processes varying with habitat quality and season. With regard to the stability of animal assemblages, this insight has major implications for biodiversity conservation.

Ecology and morphology of mouse lemurs (Microcebus spp.) in a hotspot of microendemism in northeastern Madagascar, with the description of a new species

Delimitation of cryptic species is increasingly based on genetic analyses but the integration of distributional, morphological, behavioral, and ecological data offers unique complementary insights into species diversification. We surveyed communities of nocturnal mouse lemurs (Microcebus spp.) in five different sites of northeastern Madagascar, measuring a variety of morphological parameters and assessing reproductive states for 123 individuals belonging to five different lineages. We documented two different non-sister lineages occurring in sympatry in two areas. In both cases, sympatric species pairs consisted of a locally restricted (M. macarthurii or M. sp. #3) and a more widespread lineage (M. mittermeieri or M. lehilahytsara). Estimated Extents of Occurrence (EOO) of these lineages differed remarkably with 560 and 1,500 km² versus 9,250 and 50,700 km² , respectively. Morphometric analyses distinguished unambiguously between sympatric species and detected more subtle but significant differences among sister lineages. Tail length and body size were most informative in this regard. Reproductive schedules were highly variable among lineages, most likely impacted by phylogenetic relatedness and environmental variables. While sympatric species pairs differed in their reproductive timing (M. sp. #3/M. lehilahytsara and M. macarthurii/M. mittermeieri), warmer lowland rainforests were associated with a less seasonal reproductive schedule for M. mittermeieri and M. lehilahytsara compared with populations occurring in montane forests. Distributional, morphological, and ecological data gathered in this study support the results of genomic species delimitation analyses conducted in a companion study, which identified one lineage, M. sp. #3, as meriting formal description as a new species. Consequently, a formal species description is included. Worryingly, our data also show that geographically restricted populations of M. sp. #3 and its sister species (M. macarthurii) are at high risk of local and perhaps permanent extinction from both deforestation and habitat fragmentation.

Cryptic Patterns of Speciation in Cryptic Primates: Microendemic Mouse Lemurs and the Multispecies Coalescent

Mouse lemurs (Microcebus) are a radiation of morphologically cryptic primates distributed throughout Madagascar for which the number of recognized species has exploded in the past two decades. This taxonomic revision has prompted understandable concern that there has been substantial oversplitting in the mouse lemur clade. Here, we investigate mouse lemur diversity in a region in northeastern Madagascar with high levels of microendemism and predicted habitat loss. We analyzed RADseq data with multispecies coalescent (MSC) species delimitation methods for two pairs of sister lineages that include three named species and an undescribed lineage previously identified to have divergent mtDNA. Marked differences in effective population sizes, levels of gene flow, patterns of isolation-by-distance, and species delimitation results were found among the two pairs of lineages. Whereas all tests support the recognition of the presently undescribed lineage as a separate species, the species-level distinction of two previously described species, M. mittermeieri and M. lehilahytsara is not supported-a result that is particularly striking when using the genealogical discordance index (gdi). Nonsister lineages occur sympatrically in two of the localities sampled for this study, despite an estimated divergence time of less than 1 Ma. This suggests rapid evolution of reproductive isolation in the focal lineages and in the mouse lemur clade generally. The divergence time estimates reported here are based on the MSC calibrated with pedigree-based mutation rates and are considerably more recent than previously published fossil-calibrated relaxed-clock estimates. We discuss the possible explanations for this discrepancy, noting that there are theoretical justifications for preferring the MSC estimates in this case. [Cryptic species; effective population size; microendemism; multispecies coalescent; speciation; species delimitation.].

Cryptic diversity and species boundaries within the Paragalago zanzibaricus species complex

The recently described genus Paragalago is a complex of several nocturnal and morphologically cryptic species distributed in the forests of eastern Africa. Species diversity within this genus has been mainly described using species-specific differences in their loud calls. However, molecular data are still lacking for this group and species boundaries remain unclear. In this study, we explore species diversity within the zanzibaricus-complex using a combination of mitochondrial and nuclear data and comparing multiple species delimitation methods. Our results consistently support the existence of three independent lineages, P. cocos, P. zanzibaricus, and P. granti, confirming previous hypotheses based on vocal data. We conclude that these three lineages represent valid cryptic species and we hypothesize that speciation within this complex was characterized by cycles of forest expansion and contraction in the Plio-Pleistocene.

Evolutionary significance of the variation in acoustic communication of a cryptic nocturnal primate radiation (Microcebus spp.)

Acoustic phenotypic variation is of major importance for speciation and the evolution of species diversity. Whereas selective and stochastic forces shaping the acoustic divergence of signaling systems are well studied in insects, frogs, and birds, knowledge on the processes driving acoustic phenotypic evolution in mammals is limited. We quantified the acoustic variation of a call type exchanged during agonistic encounters across eight distinct species of the smallest-bodied nocturnal primate radiation, the Malagasy mouse lemurs. The species live in two different habitats (dry forest vs. humid forest), differ in geographic distance to each other, and belong to four distinct phylogenetic clades within the genus. Genetically defined species were discriminated reliably on the phenotypic level based on their acoustic distinctiveness in a discriminant function analysis. Acoustic variation was explained by genetic distance, whereas differences in morphology, forest type, or geographic distance had no effect. The strong impact of genetics was supported by a correlation between acoustic and genetic distance and the high agreement in branching pattern between the acoustic and molecular phylogenetic trees. In sum, stochastic factors such as genetic drift best explained acoustic diversification in a social communication call of mouse lemurs.

Quantitative morphological characters of the skull suggest that Akodon oenos (Rodentia, Cricetidae, Sigmodontinae) is not a junior synonym of A. spegazzinii

Akodon oenos, a nominal form with type locality in northern Mendoza province, Argentina, has an intricate taxonomic history. Recently, it was considered a junior synonym of Akodon spegazzinii on the basis of morphological similarities of one individual captured in southern Mendoza with specimens belonging to the type series of A. oenos, plus the genetic similarity of that single specimen with representatives of A. spegazzinii. We compared specimens of the Akodon boliviensis species group recently captured in Mendoza with the type series of A. oenos and series of A. spegazzinii from northwestern Argentina. We conducted an analysis of molecular markers (cytochrome-b gene) and quantitative morphological studies (comparisons of skull measurements). Our results showed that: (1) all sequenced specimens of the A. boliviensis group from Mendoza form a monophyletic clade closely related to A. spegazzinii; (2) contrary to previous hypothesis, specimens of Akodon from Mendoza are morphologically distinctive in size and shape of the skull, from A. spegazzinii; and (3) within the A. boliviensis species group, specimens from Mendoza are distinctive. We argue that the synonymy of A. oenos under A. spegazzinii must await further studies, including molecular markers other than mitochondrial and samples of DNA sequences from topotypical specimens of A. oenos.

Genetic and morphological diversity of mouse lemurs (Microcebus spp.) in northern Madagascar: The discovery of a putative new species?

Tropical forests harbor extremely high levels of biological diversity and are quickly disappearing. Despite the increasingly recognized high rate of habitat loss, it is expected that new species will be discovered as more effort is put to document tropical biodiversity. Exploring under-studied regions is particularly urgent if we consider the rapid changes in habitat due to anthropogenic activities. Madagascar is known for its extraordinary biological diversity and endemicity. It is also threatened by habitat loss and fragmentation. It holds more than 100 endemic primate species (lemurs). Among these, Microcebus (mouse lemurs) is one of the more diverse genera. We sampled mouse lemurs from several sites across northern Madagascar, including forests never sampled before. We obtained morphological data from 99 Microcebus individuals; we extracted DNA from tissue samples of 42 individuals and amplified two mitochondrial loci (cytb and cox2) commonly used for species identification. Our findings update the distribution of three species (Microcebus tavaratra, Microcebus arnholdi, and Microcebus mamiratra), including a major increase in the distribution area of M. arnholdi. We also report the discovery of a new Microcebus lineage genetically related to M. arnholdi. Several complementary approaches suggest that the newly identified Microcebus lineage might correspond to a new putative species, to be confirmed or rejected with additional data. In addition, morphological analyses showed (a) clear phenotypic differences between M. tavaratra and M. arnholdi, but no clear differences between the new Microcebus lineage and the sister species M. arnholdi; and (b) a significant correlation between climatic variables and morphology, suggesting a possible relationship between species identity, morphology, and environment. By integrating morphological, climatic, genetic, and spatial data of two northern Microcebus species, we show that the spatial distribution of forest-dwelling species may be used as a proxy to reconstruct the past spatial changes in forest cover and vegetation type.

An empirical estimate of the generation time of mouse lemurs

The generation time of organisms drives the rate of change in populations and across evolutionary times. In long-lived species, generation time should also account for overlapping generations, and the average age of parents has been proposed as a best approximation under these conditions. This study uses this definition to estimate the generation time of a widely studied small primate, Microcebus murinus, based on parentage data generated for a free-living population over a 6-year period in northwestern Madagascar. The average age of parents was calculated separately for mothers and fathers of three different offspring cohorts that differed in the degree of demographic uncertainty. In addition, adult survival rates were calculated for males and females based on long-term capture data from the same population to estimate the possible upper limits of generation time. Adult survival was low with only 44% of adult females and 38% of adult males being recaptured at the beginning of their second breeding season. The average age of mothers was 1.56-1.91 years, pointing toward a 2-year female generation time due to the high proportion of 1-year old mothers in all three cohorts. Female generation time estimates were fairly stable across the three offspring cohorts. In contrast, the average age of fathers differed by more than 1 year from the first to the third offspring cohort (1.71-2.83 years) pointing toward a 3-year generation time, but also suggesting a higher degree of demographic uncertainty in the early years of the study. For future modeling purposes, we, therefore, propose to use the average, 2.5 years, of male and female values as new estimate for the generation time of mouse lemurs.

Unraveling patterns and processes of diversification in the South Andean-Patagonian Nassauvia subgenus Strongyloma (Asteraceae, Nassauvieae)

Congruence among different sources of data is highly desirable in phylogenetic analyses. However, plastid and nuclear DNA may record different evolutionary processes such that incongruence among results from these sources can help unravel complex evolutionary histories. That is the case of Nassauvia subgenus Strongyloma (Asteraceae), a taxon with five putative species distributed in the southern Andes and Patagonian steppe. Morphometric and phylogeographic information cast doubt on the integrity of its species, and previous molecular data even questioned the monophyly of the subgenus. We tested those questions using plastid and nuclear DNA sequences by the application of different methods such as phylogenetic trees, networks, a test of genealogical sorting, an analysis of population structure, calibration of the trees, and hybridization test, assembling non-synchronous incongruent results at subgenus and species levels in a single reconstruction. The integration of our molecular analyses and previous taxonomic, morphological, and molecular studies support subgenus Strongyloma as a monophyletic group. However, the topology of the nuclear trees and the evidence of polyploids within subgenus Nassauvia, suggest a hypothetical origin and initial radiation of Nassauvia related to an ancient hybridization event that occurred around 17-6.3 Myr ago near the Andes in west-central Patagonia. Plastid data suggest a recent diversification within subgenus Strongyloma, at most 9.8 Myr ago, towards the Patagonian steppe east of the Andes. These processes cause phylogenies to deviate from the species tree since each putative species lack exclusive ancestry. The non-monophyly of its species using both plastid and nuclear data is caused mainly by incomplete lineage sorting occurred since the Miocene. The final uplift of the Andes and Pliocene-Pleistocene glacial-interglacial and its consequences on the landscape and climate structured the genetic composition of this group of plants in the Patagonian steppe. The molecular data presented here agree with previous morphological studies, in that the five putative species typically accepted in this subgenus are not independent taxa. This study emphasizes that adding more than one sequence per species, not combining data with dissimilar inheritance patterns without first performed incongruence tests, exploring data through different methodologies, considering the timing of events, and searching for the causes of poorly resolved and/or incongruent phylogenies help to reveal complex biological underlying processes, which might otherwise remain hidden.

Challenges of next-generation sequencing in conservation management: Insights from long-term monitoring of corridor effects on the genetic diversity of mouse lemurs in a fragmented landscape

Long-term genetic monitoring of populations is essential for efforts aimed at preserving genetic diversity of endangered species. Here, we employ a framework of long-term genetic monitoring to evaluate the effects of fragmentation and the effectiveness of the establishment of corridors in restoring population connectivity and genetic diversity of mouse lemurs Microcebus ganzhorni. To this end, we supplement estimates of neutral genetic diversity with the assessment of adaptive genetic variability of the major histocompatibility complex (MHC). In addition, we address the challenges of long-term genetic monitoring of functional diversity by comparing the genotyping performance and estimates of MHC variability generated by single-stranded conformation polymorphism (SSCP)/Sanger sequencing with those obtained by high-throughput sequencing (next-generation sequencing [NGS], Illumina), an issue that is particularly relevant when previous work serves as a baseline for planning management strategies that aim to ensure the viability of a population. We report that SSCP greatly underestimates individual diversity and that discrepancies in estimates of MHC diversity attributable to the comparisons of traditional and NGS genotyping techniques can influence the conclusions drawn from conservation management scenarios. Evidence of migration among fragments in Mandena suggests that mouse lemurs are robust to the process of fragmentation and that the effect of corridors is masked by ongoing gene flow. Nonetheless, results based on a larger number of shared private alleles at neutral loci between fragment pairs found after the establishment of corridors in Mandena suggest that gene flow is augmented as a result of enhanced connectivity. Our data point out that despite low effective population size, M. ganzhorni maintains high individual heterozygosity at neutral loci and at MHC II DRB gene and that selection plays a predominant role in maintaining MHC diversity. These findings highlight the importance of long-term genetic monitoring in order to disentangle between the processes of drift and selection maintaining adaptive genetic diversity in small populations.

Diversification of the widespread neotropical frog Physalaemus cuvieri in response to Neogene-Quaternary geological events and climate dynamics

Here we reconstructed the demographical history and the dispersal dynamics of Physalemus cuvieri through the Neogene-Quaternary periods by coupling DNA regions with different mutation rates, ecological niche modelling, reconstruction of spatio-temporal lineage dispersal and coalescent simulations. Still, to test alternative diversification scenarios we used approximate Bayesian computation. Molecular phylogenetic analysis recovered four deep and strongly supported clades, which we interpret as population lineages. The ancestral location reconstruction placed the root in southcentral Amazonia, and the dispersal events indicate that spatial displacement was widespread early in the diversification of this species. The demographical scenario of "Multiple Refugia" with recent lineage admixture was the most likely hypothesis to predict the observed genetic parameters of P. cuvieri. Our results revealed that Neogene orogenic events might have played a prominent role in the early diversification of P. cuvieri. The species shows deep divergences with strong regional population structure, despite its widespread distribution. Final uplift of the central Brazilian Plateau and formation of the river basins in Central South America played an important role in the origin, diversification and the maintenance of P. cuvieri lineages.

Evidence supporting cryptic species within two sessile microinvertebrates, Limnias melicerta and L. ceratophylli (Rotifera, Gnesiotrocha)

Microorganisms, including rotifers, are thought to be capable of long distance dispersal. Therefore, they should show little population genetic structure due to high gene flow. Nevertheless, substantial genetic structure has been reported among populations of many taxa. In rotifers, genetic studies have focused on planktonic taxa leaving sessile groups largely unexplored. Here, we used COI gene and ITS region sequences to study genetic structure and delimit cryptic species in two sessile species (Limnias melicerta [32 populations]; L. ceratophylli [21 populations]). Among populations, ITS region sequences were less variable as compared to those of the COI gene (ITS; L. melicerta: 0-3.1% and L. ceratophylli: 0-4.4%; COI; L. melicerta: 0-22.7% and L. ceratophylli: 0-21.7%). Moreover, L. melicerta and L. ceratophylli were not resolved in phylogenetic analyses based on ITS sequences. Thus, we used COI sequences for species delimitation. Bayesian Species Delimitation detected nine putative cryptic species within L. melicerta and four putative cryptic species for L. ceratophylli. The genetic distance in the COI gene was 0-15.4% within cryptic species of L. melicerta and 0.5-0.6% within cryptic species of L. ceratophylli. Among cryptic species, COI genetic distance ranged 8.1-21.9% for L. melicerta and 15.1-21.2% for L. ceratophylli. The correlation between geographic and genetic distance was weak or lacking; thus geographic isolation cannot be considered a strong driver of genetic variation. In addition, geometric morphometric analyses of trophi did not show significant variation among cryptic species. In this study we used a conservative approach for species delimitation, yet we were able to show that species diversity in these sessile rotifers is underestimated.

Emerging and Enduring Issues in Primate Conservation Genetics

Conservation genetics is a branch of conservation biology that uses molecular data to assist in the conservation and management of imperiled populations, subspecies, and species. In this review, I examine conservation action plans (CAPs)—instrumental documents designed to influence conservation policy—for selected primate species. I use the information contained in CAPs as a means to guide this review. The primary genetics-based topics that are mentioned in CAPs are genetic connectivity, inbreeding, and subspecies/species delimitation. I discuss these topics as well as historical demographic inference and hybridization using examples from wild primate species to illustrate the myriad ways in which genetics can assist in conservation efforts. I also discuss some recent technological advances such as genomic capture techniques and the potential to do molecular work in remote locations.

Variation in reproduction of the smallest-bodied primate radiation, the mouse lemurs (Microcebus spp.): A synopsis

Reproduction is a fundamental trait in the life history of any species and contributes to species diversity and evolution. Here, we aim to review the barely known variation in reproductive patterns of the smallest-bodied primate radiation, the Malagasy mouse lemurs, focusing on twelve species of four phylogenetic clades. We present a new reproductive field dataset collected between May and November 1996-2016 for nine species (Microcebus murinus, M. myoxinus, M. ravelobensis, M. bongolavensis, M. danfossi, M. sambiranensis, M. margothmarshae, M. mamiratra, and M. lehilahytsara) and add published field information on three additional species. In the majority of species, the estrus of females was recorded in the period of long days (day length longer than 12 hr), whereas male testes size increased about one to three months prior to this. Reproductive schedules varied considerably between the four clades. Sympatric species-pairs of different clades differed in the timing of female and male reproduction, suggesting strong phylogenetic constraints. Populations of the same species in a different ecological setting varied in the onset of reproduction, suggesting substantial environmental plasticity. Warm temperatures and rainfall throughout the year may allow for less expressed reproductive seasonality. Our results suggest that an interplay between phylogenetic relatedness, ambient temperature (as a proxy for thermo regulatory constraints), and rainfall (as a proxy for food availability), may best explain this variation. Findings further point to a more complex control of mouse lemur reproduction than previously described and illuminate phylogenetic constraints and adaptive potentials in behavioral reaction norms of a species-rich primate radiation.

Genome-wide single-nucleotide polymorphism data reveal cryptic species within cryptic freshwater snail species-The case of the Ancylus fluviatilis species complex

DNA barcoding utilizes short standardized DNA sequences to identify species and is increasingly used in biodiversity assessments. The technique has unveiled an unforeseeably high number of morphologically cryptic species. However, if speciation has occurred relatively recently and rapidly, the use of single gene markers, and especially the exclusive use of mitochondrial markers, will presumably fail in delimitating species. Therefore, the true number of biological species might be even higher. One mechanism that can result in rapid speciation is hybridization of different species in combination with polyploidization, that is, allopolyploid speciation. In this study, we analyzed the population genetic structure of the polyploid freshwater snail Ancylus fluviatilis, for which allopolyploidization was postulated as a speciation mechanism. DNA barcoding has already revealed four cryptic species within A. fluviatilis (i.e., A. fluviatilis s. str., Ancylus sp. A-C), but early allozyme data even hint at the presence of additional cryptic lineages in Central Europe. We combined COI sequencing with high-resolution genome-wide SNP data (ddRAD data) to analyze the genetic structure of A. fluviatilis populations in a Central German low mountain range (Sauerland). The ddRAD data results indicate the presence of three cryptic species within A. fluviatilis s. str. occurring in sympatry and even syntopy, whereas mitochondrial sequence data only support the existence of one species, with shared haplotypes between species. Our study hence points to the limitations of DNA barcoding when dealing with organismal groups where speciation is assumed to have occurred rapidly, for example, through the process of allopolyploidization. We therefore emphasize that single marker DNA barcoding can underestimate the true species diversity and argue in strong favor of using genome-wide data for species delimitation in such groups.

Shared evolutionary origin of major histocompatibility complex polymorphism in sympatric lemurs

Genes of the major histocompatibility complex (MHC) play a central role in adaptive immune responses of vertebrates. They exhibit remarkable polymorphism, often crossing species boundaries with similar alleles or allelic motifs shared across species. This pattern may reflect parallel parasite-mediated selective pressures, either favouring the long maintenance of ancestral MHC allelic lineages across successive speciation events by balancing selection ("trans-species polymorphism"), or alternatively favouring the independent emergence of functionally similar alleles post-speciation via convergent evolution. Here, we investigate the origins of MHC similarity across several species of dwarf and mouse lemurs (Cheirogaleidae). We examined MHC class II variation in two highly polymorphic loci (DRB, DQB) and evaluated the overlap of gut-parasite communities in four sympatric lemurs. We tested for parasite-MHC associations across species to determine whether similar parasite pressures may select for similar MHC alleles in different species. Next, we integrated our MHC data with those previously obtained from other Cheirogaleidae to investigate the relative contribution of convergent evolution and co-ancestry to shared MHC polymorphism by contrasting patterns of codon usage at functional vs. neutral sites. Our results indicate that parasites shared across species may select for functionally similar MHC alleles, implying that the dynamics of MHC-parasite co-evolution should be envisaged at the community level. We further show that balancing selection maintaining trans-species polymorphism, rather than convergent evolution, is the primary mechanism explaining shared MHC sequence motifs between species that diverged up to 30 million years ago.

Phylogeography of the Rice Spittle Bug (Callitettix versicolor) Implies Two Long-Term Mountain Barriers in South China

South China is a region of remarkable topographic complexity. However, the impact of climate fluctuations in the Pleistocene on the local fauna and especially insects has not been extensively studied. We integrated mitochondrial DNA (mtDNA) and microsatellite data of the rice spittle bug, Callitettix versicolor, to determine the genetic structure, potential biogeographic barriers, and historical demography of this species. The mtDNA data revealed two distinct lineages (Western and Eastern) congruent with the geographically separated western and eastern sub-regions of the Hengduan Mountains. The Eastern lineage was subdivided into two sub-lineages, E1 and E2, congruent with the geographically separated northern and southern sub-regions of the Dabie Mountains. E2 was further subdivided into two sub-groups, E2-1 and E2-2, with a hybrid zone (Guizhou and Hunan Provinces) in which their areas were contiguous. The genetic structures constructed using mtDNA were corroborated by four clusters (G1-G4) of microsatellite data. The populations of each cluster were nearly consistent with a sub-lineage of the mtDNA gene tree (G1-G4 corresponded to the Western, E1, E2-2 and E2-1 lineages, respectively). The divergence time estimated between the Western and Eastern lineages was 1.17 (0.50-2.37) to 0.89 (0.39-1.78) Mya, indicating that the lineages diversified on both geographic and temporal scales. The historical demography of the Eastern lineage showed continuous population growth after the Last Interglacial (LIG) and a stable population during the Last Glacial Maximum (LGM) period. However, the Western lineage remained largely unchanged during the LIG and LGM periods. This suggests that the historical demography of C. versicolor is probably related not only to the paleoclimate of South China, but also to the geological restriction and specific habitat preferences of species.

The Mouse Lemur, a Genetic Model Organism for Primate Biology, Behavior, and Health

Systematic genetic studies of a handful of diverse organisms over the past 50 years have transformed our understanding of biology. However, many aspects of primate biology, behavior, and disease are absent or poorly modeled in any of the current genetic model organisms including mice. We surveyed the animal kingdom to find other animals with advantages similar to mice that might better exemplify primate biology, and identified mouse lemurs (Microcebus spp.) as the outstanding candidate. Mouse lemurs are prosimian primates, roughly half the genetic distance between mice and humans. They are the smallest, fastest developing, and among the most prolific and abundant primates in the world, distributed throughout the island of Madagascar, many in separate breeding populations due to habitat destruction. Their physiology, behavior, and phylogeny have been studied for decades in laboratory colonies in Europe and in field studies in Malagasy rainforests, and a high quality reference genome sequence has recently been completed. To initiate a classical genetic approach, we developed a deep phenotyping protocol and have screened hundreds of laboratory and wild mouse lemurs for interesting phenotypes and begun mapping the underlying mutations, in collaboration with leading mouse lemur biologists. We also seek to establish a mouse lemur gene "knockout" library by sequencing the genomes of thousands of mouse lemurs to identify null alleles in most genes from the large pool of natural genetic variants. As part of this effort, we have begun a citizen science project in which students across Madagascar explore the remarkable biology around their schools, including longitudinal studies of the local mouse lemurs. We hope this work spawns a new model organism and cultivates a deep genetic understanding of primate biology and health. We also hope it establishes a new and ethical method of genetics that bridges biological, behavioral, medical, and conservation disciplines, while providing an example of how hands-on science education can help transform developing countries.

Taxonomic position of Chinese voles of the tribe Arvicolini and the description of 2 new species from Xizang, China

China has 26 species in the tribe Arvicolini. The taxonomic status of these voles remains controversial despite much effort. Herein, we evaluate the taxonomic position of 22 species plus 2 unidentified taxa using mitochondrial DNA gene sequences (cytb + CO1). We also evaluate 18 species and 2 unidentified taxa using morphological data. Phylogenetic analyses of cytb resolve monophyly for the genera Alexandromys, Lasiopodomys, Microtus, Neodon, Proedromys, and Volemys with strong support. Stenocranius clusters with Chionomys but with very weak support. Analyses of concatenated cytb + CO1 resolve the same genera with strong support, but the topology of the tree differs from that of cytb in that Chionomys roots at the base of the tree independent of Stenocranius, which forms the sister-group of Lasiopodomys in a more terminal position. The matrilineal genealogy excludes the type species Arvicola amphibius from the rest of the Arvicolini. This species forms the sister-group of Ondatra with high support. Neodon includes N. irene, N. linzhiensis, N. fuscus, N. leucurus, N. sikimensis, Microtus clarkei, and 2 unidentified specimens. Alexandromys includes the former species Microtus oeconomus, M. kikuchii, M. limnophilus, M. fortis, and M. maximowiczii. Finally, Microtus has the subgenera Blanfordimys, Microtus, Mynomes, Pedomys, Pitymys, and Terricola, which includes the Chinese species M. agrestis, M. arvalis, and Blanfordimys juldaschi. General mixed Yule-coalescent species delimitation modeling demarcates 6 currently recognized species and 2 new species of Neodon. A principal component analysis of the morphological data among 7 matrilines shows that all variables have positive loadings of high magnitude on the 1st component. Canonical discriminant analysis for Neodon (including M. clarkei and 2 unidentified species) correctly classifies 93.0% of specimens. Overall, our analyses support the recognition of Alexandromys, Lasiopodomys, Microtus, Neodon, Proedromys, and Volemys as genera. Stenocranius includes Microtus gregalis, and the genealogical position of Stenocranius remains uncertain. The status of Arvicola requires further study. We assign M. clarkei to Neodon and describe 2 new species of Neodon.

Population genetics of mouse lemur vomeronasal receptors: current versus past selection and demographic inference

BACKGROUND

A major effort is underway to use population genetic approaches to identify loci involved in adaptation. One issue that has so far received limited attention is whether loci that show a phylogenetic signal of positive selection in the past also show evidence of ongoing positive selection at the population level. We address this issue using vomeronasal receptors (VRs), a diverse gene family in mammals involved in intraspecific communication and predator detection. In mouse lemurs, we previously demonstrated that both subfamilies of VRs (V1Rs and V2Rs) show a strong signal of directional selection in interspecific analyses. We predicted that ongoing sexual selection and/or co-evolution with predators may lead to current directional or balancing selection on VRs. Here, we re-sequence 17 VRs and perform a suite of selection and demographic analyses in sympatric populations of two species of mouse lemurs (Microcebus murinus and M. ravelobensis) in northwestern Madagascar.

RESULTS

M. ravelobensis had consistently higher genetic diversity at VRs than M. murinus. In general, we find little evidence for positive selection, with most loci evolving under purifying selection and one locus even showing evidence of functional loss in M. ravelobensis. However, a few loci in M. ravelobensis show potential evidence of positive selection. Using mismatch distributions and expansion models, we infer a more recent colonisation of the habitat by M. murinus than by M. ravelobensis, which most likely speciated in this region earlier on.

CONCLUSIONS

These findings suggest that the analysis of VR variation is useful in inferring demographic and phylogeographic history of mouse lemurs. In conclusion, this study reveals a substantial heterogeneity over time in selection on VR loci, suggesting that VR evolution is episodic.

Testing the adaptive radiation hypothesis for the lemurs of Madagascar

Lemurs, the diverse, endemic primates of Madagascar, are thought to represent a classic example of adaptive radiation. Based on the most complete phylogeny of living and extinct lemurs yet assembled, I tested predictions of adaptive radiation theory by estimating rates of speciation, extinction and adaptive phenotypic evolution. As predicted, lemur speciation rate exceeded that of their sister clade by nearly twofold, indicating the diversification dynamics of lemurs and mainland relatives may have been decoupled. Lemur diversification rates did not decline over time, however, as predicted by adaptive radiation theory. Optimal body masses diverged among dietary and activity pattern niches as lineages diversified into unique multidimensional ecospace. Based on these results, lemurs only partially fulfil the predictions of adaptive radiation theory, with phenotypic evolution corresponding to an 'early burst' of adaptive differentiation. The results must be interpreted with caution, however, because over the long evolutionary history of lemurs (approx. 50 million years), the 'early burst' signal of adaptive radiation may have been eroded by extinction.

Geogenetic patterns in mouse lemurs (genus Microcebus) reveal the ghosts of Madagascar's forests past

Phylogeographic analysis can be described as the study of the geological and climatological processes that have produced contemporary geographic distributions of populations and species. Here, we attempt to understand how the dynamic process of landscape change on Madagascar has shaped the distribution of a targeted clade of mouse lemurs (genus Microcebus) and, conversely, how phylogenetic and population genetic patterns in these small primates can reciprocally advance our understanding of Madagascar's prehuman environment. The degree to which human activity has impacted the natural plant communities of Madagascar is of critical and enduring interest. Today, the eastern rainforests are separated from the dry deciduous forests of the west by a large expanse of presumed anthropogenic grassland savanna, dominated by the Family Poaceae, that blankets most of the Central Highlands. Although there is firm consensus that anthropogenic activities have transformed the original vegetation through agricultural and pastoral practices, the degree to which closed-canopy forest extended from the east to the west remains debated. Phylogenetic and population genetic patterns in a five-species clade of mouse lemurs suggest that longitudinal dispersal across the island was readily achieved throughout the Pleistocene, apparently ending at ∼55 ka. By examining patterns of both inter- and intraspecific genetic diversity in mouse lemur species found in the eastern, western, and Central Highland zones, we conclude that the natural environment of the Central Highlands would have been mosaic, consisting of a matrix of wooded savanna that formed a transitional zone between the extremes of humid eastern and dry western forest types.

Integrating multiple evidences in taxonomy: species diversity and phylogeny of mustached bats (Mormoopidae: Pteronotus)

A phylogenetic systematic perspective is instrumental in recovering new species and their evolutionary relationships. The advent of new technologies for molecular and morphological data acquisition and analysis, allied to the integration of knowledge from different areas, such as ecology and population genetics, allows for the emergence of more rigorous, accurate and complete scientific hypothesis on species diversity. Mustached bats (genus Pteronotus) are a good model for the application of this integrative approach. They are a widely distributed and a morphologically homogeneous group, but comprising species with remarkable differences in their echolocation strategy and feeding behavior. The latest systematic review suggested six species with 17 subspecies in Pteronotus. Subsequent studies using discrete morphological characters supported the same arrangement. However, recent papers reported high levels of genetic divergence among conspecific taxa followed by bioacoustic and geographic agreement, suggesting an underestimated diversity in the genus. To date, no study merging genetic evidences and morphometric variation along the entire geographic range of this group has been attempted. Based on a comprehensive sampling including representatives of all current taxonomic units, we attempt to delimit species in Pteronotus through the application of multiple methodologies and hierarchically distinct datasets. The molecular approach includes six molecular markers from three genetic transmission systems; morphological investigations used 41 euclidean distances estimated through three-dimensional landmarks collected from 1628 skulls. The phylogenetic analysis reveals a greater diversity than previously reported, with a high correspondence among the genetic lineages and the currently recognized subspecies in the genus. Discriminant analysis of variables describing size and shape of cranial bones support the rising of the genetic groups to the specific status. Based on multiples evidences, we present an updated taxonomic arrangement composed by 16 extant species and a new and more robust phylogenetic hypothesis for the species included in the genus Pteronotus. Studies developed under such integrative taxonomic approach are timely for a deeper and wider comprehension of Neotropical diversity, representing the first step for answering broader questions on evolutionary and ecological aspects of Neotropical life history.

Multilocus Phylogeography and Species Delimitation in the Cumberland Plateau Salamander, Plethodon kentucki: Incongruence among Data Sets and Methods

Species are a fundamental unit of biodiversity, yet can be challenging to delimit objectively. This is particularly true of species complexes characterized by high levels of population genetic structure, hybridization between genetic groups, isolation by distance, and limited phenotypic variation. Previous work on the Cumberland Plateau Salamander, Plethodon kentucki, suggested that it might constitute a species complex despite occupying a relatively small geographic range. To examine this hypothesis, we sampled 135 individuals from 43 populations, and used four mitochondrial loci and five nuclear loci (5693 base pairs) to quantify phylogeographic structure and probe for cryptic species diversity. Rates of evolution for each locus were inferred using the multidistribute package, and time calibrated gene trees and species trees were inferred using BEAST 2 and *BEAST 2, respectively. Because the parameter space relevant for species delimitation is large and complex, and all methods make simplifying assumptions that may lead them to fail, we conducted an array of analyses. Our assumption was that strongly supported species would be congruent across methods. Putative species were first delimited using a Bayesian implementation of the GMYC model (bGMYC), Geneland, and Brownie. We then validated these species using the genealogical sorting index and BPP. We found substantial phylogeographic diversity using mtDNA, including four divergent clades and an inferred common ancestor at 14.9 myr (95% HPD: 10.8-19.7 myr). By contrast, this diversity was not corroborated by nuclear sequence data, which exhibited low levels of variation and weak phylogeographic structure. Species trees estimated a far younger root than did the mtDNA data, closer to 1.0 myr old. Mutually exclusive putative species were identified by the different approaches. Possible causes of data set discordance, and the problem of species delimitation in complexes with high levels of population structure and introgressive hybridization, are discussed.

Complete mitochondrial genome of the gray mouse lemur, Microcebus murinus (Primates, Cheirogaleidae)

We report the high-coverage complete mitochondrial genome sequence of the gray mouse lemur Microcebus murinus. The sequencing has been performed on an Illumina Hiseq 2500 platform, with a genome skimming strategy. The total length of this mitogenome is 16 963 bp, containing 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes and 1 non-coding region (D-loop region). The genome organization, nucleotide composition and codon usage are similar to those reported from other primate's mitochondrial genomes. The complete mitochondrial genome sequence reported here will be useful for comparative genomics studies in primates.

Recent advances in primate phylogenomics

The world of primate genomics is expanding rapidly in new and exciting ways owing to lowered costs and new technologies in molecular methods and bioinformatics. The primate order is composed of 78 genera and 478 species, including human. Taxonomic inferences are complex and likely a consequence of ongoing hybridization, introgression, and reticulate evolution among closely related taxa. Recently, we applied large-scale sequencing methods and extensive taxon sampling to generate a highly resolved phylogeny that affirms, reforms, and extends previous depictions of primate speciation. The next stage of research uses this phylogeny as a foundation for investigating genome content, structure, and evolution across primates. Ongoing and future applications of a robust primate phylogeny are discussed, highlighting advancements in adaptive evolution of genes and genomes, taxonomy and conservation management of endangered species, next-generation genomic technologies, and biomedicine.

Sky island diversification meets the multispecies coalescent - divergence in the spruce-fir moss spider (Microhexura montivaga, Araneae, Mygalomorphae) on the highest peaks of southern Appalachia

Microhexura montivaga is a miniature tarantula-like spider endemic to the highest peaks of the southern Appalachian mountains and is known only from six allopatric, highly disjunct montane populations. Because of severe declines in spruce-fir forest in the late 20th century, M. montivaga was formally listed as a US federally endangered species in 1995. Using DNA sequence data from one mitochondrial and seven nuclear genes, patterns of multigenic genetic divergence were assessed for six montane populations. Independent mitochondrial and nuclear discovery analyses reveal obvious genetic fragmentation both within and among montane populations, with five to seven primary genetic lineages recovered. Multispecies coalescent validation analyses [guide tree and unguided Bayesian Phylogenetics and Phylogeography (BPP), Bayes factor delimitation (BFD)] using nuclear-only data congruently recover six or seven distinct lineages; BFD analyses using combined nuclear plus mitochondrial data favour seven or eight lineages. In stark contrast to this clear genetic fragmentation, a survey of secondary sexual features for available males indicates morphological conservatism across montane populations. While it is certainly possible that morphologically cryptic speciation has occurred in this taxon, this system may alternatively represent a case where extreme population genetic structuring (but not speciation) leads to an oversplitting of lineage diversity by multispecies coalescent methods. Our results have clear conservation implications for this federally endangered taxon and illustrate a methodological issue expected to become more common as genomic-scale data sets are gathered for taxa found in naturally fragmented habitats.

Remarkable ancient divergences amongst neglected lorisiform primates

Lorisiform primates (Primates: Strepsirrhini: Lorisiformes) represent almost 10% of the living primate species and are widely distributed in sub-Saharan Africa and South/South-East Asia; however, their taxonomy, evolutionary history, and biogeography are still poorly understood. In this study we report the largest molecular phylogeny in terms of the number of represented taxa. We sequenced the complete mitochondrial cytochrome b gene for 86 lorisiform specimens, including ∼80% of all the species currently recognized. Our results support the monophyly of the Galagidae, but a common ancestry of the Lorisinae and Perodicticinae (family Lorisidae) was not recovered. These three lineages have early origins, with the Galagidae and the Lorisinae diverging in the Oligocene at about 30 Mya and the Perodicticinae emerging in the early Miocene. Our mitochondrial phylogeny agrees with recent studies based on nuclear data, and supports Euoticus as the oldest galagid lineage and the polyphyletic status of Galagoides. Moreover, we have elucidated phylogenetic relationships for several species never included before in a molecular phylogeny. The results obtained in this study suggest that lorisiform diversity remains substantially underestimated and that previously unnoticed cryptic diversity might be present within many lineages, thus urgently requiring a comprehensive taxonomic revision of this primate group. 

Agent-mediated spatial storage effect in heterogeneous habitat stabilizes competitive mouse lemur coexistence in Menabe Central, Western Madagascar

BACKGROUND

Spatio-temporal distribution patterns of species in response to natural and anthropogenic drivers provide insight into the ecological processes that determine community composition. We investigated determinants of ecological structure in a species assemblage of 4 closely related primate species of the family Cheirogaleidae (Microcebus berthae, Microcebus murinus, Cheirogaleus medius, Mirza coquereli) in western Madagascar by extensive line transect surveys across spatial and temporal heterogeneities with the specific goal of elucidating the mechanisms stabilizing competitive coexistence of the two mouse lemur species (Microcebus spp.).

RESULTS

Interspecific competition between the mouse lemurs was indicated by negative spatial associations in degraded habitat and by habitat partitioning along anthropogenic disturbance gradients during dry seasons with resource scarcity. In non-degraded habitat, intraguild predator M. coquereli, but not C. medius, was negatively associated with M. murinus on the population level, whereas its regional distribution overlapped spatially with that of M. berthae. The species' interspecific distribution pattern across spatial and temporal heterogeneities corresponded to predictions for agent-mediated coexistence and thus confirmed M. coquereli's stabilizing impact on the coexistence of mouse lemurs.

CONCLUSIONS

Interspecific interactions contribute to ecological structure in this cheirogaleid assemblage and determinants vary across spatio-temporal heterogeneities. Coexistence of Microcebus spp. is stabilized by an agent-mediated spatial storage effect: M. coquereli creates refuges from competition for M. berthae in intact habitat, whereas anthropogenic environments provide M. murinus with an escape from resource competition and intraguild predation. Species persistence in the assemblage therefore depends on the conservation of habitat content and context that stabilizing mechanisms rely on. Our large-scale population level approach did not allow for considering all potential functional and stochastic drivers of ecological structure, a key limitation that accounts for the large proportion of unexplained variance in our models.

MHC class II variation in a rare and ecological specialist mouse lemur reveals lower allelic richness and contrasting selection patterns compared to a generalist and widespread sympatric congener

The polymorphism of immunogenes of the major histocompatibility complex (MHC) is thought to influence the functional plasticity of immune responses and, consequently, the fitness of populations facing heterogeneous pathogenic pressures. Here, we evaluated MHC variation (allelic richness and divergence) and patterns of selection acting on the two highly polymorphic MHC class II loci (DRB and DQB) in the endangered primate Madame Berthe’s mouse lemur (Microcebus berthae). Using 454 pyrosequencing, we examined MHC variation in a total of 100 individuals sampled over 9 years in Kirindy Forest, Western Madagascar, and compared our findings with data obtained previously for its sympatric congener, the grey mouse lemur (Microcebus murinus). These species exhibit a contrasting ecology and demography that were expected to affect MHC variation and molecular signatures of selection. We found a lower allelic richness concordant with its low population density, but a similar level of allelic divergence and signals of historical selection in the rare feeding specialist M. berthae compared to the widespread generalist M. murinus. These findings suggest that demographic factors may exert a stronger influence than pathogen-driven selection on current levels of allelic richness in M. berthae. Despite a high sequence similarity between the two congeners, contrasting selection patterns detected at DQB suggest its potential functional divergence. This study represents a first step toward unravelling factors influencing the adaptive divergence of MHC genes between closely related but ecologically differentiated sympatric lemurs and opens new questions regarding potential functional discrepancy that would explain contrasting selection patterns detected at DQB.