Search

Cross-species introgression can have significant impacts on phylogenomic reconstruction of species divergence events. Here, we used simulations to show how the presence of even a small amount of introgression can bias divergence time estimates when gene flow is ignored in the analysis. Using advances in analytical methods under the multispecies coalescent (MSC) model, we demonstrate that by accounting for incomplete lineage sorting and introgression using large phylogenomic data sets this problem can be avoided. The multispecies-coalescent-with-introgression (MSci) model is capable of accurately estimating both divergence times and ancestral effective population sizes, even when only a single diploid individual per species is sampled. We characterize some general expectations for biases in divergence time estimation under three different scenarios: 1) introgression between sister species, 2) introgression between non-sister species, and 3) introgression from an unsampled (i.e., ghost) outgroup lineage. We also conducted simulations under the isolation-with-migration (IM) model and found that the MSci model assuming episodic gene flow was able to accurately estimate species divergence times despite high levels of continuous gene flow. We estimated divergence times under the MSC and MSci models from two published empirical datasets with previous evidence of introgression, one of 372 target-enrichment loci from baobabs (Adansonia), and another of 1000 transcriptome loci from 14 species of the tomato relative, Jaltomata. The empirical analyses not only confirm our findings from simulations, demonstrating that the MSci model can reliably estimate divergence times but also show that divergence time estimation under the MSC can be robust to the presence of small amounts of introgression in empirical datasets with extensive taxon sampling. [divergence time; gene flow; hybridization; introgression; MSci model; multispecies coalescent]

The gut microbiome is gaining recognition for its role in primate nutrition, but we stand to benefit from microbiome comparisons across diverse hosts and environmental conditions. We compared gut microbiome structure in four lemur species from four phylogenetic lineages, including 9 individual mouse lemurs (Microcebus danfossi), 6 brown lemurs (Eulemur fulvus), 20 sifakas (Propithecus coquereli), and a single sportive lemur (Lepilemur grewcockorum). In northwestern Madagascar, these species are sympatric, but use different feeding strategies to cope with environmental challenges, including relying on tree gums and insects (mouse lemurs), and some vs. significant leaf matter (brown lemurs vs. sifakas and sportive lemurs). From one fecal sample collected per lemur in the dry season in the Anjajavy Forest, we determined gut microbiome diversity, variability, and membership via 16S rRNA sequencing. The lemurs harbored strongly species-specific gut microbiomes. Brown lemurs showed more diverse and generalized consortia; mouse lemurs, sifakas, and the sportive lemur had less diverse consortia with more distinct memberships. Consistent with their fallback foods, mouse lemur microbiomes included taxa putatively associated with gum and insect digestion, whereas those of sifakas and the sportive lemur showed stronger and distinct signatures of leaf fiber and secondary compound metabolism. These results point to feeding strategy, intertwined with host phylogeny, as a driver of gut microbiome composition, but highlight real-time dietary specificity as a contributing driver of microbiome diversity. While illuminating how gut microbiomes facilitate host nutrition on challenging foods, these results help explain how ecologically diverse primates living in sympatry may differentially cope with seasonal or stochastic lean times.

Duke Lemur Centre publication #1510; International audience; Madagascar's Central Highlands are largely composed of grasslands, interspersed with patches of forest. The historical perspective was that Madagascar's grasslands had anthropogenic origins, but emerging evidence suggests that grasslands were a component of the pre-human Central Highlands vegetation. Consequently, there is now vigorous debate regarding the extent to which these grasslands have expanded due to anthropogenic pressures. Here, we shed light on the temporal dynamics of Madagascar's vegetative composition by conducting a population genomic investigation of Goodman's mouse lemur (Microcebus lehilahytsara; Cheirogaleidae). These small-bodied primates occur both in Madagascar's eastern rainforests and in the Central Highlands, making them a valuable indicator species. Population divergences among forest-dwelling mammals will reflect changes to their habitat, including fragmentation, whereas patterns of post-divergence gene flow can reveal formerly wooded migration corridors. To explore these patterns, we used RADseq data to infer population genetic structure, demographic models of post-divergence gene flow, and population size change through time. The results offer evidence that open habitats are an ancient component of the Central Highlands, and that widespread forest fragmentation occurred naturally during a period of decreased precipitation near the last glacial maximum. Models of gene flow suggest that migration across the Central Highlands has been possible from the Pleistocene through the recent Holocene via riparian corridors. Though our findings support the hypothesis that Central Highland grasslands predate human arrival, we also find evidence for human-mediated population declines. This highlights the extent to which species imminently threatened by human-mediated deforestation may already be vulnerable from paleoclimatic conditions.

DNA sequences from three mitochondrial genes and one nuclear gene were analyzed to determine the phylogeny of the Malagasy primate family Lemuridae. Whether analyzed separately or in combination, the data consistently indicate that Eulemur species comprise a clade that is sister to a Lemur catta plus Hapalemur clade. The genus Varecia is basal to both. Resolution of cladogenic events within Eulemur was found to be extremely problematic with a total of six alternative arrangements offered by various data sets and weighting regimes. We attempt to determine the best arrangement of Eulemur taxa through a variety of character and taxon sampling strategies. Because our study includes all but one Eulemur species, increased taxon sampling is probably not an option for enhancing phylogenetic accuracy. We find, however, that the combined genetic data set is more robust to changes in taxon sample than are any of the individual data sets, suggesting that increased character sampling stabilizes phylogenetic resolution. Nonetheless, due to the difficult nature of the problem, we may have to accept certain aspects of Eulemur interrelationships as uncertain.

The lemurs of Madagascar include numerous species characterized by folivory across several families. Many extant lemuriform folivores exist in sympatry in Madagascar’s remaining forests. These species avoid feeding competition by adopting different dietary strategies within folivory, reflected in behavioral, morphological, and microbiota diversity across species. These conditions make lemurs an ideal study system for understanding adaptation to leaf-eating. Most folivorous lemurs are also highly endangered. The significance of folivory for conservation outlook is complex. Though generalist folivores may be relatively well equipped to survive habitat disturbance, specialist folivores occupying narrow dietary niches may be less resilient. Characterizing the genetic bases of adaptation to folivory across species and lineages can provide insights into their differential physiology and potential to resist habitat change. We recently reported accelerated genetic change in RNASE1, a gene encoding an enzyme (RNase 1) involved in molecular adaptation in mammalian folivores, including various monkeys and sifakas (genus Propithecus; family Indriidae). Here, we sought to assess whether other lemurs, including phylogenetically and ecologically diverse folivores, might show parallel adaptive change in RNASE1 that could underlie a capacity for efficient folivory. We characterized RNASE1 in 21 lemur species representing all five families and members of the three extant folivorous lineages: (1) bamboo lemurs (family Lemuridae), (2) sportive lemurs (family Lepilemuridae), and (3) indriids (family Indriidae). We found pervasive sequence change in RNASE1 across all indriids, a dN/dS value > 3 in this clade, and evidence for shared change in isoelectric point, indicating altered enzymatic function. Sportive and bamboo lemurs, in contrast, showed more modest sequence change. The greater change in indriids may reflect a shared strategy emphasizing complex gut morphology and microbiota to facilitate folivory. This case study illustrates how genetic analysis may reveal differences in functional traits that could influence species’ ecology and, in turn, their resilience to habitat change. Moreover, our results support the body of work demonstrating that not all primate folivores are built the same and reiterate the need to avoid generalizations about dietary guild in considering conservation outlook, particularly in lemurs where such diversity in folivory has probably led to extensive specialization via niche partitioning.

