Molecular phylogeny of the major hylobatid divisions

Landscape and conservation genetics of western black crested gibbons (Nomascus concolor) in China

Despite decades of field study, very little is known about the molecular ecology of gibbons, particularly as it relates to their ability to disperse across degraded and fragmentary landscapes. The critically endangered western black crested gibbon (Nomascus concolor) has been reduced to a small, fragmented population with about 1300 individuals. In the largest population genetic study of free-ranging gibbons to date, we sampled 47 of these gibbons from 13 sites in China and generated 15 polymorphic autosomal microsatellite markers. We identify three population clusters of N. concolor in Yunnan centered in 1) the Wuliang and Ailao Mountains, 2) the Yongde Daxueshan Mountains, and 3) an isolated remnant near the border with Vietnam. Within the Wuliang Mountains, we identified four subclusters, three of which are bounded by high-altitude rhododendron forest, and one that is isolated from the main population by ~2 km of degraded forest and pasture. Least-cost path analysis and isolation by resistance modeling demonstrates that the population genetic distances among gibbons in Wuliangshan National Nature Reserve are significantly correlated with geographic paths that avoid use of high-altitude rhododendron forest in favor of evergreen broadleaf forest. Although these gibbons have likely undergone reductions in heterozygosity from recent consanguineous mating, we suggest that their active avoidance of inbreeding on the population level maintains higher than expected levels of genetic diversity. This research provides new insights into how gibbons interact with heterogeneous environments and expands our understanding of their molecular ecology and conservation genetics.

Molecular phylogeny confirms the subspecies delineation of the Malayan Siamang (Symphalangussyndactyluscontinentis) and the Sumatran Siamang (Symphalangussyndactylussyndactylus) based on the hypervariable region of mitochondrial DNA

Siamangs (Symphalangussyndactylus) are native to Peninsular Malaysia, Sumatra and southern Thailand and their taxonomical classification at subspecies level remains unclear. Morphologically, two subspecies were proposed as early as 1908 by Thomas namely Symphalangus s.syndactylus and Symphalanguss.continentis. Thus, this study aims to clarify the Siamang subspecies status, based on mtDNA D-loop sequences. Faecal samples were collected from wild Siamang populations at different localities in Peninsular Malaysia. A 600-bp sequence of the mitochondrial D-loop region was amplified from faecal DNA extracts and analysed along with GenBank sequences representing Symphalangus sp., Nomascus sp., Hylobates sp., Hoolock sp. and outgroups (Pongopygmaeus, Macacafascicularis and Papiopapio). The molecular phylogenetic analysis in this study revealed two distinct clades formed by S.s.syndactylus and S.s.continentis which supports the previous morphological delineation of the existence of two subspecies. Biogeographical analysis indicated that the Sumatran population lineage was split from the Peninsular Malaysian population lineage and a diversification occurrred in the Pliocene era (~ 3.12 MYA) through southward expansion. This postulation was supported by the molecular clock, which illustrated that the Peninsular Malaysian population (~ 1.92 MYA) diverged earlier than the Sumatran population (~ 1.85 MYA). This is the first study to use a molecular approach to validate the subspecies statuses of S.s.syndactylus and S.s.continentis. This finding will be useful for conservation management, for example, during Siamang translocation and investigations into illegal pet trade and forensics involving Malayan and Sumatran Siamangs.

Density and population size estimates of the endangered northern yellow-cheeked crested gibbon Nomascus annamensis in selectively logged Veun Sai-Siem Pang National Park in Cambodia using acoustic spatial capture-recapture methods

Many gibbon species are threatened with extinction, including the endangered northern yellow-cheeked crested gibbon, Nomascus annamensis. Assessing gibbon populations and understanding how human disturbances and environmental factors impact these populations is vital for effective conservation planning. In 2010, auditory surveys revealed that Veun Sai-Siem Pang National Park (VSSP) in Cambodia contains one of the largest known N. annamensis populations in the world, with an estimated 456 (95% CI 421-490) gibbon groups. Illegal selective logging is common in the park, but the impact of continued logging on the gibbon population has not been investigated. To determine any change in the N. annamensis population since 2010, between January and April 2019 we conducted auditory surveys at 13 sites that were at least 4 km apart. We surveyed each site for three days, each day recording the gibbon calls heard over 3.25 hours from three listening posts located 500 m apart. At the same sites, we assessed the logging intensity using transects and ecological plots. Gibbon densities can be influenced by various environmental factors such as canopy height and forest type. Therefore, in addition to investigating the relationship between the density of N. annamensis groups and logging, we included five additional environmental variables in our acoustic spatial capture-recapture models. Our best fit model with the lowest AIC value included canopy height, forest type, distance to villages, and logging. We estimate that there are 389 (95% CI 284-542) N. annamensis groups currently in VSSP. Selective logging is widespread in the park, primarily targeting four tree species. The estimated felling time of these logged trees, together with previous reports, indicate that the species most targeted in VSSP varies over time. To conserve the N. annamensis population in VSSP, it is crucial that action is taken to reduce illegal logging.

Taxonomic and genetic assessment of captive White-Handed Gibbons (Hylobateslar) in Peninsular Malaysia with implications towards conservation translocation and reintroduction programmes

Conservation translocation and reintroduction for the purpose of repopulating and reinforcing extirpated or depleted populations has been recognised as an important conservation tool, particularly for gibbon conservation in the immediate future. Feasibility assessments involving multiple factors, including taxonomic and genetic assessment of rescued and captive gibbons, are imperative prior to translocation and reintroduction programmes. In this study, we attempt to determine the subspecies and origin of captive Hylobateslar, White-handed gibbons, from Peninsular Malaysia to assist in future translocation and reintroduction programmes. A total of 12 captive and rescued H.lar samples were analysed using the control region segment of mitochondrial DNA. Sequence analyses and phylogenetic trees constructed using neighbour-joining, maximum likelihood, Bayesian inference, and network methods congruently differentiate all 12 captive individuals used in this study from other H.lar subspecies suggesting that these individuals belong to the H.larlar subspecies. In addition, two populations of H.l.lar were observed: (1) a southern population consisting of all 12 individuals from Peninsular Malaysia, and (2) a possible northern population represented by three individuals (from previous studies), which might have originated from the region between the Isthmus of Kra, Surat Thani-Krabi depression, and Kangar-Pattani. Our findings suggest that the complete control region segment can be used to determine the subspecies and origin of captive H.lar.

