DNA from extinct giant lemurs links archaeolemurids to extant indriids

Estimation of DNA Degradation in Archaeological Human Remains

The evaluation of the integrity and quantity of DNA extracted from archaeological human remains is a fundamental step before using the latest generation sequencing techniques in the study of evolutionary processes. Ancient DNA is highly fragmented and chemically modified; therefore, the present study aims to identify indices that can allow the identification of potentially amplifiable and sequenceable DNA samples, reducing failures and research costs. Ancient DNA was extracted from five human bone remains from the archaeological site of Amiternum L'Aquila, Italy dating back to the 9th-12th century and was compared with standard DNA fragmented by sonication. Given the different degradation kinetics of mitochondrial DNA compared to nuclear DNA, the mitochondrially encoded 12s RNA and 18s ribosomal RNA genes were taken into consideration; fragments of various sizes were amplified in qPCR and the size distribution was thoroughly investigated. DNA damage degree was evaluated by calculating damage frequency (λ) and the ratio between the amount of the different fragments and that of the smallest fragment (Q). The results demonstrate that both indices were found to be suitable for identifying, among the samples tested, those less damaged and suitable for post-extraction analysis; mitochondrial DNA is more damaged than nuclear, in fact, amplicons up to 152 bp and 253 bp, respectively are obtained.

Tooth chipping patterns in Archaeolemur provide insight into diet and behavior

OBJECTIVES

Archaeolemur is a recently extinct genus of lemur that is often compared to some Cercopithecidae, especially baboons. This is due in part to their derived dentition, with large anterior teeth and reduced bilophodont molars. Research involving comparative morphology, analysis of coprolites, isotopes, and enamel structure, have suggested Archaeolemur had an omnivorous diet involving mechanically challenging items. Yet, microwear analysis of posterior teeth does not necessarily support this conclusion.

MATERIALS AND METHODS

In this macroscopic study, dental chipping was recorded on permanent teeth of Archaeolemur from different localities (53 individuals; 447 permanent teeth; including both A. edwardsi and A. majori specimens). This study aimed to compare chipping patterns across the dentition of Archaeolemur with chipping in other primates.

RESULTS

The results show enamel chipping was prevalent on the anterior teeth of Archaeolemur (38.9% of anterior teeth showed at least one fracture) yet rare in posterior teeth (9%). There was a decrease in chipping frequency across the dentition, moving distally from incisors (50%; 20/40), through caniniform teeth (30%; 15/50), premolars (9.5%; 16/169), and molars (8.5%; 16/188).

DISCUSSION

The results support previous research suggesting Archaeolemur had a varied omnivorous diet in which the anterior dentition was used for extensive food processing. This likely included mechanically challenging items such as tough/hard large fruits, small vertebrates, and crustaceans. Such a high rate of chipping in the anterior dentition is uncommon in other primates, with exception of hominins.

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.

A Fast and Efficient Single-stranded Genomic Library Preparation Method Optimized for Ancient DNA

We present a protocol to prepare extracted DNA for sequencing on the Illumina sequencing platform that has been optimized for ancient and degraded DNA. Our approach, the Santa Cruz Reaction or SCR, uses directional splinted ligation of Illumina’s P5 and P7 adapters to convert natively single-stranded DNA and heat denatured double-stranded DNA into sequencing libraries in a single enzymatic reaction. To demonstrate its efficacy in converting degraded DNA molecules, we prepare 5 ancient DNA extracts into sequencing libraries using the SCR and 2 of the most commonly used approaches for preparing degraded DNA for sequencing: BEST, which targets and converts double-stranded DNA, and ssDNA2.0, which targets and converts single-stranded DNA. We then compare the efficiency with which each approach recovers unique molecules, or library complexity, given a standard amount of DNA input. We find that the SCR consistently outperforms the BEST protocol in recovering unique molecules and, despite its relative simplicity to perform and low cost per library, has similar performance to ssDNA2.0 across a wide range of DNA inputs. The SCR is a cost- and time-efficient approach that minimizes the loss of unique molecules and makes accessible a taxonomically, geographically, and a temporally broader sample of preserved remains for genomic analysis.

Mandibular symphyseal fusion in fossil primates: Insights from correlated patterns of jaw shape and masticatory function in living primates

OBJECTIVES

Variation in primate masticatory form and function has been extensively researched through both morphological and experimental studies. As a result, symphyseal fusion in different primate clades has been linked to either the recruitment of vertically directed balancing-side muscle force, the timing and recruitment of transversely directed forces, or both. This study investigates the relationship between jaw muscle activity patterns and morphology in extant primates to make inferences about masticatory function in extinct primates, with implications for understanding the evolution of symphyseal fusion.

