Molecular phylogeny and evolution of the human endogenous retrovirus HERV-W LTR family in hominoid primates

Long terminal repeats (LTRs) of human endogenous retrovirus (HERV) have contributed to the structural change or genetic variation of primate genome that are connected to speciation and evolution. Using genomic DNAs that were derived from hominoid primates (chimpanzee, gorilla, orangutan, and gibbon), we performed PCR amplification and identified thirty HERV-W LTR elements. These LTR elements showed a 82-98% sequence similarity with HERV-W LTR (AF072500). Specifically, additional sequences (GCCACCACCACTGTTT in the gorilla and TGCTGCTGACTCCCATCC in the gibbon) were noticed. Clone OR3 from the orangutan and clone GI2 from the gibbon showed a 100% sequence similarity, although they are different species. This indicates that both LTR elements were proliferated during the last 2 to 5 million years from the integration of the original LTR element. A phylogenetic tree that was obtained by the neighbor-joining method revealed a wide overlap of the LTR elements across species, suggesting that the HERV-W LTR family evolved independently during the hominoid evolution.

Comparative mapping of CDY and DAZ in higher primates

The human male specific expressed gene families CDY and DAZ are known to be repetitively clustered in the Y-specific region of the human Y chromosome. Comparative FISH-mapping of DNA clones specific for CDY and DAZ resulted in a Y-specific but diverse signal pattern within the non-recombining region of the Y-chromosomes of human and great apes. It can be concluded that the non-recombining part of the Y-chromosomes including CDY and DAZ, was exposed to species-specific amplifications, diversifications and rearrangements. Evolutionary fast fixation of any of these variations was possible as long as they did not interfere with male fertility.

Population subdivision and gene flow among wild orangutans

Genetic variability among populations of orangutans from Borneo and Sumatra was assessed using seven SSR loci. Most SSR loci were highly polymorphic and their allele frequencies exhibited substantial variation across subpopulations. While significant genetic subdivision was observed among the island populations, genetic distance did not increase with geographic distance and sufficient gene flow persists to prevent marked genetic subdivision. Since it is unlikely that the Bornean Orangutans dispersed naturally among locations separated by such formidable geographic barriers, human assistance might already have altered their genetic structure. Our data suggests that there may be at least two subspecific clades of orangutans within Borneo while Central Kalimantan animals may have become more genetically related to animals in Sumatra due to human intervention.

NK cell receptors of the orangutan (Pongo pygmaeus): a pivotal species for tracking the coevolution of killer cell Ig-like receptors with MHC-C

CD94, NKG2, Ly49, and killer cell Ig-like receptor (KIR) expressed by orangutan peripheral blood cells were examined by cloning and sequencing cDNA from a panel of individuals. Orthologs of human CD94, NKG2A, D, and F were defined. NKG2C and E are represented by one gene, Popy-NKG2CE, that is equidistant from the two human genes. Several Popy-CD94, NKG2A, and NKG2CE alleles were defined. Popy-Ly49L is expressed in cultured NK cells and has a sequence consistent with it encoding a functional receptor. Orangutan KIR corresponding to the three KIR lineages expressed in humans and chimpanzees were defined. Popy-KIR2DL4 of lineage I is the only ortholog of a human or chimpanzee KIR, but in all individuals examined, the transcripts of this gene produced premature termination, either in the D2 domain or at the beginning of the cytoplasmic domain. Ten Popy-KIR3DL and one Popy-KIR3DS of lineage II are all closely related, but represent the products of at least two genes. The two Popy-KIR2DL and four Popy-KIR2DS of lineage III also represent two genes, both being more related to KIR2DS4 than to other human and chimpanzee KIR of lineage III. The Popy-KIR2D include ones predicted to be specific for the C1 epitope of MHC-C, but none specific for C2. This correlates with the observation that all orangutan MHC-C allotypes examined have the C1 motif.

Distinctive distribution of AIM1 polymorphism among major human populations with different skin color

The genetic background for human skin color has been a major topic in human genetics; however, its molecular basis is still unclear. The gene for the AIM-1 protein (AIM1) was recently found to be responsible for the body color of medaka fish. In the search for the genes controlling human skin color variations, we have investigated genetic polymorphisms of this gene, and we have found a single-nucleotide polymorphism that has clear association with major human populations in terms of skin color.