Sensory gene families are of special interest for both what they can tell us about molecular evolution and what they imply as mediators of social communication. The vomeronasal type-1 receptors (V1Rs) have often been hypothesized as playing a fundamental role in driving or maintaining species boundaries given their likely function as mediators of intraspecific mate choice, particularly in nocturnal mammals. Here, we employ a comparative genomic approach for revealing patterns of V1R evolution within primates, with a special focus on the small-bodied nocturnal mouse and dwarf lemurs of Madagascar (genera Microcebus and Cheirogaleus, respectively). By doubling the existing genomic resources for strepsirrhine primates (i.e. the lemurs and lorises), we find that the highly speciose and morphologically cryptic mouse lemurs have experienced an elaborate proliferation of V1Rs that we argue is functionally related to their capacity for rapid lineage diversification. Contrary to a previous study that found equivalent degrees of V1R diversity in diurnal and nocturnal lemurs, our study finds a strong correlation between nocturnality and V1R elaboration, with nocturnal lemurs showing elaborate V1R repertoires and diurnal lemurs showing less diverse repertoires. Recognized subfamilies among V1Rs show unique signatures of diversifying positive selection, as might be expected if they have each evolved to respond to specific stimuli. Furthermore, a detailed syntenic comparison of mouse lemurs with mouse (genus Mus) and other mammalian outgroups shows that orthologous mammalian subfamilies, predicted to be of ancient origin, tend to cluster in a densely populated region across syntenic chromosomes that we refer to as a V1R “hotspot.”

Germline mutations are the raw material for natural selection, driving species evolution and the generation of earth's biodiversity. Without this driver of genetic diversity, life on earth would stagnate. Yet, it is a double-edged sword. An excess of mutations can have devastating effects on fitness and population viability. It is therefore one of the great challenges of molecular ecology to determine the rate and mechanisms by which these mutations accrue across the tree of life. Advances in high-throughput sequencing technologies are providing new opportunities for characterizing the rates and mutational spectra within species and populations thus informing essential evolutionary parameters such as the timing of speciation events, the intricacies of historical demography, and the degree to which lineages are subject to the burdens of mutational load. Here, we will focus on both the challenge and promise of whole-genome comparisons among parents and their offspring from known pedigrees for the detection of germline mutations as they arise in a single generation. The potential of these studies is high, but the field is still in its infancy and much uncertainty remains. Namely, the technical challenges are daunting given that pedigree-based genome comparisons are essentially searching for needles in a haystack given the very low signal to noise ratio. Despite the challenges, we predict that rapidly developing methods for whole-genome comparisons hold great promise for integrating empirically derived estimates of de novo mutation rates and mutation spectra across many molecular ecological applications.

International audience; 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 km2 versus 9,250 and 50,700 km2, 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.

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 explosion has prompted understandable concern that there has been substantial oversplitting in the mouse lemur clade. Here, we take an integrative approach to investigate species diversity in two pairs of sister lineages that occur 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 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 them. 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. mittermeieriandM. lehilahytsarais not supported – a result that is particularly striking when using the genealogical discordance index (gdi). Non-sister 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 and calibrated with pedigree-based mutation rates and are considerably more recent than previously published fossil-calibrated concatenated likelihood estimates, however. We discuss the possible explanations for this discrepancy, noting that there are theoretical justifications for preferring the MSC estimates in this case.

Hibernation is an adaptive strategy some mammals use to survive highly seasonal or unpredictable environments. We present the first investigation on the transcriptomics of hibernation in a natural population of primate hibernators: Crossley's dwarf lemurs (Cheirogaleus crossleyi). Using capture-mark-recapture techniques to track the same animals over a period of 7 months in Madagascar, we used RNA-seq to compare gene expression profiles in white adipose tissue (WAT) during three distinct physiological states. We focus on pathway analysis to assess the biological significance of transcriptional changes in dwarf lemur WAT and, by comparing and contrasting what is known in other model hibernating species, contribute to a broader understanding of genomic contributions of hibernation across Mammalia. The hibernation signature is characterized by a suppression of lipid biosynthesis, pyruvate metabolism and mitochondrial-associated functions, and an accumulation of transcripts encoding ribosomal components and iron-storage proteins. The data support a key role of pyruvate dehydrogenase kinase isoenzyme 4 (PDK4) in regulating the shift in fuel economy during periods of severe food deprivation. This pattern of PDK4 holds true across representative hibernating species from disparate mammalian groups, suggesting that the genetic underpinnings of hibernation may be ancestral to mammals.

Bamboo specialization is one of the most extreme examples of convergent herbivory, yet it is unclear how this specific high-fiber diet might selectively shape the composition of the gut microbiome compared to host phylogeny. To address these questions, we used deep sequencing to investigate the nature and comparative impact of phylogenetic and dietary selection for specific gut microbial membership in three bamboo specialists-the bamboo lemur (Hapalemur griseus, Primates: Lemuridae), giant panda (Ailuropoda melanoleuca, Carnivora: Ursidae), and red panda (Ailurus fulgens, Carnivora: Musteloideadae), as well as two phylogenetic controls-the ringtail lemur (Lemur catta) and the Asian black bear (Ursus thibetanus). We detected significantly higher Shannon diversity in the bamboo lemur (10.029) compared to both the giant panda (8.256; p = 0.0001936) and the red panda (6.484; p = 0.0000029). We also detected significantly enriched bacterial taxa that distinguished each species. Our results complement previous work in finding that phylogeny predominantly governs high-level microbiome community structure. However, we also find that 48 low-abundance OTUs are shared among bamboo specialists, compared to only 8 OTUs shared by the bamboo lemur and its sister species, the ringtail lemur (Lemur catta, a generalist). Our results suggest that deep sequencing is necessary to detect low-abundance bacterial OTUs, which may be specifically adapted to a high-fiber diet. These findings provide a more comprehensive framework for understanding the evolution and ecology of the microbiome as well as the host.

Research on animal microbiomes is increasingly aimed at determining the evolutionary and ecological factors that govern host-microbiome dynamics, which are invariably intertwined and potentially synergistic. We present three empirical studies related to this topic, each of which relies on the diversity of Malagasy lemurs (representing a total of 19 species) and the comparative approach applied across scales of analysis. In Study 1, we compare gut microbial membership across 14 species in the wild to test the relative importance of host phylogeny and feeding strategy in mediating microbiome structure. Whereas host phylogeny strongly predicted community composition, the same feeding strategies shared by distant relatives did not produce convergent microbial consortia, but rather shaped microbiomes in host lineage-specific ways, particularly in folivores. In Study 2, we compare 14 species of wild and captive folivores, frugivores, and omnivores, to highlight the importance of captive populations for advancing gut microbiome research. We show that the perturbational effect of captivity is mediated by host feeding strategy and can be mitigated, in part, by modified animal management. In Study 3, we examine various scent-gland microbiomes across three species in the wild or captivity and show them to vary by host species, sex, body site, and a proxy of social status. These rare data provide support for the bacterial fermentation hypothesis in olfactory signal production and implicate steroid hormones as mediators of microbial community structure. We conclude by discussing the role of scale in comparative microbial studies, the links between feeding strategy and host-microbiome coadaptation, the underappreciated benefits of captive populations for advancing conservation research, and the need to consider the entirety of an animal's microbiota. Ultimately, these studies will help move the field from exploratory to hypothesis-driven research.