MATERIALS AND METHODS

Three-dimensional mandibular landmark data were collected for 31 extant primates and nine fossil anthropoids and subfossil lemur species. Published electromyography (EMG) data were available for nine of the extant primate species. Partial least squares analysis and phylogenetic partial least squares analysis were used to identify relationships between EMG and jaw shape data and evaluate variation in jaw morphology.

RESULTS

Primates with partial and complete symphyseal fusion exhibit shape-function patterns associated with the wishboning motor pattern and loading regime, in contrast to shape-function patterns of primates with unfused jaws. All fossil primates examined (except Apidium) exhibit jaw morphologies suggestive of the wishboning motor pattern demonstrated in living anthropoids and indriids.

DISCUSSION

Partial fusion in Catopithecus, similar to indriids and some subfossil lemurs, may be sufficient to resist, or transfer, some amounts of transversely directed balancing-side muscle force at the symphysis, representing a transition to greater reliance on transverse jaw movement during mastication. Furthermore, possible functional convergences in physiological patterns during chewing (i.e., Archaeolemur) are identified.

Evolution and allometry of calcaneal elongation in living and extinct primates

Specialized acrobatic leaping has been recognized as a key adaptive trait tied to the origin and subsequent radiation of euprimates based on its observed frequency in extant primates and inferred frequency in extinct early euprimates. Hypothesized skeletal correlates include elongated tarsal elements, which would be expected to aid leaping by allowing for increased rates and durations of propulsive acceleration at takeoff. Alternatively, authors of a recent study argued that pronounced distal calcaneal elongation of euprimates (compared to other mammalian taxa) was related primarily to specialized pedal grasping. Testing for correlations between calcaneal elongation and leaping versus grasping is complicated by body size differences and associated allometric affects. We re-assess allometric constraints on, and the functional significance of, calcaneal elongation using phylogenetic comparative methods, and present an evolutionary hypothesis for the evolution of calcaneal elongation in primates using a Bayesian approach to ancestral state reconstruction (ASR). Results show that among all primates, logged ratios of distal calcaneal length to total calcaneal length are inversely correlated with logged body mass proxies derived from the area of the calcaneal facet for the cuboid. Results from phylogenetic ANOVA on residuals from this allometric line suggest that deviations are explained by degree of leaping specialization in prosimians, but not anthropoids. Results from ASR suggest that non-allometric increases in calcaneal elongation began in the primate stem lineage and continued independently in haplorhines and strepsirrhines. Anthropoid and lorisid lineages show stasis and decreasing elongation, respectively. Initial increases in calcaneal elongation in primate evolution may be related to either development of hallucal-grasping or a combination of grasping and more specialized leaping behaviors. As has been previously suggested, subsequent increases in calcaneal elongation are likely adaptations for more effective acrobatic leaping, highlighting the importance of this behavior in early euprimate evolution.

The lemur revolution starts now: the genomic coming of age for a non-model organism

Morris Goodman was a revolutionary. Together with a mere handful of like-minded scientists, Morris established himself as a leader in the molecular phylogenetic revolution of the 1960s. The effects of this revolution are most evident in this journal, which he founded in 1992. Happily for lemur biologists, one of Morris Goodman's primary interests was in reconstructing the phylogeny of the primates, including the tooth-combed Lorisifomes of Africa and Asia, and the Lemuriformes of Madagascar (collectively referred to as the suborder Strepsirrhini). This paper traces the development of molecular phylogenetic and evolutionary genetic trends and methods over the 50-year expanse of Morris Goodman's career, particularly as they apply to our understanding of lemuriform phylogeny, biogeography, and biology. Notably, this perspective reveals that the lemuriform genome is sufficiently rich in phylogenetic signal such that the very earliest molecular phylogenetic studies - many of which were conducted by Goodman himself - have been validated by contemporary studies that have exploited advanced computational methods applied to phylogenomic scale data; studies that were beyond imagining in the earliest days of phylogeny reconstruction. Nonetheless, the frontier still beckons. New technologies for gathering and analyzing genomic data will allow investigators to build upon what can now be considered a nearly-known phylogeny of the Lemuriformes in order to ask innovative questions about the evolutionary mechanisms that generate and maintain the extraordinary breadth and depth of biological diversity within this remarkable clade of primates.

An alu-based phylogeny of lemurs (infraorder: Lemuriformes)