A first report of separation calls in southern yellow-cheeked gibbons (Nomascus gabriellae) in captivity

The effects of social separation, including vocalization, have been studied for a very long time in non-human primates under laboratory conditions. As part of the long-term research on the vocal behaviour of Nomascus gibbons in zoos, this study provides the first record of calls of the southern yellow-cheeked gibbon (Nomascus gabriellae) in response to involuntary separation. Our study revealed that calls were also emitted by an infant (aged 1 year 8 months), and that the acoustic structure of the infant's calls was similar to that of older individuals' calls. Separation-induced calls seem to have a shorter developmental convergence than vocalizations with a stable pattern (which are specific for species and sex). The acoustic structure of the calls reported here comprised simple syllables, and differed from the sex- and species-specific vocal patterns of this species. Our findings demonstrate a novel paradigm in this genus, and provide evidence of the ability of gibbons to express distress when socially separated.

Species identification of crested gibbons (Nomascus) in captivity in China using karyotyping- and PCR-based approaches

Gibbons and siamangs (Hylobatidae) are well-known for their rapid chromosomal evolution, which has resulted in high speciation rate within the family. On the other hand, distinct karyotypes do not prevent speciation, allowing interbreeding between individuals in captivity, and the unwanted hybrids are ethically problematic as all gibbon species are endangered or critically endangered. Thus, accurate species identification is crucial for captive breeding, particularly in China where studbooks are unavailable. Identification based on external morphology is difficult, especially for hybrids, because species are usually similar in appearance. In this study, we employed G-banding karyotyping and fluorescence in situ hybridization (FISH) as well as a PCR-based approach to examine karyotypic characteristics and identify crested gibbons of the genus Nomascus from zoos and nature reserves in China. We characterized and identified five karyotypes from 21 individuals of Nomascus. Using karyotypes and mitochondrial and nuclear genes, we identified three purebred species and three hybrids, including one F2 hybrid between N. gabriellae and N. siki. Our results also supported that N. leucogenys and N. siki shared the same inversion on chromosome 7, which resolves arguments from previous studies. Our results demonstrated that both karyotyping and DNA-based approaches were suitable for identifying purebred species, though neither was ideal for hybrid identification. The advantages and disadvantages of both approaches are discussed. Our results further highlight the importance of animal ethics and welfare, which are critical for endangered species in captivity.

A computational reconstruction of Papio phylogeny using Alu insertion polymorphisms

Background

Since the completion of the human genome project, the diversity of genome sequencing data produced for non-human primates has increased exponentially. Papio baboons are well-established biological models for studying human biology and evolution. Despite substantial interest in the evolution of Papio, the systematics of these species has been widely debated, and the evolutionary history of Papio diversity is not fully understood. Alu elements are primate-specific transposable elements with a well-documented mutation/insertion mechanism and the capacity for resolving controversial phylogenetic relationships. In this study, we conducted a whole genome analysis of Alu insertion polymorphisms unique to the Papio lineage. To complete these analyses, we created a computational algorithm to identify novel Alu insertions in next-generation sequencing data.

Results

We identified 187,379 Alu insertions present in the Papio lineage, yet absent from M. mulatta [Mmul8.0.1]. These elements were characterized using genomic data sequenced from a panel of twelve Papio baboons: two from each of the six extant Papio species. These data were used to construct a whole genome Alu-based phylogeny of Papio baboons. The resulting cladogram fully-resolved relationships within Papio.

Conclusions

These data represent the most comprehensive Alu-based phylogenetic reconstruction reported to date. In addition, this study produces the first fully resolved Alu-based phylogeny of Papio baboons.

Electronic supplementary material

The online version of this article (10.1186/s13100-018-0118-3) contains supplementary material, which is available to authorized users.

The Taxonomic and Phylogenetic Affinities of Bunopithecus sericus, a Fossil Hylobatid from the Pleistocene of China

Fossil hylobatids are rare, but are known from late Miocene and Pleistocene sites throughout East Asia. The best-known fossil hylobatid from the Pleistocene of China is a left mandibular fragment with M2-3 (AMNH 18534), recovered from a pit deposit near the village of Yanjinggou in Wanzhou District, Chongqing Province. Matthew and Granger described this specimen in 1923 as a new genus and species, Bunopithecus sericus. Establishing the age of Bunopithecus has proved difficult because the Yanjinggou collection represents a mixed fauna of different ages, but it likely comes from early or middle Pleistocene deposits. Although the Bunopithecus specimen has featured prominently in discussions of hylobatid evolution and nomenclature, its systematic status has never been satisfactorily resolved. The present study reexamines the taxonomic and phylogenetic relationships of Bunopithecus by carrying out a detailed comparative morphometric study of its lower molars in relation to a large sample of modern hylobatids. Our results show that differences in M2 and M3 discriminate extant hylobatids fairly well, at least at the generic level, and that AMNH 18534 is not attributable to Hylobates, Nomascus or Symphalangus. Support for a close relationship between Bunopithecus and Hoolock is more equivocal. In most multivariate analyses, Bunopithecus presents a unique morphological pattern that falls outside the range of variation of any hylobatid taxon, although its distance from the cluster represented by extant hoolocks is relatively small. Our results support the generic distinction of Bunopithecus, which most likely represents an extinct crown hylobatid, and one that may possibly represent the sister taxon to Hoolock.

Recent advances in primate phylogenomics

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

Locational diversity of alpha satellite DNA and intergeneric hybridization aspects in the Nomascus and Hylobates genera of small apes