Population decline and fragmentation often lead to reduced genetic diversity and population differentiation. Habitat destruction throughout Madagascar has caused population decline and extinction of many endemic species. Lemur populations, including those of the largest extant lemur, Indri indri, have been fragmented into remaining forest patches. We assessed the level of genetic diversity in indri populations in three protected reserves by genotyping a total of 43 individuals at 17 microsatellite loci. Genetic diversity in terms of heterozygosity was high in all three reserves, with no differences between reserves. Population structure and F ST analyses revealed Analamazaotra Forest Station and the Torotorofotsy Conservation Area, which are separated by ca. 18 km to be genetically differentiated from each other with some admixture. Betampona Strict Nature Reserve, which is separated from the other reserves by ca. 130 km, exhibited clear population genetic differentiation, with no signs of admixture with the other reserves. Our genetic diversity estimates are similar to those for other Indridae in similar habitats and may reflect past rather than current population processes, given that populations have declined recently. Our results suggest that Betampona may be genetically isolated and that it is important to maintain gene flow between remaining populations to prevent loss of genetic diversity for the future conservation of Indri indri.

Kimura's neutral theory of molecular evolution has been essential to virtually every advance in evolutionary genetics, and by extension, is foundational to the field of conservation genetics. Conservation genetics utilizes the key concepts of neutral theory to identify species and populations at risk of losing evolutionary potential by detecting patterns of inbreeding depression and low effective population size. In turn, this information can inform the management of organisms and their habitat providing hope for the long-term preservation of both. We expand upon Avise's "inventorial" and "functional" categories of conservation genetics by proposing a third category that is linked to the coalescent and that we refer to as "process-driven." It is here that connections between Kimura's theory and conservation genetics are strongest. Process-driven conservation genetics can be especially applied to large genomic data sets to identify patterns of historical risk, such as population bottlenecks, and accordingly, yield informed intuitions for future outcomes. By examining inventorial, functional, and process-driven conservation genetics in sequence, we assess the progression from theory, to data collection and analysis, and ultimately, to the production of hypotheses that can inform conservation policies.

The rainforest biome of eastern Madagascar is renowned for its extraordinary biodiversity and restricted distribution ranges of many species, whereas the arid western region of the island is relatively species poor. We provide insight into the biogeography of western Madagascar by analyzing a multilocus phylogeographic dataset assembled for an amphibian, the widespread Malagasy bullfrog, Laliostoma labrosum. We find no cryptic species in L. labrosum (maximum 1.1% pairwise genetic distance between individuals in the 16S rRNA gene) attributable to considerable gene flow at the regional level as shown by genetic admixture in both mtDNA and three nuclear loci, especially in central Madagascar. Low breeding site fidelity, viewed as an adaptation to the unreliability of standing pools of freshwater in dry and seasonal environments, and a ubiquitous distribution within its range may underlie overall low genetic differentiation. Moreover, reductions in population size associated with periods of high aridity in western Madagascar may have purged DNA variation in this species. The mtDNA gene tree revealed seven major phylogroups within this species, five of which show mostly non-overlapping distributions. The nested positions of the northern and central mtDNA phylogroups imply a southwestern origin for all extant mtDNA lineages in L. labrosum. The current phylogeography of this species and paleo-distributions of major mtDNA lineages suggest five potential refugia in northern, western and southwestern Madagascar, likely the result of Pleistocene range fragmentation during drier and cooler climates. Lineage sorting in mtDNA and nuclear loci highlighted a main phylogeographic break between populations north and south of the Sambirano region, suggesting a role of the coastal Sambirano rainforest as a barrier to gene flow. Paleo-species distribution models and dispersal networks suggest that the persistence of some refugial populations was mainly determined by high population connectivity through space and time.

Geospatial modeling is one of the most powerful tools available to conservation biologists for estimating current species ranges of Earth's biodiversity. Now, with the advantage of predictive climate models, these methods can be deployed for understanding future impacts on threatened biota. Here, we employ predictive modeling under a conservative estimate of future climate change to examine impacts on the future abundance and geographic distributions of Malagasy lemurs. Using distribution data from the primary literature, we employed ensemble species distribution models and geospatial analyses to predict future changes in species distributions. Current species distribution models (SDMs) were created within the BIOMOD2 framework that capitalizes on ten widely used modeling techniques. Future and current SDMs were then subtracted from each other, and areas of contraction, expansion, and stability were calculated. Model overprediction is a common issue associated Malagasy taxa. Accordingly, we introduce novel methods for incorporating biological data on dispersal potential to better inform the selection of pseudo-absence points. This study predicts that 60% of the 57 species examined will experience a considerable range of reductions in the next seventy years entirely due to future climate change. Of these species, range sizes are predicted to decrease by an average of 59.6%. Nine lemur species (16%) are predicted to expand their ranges, and 13 species (22.8%) distribution sizes were predicted to be stable through time. Species ranges will experience severe shifts, typically contractions, and for the majority of lemur species, geographic distributions will be considerably altered. We identify three areas in dire need of protection, concluding that strategically managed forest corridors must be a key component of lemur and other biodiversity conservation strategies. This recommendation is all the more urgent given that the results presented here do not take into account patterns of ongoing habitat destruction relating to human activities.

DNA quantity can be a hindrance in ecological and evolutionary research programmes due to a range of factors including endangered status of target organisms, available tissue type, and the impact of field conditions on preservation methods. A potential solution to low-quantity DNA lies in whole genome amplification (WGA) techniques that can substantially increase DNA yield. To date, few studies have rigorously examined sequence bias that might result from WGA and next-generation sequencing of nonmodel taxa. To address this knowledge deficit, we use multiple displacement amplification (MDA) and double-digest RAD sequencing on the grey mouse lemur (Microcebus murinus) to quantify bias in genome coverage and SNP calls when compared to raw genomic DNA (gDNA). We focus our efforts in providing baseline estimates of potential bias by following manufacturer's recommendations for starting DNA quantities (>100 ng). Our results are strongly suggestive that MDA enrichment does not introduce systematic bias to genome characterization. SNP calling between samples when genotyping both de-novo and with a reference genome are highly congruent (>98%) when specifying a minimum threshold of 20X stack depth to call genotypes. Relative genome coverage is also similar between MDA and gDNA, and allelic dropout is not observed. SNP concordance varies based on coverage threshold, with 95% concordance reached at ~12X coverage genotyping de-novo and ~7X coverage genotyping with the reference genome. These results suggest that MDA may be a suitable solution for next-generation molecular ecological studies when DNA quantity would otherwise be a limiting factor.