LEMURS (INFRAORDER: Lemuriformes) are a radiation of strepsirrhine primates endemic to the island of Madagascar. As of 2012, 101 lemur species, divided among five families, have been described. Genetic and morphological evidence indicates all species are descended from a common ancestor that arrived in Madagascar ∼55-60 million years ago (mya). Phylogenetic relationships in this species-rich infraorder have been the subject of debate. Here we use Alu elements, a family of primate-specific Short INterspersed Elements (SINEs), to construct a phylogeny of infraorder Lemuriformes. Alu elements are particularly useful SINEs for the purpose of phylogeny reconstruction because they are identical by descent and confounding events between loci are easily resolved by sequencing. The genome of the grey mouse lemur (Microcebus murinus) was computationally assayed for synapomorphic Alu elements. Those that were identified as Lemuriformes-specific were analyzed against other available primate genomes for orthologous sequence in which to design primers for PCR (polymerase chain reaction) verification. A primate phylogenetic panel of 24 species, including 22 lemur species from all five families, was examined for the presence/absence of 138 Alu elements via PCR to establish relationships among species. Of these, 111 were phylogenetically informative. A phylogenetic tree was generated based on the results of this analysis. We demonstrate strong support for the monophyly of Lemuriformes to the exclusion of other primates, with Daubentoniidae, the aye-aye, as the basal lineage within the infraorder. Our results also suggest Lepilemuridae as a sister lineage to Cheirogaleidae, and Indriidae as sister to Lemuridae. Among the Cheirogaleidae, we show strong support for Microcebus and Mirza as sister genera, with Cheirogaleus the sister lineage to both. Our results also support the monophyly of the Lemuridae. Within Lemuridae we place Lemur and Hapalemur together to the exclusion of Eulemur and Varecia, with Varecia the sister lineage to the other three genera.

Ancient DNA studies: new perspectives on old samples

In spite of past controversies, the field of ancient DNA is now a reliable research area due to recent methodological improvements. A series of recent large-scale studies have revealed the true potential of ancient DNA samples to study the processes of evolution and to test models and assumptions commonly used to reconstruct patterns of evolution and to analyze population genetics and palaeoecological changes. Recent advances in DNA technologies, such as next-generation sequencing make it possible to recover DNA information from archaeological and paleontological remains allowing us to go back in time and study the genetic relationships between extinct organisms and their contemporary relatives. With the next-generation sequencing methodologies, DNA sequences can be retrieved even from samples (for example human remains) for which the technical pitfalls of classical methodologies required stringent criteria to guaranty the reliability of the results. In this paper, we review the methodologies applied to ancient DNA analysis and the perspectives that next-generation sequencing applications provide in this field.

A dig into the past mitochondrial diversity of Corsican goats reveals the influence of secular herding practices

The goat (Capra hircus) is one of the earliest domesticated species ca. 10,500 years ago in the Middle-East where its wild ancestor, the bezoar (Capra aegagrus), still occurs. During the Neolithic dispersal, the domestic goat was then introduced in Europe, including the main Mediterranean islands. Islands are interesting models as they maintain traces of ancient colonization, historical exchanges or of peculiar systems of husbandry. Here, we compare the mitochondrial genetic diversity of both medieval and extant goats in the Island of Corsica that presents an original and ancient model of breeding with free-ranging animals. We amplified a fragment of the Control Region for 21 medieval and 28 current goats. Most of them belonged to the A haplogroup, the most worldwide spread and frequent today, but the C haplogroup is also detected at low frequency in the current population. Present Corsican goats appeared more similar to medieval goats than to other European goat populations. Moreover, 16 out of the 26 haplotypes observed were endemic to Corsica and the inferred demographic history suggests that the population has remained constant since the Middle Ages. Implications of these results on management and conservation of endangered Corsican goats currently decimated by a disease are addressed.

"Life history space": a multivariate analysis of life history variation in extant and extinct Malagasy lemurs

Studies of primate life history variation are constrained by the fact that all large-bodied extant primates are haplorhines. However, large-bodied strepsirrhines recently existed. If we can extract life history information from their skeletons, these species can contribute to our understanding of primate life history variation. This is particularly important in light of new critiques of the classic "fast-slow continuum" as a descriptor of variation in life history profiles across mammals in general. We use established dental histological methods to estimate gestation length and age at weaning for five extinct lemur species. On the basis of these estimates, we reconstruct minimum interbirth intervals and maximum reproductive rates. We utilize principal components analysis to create a multivariate "life history space" that captures the relationships among reproductive parameters and brain and body size in extinct and extant lemurs. Our data show that, whereas large-bodied extinct lemurs can be described as "slow" in some fashion, they also varied greatly in their life history profiles. Those with relatively large brains also weaned their offspring late and had long interbirth intervals. These were not the largest of extinct lemurs. Thus, we distinguish size-related life history variation from variation that linked more strongly to ecological factors. Because all lemur species larger than 10 kg, regardless of life history profile, succumbed to extinction after humans arrived in Madagascar, we argue that large body size increased the probability of extinction independently of reproductive rate. We also provide some evidence that, among lemurs, brain size predicts reproductive rate better than body size.