Intra- and interspecific variation in primate gene expression patterns

Although humans and their closest evolutionary relatives, the chimpanzees, are 98.7% identical in their genomic DNA sequences, they differ in many morphological, behavioral, and cognitive aspects. The underlying genetic basis of many of these differences may be altered gene expression. We have compared the transcriptome in blood leukocytes, liver, and brain of humans, chimpanzees, orangutans, and macaques using microarrays, as well as protein expression patterns of humans and chimpanzees using two-dimensional gel electrophoresis. We also studied three mouse species that are approximately as related to each other as are humans, chimpanzees, and orangutans. We identified species-specific gene expression patterns indicating that changes in protein and gene expression have been particularly pronounced in the human brain.

Species-specific evolution of repeated DNA sequences in great apes

DNA sequencing reveals that the genomes of the human, gorilla and chimpanzee share more than 98% homology. Comparative chromosome painting and gene mapping have demonstrated that only a few rearrangements of a putative ancestral mammalian genome occurred during great ape and human evolution. However, interspecies representational difference analysis (RDA) of the gorilla between human and gorilla revealed gorilla-specific DNA sequences. Cloning and sequencing of gorilla-specific DNA sequences indicate that there are repetitive elements. Gorilla-specific DNA sequences were mapped by fluorescence in-situ hybridization (FISH) to the subcentromeric/centromeric regions of three pairs of gorilla submetacentric chromosomes. These sequences could represent either ancient sequences that got lost in other species, such as human and orang-utan, or, more likely, recent sequences which evolved or originated specifically in the gorilla genome.

Telomere biology and cellular aging in nonhuman primate cells

To determine how cellular aging is conserved among primates, we analyzed the replicative potential and telomere shortening in skin fibroblasts of anthropoids and prosimians. The average telomere length of the New World primates Ateles geoffroyi (spider monkey) and Saimiri sciureus (squirrel monkey) and the Old World primates Macaca mulatta (rhesus monkey), Pongo pygmaeus (orangutan), and Pan paniscus (pigmy chimpanzee) ranged from 4 to 16 kb. We found that telomere shortening limits the replicative capacity of anthropoid fibroblasts and that the expression of human telomerase produced telomere elongation and the extension of their in vitro life span. In contrast the prosimian Lemur catta (ring-tailed lemur) had both long and short telomeres and telomere shortening did not provide an absolute barrier to immortalization. Following a transient growth arrest a subset of cells showing a reduced number of chromosomes overgrew the cultures without activation of telomerase. Here we show that the presence of continuous TTAGGG repeats at telomeres and rigorous control of replicative aging by telomere shortening appear to be conserved among anthropoid primates but is less effective in prosimian lemurs.

Characterization of the MHC class I-related MR1 locus in nonhuman primates

We characterized the MHC-related 1 ( MR1) locus in two nonhuman primates species, Pongo pygmaeus and Pan troglodytes. MR1 cDNA sequences encoding several isoforms generated through alternative splicing were observed in both species. Amino acid alignment between the five species in which MR1 has been characterized to date - human, chimpanzee, orangutan, mouse, and rat - reveals a very high degree of conservation specially in the alpha1 and alpha2 domains of the molecule. The main differences concentrate in the transmembrane and cytoplasmic domains. In the three primates species there is a lysine residue inside the putative transmembrane domain which is not present in rodents. Furthermore, the MR1 cytoplasmic region is longer in rodents, with a conserved serine-containing motif that could be involved in endocytosis; remarkably, this motif is absent in the three primate species. We also describe the presence in the chimpanzee of a sequence homologous to the MR1P1 pseudogene previously found in humans.

Alu-mediated inactivation of the human CMP- N-acetylneuraminic acid hydroxylase gene

Inactivation of the CMP-N-acetylneuraminic acid hydroxylase gene has provided an example of human-specific genomic mutation that results in a widespread biochemical difference between human and nonhuman primates. We have found that, although a region containing a 92-bp exon and an AluSq element in the hydroxylase gene is intact in all nonhuman primates examined, the same region in the human genome is replaced by an AluY element that was disseminated at least one million years ago. We propose a mechanistic model for this Alu-mediated replacement event, which deleted the 92-bp exon and thus inactivated the human hydroxylase gene. It is suggested that Alu elements have played potentially important roles in genotypic and phenotypic evolution in the hominid lineage.