Recently, we discovered that alpha satellite DNA has unique and genus-specific localizations on the chromosomes of small apes. This study describes the details of alpha satellite localization in the genera Nomascus and Hylobates and explores their usefulness in distinguishing parental genome sets in hybrids between these genera. Fluorescence in situ hybridization was used to establish diagnostic criteria of alpha satellite DNA markers in discriminating small ape genomes. In particular we established the genus specificity of alpha satellite distribution in three species of light-cheeked gibbons (Nomascus leucogenys, N. siki, and N. gabriellae) in comparison to that of Hylobates lar. Then we determined the localization of alpha satellite DNA in a hybrid individual which resulted from a cross between these two genera. In Nomascus the alpha satellite DNA blocks were located at the centromere, telomere, and four interstitial regions. In Hylobates detectable amounts of alpha satellite DNA were seen only at centromeric regions. The differences in alpha satellite DNA locations between Nomascus and Hylobates allowed us to easily distinguish the parental chromosomal sets in the genome of intergeneric hybrid individuals found in Thai and Japanese zoos. Our study illustrates how molecular cytogenetic markers can serve as diagnostic tools to identify the origin of individuals. These molecular tools can aid zoos, captive breeding programs and conservation efforts in managing small apes species. Discovering more information on alpha satellite distribution is also an opportunity to examine phylogenetic and evolutionary questions that are still controversial in small apes.

Evolution of life history and behavior in Hominidae: towards phylogenetic reconstruction of the chimpanzee-human last common ancestor

The origin of the fundamental behavioral differences between humans and our closest living relatives is one of the central issues of evolutionary anthropology. The prominent, chimpanzee-based referential model of early hominin behavior has recently been challenged on the basis of broad multispecies comparisons and newly discovered fossil evidence. Here, we argue that while behavioral data on extant great apes are extremely relevant for reconstruction of ancestral behaviors, these behaviors should be reconstructed trait by trait using formal phylogenetic methods. Using the widely accepted hominoid phylogenetic tree, we perform a series of character optimization analyses using 65 selected life-history and behavioral characters for all extant hominid species. This analysis allows us to reconstruct the character states of the last common ancestors of Hominoidea, Hominidae, and the chimpanzee-human last common ancestor. Our analyses demonstrate that many fundamental behavioral and life-history attributes of hominids (including humans) are evidently ancient and likely inherited from the common ancestor of all hominids. However, numerous behaviors present in extant great apes represent their own terminal autapomorphies (both uniquely derived and homoplastic). Any evolutionary model that uses a single extant species to explain behavioral evolution of early hominins is therefore of limited use. In contrast, phylogenetic reconstruction of ancestral states is able to provide a detailed suite of behavioral, ecological and life-history characters for each hypothetical ancestor. The living great apes therefore play an important role for the confident identification of the traits found in the chimpanzee-human last common ancestor, some of which are likely to represent behaviors of the fossil hominins.

A mitogenomic phylogeny of living primates

Primates, the mammalian order including our own species, comprise 480 species in 78 genera. Thus, they represent the third largest of the 18 orders of eutherian mammals. Although recent phylogenetic studies on primates are increasingly built on molecular datasets, most of these studies have focused on taxonomic subgroups within the order. Complete mitochondrial (mt) genomes have proven to be extremely useful in deciphering within-order relationships even up to deep nodes. Using 454 sequencing, we sequenced 32 new complete mt genomes adding 20 previously not represented genera to the phylogenetic reconstruction of the primate tree. With 13 new sequences, the number of complete mt genomes within the parvorder Platyrrhini was widely extended, resulting in a largely resolved branching pattern among New World monkey families. We added 10 new Strepsirrhini mt genomes to the 15 previously available ones, thus almost doubling the number of mt genomes within this clade. Our data allow precise date estimates of all nodes and offer new insights into primate evolution. One major result is a relatively young date for the most recent common ancestor of all living primates which was estimated to 66-69 million years ago, suggesting that the divergence of extant primates started close to the K/T-boundary. Although some relationships remain unclear, the large number of mt genomes used allowed us to reconstruct a robust primate phylogeny which is largely in agreement with previous publications. Finally, we show that mt genomes are a useful tool for resolving primate phylogenetic relationships on various taxonomic levels.

Inferring the evolutionary histories of divergences in Hylobates and Nomascus gibbons through multilocus sequence data

BACKGROUND

Gibbons (Hylobatidae) are the most diverse group of living apes. They exist as geographically-contiguous species which diverged more rapidly than did their close relatives, the great apes (Hominidae). Of the four extant gibbon genera, the evolutionary histories of two polyspecific genera, Hylobates and Nomascus, have been the particular focus of research but the DNA sequence data used was largely derived from the maternally inherited mitochondrial DNA (mtDNA) locus.

RESULTS

To investigate the evolutionary relationships and divergence processes of gibbon species, particularly those of the Hylobates genus, we produced and analyzed a total of 11.5 kb DNA of sequence at 14 biparentally inherited autosomal loci. We find that on average gibbon genera have a high average sequence diversity but a lower degree of genetic differentiation as compared to great ape genera. Our multilocus species tree features H. pileatus in a basal position and a grouping of the four Sundaic island species (H. agilis, H. klossii, H. moloch and H. muelleri). We conducted pairwise comparisons based on an isolation-with-migration (IM) model and detect signals of asymmetric gene flow between H. lar and H. moloch, between H. agilis and H. muelleri, and between N. leucogenys and N. siki.

CONCLUSIONS

Our multilocus analyses provide inferences of gibbon evolutionary histories complementary to those based on single gene data. The results of IM analyses suggest that the divergence processes of gibbons may be accompanied by gene flow. Future studies using analyses of multi-population model with samples of known provenance for Hylobates and Nomascus species would expand the understanding of histories of gene flow during divergences for these two gibbon genera.

Incomplete lineage sorting is common in extant gibbon genera