Humans first arrived on Madagascar only a few thousand years ago. Subsequent habitat destruction and hunting activities have had significant impacts on the island's biodiversity, including the extinction of megafauna. For example, we know of 17 recently extinct ‘subfossil’ lemur species, all of which were substantially larger (body mass ∼11–160 kg) than any living population of the ∼100 extant lemur species (largest body mass ∼6.8 kg). We used ancient DNA and genomic methods to study subfossil lemur extinction biology and update our understanding of extant lemur conservation risk factors by i) reconstructing a comprehensive phylogeny of extinct and extant lemurs, and ii) testing whether low genetic diversity is associated with body size and extinction risk. We recovered complete or near-complete mitochondrial genomes from five subfossil lemur taxa, and generated sequence data from population samples of two extinct and eight extant lemur species. Phylogenetic comparisons resolved prior taxonomic uncertainties and confirmed that the extinct subfossil species did not comprise a single clade. Genetic diversity estimates for the two sampled extinct species were relatively low, suggesting small historical population sizes. Low genetic diversity and small population sizes are both risk factors that would have rendered giant lemurs especially susceptible to extinction. Surprisingly, among the extant lemurs, we did not observe a relationship between body size and genetic diversity. The decoupling of these variables suggests that risk factors other than body size may have as much or more meaning for establishing future lemur conservation priorities.

Background The de novo assembly of repeat-rich mammalian genomes using only high-throughput short read sequencing data typically results in highly fragmented genome assemblies that limit downstream applications. Here, we present an iterative approach to hybrid de novo genome assembly that incorporates datasets stemming from multiple genomic technologies and methods. We used this approach to improve the gray mouse lemur (Microcebus murinus) genome from early draft status to a near chromosome-scale assembly. Methods We used a combination of advanced genomic technologies to iteratively resolve conflicts and super-scaffold the M. murinus genome. Results We improved the M. murinus genome assembly to a scaffold N50 of 93.32 Mb. Whole genome alignments between our primary super-scaffolds and 23 human chromosomes revealed patterns that are congruent with historical comparative cytogenetic data, thus demonstrating the accuracy of our de novo scaffolding approach and allowing assignment of scaffolds to M. murinus chromosomes. Moreover, we utilized our independent datasets to discover and characterize sequences associated with centromeres across the mouse lemur genome. Quality assessment of the final assembly found 96% of mouse lemur canonical transcripts nearly complete, comparable to other published high-quality reference genome assemblies. Conclusions We describe a new assembly of the gray mouse lemur (Microcebus murinus) genome with chromosome-scale scaffolds produced using a hybrid bioinformatic and sequencing approach. The approach is cost effective and produces superior results based on metrics of contiguity and completeness. Our results show that emerging genomic technologies can be used in combination to characterize centromeres of non-model species and to produce accurate de novo chromosome-scale genome assemblies of complex mammalian genomes. Electronic supplementary material The online version of this article (doi:10.1186/s12915-017-0439-6) contains supplementary material, which is available to authorized users.

Microbiologists often evaluate microbial community dynamics by formulating functional hypotheses based on ecological processes. Indeed, many of the methods and terms currently used to describe animal microbiomes derive from ecology and evolutionary biology. As our understanding of the composition and functional dynamics of "the microbiome" grows, we increasingly refer to the host as an ecosystem within which microbial processes play out. Even so, an ecosystem service framework that extends to the context of the host has thus far been lacking. Here, we argue that ecosystem services are a useful framework with which to consider the value of microbes to their hosts. We discuss those "microbiome services" in the specific context of the mammalian gut, providing a context from which to develop new hypotheses and to evaluate microbial functions in future studies and novel systems.

Primates have long been a test case for the development of phylogenetic methods for divergence time estimation. Despite a large number of studies, however, the timing of origination of crown Primates relative to the K-Pg boundary and the timing of diversification of the main crown groups remain controversial. Here we analysed a dataset of 372 taxa (367 Primates and 5 outgroups, 61 thousand base pairs) that includes nine complete primate genomes (3.4 million base pairs). We systematically explore the effect of different interpretations of fossil calibrations and molecular clock models on primate divergence time estimates. We find that even small differences in the construction of fossil calibrations can have a noticeable impact on estimated divergence times, especially for the oldest nodes in the tree. Notably, choice of molecular rate model (auto-correlated or independently distributed rates) has an especially strong effect on estimated times, with the independent rates model producing considerably more ancient estimates for the deeper nodes in the phylogeny. We implement thermodynamic integration, combined with Gaussian quadrature, in the program MCMCTree, and use it to calculate Bayes factors for clock models. Bayesian model selection indicates that the auto-correlated rates model fits the primate data substantially better, and we conclude that time estimates under this model should be preferred. We show that for eight core nodes in the phylogeny, uncertainty in time estimates is close to the theoretical limit imposed by fossil uncertainties. Thus, these estimates are unlikely to be improved by collecting additional molecular sequence data. All analyses place the origin of Primates close to the K-Pg boundary, either in the Cretaceous or straddling the boundary into the Palaeogene.

Vomeronasal receptor genes have frequently been invoked as integral to the establishment and maintenance of species boundaries among mammals due to the elaborate one-to-one correspondence between semiochemical signals and neuronal sensory inputs. Here, we report the most extensive sample of vomeronasal receptor class 1 (V1R) sequences ever generated for a diverse yet phylogenetically coherent group of mammals, the tooth-combed primates (suborder Strepsirrhini). Phylogenetic analysis confirms our intensive sampling from a single V1R subfamily, apparently unique to the strepsirrhine primates. We designate this subfamily as V1Rstrep. The subfamily retains extensive repertoires of gene copies that descend from an ancestral gene duplication that appears to have occurred prior to the diversification of all lemuriform primates excluding the basal genus Daubentonia (the aye-aye). We refer to the descendent clades as V1Rstrep-α and V1Rstrep-β. Comparison of the two clades reveals different amino acid compositions corresponding to the predicted ligand-binding site and thus potentially to altered functional profiles between the two. In agreement with previous studies of the mouse lemur (genus, Microcebus), the majority of V1Rstrep gene copies appear to be intact and under strong positive selection, particularly within transmembrane regions. Finally, despite the surprisingly high number of gene copies identified in this study, it is nonetheless probable that V1R diversity remains underestimated in these nonmodel primates and that complete characterization will be limited until high-coverage assembled genomes are available. peerReviewed

Numerous hypotheses have been proposed for the historical processes governing the rich endemism of Madagas- car's biodiversity. The 'watershed model' suggests that drier climates in the recent geological past have resulted in the contraction of forests around major watersheds, thereby defining areas of endemism. We test whether this hypothesis explains phylogeographical patterns in a dry forest-dependent rodent, Eliurus myoxinus, an endemic species widely distributed through western Madagascar. We sequenced the mitochondrial cytochrome b locus and nuclear introns of the β-fibrinogen and the growth hormone receptor genes for E. myoxinus. Using a parametric bootstrapping approach, we tested whether the mitochondrial gene tree data fit expectations of local differentiation given the watershed model. We additionally estimated population differentiation and historical demographic parameters, and reconstructed the spatial history of E. myoxinus to highlight spatial and temporal patterns of differentiation. The data do not support the watershed model as a clear explanation for the genetic patterns of diversity within extant E. myoxinus populations. We find striking patterns of latitudinal genetic structure within western Madagascar, and indicate possible roles for environmental and ecological gradients along this axis in generating phylogeographical diversity. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 110, 500-517.

