Mitochondrial evidence for multiple radiations in the evolutionary history of small apes
© Thinh et al; licensee BioMed Central Ltd. 2010
Received: 20 October 2009
Accepted: 12 March 2010
Published: 12 March 2010
Gibbons or small apes inhabit tropical and subtropical rain forests in Southeast Asia and adjacent regions, and are, next to great apes, our closest living relatives. With up to 16 species, gibbons form the most diverse group of living hominoids, but the number of taxa, their phylogenetic relationships and their phylogeography is controversial. To further the discussion of these issues we analyzed the complete mitochondrial cytochrome b gene from 85 individuals representing all gibbon species, including most subspecies.
Based on phylogenetic tree reconstructions, several monophyletic clades were detected, corresponding to genera, species and subspecies. A significantly supported branching pattern was obtained for members of the genus Nomascus but not for the genus Hylobates. The phylogenetic relationships among the four genera were also not well resolved. Nevertheless, the new data permitted the estimation of divergence ages for all taxa for the first time and showed that most lineages emerged during four short time periods. In the first, between ~6.7 and ~8.3 mya, the four gibbon genera diverged from each other. In the second (~3.0 - ~3.9 mya) and in the third period (~1.3 - ~1.8 mya), Hylobates and Hoolock differentiated. Finally, between ~0.5 and ~1.1 mya, Hylobates lar diverged into subspecies. In contrast, differentiation of Nomascus into species and subspecies was a continuous and prolonged process lasting from ~4.2 until ~0.4 mya.
Although relationships among gibbon taxa on various levels remain unresolved, the present study provides a more complete view of the evolutionary and biogeographic history of the hylobatid family, and a more solid genetic basis for the taxonomic classification of the surviving taxa. We also show that mtDNA constitutes a useful marker for the accurate identification of individual gibbons, a tool which is urgently required to locate hunting hotspots and select individuals for captive breeding programs. Further studies including nuclear sequence data are necessary to completely understand the phylogeny and phylogeography of gibbons.
Although in several aspects unique among primates and with up to 16 species the most diverse group of apes, gibbons are still in the shadow of great apes in respect of scientific studies, conservation efforts and public awareness. However, many gibbon species are on the brink of extinction and most of them are classified as "Endangered" or even "Critically Endangered" . With approximately 20 individuals left in its native habitat, the Hainan gibbon (Nomascus hainanus) is the rarest primate in the world [6, 13, 14]. Responsible for this critical situation is habitat loss and hunting, which both have seriously reduced gibbon populations throughout their range [15, 16]. Hence, much more attention has to be drawn on the gibbons' situation and extensive conservation actions are urgently required to save them from extinction .
While gibbons are widely considered to form a monophyletic clade, there is no consensus about the phylogeny and taxonomy within the family. Although various studies based on morphology, behaviour, vocalisation, protein electrophoresis, karyotyping and DNA sequencing were conducted [3–5, 7, 17–35], neither a congruent phylogeny nor a consistent taxonomic classification was obtained. Moreover, incomplete taxon sampling as well as misidentified specimens resulted in only fragmentary or even false conclusions. Accordingly, the classification of gibbon taxa at various taxonomic levels as well as their phylogenetic relationships remain disputed and a consensus is far from being available.
Common names, IUCN classification and proposed classification of gibbons.
IUCN classification 
Kloss' s gibbon
Eastern Müller's Bornean gibbon
Hylobates muelleri muelleri
Northern Müller's Bornean gibbon
Hylobates muelleri funereus
Abbott's Müller's Bornean gibbon
Hylobates muelleri abbotti
Bornean white-bearded gibbon
Malayan lar gibbon
Hylobates lar lar
Hylobates lar lar*
Sumatran lar gibbon
Hylobates lar vestitus
Hylobates lar vestitus*
Mainland lar gibbon
Hylobates lar entelloides
Hylobates lar entelloides*
Carpenter's lar gibbon
Hylobates lar carpenteri
Hylobates lar carpenteri*
Yunnan lar gibbon
Hylobates lar yunnanensis
Hylobates lar yunnanensis*
Silvery Javan gibbon
Western hoolock gibbon
Eastern hoolock gibbon
Cao-vit crested gibbon
Black crested gibbon
Nomascus concolor concolor
Nomascus concolor concolor*
West Yunnan black crested gibbon
Nomascus concolor furvogaster
Nomascus concolor concolor*
Central Yunnan black crested gibbon
Nomascus concolor jingdongensis
Nomascus concolor concolor*
Laotian black crested gibbon
Nomascus concolor lu
Nomascus concolor lu*
Northern white-cheeked gibbon
Southern white-cheeked gibbon
16 species, 12 subspecies
18 species, 7 subspecies
In the present study, we analyse the complete mitochondrial cytochrome b (cytb) gene from 85 individuals, which represent all gibbon genera and species, and most subspecies. Based on our data, we are able to 1) provide the most complete phylogeny of gibbons on all taxonomic levels, 2) estimate divergence times between lineages, 3) establish a reliable classification, 4) elucidate gibbon phylogeography, and 5) provide a tool for the species identification of gibbon individuals.
From all 85 gibbons, we successfully generated sequences of the complete mitochondrial cytb gene (1,140 bp). A contamination of our dataset with nuclear pseudogenes (numts) can be regarded as minimal, because no multiple amplifications of different copies were detected by direct sequencing. All sequences were correctly transcribed, and identical sequences were obtained for the same individual in cases where different material types were available. Moreover, no inconsistent positions were detected in alignments, which were assembled from overlapping sequences. Cross-contamination between individuals can be excluded as well, since all negative controls revealed no amplifications and randomly repeated PCRs for the same individual produced identical sequences.
Among the 85 individual gibbons studied, no identical haplotypes were detected. The cytb alignment comprising solely gibbons was characterized by 429 variable sites, of which 374 were parsimony-informative. In the complete alignment, which additionally contained great ape, human and hamadryas baboon representatives, we observed 565 variable sites, of which 462 were parsimony-informative.
Support values and Bayesian divergence date estimates (in mya)*.
Mean (95% CI)
Papio - Hominoidea
Hylobatidae - Hominidae
Pongo - Gorilla/Pan/Homo
Gorilla - Pan/Homo
Pan - Homo
Pan troglodytes - P. paniscus
Pongo pygmaeus - P. abelii
Nomascus - Symphalangus/Hoolock/Hylobates
Symphalangus - Hoolock/Hylobates
Hoolock - Hylobates
Hylobates klossii - H. pileatus/H. moloch/H. agilis/H. albibarbis/H. lar/H. muelleri
H. pileatus/H. moloch - H. agilis/H. albibarbis/H. lar/H. muelleri
H. muelleri - H. agilis/H. albibarbis/H. lar
H. agilis/H. albibarbis - H. lar
H. agilis - H. albibarbis
H. pileatus - H. moloch
H. muelleri funereus - H. m. abbotti/H. m. muelleri
H. muelleri abbotti - H. m. muelleri
H. agilis agilis - H. a. unko
H. lar vestitus - H. l. lar/H. l. entelloides/H. l. carpenteri/H. l. yunnanensis
H. l. lar - H. entelloides/H. l. carpenteri/H. l. yunnanensis
H. l. entelloides - H. l. carpenteri/H. l. yunnanensis
H. l. carpenteri - H. l. yunnanensis
MRCA H. klossii
MRCA H. muelleri muelleri
MRCA H. albibarbis
MRCA H. agilis unko
MRCA H. agilis agilis
MRCA H. lar carpenteri
MRCA H. lar entelloides
MRCA H. pileatus
MRCA H. moloch
Hoolock hoolock - H. leuconedys
MRCA H. leuconedys
MRCA H. hoolock
MRCA Symphalangus syndactylus
Nomascus hainanus/N. nasutus - N. concolor/N. gabriellae/N. leucogenys/N. siki
N. hainanus - N. nasutus
N. concolor - N. gabriellae/N. leucogenys/N. siki
N. gabriellae - N. leucogenys/N. siki
N. leucogenys - N. siki
N. concolor lu - N. c. concolor/N. c. furvogaster/N. c. jingdongensis
MRCA N. nasutus
MRCA N. concolor lu
MRCA N. concolor concolor/N. c. furvogaster/N. jingdongensis
MRCA N. gabriellae
MRCA N. leucogenys
MRCA N. siki
Based on divergence age estimates, gibbons separated from great apes and humans 16.26 million years ago (mya) (for 95% credibility intervals see Table 2). Within hominids, Pongo branched off first (13.83 mya), followed by Gorilla (8.90 mya), before finally Homo and Pan diverged from each other (6.56 mya). The differentiation of Pongo and Pan into species occurred 4.12 and 2.74 mya, respectively. In an initial radiation, gibbons diverged within a relative short time period of only 1.65 million years (6.69-8.34 mya) into four genera. Within Hylobates, most species diverged from each other between 3.02 and 3.90 mya. The only exception was the separation of H. albibarbis from H. agilis 1.56 mya, which was in the time frame of subspecies splits within H. muelleri (1.42-1.78 mya). Differentiation of H. lar into subspecies occurred even later (0.52-1.05 mya). The two Hoolock species diverged 1.42 mya from each other. In Nomascus, differentiation into species took place over a longer time period, lasting from 4.24 until 0.55 mya. The most recent species divergence within Nomascus occurred between N. siki and N. leucogenys (0.55 mya), which was in a similar range as the separation of N. concolor lu from the other N. concolor subspecies (0.43 mya).
By analysing all species and most subspecies, the present study provides the most complete view into the evolutionary history of the gibbon family. However, as in earlier molecular studies on gibbons [26–35], relationships on various taxonomic levels are less resolved and partially contradict earlier findings. While the herein depicted branching pattern among genera is identical with that found in earlier studies using also cytb  or D-loop  sequences, it differs from another cytb-based study  in placing Nomascus and not Symphalangus as most basal genus. Studies based on mitochondrial ND3-ND4 sequences  or chromosomal rearrangements  suggest Hoolock as most ancestral lineage, and Nomascus together with either Hylobates  or Symphalangus  as the most recently diverged genera. For Hylobates, our data indicate a basal position of H. klossii, and a further division into a clade consisting of H. lar, H. muelleri, H. agilis and H. albibarbis, and another one with H. moloch and H. pileatus. Various branching patterns among Hylobates species are proposed [27, 31, 32, 35], which all differ from our one, but respective support values are similarly low as in our study. In contrast, the relationships found among species of the genus Nomascus are well resolved and identical with that suggested by [30, 31, 33, 34].
According to our and earlier data, relationships among gibbon genera and Hylobates species remain disputed, which most likely can be explained by the separation of respective lineages within relative short time periods. This becomes even more obvious when considering estimated divergence ages, which fall into four temporal windows. In the first, between ~6.7 and ~8.3 mya, the four gibbon genera originated. In a second radiation, between ~3.0 and ~3.9 mya, Hylobates split into various species, and in a third burst, between ~1.3 and ~1.8 mya, H. muelleri, the H. agilis + H. albibarbis clade and Hoolock further differentiated. Finally, in a fourth radiation, between ~0.5 and ~1.1 mya, H. lar diverged into subspecies. In contrast, speciation in Nomascus was a continuous process, lasting from 4.24 until 0.55 mya.
Our data show that mitochondrial DNA (mtDNA) provides a powerful tool for the identification and taxonomic classification of gibbons, because taxa form strongly supported monophyletic clades, or at least appear to form distinct lineages in those cases where only one individual per taxon was tested. Moreover, most differentiation events fall into four temporal periods, which allow a hierarchical ranking as proposed by Goodman et al. , though the threshold for the recognition of a certain taxonomic unit whether genus, species, or subspecies remains disputed. Hence, to provide a more reliable classification, we compare divergence ages among gibbon lineages with those among other Asian primates and hominids.
Accordingly and concordant with recent classifications [4, 12, 16, 29, 34, 37, 38, 41], the four major gibbon lineages are proposed as distinct genera (Table 1), since they split from each other in a similar time range as did colobine genera [[42, 43], Roos C, Zinner D, Schwarz C, Nash SD, Xing J, Batzer MA, Leendertz FH, Ziegler T, Perwitasari-Farajallah D, Nadler T, Walter L, Osterholz M: Nuclear versus mitochondrial DNA: evidence for hybridization in colobine monkeys, submitted] or African great apes and human [40, 42]. Most species of Hylobates and Nomascus emerged in or around the second radiation, which is on the same time scale as species splits within Pongo and Pan, and the separation of species groups within Macaca [44, 45] and Trachypithecus . Thus, taxa originating in this time period should be recognized as distinct species (H. moloch, H. pileatus, H. klossii, H. lar, H. muelleri, H. agilis/H. albibarbis, H. hoolock/H. leuconedys, N. nasutus, N. hainanus, N. concolor, N. gabriellae/N. leucogenys/N. siki), and might be even classified as species groups. Further differentiation events among gibbons occurred in the third time period, which is in a similar window as several speciation events within macaques [44, 45]. Accordingly, H. leuconedys and H. albibarbis should be separated from H. hoolock and H. agilis on species level, respectively, and the three subspecies of H. muelleri could be considered for elevation to species level. Moreover, H. agilis is divided into two clades, which refer to individuals identified by pelage coloration as H. agilis agilis and H. agilis unko. However, in a recent work based on a larger number of individuals a reciprocal monophyly of both lineages is doubted , and, hence, we provisionally recognize H. agilis as monotypic. For H. lar, only a few unambiguously identified specimens were available for our study, but these represent at least four of the five recognized subspecies, while the identity of the putative H. lar yunnanensis individual remains uncertain. Based on our data, H. lar subspecies form distinct lineages, which diverged relative recently. We provisionally accept all five subspecies, though ongoing studies might reject some or all of them. For N. concolor, our data indicate a separation of N. concolor lu from the remaining subspecies, which form a clade without further subdivision into taxa. Hence and concordant with Monda et al.  and Roos et al. , we provisionally classify N. concolor furvogaster and N. concolor jingdongensis as synonyms of N. concolor concolor, while we feel N. concolor lu is a separate subspecies. We further separate N. gabriellae from N. siki/N. leucogenys on species level, while it is questionable whether the latter two should be recognized as species or subspecies. Our study reveals a split between both taxa just 0.55 mya, which is in a similar range as the subspecies differentiation within H. lar or N. concolor. Hence, a separation of both taxa only on subspecies level would be indicated. However, both taxa show slight differences in vocalisation and facial colouration [4, 5, 15], and Carbone et al.  found a chromosomal inversion unique to N. leucogenys. Accordingly, we follow here the current view and recognize N. leucogenys and N. siki as distinct species. In summary, we recognize four gibbon genera with 18 species and seven subspecies (Table 1).
Multiple radiations in the evolutionary history of gibbons suggest a complicated biogeographic pattern leading to the current distribution of gibbon taxa. Since gibbons are arboreal [7, 39], radiations most likely were correlated with expanding forest habitats. In fact, the complete range of gibbons experienced complex geographical and environmental changes during the last ten million years. Notably, in the late Miocene as well as in the Plio- and Pleistocene, a series of dramatic climatic changes influenced the geography and vegetation in the region, leading to shifts in the extension and distribution of different habitat types [49–54]. In particular, periods of maximum glaciation might have reduced rainforest cover, resulting in the appearance of more open and deciduous vegetation types in many parts of the region [[52–57], but see ]. Moreover, due to the alternately falling and rising sea water levels during the several glacial and interglacial periods [59–64], connections and separations of landmasses were common, and repeated migration between islands and today's mainland was possible [65–67].
By combining the available information, we develop the following dispersal scenario for gibbons, which is in general agreement with that proposed by Chatterjee [32, 68], Harrison et al. , and Jablonski and Chaplin , but which differs substantially from them in some aspects. Accordingly, gibbons most likely originated on the Asian mainland, because all four gibbon genera occur there. Specifically, the Hengduan mountains in the border region of today's Burma, India and China might have been a possible diversification hotspot [71, 72]. In the region, all the larger Southeast Asian rivers (Mekong, Salween, Yangtze) rise, which are all well known as barriers for arboreal primates . Although these rivers changed their courses several times, their upper reaches in the Hengduan mountains exist at least since the early Miocene . Recently, the Hengduan mountains were also proposed as a region of differentiation for colobine monkeys, and, most interestingly, respective splitting events occurred on a similar time scale as in gibbons [Roos C, Zinner D, Schwarz C, Nash SD, Xing J, Batzer MA, Leendertz FH, Ziegler T, Perwitasari-Farajallah D, Nadler T, Walter L, Osterholz M: Nuclear versus mitochondrial DNA: evidence for hybridization in colobine monkeys, submitted]. In fact, in the late Miocene, widely distributed rain forest habitats promoted range extension for arboreal primates [50, 54]. Accordingly, in the late Miocene, Nomascus invaded the region east of the Mekong, Hoolock entered the region west of the Salween, and Hylobates and Symphalangus migrated into the area in-between and later on into Sundaland.
Hylobates successfully colonized large parts of Sundaland, but also survived on the Asian mainland. Shortly after its arrival in Sundaland in the Pliocene, populations on the Asian mainland, the Malay peninsula, Sumatra, Borneo, Java and the Mentawai archipelago became isolated. At the same time, various species groups of the genera Macaca and Trachypithecus diverged [44–46], indicating dramatic environmental changes. In fact, this time period was characterized by global warming and sea levels similar to today [54, 61–63], which prevented migration between landmasses and, thus, promoted speciation due to vicariance. Whether Symphalangus experienced a similar range expansion in Sundaland like Hylobates, remains questionable. Today the genus appears only on Sumatra and the Malay peninsula, and fossil data provide only evidence for its historical occurrence on Java and Sumatra . In the early Pleistocene, further differentiation in Hylobates occurred on Borneo and Sumatra, and in Hoolock on the mainland which is on a similar time scale when macaque species diverged [44, 45], and which might has been triggered by the shrinking of forest habitats due to cold phases [, but see ]. Notably, H. albibarbis is mitochondrially closer related to Sumatran H. agilis than to the other Bornean gibbons, and acoustic, morphological and chromosomal data suggest an intermediate position [2, 5, 47, 75]. Accordingly, H. albibarbis might be the product of an ancient hybridization event, in which proto-H. agilis invaded Borneo during sea level lowstands [61–64], and successfully reproduced with proto-H. muelleri. As we find mtDNA of proto-H. agilis in H. albibarbis, female introgression is the most likely hybridization scenario, which is in agreement with recent findings, that gibbon females disperse over longer distances than males . Finally, in a last range expansion in the early to middle Pleistocene, H. lar colonized, starting from its Sumatran refuge, the Malaysian peninsular and mainland Southeast Asia [see also ].
In contrast to the biogeographic pattern found in Hylobates and to the scenario proposed by Chatterjee [32, 68], for Nomascus not a radiation but a successive migration from North to South over a long time period becomes evident. Based on our data, Nomascus originated in the border region of Vietnam and China in the early Pliocene and it took to the early Pleistocene until the genus reached the southern extend of its current distribution in southern Vietnam and Cambodia.
All gibbon species are on the brink of extinction and, with the exception of H. leuconedys (Vulnerable), are classified as "Endangered" or even "Critically Endangered" [12, 16]. With approximately 20 individuals left in its native habitat, the Hainan gibbon (N. hainanus) is the rarest primate in the world [6, 13, 14], and the situation for its closest relative, the Cao-vit crested gibbon (N. nasutus) with approximately 100 individuals left [12, 77], as well as for other gibbon species, the situation is alarming. Reasons for the decline of gibbons are manifold, but habitat loss due to forest clearance for agricultural use, oil palm or rubber plantations, gold mining, or charcoal and timber production, as well as illegal hunting for food and sport, and the trade for pets or medicine are major threats to wild gibbon populations [15, 16].
To save gibbons from extinction, urgent actions are required to prevent ongoing habitat destruction and hunting, and to build up a viable gene pool in captivity for later release purposes. Specifically, to prevent or at least reduce hunting, hunting hotspots have to be identified. Therefore, it is crucial to confirm the taxon identity and if possible the geographical origin of confiscated gibbons or their remains. Similarly, to avoid artificial hybrids, only gibbons with clear taxon identity should be considered for reproduction in zoos or rescue centres. Finally, if captive gibbons are reintroduced into the wild, it has to be ascertained that these gibbons are pure individuals and of the same taxon as those, which naturally occur in the area they are to be released.
An accurate taxonomic identification of gibbons based on vocal data or pelage colouration is sometimes complicated [4, 5]. In this respect, mtDNA analysis might be a promising tool. As shown in our study, gibbon taxa can be diagnosed through mtDNA, and, hence, a secure identification can easily be obtained. Yet since mtDNA is only maternally inherited, possible hybrids will not be detected in such analysis, so that additional markers should be studied as well.
Due to a nearly complete taxon sampling, the present study provides the most comprehensive insights into the evolutionary and biogeographic history of the hylobatid family. Based on estimated divergence ages and unresolved relationships among gibbon taxa on various levels, several radiation-like splitting events are indicated, which suggest a complex biogeographic history of gibbons. Presumably, most of these differentiation events occurred in wave-like range expansions in Sundaland and the Asian mainland followed by vicariance effects, most likely caused by alternately shrinking and expanding rain forest habitats and by repeated separations and connections of landmasses. In contrast, in the region east of the Mekong river gibbons underwent a successive North-to-South migration. Our study also shows that mtDNA provides a solid platform for the taxonomic classification of gibbons and that mtDNA can be successfully applied to accurately identify the species affiliation of gibbon individuals, which is urgently required for conservation purposes. However, to completely understand the phylogeny and phylogeography of gibbons, to identify hybrids in captivity, or to trace possible ancient hybridization events as it might be indicated for H. albibarbis, further studies including extended mitochondrial as well as autosomal, X and Y chromosomal sequence data, are necessary.
A total of 85 specimens representing all species and most subspecies of hylobatids were included in our study. Blood, tissue, faecal or hair samples were collected during field surveys, in zoos or rescue centres, or from museum specimens between 1995 and 2008 (Additional file 1). Blood and hair samples were taken during routine health checks by veterinarians. Tissue samples were obtained only from deceased animals. Taxon identity of individuals was ascertained by pelage coloration, morphology and if possible by vocalization and geographic origin. With the exception of some H. lar individuals for which subspecies identity could not be traced, only clearly identified specimens were included in our study. Fresh tissue or faecal samples were preserved in 80-90% ethanol and dry samples (tissue, museum skins, hair samples) were placed in plastic bags without any additive. Samples were stored at ambient temperature for up to six months before further processing.
Total genomic DNA was extracted with the DNeasy Blood & Tissue and QIAamp DNA Stool Mini kits from Qiagen. When hair follicle cells were used, up to three hairs were directly implemented into the PCR reaction. From high-quality DNA, the complete mitochondrial cytb gene was PCR-amplified in a single fragment with the primers 5'-AATGATATGAAAAACCATCGTTGTA-3' and 5'-TTCATTTCCGGCTTACAAGAC-3'. For low-quality DNA, extracted from faeces or museums material, two to seven overlapping PCR products were amplified with primers constructed on the basis of sequences from conspecifics (respective primers are available from the authors upon request). For all amplifications, wax-mediated hot-start PCRs were performed for 40 cycles, each with a denaturation step at 92°C for 1 min, annealing at 60°C for 1 min, and extension at 72°C for 0.5-1.5 min, followed by a final extension step at 72°C for 5 min. The results of the PCR amplifications were checked on 1% agarose gels. Subsequently, PCR products were cleaned with the Qiagen Gel Extraction kit and sequenced on an ABI 3130xl sequencer using the BigDye Cycle Sequencing kit. Sequences were assembled with Geneious v4.6.1  and checked for their potential to be correctly transcribed. Gibbon haplotypes were deposited at GenBank and are available under the accession numbers GU321245-GU321329 (see also Additional file 1).
To prevent cross-species contaminations, laboratory procedures followed described standards . To exclude contaminations of the dataset with numts, we mainly used material in which nuclear DNA is highly degraded (faeces, museum tissue) [79, 80]. Moreover, the applied primers are known to amplify solely the mitochondrial copy of the gene in hylobatids , and for cross-validation purposes, for some specimens, sequences were generated using different material types (blood, faeces).
For phylogenetic reconstructions, we expanded our dataset with orthologous sequences from various hominids (Homo, Pan, Gorilla, Pongo) and Papio hamadryas, which was used as outgroup. Phylogenetic trees were constructed with MP and NJ algorithms as implemented in PAUP v4.0b10  as well as with ML and Bayesian algorithms, using the programs GARLI v0.951  and MrBayes v3.1.2 [83, 84]. For MP analysis, all characters were treated as unordered and equally weighted throughout. A heuristic search was performed with the maximum number of trees set to 100. For NJ and ML reconstructions, the optimal nucleotide substitution model (GTR + Γ) was chosen using Akaike information criterion (AIC) as implemented in MODELTEST v3.7 . Relative support of internal nodes was performed by bootstrap analyses with 10,000 (MP, NJ) or 500 replications (ML). In GARLI, only the model specification settings were adjusted according to the dataset, while all other settings were left at their default value. ML majority-rule consensus trees were calculated in PAUP. For Bayesian reconstructions, the dataset was partitioned into codon positions, each with its own substitution model. We used four Markov Chain Monte Carlo (MCMC) chains with the default temperature of 0.1. Four repetitions were run for 10,000,000 generations with tree and parameter sampling occurring every 100 generations. The first 25% of samples were discarded as burnin, leaving 75,001 trees per run. Posterior probabilities for each split and a phylogram with mean branch lengths were calculated from the posterior density of trees.
To estimate divergence times, a Bayesian MCMC method, which employs a relaxed molecular clock approach , as implemented in BEAST v1.4.8 , was used. A relaxed lognormal model of lineage variation and a Yule prior for branching rates was assumed. The alignment was partitioned into codon positions, and the substitution model, rate heterogeneity and base frequencies were unlinked across codon positions. Optimal nucleotide substitution models were chosen using AIC in MODELTEST.
For calibrations we used the fossil-based divergence between Homo and Pan, which was dated at 6 - 7 mya [88–90], the separation of Pongo from the Homo/Pan lineage ~14 mya , and the divergence of hominoids and cercopithecoids ~23 mya [92, 93]. Instead of hardbounded calibration points, we used the published dates as a normal distribution prior for the respective node. For the Homo - Pan divergence, this translates into a normal distribution with a mean of 6.5 mya and a standard deviation (SD) of 0.5 mya, for the separation of Pongo from the Homo/Pan clade into a mean of 14.0 mya and a SD of 1.0 mya, and for the hominoid - cercopithecoid divergence into a mean of 23 mya and a SD of 2 mya.
Since the estimation of phylogenetic relationships was not the main aim of this analysis, for the calculation an a-priori fixed tree topology as obtained from NJ reconstructions using the GTR + Γ model (Figure 2) was implemented. Four replicates were run for 10,000,000 generations with tree and parameter sampling occurring every 100 generations. The adequacy of a 10% burnin and convergence of all parameters were assessed by visual inspection of the trace of the parameters across generations using TRACER v1.4.1 . Subsequently, the sampling distributions were combined (25% burnin) using the software LogCombiner v1.4.8, and a consensus chronogram with node height distribution was generated and visualized with TreeAnnotator v1.4.8 and FigTree v1.2.2 .
We are grateful to the following colleagues, zoos, and institutions for providing permits or valuable gibbon materials: Claudia Barelli, Gareth Goldthorpe, Andrew Kitchener, Nicolas Lormée, Annette Schrod, Chris Smeenk, Berlin Zoo, Duisburg Zoo, Leipzig Zoo, Munich Zoo, Nuremberg Zoo, Rostock Zoo, Schwerin Zoo, Wuppertal Zoo, Zurich Zoo, Besancon Zoo, Plock Zoo, Banham Zoo, Bristol Zoo, Howletts Wild Animal Park, Paignton Zoo, Twycross Zoo, Beijing Zoo, Dhaka Zoo, Jakarta Zoo, Taman Safari, Singapore Zoo, Perth Zoo, Louisiana Purchase Gardens and Zoo, Bogor Agricultural University, the Indonesian Institute for Science, Museum für Naturkunde Berlin, National Museums Scotland Edinburgh, Natural History Museum Leiden, National Museum of Natural History Washington, Institute of Zoology of the Chinese Academy of Sciences Beijing, Bawangling National Nature Reserve, Bokeo Nature Reserve, Cat Tien National Park, Phong Nha-Ke Bang National Park and Khao Yai National Park. No international or national rules and regulations have been violated during sampling, and shipping. Many thanks also to Christiane Schwarz for her excellent laboratory work, and to Colin Groves and two anonymous reviewers for valuable comments on an earlier version of the manuscript. This study was financially supported by the German Primate Center, and the Biodiversitäts-Pakt of the Wissenschaftsgemeinschaft Gottfried-Wilhelm Leibniz.
- Fleagle JG: Primate Adaptations and Evolution. 1999, San Diego: Academic PressGoogle Scholar
- Groves CP: Primate Taxonomy. 2001, Washington: Smithsonian Institution PressGoogle Scholar
- Geissmann T: Taxonomy and evolution of gibbons. Evol Anthropol. 2002, 11: 28-31. 10.1002/evan.10047.View ArticleGoogle Scholar
- Mootnick AR: Gibbon (Hylobatidae) species identification recommended for rescue or breeding centers. Primate Conserv. 2006, 21: 103-138. 10.1896/0898-6188.8.131.52.View ArticleGoogle Scholar
- Geissmann T: Gibbon systematics and species identification. Int Zoo News. 1995, 42: 65-77.Google Scholar
- Cunningham C, Mootnick AR: Gibbons. Curr Biol. 2009, 19: R543-R544. 10.1016/j.cub.2009.05.013.View ArticlePubMedGoogle Scholar
- Geissmann T: Duet-splitting and the evolution of gibbon songs. Biol Rev. 2002, 77: 57-76.View ArticlePubMedGoogle Scholar
- Müller S, Hollatz M, Wienberg J: Chromosomal phylogeny and evolution of gibbons (Hylobatidae). Hum Genet. 2003, 113: 493-501. 10.1007/s00439-003-0997-2.View ArticlePubMedGoogle Scholar
- Roberto R, Capozzi O, Wilson RK, Mardis ER, Lomiento M, Tuzun E, Cheng Z, Mootnick AR, Archidiacono N, Rocchi M, Eichler EE: Molecular refinement of gibbon genome rearrangements. Genome Res. 2007, 17: 249-257. 10.1101/gr.6052507.PubMed CentralView ArticlePubMedGoogle Scholar
- Misceo D, Capozzi O, Roberto R, Dell'Oglio MP, Rocchi M, Stanyon R, Archidiancono N: Tracking the complex flow of chromosome rearrangements from the Hominoidea ancestor to extant Hylobates and Nomascus gibbons by high-resolution synteny mapping. Genome Res. 2008, 18: 1530-1537. 10.1101/gr.078295.108.PubMed CentralView ArticlePubMedGoogle Scholar
- Giriajan S, Chen L, Graves T, Marques-Bonet T, Ventura M, Fronick C, Fulton L, Rocchi M, Fulton RS, Wilson RK, Mardis ER, Eichler EE: Sequencing human-gibbon breakpoints of synteny reveals mosaic new insertions at rearrangement sites. Genome Res. 2009, 19: 178-190. 10.1101/gr.086041.108.View ArticleGoogle Scholar
- IUCN Red List of Threatened Species 2009.2. 2009, [http://www.iucnredlist.org]
- Chan BPL, Fellowes JR, Geissmann T, Zhang J: Hainan Gibbon Status Survey and Conservation Action Plan. Version 1. 2005, Hong Kong: Kadoorie Farm & Botanic GardenGoogle Scholar
- Mootnick AR, Xiaoming W, Moisson P, Chan BPL, Fellowes JR, Nadler T: Hainan Gibbon, Nomascus hainanus (Thomas, 1892). Primates in Peril: The World's 25 Most Endangered Primates, 2006-2008. Edited by: Mittermeier RA, Ratsimbazafy J, Rylands AB, Williamson L, Oates JF, Mbora D, Ganzhorn JU, Rodriguez-Luna E, Palacios E, Heymann EW, Cecilia M, Kierulff M, Yongcheng L, Supriatna J, Roos C, Walker S, Aguiar JM. 2007, Arlington: IUCN/SSC Primate Specialist Group (PSG), International Primatological Society (IPS), and Conservation International (CI), 16-17.Google Scholar
- Geissmann T, Dang Xuan Nguyen, Lormée N, Momberg F: Vietnam Primate Conservation Status Review. Part 1: Gibbon. 2000, Hanoi: Fauna & Flora International, Indochina ProgrammeGoogle Scholar
- Geissmann T: Status reassessment of the gibbons: results of the Asian primate red list workshop 2006. Gibbon J. 2007, 3: 5-15.Google Scholar
- Napier JR, Napier PH: A Handbook of Living Primates. 1967, London: Academic PressGoogle Scholar
- Groves CP: Systematics and phylogeny of gibbons. Gibbon and Siamang. Edited by: Rumbaugh DM. 1972, Basel: Karger, 1-89.Google Scholar
- Bruce EJ, Ayala FJ: Phylogenetic relationships between man and the apes: eletrophoretic evidence. Evolution. 1979, 33: 1040-1056. 10.2307/2407465.View ArticleGoogle Scholar
- Haimoff EH, Chivers DJ, Gittins SP, Whitten AJ: A phylogeny of gibbons (Hylobates spp.) based on morphological and behavioural characters. Folia Primatol. 1982, 39: 213-237. 10.1159/000156079.View ArticlePubMedGoogle Scholar
- Prouty LA, Buchanan PD, Pollitzer WS, Mootnick AR: Bunopithecus: A genus-level taxon for the hoolock gibbon (Hylobates hoolock). Am J Primatol. 1983, 5: 83-87. 10.1002/ajp.1350050110.View ArticleGoogle Scholar
- Creel N, Preuschoft H: Systematics of the lesser apes: A quantitative taxonomic analysis of craniometric and other variables. The Lesser Apes. Evolutionary and Behavioural Biology. Edited by: Preuschoft H, Chivers DJ, Brockelman WY, Creel N. 1984, Edinburgh: Edinburgh University Press, 562-613.Google Scholar
- Marshall JT, Sugardjito J: Gibbon systematics. Comparative Primate Biology, Systematics, Evolution and Anatomy. Edited by: Swindler DR, Erwin J. 1986, New York: Liss, 1: 137-185.Google Scholar
- Shafer DA: Evolutionary cytogenetics of the siabon (gibbon-siamang) hybrid apes. Current Perspectives in Primate Biology. Edited by: Taub DM, King FA. 1986, New York: Van Nostrand Reinhold Co, 226-239.Google Scholar
- Liu R, Shi L, Chen Y: A study on the chromosomes of white-browed gibbon (Hylobates hoolock leuconedys). Acta Theriol Sinica. 1987, 7: 1-7.Google Scholar
- Garza JC, Woodruff DS: A phylogenetic study of the gibbons (Hylobates) using DNA obtained non-invasively from hair. Mol Phylogenet Evol. 1992, 1: 202-210. 10.1016/1055-7903(92)90016-A.View ArticlePubMedGoogle Scholar
- Hayashi S, Hayasaka K, Takenaka O, Horai S: Molecular phylogeny of gibbons inferred from mitochondrial DNA sequences: preliminary report. J Mol Evol. 1995, 41: 359-365. 10.1007/BF01215183.View ArticlePubMedGoogle Scholar
- Hall LM, Jones DS, Wood BA: Evolution of the gibbon subgenera inferred from cytochrome b DNA sequence data. Mol Phylogenet Evol. 1998, 10: 281-286. 10.1006/mpev.1998.0539.View ArticlePubMedGoogle Scholar
- Roos C, Geissmann T: Molecular phylogeny of the major hylobatid divisions. Mol Phylogenet Evol. 2001, 19: 486-494. 10.1006/mpev.2001.0939.View ArticlePubMedGoogle Scholar
- Roos C: Molecular evolution and systematics of Vietnamese primates. Conservation of Primates in Vietnam. Edited by: Nadler T, Streicher U, Ha Thang Long. 2004, Hanoi: Endangered Primate Rescue Center, Frankfurt Zoological Society, 23-28.Google Scholar
- Takacs Z, Morales JC, Geissmann T, Melnick DJ: A complete species-level phylogeny of the Hylobatidae based on mitochodrial ND3-ND4 gene sequence. Mol Phylogenet Evol. 2005, 36: 456-467. 10.1016/j.ympev.2005.03.032.View ArticlePubMedGoogle Scholar
- Chatterjee HJ: Phylogeny and biogeography of gibbons: a dispersal-vicariance analysis. Int J Primatol. 2006, 27: 699-712. 10.1007/s10764-006-9044-1.View ArticleGoogle Scholar
- Monda K, Simmons RE, Kressirer P, Su B, Woodruff DS: Mitochondrial DNA hypervariable region-1 sequence variation and phylogeny of the concolor gibbons, Nomascus. Am J Primatol. 2007, 69: 1-22. 10.1002/ajp.20439.View ArticleGoogle Scholar
- Roos C, Thanh Vu Ngoc, Walter L, Nadler T: Molecular systematics of Indochinese primates. Vietn J Primatol. 2007, 1 (1): 41-53.Google Scholar
- Whittaker DJ, Morales JC, Melnick DJ: Resolution of the Hylobates phylogeny: congruence of mitochondrial D-loop sequences with molecular, behavioral, and morphological datasets. Mol Phylogenet Evol. 2007, 45: 620-628. 10.1016/j.ympev.2007.08.009.View ArticlePubMedGoogle Scholar
- Schultz AH: Observations on the growth, classification and evolutionary specialization of gibbons and siamangs. Human Biol. 1933, 5: 212-255. 385-428Google Scholar
- Brandon-Jones D, Eudey AA, Geissmann T, Groves CP, Melnick DJ, Morales JC, Shekelle M, Stewart CB: Asian primate classification. Int J Primatol. 2004, 25: 97-164. 10.1023/B:IJOP.0000014647.18720.32.View ArticleGoogle Scholar
- Mootnick AR, Groves CP: A new generic name for the hoolock gibbon (Hylobatidae). Int J Primatol. 2005, 26: 971-976. 10.1007/s10764-005-5332-4.View ArticleGoogle Scholar
- Chivers DJ: The lesser apes. Primate Conservation. Edited by: Prince Rainier III of Monaco, Bourne GH. 1977, London & New York: Academic Press, 539-598.Google Scholar
- Goodman M, Porter CA, Czelusniak J, Page SL, Schneider H, Shoshani J, Gunnell G, Groves CP: Toward a phylogenetic classification of primates based on DNA evidence complemented by fossil evidence. Mol Phylogenet Evol. 1998, 9: 585-598. 10.1006/mpev.1998.0495.View ArticlePubMedGoogle Scholar
- The Gibbon Research Lab. [http://www.gibbons.de]
- Raaum RL, Sterner KN, Noviello CM, Stewart CB, Disotell TR: Catarrhine primate divergence dates estimated from complete mitochondrial genomes: concordance with fossil and nuclear DNA evidence. J Hum Evol. 2005, 48: 237-257. 10.1016/j.jhevol.2004.11.007.View ArticlePubMedGoogle Scholar
- Sterner KN, Raaum RL, Zhang YP, Stewart CB, Disotell TR: Mitochondrial data support an odd-nosed colobine clade. Mol Phylogenet Evol. 2006, 40: 1-7. 10.1016/j.ympev.2006.01.017.View ArticlePubMedGoogle Scholar
- Tosi AJ, Morales JC, Melnick DJ: Paternal, maternal, and biparental molecular markers provide unique windows onto the evolutionary history of macaque monkeys. Evolution. 2003, 57: 1419-1435.View ArticlePubMedGoogle Scholar
- Ziegler T, Abegg C, Meijaard E, Perwitasari-Farajallah D, Walter L, Hodges JK, Roos C: Molecular phylogeny and evolutionary history of Southeast Asian macaques forming the M. silenus group. Mol Phylogenet Evol. 2007, 42: 807-816. 10.1016/j.ympev.2006.11.015.View ArticlePubMedGoogle Scholar
- Roos C, Nadler T, Walter L: Mitochondrial phylogeny, taxonomy and biogeography of the silvered langur species group (Trachypithecus cristatus). Mol Phylogenet Evol. 2008, 47: 629-636. 10.1016/j.ympev.2008.03.006.View ArticlePubMedGoogle Scholar
- Hirai H, Hayano A, Tanaka H, Mootnick AR, Wijayanto H, Perwitasari-Farajallah D: Genetic differentiation of agile gibbons between Sumatra and Kalimantan in Indonesia. The Gibbons. New Perspectives on Small Ape Socioecology and Population Biology. Edited by: Lappan S, Whittaker DJ. 2009, New York: Springer, 37-49.Google Scholar
- Carbone L, Mootnick A, Nadler T, Moisson P, Ryder O, Roos C, de Jong PJ: A chromosomal inversion unique to the northern white-cheeked gibbon. PLoS ONE. 2009, 4: e4999-10.1371/journal.pone.0004999.PubMed CentralView ArticlePubMedGoogle Scholar
- Eudey AA: Pleistocene glacial phenomena and the evolution of Asian macaques. The Macaques: Studies in Ecology, Behavior and Evolution. Edited by: Lindburg DG. 1980, New York: Van Nostrand Rheinhold, 52-83.Google Scholar
- Morley RJ, Flenley JR: Late Cainozoic vegetational and environmental changes in the Malay Archipelago. Biogeographical Evolution of the Malay Archipelago. Edited by: Whitmore TC. 1987, Oxford: Oxford Scientific Publications, 50-59.Google Scholar
- Morley RJ: Origin and Evolution of Tropical Rain Forests. 2000, West Sussex: John Wiley & SonsGoogle Scholar
- Bird MI, Taylor D, Hunt C: Palaeoenvironments of insular Southeast Asia during the last glacial period: a savanna corridor in Sundaland?. Quaternary Sci Rev. 2005, 24: 2228-2242. 10.1016/j.quascirev.2005.04.004.View ArticleGoogle Scholar
- Meijaard E: Solving mammalian riddles. A reconstruction of the Tertiary and Quaternary distribution of mammals and their palaeoenvironments in island South-East Asia. PhD thesis. 2004, Australian National University, Anthropology and Archaeology DepartmentGoogle Scholar
- Meijaard E, Groves CP: The geography of mammals and rivers in mainland Southeast Asia. Primate Biogeography. Edited by: Lehman SM, Fleagle JG. 2006, New York: Springer, 305-329. full_text.View ArticleGoogle Scholar
- Heaney LR: A synopsis of climatic and vegetational change in Southeast Asia. Climatic Change. 1991, 19: 53-61. 10.1007/BF00142213.View ArticleGoogle Scholar
- Urushibara-Yoshino K, Yoshino M: Palaeoenvironmental change in Java island and its surrounding areas. J Quaternary Sci. 1997, 12: 435-442. 10.1002/(SICI)1099-1417(199709/10)12:5<435::AID-JQS325>3.0.CO;2-E.View ArticleGoogle Scholar
- Kaars van der S: Pollen distribution in marine sediments from the south-eastern Indonesian waters. Palaeogeography, Palaeoclimatology, Palaeoecology. 2001, 171: 341-361. 10.1016/S0031-0182(01)00253-X.View ArticleGoogle Scholar
- Cannon CH, Morley RJ, Bush ABG: The current refugial rainforests of Sundaland are unrepresentative of their biogeographic past and highly vulnerable to disturbance. Proc Natl Acad Sci USA. 2009, 106: 11188-11193. 10.1073/pnas.0809865106.PubMed CentralView ArticlePubMedGoogle Scholar
- Jablonski NG, Whitfort MJ: Environmental changes during the Quaternary in East Asia and its consequences for mammals. Records of the Western Australian Museum. 1999, 57: 307-315.Google Scholar
- Meijaard E: Mammals of South-East Asian islands and their late Pleistocene environments. J Biogeography. 2003, 30: 1245-1257. 10.1046/j.1365-2699.2003.00890.x.View ArticleGoogle Scholar
- Miller KG, Kominz MA, Browning JV, Wright JD, Mountain GS, Katz ME, Sugarman PJ, Cramer BS, Christie-Blick N, Pekar SF: The phanerozoic record of global sea-level change. Science. 2005, 310: 1293-1298. 10.1126/science.1116412.View ArticlePubMedGoogle Scholar
- Lisiecki LE, Raymo ME: A Pliocene-Pleistocene stack of 57 globally distributed benthic δ18 O records. Paleoceanography. 2005, 20: PA1003-10.1029/2004PA001071.Google Scholar
- Naish TR, Wilson GS: Constraints on the amplitude of Mid-Pliocene (3.6-2.4?Ma) eustatic sea-level fluctuations from the New Zealand shallow-marine sediment record. Phil Trans R Soc A. 2009, 367: 169-187. 10.1098/rsta.2008.0223.View ArticlePubMedGoogle Scholar
- Woodruff DS: Biogeography and conservation in Southeast Asia: how 2.7 million years of repeated environmental fluctuations affect today's patterns and the future of the remaining refugial-based biodiversity. Biodivers Conserv.
- Verstappen HT: On paleo-climates and landform development in Malesia. Modern Quaternary Research in Southeast Asia. Edited by: Barstra GJ, Casparie WA. 1975, Rotterdam: Balkema, 3-35.Google Scholar
- Tougard C: Biogeography and migration routes of large mammal faunas in South-East Asia during the late middle Pleistocene: focus on the fossil and extant faunas from Thailand. Palaeogeography, Palaeoclimatology, Palaeoecology. 2001, 168: 337-358. 10.1016/S0031-0182(00)00243-1.View ArticleGoogle Scholar
- Woodruff DS, Turner LM: The Indochinese-Sundaic zoogeographic transition: a description and analysis of terrestrial mammal species distributions. J Biogeography. 2009, 36: 803-821. 10.1111/j.1365-2699.2008.02071.x.View ArticleGoogle Scholar
- Chatterjee HJ: Evolutionary relationships among the gibbons: a biogeographic perspective. Gibbons: New Perspectives on Small Ape Socioecology and Population Biology. Edited by: Lappan S, Whittaker DJ. 2009, New York: Springer, 13-36.View ArticleGoogle Scholar
- Harrison T, Krigbaum J, Manser J: Primate biogeography and ecology on the Sunda Shelf islands: a paleontological and zooarchaeological perspective. Primate Biogeography. Edited by: Lehman SM, Fleagle JG. 2006, New York: Springer, 331-374. full_text.View ArticleGoogle Scholar
- Jablonski NG, Chaplin G: The fossil record of gibbons. Gibbons: New Perspectives on Small Ape Socioecology and Population Biology. Edited by: Lappan S, Whittaker DJ. 2009, New York: Springer, 111-130.View ArticleGoogle Scholar
- Peng YZ, Pan RL, Jablonski NG: Classification and evolution of Asian colobines. Folia Primatol. 1993, 60: 106-117. 10.1159/000156680.View ArticlePubMedGoogle Scholar
- Jablonski NG: Natural History of the Doucs and Snub-nosed Monkeys. 1998, New Jersey: World Scientific Publishing CompanyView ArticleGoogle Scholar
- Hallet B, Molnar P: Distorted drainage basins as markers of crustal strain east of the Himalayas. J Geophysical Res. 2001, 106: 13697-13709. 10.1029/2000JB900335.View ArticleGoogle Scholar
- Singh AD, Srinivasan MS: Quaternary climate changes indicated by planktonic forminifera of northern Indian ocean. Curr Sci. 1993, 64: 908-915.Google Scholar
- Hirai H, Wijayanto H, Tanaka H, Mootnick AR, Hayano A, Perwitasari-Farajallah D, Iskandriati D, Sajuthi D: A whole-arm translocation (WAT8/9) separating Sumatran and Bornean agile gibbons, and its evolutionary features. Chromosome Res. 2005, 13: 123-133. 10.1007/s10577-005-7475-8.View ArticlePubMedGoogle Scholar
- Lappan S: Patterns of dispersal in Sumatran siamangs (Symphalangus syndactylus): preliminary mtDNA evidence suggests more frequent male than female dispersal to adjacent groups. Am J Primatol. 2007, 69: 692-698. 10.1002/ajp.20382.View ArticlePubMedGoogle Scholar
- Long Y, Nadler T: Eastern black crested gibbon Nomascus nasutus (Kunkel d'Herculais, 1884). Primates in Peril: The World's 25 Most Endangered Primates 2008-2010. Edited by: Mittermeier RA, Wallis J, Rylands AB, Ganzhorn JU, Oates JF, Williamson EA, Palacios E, Heymann EW, Kierulff MCM, Long Y, Supriatna J, Roos C, Walker S, Cortés-Ortiz L, Schwitzer C. 2009, Arlington: IUCN/SSC Primate Specialist Group (PSG), International Primatological Society (IPS), and Conservation International (CI), 60-61.Google Scholar
- Drummond AJ, Kearse M, Heled J, Moir R, Thierer T, Ashton B, Wilson A, Stones-Havas S: Geneious, version 4.6.1. 2008, [http://www.geneious.com]Google Scholar
- Hofreiter M, Siedel H, Van Neer W, Vigilant L: Mitochondrial DNA sequence from an enigmatic gorilla population (Gorilla gorilla uellensis). Am J Phys Anthropol. 2003, 121: 361-368. 10.1002/ajpa.10186.View ArticlePubMedGoogle Scholar
- Thalmann O, Hebler J, Poinar HN, Pääbo S, Vigilant L: Unreliable mtDNA data due to nuclear insertions: a cautionary tale from analysis of humans and other great apes. Mol Ecol. 2004, 13: 321-335. 10.1046/j.1365-294X.2003.02070.x.View ArticlePubMedGoogle Scholar
- Swofford DL: PAUP*: Phylogenetic analysis using parsimony (*and other methods), Version 4.0b10. 2003, Sunderland: Sinauer AssociatesGoogle Scholar
- Zwickl DJ: Genetic algorithm approaches for the phylogenetic analysis of large biological sequence data sets under the maximum likelihood criterion. PhD thesis. 2006, Texas University, AustinGoogle Scholar
- Huelsenbeck JP, Ronquist F, Nielsen R, Bollback JP: Bayesian inference of phylogeny and its impact on evolutionary biology. Science. 2001, 294: 2310-2314. 10.1126/science.1065889.View ArticlePubMedGoogle Scholar
- Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.View ArticlePubMedGoogle Scholar
- Posada D, Crandall KA: Modeltest: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.View ArticlePubMedGoogle Scholar
- Drummond AJ, Ho SYW, Phillips MJ, Rambaut A: Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006, 4: e88-10.1371/journal.pbio.0040088.PubMed CentralView ArticlePubMedGoogle Scholar
- Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: e214-10.1186/1471-2148-7-214.View ArticleGoogle Scholar
- Vignaud P, Duringer P, Mackaye HT, Likius A, Blondel C, Boisserie JR, De Bonis L, Eisenmann V, Etienne ME, Geraads D, Guy F, Lehmann T, Lihoreau F, Lopez-Martinez N, Mourer-Chauvire C, Otero O, Rage JC, Schuster M, Viriot L, Zazzo A, Brunet M: Geology and palaeontology of the upper Miocene Toros-Menalla hominid locality, Chad. Nature. 2002, 418: 152-155. 10.1038/nature00880.View ArticlePubMedGoogle Scholar
- Brunet M, Guy F, Pilbeam D, Lieberman DE, Likius A, Mackaye HT, Ponce de León MS, Zollikofer CP, Vignaud P: New material of the earliest hominid from the upper Miocene of Chad. Nature. 2005, 434: 752-755. 10.1038/nature03392.View ArticlePubMedGoogle Scholar
- Lebatard AE, Bourles DL, Duringer P, Jolivet M, Braucher R, Carcaillet J, Schuster M, Arnaud N, Monie P, Lihoreau F, Likius A, Mackaye HT, Vignaud P, Brunet M: Cosmogenic nuclide dating of Sahelanthropus tchadensis and Australopithecus bahrelghazali: Mio-Pliocene hominids from Chad. Proc Natl Acad Sci USA. 2008, 105: 3226-3231. 10.1073/pnas.0708015105.PubMed CentralView ArticlePubMedGoogle Scholar
- Kelley J: The hominoid radiation in Asia. The Primate Fossil Record. Edited by: Hartwig WC. 2002, Cambridge: Cambridge University Press, 369-384.Google Scholar
- Benefit BR, McCrossin ML: The Victoriapithecidae, Cercopithecoidea. The Primate Fossil Record. Edited by: Hartwig WC. 2002, Cambridge: Cambridge University Press, 241-253.Google Scholar
- Young NM, MacLatchy L: The phylogenetic position of Mortopithecus. J Hum Evol. 2004, 46: 163-184. 10.1016/j.jhevol.2003.11.002.View ArticlePubMedGoogle Scholar
- Rambaut A, Drummond AJ: Tracer v1.4.1: MCMC trace analysis tool. 2007, Institute of Evolutionary Biology, University of EdinburghGoogle Scholar
- Rambaut A: FigTree: Tree figure drawing tool, version 1.2.2. 2008, Institute of Evolutionary Biology, University of EdinburghGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.