We sequenced reduced representation libraries by means of Illumina technology to generate over 1.5 Mb of orthologous sequence from a representative of each of the four extant gibbon genera (Nomascus, Hylobates, Symphalangus, and Hoolock). We used these data to assess the evolutionary relationships between the genera by evaluating the likelihoods of all possible bifurcating trees involving the four taxa. Our analyses provide weak support for a tree with Nomascus and Hylobates as sister taxa and with Hoolock and Symphalangus as sister taxa, though bootstrap resampling suggests that other phylogenetic scenarios are also possible. This uncertainty is due to short internal branch lengths and extensive incomplete lineage sorting across taxa. The true phylogenetic relationships among gibbon genera will likely require a more extensive whole-genome sequence analysis.

A comparative analysis of Y chromosome and mtDNA phylogenies of the Hylobates gibbons

Background

The evolutionary relationships of closely related species have long been of interest to biologists since these species experienced different evolutionary processes in a relatively short period of time. Comparison of phylogenies inferred from DNA sequences with differing inheritance patterns, such as mitochondrial, autosomal, and X and Y chromosomal loci, can provide more comprehensive inferences of the evolutionary histories of species. Gibbons, especially the genus Hylobates, are particularly intriguing as they consist of multiple closely related species which emerged rapidly and live in close geographic proximity. Our current understanding of relationships among Hylobates species is largely based on data from the maternally-inherited mitochondrial DNAs (mtDNAs).

Results

To infer the paternal histories of gibbon taxa, we sequenced multiple Y chromosomal loci from 26 gibbons representing 10 species. As expected, we find levels of sequence variation some five times lower than observed for the mitochondrial genome (mtgenome). Although our Y chromosome phylogenetic tree shows relatively low resolution compared to the mtgenome tree, our results are consistent with the monophyly of gibbon genera suggested by the mtgenome tree. In a comparison of the molecular dating of divergences and on the branching patterns of phylogeny trees between mtgenome and Y chromosome data, we found: 1) the inferred divergence estimates were more recent for the Y chromosome than for the mtgenome, 2) the species H. lar and H. pileatus are monophyletic in the mtgenome phylogeny, respectively, but a H. pileatus individual falls into the H. lar Y chromosome clade.

Conclusions

Based on the ~6.4 kb of Y chromosomal DNA sequence data generated for each of the 26 individuals in this study, we provide molecular inferences on gibbon and particularly on Hylobates evolution complementary to those from mtDNA data. Overall, our results illustrate the utility of comparative studies of loci with different inheritance patterns for investigating potential sex specific processes on the evolutionary histories of closely related taxa, and emphasize the need for further sampling of gibbons of known provenance.

Phylogenetic separation in limb use in captive gibbons (Hylobatidae): a comparison across the primate order

Although there have been few studies of self-scratching in primates, some have reported distinct differences in whether hands or feet are used, and these variations seem to reflect the evolutionary history of the Order. Monkeys and prosimians use both hands and feet to self-scratch while African great apes use hands almost exclusively. Gibbons represent an evolutionary divergence between monkeys and great apes and incidental observations at the Gibbon Conservation Center pointed to a difference in self-scratching among the four extant gibbon genera (Hoolock, Nomascus, Symphalangus, and Hylobates). To validate and further explore these preliminary observations, we collected systematic data on self-scratching from 32 gibbons, including nine species and all four genera. To supplement gibbon data, we also collected self-scratching information from 18 great apes (four species), five prosimians (two species), 26 New World Monkeys (nine species) and 20 Old World Monkeys (seven species). All monkeys and some prosimians used both hands and feet to self-scratch, whereas one prosimian species used only feet. All African great apes used hands exclusively (orangutans were an exception displaying occasional foot-use). This appears to represent a fundamental difference between monkeys and great apes in limb use. Interestingly, there was a clear difference in self-scratching between the four gibbon genera. Hylobates and Symphalangus self-scratched only with hands (like all African great apes), while Hoolock and Nomascus self-scratched with both hands and feet (like monkeys and prosimians). This difference in gibbon behavior may reflect the evolutionary history of gibbons as Hoolock and Nomascus are thought to have evolved before both Hylobates and Symphalangus. What evolutionary pressures led to this divergent pattern is currently opaque; however, this shift in limb preference may result from niche separation across the order facilitating differences in the behavioral repertoire associated with hind and forelimbs.

An Alu-based phylogeny of gibbons (hylobatidae)

Gibbons (Hylobatidae) are small, arboreal apes indigenous to Southeast Asia that diverged from other apes ∼15-18 Ma. Extant lineages radiated rapidly 6-10 Ma and are organized into four genera (Hylobates, Hoolock, Symphalangus, and Nomascus) consisting of 12-19 species. The use of short interspersed elements (SINEs) as phylogenetic markers has seen recent popularity due to several desirable characteristics: the ancestral state of a locus is known to be the absence of an element, rare potentially homoplasious events are relatively easy to resolve, and samples can be quickly and inexpensively genotyped. During radiation of primates, one particular family of SINEs, the Alu family, has proliferated in primate genomes. Nomascus leucogenys (northern white-cheeked gibbon) sequences were analyzed for repetitive content with RepeatMasker using a custom library. The sequences containing Alu elements identified as members of a gibbon-specific subfamily were then compared with orthologous positions in other primate genomes. A primate phylogenetic panel consisting of 18 primate species, including 13 gibbon species representing all four extant genera, was assayed for all loci, and a total of 125 gibbon-specific Alu insertions were identified. The resulting amplification patterns were used to generate a phylogenetic tree. We demonstrate significant support for Symphalangus as the most basal lineage within the family. Our findings also place Nomascus as a derived lineage, sister to Hoolock, with the Nomascus-Hoolock clade sister to Hylobates. Further, our analysis groups N. leucogenys and Nomascus siki as sister taxa to the exclusion of the other Nomascus species assayed. This study represents the first use of SINEs to determine the genus level phylogenetic relationships within the family Hylobatidae. These relationships have been resolved with robust support at most internal nodes, demonstrating the utility of SINE-based phylogenetic analysis. We postulate that hybridization and rapid radiation may have contributed to the complex and contradictory findings of the previous studies. Our findings will aid in the conservation of these threatened primates and inform future studies of the biogeographical history and distribution of modern gibbon species.