Sequence variation in humans and other primates at six short tandem repeat loci used in forensic identity testing

A large number of alleles from the six different short tandem repeat (STR) loci FGA, D3S1358, vWA, CSF1PO, TPOX and TH01, used in human identity testing were sequenced to provide support for the robustness of fluorescent STR DNA typing by allele size. Sequence information for some of these loci (FGA, vWA, TH01) is an extension of published work, whereas no extensive sequence information is available with respect to the D3S1358, CSF1PO, and TPOX loci. Sequencing of alleles at each locus has provided quantitative data with respect to the true nucleotide length of common alleles, and of alleles that vary in length from the common alleles. All alleles that were identified as "off-ladder" alleles through fluorescent typing at these STR loci have proven to be true length variant alleles. Sequencing at the D3S1358 and CSF1PO loci allowed for the establishment of a common nomenclature for these loci. A correlation between percent stutter and the length of the core tandem repeat is demonstrated at the FGA locus. Alleles in which the core tandem repeat is interrupted by a repeat unit of different sequence have a reduced percent stutter. DNA samples from three non-human primates (chimpanzee, orangutan, and gorilla) were compared to the human sequences, and shown to differ markedly across loci with respect to their homology. The effects of primer binding site mutations on the amplification efficiency at a particular locus, and methods used to interpret amplification imbalance of heterozygous alleles at a locus is also addressed.

"Bar-coding" primate chromosomes: molecular cytogenetic screening for the ancestral hominoid karyotype

Two recently introduced multicolor FISH approaches, cross-species color banding (also termed Rx-FISH) and multiplex FISH using painting probes derived from somatic cell hybrids retaining fragments of human chromosomes, were applied in a comparative molecular cytogenetic study of higher primates. We analyzed these "chromosome bar code" patterns to obtain an overview of chromosomal rearrangements that occurred during higher primate evolution. The objective was to reconstruct the ancestral genome organization of hominoids using the macaque as outgroup species. Approximately 160 individual and discernible molecular cytogenetic markers were assigned in these species. Resulting comparative maps allowed us to identify numerous intra-chromosomal rearrangements, to discriminate them from previous contradicting chromosome banding interpretations and to propose an ancestral karyotype for hominoids. From 25 different chromosome forms in an ancestral karyotype for all hominoids of 2N=48 we propose 21. Probes for chromosomes 2p, 4, 9 and Y were not informative in the present experiments. The orangutan karyotype was very similar to the proposed ancestral organization and conserved 19 of the 21 ancestral forms; thus most chromosomes were already present in early hominoid evolution, while African apes and human show various derived changes.

Genetic divergence of orangutan subspecies (Pongo pygmaeus)

Microsatellites and mitochondrial DNA sequences were studied for the two subspecies of orangutans (Pongo pygmaeus), which are located in Borneo (P. p. pygmaeus) and Sumatra (P. p. abelii), respectively. Both subspecies possess marked genetic diversity. Genetic subdivision was identified within the Sumatran orangutans. The genetic differentiation between the two subspecies is highly significant for ND5 region but not significant for 16s rRNA or microsatellite data by exact tests, although FST estimates are highly significant for these markers. Divergence time between the two subspecies is approximately 2.3 +/- 0.5 million years ago (MYA) estimated from our data, much earlier than the isolation of their geological distribution. Neither subspecies underwent a recent bottleneck, though the Sumatran subspecies might have experienced expansion approximately 82,000 years ago. The estimated effective population sizes for both subspecies are on the order of 104. Our results contribute additional information that may be interpreted in the context of orangutan conservation efforts.

A population genetics study of single nucleotide polymorphisms in the interleukin 4 receptor alpha (IL4RA) gene

Interleukin 4 (IL4) plays a critical role in T helper 2 (Th2) immune responses. Here we report a population genetics study of variation in the gene encoding the alpha-chain of the IL4 receptor (IL4RA) in three ethnic groups: African Americans, European Americans and East Asians. A 2941-bp region spanning exon 12 of IL4RA gene was sequenced in 12 individuals from each group. A total of 24 single nucleotide polymorphisms (SNPs) were identified in the combined sample. The genetic variation of the coding region of exon 12 is two to three times higher than in other reported genes. A significant departure from the expectation of evolutionary neutrality was observed, suggesting that natural selection may have influenced the evolution of this gene. We propose a model in which past selection by pathogens contributed to the increasing prevalence of atopic disorders in Western societies.