The number of newly discovered Malagasy vertebrate taxa has multiplied in recent years, emphasizing the importance of complete taxon sampling for phylogenetics, biogeography, functional ecology, and conservation biology, especially in such a biodiversity hotspot. In particular, the diversity of extant lemurs is much higher than previously thought, and we have yet to comprehend fully the full extent of lemuriform biodiversity. A recent genetic analysis of mtDNA and nDNA sequence data in Malagasy mouse lemurs revealed the existence of several novel mtDNA clades based on new field sampling. These geographically defined and previously unrecognized mtDNA clades corresponded precisely to patterns of population structure revealed in the analysis of the nDNA data, thus confirming their evolutionary divergence from other mouse lemur clades. Two of these independently evolving lineages correspond to specimens that were collected by us in the Marolambo and Manantantely/Ivorona regions. Here we summarize the genetic evidence and report on the morphometric and external characteristics of these animals, formally describing them as new species. This report thus brings the number of currently recognized and described mouse lemur species to 20. The forests in which these mouse lemurs were discovered have been heavily degraded in the past decade, prompting the classification of one of the new species as Endangered by the IUCN, even before its formal description. As with several other newly described lemur species, immediate field studies and appropriate conservation actions are therefore urgent.

Deforestation and a changing climate threaten the health and survival of lemurs in Madagascar. An important component of lemur health, parasite infection can reduce fitness and survival outcomes. Future lemur parasite richness, abundance and distribution may be highly influenced by climate change. Current knowledge of lemur parasites is narrow in geographic and temporal scope, with sampling at a limited number of sites, and thus far, there have been no attempts to assess the effects of climate change on lemur parasite distributions. We used geospatial tools to predict the distributions of six lemur parasites of high frequency and pathogenic potential. We then assessed how anticipated climate shifts in Madagascar may alter the distributions of these lemur parasites in the future. Under current climate conditions, we found that the focal parasites exhibited widespread potential distributions across Madagascar, covering 12–26% of surface land area and 40–86% of forested area. Our analyses also showed that parasites responded differently to projected climate changes, with shifts ranging from a contraction of current distributions by 7% to an expansion of 60%. A predicted net expansion in parasite distribution may expose naive lemur hosts to new parasites, which could have a profound effect on lemur health. Those parasites with the greatest potential for harmful effects are predicted to experience the largest expansion in range. Predicting these changing distributions will be critical for assessing population health, improving protected area design, preparing for reintroduction efforts and addressing potential parasite risk in lemurs, humans and domestic animals.

— Resolving phylogenetic relationships within the economically important genus Coffea (Rubiaceae) has proven to be difficult due to low levels of plastid and nuclear ITS sequence divergence. The recent identification of a large number of conserved, single-copy, putatively orthologous genes (COSII) provides a unique opportunity to apply regions of the nuclear genome to phylogenetic studies of Euasterid species. We expand a previously published plastid data set of Coffea with the addition of three new plastid regions and a single COSII nuclear locus. Significant phylogenetic incongruence is identified between the plastid and nuclear COSII data sets, and a combined analysis is performed after removal of incongruent taxa. Phylogenetic results from plastid, nuclear, and combined plastid/nuclear data sets allow the novel identification of early diverging lineages in Coffea. Specifically, the data show that a group of predominately Lower-Guinea/Congolian Coffea species form a clade that is sister to the ...

Understanding the remarkably high species diversity and levels of endemism found among Madagascar’s flora and fauna has been the focus of many studies. One hypothesis that has received much attention proposes that Quaternary climate fluctuations spurred diversification. While spatial patterns of distribution and phylogenetic relationships can provide support for biogeographic predictions, temporal estimates of divergence are required to determine the fit of these geospatial patterns to climatic or biogeographic mechanisms. We use multilocus DNA sequence data to test whether divergence times among Malagasy iguanid lizards of the subfamily Oplurinae are compatible with a hypotheses of Pliocene‐Pleistocene diversification. We estimate the oplurine species tree and associated divergence times under a relaxed-clock model. In addition, we examine the phylogeographic structure and population divergence times within two sister species of Oplurus primarily distributed in the north-west and south-west of Madagascar (Oplurus cuvieriandOpluruscyclurus,respectively).Wefindthatdivergenceeventsamongoplurine lineages occurred in the Oligocene and Miocene and are thus far older and incompatible with the hypothesis that recent climate fluctuations are related to current species diversity. However,thetimingofintraspecificdivergencesandspatialpatternsofpopulationgenetic structure within O. cuvieri and O. cyclurus suggest a role for both intrinsic barriers and recentclimate fluctuationsat population-leveldivergences.Integratinginformationacross spatial and temporal scales allows us to identify and better understand the mechanisms generating patterns diversity.

The systematics and speciation literature is rich with discussion relating to the potential for gene tree/species tree discordance. Numerous mechanisms have been proposed to generate discordance, including differential selection, longbranch attraction, gene duplication, genetic introgression, and/or incomplete lineage sorting. For speciose clades in which divergence has occurred recently and rapidly, recovering the true species tree can be particularly problematic due to incomplete lineage sorting. Unfortunately, the availability of multilocus or ‘‘phylogenomic’’ data sets does not simply solve the problem, particularly when the data are analyzed with standard concatenation techniques. In our study, we conduct a phylogenetic study for a nearly complete species sample of the dwarf and mouse lemur clade, Cheirogaleidae. Mouse lemurs (genus, Microcebus) have been intensively studied over the past decade for reasons relating to their high level of cryptic species diversity, and although there has been emerging consensus regarding the evolutionary diversity contained within the genus, there is no agreement as to the inter-specific relationships within the group. We attempt to resolve cheirogaleid phylogeny, focusing especially on the mouse lemurs, by employing a large multilocus data set. We compare the results of Bayesian concordance methods with those of standard gene concatenation, finding that though concatenation yields the strongest results as measured by statistical support, these results are found to be highly misleading. By employing an approach where individual alleles are treated as operational taxonomic units, we show that phylogenetic results are substantially influenced by the selection of alleles in the concatenation process.

We investigate the genetic relationships between purported island species of Pteropus fruit bat (Megachiroptera) from the western Indian Ocean islands using mitochondrial DNA sequencing in order to infer the pattern of colonisation of this biogeographic region. Most significantly, our genetic data questions the current taxonomic assignment based on morphology of many of the island species and subspecies, suggesting instead that many of the western Indian Ocean islands harbour 'races' of P. giganteus from mainland India. Our results strongly argue against a single colonisation event from mainland Asia. Evidence is presented for three colonisation events; the first to the western-most extremity of their range (Comoros and Pemba Island), the second to Rodrigues Island; and a third giving rise to the remaining extant island taxa, the latter two events occurring relatively recently and rapidly.