Gene conversion and purifying selection shape nucleotide variation in gibbon L/M opsin genes

Background

Routine trichromatic color vision is a characteristic feature of catarrhines (humans, apes and Old World monkeys). This is enabled by L and M opsin genes arrayed on the X chromosome and an autosomal S opsin gene. In non-human catarrhines, genetic variation affecting the color vision phenotype is reported to be absent or rare in both L and M opsin genes, despite the suggestion that gene conversion has homogenized the two genes. However, nucleotide variation of both introns and exons among catarrhines has only been examined in detail for the L opsin gene of humans and chimpanzees. In the present study, we examined the nucleotide variation of gibbon (Catarrhini, Hylobatidae) L and M opsin genes. Specifically, we focused on the 3.6~3.9-kb region that encompasses the centrally located exon 3 through exon 5, which encode the amino acid sites functional for the spectral tuning of the genes.

Results

Among 152 individuals representing three genera (Hylobates, Nomascus and Symphalangus), all had both L and M opsin genes and no L/M hybrid genes. Among 94 individuals subjected to the detailed DNA sequencing, the nucleotide divergence between L and M opsin genes in the exons was significantly higher than the divergence in introns in each species. The ratio of the inter-LM divergence to the intra-L/M polymorphism was significantly lower in the introns than that in synonymous sites. When we reconstructed the phylogenetic tree using the exon sequences, the L/M gene duplication was placed in the common ancestor of catarrhines, whereas when intron sequences were used, the gene duplications appeared multiple times in different species. Using the GENECONV program, we also detected that tracts of gene conversions between L and M opsin genes occurred mostly within the intron regions.

Conclusions

These results indicate the historical accumulation of gene conversions between L and M opsin genes in the introns in gibbons. Our study provides further support for the homogenizing role of gene conversion between the L and M opsin genes and for the purifying selection against such homogenization in the central exons to maintain the spectral difference between L and M opsins in non-human catarrhines.

Anatomical Contributions to Hylobatid Taxonomy and Adaptation

Compared with the great apes, the small-bodied hylobatids were treated historically as a relatively uniform group with 2 genera, Hylobates and the larger-bodied Symphalangus. Four genera are now recognized, each with a different chromosome number: Hoolock (hoolock) (38), Hylobates (44), Nomascus (crested gibbon) (52), and Symphalangus (siamang) (50). Previous morphological studies based on relative bone lengths, e.g., intermembral indices; molar tooth sizes; and body masses did not distinguish the 4 genera from each other. We applied quantitative anatomical methods to test the hypothesis that each genus can be differentiated from the others using the relative distribution of body mass to the forelimbs and hind limbs. Based on dissections of 13 hylobatids from captive facilities, our findings demonstrate that each of the 4 genera has a distinct pattern of body mass distribution. For example, the adult Hoolock has limb proportions of nearly equal mass, a pattern that differentiates it from species in the genus Hylobates, e.g., H. lar (lar gibbon), H. moloch (Javan gibbon), H. pileatus (pileated gibbon), Nomascus, and Symphalangus. Hylobates is distinct in having heavy hind limbs. Although Symphalangus has been treated as a scaled up version of Hylobates, its forelimb exceeds its hind limb mass, an unusual primate pattern otherwise found only in orangutans. This research provides new information on whole body anatomy and adds to the genetic, ecological, and behavioral evidence for clarifying the taxonomy of the hylobatids. The research also underscores the important contribution of studies on rare species in captivity.

A molecular phylogeny of living primates

Comparative genomic analyses of primates offer considerable potential to define and understand the processes that mold, shape, and transform the human genome. However, primate taxonomy is both complex and controversial, with marginal unifying consensus of the evolutionary hierarchy of extant primate species. Here we provide new genomic sequence (∼8 Mb) from 186 primates representing 61 (∼90%) of the described genera, and we include outgroup species from Dermoptera, Scandentia, and Lagomorpha. The resultant phylogeny is exceptionally robust and illuminates events in primate evolution from ancient to recent, clarifying numerous taxonomic controversies and providing new data on human evolution. Ongoing speciation, reticulate evolution, ancient relic lineages, unequal rates of evolution, and disparate distributions of insertions/deletions among the reconstructed primate lineages are uncovered. Our resolution of the primate phylogeny provides an essential evolutionary framework with far-reaching applications including: human selection and adaptation, global emergence of zoonotic diseases, mammalian comparative genomics, primate taxonomy, and conservation of endangered species.

Advances in human biomedicine, including those focused on changes in genes triggered or disrupted in development, resistance/susceptibility to infectious disease, cancers, mechanisms of recombination, and genome plasticity, cannot be adequately interpreted in the absence of a precise evolutionary context or hierarchy. However, little is known about the genomes of other primate species, a situation exacerbated by a paucity of nuclear molecular sequence data necessary to resolve the complexities of primate divergence over time. We overcome this deficiency by sequencing 54 nuclear gene regions from DNA samples representing ∼90% of the diversity present in living primates. We conduct a phylogenetic analysis to determine the origin, evolution, patterns of speciation, and unique features in genome divergence among primate lineages. The resultant phylogenetic tree is remarkably robust and unambiguously resolves many long-standing issues in primate taxonomy. Our data provide a strong foundation for illuminating those genomic differences that are uniquely human and provide new insights on the breadth and richness of gene evolution across all primate lineages.

Phylogeny and distribution of crested gibbons (genus Nomascus) based on mitochondrial cytochrome b gene sequence data

Crested gibbons, genus Nomascus, are endemic to the Indochinese bioregion and occur only in Vietnam, Laos, Cambodia, and southern China. However, knowledge about the number of species to be recognized and their exact geographical distributions is still limited. To further elucidate the evolutionary history of crested gibbon species and to settle their distribution ranges, we analyzed the complete mitochondrial cytochrome b gene from 79 crested gibbon individuals from known locations. Based on our findings, crested gibbons should be classified into seven species. Within N. concolor, we recognize two subspecies, N. concolor concolor and N. concolor lu. Phylogenetic reconstructions indicate that the northernmost species, N. hainanus, N. nasutus, and N. concolor branched off first, suggesting that the genus originated in the north and successively migrated to the south. The most recent splits within Nomascus occurred between N. leucogenys and N. siki, and between Nomascus sp. and N. gabriellae. Based on our data, the currently postulated distributions of the latter four species have to be revised. Our study shows that molecular methods are a useful tool to elucidate phylogenetic relationships among crested gibbons and to determine species boundaries.

