Search

Nomascus leucogenys is a critically endangered species of small apes. Here, we sequenced and assembled the male genome of N. leucogenys, using PacBio and Hi-C datasets, with a particular focus on its Y-chromosome. The resulting high-quality haplotype-phased assemblies are at chromosome-scale, with scaffold/contig N50 values of 124.2/102.2 Mb for Haplotype 1 and 121.2/85.67 Mb for Haplotype 2. The assembled Y-chromosome spans 16.06 Mb. BUSCO assessment indicated completeness scores exceeding 95%. We predicted 18,925 protein-coding genes (23,783 mRNAs), including 58 genes on the Y-chromosome. Approximately 50% of the genome comprises repetitive elements. These comprehensive genome datasets will serve as a valuable resource for future studies on the genetics and protection of gibbons and improve our understanding on the evolution of Y-chromosome-related genes in primates., Competing Interests: Competing interests: The authors declare no competing interests., (© 2024. The Author(s).)

Human and non-human primates have strikingly similar genomes, but they strongly differ in many brain-based processes (e.g., behaviour and cognition). While the functions of protein-coding genes have been extensively studied, rather little is known about the role of non-coding RNAs such as long non-coding RNAs (lncRNAs). Here, we predicted lncRNAs and analysed their expression pattern across different brain regions of human and non-human primates (chimpanzee, gorilla, and gibbon). Our analysis identified shared orthologous and non-orthologous lncRNAs, showing striking differences in the genomic features. Differential expression analysis of the shared orthologous lncRNAs from humans and chimpanzees revealed distinct expression patterns in subcortical regions (striatum, hippocampus) and neocortical areas while retaining a homogeneous expression in the cerebellum. Co-expression analysis of lncRNAs and protein-coding genes revealed massive proportions of co-expressed pairs in neocortical regions of humans compared to chimpanzees. Network analysis of co-expressed pairs revealed the distinctive role of the hub-acting orthologous lncRNAs in a region- and species-specific manner. Overall, our study provides novel insight into lncRNA driven gene regulatory landscape, neural regulation, brain evolution, and constitutes a resource for primate's brain lncRNAs., (© 2024. The Author(s).)

A- and B-antigens are present on red blood cells (RBCs) as well as other cells and secretions in Hominoidea including humans and apes such as chimpanzees and gibbons, whereas expression of these antigens on RBCs is subtle in monkeys such as Japanese macaques. Previous studies have indicated that H-antigen expression has not completely developed on RBCs in monkeys. Such antigen expression requires the presence of H-antigen and A- or B-transferase expression in cells of erythroid lineage, although whether or not ABO gene regulation is associated with the difference of A- or B-antigen expression between Hominoidea and monkeys has not been examined. Since it has been suggested that ABO expression on human erythrocytes is dependent upon an erythroid cell-specific regulatory region or the + 5.8-kb site in intron 1, we compared the sequences of ABO intron 1 among non-human primates, and demonstrated the presence of sites orthologous to the + 5.8-kb site in chimpanzees and gibbons, and their absence in Japanese macaques. In addition, luciferase assays revealed that the former orthologues enhanced promoter activity, whereas the corresponding site in the latter did not. These results suggested that the A- or B-antigens on RBCs might be ascribed to emergence of the + 5.8-kb site or the corresponding regions in ABO through genetic evolution., (© 2023. The Author(s).)

The Javan gibbon, Hylobates moloch, is an endangered gibbon species restricted to the forest remnants of western and central Java, Indonesia, and one of the rarest of the Hylobatidae family. Hylobatids consist of 4 genera (Holoock, Hylobates, Symphalangus, and Nomascus) that are characterized by different numbers of chromosomes, ranging from 38 to 52. The underlying cause of this karyotype plasticity is not entirely understood, at least in part, due to the limited availability of genomic data. Here we present the first scaffold-level assembly for H. moloch using a combination of whole-genome Illumina short reads, 10X Chromium linked reads, PacBio, and Oxford Nanopore long reads and proximity-ligation data. This Hylobates genome represents a valuable new resource for comparative genomics studies in primates., (© The American Genetic Association. 2022.)

