Comparative and population mitogenomic analyses of Madagascar's extinct, giant ‘subfossil’ lemurs

Twin vocal folds as a novel evolutionary adaptation for vocal communications in lemurs

Primates have varied vocal repertoires to communicate with conspecifics and sometimes other species. The larynx has a central role in vocal source generation, where a pair of vocal folds vibrates to modify the air flow. Here, we show that Madagascan lemurs have a unique additional pair of folds in the vestibular region, parallel to the vocal folds. The additional fold has a rigid body of a vocal muscle branch and it is covered by a stratified squamous epithelium, equal to those of the vocal fold. Such anatomical features support the hypothesis that it also vibrates in a manner like the vibrations that occur in the vocal folds. To examine the acoustic function of the two pairs of folds, we made a silicone compound model to demonstrate that they can simultaneously vibrate to lower the fundamental frequency and increase vocal efficiency. Similar acoustic effects are achieved using different features of the larynx for the other primates, e.g., by vibrating multiple sets of ventricular folds in several species and further by an evolutionary modification of enlarged larynx in howler monkeys. Our multidisciplinary approaches found that these functions were acquired through a unique evolutionary adaptation of the twin vocal folds in Madagascan lemurs.

TeloBase: a community-curated database of telomere sequences across the tree of life

Discoveries over the recent decade have demonstrated the unexpected diversity of telomere DNA motifs in nature. However, currently available resources, 'Telomerase database' and 'Plant rDNA database', contain just fragments of all relevant literature published over decades of telomere research as they have a different primary focus and limited updates. To fill this gap, we gathered data about telomere DNA sequences from a thorough literature screen as well as by analysing publicly available NGS data, and we created TeloBase (http://cfb.ceitec.muni.cz/telobase/) as a comprehensive database of information about telomere motif diversity. TeloBase is supplemented by internal taxonomy utilizing popular on-line taxonomic resources that enables in-house data filtration and graphical visualisation of telomere DNA evolutionary dynamics in the form of heat tree plots. TeloBase avoids overreliance on administrators for future data updates by having a simple form and community-curation system for application and approval, respectively, of new telomere sequences by users, which should ensure timeliness of the database and topicality. To demonstrate TeloBase utility, we examined telomere motif diversity in species from the fungal genus Aspergillus, and discovered (TTTATTAGGG)n sequence as a putative telomere motif in the plant family Chrysobalanaceae. This was bioinformatically confirmed by analysing template regions of identified telomerase RNAs.

A comprehensive phylogenomic study unveils evolutionary patterns and challenges in the mitochondrial genomes of Carcharhiniformes: A focus on Triakidae

The complex evolutionary patterns in the mitochondrial genome (mitogenome) of the most species-rich shark order, the Carcharhiniformes (ground sharks) has led to challenges in the phylogenomic reconstruction of the families and genera belonging to the order, particularly the family Triakidae (houndsharks). The current state of Triakidae phylogeny remains controversial, with arguments for both monophyly and paraphyly within the family. We hypothesize that this variability is triggered by the selection of different a priori partitioning schemes to account for site and gene heterogeneity within the mitogenome. Here we used an extensive statistical framework to select the a priori partitioning scheme for inference of the mitochondrial phylogenomic relationships within Carcharhiniformes, tested site heterogeneous CAT + GTR + G4 models and incorporated the multi-species coalescent model (MSCM) into our analyses to account for the influence of gene tree discordance on species tree inference. We included five newly assembled houndshark mitogenomes to increase resolution of Triakidae. During the assembly procedure, we uncovered a 714 bp-duplication in the mitogenome of Galeorhinus galeus. Phylogenetic reconstruction confirmed monophyly within Triakidae and the existence of two distinct clades of the expanded Mustelus genus. The latter alludes to potential evolutionary reversal of reproductive mode from placental to aplacental, suggesting that reproductive mode has played a role in the trajectory of adaptive divergence. These new sequences have the potential to contribute to population genomic investigations, species phylogeography delineation, environmental DNA metabarcoding databases and, ultimately, improved conservation strategies for these ecologically and economically important species.

Deciphering the mandibular shape variation in a group of Malagasy primates using Fourier outline analysis

OBJECTIVES

Among living Malagasy primates, the family Lemuridae has previously been recognized as presenting a higher mandibular morphological variation than other families. We conducted a quantitative analysis of mandibular size and shape within the five genera (Lemur, Eulemur, Hapalemur, Prolemur, and Varecia) associated with a set of covariables that could explain this variation.

MATERIALS AND METHODS

We used Fourier outline analysis on the left hemimandible of 182 specimens covering the Lemuridae family. The influence of the phylogeny but also seven covariables (genus, diet, sex, sexual behavior, mating system, ecoregion, and forest type) on mandibular variation was examined using multivariate statistics and model selection.

RESULTS

Our results indicate that the high level of morphological variation within the family, associated with a phylogenetic effect and differences in diet, is due to a strong distinction between the genera Prolemur and Hapalemur and the other genera of the family. A second analysis, correcting this strong effect, indicates that mandibular shape variation is influenced not only by the phylogeny and the diet but by a combination of all the covariables.

DISCUSSION

The analysis of morphological variation is a powerful tool with major applications, both for the estimation of biological diversity and for the understanding of the fundamental parameters of species' ecology. Our work indicates that, if mandibular shape variation is mainly driven by dietary adaptation, other variables describing ecology and habitat should be considered and taken into account for an integrative understanding of species resources and the establishment of conservation measures.

The colonisation of Madagascar by land-bound vertebrates

Despite discussions extending back almost 160 years, the means by which Madagascar's iconic land vertebrates arrived on the island remains the focus of active debate. Three options have been considered: vicariance, range expansion across land bridges, and dispersal over water. The first assumes that a group (clade/lineage) occupied the island when it was connected with the other Gondwana landmasses in the Mesozoic. Causeways to Africa do not exist today, but have been proposed by some researchers for various times in the Cenozoic. Over-water dispersal could be from rafting on floating vegetation (flotsam) or by swimming/drifting. A recent appraisal of the geological data supported the idea of vicariance, but found nothing to justify the notion of past causeways. Here we review the biological evidence for the mechanisms that explain the origins of 28 of Madagascar's land vertebrate clades [two other lineages (the geckos Geckolepis and Paragehyra) could not be included in the analysis due to phylogenetic uncertainties]. The podocnemid turtles and typhlopoid snakes are conspicuous for they appear to have arisen through a deep-time vicariance event. The two options for the remaining 26 (16 reptile, five land-bound-mammal, and five amphibian), which arrived between the latest Cretaceous and the present, are dispersal across land bridges or over water. As these would produce very different temporal influx patterns, we assembled and analysed published arrival times for each of the groups. For all, a 'colonisation interval' was generated that was bracketed by its 'stem-old' and 'crown-young' tree-node ages; in two instances, the ranges were refined using palaeontological data. The synthesis of these intervals for all clades, which we term a colonisation profile, has a distinctive shape that can be compared, statistically, to various models, including those that assume the arrivals were focused in time. The analysis leads us to reject the various land bridge models (which would show temporal concentrations) and instead supports the idea of dispersal over water (temporally random). Therefore, the biological evidence is now in agreement with the geological evidence, as well as the filtered taxonomic composition of the fauna, in supporting over-water dispersal as the mechanism that explains all but two of Madagascar's land-vertebrate groups.

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.

Is species richness mediated by functional and genetic divergence? A global analysis in birds

Unravelling why species richness shows such dramatic spatial variation is an ongoing challenge. Common to many theories is that increasing species richness (e.g. with latitude) requires a compensatory trade-off on an axis of species' ecology. Spatial variation in species richness may also affect genetic diversity if large numbers of coexisting, related species result in smaller population sizes.Here, we test whether increasing species richness results in differential occupation of morphospace by the constituent species, or decreases species' genetic diversity. We test for two potential mechanisms of morphological accommodation: denser packing in ecomorphological space, and expansion of the space. We then test whether species differ in their nucleotide diversity depending on allopatry or sympatry with relatives, indicative of potential genetic consequences of coexistence that would reduce genetic diversity in sympatry. We ask these questions in a spatially explicit framework, using a global database of avian functional trait measurements in combination with >120,000 sequences downloaded from GenBank.We find that higher species richness within families is not systematically correlated with either packing in morphological space or overdispersion but, at the Class level, we find a general positive relationship between packing and species richness, but that points sampled in the tropics have comparatively greater packing than temperate ones relative to their species richness. We find limited evidence that geographical co-occurrence with closely related species or tropical distributions decreases nucleotide diversity of nuclear genes; however, this requires further analysis.Our results suggest that avian families can accumulate species regionally with minimal tradeoffs or cost, implying that external biotic factors do not limit species richness. Read the free Plain Language Summary for this article on the Journal blog.

The evolution of primate malaria parasites: A study on the origin and diversification of Plasmodium in lemurs

Among the primate malaria parasites, those found in lemurs have been neglected. Here, six Plasmodium lineages were detected in 169 lemurs. Nearly complete mitochondrial genomes (mtDNA, ≈6Kb) and apicoplast loci (≈6Kb) were obtained from these parasites and other Haemosporida species. Plasmodium spp. in lemurs are a diverse clade that shares a common ancestor with other primate parasites from continental Africa. Time-trees for the mtDNA were estimated under different scenarios, and the origin of the lemur clade coincides with the proposed time of their host species' most recent common ancestor (Lemuridae-Indriidae). A time tree with fewer taxa was estimated with mtDNA + Apicoplast loci. Those time estimates overlapped but were younger and had narrower credibility intervals than those from mtDNA alone. Importantly, the mtDNA + Apicoplast estimates that the clade including the most lethal malaria parasite in humans, Plasmodium falciparum, may have originated with Homininae (African apes). Finally, the phylogenetic congruence of the lemurs and their parasites was explored. A statistically significant scenario identified four cospeciation, two duplications, four transfer (host-switches), and zero loss events. Thus, the parasite species sampled in lemurs seem to be radiating with their hosts.

Nuclear and Mitochondrial Phylogenomics of the Sifakas Reveal Cryptic Variation in the Diademed Sifaka

The most comprehensive phylogenomic reconstruction to date was generated on all nominal taxa within the lemur genus Propithecus. Over 200 wild-caught individuals were included in this study to evaluate the intra and interspecific relationships across this genus. Ultraconserved Elements (UCEs) resulted in well-supported phylogenomic trees. Complete mitochondrial genomes (CMGs) largely agreed with the UCEs, except where a mitochondrial introgression was detected between one clade of the diademed sifaka (Propithecus diadema) and the Milne-Edwards sifaka (P. edwardsi). Additionally, the crowned (P. coronatus) and Von der Decken’s (P. deckeni) sifakas belonged to a single admixed lineage from UCEs. Further sampling across these two species is warranted to determine if our sampling represents a hybrid zone. P. diadema recovered two well-supported clades, which were dated and estimated as being ancient as the split between the Perrier’s (P. perrierii) and silky (P. candidus) sifakas. The reconstructed demographic history of the two clades also varied over time. We then modeled the modern ecological niches of the two cryptic P. diadema clades and found that they were significantly diverged (p < 0.01). These ecological differences result in a very limited zone of geographic overlap for the P. diadema clades (<60 km2). Niche models also revealed that the Onive River acts as a potential barrier to dispersal between P. diadema and P. edwardsi. Further taxonomic work is required on P. diadema to determine if its taxonomic status should be revised. This first genomic evaluation of the genus resolved the relationships between the taxa and the recovered cryptic diversity within one species.

Rocks and clocks revised: New promises and challenges in dating the primate tree of life

In recent years, multiple technological and methodological advances have increased our ability to estimate phylogenies, leading to more accurate dating of the primate tree of life. Here we provide an overview of the limitations and potentials of some of these advancements and discuss how dated phylogenies provide the crucial temporal scale required to understand primate evolution. First, we review new methods, such as the total-evidence dating approach, that promise a better integration between the fossil record and molecular data. We then explore how the ever-increasing availability of genomic-level data for more primate species can impact our ability to accurately estimate timetrees. Finally, we discuss more recent applications of mutation rates to date divergence times. We highlight example studies that have applied these approaches to estimate divergence dates within primates. Our goal is to provide a critical overview of these new developments and explore the promises and challenges of their application in evolutionary anthropology.

Comparative Analysis of the Vocal Repertoires of the Indri (Indri indri) and the Diademed Sifaka (Propithecus diadema)

Strepsirrhine vocalisations are extraordinarily diverse and cross-species comparisons are needed to explore how this variability evolved. We contributed to the investigation of primate acoustic diversity by comparing the vocal repertoire of two sympatric lemur species, Propithecus diadema and Indri indri. These diurnal species belong to the same taxonomic family and have similar activity patterns but different social structures. These features make them excellent candidates for an investigation of the phylogenetic, environmental, and social influence on primate vocal behavior. We recorded 3 P. diadema groups in 2014 and 2016. From 1,872 recordings we selected and assigned 3814 calls to 9 a priori call types, on the basis of their acoustic structure. We implemented a reproducible technique performing an acoustic feature extraction relying on frequency bins, t-SNE data reduction, and a hard-clustering analysis. We first quantified the vocal repertoire of P. diadema, finding consistent results for the 9 putatively identified call types. When comparing this repertoire with a previously published repertoire of I. indri, we found highly species-specific repertoires, with only 2% of the calls misclassified by species identity. The loud calls of the two species were very distinct, while the low-frequency calls were more similar. Our results pinpoint the role of phylogenetic history, social and environmental features on the evolution of communicative systems and contribute to a deeper understanding of the evolutionary roots of primate vocal differentiation. We conclude by arguing that standardized and reproducible techniques, like the one we employed, allow robust comparisons and should be prioritized in the future.

Habitat use by the island lemurs of Nosy Be, Madagascar

Madagascar's lemurs are threatened by forest loss, fragmentation, and degradation. Many species use flexible behaviors to survive in degraded habitat, but their ability to persist in very small areas may be limited. Insular lemurs, like those found on Nosy Be, an island off the northwestern coast of Madagascar, are at heightened risk of sudden population declines and extirpation. Nosy Be is home to two Critically Endangered species-the endemic Nosy Be sportive lemur (Lepilemur tymerlachsoni) and Claire's mouse lemur (Microcebus mamiratra)-as well as the Endangered black lemur (Eulemur macaco). Most of the remaining forest on Nosy Be is protected by the 862-ha Lokobe National Park. To document how Nosy Be lemurs use their restricted habitat, we conducted vegetation and reconnaissance surveys on 53 transects in and around Lokobe. We collected data on tree size, canopy cover, understory visibility, and elevation for 248 lemur sightings. We used a spatially explicit, multi-species occupancy model to investigate which forest-structure variables are important to lemurs. Our results represent some of the first data on habitat use by insular lemurs. Black lemurs preferred significantly larger trees and areas with less dense understory. They also occurred significantly less outside of Lokobe National Park, even when accounting for sampling effort and geography. The distributions of the sportive and mouse lemurs were not related to the forest structure variables we documented, but they did negatively predict each other-perhaps because their habitat requirements differ. These results also underscore the importance of the national park to protecting the black lemur population on Nosy Be and raise questions about what factors do influence the distribution of Nosy Be's smaller lemurs. Close monitoring is needed to prevent these populations and the ecosystem services they provide from disappearing, as have other island lemurs.

Molecular Adaptation to Folivory and the Conservation Implications for Madagascar’s Lemurs

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 &gt; 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.

Integrative approaches to dispersing science: A case study of March Mammal Madness

OBJECTIVES

Public engagement is increasingly viewed as an important pillar of scientific scholarship. For early career and established scholars, navigating the mosaic landscape of public education and science communication, noted for rapid "ecological" succession, can be daunting. Moreover, academics are characterized by diverse skills, motivations, values, positionalities, and temperaments that may differentially incline individuals to particular public translation activities.

METHODS

Here we briefly contextualize engagement activities within a scholarly portfolio, describe the use of one public education program-March Mammal Madness (MMM)- to highlight approaches to science communication, and explore essential elements and practical considerations for creating and sustaining outreach pursuits in tandem with other scholarly activities.

RESULTS

MMM, an annual simulated tournament of living and fossil animal taxa, has reached hundreds of thousands of learners since 2013. This program has provided a platform to communicate research findings from biology and anthropology and showcase numerous scholars in these fields. MMM has leveraged tournament devices to intentionally address topics of climate change, capitalist environmental degradation, academic sexism, and racist settler-colonialism. The tournament, however, has also perpetuated implicit biases that need disrupting.

CONCLUSIONS

By embracing reflexive, self-interrogative, and growth attitudes, the tournament organizers iteratively refine and improve this public science education program to better align our activities with our values and goals. Our experiences with MMM suggest that dispersing science is most sustainable when we combine ancestral adaptations for cooperation, community, and storytelling with good-natured competition in the context of shared experiences and shared values.

Evolutionary and phylogenetic insights from a nuclear genome sequence of the extinct, giant, "subfossil" koala lemur Megaladapis edwardsi

No endemic Madagascar animal with body mass >10 kg survived a relatively recent wave of extinction on the island. From morphological and isotopic analyses of skeletal "subfossil" remains we can reconstruct some of the biology and behavioral ecology of giant lemurs (primates; up to ∼160 kg) and other extraordinary Malagasy megafauna that survived into the past millennium. Yet, much about the evolutionary biology of these now-extinct species remains unknown, along with persistent phylogenetic uncertainty in some cases. Thankfully, despite the challenges of DNA preservation in tropical and subtropical environments, technical advances have enabled the recovery of ancient DNA from some Malagasy subfossil specimens. Here, we present a nuclear genome sequence (∼2× coverage) for one of the largest extinct lemurs, the koala lemur Megaladapis edwardsi (∼85 kg). To support the testing of key phylogenetic and evolutionary hypotheses, we also generated high-coverage nuclear genomes for two extant lemurs, Eulemur rufifrons and Lepilemur mustelinus, and we aligned these sequences with previously published genomes for three other extant lemurs and 47 nonlemur vertebrates. Our phylogenetic results confirm that Megaladapis is most closely related to the extant Lemuridae (typified in our analysis by E. rufifrons) to the exclusion of L. mustelinus, which contradicts morphology-based phylogenies. Our evolutionary analyses identified significant convergent evolution between M. edwardsi and an extant folivore (a colobine monkey) and an herbivore (horse) in genes encoding proteins that function in plant toxin biodegradation and nutrient absorption. These results suggest that koala lemurs were highly adapted to a leaf-based diet, which may also explain their convergent craniodental morphology with the small-bodied folivore Lepilemur.

Paleogenomics illuminates the evolutionary history of the extinct Holocene "horned" crocodile of Madagascar, Voay robustus

Ancient DNA is transforming our ability to reconstruct historical patterns and mechanisms shaping modern diversity and distributions. In particular, molecular data from extinct Holocene island faunas have revealed surprising biogeographic scenarios. Here, we recovered partial mitochondrial (mt) genomes for 1300-1400 year old specimens (n = 2) of the extinct "horned" crocodile, Voay robustus, collected from Holocene deposits in southwestern Madagascar. Phylogenetic analyses of partial mt genomes and tip-dated timetrees based on molecular, fossil, and stratigraphic data favor a sister group relationship between Voay and Crocodylus (true crocodiles). These well supported trees conflict with recent morphological systematic work that has consistently placed Voay within Osteolaeminae (dwarf crocodiles and kin) and provide evidence for likely homoplasy in crocodylian cranial anatomy and snout shape. The close relationship between Voay and Crocodylus lends additional context for understanding the biogeographic origins of these genera and refines competing hypotheses for the recent extinction of Voay from Madagascar.

Insights from macroevolutionary modelling and ancestral state reconstruction into the radiation and historical dietary ecology of Lemuriformes (Primates, Mammalia)

Background

Lemurs once rivalled the diversity of rest of the primate order despite thier confinement to the island of Madagascar. We test the adaptive radiation model of Malagasy lemur diversity using a novel combination of phylogenetic comparative methods and geometric methods for quantifying tooth shape.

Results

We apply macroevolutionary model fitting approaches and disparity through time analysis to dental topography metrics associated with dietary adaptation, an aspect of mammalian ecology which appears to be closely related to diversification in many clades. Metrics were also reconstructed at internal nodes of the lemur tree and these reconstructions were combined to generate dietary classification probabilities at internal nodes using discriminant function analysis. We used these reconstructions to calculate rates of transition toward folivory per million-year intervals. Finally, lower second molar shape was reconstructed at internal nodes by modelling the change in shape of 3D meshes using squared change parsimony along the branches of the lemur tree. Our analyses of dental topography metrics do not recover an early burst in rates of change or a pattern of early partitioning of subclade disparity. However, rates of change in adaptations for folivory were highest during the Oligocene, an interval of possible forest expansion on the island.

Conclusions

There was no clear phylogenetic signal of bursts of morphological evolution early in lemur history. Reconstruction of the molar morphologies corresponding to the ancestral nodes of the lemur tree suggest that this may have been driven by a shift toward defended plant resources, however. This suggests a response to the ecological opportunity offered by expanding forests, but not necessarily a classic adaptive radiation initiated by dispersal to Madagascar.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01793-x.

Primate phylogenomics uncovers multiple rapid radiations and ancient interspecific introgression

Our understanding of the evolutionary history of primates is undergoing continual revision due to ongoing genome sequencing efforts. Bolstered by growing fossil evidence, these data have led to increased acceptance of once controversial hypotheses regarding phylogenetic relationships, hybridization and introgression, and the biogeographical history of primate groups. Among these findings is a pattern of recent introgression between species within all major primate groups examined to date, though little is known about introgression deeper in time. To address this and other phylogenetic questions, here, we present new reference genome assemblies for 3 Old World monkey (OWM) species: Colobus angolensis ssp. palliatus (the black and white colobus), Macaca nemestrina (southern pig-tailed macaque), and Mandrillus leucophaeus (the drill). We combine these data with 23 additional primate genomes to estimate both the species tree and individual gene trees using thousands of loci. While our species tree is largely consistent with previous phylogenetic hypotheses, the gene trees reveal high levels of genealogical discordance associated with multiple primate radiations. We use strongly asymmetric patterns of gene tree discordance around specific branches to identify multiple instances of introgression between ancestral primate lineages. In addition, we exploit recent fossil evidence to perform fossil-calibrated molecular dating analyses across the tree. Taken together, our genome-wide data help to resolve multiple contentious sets of relationships among primates, while also providing insight into the biological processes and technical artifacts that led to the disagreements in the first place.

Combining three newly sequenced primate genomes with other published genomes, this study adapts a little-known method for detecting ancient introgression to genome-scale data, revealing multiple previously unknown examples of hybridization between primate species.

Mid-Cenozoic climate change, extinction, and faunal turnover in Madagascar, and their bearing on the evolution of lemurs

Background

Was there a mid-Cenozoic vertebrate extinction and recovery event in Madagascar and, if so, what are its implications for the evolution of lemurs? The near lack of an early and mid-Cenozoic fossil record on Madagascar has inhibited direct testing of any such hypotheses. We compare the terrestrial vertebrate fauna of Madagascar in the Holocene to that of early Cenozoic continental Africa to shed light on the probability of a major mid-Cenozoic lemur extinction event, followed by an “adaptive radiation” or recovery. We also use multiple analytic approaches to test competing models of lemur diversification and the null hypothesis that no unusual mid-Cenozoic extinction of lemurs occurred.

Results

Comparisons of the terrestrial vertebrate faunas of the early Cenozoic on continental Africa and Holocene on Madagascar support the inference that Madagascar suffered a major mid-Cenozoic extinction event. Evolutionary modeling offers some corroboration, although the level of support varies by phylogeny and model used. Using the lemur phylogeny and divergence dates generated by Kistler and colleagues, RPANDA and TESS offer moderate support for the occurrence of unusual extinction at or near the Eocene-Oligocene (E-O) boundary (34 Ma). TreePar, operating under the condition of obligate mass extinction, found peak diversification at 31 Ma, and low probability of survival of prior lineages. Extinction at the E-O boundary received greater support than other candidate extinctions or the null hypothesis of no major extinction. Using the lemur phylogeny and divergence dates generated by Herrera & Dàvalos, evidence for large-scale extinction diminishes and its most likely timing shifts to before 40 Ma, which fails to conform to global expectations.

Conclusions

While support for large-scale mid-Cenozoic lemur extinction on Madagascar based on phylogenetic modeling is inconclusive, the African fossil record does provide indirect support. Furthermore, a major extinction and recovery of lemuriforms during the Eocene-Oligocene transition (EOT) would coincide with other major vertebrate extinctions in North America, Europe, and Africa. It would suggest that Madagascar’s lemurs were impacted by the climate shift from “greenhouse” to “ice-house” conditions that occurred at that time. This could, in turn, help to explain some of the peculiar characteristics of the lemuriform clade.

Genetic population structure of endangered ring-tailed lemurs (Lemur catta) from nine sites in southern Madagascar

Madagascar's ring-tailed lemurs (Lemur catta) are experiencing rapid population declines due to ongoing habitat loss and fragmentation, as well as increasing exploitation for bushmeat and the illegal pet trade. Despite being the focus of extensive and ongoing behavioral studies, there is comparatively little known about the genetic population structuring of the species. Here, we present the most comprehensive population genetic analysis of ring-tailed lemurs to date from across their likely remaining geographic range. We assessed levels of genetic diversity and population genetic structure using multilocus genotypes for 106 adult individuals from nine geographically representative localities. Population structure and F ST analyses revealed moderate genetic differentiation with localities being geographically partitioned into northern, southern, western and also potentially central clusters. Overall genetic diversity, in terms of allelic richness and observed heterozygosity, was high in the species (AR = 4.74, H O = 0.811). In fact, it is the highest among all published lemur estimates to date. While these results are encouraging, ring-tailed lemurs are currently affected by ongoing habitat fragmentation and occur at lower densities in poorer quality habitats. The effects of continued isolation and fragmentation, coupled with climate-driven environmental instability, will therefore likely impede the long-term viability of the species.

The evolution of palate shape in the Lepilemur-Cheirogaleidae clade (Primates: Strepsirrhini)

OBJECTIVES

Phylogenies consistently group the folivorous Lepilemur species with the small-bodied insectivorous-frugivorous cheirogaleids. Juvenile lepilemurs and adult cheirogaleids share allometries in most aspects of skull morphology, except the palate. We investigated potential influences on palate shape in these taxa and several outgroups using geometric morphometrics.

MATERIALS AND METHODS

Our sample included representatives of four extant strepsirrhine families, Cheirogaleidae (including Lepilemurinae), Lemuridae, Indriidae, and Galagidae, and one subfossil Megaladapis. Our dataset comprised 32 landmarks collected from 397 specimens representing 15 genera and 28 species, and was analyzed using generalized procrustes analyses and between group principal component analysis. We explored the influence of size, phylogeny, diet, and the propagation of loud vocalizations on palate shape.

RESULTS

While congeneric species clustered within the morphospace, the phylomorphospace did not mirror molecular phylogenetic hypotheses of higher-order relationships. Four palate forms were distinguished within the Cheirogaleidae. Diet, strongly linked to body size, had the single greatest influence on palate shape. The production of long-distance advertisement calls was most often associated with positive scores on the PC1 axis.

DISCUSSION

Our results suggest that the extensive variation in palate shape among Cheirogaleidae is related to dietary shifts that accompanied changes in body size during the clade's radiation. Molecular phylogenies indicate that cheirogaleid diversification involved repeated dwarfing events, which in turn drove dietary shifts from ancestral folivory-frugivory to frugivory, gummivory, and faunivory in the descendant species. The elongated Lepilemur palate is probably related to accelerated eruption of the cheek teeth to render juveniles competent to shear leaves upon weaning.

A comparative study of gut microbiomes in captive nocturnal strepsirrhines

Feeding strategy and diet are increasingly recognized for their roles in governing primate gut microbiome (GMB) composition. Whereas feeding strategy reflects evolutionary adaptations to a host's environment, diet is a more proximate measure of food intake. Host phylogeny, which is intertwined with feeding strategy, is an additional, and often confounding factor that shapes GMBs across host lineages. Nocturnal strepsirrhines are an intriguing and underutilized group in which to examine the links between these three factors and GMB composition. Here, we compare GMB composition in four species of captive, nocturnal strepsirrhines with varying feeding strategies and phylogenetic relationships, but nearly identical diets. We use 16S rRNA sequences to determine gut bacterial composition. Despite similar husbandry conditions, including diet, we find that GMB composition varies significantly across host species and is linked to host feeding strategy and phylogeny. The GMBs of the omnivorous and the frugivorous species were significantly more diverse than were those of the insectivorous and exudativorous species. Across all hosts, GMBs were enriched for bacterial taxa associated with the macronutrient resources linked to the host's respective feeding strategy. Ultimately, the reported variation in microbiome composition suggests that the impacts of captivity and concurrent diet do not overshadow patterns of feeding strategy and phylogeny. As our understanding of primate GMBs progresses, populations of captive primates can provide insight into the evolution of host-microbe relationships, as well as inform future captive management protocols that enhance primate health and conservation.

Resolving the phylogenetic position of Darwin's extinct ground sloth (Mylodon darwinii) using mitogenomic and nuclear exon data

Mylodon darwinii is the extinct giant ground sloth named after Charles Darwin, who first collected its remains in South America. We have successfully obtained a high-quality mitochondrial genome at 99-fold coverage using an Illumina shotgun sequencing of a 12 880-year-old bone fragment from Mylodon Cave in Chile. Low level of DNA damage showed that this sample was exceptionally well preserved for an ancient subfossil, probably the result of the dry and cold conditions prevailing within the cave. Accordingly, taxonomic assessment of our shotgun metagenomic data showed a very high percentage of endogenous DNA with 22% of the assembled metagenomic contigs assigned to Xenarthra. Additionally, we enriched over 15 kb of sequence data from seven nuclear exons, using target sequence capture designed against a wide xenarthran dataset. Phylogenetic and dating analyses of the mitogenomic dataset including all extant species of xenarthrans and the assembled nuclear supermatrix unambiguously place Mylodon darwinii as the sister-group of modern two-fingered sloths, from which it diverged around 22 million years ago. These congruent results from both the mitochondrial and nuclear data support the diphyly of the two modern sloth lineages, implying the convergent evolution of their unique suspensory behaviour as an adaption to arboreality. Our results offer promising perspectives for whole-genome sequencing of this emblematic extinct taxon.

Primate tails: Ancestral state reconstruction and determinants of interspecific variation in primate tail length

OBJECTIVE

Living primates vary considerably in tail length-body size relation, ranging from tailless species to those where the tail is more than twice as long as the body. Because the general pattern and determinants of tail evolution remain incompletely known, we reconstructed evolutionary changes in relative tail length across all primates and sought to explain interspecific variation in this trait.

METHODS

We combined data on tail length, head-body length, intermembral index (IMI), habitat use, locomotion type, and range latitude for 340 species from published sources. We reconstructed the evolution of relative tail length to identify all independent cases of regime shifts on a primate phylogeny, using several methods based on Ornstein-Uhlenbeck (OU) models. Accounting for phylogeny, we also examined the effects of habitat, locomotion type, distance from the equator and IMI on interspecific variation in tail length-body size relation.

RESULTS

Primate tail length is not sexually dimorphic. A phylogenetic reconstruction allowing multiple optima explains the observed regime shifts best. During the evolutionary history of primates, relative tail length changed 50 times under an OU model. Specifically, relative tail length increased 26 and decreased 24 times. Most of these changes occurred among Old World primates. Among the variables tested here, interspecific variation in IMI and the difference between leaping and non-leaping locomotion explained interspecific variation in relative tail length: Evolutionary decreases in relative tail length are generally associated with an increase in IMI and an absence of leaping behavior.

CONCLUSIONS

Regime shifts for relative tail length in living primates occurred in concert with fundamental changes in IMI and a change from leaping to non-leaping locomotion, or vice versa. Exceptions from this general pattern are linked to the presence of a prehensile tail or specialized foraging strategies. Thus, the primate tail appears to have evolved in functional coordination with limb proportions, presumably to assist body balance.

An environmental DNA sampling method for aye-ayes from their feeding traces

Noninvasive sampling is an important development in population genetic monitoring of wild animals. Particularly, the collection of environmental DNA (eDNA) which can be collected without needing to encounter the target animal facilitates the genetic analysis of endangered species. One method that has been applied to these sample types is target capture and enrichment which overcomes the issue of high proportions of exogenous (nonhost) DNA from these lower quality samples. We tested whether target capture of mitochondrial DNA from sampled feeding traces of the aye-aye, an endangered lemur species would yield mitochondrial DNA sequences for population genetic monitoring. We sampled gnawed wood where aye-ayes excavate wood-boring insect larvae from trees. We designed RNA probes complementary to the aye-aye's mitochondrial genome and used these to isolate aye-aye DNA from other nontarget DNA in these samples. We successfully retrieved six near-complete mitochondrial genomes from two sites within the aye-aye's geographic range that had not been sampled previously. Our method demonstrates the application of next-generation molecular techniques to species of conservation concern. This method can likely be applied to alternative foraged remains to sample endangered species other than aye-ayes.

Fossil lemurs from Egypt and Kenya suggest an African origin for Madagascar's aye-aye

In 1967 G.G. Simpson described three partial mandibles from early Miocene deposits in Kenya that he interpreted as belonging to a new strepsirrhine primate, Propotto. This interpretation was quickly challenged, with the assertion that Propotto was not a primate, but rather a pteropodid fruit bat. The latter interpretation has not been questioned for almost half a century. Here we re-evaluate the affinities of Propotto, drawing upon diverse lines of evidence to establish that this strange mammal is a strepsirrhine primate as originally suggested by Simpson. Moreover, our phylogenetic analyses support the recognition of Propotto, together with late Eocene Plesiopithecus from Egypt, as African stem chiromyiform lemurs that are exclusively related to the extant aye-aye (Daubentonia) from Madagascar. Our results challenge the long-held view that all lemurs are descended from a single ancient colonization of Madagascar, and present an intriguing alternative scenario in which two lemur lineages dispersed from Africa to Madagascar independently, possibly during the later Cenozoic.

Phylogenomic Reconstruction of Sportive Lemurs (genus Lepilemur) Recovered from Mitogenomes with Inferences for Madagascar Biogeography

The family Lepilemuridae includes 26 species of sportive lemurs, most of which were recently described. The cryptic morphological differences confounded taxonomy until recent molecular studies; however, some species’ boundaries remain uncertain. To better understand the genus Lepilemur, we analyzed 35 complete mitochondrial genomes representing all recognized 26 sportive lemur taxa and estimated divergence dates. With our dataset we recovered 25 reciprocally monophyletic lineages, as well as an admixed clade containing Lepilemur mittermeieri and Lepilemur dorsalis. Using modern distribution data, an ancestral area reconstruction and an ecological vicariance analysis were performed to trace the history of diversification and to test biogeographic hypotheses. We estimated the initial split between the eastern and western Lepilemur clades to have occurred in the Miocene. Divergence of most species occurred from the Pliocene to the Pleistocene. The biogeographic patterns recovered in this study were better addressed with a combinatorial approach including climate, watersheds, and rivers. Generally, current climate and watershed hypotheses performed better for western and eastern clades, while speciation of northern clades was not adequately supported using the ecological factors incorporated in this study. Thus, multiple mechanisms likely contributed to the speciation and distribution patterns in Lepilemur.

Novel opsin gene variation in large-bodied, diurnal lemurs

Some primate populations include both trichromatic and dichromatic (red-green colour blind) individuals due to allelic variation at the X-linked opsin locus. This polymorphic trichromacy is well described in day-active New World monkeys. Less is known about colour vision in Malagasy lemurs, but, unlike New World monkeys, only some day-active lemurs are polymorphic, while others are dichromatic. The evolutionary pressures underlying these differences in lemurs are unknown, but aspects of species ecology, including variation in activity pattern, are hypothesized to play a role. Limited data on X-linked opsin variation in lemurs make such hypotheses difficult to evaluate. We provide the first detailed examination of X-linked opsin variation across a lemur clade (Indriidae). We sequenced the X-linked opsin in the most strictly diurnal and largest extant lemur, Indri indri, and nine species of smaller, generally diurnal indriids (Propithecus). Although nocturnal Avahi (sister taxon to Propithecus) lacks a polymorphism, at least eight species of diurnal indriids have two or more X-linked opsin alleles. Four rainforest-living taxa-I. indri and the three largest Propithecus species-have alleles not previously documented in lemurs. Moreover, we identified at least three opsin alleles in Indri with peak spectral sensitivities similar to some New World monkeys.

Evolution of Acidic Mammalian Chitinase Genes (CHIA) Is Related to Body Mass and Insectivory in Primates

Insects are an important food resource for many primates, but the chitinous exoskeletons of arthropods have long been considered to be indigestible by the digestive enzymes of most mammals. However, recently mice and insectivorous bats were found to produce the enzyme acidic mammalian chitinase (AMCase) to digest insect exoskeletons. Here, we report on the gene CHIA and its paralogs, which encode AMCase, in a comparative sample of nonhuman primates. Our results show that early primates likely had three CHIA genes, suggesting that insects were an important component of the ancestral primate diet. With some exceptions, most extant primate species retain only one functional CHIA paralog. The exceptions include two colobine species, in which all CHIA genes have premature stop codons, and several New World monkey species that retain two functional genes. The most insectivorous species in our sample also have the largest number of functional CHIA genes. Tupaia chinensis and Otolemur garnettii retain three functional CHIA paralogs, whereas Tarsius syrichta has a total of five, two of which may be duplications specific to the tarsier lineage. Selection analyses indicate that CHIA genes are under more intense selection in species with higher insect consumption, as well as in smaller-bodied species (<500 g), providing molecular support for Kay's Threshold, a well-established component of primatological theory which proposes that only small primates can be primarily insectivorous. These findings suggest that primates, like mice and insectivorous bats, may use the enzyme AMCase to digest the chitin in insect exoskeletons, providing potentially significant nutritional benefits.

A massively parallel strategy for STR marker development, capture, and genotyping

Short tandem repeat (STR) variants are highly polymorphic markers that facilitate powerful population genetic analyses. STRs are especially valuable in conservation and ecological genetic research, yielding detailed information on population structure and short-term demographic fluctuations. Massively parallel sequencing has not previously been leveraged for scalable, efficient STR recovery. Here, we present a pipeline for developing STR markers directly from high-throughput shotgun sequencing data without a reference genome, and an approach for highly parallel target STR recovery. We employed our approach to capture a panel of 5000 STRs from a test group of diademed sifakas (Propithecus diadema, n = 3), endangered Malagasy rainforest lemurs, and we report extremely efficient recovery of targeted loci—97.3–99.6% of STRs characterized with ≥10x non-redundant sequence coverage. We then tested our STR capture strategy on P. diadema fecal DNA, and report robust initial results and suggestions for future implementations. In addition to STR targets, this approach also generates large, genome-wide single nucleotide polymorphism (SNP) panels from flanking regions. Our method provides a cost-effective and scalable solution for rapid recovery of large STR and SNP datasets in any species without needing a reference genome, and can be used even with suboptimal DNA more easily acquired in conservation and ecological studies.

The contribution of admixture to primate evolution

Genome-wide data on genetic variation are now available for multiple primate species and populations, facilitating analyses of evolutionary history within and across taxa. One emerging theme from these studies involves the central role of admixture. Genomic data sets indicate that both ancient gene flow following initial taxonomic divergence and ongoing gene flow at current species boundaries are common. These findings are of particular interest given evidence for a complex history of admixture in our own lineage, including examples of ecologically driven adaptive introgression. Like other aspects of human biology, studies of nonhuman primates thus provide both comparative context and a living model for understanding admixture dynamics in hominins. We highlight several open questions that could be addressed in future work.

A new model for ancient DNA decay based on paleogenomic meta-analysis

The persistence of DNA over archaeological and paleontological timescales in diverse environments has led to a revolutionary body of paleogenomic research, yet the dynamics of DNA degradation are still poorly understood. We analyzed 185 paleogenomic datasets and compared DNA survival with environmental variables and sample ages. We find cytosine deamination follows a conventional thermal age model, but we find no correlation between DNA fragmentation and sample age over the timespans analyzed, even when controlling for environmental variables. We propose a model for ancient DNA decay wherein fragmentation rapidly reaches a threshold, then subsequently slows. The observed loss of DNA over time may be due to a bulk diffusion process in many cases, highlighting the importance of tissues and environments creating effectively closed systems for DNA preservation. This model of DNA degradation is largely based on mammal bone samples due to published genomic dataset availability. Continued refinement to the model to reflect diverse biological systems and tissue types will further improve our understanding of ancient DNA breakdown dynamics.

Disease-associated mitochondrial mutations and the evolution of primate mitogenomes

Several human diseases have been associated with mutations in mitochondrial genes comprising a set of confirmed and reported mutations according to the MITOMAP database. An analysis of complete mitogenomes across 139 primate species showed that most confirmed disease-associated mutations occurred in aligned codon positions and gene regions under strong purifying selection resulting in a strong evolutionary conservation. Only two confirmed variants (7.1%), coding for the same amino acids accounting for severe human diseases, were identified without apparent pathogenicity in non-human primates, like the closely related Bornean orangutan. Conversely, reported disease-associated mutations were not especially concentrated in conserved codon positions, and a large fraction of them occurred in highly variable ones. Additionally, 88 (45.8%) of reported mutations showed similar variants in several non-human primates and some of them have been present in extinct species of the genus Homo. Considering that recurrent mutations leading to persistent variants throughout the evolutionary diversification of primates are less likely to be severely damaging to fitness, we suggest that these 88 mutations are less likely to be pathogenic. Conversely, 69 (35.9%) of reported disease-associated mutations occurred in extremely conserved aligned codon positions which makes them more likely to damage the primate mitochondrial physiology.

Middle Pleistocene protein sequences from the rhinoceros genus Stephanorhinus and the phylogeny of extant and extinct Middle/Late Pleistocene Rhinocerotidae

Background

Ancient protein sequences are increasingly used to elucidate the phylogenetic relationships between extinct and extant mammalian taxa. Here, we apply these recent developments to Middle Pleistocene bone specimens of the rhinoceros genus Stephanorhinus. No biomolecular sequence data is currently available for this genus, leaving phylogenetic hypotheses on its evolutionary relationships to extant and extinct rhinoceroses untested. Furthermore, recent phylogenies based on Rhinocerotidae (partial or complete) mitochondrial DNA sequences differ in the placement of the Sumatran rhinoceros (Dicerorhinus sumatrensis). Therefore, studies utilising ancient protein sequences from Middle Pleistocene contexts have the potential to provide further insights into the phylogenetic relationships between extant and extinct species, including Stephanorhinus and Dicerorhinus.

Methods

ZooMS screening (zooarchaeology by mass spectrometry) was performed on several Late and Middle Pleistocene specimens from the genus Stephanorhinus, subsequently followed by liquid chromatography-tandem mass spectrometry (LC-MS/MS) to obtain ancient protein sequences from a Middle Pleistocene Stephanorhinus specimen. We performed parallel analysis on a Late Pleistocene woolly rhinoceros specimen and extant species of rhinoceroses, resulting in the availability of protein sequence data for five extant species and two extinct genera. Phylogenetic analysis additionally included all extant Perissodactyla genera (Equus, Tapirus), and was conducted using Bayesian (MrBayes) and maximum-likelihood (RAxML) methods.

Results

Various ancient proteins were identified in both the Middle and Late Pleistocene rhinoceros samples. Protein degradation and proteome complexity are consistent with an endogenous origin of the identified proteins. Phylogenetic analysis of informative proteins resolved the Perissodactyla phylogeny in agreement with previous studies in regards to the placement of the families Equidae, Tapiridae, and Rhinocerotidae. Stephanorhinus is shown to be most closely related to the genera Coelodonta and Dicerorhinus. The protein sequence data further places the Sumatran rhino in a clade together with the genus Rhinoceros, opposed to forming a clade with the black and white rhinoceros species.

Discussion

The first biomolecular dataset available for Stephanorhinus places this genus together with the extinct genus Coelodonta and the extant genus Dicerorhinus. This is in agreement with morphological studies, although we are unable to resolve the order of divergence between these genera based on the protein sequences available. Our data supports the placement of the genus Dicerorhinus in a clade together with extant Rhinoceros species. Finally, the availability of protein sequence data for both extinct European rhinoceros genera allows future investigations into their geographic distribution and extinction chronologies.

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.

Implications of lemuriform extinctions for the Malagasy flora

Madagascar's lemurs display a diverse array of feeding strategies with complex relationships to seed dispersal mechanisms in Malagasy plants. Although these relationships have been explored previously on a case-by-case basis, we present here the first comprehensive analysis of lemuriform feeding, to our knowledge, and its hypothesized effects on seed dispersal and the long-term survival of Malagasy plant lineages. We used a molecular phylogenetic framework to examine the mode and tempo of diet evolution, and to quantify the associated morphological space occupied by Madagascar's lemurs, both extinct and extant. Using statistical models and morphometric analyses, we demonstrate that the extinction of large-bodied lemurs resulted in a significant reduction in functional morphological space associated with seed dispersal ability. These reductions carry potentially far-reaching consequences for Malagasy ecosystems, and we highlight large-seeded Malagasy plants that appear to be without extant animal dispersers. We also identify living lemurs that are endangered yet occupy unique and essential dispersal niches defined by our morphometric analyses.

Bayesian molecular clock dating of species divergences in the genomics era

Five decades have passed since the proposal of the molecular clock hypothesis, which states that the rate of evolution at the molecular level is constant through time and among species. This hypothesis has become a powerful tool in evolutionary biology, making it possible to use molecular sequences to estimate the geological ages of species divergence events. With recent advances in Bayesian clock dating methodology and the explosive accumulation of genetic sequence data, molecular clock dating has found widespread applications, from tracking virus pandemics and studying the macroevolutionary process of speciation and extinction to estimating a timescale for life on Earth.

What did Hadropithecus eat, and why should paleoanthropologists care?

Over 40 years ago, Clifford Jolly noted different ways in which Hadropithecus stenognathus converged in its craniodental anatomy with basal hominins and with geladas. The Malagasy subfossil lemur Hadropithecus departs from its sister taxon, Archaeolemur, in that it displays comparatively large molars, reduced incisors and canines, a shortened rostrum, and thickened mandibular corpus. Its molars, however, look nothing like those of basal hominins; rather, they much more closely resemble molars of grazers such as Theropithecus. A number of tools have been used to interpret these traits, including dental microwear and texture analysis, molar internal and external morphology, and finite element analysis of crania. These tools, however, have failed to provide support for a simple dietary interpretation; whereas there is some consistency in the inferences they support, dietary inferences (e.g., that it was graminivorous, or that it specialized on hard objects) have been downright contradictory. Cranial shape may correlate poorly with diet. But a fundamental question remains unresolved: why do the various cranial and dental convergences exemplified by Hadropithecus, basal hominins, and Theropithecus exist? In this paper we review prior hypotheses regarding the diet of Hadropithecus. We then use stable carbon and nitrogen isotope data to elucidate this species' diet, summarizing earlier stable isotope analyses and presenting new data for lemurs from the central highlands of Madagascar, where Hadropithecus exhibits an isotopic signature strikingly different from that seen in other parts of the island. We offer a dietary explanation for these differences. Hadropithecus likely specialized neither on grasses nor hard objects; its staples were probably the succulent leaves of CAM plants. Nevertheless, aspects of prior hypotheses regarding the ecological significance of its morphology can be supported. Am. J. Primatol. 78:1098-1112, 2016. © 2015 Wiley Periodicals, Inc.

Sequence characterization and phylogeny analysis of the complete mitochondrial genome of verreaux's sifaka, Propithecus verreauxi (primates: indriidae)

The Verreaux's sifaka, Propithecus verreauxi, is a medium-sized lemur that inhabits in tropical dry lowland and montane forest. Here, we reported the complete mitochondrial genome sequence of this species. The mitogenome of P. verreauxi is 17 106 bp in length and composed of 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, 22 transfer RNA genes and a control region. The structure about gene order and composition is identical to that of the other lemur species and related genera. The overall base composition of the heavy strand in descending order is A (32.84%), T (27.13%), C (26.95%), and G (13.08%). Most of the genes are encoded on the heavy strand except for the NADH dehydrogenase subunit 6 (ND6) and eight tRNA genes. With the exception of COX3, which terminates with an incomplete stop codon (T-), all the other PCGs initiates with a traditional ATN start codon and ends with the typical mitochondrial stop codon (TAG/TAA/AGA). The phylogenetic tree constructed using the complete mitochondrial genome sequences of P. verreauxi together with 15 other closely related species with Neighbor-Joining (NJ) method shows that P. verreauxi is closer to P. coquereli in the phylogenetic relationship.