Mitochondrial genome sequences effectively reveal the phylogeny of Hylobates gibbons

Background

Uniquely among hominoids, gibbons exist as multiple geographically contiguous taxa exhibiting distinctive behavioral, morphological, and karyotypic characteristics. However, our understanding of the evolutionary relationships of the various gibbons, especially among Hylobates species, is still limited because previous studies used limited taxon sampling or short mitochondrial DNA (mtDNA) sequences. Here we use mtDNA genome sequences to reconstruct gibbon phylogenetic relationships and reveal the pattern and timing of divergence events in gibbon evolutionary history.

Methodology/Principal Findings

We sequenced the mitochondrial genomes of 51 individuals representing 11 species belonging to three genera (Hylobates, Nomascus and Symphalangus) using the high-throughput 454 sequencing system with the parallel tagged sequencing approach. Three phylogenetic analyses (maximum likelihood, Bayesian analysis and neighbor-joining) depicted the gibbon phylogenetic relationships congruently and with strong support values. Most notably, we recover a well-supported phylogeny of the Hylobates gibbons. The estimation of divergence times using Bayesian analysis with relaxed clock model suggests a much more rapid speciation process in Hylobates than in Nomascus.

Conclusions/Significance

Use of more than 15 kb sequences of the mitochondrial genome provided more informative and robust data than previous studies of short mitochondrial segments (e.g., control region or cytochrome b) as shown by the reliable reconstruction of divergence patterns among Hylobates gibbons. Moreover, molecular dating of the mitogenomic divergence times implied that biogeographic change during the last five million years may be a factor promoting the speciation of Sundaland animals, including Hylobates species.

Unresolved molecular phylogenies of gibbons and siamangs (Family: Hylobatidae) based on mitochondrial, Y-linked, and X-linked loci indicate a rapid Miocene radiation or sudden vicariance event

According to recent taxonomic reclassification, the primate family Hylobatidae contains four genera (Hoolock, Nomascus, Symphalangus, and Hylobates) and between 14 and 18 species, making it by far the most species-rich group of extant hominoids. Known as the "small apes", these small arboreal primates are distributed throughout Southeast, South and East Asia. Considerable uncertainty surrounds the phylogeny of extant hylobatids, particularly the relationships among the genera and the species within the Hylobates genus. In this paper we use parsimony, likelihood, and Bayesian methods to analyze a dataset containing nearly 14 kilobase pairs, which includes newly collected sequences from X-linked, Y-linked, and mitochondrial loci together with data from previous mitochondrial studies. Parsimony, likelihood, and Bayesian analyses largely failed to find a significant difference among phylogenies with any of the four genera as the most basal taxon. All analyses, however, support a tree with Hylobates and Symphalangus as most closely related genera. One strongly supported phylogenetic result within the Hylobates genus is that Hylobates pileatus is the most basal taxon. Multiple analyses failed to find significant support for any singular genus-level phylogeny. While it is natural to suspect that there might not be sufficient data for phylogenetic resolution (whenever that situation occurs), an alternative hypothesis relating to the nature of gibbon speciation exists. This lack of resolution may be the result of a rapid radiation or a sudden vicariance event of the hylobatid genera, and it is likely that a similarly rapid radiation occurred within the Hylobates genus. Additional molecular and paleontological evidence are necessary to better test among these, and other, hypotheses of hylobatid evolution.

A comparative analysis of internal cranial anatomy in the hylobatidae

Craniometric studies on the hylobatids using external metrics (Creel and Preuschoft, 1976, 1984) sorted hylobatid populations into primary species groupings which are in accordance with the four currently recognized generic-level groupings. The goal of the current study was to assess the relative orientations of the orbits, palate, and basioccipital clivus among the hylobatid genera in an effort to further clarify whether the lesser apes differ significantly in these internal cranial features and how that variation patterns across the groups. Nine angular variables quantifying orbital, palatal, and basioccipital clivus orientations were measured on lateral view radiographs of adults representing three of the four hylobatid genera: Hylobates; Nomascus; and, Symphalangus. The interspecific adult hylobatid means for the angular variables were analyzed using t-test contrasts. The total sample was further subjected to discriminant function analysis (DFA) to test for the ability of craniofacial angular variables to distinguish the hylobatid genera from one another. The three hylobatid genera displayed significant morphological differentiation in orbital, palatal, and posterior skull base orientations. Normal, jackknifed, and cross-validation DFA procedures correctly identified the hylobatids 50-100% of the time. The observed morphological patterns generally mapped onto the findings of earlier external craniometric hylobatid studies and suggest concordance between specific internal and external cranial features. This article is the first comprehensive study of variation in internal cranial anatomy of the Hylobatidae and includes the first published craniofacial angular data for Nomascus.

Mitochondrial evidence for multiple radiations in the evolutionary history of small apes

BACKGROUND

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.

RESULTS

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 approximately 6.7 and approximately 8.3 mya, the four gibbon genera diverged from each other. In the second (approximately 3.0 - approximately 3.9 mya) and in the third period (approximately 1.3 - approximately 1.8 mya), Hylobates and Hoolock differentiated. Finally, between approximately 0.5 and approximately 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 approximately 4.2 until approximately 0.4 mya.

CONCLUSIONS

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.

A small, variable, and irregular killer cell Ig-like receptor locus accompanies the absence of MHC-C and MHC-G in gibbons