Speciation and intrasubspecific variation of Bornean orangutans, Pongo pygmaeus pygmaeus

Mitochondrial DNA control region sequences of orangutans (Pongo pygmaeus) from six different populations on the island of Borneo were determined and analyzed for evidence of regional diversity and were compared separately with orangutans from the island of Sumatra. Within the Bornean population, four distinct subpopulations were identified. Furthermore, the results of this study revealed marked divergence, supportive evidence of speciation between Sumatran and Bornean orangutans. This study demonstrates that, as an entire population, Bornean orangutans have not experienced a serious genetic bottleneck, which has been suggested as the cause of low diversity in humans and east African chimpanzees. Based on these new data, it is estimated that Bornean and Sumatran orangutans diverged approximately 1.1 MYA and that the four distinct Bornean populations diverged 860,000 years ago. These findings have important implications for management, breeding, and reintroduction practices in orangutan conservation efforts.

Analyses of cross species polymerase chain reaction products to infer the ancestral state of human polymorphisms

In numerous population genetic and disease association studies decisions about the ancestry of polymorphic alleles are often made based on the relative frequency of the alleles in the extant populations with the most frequent allele being deemed as ancestral. However, the frequency of an allele in a population is generally not a perfect indicator of its ancestral status. A more accurate method to assess ancestral/derived status of polymorphic alleles involves identification of shared alleles between species. We used this strategy to examine genomic regions homologous to several human polymorphisms in four species of non-human primates. Cross species polymerase chain reaction (CS-PCR), with primers designed from human sequence, was used to investigate regions of interest. Nineteen polymorphisms at six loci (DRD2, HOXB@, PAH, D4S10, RBP3, and RET) were examined either by restriction fragment length analysis of PCR products (PCR-RFLP) or by direct sequencing. At seventeen of the eighteen PCR-RFLPs, non-human primates were monomorphic and identical to each other for either lack of restriction enzyme site or presence of the site. Thus, at these seventeen polymorphic sites the shared alleles are most likely to be the ancestral ones in humans. In several cases we have used sequence data to further demonstrate that the nucleotide at the site of the polymorphism is conserved between species confirming the hypothesis of a single ancestral allele. However, not all human alleles can be simply resolved into ancestral and derived; sequence data from one PCR-RFLP (in an intron of the PAH locus) and a single strand conformational polymorphism (SSCP) in the 3' untranslated region (UTR) of the DRD2 gene illustrate this point.

Evolution of a D2 dopamine receptor intron within the great apes and humans

Although direct DNA sequencing may allow rapid and high quality comparative phylogenetic analyses among species, such an approach may not be the most efficient method by which to make a large number of cross-species comparisons. We illustrate the use of Denaturing Gradient Gel Electrophoresis (DGGE) to screen a D2 Dopamine Receptor intron for DNA sequence variation, both within and between closely related species, in order to infer their evolutionary relationships. Our results suggest that: a) humans have less genetic variation than the great apes; b) pygmy chimpanzees have less genetic variation than common chimpanzees; and c) DNA sequence comparative analyses of primates require adequate sampling, both in number and in geographical range.

Attempts to detect retrotransposition and de novo deletion of Alus and other dispersed repeats at specific loci in the human genome

Dispersed repeat elements contribute to genome instability by de novo insertion and unequal recombination between repeats. To study the dynamics of these processes, we have developed single DNA molecule approaches to detect de novo insertions at a single locus and Alu-mediated deletions at two different loci in human genomic DNA. Validation experiments showed these approaches could detect insertions and deletions at frequencies below 10(-6) per cell. However, bulk analysis of germline (sperm) and somatic DNA showed no evidence for genuine mutant molecules, placing an upper limit of insertion and deletion rates of 2 x 10(-7) and 3 x 10(-7), respectively, in the individuals tested. Such re-arrangements at these loci therefore occur at a rate lower than that detectable by the most sensitive methods currently available.

Genomic divergences between humans and other hominoids and the effective population size of the common ancestor of humans and chimpanzees