The chiropteran family Myzopodidae is endemic to Madagascar and is characterized by several unique morphologies, such as sessile adhesive discs on the thumb and sole. A new species, Myzopoda schliemanni, was recently described from western Madagascar that is morphologically distinct and geographically disjunct from the eastern species, M. aurita, the only other member of this family. Geographic variation within Myzopoda has only recently been studied at the morphological level and has never been addressed at the genetic level. We used a combination of phylogenetic, coalescent, and population genetic analyses to characterize the speciation history of Myzopoda and to clarify current and former patterns of gene flow within and between Myzopoda. Mitochondrial DNA sequences were used to determine whether genetic data support the morphologically distinct species M. schliemanni, to infer the distribution of the common ancestor of extant Myzopoda, to estimate effective population sizes (Ne) and levels of migration between species, and to determine patterns of population structure within species. Phylogenetic and network analyses revealed the existence of 4 well-supported clades in Myzopoda, but could not resolve relationships among those clades. Divergent haplotypes within species may result from either recent gene flow between the 2 species or more likely from incomplete lineage sorting. Multiple coalescent-based methodologies produced concordant estimates of Ne for Myzopoda, but conflicting signals for migration between the species, probably reflecting differences in the underlying models used by the methods. We found significant genetic structure within M. aurita, but no correlation with geography. This pattern may result from recent gene flow facilitated by expansion of Ravenala stands, an important day-roost tree for Myzopoda, associated with anthropogenic deforestation and the opening up of new habitat for members of this genus.

One of the most obvious physiological changes accompanying seasonal heterothermy in mammals is a fattening stage preceding periods of resource scarcity. This phenomenon reflects the interplay of both diet and physiology. Though the accrual of fat stores is known to be essential for overwintering in some species, the influence of diet on the physiology of torpor is not fully understood. Results from captive studies in heterothermic rodents and marsupials have indicated that when autumn diets are enriched with polyunsaturated fatty acids (PUFAs), animals receiving these diets experience deeper and more frequent torpor bouts than their counterparts receiving a control diet. Our study investigates this potential effect of dietary composition in animals that use daily torpor rather than prolonged torpor (i.e., hibernation). In so doing, we investigate the degree to which dietary effects on torpor are restricted to cold-adapted rodents and marsupials, or are a more general feature of mammalian heterothermy. We examined the effects of a PUFA diet and a control diet on the thermoregulation of one of the few species of primates known to use daily torpor: the grey mouse lemur (Microcebus murinus). Though the results of this study are largely inconclusive regarding the impact of dietary manipulations on torpor frequency and duration, we nonetheless find that the propensity of animals to enter torpor is directly influenced by age and seasonal changes in body mass, and thus reflect important physiological aspects of flexible thermoregulatory responses.

Hibernation is a complex physiological response that some mammalian species employ to evade energetic demands. Previous work in mammalian hibernators suggests that hibernation is activated not by a set of genes unique to hibernators, but by differential expression of genes that are present in all mammals. This question of universal genetic mechanisms requires further investigation and can only be tested through additional investigations of phylogenetically dispersed species. To explore this question, we use RNA-Seq to investigate gene expression dynamics as they relate to the varying physiological states experienced throughout the year in a group of primate hibernators—Madagascar’s dwarf lemurs (genus Cheirogaleus). In a novel experimental approach, we use longitudinal sampling of biological tissues as a method for capturing gene expression profiles from the same individuals throughout their annual hibernation cycle. We identify 90 candidate genes that have variable expression patterns when comparing two active states (Active 1 and Active 2) with a torpor state. These include genes that are involved in metabolic pathways, feeding behavior, and circadian rhythms, as might be expected to correlate with seasonal physiological state changes. The identified genes appear to be critical for maintaining the health of an animal that undergoes prolonged periods of metabolic depression concurrent with the hibernation phenotype. By focusing on these differentially expressed genes in dwarf lemurs, we compare gene expression patterns in previously studied mammalian hibernators. Additionally, by employing evolutionary rate analysis, we find that hibernation-related genes do not evolve under positive selection in hibernating species relative to nonhibernators., This work was supported by grants from the Society of Integrative and Comparative Biology and Sigma Xi to S.L.F., a National Science Foundation Doctoral Research Improvement Grant (NSF-BCS 1455809 to A.D.Y. and S.L.F.), a grant from the Duke University Institute for Genome Science and Policy to A.D.Y., and funding from the Ministerio de Economía y Competitividad of the Spanish Government (BIO2009-08160 to J.L.V.C. and BFU2012-36820 to J.L.V.C. and M.M.A.); and Fundació ICREA to M.M.A.

New applications of genetic data to questions of historical biogeography have revolutionized our understanding of how organisms have come to occupy their present distributions. Phylogenetic methods in combination with divergence time estimation can reveal biogeographical centres of origin, differentiate between hypotheses of vicariance and dispersal, and reveal the directionality of dispersal events. Despite their power, however, phylogenetic methods can sometimes yield patterns that are compatible with multiple, equally well-supported biogeographical hypotheses. In such cases, additional approaches must be integrated to differentiate among conflicting dispersal hypotheses. Here, we use a synthetic approach that draws upon the analytical strengths of coalescent and population genetic methods to augment phylogenetic analyses in order to assess the biogeographical history of Madagascar's Triaenops bats (Chiroptera: Hipposideridae). Phylogenetic analyses of mitochondrial DNA sequence data for Malagasy and east African Triaenops reveal a pattern that equally supports two competing hypotheses. While the phylogeny cannot determine whether Africa or Madagascar was the centre of origin for the species investigated, it serves as the essential backbone for the application of coalescent and population genetic methods. From the application of these methods, we conclude that a hypothesis of two independent but unidirectional dispersal events from Africa to Madagascar is best supported by the data.

Implementation of the coalescent model in a Bayesian framework is an emerging strength in genetically based species delimitation studies. By providing an objective measure of species diagnosis, these methods represent a quantitative enhancement to the analysis of multilocus data, and complement more traditional methods based on phenotypic and ecological characteristics. Recognized as two species 20 years ago, mouse lemurs (genus Microcebus) now comprise more than 20 species, largely diagnosed from mtDNA sequence data. With each new species description, enthusiasm has been tempered with scientific scepticism. Here, we present a statistically justified and unbiased Bayesian approach towards mouse lemur species delimitation. We perform validation tests using multilocus sequence data and two methodologies: (i) reverse-jump Markov chain Monte Carlo sampling to assess the likelihood of different models defined a priori by a guide tree, and (ii) a Bayes factor delimitation test that compares different species-tree models without a guide tree. We assess the sensitivity of these methods using randomized individual assignments, which has been used in bpp studies, but not with Bayes factor delimitation tests. Our results validate previously diagnosed taxa, as well as new species hypotheses, resulting in support for three new mouse lemur species. As the challenge of multiple researchers using differing criteria to describe diversity is not unique to Microcebus, the methods used here have significant potential for clarifying diversity in other taxonomic groups. We echo previous studies in advocating that multiple lines of evidence, including use of the coalescent model, should be trusted to delimit new species.