Stable carbon and nitrogen isotope enrichment in primate tissues

Isotopic studies of wild primates have used a wide range of tissues to infer diet and model the foraging ecologies of extinct species. The use of mismatched tissues for such comparisons can be problematic because differences in amino acid compositions can lead to small isotopic differences between tissues. Additionally, physiological and dietary differences among primate species could lead to variable offsets between apatite carbonate and collagen. To improve our understanding of the isotopic chemistry of primates, we explored the apparent enrichment (ε*) between bone collagen and muscle, collagen and fur or hair keratin, muscle and keratin, and collagen and bone carbonate across the primate order. We found that the mean ε* values of proteinaceous tissues were small (≤1‰), and uncorrelated with body size or phylogenetic relatedness. Additionally, ε* values did not vary by habitat, sex, age, or manner of death. The mean ε* value between bone carbonate and collagen (5.6 ± 1.2‰) was consistent with values reported for omnivorous mammals consuming monoisotopic diets. These primate-specific apparent enrichment values will be a valuable tool for cross-species comparisons. Additionally, they will facilitate dietary comparisons between living and fossil primates.

The climatic niche diversity of malagasy primates: a phylogenetic perspective

BACKGROUND

Numerous researchers have posited that there should be a strong negative relationship between the evolutionary distance among species and their ecological similarity. Alternative evidence suggests that members of adaptive radiations should display no relationship between divergence time and ecological similarity because rapid evolution results in near-simultaneous speciation early in the clade's history. In this paper, we performed the first investigation of ecological diversity in a phylogenetic context using a mammalian adaptive radiation, the Malagasy primates.

METHODOLOGY/PRINCIPAL FINDINGS

We collected data for 43 extant species including: 1) 1064 species by locality samples, 2) GIS climate data for each sampling locality, and 3) the phylogenetic relationships of the species. We calculated the niche space of each species by summarizing the climatic variation at localities of known occurrence. Climate data from all species occurrences at all sites were entered into a principal components analysis. We calculated the mean value of the first two PCA axes, representing rainfall and temperature diversity, for each species. We calculated the K statistic using the Physig program for Matlab to examine how well the climatic niche space of species was correlated with phylogeny.

CONCLUSIONS/SIGNIFICANCE

We found that there was little relationship between the phylogenetic distance of Malagasy primates and their rainfall and temperature niche space, i.e., closely related species tend to occupy different climatic niches. Furthermore, several species from different genera converged on a similar climatic niche. These results have important implications for the evolution of ecological diversity, and the long-term survival of these endangered species.

Patterns of macroevolution among Primates inferred from a supermatrix of mitochondrial and nuclear DNA

Here, we present a new primate phylogeny inferred from molecular supermatrix analyses of size 42 kb containing 70% of missing data, and representing 75% of primate species diversity. The supermatrix was analysed using a gene-partitioned maximum likelihood approach to obtain an exhaustive molecular phylogenetic framework. All clades recovered from recent molecular works were upheld in our analyses demonstrating that the presence of missing data did not bias our supermatrix inference. The resulting phylogenetic tree was subsequently dated with a molecular dating method to provide a timescale for speciation events. Results obtained from our relaxed molecular clock analyses concurred with previous works based on the same fossil constraints. The resulting dated tree allowed to infer of macroevolutionary processes among the primates. Shifts in diversification rate and speciation rates were determined using the SymmeTREE method and a birthdeath process. No significant asymmetry was detected for the primate clade, but significant shifts in diversification rate were identified for seven clades: Anthropoidea, Lemuriformes, Lemuridae, Galagidae, Callithrix genus, the Cercopithecinae and Asian Macaca. Comparisons with previous primate supertree results reveal that (i) there was a diversification event at the root of the Lemuriformes, (ii) a higher diversification rate is detected for Cercopithecidae and Anthropoidea and (iii) a shift in diversification is always recovered for Macaca genus. Macroevolutionary inferences and primate divergence dates show that major primate diversification events occurred after the Paleogene, suggesting the extinction of ancient primate lineages.