To study the genomic divergences among hominoids and to estimate the effective population size of the common ancestor of humans and chimpanzees, we selected 53 autosomal intergenic nonrepetitive DNA segments from the human genome and sequenced them in a human, a chimpanzee, a gorilla, and an orangutan. The average sequence divergence was only 1.24% +/- 0.07% for the human-chimpanzee pair, 1.62% +/- 0.08% for the human-gorilla pair, and 1.63% +/- 0.08% for the chimpanzee-gorilla pair. These estimates, which were confirmed by additional data from GenBank, are substantially lower than previous ones, which included repetitive sequences and might have been based on less-accurate sequence data. The average sequence divergences between orangutans and humans, chimpanzees, and gorillas were 3.08% +/- 0.11%, 3.12% +/- 0.11%, and 3.09% +/- 0.11%, respectively, which also are substantially lower than previous estimates. The sequence divergences in other regions between hominoids were estimated from extensive data in GenBank and the literature, and Alus showed the highest divergence, followed in order by Y-linked noncoding regions, pseudogenes, autosomal intergenic regions, X-linked noncoding regions, synonymous sites, introns, and nonsynonymous sites. The neighbor-joining tree derived from the concatenated sequence of the 53 segments--24,234 bp in length--supports the Homo-Pan clade with a 100% bootstrap value. However, when each segment is analyzed separately, 22 of the 53 segments (approximately 42%) give a tree that is incongruent with the species tree, suggesting a large effective population size (N(e)) of the common ancestor of Homo and Pan. Indeed, a parsimony analysis of the 53 segments and 37 protein-coding genes leads to an estimate of N(e) = 52,000 to 96,000. As this estimate is 5 to 9 times larger than the long-term effective population size of humans (approximately 10,000) estimated from various genetic polymorphism data, the human lineage apparently had experienced a large reduction in effective population size after its separation from the chimpanzee lineage. Our analysis assumes a molecular clock, which is in fact supported by the sequence data used. Taking the orangutan speciation date as 12 to 16 million years ago, we obtain an estimate of 4.6 to 6.2 million years for the Homo-Pan divergence and an estimate of 6.2 to 8.4 million years for the gorilla speciation date, suggesting that the gorilla lineage branched off 1.6 to 2.2 million years earlier than did the human-chimpanzee divergence.

Great ape DNA sequences reveal a reduced diversity and an expansion in humans

The extent of DNA sequence variation of chimpanzees is several-fold greater than that of humans. It is unclear, however, if humans or chimpanzees are exceptional among primates in having low and high amounts of DNA sequence diversity, respectively. To address this, we have determined approximately 10,000 bp of noncoding DNA sequences at Xq13.3 (which has been extensively studied in both humans and chimpanzees) from 10 western lowland gorillas (Gorilla gorilla gorilla) and 1 mountain gorilla (Gorilla gorilla beringei; that is, from 2 of the 3 currently recognized gorilla subspecies), as well as 8 Bornean (Pongo pygmaeus pygmaeus) and 6 Sumatran (Pongo pygmaeus abelii) orang-utans, representing both currently recognized orang-utan subspecies. We show that humans differ from the great apes in having a low level of genetic variation and a signal of population expansion.

Global patterns of human DNA sequence variation in a 10-kb region on chromosome 1

Human DNA variation is currently a subject of intense research because of its importance for studying human origins, evolution, and demographic history and for association studies of complex diseases. A approximately 10-kb region on chromosome 1, which contains only four small exons (each <155 bp), was sequenced for 61 humans (20 Africans, 20 Asians, and 21 Europeans) and for 1 chimpanzee, 1 gorilla, and 1 orangutan. We found 52 polymorphic sites among the 122 human sequences and 382 variant sites among the human, chimpanzee, gorilla, and orangutan sequences. For the introns sequenced (8,991 bp), the nucleotide diversity (pi) was 0.058% among all sequences, 0.076% among the African sequences, 0.047% among the Asian sequences, and 0.045% among the European sequences. A compilation of data revealed that autosomal regions have, on average, the highest pi value (0.091%), X-linked regions have a somewhat lower pi value (0.079%), and Y-linked regions have a very low pi value (0.008%). The lower polymorphism in the present region may be due to a lower mutation rate and/or selection in the gene containing these introns or in genes linked to this region. The present region and two other 10-kb noncoding regions all show a strong excess of low-frequency variants, indicating a relatively recent population expansion. This region has a low mutation rate, which was estimated to be 0.74 x 10 per nucleotide per year. An average estimate of approximately 12,600 for the long-term effective population size was obtained using various methods; the estimate was not far from the commonly used value of 10,000. Fu and Li's tests rejected the assumption of an equilibrium neutral Wright-Fisher population, largely owing to the high proportion of low-frequency variants. The age of the most recent common ancestor of the sequences in our sample was estimated to be more than 1 Myr. Allowing for some unrealistic assumptions in the model, this estimate would still suggest an age of more than 500,000 years, providing further evidence for a genetic history of humans much more ancient than the emergence of modern humans. The fact that many unique variants exist in Europe and Asia also suggests a fairly long genetic history outside of Africa and argues against a complete replacement of all indigenous populations in Europe and Asia by a small Africa stock. Moreover, the ancient genetic history of humans indicates no severe bottleneck during the evolution of humans in the last half million years; otherwise, much of the ancient genetic history would have been lost during a severe bottleneck. We suggest that both the "Out of Africa" and the multiregional models are too simple to explain the evolution of modern humans.