The killer cell Ig-like receptors (KIRs) of NK cells recognize MHC class I ligands and function in placental reproduction and immune defense against pathogens. During the evolution of monkeys, great apes, and humans, an ancestral KIR3DL gene expanded to become a diverse and rapidly evolving gene family of four KIR lineages. Characterizing the KIR locus are three framework regions, defining two intervals of variable gene content. By analysis of four KIR haplotypes from two species of gibbon, we find that the smaller apes do not conform to these rules. Although diverse and irregular in structure, the gibbon haplotypes are unusually small, containing only two to five functional genes. Comparison with the predicted ancestral hominoid KIR haplotype indicates that modern gibbon KIR haplotypes were formed by a series of deletion events, which created new hybrid genes as well as eliminating ancestral genes. Of the three framework regions, only KIR3DL3 (lineage V), defining the 5' end of the KIR locus, is present and intact on all gibbon KIR haplotypes. KIR2DL4 (lineage I) defining the central framework region has been a major target for elimination or inactivation, correlating with the absence of its putative ligand, MHC-G, in gibbons. Similarly, the MHC-C-driven expansion of lineage III KIR genes in great apes has not occurred in gibbons because they lack MHC-C. Our results indicate that the selective forces shaping the size and organization of the gibbon KIR locus differed from those acting upon the KIR of other hominoid species.

Visual preference in a human-reared agile gibbon (Hylobates agilis)

Visual preference was evaluated in a male agile gibbon. The subject was raised by humans immediately after birth, but lived with his biological family from one year of age. Visual preference was assessed using a free-choice task in which five or six photographs of different primate species, including humans, were presented on a touch-sensitive screen. The subject touched one of them. Food rewards were delivered irrespective of the subject's responses. We prepared two types of stimulus sets. With set 1, the subject touched photographs of humans more frequently than those of other species, recalling previous findings in human-reared chimpanzees. With set 2, photographs of nine species of gibbons were presented. Chimpanzees touched photographs of white-handed gibbons more than those of other gibbon species. The gibbon subject initially touched photographs of agile gibbons more than white-handed gibbons, but after one and two years his choice patterns resembled the chimpanzees'. The results suggest that, as in chimpanzees, visual preferences of agile gibbons are not genetically programmed but develop through social experience during infancy.

Estimating the phylogeny and divergence times of primates using a supermatrix approach

Background

The primates are among the most broadly studied mammalian orders, with the published literature containing extensive analyses of their behavior, physiology, genetics and ecology. The importance of this group in medical and biological research is well appreciated, and explains the numerous molecular phylogenies that have been proposed for most primate families and genera. Composite estimates for the entire order have been infrequently attempted, with the last phylogenetic reconstruction spanning the full range of primate evolutionary relationships having been conducted over a decade ago.

Results

To estimate the structure and tempo of primate evolutionary history, we employed Bayesian phylogenetic methods to analyze data supermatrices comprising 7 mitochondrial genes (6,138 nucleotides) from 219 species across 67 genera and 3 nuclear genes (2,157 nucleotides) from 26 genera. Many taxa were only partially represented, with an average of 3.95 and 5.43 mitochondrial genes per species and per genus, respectively, and 2.23 nuclear genes per genus. Our analyses of mitochondrial DNA place Tarsiiformes as the sister group of Strepsirrhini. Within Haplorrhini, we find support for the primary divergence of Pitheciidae in Platyrrhini, and our results suggest a sister grouping of African and non-African colobines within Colobinae and of Cercopithecini and Papionini within Cercopthecinae. Date estimates for nodes within each family and genus are presented, with estimates for key splits including: Strepsirrhini-Haplorrhini 64 million years ago (MYA), Lemuriformes-Lorisiformes 52 MYA, Platyrrhini-Catarrhini 43 MYA and Cercopithecoidea-Hominoidea 29 MYA.

Conclusion

We present an up-to-date, comprehensive estimate of the structure and tempo of primate evolutionary history. Although considerable gaps remain in our knowledge of the primate phylogeny, increased data sampling, particularly from nuclear loci, will be able to provide further resolution.

A chromosomal inversion unique to the northern white-cheeked gibbon

The gibbon family belongs to the superfamily Hominoidea and includes 15 species divided into four genera. Each genus possesses a distinct karyotype with chromosome numbers varying from 38 to 52. This diversity is the result of numerous chromosomal changes that have accumulated during the evolution of the gibbon lineage, a quite unique feature in comparison with other hominoids and most of the other primates. Some gibbon species and subspecies rank among the most endangered primates in the world. Breeding programs can be extremely challenging and hybridization plays an important role within the factors responsible for the decline of captive gibbons. With less than 500 individuals left in the wild, the northern white-cheeked gibbon (Nomascus leucogenys leucogenys, NLE) is the most endangered primate in a successful captive breeding program. We present here the analysis of an inversion that we show being specific for the northern white-cheeked gibbon and can be used as one of the criteria to distinguish this subspecies from other gibbon taxa. The availability of the sequence spanning for one of the breakpoints of the inversion allows detecting it by a simple PCR test also on low quality DNA. Our results demonstrate the important role of genomics in providing tools for conservation efforts.

Tracking the complex flow of chromosome rearrangements from the Hominoidea Ancestor to extant Hylobates and Nomascus Gibbons by high-resolution synteny mapping

In this study we characterized the extension, reciprocal arrangement, and orientation of syntenic chromosomal segments in the lar gibbon (Hylobates lar, HLA) by hybridization of a panel of approximately 1000 human BAC clones. Each lar gibbon rearrangement was defined by a splitting BAC clone or by two overlapping clones flanking the breakpoint. A reconstruction of the synteny arrangement of the last common ancestor of all living lesser apes was made by combining these data with previous results in Nomascus leucogenys, Hoolock hoolock, and Symphalangus syndactylus. The definition of the ancestral synteny organization facilitated tracking the cascade of chromosomal changes from the Hominoidea ancestor to the present day karyotype of Hylobates and Nomascus. Each chromosomal rearrangement could be placed within an approximate phylogenetic and temporal framework. We identified 12 lar-specific rearrangements and five previously undescribed rearrangements that occurred in the Hylobatidae ancestor. The majority of the chromosomal differences between lar gibbons and humans are due to rearrangements that occurred in the Hylobatidae ancestor (38 events), consistent with the hypothesis that the genus Hylobates is the most recently evolved lesser ape genus. The rates of rearrangements in gibbons are 10 to 20 times higher than the mammalian default rate. Segmental duplication may be a driving force in gibbon chromosome evolution, because a consistent number of rearrangements involves pericentromeric regions (10 events) and centromere inactivation (seven events). Both phenomena can be reasonably supposed to have strongly contributed to the euchromatic dispersal of segmental duplications typical of pericentromeric regions. This hypothesis can be more fully tested when the sequence of this gibbon species becomes available. The detailed synteny map provided here will, in turn, substantially facilitate sequence assembly efforts.