Madagascar is one of the world's hottest biodiversity hot spots due to its diverse, endemic, and highly threatened biota. This biota shows a distinct signature of evolution in isolation, both in the high levels of diversity within lineages and in the imbalance of lineages that are represented. For example, chameleon diversity is the highest of any place on Earth, yet there are no salamanders. These biotic enigmas have inspired centuries of speculation relating to the mechanisms by which Madagascar's biota came to reside there. The two most probable causal factors are Gondwanan vicariance and/or Cenozoic dispersal. By reviewing a comprehensive sample of phylogenetic studies of Malagasy biota, we find that the predominant pattern is one of sister group relationships to African taxa. For those studies that include divergence time analysis, we find an overwhelming indication of Cenozoic origins for most Malagasy clades. We conclude that most of the present-day biota of Madagascar is comprised of the descendents of Cenozoic dispersers, predominantly with African origins.

The living and extinct giant tortoises of the western Indian Ocean islands have been extensively studied in recent times (Bourn et al., 1999; Austin and Arnold, 2001; Palkovacs et al., 2002; Austin et al., 2003; Palkovacs et al., 2003). Currently, indigenous wild populations of Indian Ocean giant tortoises (genus Dipsochelys) occur only on the island of Aldabra, northwest of Madagascar, with most wild tortoises having become extinct on the Seychelles as early as 1800, largely through overexploitation by European sailors and settlers (Stoddardt et al., 1979; Austin et al., 2003). Recent molecular genetic studies on Dipsochelys suggest that there is very little genetic variation among extant Indian Ocean tortoises. Palkovacs et al. (2002) analyzed three mitochondrial genes, 12S and 16S ribosomal RNA (rRNA) and cytochrome b (cyt-b), and found no variation among captive tortoises assigned to three Dipsochelys species. Austin et al. (2003) looked at cyt-b sequences from modern and museum skins of Dipsochelys from the islands of Aldabra and the Seychelles and found ex

The increasing use of mitochondrial DNA (mtDNA) to explore and test species limits among morphologically similar species is potentially compromised by phenomena poorly reflective of organismal history and speciation, including (but not limited to) stochastic lineage sorting and gene flow. In situations where molecular data are only available from a single gene or linkage partition (e.g. mtDNA), corroboration of suspected species boundaries should be sought from independent lines of evidence, such as morphology. Recent attempts to delimit species using mtDNA and morphology have either implicitly or explicitly ignored the possibility that distinct species can occur in direct sympatry throughout much of their range, presumably because such situations are believed to be rare. We examined phylogenetic relationships within the long-tailed shrew tenrecs (Mammalia: Tenrecidae; Microgale spp.) from Madagascar. Current taxonomy recognizes two broadly sympatric species, though as many as six have been described. Given that alpha taxonomy within shrew tenrecs has been controversial, and that patterns of morphological variation can be especially difficult to assess for this group, some authors have suggested that additional cryptic species may exist. To examine this possibility, we conducted a phylogenetic study using the mitochondrial NADH dehydrogenase subunit 2 gene and a morphometric analysis of 29 craniodental, postcranial, and external measurements from a broad geographical sample of long-tailed shrew tenrecs. The two data sets were nearly perfectly congruent in identifying four groups that can be classified as species, thereby doubling the currently recognized number of species. We present previously unrecognized distributional evidence consistent with our conclusions and provide an empirical example of how a revised understanding of species limits alters inferences of geographic variation and species coexistence, particularly with respect to fine-scale habitat partitioning. The results of this study suggest that certain species pairs, previously assumed to be single species occupying broad elevational ranges, are actually reproductively isolated units that are partitioning their environment along elevational lines. © 2004 The Linnean Society of London, Biological Journal of the Linnean Society, 2004, 83, 1‐22. ADDITIONAL KEYWORDS: Madagascar ‐ mitochondrial DNA ‐ morphometry ‐ phylogeography ‐ species limits.

Divergence time and substitution rate are seriously confounded in phylogenetic analysis, making it difficult to estimate divergence times when the molecular clock (rate constancy among lineages) is violated. This problem can be alleviated to some extent by analyzing multiple gene loci simultaneously and by using multiple calibration points. While different genes may have different patterns of evolutionary rate change, they share the same divergence times. Indeed, the fact that each gene may violate the molecular clock differently leads to the advantage of simultaneous analysis of multiple loci. Multiple calibration points provide the means for characterizing the local evolutionary rates on the phylogeny. In this paper, we extend previous likelihood models of local molecular clock for estimating species divergence times to accommodate multiple calibration points and multiple genes. Heterogeneity among different genes in evolutionary rate and in substitution process is accounted for by the models. We apply the likelihood models to analyze two mitochondrial protein-coding genes, cytochrome oxidase II and cytochrome b, to estimate divergence times of Malagasy mouse lemurs and related outgroups. The likelihood method is compared with the Bayes method of Thorne et al. (1998, Mol. Biol. Evol. 15:1647-1657), which uses a probabilistic model to describe the change in evolutionary rate over time and uses the Markov chain Monte Carlo procedure to derive the posterior distribution of rates and times. Our likelihood implementation has the drawbacks of failing to accommodate uncertainties in fossil calibrations and of requiring the researcher to classify branches on the tree into different rate groups. Both problems are avoided in the Bayes method. Despite the differences in the two methods, however, data partitions and model assumptions had the greatest impact on date estimation. The three codon positions have very different substitution rates and evolutionary dynamics, and assumptions in the substitution model affect date estimation in both likelihood and Bayes analyses. The results demonstrate that the separate analysis is unreliable, with dates variable among codon positions and between methods, and that the combined analysis is much more reliable. When the three codon positions were analyzed simultaneously under the most realistic models using all available calibration information, the two methods produced similar results. The divergence of the mouse lemurs is dated to be around 7-10 million years ago, indicating a surprisingly early species radiation for such a morphologically uniform group of primates.