Regions of low single-nucleotide polymorphism incidence in human and orangutan xq: deserts and recent coalescences

While scanning for single-nucleotide polymorphisms (SNPs) in the human Xq25-q28 region of CEPH families, we found six long "deserts" of low SNP incidence representing 28% of the investigated genome. One was 1.66 Mb in length. To determine whether these SNP deserts were due to reduced input of mutations or to recent coalescent events such as bottlenecks or selective sweeps, comparative sequence was determined from a female orangutan. The mean divergence was 2.9% and was not reduced in deserts compared with nondesert regions. Thus, the best explanation for the SNP deserts is recent coalescent events in humans. These events are the cause of substantial variation in human noncoding SNP incidence. In addition, the mutational spectrum in humans and orangutans was estimated as 63% AG (and CT), 17% AC (and GT), 8% CG, 4% AT, and 8% insertion/deletions. The average lifetime of a SNP destined to become fixed for a new allele between these species was estimated as 284,000 years.

Integrated analysis of sequence evolution and population history using hypervariable compound haplotypes

We have examined compound haplotypes from a highly informative region of human chromosome 16, in which information from the rapid evolution of a highly unstable minisatellite is integrated with data on the longer-term evolution of this segment from 10 flanking substitutional polymorphisms. Combined with sequence data from non-human primates, analysis of relationships between these compound haplotypes allows the reconstruction of a rooted network of the evolutionary pathways between them. Most relationships can be explained via simple substitutional mutations, although the origins of some haplotypes involve recurrent events at a hotspot for substitutional mutation and/or gene conversion. For compound haplotypes including the minisatellite array, the network found in a range of world-wide populations constitutes a highly informative data set for the analysis of population history (437 different compound haplotypes were discriminated among 658 studied). Since the mutation rates and processes of the minisatellite array are known from direct studies, ages for individual lineages have been estimated using associated minisatellite diversity. These analyses suggest that the higher information content and sampling depth of these compound haplotypes may allow more precise calibration of lineage ages than is possible using coalescent analysis of DNA sequence. Using this method we have dated the oldest Eurasian lineage as 52,000-66,000 years and the oldest European specific lineage as 37,600-56,200 years.

mtDNA sequence diversity of orangutans from the islands of Borneo and Sumatra

A comparison of mitochondrial DNA sequences was undertaken for two genes among orangutans from Borneo and Sumatra. The distribution of haplotypes among 42 individuals for NADH dehydrogenease subunit 3 and 39 individuals for cytochrome B was used to infer population structure. The haplotypes among all Bornean orangutans form a cluster of closely related individuals for both genes, with two distinct haplotypes occupying different regions of the island. Sumatran haplotypes fall into three distinct, and highly diverged, groups. Strikingly, one of the Sumatran haplotypes shares sequence identity with the most widespread Bornean haplotype. This haplotype distribution is considered in the context of the highly complex geological history for the area around the Malay Archipelago. Alternating periods of geographic isolation and reunion, resulting from glacially induced land bridge formation, presented substantial opportunity for population dispersal between periodically isolated demes. We present a paleodispersal model that is consistent with genetic, geological, paleoecological, and fossil data. The disparity of sequences present in orangutan populations argues against a simple Sumatra-Borneo dichotomy. Our evidence, and that of others, suggests that the Sumatran population alone contains the remnants of at least three separate lineages.