The growth hormone (GH) locus has experienced a dramatic evolution in primates, becoming multigenic and diverse in anthropoids. Despite sequence information from a vast number of primate species, it has remained unclear how the multigene family was favored. We compared the structure and composition of apes' GH loci as a prerequisite to understanding their origin and possible evolutionary role. These thorough analyses of the GH loci of the chimpanzee, gorilla, and orangutan were done by resorting to previously sequenced bacterial artificial chromosomes (BACs) harboring them, as well as to their respective genome projects data available in GenBank. The GH loci of modern man, Neanderthal, gibbon, and wild boar were retrieved from GenBank. Coding regions, regulatory elements, and repetitive sequences were identified and compared among species. The GH loci of all the analyzed species are flanked by the genes CD79B (5') and ICAM-1 (3'). In man, Neanderthal, and chimpanzee, the loci were integrated by five almost indistinguishable genes; however, in the former two, they rendered three different hormones, and in the latter, four different proteins were derived. Gorilla exhibited six genes, gibbon seven, and orangutan four. The sequences of the proximal promoters, enhancers, P-elements, and a locus control region (LCR) were highly conserved. The locus evolution might have implicated duplications of the ancestral pituitary gene ( GH-N ) and subsequent diversification of the copies, leading to the placental single GH-V gene and the multiple CSH genes.

Gibbons are the most speciose family of living apes, characterized by a diverse chromosome number and rapid rate of large-scale rearrangements. Here we performed single-cell template strand sequencing (Strand-seq), molecular cytogenetics, and deep in silico analysis of a southern white-cheeked gibbon genome, providing the first comprehensive map of 238 previously hidden small-scale inversions. We determined that more than half are gibbon specific, at least fivefold higher than shown for other primate lineage-specific inversions, with a significantly high number of small heterozygous inversions, suggesting that accelerated evolution of inversions may have played a role in the high sympatric diversity of gibbons. Although the precise mechanisms underlying these inversions are not yet understood, it is clear that segmental duplication-mediated NAHR only accounts for a small fraction of events. Several genomic features, including gene density and repeat (e.g., LINE-1) content, might render these regions more break-prone and susceptible to inversion formation. In the attempt to characterize interspecific variation between southern and northern white-cheeked gibbons, we identify several large assembly errors in the current GGSC Nleu3.0/nomLeu3 reference genome comprising more than 49 megabases of DNA. Finally, we provide a list of 182 candidate genes potentially involved in gibbon diversification and speciation., (© 2022 Mercuri et al.; Published by Cold Spring Harbor Laboratory Press.)

We expand upon prior work (Gibbons et al., ) relating childhood stressor effects, particularly harsh childhood environments, to risky behavior and ultimately physical health by adding longer-term outcomes - deoxyribonucleic acid (DNA) methylation-based measures of accelerated aging (DNA m -aging). Further, following work on the effects of early exposure to danger (McLaughlin et al., ), we also identify an additional pathway from harsh childhood environments to DNA m -aging that we label the danger/FKBP5 pathway, which includes early exposure to dangerous community conditions that are thought to impact glucocorticoid regulation and pro-inflammatory mechanisms. Because different DNA m -aging indices provide different windows on accelerated aging, we contrast effects on early indices of DNA m -aging based on chronological age with later indices that focused on predicting biological outcomes. We utilize data from Family and Community Health Study participants ( N = 449) from age 10 to 29. We find that harshness influences parenting, which, in turn, influences accelerated DNA m -aging through the risky cognitions and substance use (i.e., behavioral) pathway outlined by Gibbons et al. (). Harshness is also associated with increased exposure to threat/danger, which, in turn, leads to accelerated DNA m -aging through effects on FKBP5 activity and enhanced pro-inflammatory tendencies (i.e., the danger/FKBP5 pathway).

Gibbons (Family Hylobatidae) are a suitable model for exploring hybridization in pair-living primates as several species form hybrid zones. In Khao Yai National Park, Thailand, white-handed gibbons (Hylobates lar) and pileated gibbons (Hylobates pileatus) are distributed parapatrically and hybridize in a narrow zone. Their phenotypic characteristics suggest limited inter-species gene flow, although this has never been assessed. To uncover the history and degree of gene flow between the two species, we studied the genetic structure of gibbons in the hybrid zone by analyzing fecal DNA samples, phenotypic characteristics, vocalizations and individuals' social status. We determined eight autosomal single nucleotide variant (SNV) loci, and mitochondrial DNA (mtDNA) and Y-chromosomal haplotypes of 72 gibbons. We compared these markers with reference types of wild pureblood white-handed gibbons (n = 12) in Kaeng Krachan National Park and pureblood pileated gibbons (n = 4) in Khao Soi Dao Wildlife Sanctuary. Autosomal genotypic analyses confirmed the various levels of mixed ancestry for several adult gibbons with or without atypical phenotypic traits in Khao Yai National Park. In some other adult gibbons, the mixed ancestry was not detected in either autosomal SNVs or their phenotypic traits but the mtDNA. Both male and female adult hybrids formed reproductive units mainly with a phenotypic pureblood partner and many of them produced offspring. Taken together, our results suggest that once hybridization occurs, white-handed-pileated-gibbon hybrids can reproduce with either parental species and that the backcrossing and thus introgression may occur in successive generations, with no drastic changes in phenotypic appearance., Competing Interests: The authors have declared that no competing interests exist.

Natural hybridization has played various roles in the evolutionary history of primates. Its consequences range from genetic introgression between taxa, formation of hybrid zones, and formation of new lineages. Hylobates lar, the white-handed gibbon, and Hylobates pileatus, the pileated gibbon, are largely allopatric species in Southeast Asia with a narrow contact zone in Khao Yai National Park, Thailand, which contains both parental types and hybrids. Hybrid individuals in the zone are recognizable by their intermediate pelage and vocal patterns, but have not been analyzed genetically. We analyzed mitochondrial and microsatellite DNA of 52 individuals to estimate the relative genetic contributions of the parental species to each individual, and the amount of introgression into the parental species. We obtained fecal samples from 33 H. lar, 15 H. pileatus and four phenotypically intermediate individuals in the contact zone. Both mitochondrial and microsatellite markers confirmed distinct differences between these taxa. Both H. lar and H. pileatus contributed to the maternal lineages of the hybrids based on mitochondrial analysis; hybrids were viable and present in socially normal reproductive pairs. The microsatellite analysis identified ten admixed individuals, four F1 hybrids, which corresponded to phenotypic hybrids, and six H. lar-like backcrosses. All 15 H. pileatus samples were identified as originating from genetically H. pileatus individuals with no H. lar admixture; hence, backcrossing is biased toward H. lar. A relatively low number of phenotypic hybrids and backcrossed individuals along with a high number of parental types indicates a bimodal hybrid zone, which suggests relatively strong bias in mate selection between the species., (© 2021. Japan Monkey Centre.)

Gibbons of the genus Hylobates, which inhabit Southeast Asia, show great diversity and comprise seven to nine species. Natural hybridisation has been observed in several species contact zones, but the history and extent of hybridisation and introgression in possibly historical and the current contact zones remain unclear. To uncover Hylobates species phylogeny and the extent of introgression in their evolution, genotyping by random amplicon sequencing-direct (GRAS-Di) was applied to 47 gibbons, representing seven Hylobates species/subspecies and two outgroup gibbon species. Over 200,000 autosomal single-nucleotide variant sites were identified. The autosomal phylogeny supported that divergence from the mainland species began ~3.5 million years ago, and subsequently occurred among the Sundaic island species. Significant introgression signals were detected between H. lar and H. pileatus, H. lar and H. agilis and H. albibarbis and H. muelleri, which all are parapatric and form ongoing hybrid zones. Furthermore, the introgression signals were detected in every analysed individual of these species, indicating a relatively long history of hybridisation, which might have affected the entire gene pool. By contrast, signals of introgression were either not detected or doubtful in other species pairs living on different islands, indicating the rarity of hybridisation and introgression, even though the Sundaic islands were connected during the Pliocene and Pleistocene glacial events., (© 2021. The Author(s), under exclusive licence to The Genetics Society.)

Centromeres are functionally conserved chromosomal loci essential for proper chromosome segregation during cell division, yet they show high sequence diversity across species. Despite their variation, a near universal feature of centromeres is the presence of repetitive sequences, such as DNA satellites and transposable elements (TEs). Because of their rapidly evolving karyotypes, gibbons represent a compelling model to investigate divergence of functional centromere sequences across short evolutionary timescales. In this study, we use ChIP-seq, RNA-seq, and fluorescence in situ hybridization to comprehensively investigate the centromeric repeat content of the four extant gibbon genera (Hoolock, Hylobates, Nomascus, and Siamang). In all gibbon genera, we find that CENP-A nucleosomes and the DNA-proteins that interface with the inner kinetochore preferentially bind retroelements of broad classes rather than satellite DNA. A previously identified gibbon-specific composite retrotransposon, LAVA, known to be expanded within the centromere regions of one gibbon genus (Hoolock), displays centromere- and species-specific sequence differences, potentially as a result of its co-option to a centromeric function. When dissecting centromere satellite composition, we discovered the presence of the retroelement-derived macrosatellite SST1 in multiple centromeres of Hoolock, whereas alpha-satellites represent the predominate satellite in the other genera, further suggesting an independent evolutionary trajectory for Hoolock centromeres. Finally, using de novo assembly of centromere sequences, we determined that transcripts originating from gibbon centromeres recapitulate the species-specific TE composition. Combined, our data reveal dynamic shifts in the repeat content that define gibbon centromeres and coincide with the extensive karyotypic diversity within this lineage., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Co-option of transposable elements (TEs) to become part of existing or new enhancers is an important mechanism for evolution of gene regulation. However, contributions of lineage-specific TE insertions to recent regulatory adaptations remain poorly understood. Gibbons present a suitable model to study these contributions as they have evolved a lineage-specific TE called LAVA (LINE- Alu Sz-VNTR- Alu LIKE ), which is still active in the gibbon genome. The LAVA retrotransposon is thought to have played a role in the emergence of the highly rearranged structure of the gibbon genome by disrupting transcription of cell cycle genes. In this study, we investigated whether LAVA may have also contributed to the evolution of gene regulation by adopting enhancer function. We characterized fixed and polymorphic LAVA insertions across multiple gibbons and found 96 LAVA elements overlapping enhancer chromatin states. Moreover, LAVA was enriched in multiple transcription factor binding motifs, was bound by an important transcription factor (PU.1), and was associated with higher levels of gene expression in cis We found gibbon-specific signatures of purifying/positive selection at 27 LAVA insertions. Two of these insertions were fixed in the gibbon lineage and overlapped with enhancer chromatin states, representing putative co-opted LAVA enhancers. These putative enhancers were located within genes encoding SETD2 and RAD9A, two proteins that facilitate accurate repair of DNA double-strand breaks and prevent chromosomal rearrangement mutations. Co-option of LAVA in these genes may have influenced regulation of processes that preserve genome integrity. Our findings highlight the importance of considering lineage-specific TEs in studying evolution of gene regulatory elements., Competing Interests: The authors declare no competing interest.

A novel gamma-retroviral sequence (7912 bp), inclusive of both partial 5' and 3' long terminal repeat regions, was identified from the brain of a black flying-fox (Pteropus alecto), Queensland, Australia. The sequence was distinct from other retroviral sequences identified in bats and showed greater identity to Koala, Gibbon ape leukaemia, Melomys burtoni and Woolly monkey retroviruses, forming their own phylogenetic clade. This finding suggests that these retroviruses may have an unknown common ancestor and that further investigation into the diversity of gamma-retroviruses in Australian Pteropus species may elucidate their evolutionary origins.

Objectives: Semenogelin 1 and 2 (SEMG1 and SEMG2) are known as semen coagulating proteins in primates with a repetitive structure of 60-amino acids. The number of repeats varies among species and is hypothesized to be related to the level of primate sperm competition. Gibbons until recently were thought to be monogamous primates, but it is now known that gibbon social structure is flexible. Thus, hypotheses of the relationship between the SEMGs evolution and mating systems were tested., Materials and Methods: The sequences of the exon 2 of the SEMG1 and SEMG2 were obtained from 50 captive gibbons comprising six species belonging to three genera (Hylobates, Symphalangus, and Nomascus). Then we quantified the levels of polymorphism and estimated rates of protein evolution by calculating d N /d S ratio., Results: Several mutations that create a premature stop codon in the SEMG1 and a reduction in the repeats of the SEMG2 in the genus Hylobates were observed and may alter the coding properties for these proteins. We also found different level of nucleotide diversity in each gene and between genera. Strikingly, in Nomascus leucogenys we discovered a high d N /d S ratio in the SEMG1 and SEMG2. The Nomascus SEMG2 also showed significantly lower nucleotide diversity than the other two genera., Discussion: These results are consistent with the presence of a strong positive selection in the Nomascus lineage even if the exact selective forces acting on these genes are not yet conclusively known. We were not able to demonstrate, among gibbons, unambiguous relationships between the SEMGs evolution and mating systems., (© 2018 Wiley Periodicals, Inc.)

Retrogenes are a class of functional genes derived from the mRNA of various intron-containing genes. PIPSL was created through a unique mechanism, whereby distinct genes were assembled at the RNA level, and the resulting chimera was then reverse transcribed and integrated into the genome by the L1 retrotransposon. Expression of PIPSL RNA via its transcription start sites (TSSs) has been confirmed in the testes of humans and chimpanzee. Here, we demonstrated that PIPSL RNA is expressed in the testis of the white-handed gibbon. The 5'-end positions of gibbon RNAs were confined to a narrow range upstream of the PIPSL start codon and overlapped with those of orangutan and human, suggesting that PIPSL TSSs are similar among hominoid species. Reporter assays using a luciferase gene and the flanking sequences of human PIPSL showed that an upstream sequence exhibits weak promoter activity in human cells. Our findings suggest that PIPSL might have acquired a promoter at an early stage of hominoid evolution before the divergence of gibbons and ultimately retained similar TSSs in all of the lineages. Moreover, the upstream sequence derived from the phosphatidylinositol-4-phosphate 5-kinase, type I, alpha 5' untranslated region and/or neighboring repetitive sequences in the genome possibly exhibits promoter activity. Furthermore, we observed that a TATA-box-like sequence has emerged by nucleotide substitution in a lineage leading to humans, with this possibly responsible for a broader distribution of the human PIPSL TSSs., (Copyright © 2018 Elsevier B.V. All rights reserved.)

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.

Objectives: It has long been recognized that in gibbons both sexes disperse from the natal group. However, the fate of dispersed individuals was rarely documented. Here we provide the first detailed information on sex differences in dispersal patterns by analyzing the spatial genetic structure of a well-known white-handed gibbon (Hylobates lar) population., Materials and Methods: Mitochondrial DNA (mtDNA) and Y-chromosomal haplotypes, and autosomal microsatellite genotypes were determined for individuals of the Mo Singto study site, Khao Yai National Park, Thailand. Mantel tests for the three genetic marker types were performed for 17 gibbon groups comprising 23 adult males and 18 adult females., Results: Significant positive Mantel correlations were observed for spatial distance and both autosomal microsatellite-based as well as Y-chromosomal haplotype-based genetic distance among adult males. Neighboring adult males tended to be genetically related and share Y-chromosomal haplotypes. Conversely, no significant Mantel correlations were observed either in autosomal microsatellites or mtDNA among adult females., Discussion: Our results confirm, at a genetic level, hypotheses from long-term demographic observations that white-handed gibbon males of the Mo Singto population primarily disperse into adjacent groups. Instead, females disperse more opportunistically either to adjacent or more distant groups. This sex-specific difference reflects an apparent greater tolerance between males than between females. The higher tolerance of adult males allows the formation of stable multimale groups and facilitates male dispersal into an adjacent group. Stable multifemale groups have never been documented for white-handed gibbons probably due to feeding competition between females., (© 2018 Wiley Periodicals, Inc.)

The phylogenetic relationships among extant gibbon species remain unresolved despite numerous efforts using morphological, behavorial, and genetic data and the sequencing of whole genomes. A major challenge in reconstructing the gibbon phylogeny is the radiative speciation process, which resulted in extremely short internal branches in the species phylogeny and extensive incomplete lineage sorting with extensive gene-tree heterogeneity across the genome. Here, we analyze two genomic-scale data sets, with ∼10,000 putative noncoding and exonic loci, respectively, to estimate the species tree for the major groups of gibbons. We used the Bayesian full-likelihood method bpp under the multispecies coalescent model, which naturally accommodates incomplete lineage sorting and uncertainties in the gene trees. For comparison, we included three heuristic coalescent-based methods (mp-est, SVDQuartets, and astral) as well as concatenation. From both data sets, we infer the phylogeny for the four extant gibbon genera to be (Hylobates, (Nomascus, (Hoolock, Symphalangus))). We used simulation guided by the real data to evaluate the accuracy of the methods used. Astral, while not as efficient as bpp, performed well in estimation of the species tree even in presence of excessive incomplete lineage sorting. Concatenation, mp-est and SVDQuartets were unreliable when the species tree contains very short internal branches. Likelihood ratio test of gene flow suggests a small amount of migration from Hylobates moloch to H. pileatus, while cross-genera migration is absent or rare. Our results highlight the utility of coalescent-based methods in addressing challenging species tree problems characterized by short internal branches and rampant gene tree-species tree discordance., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Superfamily Hominoidea, which consists of Hominidae (humans and great apes) and Hylobatidae (gibbons), is well-known for sharing human-like characteristics, however, the genomic origins of these shared unique phenotypes have mainly remained elusive. To decipher the underlying genomic basis of Hominoidea-restricted phenotypes, we identified and characterized Hominoidea-restricted highly conserved noncoding sequences (HCNSs) that are a class of potential regulatory elements which may be involved in evolution of lineage-specific phenotypes. We discovered 679 such HCNSs from human, chimpanzee, gorilla, orangutan and gibbon genomes. These HCNSs were demonstrated to be under purifying selection but with lineage-restricted characteristics different from old CNSs. A significant proportion of their ancestral sequences had accelerated rates of nucleotide substitutions, insertions and deletions during the evolution of common ancestor of Hominoidea, suggesting the intervention of positive Darwinian selection for creating those HCNSs. In contrary to enhancer elements and similar to silencer sequences, these Hominoidea-restricted HCNSs are located in close proximity of transcription start sites. Their target genes are enriched in the nervous system, development and transcription, and they tend to be remotely located from the nearest coding gene. Chip-seq signals and gene expression patterns suggest that Hominoidea-restricted HCNSs are likely to be functional regulatory elements by imposing silencing effects on their target genes in a tissue-restricted manner during fetal brain development. These HCNSs, emerged through adaptive evolution and conserved through purifying selection, represent a set of promising targets for future functional studies of the evolution of Hominoidea-restricted phenotypes., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Is the origin of gibbon ape leukemia virus (GALV) human after all? When GALV was discovered and found to cause neoplastic disease in gibbons, it stimulated a great deal of research including investigations into the origins of this virus. A number of publications have suggested that the GALV progenitor was a retrovirus present in one of several species of South East Asian rodents that had close contact with captive gibbons. However, there are no published retroviral sequences from any South East Asian species to support this view. Here we present an alternative hypothesis that the origin of GALV is a virus closely related to Melomys burtoni retrovirus, and that this virus infected human patients in Papua New Guinea from whom biological material was obtained or in some way contaminated these samples. This material we propose contained infectious MbRV-related virus that was then unwittingly introduced into gibbons which subsequently developed GALV infections.

Urate oxidase (Uox) catalyses the first reaction of oxidative uricolysis, a three-step enzymatic pathway that allows some animals to eliminate purine nitrogen through a water-soluble compound. Inactivation of the pathway in hominoids leads to elevated levels of sparingly soluble urate and puts humans at risk of hyperuricemia and gout. The uricolytic activities lost during evolution can be replaced by enzyme therapy. Here we report on the functional and structural characterization of Uox from zebrafish and the effects on the enzyme of the missense mutation (F216S) that preceded Uox pseudogenization in hominoids. Using a kinetic assay based on the enzymatic suppression of the spectroscopic interference of the Uox reaction product, we found that the F216S mutant has the same turnover number of the wild-type enzyme but a much-reduced affinity for the urate substrate and xanthine inhibitor. Our results indicate that the last functioning Uox in hominoid evolution had an increased Michaelis constant, possibly near to upper end of the normal range of urate in the human serum (~300 μM). Changes in the renal handling of urate during primate evolution can explain the genetic modification of uricolytic activities in the hominoid lineage without the need of assuming fixation of deleterious mutations.

LINE-Alu-VNTR-Alu-like (LAVA) elements comprise a family of non-autonomous, composite, non-LTR retrotransposons specific to gibbons and may have played a role in the evolution of this lineage. A full-length LAVA element consists of portions of repeats found in most primate genomes: CT-rich, Alu-like, and VNTR regions from the SVA retrotransposon, and portions of the AluSz and L1ME5 elements. To evaluate whether the gibbon genome currently harbors functional LAVA elements capable of mobilization by the endogenous LINE-1 (L1) protein machinery and which LAVA components are important for retrotransposition, we established a trans-mobilization assay in HeLa cells. Specifically, we tested if a full-length member of the older LAVA subfamily C that was isolated from the gibbon genome and named LAVA C , or its components, can be mobilized in the presence of the human L1 protein machinery. We show that L1 proteins mobilize the LAVA C element at frequencies exceeding processed pseudogene formation and human SVA E retrotransposition by > 100-fold and ≥3-fold, respectively. We find that only the SVA-derived portions confer activity, and truncation of the 3' L1ME5 portion increases retrotransposition rates by at least 100%. Tagged de novo insertions integrated into intronic regions in cell culture, recapitulating findings in the gibbon genome. Finally, we present alternative models for the rise of the LAVA retrotransposon in the gibbon lineage., (© The Author 2016. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Evidence-based conservation planning is crucial for informing management decisions for species of extreme rarity, but collection of robust data on genetic status or other parameters can be extremely challenging for such species. The Hainan gibbon, possibly the world's rarest mammal, consists of a single population of ~25 individuals restricted to one protected area on Hainan Island, China, and has persisted for over 30 years at exceptionally low population size. Analysis of genotypes at 11 microsatellite loci from faecal samples for 36% of the current global population and tissue samples from 62% of existing historical museum specimens demonstrates limited current genetic diversity (Na = 2.27, Ar = 2.24, He  = 0.43); diversity has declined since the 19th century and even further within the last 30 years, representing declines of ~30% from historical levels (Na = 3.36, Ar = 3.29, He  = 0.63). Significant differentiation is seen between current and historical samples (FST  = 0.156, P = 0.0315), and the current population exhibits extremely small Ne (current Ne  = 2.16). There is evidence for both a recent population bottleneck and an earlier bottleneck, with population size already reasonably low by the late 19th century (historical Ne  = 1162.96). Individuals in the current population are related at the level of half- to full-siblings between social groups, and full-siblings or parent-offspring within a social group, suggesting that inbreeding is likely to increase in the future. The species' current reduced genetic diversity must be considered during conservation planning, particularly for expectations of likely population recovery, indicating that intensive, carefully planned management is essential., (© 2016 John Wiley & Sons Ltd.)

Background: For conservation of highly threatened species to be effective, it is crucial to differentiate natural population parameters from atypical behavioural, ecological and demographic characteristics associated with human disturbance and habitat degradation, which can constrain population growth and recovery. Unfortunately, these parameters can be very hard to determine for species of extreme rarity. The Hainan gibbon (Nomascus hainanus), the world's rarest ape, consists of a single population of c.25 individuals, but intensive management is constrained by a limited understanding of the species' expected population characteristics and environmental requirements. In order to generate a more robust evidence-base for Hainan gibbon conservation, we employed a comparative approach to identify intrinsic and extrinsic drivers of variation in key ecological and behavioural traits (home range size, social group size, mating system) across the Hylobatidae while controlling for phylogenetic non-independence., Results: All three studied traits show strong phylogenetic signals across the Hylobatidae. Although the Hainan gibbon and some closely related species have large reported group sizes, no observed gibbon group size is significantly different from the values expected on the basis of phylogenetic relationship alone. However, the Hainan gibbon and two other Nomascus species (N. concolor, N. nasutus) show home range values that are higher than expected relative to all other gibbon species. Predictive models incorporating intraspecific trait variation but controlling for covariance between population samples due to phylogenetic relatedness reveal additional environmental and biological determinants of variation in gibbon ranging requirements and social structure, but not those immediately associated with recent habitat degradation., Conclusions: Our study represents the first systematic assessment of behavioural and ecological trait patterns across the Hylobatidae using recent approaches in comparative analysis. By formally contextualising the Hainan gibbon's observed behavioural and ecological characteristics within family-wide variation in gibbons, we are able to determine natural population parameters expected for this Critically Endangered species, as well as wider correlates of variation for key population characteristics across the Hylobatidae. This approach reveals key insights with a direct impact on future Hainan gibbon conservation planning, and demonstrates the usefulness of the comparative approach for informing management of species of conservation concern.

Gibbons are believed to have diverged from the larger great apes ∼16.8 MYA and today reside in the rainforests of Southeast Asia. Based on their diploid chromosome number, the family Hylobatidae is divided into four genera, Nomascus, Symphalangus, Hoolock, and Hylobates. Genetic studies attempting to elucidate the phylogenetic relationships among gibbons using karyotypes, mitochondrial DNA (mtDNA), the Y chromosome, and short autosomal sequences have been inconclusive . To examine the relationships among gibbon genera in more depth, we performed second-generation whole genome sequencing (WGS) to a mean of ∼15× coverage in two individuals from each genus. We developed a coalescent-based approximate Bayesian computation (ABC) method incorporating a model of sequencing error generated by high coverage exome validation to infer the branching order, divergence times, and effective population sizes of gibbon taxa. Although Hoolock and Symphalangus are likely sister taxa, we could not confidently resolve a single bifurcating tree despite the large amount of data analyzed. Instead, our results support the hypothesis that all four gibbon genera diverged at approximately the same time. Assuming an autosomal mutation rate of 1 × 10(-9)/site/year this speciation process occurred ∼5 MYA during a period in the Early Pliocene characterized by climatic shifts and fragmentation of the Sunda shelf forests. Whole genome sequencing of additional individuals will be vital for inferring the extent of gene flow among species after the separation of the gibbon genera., (Copyright © 2015 by the Genetics Society of America.)

Amino acid usage varies from species to species. A previous study has found a universal trend in amino acid gain and loss in many taxa and a one-way model of amino acid evolution in which the number of new amino acids increases as the number of old amino acids decreases was proposed. Later studies showed that this pattern of amino acid gain and loss is likely to be compatible with the neutral theory. The present work aimed to further study this problem by investigating the evolutionary patterns of amino acids in 8 primates (the nucleotide and protein alignments are available online http://gattaca.nju.edu.cn/pub_data.html). First, the number of amino acids gained and lost was calculated and the evolution trend of each amino acid was inferred. These values were found to be closely related to the usage of each amino acid. Then we analyzed the mutational trend of amino acid substitution in human using SNPs, this trend is highly correlated with fixation trend only with greater variance. Finally, the trends in the evolution of 20 amino acids were evaluated in human on different time scales, and the increasing rate of 5 significantly increasing amino acids was found to decrease as a function of time elapsed since divergence, and the dS/dN ratio also found to increase as a function of time elapsed since divergence. These results suggested that the observed amino acid substitution pattern is influenced by mutation and purifying selection. In conclusion, the present study shows that usage of amino acids is an important factor capable of influencing the observed pattern of amino acid evolution, and also presented evidences suggesting that the observed universal trend of amino acid gain and loss is compatible with neutral evolution., (Copyright © 2014 Elsevier B.V. All rights reserved.)

SVA (SINE-R-VNTR-Alu) elements constitute the youngest family of composite non-LTR retrotransposons in hominoid primates. The sequence of their assembly, however, remains unclear. Recently, a second family of VNTR-containing composites, LAVA (L1-Alu-VNTR-Alu), has been identified in gibbons. We now report the existence of two additional VNTR composite families, PVA (PTGR2-VNTR-Alu) and FVA (FRAM-VNTR-Alu), in the genome of Nomascus leucogenys. Like LAVA, they share the 5'-Alu-like region and VNTR with SVA, but differ at their 3'-ends. The 3'-end of PVA comprises part of the PTGR2 gene, whereas FVA is characterized by the presence of a partial FRAM element in its 3'-domain. Splicing could be identified as the mechanism of acquisition of the variant 3'-ends in all four families of VNTR composites. SVAs have been shown to be mobilized by the L1 protein machinery in trans. A critical role in this process has been ascribed to their 5'-hexameric repeat/ Alu-like region. The Alu-like region displays specific features in each of the VNTR composite families/subfamilies with characteristic deletions found in the evolutionary younger subfamilies. Using reciprocal exchanges between SVA_E and PVA/FVA elements, we demonstrate that the structure, not the presence of the (CCCTCT)n/ Alu-like region determines mobilization capacity. Combination of LAVA and SVA_E domains does not yield any active elements-suggesting the use of different combinations of host factors for the two major groups of VNTR composites. Finally, we demonstrate that the LAVA 3'-L1ME5 fragment attenuates mobilization capacity., (© The Author 2014. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

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.

Gibbons are small arboreal apes that display an accelerated rate of evolutionary chromosomal rearrangement and occupy a key node in the primate phylogeny between Old World monkeys and great apes. Here we present the assembly and analysis of a northern white-cheeked gibbon (Nomascus leucogenys) genome. We describe the propensity for a gibbon-specific retrotransposon (LAVA) to insert into chromosome segregation genes and alter transcription by providing a premature termination site, suggesting a possible molecular mechanism for the genome plasticity of the gibbon lineage. We further show that the gibbon genera (Nomascus, Hylobates, Hoolock and Symphalangus) experienced a near-instantaneous radiation ∼5 million years ago, coincident with major geographical changes in southeast Asia that caused cycles of habitat compression and expansion. Finally, we identify signatures of positive selection in genes important for forelimb development (TBX5) and connective tissues (COL1A1) that may have been involved in the adaptation of gibbons to their arboreal habitat.

Variable number of tandem repeats (VNTRs) are scattered throughout the primate genome, and genetic variation of these VNTRs have been accumulated during primate radiation. Here, we analyzed VNTRs upstream of the monoamine oxidase A (MAOA) gene in 11 different gibbon species. An abundance of truncated VNTR sequences and copy number differences were observed compared to those of human VNTR sequences. To better understand the biological role of these VNTRs, a luciferase activity assay was conducted and results indicated that selected VNTR sequences of the MAOA gene from human and three different gibbon species (Hylobates klossii, Hylobates lar, and Nomascus concolor) showed silencing ability. Together, these data could be useful for understanding the evolutionary history and functional significance of MAOA VNTR sequences in gibbon species.

Tandem repeats are common in eukaryotic genomes, but due to difficulties in assaying them remain poorly studied. Here, we demonstrate the utility of Nanostring technology as a targeted approach to perform accurate measurement of tandem repeats even at extremely high copy number, and apply this technology to genotype 165 HapMap samples from three different populations and five species of non-human primates. We observed extreme variability in copy number of tandemly repeated genes, with many loci showing 5-10 fold variation in copy number among humans. Many of these loci show hallmarks of genome assembly errors, and the true copy number of many large tandem repeats is significantly under-represented even in the high quality 'finished' human reference assembly. Importantly, we demonstrate that most large tandem repeat variations are not tagged by nearby SNPs, and are therefore essentially invisible to SNP-based GWAS approaches. Using association analysis we identify many cis correlations of large tandem repeat variants with nearby gene expression and DNA methylation levels, indicating that variations of tandem repeat length are associated with functional effects on the local genomic environment. This includes an example where expansion of a macrosatellite repeat is associated with increased DNA methylation and suppression of nearby gene expression, suggesting a mechanism termed "repeat induced gene silencing", which has previously been observed only in transgenic organisms. We also observed multiple signatures consistent with altered selective pressures at tandemly repeated loci, suggesting important biological functions. Our studies show that tandemly repeated loci represent a highly variable fraction of the genome that have been systematically ignored by most previous studies, copy number variation of which can exert functionally significant effects. We suggest that future studies of tandem repeat loci will lead to many novel insights into their role in modulating both genomic and phenotypic diversity.

Knowledge of the genetic mating system of animal species is essential for our understanding of the evolution of social systems and individual reproductive strategies. In recent years, genetic methods have uncovered an unexpected diversity of paternal genetic contributions across diverse animal social mating systems, but particularly in pair-living species. In most pair-living birds, for example, genetic and behavioral observations have confirmed a previously unknown significance of extra-pair copulations (EPCs) and extra-pair paternity. Among mammals, white-handed gibbons (Hylobates lar) are also known to live in pairs and are traditionally believed to be single-male single-female breeders. However, at Khao Yai National Park, Thailand, behavioral observations have confirmed the occurrence of both EPCs and functional multi-male grouping, but knowledge about the genetic mating system is still unavailable. In this study, we genotyped 89 white-handed gibbons of the Khao Yai population based on fecal samplings and were able to determine paternity for 41 offspring through short tandem repeat analysis. We found that females' stable social partners sired the majority (90.5%) of offspring (N = 38), while only a few (7.1%) offspring (n = 2 confirmed cases; n = 1 inferred case) were conceived with extra-pair partners. The paternity of one offspring remained inconclusive (2.4%), because the offspring's genotype did not mismatch with the genotypes of two potential sires. Like other predominantly pair-living species, gibbons appear to follow a mixed-reproductive strategy. The genetic mating system of wild white-handed gibbons is best described as flexible, primarily monogamous and opportunistically promiscuous. Inc., (© 2013 Wiley Periodicals, Inc.)

Alpha satellite DNA (AS), a major DNA component of primate centromeres, is composed of a tandem array of repeat units of approximately 170 bp. The AS of hominids (family Hominidae; humans and great apes) includes sequences organized into higher-order repeat (HOR) structures, with a periodic appearance of multiple copies of the basic repeat units. Here, we identified an HOR in AS of the siamang, a small ape phylogenetically distinct from hominids but included in hominoids (superfamily Hominoidea). We sequenced long stretches of genomic DNA, and found a repetition of blocks consisting of six and four basic repeat units. Thus, AS organization into HOR is an attribute of hominoids, rather than, as currently postulated, hominids. In addition to centromeres, siamangs carry AS in terminal heterochromatin blocks, and it cannot be determined at present whether these HOR-containing AS sequences originate from the centromere or from the terminal heterochromatin. Even if the latter is the case, these sequences might affect the composition of centromeric AS by being transferred to the centromere.

Hoolock gibbons (genus Hoolock) are a group of very endangered primate species that belong to the small ape family (family Hylobatidae). The entire population that is distributed in the northeast and southeast of Bangladesh is estimated to include only around 350 individuals. A conservation program is thus necessary as soon as possible. Genetic markers are significant tools for planning such programs. In this study, we examined chromosomal characteristics of two western hoolock gibbons that were captured in a Bangladesh forest. During chromosome analysis, we encountered two chromosome variations that were observed for the first time in the wild-born western hoolock gibbons (Hoolock hoolock). The first one was a nonhomologous centromere position in chromosome 8 that was observed in the two examined individuals. The alteration was identical in the two individuals, which were examined by G-band and DAPI-band analyses. Chromosome paint analyses revealed that the difference in the centromere position was due to a single small pericentric inversion. The second variation was a heterozygous elongation in chromosome 9. Analysis by sequential techniques of fluorescence in situ hybridization with 18S rDNA and silver nitrate staining revealed a single and an inverted tandem duplication, respectively, of the nucleolus organizer region in two individuals. These chromosome variations provide useful information for the next steps to consider the evolution and conservation of the hoolock gibbon.

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.

Recombinations are known to disrupt bifurcating tree structure of gene genealogies. Although recently occurred recombinations are easily detectable by using conventional methods, recombinations may have occurred at any time. We devised a new method for detecting ancient recombinations through phylogenetic network analysis, and detected five ancient recombinations in gibbon ABO blood group genes [Kitano et al., 2009. Mol. Phylogenet. Evol., 51, 465-471]. We present applications of this method, now named as "PNarec", to various virus sequences as well as HLA genes., (Copyright © 2012 Elsevier Inc. All rights reserved.)

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.

Inter-specific hybrid zones for Hylobates gibbons are known in Southeast Asia. Among these, one hybrid zone between Hylobates lar and H. pileatus is located in Khao Yai National Park, Thailand. To find molecular evidence for the natural hybridization of the gibbons in this region, we studied mitochondrial DNA (mtDNA) of 68 gibbons of the H. lar phenotype living adjacent to the hybrid zone. Nucleotide sequencing of a fragment of mtDNA spanning hyper variable segment I showed that nine gibbons had an mtDNA haplotype of H. pileatus, and that seven of these nine gibbons belonged to a single maternal lineage over three generations. It is thus confirmed that introgression between H. lar and H. pileatus exists and the initial hybridization took place ages ago.

Each genus of small apes has a highly distinctive karyotype (karyomorph) at every level of cytogenetic analysis. Early workers using classical staining and banding had problems integrating the karyolocial data with that of other primates. Chromosome painting allowed syntenic homology maps to be constructed for each of the four karyomorphs (2n = 38, 44, 50 and 52). They revealed that the great apes and Old World monkeys had strongly conserved karyotypes while those of small apes were highly rearranged. However, they provided contradictory phylogenetic results to other bio-molecular tree of small ape evolution. More recently BAC-FISH investigations using a panel of about 900 BACs defined each breakpoint by spanning or flanking BAC clones The syntenic map was refined and now includes small segments of homology which had previously gone undected, marker order (synteny block orientation) and the location of ancestral and Evolutionarily New Centromeres. However, the BAC-FISH data similar to other biomolecular methods used up to now could not resolve the phylogenetic tree of hylobatids. These difficulties may be explained by the rapid divergence of crown hylobatids, reticulate evolution and incomplete lineage sorting. The lack of significant cytogenetic landmarks at the nodes of the gibbon tree could indicate that chromosomal rearrangements did not play a primary role in hylobatid speciation.

Chromosome rearrangements in small apes are up to 20 times more frequent than in most mammals. Because of their complexity, the full extent of chromosome evolution in these hominoids is not yet fully documented. However, previous work with array painting, BAC-FISH, and selective sequencing in two of the four karyomorphs has shown that high-resolution methods can precisely define chromosome breakpoints and map the complex flow of evolutionary chromosome rearrangements. Here we use these tools to precisely define the rearrangements that have occurred in the remaining two karyomorphs, genera Symphalangus (2n = 50) and Hoolock (2n = 38). This research provides the most comprehensive insight into the evolutionary origins of chromosome rearrangements involved in transforming small apes genome. Bioinformatics analyses of the human-gibbon synteny breakpoints revealed association with transposable elements and segmental duplications, providing some insight into the mechanisms that might have promoted rearrangements in small apes. In the near future, the comparison of gibbon genome sequences will provide novel insights to test hypotheses concerning the mechanisms of chromosome evolution. The precise definition of synteny block boundaries and orientation, chromosomal fusions, and centromere repositioning events presented here will facilitate genome sequence assembly for these close relatives of humans.

The siamang (Symphalangus syndactylus), a species of the family Hylobatidae (gibbons), carries large blocks of constitutive heterochromatin in the telomere region of chromosomes. We recently found that alpha satellite DNA constitutes these heterochromatin blocks as a main component. Alpha satellite DNA, tandem repeat sequences of 171-bp repeat units, is a major component of centromeres in primates. In addition to the siamang, the white-cheeked gibbon (Nomascus leucogenys) was previously found to carry the alpha satellite DNA in the telomere region, although not as large a scale as the siamang. Gibbons comprise four genera: Hoolock, Hylobates, Nomascus, and Symphalangus. Here, we report that the amplification of alpha satellite DNA in the telomere region is probably confined to two genera: Nomascus and Symphalangus. We examined one species of Hoolock and four species of Hylobates and obtained evidence against such an amplification event in these species. The phylogenetic relationship of the four gibbon genera remains unclear. One simple explanation for the current distribution of the telomere region alpha satellite DNA would be that Nomascus and Symphalangus are relatively closely related and the amplification occurred in their common ancestor.

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.

Chromosomes of the siamang Symphalangus syndactylus (a small ape) carry large-scale heterochromatic structures at their ends. These structures look similar, by chromosome C-banding, to chromosome-end heterochromatin found in chimpanzee, bonobo and gorilla (African great apes), of which a major component is tandem repeats of 32-bp-long, AT-rich units. In the present study, we identified repetitive sequences that are a major component of the siamang heterochromatin. Their repeat units are 171 bp in length, and exhibit sequence similarity to alpha satellite DNA, a major component of the centromeres in primates. Thus, the large-scale heterochromatic structures have different origins between the great apes and the small ape. The presence of alpha satellite DNA in the telomere region has previously been reported in the white-cheeked gibbon Nomascus leucogenys, another small ape species. There is, however, a difference in the size of the telomere-region alpha satellite DNA, which is far larger in the siamang. It is not known whether the sequences of these two species (of different genera) have a common origin because the phylogenetic relationship of genera within the small ape family is still not clear. Possible evolutionary scenarios are discussed.

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.

Gibbons (Hylobatidae) shared a common ancestor with the other hominoids only 15-18 million years ago. Nevertheless, gibbons show very distinctive features that include heavily rearranged chromosomes. Previous observations indicate that this phenomenon may be linked to the attenuated epigenetic repression of transposable elements (TEs) in gibbon species. Here we describe the massive expansion of a repeat in almost all the centromeres of the eastern hoolock gibbon (Hoolock leuconedys). We discovered that this repeat is a new composite TE originating from the combination of portions of three other elements (L1ME5, AluSz6, and SVA_A) and thus named it LAVA. We determined that this repeat is found in all the gibbons but does not occur in other hominoids. Detailed investigation of 46 different LAVA elements revealed that the majority of them have target site duplications (TSDs) and a poly-A tail, suggesting that they have been retrotransposing in the gibbon genome. Although we did not find a direct correlation between the emergence of LAVA elements and human-gibbon synteny breakpoints, this new composite transposable element is another mark of the great plasticity of the gibbon genome. Moreover, the centromeric expansion of LAVA insertions in the hoolock closely resembles the massive centromeric expansion of the KERV-1 retroelement reported for wallaby (marsupial) interspecific hybrids. The similarity between the two phenomena is consistent with the hypothesis that evolution of the gibbons is characterized by defects in epigenetic repression of TEs, perhaps triggered by interspecific hybridization.

We investigated the ecology and interspecific interactions of the two gibbon species (Hylobates lar and H. pileatus) that overlap in distribution within a narrow zone of contact in the headwaters of the Takhong River at Khao Yai National Park in central Thailand. The zone is about 10-km wide, with phenotypic hybrids comprising 6.5% of the adult population. We compared species with respect to diet, territory size, intra- and interspecific group encounters, and territory quality. The two gibbon species exploited the same types of resources within their territories despite variation in the relative abundance of food-plant species between territories. The gibbons were interspecifically territorial, and males of both species displayed aggressive behaviors at both intraspecific and interspecific territorial boundaries. There was no marked difference in the amount of overlap between territories of conspecific and heterospecific pairs of groups. Although the habitat was not homogeneous, territory quality did not vary significantly between species. The species have not diverged in habitat preference or in diet. Neither species dominated in interspecific encounters, and both were reproducing well in the contact zone. We analyzed the potential consequences of several types of interspecific interactions on individual dispersal options and on the structure of the contact zone. Interference competition through interspecific territoriality affects the dispersal of individuals into the range of the other species. In general, territorial competition coupled with limited hybridization leads to predictions of a narrow contact zone or parapatry between species; thus, behavioral and ecological interactions between species need to be considered as potential factors in explaining range borders of primate species.

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.

Gibbons are small, arboreal, highly endangered apes that are understudied compared with other hominoids. At present, there are four recognized genera and approximately 17 species, all likely to have diverged from each other within the last 5-6 My. Although the gibbon phylogeny has been investigated using various approaches (i.e., vocalization, morphology, mitochondrial DNA, karyotype, etc.), the precise taxonomic relationships are still highly debated. Here, we present the first survey of nuclear sequence variation within and between gibbon species with the goal of estimating basic population genetic parameters. We gathered ~60 kb of sequence data from a panel of 19 gibbons representing nine species and all four genera. We observe high levels of nucleotide diversity within species, indicative of large historical population sizes. In addition, we find low levels of genetic differentiation between species within a genus comparable to what has been estimated for human populations. This is likely due to ongoing or episodic gene flow between species, and we estimate a migration rate between Nomascus leucogenys and N. gabriellae of roughly one migrant every two generations. Together, our findings suggest that gibbons have had a complex demographic history involving hybridization or mixing between diverged populations.