The phylogeny of Rosoideae was investigated using 44 species. Here we report new sequence data from the chloroplast trnL/F region as well as an increased sample of species. The analysis of these new data, along with previously used data from the nuclear ribosomal internal transcribed spacers (ITS), significantly increased resolution as well as confidence for Rosoideae phylogeny. Using both Bayesian inference and parsimony methods, we conducted analyses on the data sets separately and in combination. The resulting phylogenies are congruent with all well‐supported clades of Rosoideae found in previous analyses of ITS or rbcL data. The support for these and other clades is improved, and we consider several clades to be supported well enough to be named. The following clades are given phylogenetic definitions: Sanguisorbeae and its subclades Agrimoniinae and Sanguisorbinae, Potentilleae and its subclades Fragariinae and Potentilla, Roperculina (Rosa + Sanguisorbeae + Potentilleae), and Sanpotina (Sanguisorbea...

Recent morphological and molecular phylogenetic studies of mouse lemurs (Microcebus) living in the western and southern regions of Madagascar have shown that specific diversity had been considerably underestimated. In large part, this underestimate was due to the lack of sufficient specimens from given localities to assess properly the level of phenotypic variation within and between populations. The accurate delineation of specific boundaries has no doubt been confounded by the diminutive size, nocturnal habits, and subtle morphological variation characteristic of mouse lemurs, which can make field identification of individuals problematic. We illustrate the use of molecular phylogenetic analysis to reveal reproductive isolation in two sympatric mouse lemur species, Microcebus murinus and M. griseorufus. Their documentation in the Berenty Private Reserve in the extreme south of Madagascar verifies the historically-broad distribution of Microcebus griseorufus, a species recently resurrected from synonomy.

Although numerous studies have attempted to find single unifying mechanisms for generating Madagascar's unique flora and fauna, little consensus has been reached regarding the relative importance of climatic, geologic and ecological processes as catalysts of diversification of the region's unique biota. Rather, recent work has shown that both biological and physical drivers of diversification are best analysed in a case-by-case setting with attention focused on the ecological and life-history requirements of the specific phylogenetic lineage under investigation. Here, we utilize a comprehensive analytical approach to examine evolutionary drivers and elucidate the biogeographic history of Malagasy plated lizards (Zonosaurinae). Data from three genes are combined with fossil information to construct time-calibrated species trees for zonosaurines and their African relatives, which are used to test alternative diversification hypotheses. Methods are utilized for explicitly incorporating phylogenetic uncertainty into downstream analyses. Species distribution models are created for 14 of 19 currently recognized species, which are then used to estimate spatial patterns of species richness and endemicity. Spatially explicit analyses are employed to correlate patterns of diversity with both topographic heterogeneity and climatic stability through geologic time. We then use inferred geographic ranges to estimate the biogeographic history of zonosaurines within each of Madagascar's major biomes. Results suggest constant Neogene and Quaternary speciation with divergence from the African most recent common ancestor ~30 million years ago when oceanic currents and African rivers facilitated dispersal. Spatial patterns of diversity appear concentrated along coastal regions of northern and southern Madagascar. We find no relationship between either topographic heterogeneity or climatic stability and patterns of diversity. Ancestral state reconstructions suggest that western dry forests were important centres of origin with recent invasion into spiny and rain forest. These data highlight the power of combining multilocus phylogenetic and spatially explicit analyses for testing alternative diversification hypotheses within Madagascar's unique biota and more generally, particularly as applied to phylogenetically and biologically constrained systems.

A recent survey of the high-mountain zone of the Madagascar Parc National (PN) d'Andringitra revealed the presence of an apparently isolated troop of the ring-tailed lemur Lemur catta. These animals display phenotypic and ecological characteristics that are unusual for the monotypic genus Lemur, thus raising the possibility that they are members of a different undescribed species. We present analyses of two mitochondrial genes to test the hypothesis that L. catta from Andringitra should be considered a distinct species. The results indicate that taxonomic revision is not warranted under the expectations of the phylogenetic, coalescent, or biological species concepts. Rather, the genetic patterns observed among the Andringitra and lowland mitochondrial haplotypes are consistent with those expected for a single species.

High-throughput sequencing platforms are generating massive amounts of genomic data from nonmodel species, and these data sets are valuable resources that can be mined to advance a number of research areas. An example is the growing amount of transcriptome data that allow for examination of gene expression in nonmodel species. Here, we show how publicly available transcriptome data from nonmodel primates can be used to design novel research focused on immunogenomics. We mined transcriptome data from the world's most endangered group of primates, the lemurs of Madagascar, for sequences corresponding to immunoglobulins. Our results confirmed homology between strepsirrhine and haplorrhine primate immunoglobulins and allowed for high-throughput sequencing of expressed antibodies (Ig-seq) in Coquerel's sifaka (Propithecus coquereli). Using both Pacific Biosciences RS and Ion Torrent PGM sequencing, we performed Ig-seq on two individuals of Coquerel's sifaka. We generated over 150 000 sequences of expressed antibodies, allowing for molecular characterization of the antigen-binding region. Our analyses suggest that similar VDJ expression patterns exist across all primates, with sequences closely related to the human VH 3 immunoglobulin family being heavily represented in sifaka antibodies. Moreover, the antigen-binding region of sifaka antibodies exhibited similar amino acid variation with respect to haplorrhine primates. Our study represents the first attempt to characterize sequence diversity of the expressed antibody repertoire in a species of lemur. We anticipate that methods similar to ours will provide the framework for investigating the adaptive immune response in wild populations of other nonmodel organisms and can be used to advance the burgeoning field of eco-immunology.

DNA sequences of the complete cytochrome b gene are shown to contain robust phylogenetic signal for the strepsirrhine primates (i.e., lemurs and lorises). The phylogeny derived from these data conforms to other molecular studies of strepsirrhine relationships despite the fact that uncorrected nucleotide distances are high for nearly all intrastrepsirrhine comparisons, with most in the 15%-20% range. Cytochrome b sequences support the hypothesis that Malagasy lemuriforms and Afro-Asian lorisiforms each comprise clades that share a sister-group relationship. A study (Adkins and Honeycutt 1994) of the cytochrome c oxidase subunit II (COII) gene placed one Malagasy primate (Daubentonia) at the base of the strepsirrhine clade, thereby suggesting a diphyletic Lemuriformes. The reanalysis of COII third-position transversions, either alone or in combination with cytochrome b third-position transversions, however, yields a tree that is congruent with phylogenetic hypotheses derived from cytochrome b and other genetic data sets.

We report the development of 18, single-copy, anonymous nuclear loci from the Malagasy plated lizard Zonosaurus madagascariensis. More than 140 clones from a genomic library were examined and 38 potential loci tested across both closely and distantly related lizards. Of the 18 loci reported here, more than half (10) work in closely related zonosaurines although only one successfully amplified a homologous fragment in the distantly related iguanid (Oplurus). Sequences of these loci revealed a high frequency of single nucleotide polymorphisms, supporting previous reports of high levels of intraspecific variation in lizards.

Anthropogenic habitat disturbance impairs ecosystem health by fragmenting forested areas, introducing environmental contamination, and reducing the quality of habitat resources. The effect of this disturbance on wildlife health is of particular concern in Madagascar, one of the world's biodiversity hotspots, where anthropogenic pressures on the environment remain high. Despite the conservation importance of threatened lemur populations in Madagascar, few data exist on the effects of anthropogenic disturbance on lemur health. To examine these impacts, indri (Indri indri) populations were evaluated from two forest reserves that differ in their exposure to anthropogenic disturbance. We compared the health status of 36 indri individuals from two sites: one population from a protected, undisturbed area of lowland evergreen humid forest and the other population from a reserve exposed to frequent tourism and forest degradation. Comparison of indri health parameters between sites suggests an impact of anthropogenic disturbance, including significant differences in leukocyte count and differential, 12 serum parameters, 6 trace minerals, and a higher diversity of parasites, with a significant difference in the presence of the louse, Trichophilopterus babakotophilus. These data suggest that indri living in disturbed forests may experience physiological changes and increased susceptibility to parasitism, which may ultimately impair reproductive success and survival. Am. J. Primatol. 73:632–642, 2011. © 2011 Wiley-Liss, Inc.