Patterns of dispersal in Sumatran siamangs (Symphalangus syndactylus): preliminary mtDNA evidence suggests more frequent male than female dispersal to adjacent groups

Gibbons of both sexes have been observed emigrating from their natal groups, but the consequences of dispersal in gibbons are poorly understood, and it is unclear whether these are the same for both sexes. I sequenced a 350-bp fragment of mitochondrial DNA from 18 adults in seven siamang (Symphalangus syndactylus) groups at the Way Canguk Research Station in southern Sumatra to assess patterns of matrilineal relatedness among and within siamang groups, and to assess their fit with different patterns of sex-specific dispersal. A total of 11 haplotypes were identified in the seven study groups; 50% of adult males in five contiguous groups shared a haplotype with a member of an immediately adjacent group, whereas only 16.7% of females shared a haplotype with a neighbor. The apparent difference persisted if only same-sex individuals were considered (37.5% of males vs. 0% of females). Four of the seven study groups contained two adult males and a single adult female. In three multimale groups, the three adults all had different haplotypes, suggesting that neither male was the retained adult offspring of the female, whereas in the fourth group, the haplotype of one male was identical with that of the adult female. The high diversity of haplotypes and the absence of clustering among female haplotypes in the study neighborhood suggest that female dispersal to territories adjacent to the natal group may be relatively rare. The presence of some clustering of male haplotypes suggests that shorter dispersal distances may be more common in males.

Mitochondrial DNA hypervariable region-1 sequence variation and phylogeny of the concolor gibbons, Nomascus

The still little known concolor gibbons are represented by 14 taxa (five species, nine subspecies) distributed parapatrically in China, Myanmar, Vietnam, Laos and Cambodia. To set the stage for a phylogeographic study of the genus we examined DNA sequences from the highly variable mitochondrial hypervariable region-1 (HVR-1 or control region) in 51 animals, mostly of unknown geographic provenance. We developed gibbon-specific primers to amplify mtDNA noninvasively and obtained >477 bp sequences from 38 gibbons in North American and European zoos and >159 bp sequences from ten Chinese museum skins. In hindsight, we believe these animals represent eight of the nine nominal subspecies and four of the five nominal species. Bayesian, maximum likelihood and maximum parsimony haplotype network analyses gave concordant results and show Nomascus to be monophyletic. Significant intraspecific variation within N. leucogenys (17 haplotypes) is comparable with that reported earlier in Hylobates lar and less than half the known interspecific pairwise distances in gibbons. Sequence data support the recognition of five species (concolor, leucogenys, nasutus, gabriellae and probably hainanus) and suggest that nasutus is the oldest and leucogenys, the youngest taxon. In contrast, the subspecies N. c. furvogaster, N. c. jingdongensis, and N. leucogenys siki, are not recognizable at this otherwise informative genetic locus. These results show that HVR-1 sequence is variable enough to define evolutionarily significant units in Nomascus and, if coupled with multilocus microsatellite or SNP genotyping, more than adequate to characterize their phylogeographic history. There is an urgent need to obtain DNA from gibbons of known geographic provenance before they are extirpated to facilitate the conservation genetic management of the surviving animals.

Resolution of the Hylobates phylogeny: congruence of mitochondrial D-loop sequences with molecular, behavioral, and morphological data sets

Gibbons of the genus Hylobates likely speciated very rapidly following isolation by rising sea levels during the Pleistocene. We sequenced the hypervariable region I (HV-I) of the mitochondrial D-loop to reconstruct the phylogeny of this group. Although the results clearly supported monophyly of each of the six species, the relationships among them were not clearly resolved by these data alone. A homogeneity test against published data sets of a coding mitochondrial locus (ND3-ND4 region), behavioral characters (vocalizations), and morphological traits (including skeletal and soft tissue anatomy) revealed no significant incongruence, and combining them resulted in a phylogenetic tree with much stronger support. The Kloss's gibbon (H. klossii), long considered a primitive taxon based on morphology, shares many molecular and vocal characteristics with the Javan gibbon (H. moloch), and appear as the most recently derived species. The northernmost species (H. lar and H. pileatus) are the most basal taxa. These data suggest that ancestral gibbons radiated from north to south. Unlike other markers, the HV-I region can accurately identify members of different gibbon species much like a DNA barcode, with potential applications to conservation.

A most distant intergeneric hybrid offspring (Larcon) of lesser apes, Nomascus leucogenys and Hylobates lar

Unlike humans, which are the sole remaining representatives of a once larger group of bipedal apes (hominins), the "lesser apes" (hylobatids) are a diverse radiation with numerous extant species. Consequently, the lesser apes can provide a valuable evolutionary window onto the possible interactions (e.g., interbreeding) of hominin lineages coexisting in the same time and place. In the present work, we employ chromosomal analyses to verify the hybrid ancestry of an individual (Larcon) produced by two of the most distant genera of lesser apes, Hylobates (lar-group gibbons) and Nomascus (concolor-group gibbons). In addition to a mixed pelage pattern, the hybrid animal carries a 48-chromosome karyotype that consists of the haploid complements of each parental species: Hylobates lar (n = 22) and Nomascus leucogenys leucogenys (n = 26). Studies of this animal's karyotype shed light onto the processes of speciation and genus-level divergence in the lesser apes and, by extension, across the Hominoidea.